[hiep.nguyentrong]

When i try to run rhoSimpleFoam or steadyCompressibleFoam (foam-extend-4.1) with turbulence models (kEpsilon or kOmegaSST), i get this error.
My case is transonic flow over NACA0012. It run fine with laminar or other version like OpenFOAM-v2112 or OpenFOAM-dev.

How can i solve this? Is it a bug?

Selecting RAS turbulence model kEpsilon

Program received signal SIGFPE, Arithmetic exception.
0x00007ffff721d84f in Foam::compressible::RASModels::mutkWallFunctionFvPatchScalarField::calcMut (this=0x555555dd1760) at derivedFvPatchFields/wallFunctions/mutWallFunctions/mutkWallFunction/mutkWallFunctionFvPatchScalarField.C:69
69	            Cmu25*y[faceI]*sqrt(k[faceCellI])/(muw[faceI]/rhow[faceI]);

[mAlletto]

Seems like a division by zero -> rhok is probably zero some where

[hiep.nguyentrong]

Thank you for reply.
I tried to initiate all turbulence variable with difference value but still get floating point at start

[mAlletto]

Usually one adds a small value to the denominator to avoid divisions by zero

[hiep.nguyentrong]

i tried add 0.1 value at wall and internalField, but still get the error. I did same with other field and nothing change

    "(WALL3|WALL4)"
    {
        type            mutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.1;
    }

[mAlletto]

The small value like 1.0e-04 should go in the denominator in the source code

[hiep.nguyentrong]

Thank you but i did not solve the problem. Here is my 0 file:

alphat:

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

internalField   uniform 1e-4;

boundaryField
{
    "(INLE1|PRES2)"
    {
        type            calculated;
        value           uniform 1e-4;
    }

    "(WALL3|WALL4)"
    {
        type            alphatWallFunction;
        value           uniform 1e-4;
    }

}

k:

boundaryField
{
    "(INLE1|PRES2)"
    {
        type            inletOutlet;
        inletValue      uniform $kInlet;
        value           uniform $kInlet;
    }

    "(WALL3|WALL4)"
    {
        type            kqRWallFunction;
        value           uniform $kInlet;
    }

    //#includeEtc "caseDicts/setConstraintTypes"
}

mut:

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

internalField   uniform 1e-4;

boundaryField
{
    "(WALL3|WALL4)"
    {
        type            mutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 1e-4;
    }
    "(INLE1|PRES2)"
    {
        type            calculated;
        value           uniform 1e-4;
    }
    frontAndBackPlanes
    {
        type            empty;
    }
}

omega

//omegaInlet      10;
omegaInlet 		4000;

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

internalField   uniform $omegaInlet;

boundaryField
{
    "(INLE1|PRES2)"
    {
        type            inletOutlet;
        inletValue      uniform $omegaInlet;
        value           uniform $omegaInlet;
    }

    "(WALL3|WALL4)"
    {
        type            omegaWallFunction;
        value           uniform $omegaInlet;
    }

    //#includeEtc "caseDicts/setConstraintTypes"
}

[mAlletto]

The keyword value is just a placeholder since the boundary conditions are derived from a fixedValue b.c. you have to modify the source code to avoid a floating point exception