[xye@localhost fluid]$ fe31
[xye@localhost fluid]$ fsiFoam.rheo
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.1-1dd681f6e943
Exec : fsiFoam.rheo
Date : Jun 19 2019
Time : 14:30:54
Host : localhost.localdomain
PID : 9607
CtrlDict : /home/xye/foam/foam-extend-3.1/etc/controlDict
Case : /var/tmp/tube.gk3/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: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading stress mesh
Create dynamic mesh2 for time = 0
--> 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: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}
Creating traction displacement boundary conditions
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: 13602
History point coordinates: (-0.010830839 0.00043966136 0.05956)
Reference point coordinates: (-0.011 0 0.06)
Time = 1e-05
Create extended GGI zone-to-zone interpolator
Checking fluid-to-solid face interpolator
Fluid-to-solid face interpolation error: 9.2188672e-05
Checking solid-to-fluid point interpolator (GGI)
Extended GGI, master point distance, max: 1.0600422e-05, avg: 2.6501053e-06, min: 0
Extended GGI, master point orientation (<0), max: -0.99813066, min: -1, nIncorrectPoints: 0/6480
Solid-to-fluid point interpolation error (GGI): 1.0600422e-05
Number of uncovered master faces: 0
Number of uncovered slave faces: 0
Create extended GGI zone-to-zone interpolator
Checking fluid2-to-solid face interpolator
Fluid2-to-solid face interpolation error: 0.00026333936
Checking solid-to-fluid2 point interpolator (GGI)
Extended GGI, master point distance, max: -4.9391652e-05, avg: -6.3485005e-05, min: -6.9931993e-05
Extended GGI, master point orientation (<0), max: -0.99786207, min: -1, nIncorrectPoints: 0/3240
Solid-to-fluid2 point interpolation error (GGI): 6.9931993e-05
Number of uncovered master faces: 0
Number of uncovered slave faces: 0
Time = 1e-05, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.01
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 cellMotionUz, 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
GAMG: Solving for cellMotionUz, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving flow model
Courant Number mean: 0.0023783633 max: 0.036268386 velocity magnitude: 3
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.1156212e-07, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 2.8372056e-07, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 7.9681438e-07, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 8.3190488e-07, No Iterations 26
GAMG: Solving for p, Initial residual = 0.0019855647, Final residual = 4.0332777e-07, No Iterations 11
time step continuity errors : sum local = 1.3145313e-08, global = 6.5857444e-10
GAMG: Solving for p, Initial residual = 0.018987894, Final residual = 9.9901421e-07, No Iterations 19
GAMG: Solving for p, Initial residual = 0.0015550554, Final residual = 6.5120828e-07, No Iterations 10
time step continuity errors : sum local = 1.6829818e-08, global = -6.6320371e-09
GAMG: Solving for p, Initial residual = 0.0034704502, Final residual = 8.5792427e-07, No Iterations 11
GAMG: Solving for p, Initial residual = 0.00042144851, Final residual = 9.4952659e-07, No Iterations 6
time step continuity errors : sum local = 2.4210026e-08, global = -1.0611199e-08
bounding k, min: 0 max: 0.03375 average: 0.03375
bounding epsilon, min: 0 max: 0.03421 average: 0.03421
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch downup. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch downup. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch downup. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch middleside. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch middleside. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch middleside. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch updown. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch updown. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch updown. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch wall. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch wall. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch wall. Evaluating as zeroGradient
DILUPBiCG: Solving for epsilon, Initial residual = 0.092836262, Final residual = 2.2445764e-07, No Iterations 3
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 1.8421419e-07, No Iterations 4
Setting traction on solid patch
Total force (fluid) = (0.03831704 -0.00094789465 0.13955462)
Total force (solid) = (1956.2512 -4007.049 -0.079763461)
Time = 1e-05, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.01
Creating accumulated fluid2 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 cellMotionUz, 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
GAMG: Solving for cellMotionUz, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving flow model
Courant Number mean: 0.0047042481 max: 0.055231699 velocity magnitude: 5.14
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 2.2245216e-07, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 2.1000714e-07, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.6929291e-08, No Iterations 4
GAMG: Solving for p, Initial residual = 1, Final residual = 5.8415859e-07, No Iterations 18
GAMG: Solving for p, Initial residual = 0.0042381369, Final residual = 5.9645218e-07, No Iterations 8
time step continuity errors : sum local = 5.1677233e-09, global = -9.5521909e-10
GAMG: Solving for p, Initial residual = 0.23149673, Final residual = 6.4477968e-07, No Iterations 20
GAMG: Solving for p, Initial residual = 0.064111542, Final residual = 4.0324494e-07, No Iterations 12
time step continuity errors : sum local = 2.1015578e-10, global = 2.5631328e-11
GAMG: Solving for p, Initial residual = 0.81432421, Final residual = 8.9817477e-07, No Iterations 19
GAMG: Solving for p, Initial residual = 0.037963488, Final residual = 7.9238394e-07, No Iterations 15
time step continuity errors : sum local = 8.4133791e-11, global = -2.6803614e-11
bounding k, min: 0 max: 0.09907 average: 0.09907
bounding epsilon, min: 0 max: 0.172 average: 0.172
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch interfaceWater. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch interfaceWater. Evaluating as zeroGradient
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch interfaceWater. Evaluating as zeroGradient
DILUPBiCG: Solving for epsilon, Initial residual = 0.095675524, Final residual = 2.7053405e-08, No Iterations 3
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 1.7321655e-07, No Iterations 4
Setting traction on solid patch
Total force (fluid2) = (9.5013626e-06 -2.1033763e-05 0.27801767)
Total force (solid) = (-9.4576616e-06 2.0935566e-05 -0.27672832)
Evolving stress model
DICPCG: Solving for D, Initial residula = 1, Final residual = 0.0083094099, No outer iterations = 98
Max relative residual = 1, Relative residual = 0.00099513769, enforceLinear = 0
Segmentation fault (core dumped)
又是这个Segmentation fault (core dumped)，这回我把VTK删除重新建立耦合面也不行