[hxaxtma]

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:

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"

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

[hxaxtma]

no one has any idea?

[hxaxtma]

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?

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

[stephie]

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

best regards,
Stephie

[songchang]

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.

[thegauravonline]

Originally Posted by songchang

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.

Hi song

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

[popa]

Originally Posted by thegauravonline

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!

[Hgholami]

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

[popa]

Originally Posted by Hgholami

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.

[Hgholami]

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

[popa]

Originally Posted by Hgholami

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.

[Hgholami]

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

[chandra shekhar pant]

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!