[jason_huang]

Hi everyone, I am simulating flow around DARPA SUBOFF using steady simpleFOAM k-omega SST model, but the metioned error occur.

boundary condition descreption:
1. wall - darpa suboff
2. in - inlet
3. out - outlet
4. sym - symmetry plane

error:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> cons, 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<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:?
#9  ? at ??:?
#10  __libc_start_main in /lib64/libc.so.6
#11  ? at ??:?
Floating point exception (core dumped)

checkMesh output: OK

Mesh stats
    points:           228080
    faces:            663159
    internal faces:   642561
    cells:            217620
    faces per cell:   6
    boundary patches: 4
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     217620
    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
    in                  961      1024     ok (non-closed singly connected)
    out                 961      1024     ok (non-closed singly connected)
    sym                 13096    13340    ok (non-closed singly connected)
    wall                5580     5702     ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking geometry...
    Overall domain bounding box (-185.03937 -196.85039 -8.9412195e-10) (346.45669 196.85039 196.85039)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (1.4344067e-16 -3.0733581e-16 -5.3925414e-16) OK.
    Max cell openness = 3.4660301e-16 OK.
    Max aspect ratio = 35.858454 OK.
    Minimum face area = 0.0094103014. Maximum face area = 445.01098.  Face area magnitudes OK.
    Min volume = 0.00098242557. Max volume = 4140.3589.  Total volume = 24304384.  Cell volumes OK.
    Mesh non-orthogonality Max: 66.862872 average: 14.839491
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.98146085 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

fvSchemes

ddtSchemes
{
    default steadyState;
}

gradSchemes
{
    default Gauss linear;

}

divSchemes
{
    default Gauss linear;

}

laplacianSchemes
{
    default Gauss linear corrected;
}

interpolationSchemes
{
    default linear;
}

snGradSchemes
{
    default corrected;
}

wallDist
{
    method meshWave;
}

fvSolution:

solvers
{
    p
    {
        solver GAMG;
        smoother GaussSeidel;
        tolerance 1e-6;
        relTol 0.1;
    }
    Phi
    {
        $p;
    }
    "(U|k|omega)"
    {
       solver smoothSolver;
       smoother GaussSeidel;
       tolerance 1e-8;
       relTol 0.1;
       nSweeps 1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    pRefCell 0;
    pRefvalue 0;
    consistent yes;
}

relaxationFactors
{
    fields
    {
        p 0.3;
    }
    equations
    {
        U 0.9;
        ".*" 0.9;
    }
}

controlDict:

application simpleFoam;
startFrom  startTime;
startTime 0;
stopAt endTime;
endTime 200;
deltaT 1;

writeControl timeStep;
writeInterval 20;
purgeWrite 0;
writeFormat ascii;
writePrecision 8;
writeConpression off;
timeFormat fixed;
timePrecision 4;
graphFormat gnuplot;

runTimeModifiable false;

initial U:

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (9.25485556 0 0);

boundaryField
{
    in
    {
        type            freestream;
        freestreamValue $internalField;
    }

    out
    {
         type            freestream;
        freestreamValue  $internalField;
    }

    sym
    {
        type            symmetry;
    }

    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

}

initial p:

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    in
    {
      type freestreamPressure;
      freestreamValue $internalField;

    out
    {
      type freestreamPressure;
      freestreamValue $internalField;
    }

    sym
    {
        type            symmetry;
    }

    wall
    {
        type            zeroGradient;
    }
}

initial nut

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 1.25e-7;

boundaryField
{
    in
    {
        type            freestream;
        freestreamValue           $internalField;
    }

    out
    {
              type            freestream;
        freestreamValue           $internalField;

    }

    wall
    {
        type            freestream;
        freestreamValue           $internalField;
    }

    sym
    {
        type            symmetry;
    }

}

initial k

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform  1.5e-4; // 1.08e-3 = 1.125*(69.4)^2/5e6 // 11=3.9% 18.1=5%;

boundaryField
{
    in
    {
        type            fixedValue;
        value           $internalField;
    }

    out
    {
        type            zeroGradient;
    }

    wall
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    sym
    {
        type            symmetry;
    }

}

initial omega

dimensions      [0 0 -1 0 0 0 0];

internalField   uniform 2.236;

boundaryField
{
    in
    {
        type            fixedValue;
        value           $internalField;
    }

    out
    {
        type            zeroGradient;
    }

    sym
    {
        type            symmetry;
    }

    wall
    {
        type            omegaWallFunction;
        value           $internalField;
    }

}

Thank you all

my case if needed
https://drive.google.com/file/d/1YCN...ew?usp=sharing