[nayem.cosmic]

Hello everyone
I am trying to calculate resistance of wigley hull with interFoam solver in OpenFOAM 9. Used following base case.
tutorials/multiphase/interFoam/RAS/DTCHull
Just changed the geometry stl file and changed relevant values in some files.

Run the solver in single core. Solver is giving the following error after running 80-90 time steps.

#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::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#4  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::kinematicTransportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::kinematicTransportModel> >::correct() at ??:?
#5  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
Floating point exception (core dumped)

Have checked the mesh. It prints OK.

Mesh stats
    points:           1826029
    faces:            5286320
    internal faces:   5205924
    cells:            1731246
    faces per cell:   6.06052
    boundary patches: 7
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     1638870
    prisms:        7727
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     84649
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            4   1403
            5   1230
            6   1458
            7   62037
            8   1069
            9   17452

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
              atmosphere      798      860  ok (non-closed singly connected)
                   inlet     3382     3580  ok (non-closed singly connected)
                  outlet     3382     3580  ok (non-closed singly connected)
                  bottom      798      860  ok (non-closed singly connected)
                    side     7476     7697  ok (non-closed singly connected)
                midPlane    29595    30776  ok (non-closed singly connected)
                    hull    34965    36843  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-26 -19 -16) (16 3.51072e-12 4)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (3.66952e-16 1.95783e-15 -1.02697e-16) OK.
    Max cell openness = 4.12411e-16 OK.
    Max aspect ratio = 71.827 OK.
    Minimum face area = 1.26253e-05. Maximum face area = 1.00711.  Face area magnitudes OK.
    Min volume = 7.66276e-08. Max volume = 0.934091.  Total volume = 15957.3.  Cell volumes OK.
    Mesh non-orthogonality Max: 69.9801 average: 5.7433
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.99295 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

fvScheme file:

ddtSchemes
{
    default         localEuler;
}

gradSchemes
{
    default         Gauss linear;
    limitedGrad     cellLimited Gauss linear 1;
}

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss interfaceCompression vanLeer 1;
    div(phi,k)      Gauss linearUpwind limitedGrad;
    div(phi,omega)  Gauss linearUpwind limitedGrad;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

fvSolution file:

solvers
{
    "alpha.water.*"
    {
        nAlphaCorr      2;
        nAlphaSubCycles 1;

        MULESCorr       yes;
        nLimiterIter    10;
        alphaApplyPrevCorr  yes;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          0;
        minIter         1;
    }

    "pcorr.*"
    {
        solver          PCG;

        preconditioner
        {
            preconditioner  GAMG;

            smoother        GaussSeidel;

            tolerance       1e-5;
            relTol          0;
        };

        tolerance       1e-5;
        relTol          0;
    };

    p_rgh
    {
        solver          GAMG;

        smoother        DIC;

        tolerance       1e-7;
        relTol          0.01;
    };

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "(U|k|omega).*"
    {
        solver          smoothSolver;

        smoother        symGaussSeidel;
        nSweeps         1;

        tolerance       1e-7;
        relTol          0.1;
        minIter         1;
    };
}

PIMPLE
{
    momentumPredictor   no;

    nOuterCorrectors    1;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;

    maxCo               10;
    maxAlphaCo          5;

    rDeltaTSmoothingCoeff 0.05;
    rDeltaTDampingCoeff 0.5;
    nAlphaSpreadIter    0;
    nAlphaSweepIter     0;
    maxDeltaT           1;
}

relaxationFactors
{
    equations
    {
        ".*" 1;
    }
}

cache
{
    grad(U);
}

Please find interFoam log file in attachment.
Also find the case directory here.

I am new to OpenFOAM. Can you help me?
Thanks.

[nayem.cosmic]

Update:
Solved my issue.
Changed relaxation factor to below 1 in fvSolution file.

relaxationFactors
{
    equations
    {
        ".*" 0.7;
    }
}

[chandrana70]

Originally Posted by nayem.cosmic

Update:
Solved my issue.
Changed relaxation factor to below 1 in fvSolution file.

relaxationFactors
{
    equations
    {
        ".*" 0.7;
    }
}

I'm just wondering, why should be it at 0.7? I've tried my submarine simulation, and the last solver, it gives the same error. could you tell me?

[nayem.cosmic]

There is no fixed value as far as I know. You might lower the value and try.

[chandrana70]

Originally Posted by nayem.cosmic

There is no fixed value as far as I know. You might lower the value and try.

Ahh I see. Thanks for your enlightenment