This solver is only a demo. It should not be used to simulate the real world
Ah. And by the way: it does not work in parallel
At each time-step the following things are done:
- the equations in the fluid region are solved
- fluid pressures at the fluid/solid-interface are written as boundary equations to the solid
- the equations for the solid are solved
- displacements for the interface-nodes of the solid region are written to the interface-nodes of the fluid region
- the grid of the fluid region is deformed
In order for blockMesh to work on the grid regions you've got to fool it:
ln -s region1/polyMesh/ icoStructFoamTest/constant/polyMesh blockMesh . icoStructFoamTest rm icoStructFoamTest/constant/polyMesh ln -s region2/polyMesh/ icoStructFoamTest/constant/polyMesh blockMesh . icoStructFoamTest rm icoStructFoamTest/constant/polyMesh
foamToVTK . icoStructFoamTest -mesh region1 foamToVTK . icoStructFoamTest -mesh region2
An alternative approach is described in How to postprocess multiple regions.
The fluid mesh is deformed. The solid mesh has to be deformed while post-processing (with the Warp (vector) filter)
4 Example cases
Please note that none of these cases resembles anything from the real world and only server for demonstration purposes
4.1 Deformable channel
This case simulates flow in a channel. One wall of the channel is coupled to a deformable solid.
The two regions are:
- Lower Region: channel with velocity inlet on the left and pressure outlet on the right
- Upper Region: solid which is fixed at the top and free at the left and the right
In the beginning both of them are of equal size and square.
Due to the higher pressure at the inlet the solid gets deformed.
Both pictures are from the end of the simulation. For the solid the strain is shown. For the fluid the left picture the colors show the pressure, the right picture shows the velocity.
4.2 Soft thing on a stick
This case simulates the flow through a channel with an obstacle that consists of two parts: a solid, fixed part around which a deformable material is wrapped.
The two regions are shown in the picture on the left:
- Blue: the fluid. On the left-side a velocity inlet, on the right side a pressure outlet. Top is a wall, botto, a symmetry boundary.
- Red: the solid. The bottom edges are symmetry boundaries. The inner edges are fixed (the stick in the title)
The picture on the left is the undeformed geometry at the beginning. The picture on the right is the geometry at the end of the simulation. The solid is colored with the strain, the fluid with the pressure.
In certain parts of Austria this is also known as the Lentos-case.
5.1 OpenFOAM 1.3
These files are meant to be used with :
5.2 OpenFOAM 1.5-dev
It's available at the following SVN repository folder: trunk/Breeder_1.5/solvers/other/IcoStructFoam @ sourceforge.net/openfoam-extend
Instructions on how to use it are given in the README file.
- Download can be done by using SVN:
svn checkout svn://svn.code.sf.net/p/openfoam-extend/svn/trunk/Breeder_1.5/solvers/other/IcoStructFoam IcoStructFoam
- Download using Git from the Unofficial-Extend-Project-Mirror:
git clone https://github.com/Unofficial-Extend-Project-Mirror/openfoam-extend-Breeder1.5-solvers-other-IcoStructFoam.git IcoStructFoam
- 11 Nov 2005: Initial upload
- 10 Nov 2006: Changed Version for 1.3
- 7 Oct 2008: Added version for 1.5-dev, but only on the SVN repository
- 11 Oct 2014: Added instructions for getting version for 1.5-dev from the SVN repository and unofficial Git mirror