各位大佬，我参考文献 3rd High-Lift Workshop Summary Paper - OpenFOAM, STAR-CCM+ & LAVA simulations on Unstructured Grids链接https://arc.aiaa.org/doi/10.2514/6.2018-1253。用rhoPimpleFoam计算飞机，但是算了几步就发散了，希望各位大佬帮忙看看是fvSchemes和fvSolution的问题嘛，还是nut、nuTilda设置的问题。
目前用的是作者上传到会议网站的网格，按照文章里说的，采用的是SA模型，以及动量和湍流项采用的是二阶迎风。
Courant Number mean: 1.4881475e-08 max: 1.8237221
Time = 3.2e-09
*
* PIMPLE: Iteration 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 2.421018e-11, Final residual = 2.421018e-11, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 2.1538356e-10, Final residual = 2.1538356e-10, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 2.0979924e-10, Final residual = 2.0979924e-10, No Iterations 0
DILUPBiCGStab: Solving for e, Initial residual = 4.1760218e-08, Final residual = 2.1963671e-13, No Iterations 1
GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.0250634e-11, global = 5.0210593e-11, cumulative = 3.7505015e-09
GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.0250633e-11, global = 5.0210593e-11, cumulative = 3.8007121e-09
PIMPLE: Iteration 2
DILUPBiCGStab: Solving for Ux, Initial residual = 1.3321419e-08, Final residual = 1.002751e-14, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 6.2117311e-09, Final residual = 6.2117311e-09, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 1.4225061e-10, Final residual = 1.4225061e-10, No Iterations 0
DILUPBiCGStab: Solving for e, Initial residual = 1.5129023e-07, Final residual = 2.4037934e-13, No Iterations 1
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 3.900971e-09
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 4.00123e-09
PIMPLE: Iteration 3
DILUPBiCGStab: Solving for Ux, Initial residual = 2.4133e-08, Final residual = 2.8587861e-14, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1.1420948e-08, Final residual = 1.9949995e-14, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1.4399722e-10, Final residual = 1.4399722e-10, No Iterations 0
DILUPBiCGStab: Solving for e, Initial residual = 1.3956424e-07, Final residual = 2.3199397e-13, No Iterations 1
[73] #0 Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found
addr2line failed
[73] #1 Foam::sigFpe::sigHandler(int)sh: addr2line: command not found
addr2line failed
[73] #2 ?sh: addr2line: command not found
addr2line failed
[73] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate()sh: addr2line: command not found
addr2line failed
[73] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: addr2line: command not found
addr2line failed
[73] #5 ?sh: addr2line: command not found
addr2line failed
[73] #6 __libc_start_mainsh: addr2line: command not found
addr2line failed
[73] #7 ?sh: addr2line: command not found
addr2line failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x3c75cd0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=22, req_array=0x3966230, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=52, req_array=0x304c450, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=46, req_array=0x355b680, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=36, req_array=0x3dfabb0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x349e540, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=26, req_array=0x3bd5ce0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=12, req_array=0x3f18790, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=18, req_array=0x3f97b30, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=14, req_array=0x345ac20, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=32, req_array=0x3df8960, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x4316400, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x2fe0f30, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=50, req_array=0x3a36c60, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x3f6c5f0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
Fatal error in PMPI_Waitall: Unknown error class, error stack:
PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x520edf0, status_array=0x1) failed
MPIR_Waitall_impl(181)................: fail failed
MPIDU_Complete_posted_with_error(1710): Process failed
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
fvSchmes
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
limited cellLimited Gauss linear 1;
grad(U) $limited;
grad(k) $limited;
grad(omega) $limited;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);//Gauss linearUpwind limited;
turbulence Gauss linearUpwind default;//Gauss linearUpwind limited;
energy Gauss linearUpwind default;//Gauss linearUpwind limited;
div(phi,k) $turbulence;
div(phi,omega) $turbulence;
div(phi,e) $energy;
div(phi,K) $energy;
div(phi,Ekp) $energy;
div(phi,nuTilda) $energy;
div(phiv,p) Gauss upwind;
div((phi|interpolate(rho)),p) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited 0.5;//Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.5;//corrected;
}
wallDist
{
method meshWave;
}
fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
smoother GaussSeidel;
tolerance 1e-5;
relTol 0.01;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
ncellsInCoarsetLevel 1000;
mergeLevels 1;
}
pFinal
{
$p;
relTol 0;
}
"(rho|U|k|omega|e)"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-8;
relTol 0.0;
}
"(rho|U|k|omega|e|nuTilda)Final"
{
$U;
relTol 0;
}
}
SIMPLE
{
residualControl
{
p 1e-5;
U 1e-5;
"(k|omega|e)" 1e-5;
}
nNonOrthogonalCorrectors 3;
pMinFactor 0.1;
pMaxFactor 2;
}
PIMPLE
{
nCorrectors 2;
nNonOrthogonalCorrectors 3;
nOuterCorrectors 3;
pMinFactor 0.1;
pMaxFactor 2;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
omega 0.7;
}
// ************************************************************************* //
nut
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
freestream
{
type calculated;
value uniform 0;
}
wall
{
type nutkWallFunction;
type nutkWallFunction;
value uniform 0;
}
#includeEtc "caseDicts/setConstraintTypes"
}
// ************************************************************************* //
nuTilda
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
freestream
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}
wall
{
type fixedValue;
value uniform 0;
}
symmetry
{
type symmetry;
}
#includeEtc "caseDicts/setConstraintTypes"
}
// ************************************************************************* //
