pointwise输出二维网格到openfoam无法计算

[chengan.wang]

大家好，我用pointwise输出了一个翼型的网格，使用checkMesh检查没有错误

Mesh stats
    points:           28336
    internal points:  0
    faces:            56210
    internal faces:   27874
    cells:            14014
    faces per cell:   6
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     14014
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    airFoil2D           154      308      ok (non-closed singly connected)  
    farfield            154      308      ok (non-closed singly connected)  
    frontAndBack        28028    28336    ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-4.99867 -5.67032 0) (6.78888 5.66637 3)
    Mesh (non-empty, non-wedge) directions (1 1 0)
    Mesh (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (-1.85094e-17 2.46792e-17 -2.22841e-16) OK.
    Max cell openness = 4.12267e-15 OK.
    Max aspect ratio = 170.167 OK.
    Minimum face area = 4.43442e-08. Maximum face area = 1.67744.  Face area magnitudes OK.
    Min volume = 1.33032e-07. Max volume = 0.378026.  Total volume = 313.461.  Cell volumes OK.
    Mesh non-orthogonality Max: 17.8429 average: 1.07748
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.9585 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

我用paraFoam打开也没有问题

但使用simpleFoam运算就会出现问题：

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 156.2, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 4.72215, No Iterations 1000
GAMG:  Solving for p, Initial residual = 1, Final residual = 2.66928e+147, No Iterations 1000
time step continuity errors : sum local = 1.37947e+146, global = 9.17194e+129, cumulative = 9.17194e+129
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve() at ??:?
#9  Foam::incompressible::RASModels::SpalartAllmaras::correct() at ??:?
#10  ? at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  ? at ??:?
Floating point exception (core dumped)

我用Construct2D或者snappyHexmesh生成的网格却没有问题，说明求解器设置没有问题，所以错误出自哪里呢？
谢谢各位！

[李东岳]

这个边界条件设置的有问题，速度压力都没收敛

[zzcfd]

最后解决了吗。我现在还是会遇到这个问题，改变了边界条件还是这样的计算结果。
Create time

Create mesh for time = 0

SIMPLE: convergence criteria
field p tolerance 1e-05
field U tolerance 1e-05
field nuTilda tolerance 1e-05

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model SpalartAllmaras
Selecting patchDistMethod meshWave
RAS
{
RASModel SpalartAllmaras;
turbulence on;
printCoeffs on;
sigmaNut 0.666667;
kappa 0.41;
Cb1 0.1355;
Cb2 0.622;
Cw2 0.3;
Cw3 2;
Cv1 7.1;
Cs 0.3;
}

No MRF models present

No finite volume options present

Starting time loop

momentumError momErr:
turbulenceFields turbulenceFields1: storing fields:
turbulenceProperties:nut
turbulenceProperties:nuEff
turbulenceProperties:k
turbulenceProperties:R
turbulenceProperties:I
turbulenceProperties:L
turbulenceProperties:devReff
turbulenceProperties:epsilon

Time = 0.01

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00113906, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.402392, No Iterations 1000
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrListFoam::lduMatrix::smoother const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9 Foam::fvMatrix<double>::solve() in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#10 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#11 ? in /lib/x86_64-linux-gnu/libc.so.6
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam
浮点异常 (核心已转储)