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.
- OpenFOAM (pimpleFoam). In case you are using a very old OpenFOAM version, you will need to adjust the solver to
Fluid/system/controlDictfile. For more information, have a look at the OpenFOAM adapter documentation.
CalculiX. For more information, have a look at the CalculiX adapter documentation.
FEniCS. The structural model is currently limited to linear elasticity. Currently 3D functionality is experimental in the FEniCS adapter and more details can be found here For more information, have a look at the FeniCS adapter documentation.
Running the simulation
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: