CFD Online Discussion Forums - icoFsiElasticNonLinULSolidFoam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - icoFsiElasticNonLinULSolidFoam - simulation does not start (https://www.cfd-online.com/Forums/openfoam-solving/130508-icofsielasticnonlinulsolidfoam-simulation-does-not-start.html)

icoFsiElasticNonLinULSolidFoam - simulation does not start

Hi guys,

I am trying to setup an icoElasticFsi...simulation. Therefore I copied the HronTurek TestCase from foam-extend3.0 and changed the settings to my 3D Mesh settings.

My testcase is a simple finite flap in an couette channel.

My fluid boundaries are:
in, out (cyclicGgi), top (moving wall), bottom (no slip wall) and the both channelsides (slip walls). The flap itself ist called "cyl"

My solid boundaries are:
cyl_fixed and cyl_solid -> fixed is the ground patch of the flap and solid is the moving part.

When I run icoFsiEl....I got following error immediatly:

Code:

Create time





Create dynamic mesh for time = 0



 --> FOAM Warning : 

    From function dlLibraryTable::open(const fileName& functionLibName)

    in file db/dlLibraryTable/dlLibraryTable.C at line 85

    could not load libmesquiteMotionSolver.so: cannot open shared object file: No such file or directory

--> FOAM Warning : 

    From function dynamicFvMesh::New(const IOobject& io)

    in file dynamicFvMesh/newDynamicFvMesh.C at line 42

    Loading of dynamic mesh library libmesquiteMotionSolver.so unsuccesful. Some dynamic mesh  methods may not be  available

Selecting dynamicFvMesh dynamicMotionSolverFvMesh

Initializing the GGI interpolator between master/shadow patches: in/out

Selecting motion solver: laplace

Selecting motion diffusivity: quadratic





Reading transportProperties





Reading field p





Reading field U





Reading/calculating face flux field phi





Reading incremental displacement field DU





Patch cyl_solid    Traction boundary field: DU

    nonLinear set to updated Lagrangian

--> FOAM Warning : The gradScheme for DU should be "extendedLeastSquares 0" for the boundary non-orthogonal correction to be right

Reading incremental displacement field DV





Reading accumulated velocity field V





Reading accumulated stress field sigma





Reading incremental stress field DSigma





Selecting rheology model linearElastic

Creating constitutive model

Floating point exception

All my testfiles are included in the attachement (also boundary file, except faces, zones, points, sets, eg etc...)

Maybe you can help me, cause the error log is not given that much information.

In the original Turek Totorial would be the next step after "Creating constitutive model"

Code:

Reading coupling properties

Create fluid-to-solid and solid-to-fluid interpolators

Check fluid-to-solid and solid-to-fluid interpolators

Fluid-to-solid face interpolation error: 1.24127e-16

Solid-to-fluid face interpolation error: 1.24127e-16

So there must be somehting wrong with my coupling, but I cannot find the failure.

Thank you for your help

no one has any idea?

Ok simple Error, but I forgot to set TwoDMotion to no: Well now I got immediatly a floating point error after 1 interation. Any idea why?

Code:

Create time



Create dynamic mesh for time = 0



--> FOAM Warning : 

    From function dlLibraryTable::open(const fileName& functionLibName)

    in file db/dlLibraryTable/dlLibraryTable.C at line 85

    could not load libmesquiteMotionSolver.so: cannot open shared object file: No such file or directory

--> FOAM Warning : 

    From function dynamicFvMesh::New(const IOobject& io)

    in file dynamicFvMesh/newDynamicFvMesh.C at line 42

    Loading of dynamic mesh library libmesquiteMotionSolver.so unsuccesful. Some dynamic mesh  methods may not be  available

Selecting dynamicFvMesh dynamicMotionSolverFvMesh

Selecting motion solver: laplace

Selecting motion diffusivity: quadratic



Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Reading incremental displacement field DU

Patch plate    Traction boundary field: DU

    nonLinear set to updated Lagrangian

--> FOAM Warning : The gradScheme for DU should be "extendedLeastSquares 0" for the boundary non-orthogonal correction to be right

Reading incremental displacement field DV

Reading accumulated velocity field V

Reading accumulated stress field sigma

Reading incremental stress field DSigma



Selecting rheology model linearElastic

Creating constitutive model

Reading coupling properties

Create fluid-to-solid and solid-to-fluid interpolators

Check fluid-to-solid and solid-to-fluid interpolators

Fluid-to-solid face interpolation error: 1.78415e-15

Solid-to-fluid face interpolation error: 1.78415e-15



Starting time loop

Time = 1e-05

Selecting coupling scheme Aitken

Time = 1e-05, iteration: 1

Current fsi under-relaxation factor: 0.01

Maximal accumulated displacement of interface points: 0

Courant Number mean: 0 max: 0 velocity magnitude: 0

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 8.89011e-14, No Iterations 1

DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 4.03802e-11, No Iterations 1

DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 1.98383e-11, No Iterations 1

Floating point exception

Hey,
I get now the same mistake. After 20 seconds the simultaions allways interrupts.
Did you find any solution?

best regards,
Stephie

How to use icoFsiElasticNonLinULSolidFoam? Looking forward to a testcase.

Hey,

I am trying to setup an icoElasticFsi...simulation. Glad to hear that the solver icoFsiElasticNonLinULSolidFoam can make it. Anybody can give me detailed guidance on the HronTurek TestCase using openfoam 4.0？ The necessary code file and Instruction about the entire process will be more useful.

Looking forward to your reply.

Quote:

Originally Posted by songchang (Post 663177)

Hi song

I installed Foam extend 3.2 in order to work on icoFsiElasticNonLinULSolidFoam solver.

Quote:

Originally Posted by thegauravonline (Post 664933)

Hi song

I installed Foam extend 3.2 in order to work on icoFsiElasticNonLinULSolidFoam solver.

hello,can you get the result that the solid can swing up and down? I run the case which default,but it stay static.why!

I change the materials properties as the HronTurek's study (FSI2) and the structure start to oscillate.

Quote:

Originally Posted by Hgholami (Post 722322)

I change the materials properties as the HronTurek's study (FSI2) and the structure start to oscillate.

Yeah,i use foam-extend4.0,and i also change the materials properties(E),but i can't get a stable oscillate,do you have the same situation? or

can you provided what properties did you change? Thank you.

Dear Popa
I uses FSI2 of hronTurek study. the stable simulation achieve with fsiFoam

Quote:

Originally Posted by Hgholami (Post 728410)

Dear Popa
I uses FSI2 of hronTurek study. the stable simulation achieve with fsiFoam

What version of openFoam do you use?I haven't the solver fsiFoam in form_extend4.0,and once,i did use a tool(http://openfoamwiki.net/index.php/Ex...re_interaction) and get the result and its solver is just the fsiFoam. Do you do so or other ways? Thank you.

YES. I uses toolkit of openFoam-ext4.0 and then fsiFoam

Dear All,

In this regard I have a query, I already have OpenFoam 19.06, do I need to install OpenFoam-extend-4.0? If I install OpenFoam-extend-4.0 in the same workstation, will the mess up ? My ultimate aim is to install OpenFPCI, which essentially requires "both OpenFOAM-Extend and ParaFEM are compiled on the system, using the Open-MPI available on the system".
Any help/suggestion is welcomed.

Thanks!