新流固流耦合算例报错

[[email protected] fluid]$ fe31
[[email protected] 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删除重新建立耦合面也不行

可以考虑下是不是耦合面的问题?提示的是

Segmentation fault (core dumped)

有可能是指针不对,或是数组超限了