Tutorial for an FSI simulation of a three-dimensional expanding tube scenario

## Setup

The expanding tube test case involves a cylindrical fluid domain surrounded by a solid domain. A pressure inlet boundary condition is applied at the inlet for 3 milliseconds, and then 0 set to zero for a further 7 millisecond. The pressure of the fluid expands the tube which then relaxes once the pressure decreases.

The expanding tube test case comes with the interface surface mesh connectivity of the solid domain. This allows the use of nearest-projection mapping of the displacements of the solid domain. In order to run the example with nearest projection mapping, the “node-mesh-with-connectivity” has been specified in the solid-calculix/config.yml file. More details can be found in the CalculiX configuration description.

## Available solvers

Fluid participant:

• OpenFOAM (pimpleFoam). In case you are using a very old OpenFOAM version, you will need to adjust the solver to pimpleDyMFoam in the Fluid/system/controlDict file. For more information, have a look at the OpenFOAM adapter documentation.

Solid participant:

You can start the simulation by running the script ./run.sh located in each participant directory. OpenFOAM can be executed in parallel using run.sh -parallel. The default setting uses 4 MPI ranks.
You can visualize the results using paraView or cgx(for native CalculiX result files), as usual. The total deformation is rather small. Multiplying the deformation by factor of 10 (warp by vector filter in paraView) and visualizing the fluid domain at t=0.005s looks as follows:
You can also plot the displacement of the midpoint of the tube by running sh plot-displacement.sh <filename>. The displacement plot for each solver combination looks like: