Difference between revisions of "Contrib icoStructFoam"

From OpenFOAMWiki
m
 
m
Line 38: Line 38:
  
 
=== Deformable channel ===
 
=== Deformable channel ===
 +
 +
[[Image:icoStructFoam_testDeformed.png|left|thumb|Pressure and strain]]
 +
 +
[[Image:icoStructFoam_testDeformed.png|right|thumb|Velocity and strain]]
  
 
=== Soft thing on a stick ===
 
=== Soft thing on a stick ===
 +
 +
[[Image:icoStructFoam_thingOnStickGeometry.png|left|thumb|The two regions]]
 +
 +
[[Image:icoStructFoam_thingOnStickDeformed.png|right|thumb|Pressure and strain]]
  
 
== Download ==
 
== Download ==
 +
 +
[[Media:icoStructFoam.tar.gz|The solver]]
 +
 +
[[Media:icoStructFoam_icoStructFoamTest.tar.gz|Case of the deformable channel]]
 +
 +
[[Media:icoStructFoam_thingOnStick.tar.gz|Case of the soft thing on a stick]]
  
 
--[[User:Bgschaid|Bgschaid]] 13:34, 11 Nov 2005 (CET)
 
--[[User:Bgschaid|Bgschaid]] 13:34, 11 Nov 2005 (CET)

Revision as of 13:26, 11 November 2005

UNDER CONSTRUCTION

1 Description

This solver is a hybrid of icoFoam and stressedFoam. One grid-region is an imcompressible fluid, the other is a solid. The pressure of the fluid is used to deform the solid.

2 Algorithm

At each time-step the following things are done:

  1. the equations in the fluid region are solved
  2. fluid pressures at the fluid/solid-interface are written as boundary equations to the solid
  3. the equations for the solid are solved
  4. displacements for the interface-nodes of the solid region are written to the interface-nodes of the fluid region
  5. the grid of the fluid region is deformed

3 Usage

Due to the two grids the standard Utilities for pre- and post-processing (FoamX and paraFoam) can't be used. Everthing is setup with the help of a text-editor.

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

To write the data of the two regions for paraview just use foamToVTK:

foamToVTK . icoStructFoamTest -mesh region1
foamToVTK . icoStructFoamTest -mesh region2

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

Pressure and strain
Velocity and strain

4.2 Soft thing on a stick

The two regions
Pressure and strain

5 Download

The solver

Case of the deformable channel

Case of the soft thing on a stick

--Bgschaid 13:34, 11 Nov 2005 (CET)