Error for solving FSI problem by using fsiFoam, thanks.
 

Hi,

Currently I am using fsiFoam of extend4.0 to run a standard turek&horn cylinder-foil fluid structures interaction case. I compiled whole extend and Toolkits step by step follow http://www.openfoamwiki.net/index.ph...re_interaction.

The solidFoam can run, but fsiFoam was not able to run. Error information is as follows:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 4.0-10298020aac7
Exec : fsiFoam
Date : Mar 23 2018
Time : 10:36:11
Host : workstation
PID : 29565
CtrlDict : "/home/liyanmiho/foam/liyanmiho-4.0/FluidSolidInteraction/run/fsiFoam/HronTurekFsi3/fluid/system/controlDict"
Case : /home/liyanmiho/foam/liyanmiho-4.0/FluidSolidInteraction/run/fsiFoam/HronTurekFsi3/fluid
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: refVelocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading stress mesh
Creating fixed displacement boundary condition
Creating traction displacement boundary conditions
Selecting rheology model linearElastic
Creating constitutive model
Selecting coupling scheme IQN-ILS

Starting time loop

Creating pointHistory function object.
History point ID: 1165
History point coordinates: (0.6 0.2 0.015)
Reference point coordinates: (0.6 0.2 0.025334)
Creating hronTurekReport function object.
Time = 0.001 (dt = 0.001)

Create extended GGI zone-to-zone interpolator
Checking fluid-to-solid face interpolator
Fluid-to-solid face interpolation error: 0.00333563
Checking solid-to-fluid point interpolator (GGI)
calcMasterPointAddressing() const
Extended GGI, master point distance, max: 0, avg: 0, min: 0
Solid-to-fluid point interpolation error (GGI): 1.24175e-16
Number of uncovered master faces: 0
Number of uncovered slave faces: 0
calcSlavePointAddressing() const
Extended GGI, slave point distance, max: 0, avg: 0, min: 0
Extended GGI, slave point orientation (<0), max: -0.707107, min: -1, nIncorrectPoints: 0/434
Setting traction on solid patch
Total force (fluid) = (0 0 0)
Total force (solid) = (0 0 0)
Evolving solid solver: unsTotalLagrangianSolid
fvMatrix<Type>::solve: Solving for D, Initial residula = 0, Final residual = 0, No outer iterations = 1
Max relative residual = 0, Relative residual = 0, enforceLinear = 0
Current fsi relative residual norm: 0
Alternative fsi residual: 0

Time = 0.001, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.05
Creating accumulated fluid interface displacement
Maximal accumulated displacement of interface points: 0
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving fluid solver: consistentIcoFluid
Courant Number mean: 0 max: 0 velocity magnitude: 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.2281e-08, No Iterations 2
DILUPBiCG: Solving for Uy: solution singularity
GAMG: Solving for p, Initial residual = 1, Final residual = 8.62349e-07, No Iterations 74
GAMG: Solving for p, Initial residual = 0.035702, Final residual = 9.96795e-07, No Iterations 19


--> FOAM FATAL ERROR:
flux requested but p not specified in the fluxRequired sub-dictionary of fvSchemes.

From function fvMatrix<Type>::flux()
in file lnInclude/fvMatrix.C at line 1046.

FOAM aborting

Aborted (core dumped)

Could you please give me valuable suggestion or information to fix this problem?

Thank you very much!

Liyan

Hi Liyan,

I realize this may be two months too late, but to fix the flux required problem, you need to add

to the flow solver (eg consistentIcoFluid). I added it right after mesh is defined in the evolve() function. This should fix the issue!

Hi amanbearpig
I run HronTurek test with fsiFoam (FSI2). the problem is when deflection increase, the sulotion diverge. Do you any recommend for fvScheme ond fvSolution?