[phanina]

Dear Philip and all,

I'm trying to modify and model a fluid-structure interaction case using solids4foam toolbox.I took Turek Hron FSI tutorial case and tried to modify it to conform to my probelm(see attached).I have circular solid domain which is supposed to deform/move as a result of applied fluid velocity and pressure.However, in my simulation the solid domain moves for some time and suddenly stops.I'm not sure what could be the cause of this problem and how to diagnose it.

I'd really appreciate any suggestion in finding solution for the problem

below is the output of log file just before the simulation stops:





--> FOAM Warning :
From function void neoHookeanElastic::correct(surfaceSymmTensorField& sigma)
in file materialModels/mechanicalModel/mechanicalLaws/mechanicalLaw/mechanicalLaw.C at line 549
Material linearity enforced for stability!
--> FOAM Warning :
From function void neoHookeanElastic::correct(surfaceSymmTensorField& sigma)
in file materialModels/mechanicalModel/mechanicalLaws/mechanicalLaw/mechanicalLaw.C at line 549
Material linearity enforced for stability!
--> FOAM Warning :
From function void neoHookeanElastic::correct(surfaceSymmTensorField& sigma)
in file materialModels/mechanicalModel/mechanicalLaws/mechanicalLaw/mechanicalLaw.C at line 549
Material linearity enforced for stability!
--> FOAM Warning :
From function void neoHookeanElastic::correct(surfaceSymmTensorField& sigma)
in file materialModels/mechanicalModel/mechanicalLaws/mechanicalLaw/mechanicalLaw.C at line 549
Material linearity enforced for stability!
--> FOAM Warning :
From function void neoHookeanElastic::correct(volSymmTensorField& sigma)
in file materialModels/mechanicalModel/mechanicalLaws/mechanicalLaw/mechanicalLaw.C at line 437
Material linearity enforced for stability!
FDICPCG: Solving for D, Initial residual = 1, Final residual = 8.86576e-06, No outer iterations = 1000
Max relative residual = 1.54678, Relative residual = 2.82707e-05, enforceLinear = true
Interpolating point values using GGI
Interpolating point values using GGI
FSI relative residual1 norm for interface 0: 1
FSI residual2 norm for interface 0: 1

Time = 2.715, iteration: 27
Maximal accumulated displacement of interface 0: 2.36062e+13
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 5.81041e-07, No Iterations 5
GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 7.28632e-07, No Iterations 5
GAMG: Solving for cellMotionUx, Initial residual = 0.0438151, Final residual = 3.06917e-07, No Iterations 5
GAMG: Solving for cellMotionUy, Initial residual = 0.0400751, Final residual = 3.38384e-07, No Iterations 5
Evolving fluid model: icoFluid
Courant Number mean: 291.129 max: 35487.5 velocity magnitude: 5.96034e+14
DILUPBiCG: Solving for Ux, Initial residual = 0.999993, Final residual = 5.64252e-07, No Iterations 201
DILUPBiCG: Solving for Uy, Initial residual = 0.999997, Final residual = 4.28668e-07, No Iterations 221