计算情况大致如下，管内流体设置为湍流，不过速度 k 等参数都为0
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: 39056
History point coordinates: (0.044069226 0.00043439053 0.024247491)
Reference point coordinates: (0 0.005 0.025)
Time = 0.00025
Create extended GGI zone-to-zone interpolator
Checking fluid-to-solid face interpolator
Fluid-to-solid face interpolation error: 0.0010269982
Checking solid-to-fluid point interpolator (GGI)
Extended GGI, master point distance, max: 2.1053987e-05, avg: 1.3178253e-05, min: -8.2349703e-07
Extended GGI, master point orientation (<0), max: -0.99512002, min: -0.99999995, nIncorrectPoints: 0/2880
Solid-to-fluid point interpolation error (GGI): 2.1053987e-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: 5.1311848e-06
Checking solid-to-fluid2 point interpolator (GGI)
Extended GGI, master point distance, max: 7.9914926e-10, avg: -2.9176373e-06, min: -5.286356e-06
Extended GGI, master point orientation (<0), max: -0.99600636, min: -0.99720899, nIncorrectPoints: 0/10800
Solid-to-fluid2 point interpolation error (GGI): 5.286356e-06
Number of uncovered master faces: 0
Number of uncovered slave faces: 0
Time = 0.00025, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.6
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.080258043 max: 0.44144684 velocity magnitude: 3.0178771
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 6.4947932e-08, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 7.0788482e-08, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 9.6953364e-08, No Iterations 5
GAMG: Solving for p, Initial residual = 1, Final residual = 9.6088419e-07, No Iterations 15
GAMG: Solving for p, Initial residual = 0.10051973, Final residual = 8.7262466e-07, No Iterations 10
time step continuity errors : sum local = 3.3363082e-08, global = 1.9018035e-09
GAMG: Solving for p, Initial residual = 0.15929902, Final residual = 8.8099237e-07, No Iterations 14
GAMG: Solving for p, Initial residual = 0.039798, Final residual = 8.6682036e-07, No Iterations 9
time step continuity errors : sum local = 2.720883e-08, global = -4.1907837e-09
GAMG: Solving for p, Initial residual = 0.020351497, Final residual = 9.4198148e-07, No Iterations 9
GAMG: Solving for p, Initial residual = 0.0017079998, Final residual = 7.0146151e-07, No Iterations 6
time step continuity errors : sum local = 2.0868976e-08, global = -1.182766e-09
bounding k, min: 0 max: 0.03375 average: 0.03375
bounding epsilon, min: 0 max: 0.02547 average: 0.02547
From function void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
in file derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C at line 174
Cannot access RASModel::G field for patch INTERFACE. 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 INTERFACE. 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 INTERFACE. Evaluating as zeroGradient
DILUPBiCG: Solving for epsilon, Initial residual = 0.53181342, Final residual = 7.5558784e-08, No Iterations 4
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 3.4610577e-07, No Iterations 4
Setting traction on solid patch
Total force (fluid) = (-0.17594396 557.35816 0.00022212502)
Total force (solid) = (0.18797305 -559.54374 -0.00022353041)
Time = 0.00025, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.6
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 max: 0 velocity magnitude: 0
DILUPBiCG: Solving for Ux, Initial residual = 0.98174913, Final residual = 3.7927751e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.9723697, Final residual = 4.0377006e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.99617218, Final residual = 3.8146891e-07, No Iterations 1
GAMG: Solving for p, Initial residual = 0.99810585, Final residual = 9.2691854e-07, No Iterations 83
GAMG: Solving for p, Initial residual = 0.041756438, Final residual = 7.1080034e-07, No Iterations 8
time step continuity errors : sum local = 5.5955825e-23, global = -1.1884079e-23
GAMG: Solving for p, Initial residual = 0.0036700866, Final residual = 6.647089e-07, No Iterations 7
GAMG: Solving for p, Initial residual = 0.00067284946, Final residual = 5.8457543e-07, No Iterations 5
time step continuity errors : sum local = 4.7407162e-23, global = 2.6456577e-24
GAMG: Solving for p, Initial residual = 0.00032557593, Final residual = 7.899604e-07, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00028101735, Final residual = 7.7354747e-07, No Iterations 4
time step continuity errors : sum local = 6.2751512e-23, global = 4.6937306e-25
bounding k, min: 0 max: 0 average: 0
bounding epsilon, min: 0 max: 0 average: 0
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
Floating point exception (core dumped)
不过，把内部流体的湍流关闭后，却又可以计算了