CFD Online Discussion Forums - Foam::error::printStack(Foam::Ostream&) with simpleFoam -parallel

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Foam::error::printStack(Foam::Ostream&) with simpleFoam -parallel (https://www.cfd-online.com/Forums/openfoam-solving/85341-foam-error-printstack-foam-ostream-simplefoam-parallel.html)

Foam::error::printStack(Foam::Ostream&)

I solved the problem changing from steadyState to Euler (on fvSchemes).

Now it is:

ddtSchemes
{
default Euler;
}

Foam::error::printStack(Foam::Ostream&)

Hi, I'm getting similar error. I've tried a lot of things but i didn't solve my problem. My solver is interFoam. Just look at deltaT and alpha.water values in error file. How can i fix this? Any suggestion?
(I used SolidWorks 2016, salome 8.4.0, openfoam-7, ubuntu 18.04 and i customized classic laminar/damBreak case.)

boundary

Code:

{

    version     2.0;

    format      ascii;

    class       polyBoundaryMesh;

    location    "constant/polyMesh";

    object      boundary;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



5

(

    leftWall

    {

        type            wall;

        inGroups        List<word> 1(wall);

        nFaces          116;

        startFace       26049;

    }

    rightWall

    {

        type            wall;

        inGroups        List<word> 1(wall);

        nFaces          80;

        startFace       26165;

    }

    lowerWall

    {

        type            wall;

        inGroups        List<word> 1(wall);

        nFaces          136;

        startFace       26245;

    }

    atmosphere

    {

        type            patch;

        nFaces          120;

        startFace       26381;

    }

    defaultFaces

    {

        type            empty;

        inGroups        List<word> 1(empty);

        nFaces          9918;

        startFace       26501;

    }

)



// ************************************************************************* //

alpha.water

Code:

FoamFile

{

    version     2.0;

    format      ascii;

    class       volScalarField;

    object      alpha.water;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



dimensions      [0 0 0 0 0 0 0];



internalField   uniform 0;



boundaryField

{

    leftWall

    {

        type            zeroGradient;

    }



    rightWall

    {

        type            zeroGradient;

    }



    lowerWall

    {

        type            zeroGradient;

    }



    atmosphere

    {

        type            inletOutlet;

        inletValue      uniform 0;

        value           uniform 0;

    }



    defaultFaces

    {

        type            empty;

    }

}

p_rgh

Code:

FoamFile

{

    version     2.0;

    format      ascii;

    class       volScalarField;

    object      p_rgh;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



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



internalField   uniform 0;



boundaryField

{

    leftWall

    {

        type            fixedFluxPressure;

        value           uniform 0;

    }



    rightWall

    {

        type            fixedFluxPressure;

        value           uniform 0;

    }



    lowerWall

    {

        type            fixedFluxPressure;

        value           uniform 0;

    }



    atmosphere

    {

        type            totalPressure;

        p0              uniform 0;

    }



    defaultFaces

    {

        type            empty;

    }

}

Code:

FoamFile

{

    version     2.0;

    format      ascii;

    class       volVectorField;

    location    "0";

    object      U;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



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



internalField   uniform (0 0 0);



boundaryField

{

    leftWall

    {

        type            noSlip;

    }

    rightWall

    {

        type            noSlip;

    }

    lowerWall

    {

        type            noSlip;

    }

    atmosphere

    {

        type            pressureInletOutletVelocity;

        value           uniform (0 0 0);

    }

    defaultFaces

    {

        type            empty;

    }

}

HERE MY ERROR

Code:

smoothSolver:  Solving for alpha.water, Initial residual = 5.98474e-05, Final residual = 3.88048e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.44338e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -3.52175e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.781862, Final residual = 0.0317406, No Iterations 8

time step continuity errors : sum local = 0.000492867, global = -3.45236e-08, cumulative = 0.00744918

DICPCG:  Solving for p_rgh, Initial residual = 0.608476, Final residual = 0.0248548, No Iterations 8

time step continuity errors : sum local = 0.00162432, global = 1.66769e-08, cumulative = 0.0074492

DICPCG:  Solving for p_rgh, Initial residual = 0.618346, Final residual = 9.228e-08, No Iterations 362

time step continuity errors : sum local = 2.749e-08, global = -8.68054e-10, cumulative = 0.0074492

ExecutionTime = 35.66 s  ClockTime = 36 s



Courant Number mean: 0.0396433 max: 17.064

Interface Courant Number mean: 0.00390682 max: 5.60411

deltaT = 7.77403e-100

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 5.84367e-05, Final residual = 3.60982e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.44162e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -5.0219e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.779072, Final residual = 0.0316353, No Iterations 8

time step continuity errors : sum local = 7.53892e-05, global = -4.91121e-09, cumulative = 0.00744919

DICPCG:  Solving for p_rgh, Initial residual = 0.607604, Final residual = 0.0248692, No Iterations 8

time step continuity errors : sum local = 0.000248544, global = 2.35874e-09, cumulative = 0.00744919

DICPCG:  Solving for p_rgh, Initial residual = 0.616997, Final residual = 9.84716e-08, No Iterations 362

time step continuity errors : sum local = 4.46598e-09, global = -1.34046e-10, cumulative = 0.00744919

ExecutionTime = 36.11 s  ClockTime = 37 s



Courant Number mean: 0.00599844 max: 2.59439

Interface Courant Number mean: 0.000588719 max: 0.891569

deltaT = 1.49824e-100

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 5.73287e-05, Final residual = 3.53128e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.43785e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -6.81171e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.776517, Final residual = 0.0313868, No Iterations 8

time step continuity errors : sum local = 0.00048802, global = -2.93771e-08, cumulative = 0.00744916

DICPCG:  Solving for p_rgh, Initial residual = 0.606532, Final residual = 0.0247253, No Iterations 8

time step continuity errors : sum local = 0.00160604, global = 1.3903e-08, cumulative = 0.00744918

DICPCG:  Solving for p_rgh, Initial residual = 0.615519, Final residual = 9.41518e-08, No Iterations 340

time step continuity errors : sum local = 2.75985e-08, global = -1.00355e-09, cumulative = 0.00744918

ExecutionTime = 36.54 s  ClockTime = 37 s



Courant Number mean: 0.0385174 max: 16.7802

Interface Courant Number mean: 0.00380786 max: 6.09569

deltaT = 4.4643e-102

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 5.61634e-05, Final residual = 3.3963e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.42924e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.35845e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.773262, Final residual = 0.0307347, No Iterations 8

time step continuity errors : sum local = 7.27895e-05, global = -3.94879e-09, cumulative = 0.00744917

DICPCG:  Solving for p_rgh, Initial residual = 0.604857, Final residual = 0.024329, No Iterations 8

time step continuity errors : sum local = 0.000239488, global = 2.02513e-09, cumulative = 0.00744917

DICPCG:  Solving for p_rgh, Initial residual = 0.613341, Final residual = 9.64352e-08, No Iterations 362

time step continuity errors : sum local = 4.2589e-09, global = -1.19707e-10, cumulative = 0.00744917

ExecutionTime = 36.98 s  ClockTime = 38 s



Courant Number mean: 0.00578662 max: 2.56257

Interface Courant Number mean: 0.000578441 max: 0.973182

deltaT = 8.71061e-103

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 5.49118e-05, Final residual = 3.2171e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.4112e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -3.96341e-16  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.770017, Final residual = 0.0299674, No Iterations 8

time step continuity errors : sum local = 0.000454079, global = -2.18559e-08, cumulative = 0.00744915

DICPCG:  Solving for p_rgh, Initial residual = 0.602994, Final residual = 0.0293486, No Iterations 5

time step continuity errors : sum local = 0.00183973, global = 9.67059e-08, cumulative = 0.00744925

DICPCG:  Solving for p_rgh, Initial residual = 0.620295, Final residual = 9.23708e-08, No Iterations 342

time step continuity errors : sum local = 2.46046e-08, global = -9.39872e-11, cumulative = 0.00744925

ExecutionTime = 37.41 s  ClockTime = 38 s



Courant Number mean: 0.0322515 max: 17.5305

Interface Courant Number mean: 0.00359975 max: 6.76509

deltaT = 2.48441e-104

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 4.83567e-05, Final residual = 1.91459e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -1.04402e-46  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -2.51392e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.742106, Final residual = 0.0363213, No Iterations 7

time step continuity errors : sum local = 6.49571e-05, global = -5.59286e-11, cumulative = 0.00744925

DICPCG:  Solving for p_rgh, Initial residual = 0.606875, Final residual = 0.0274085, No Iterations 5

time step continuity errors : sum local = 0.000207145, global = 4.45286e-09, cumulative = 0.00744925

DICPCG:  Solving for p_rgh, Initial residual = 0.619751, Final residual = 9.17399e-08, No Iterations 362

time step continuity errors : sum local = 3.04518e-09, global = -6.90435e-11, cumulative = 0.00744925

ExecutionTime = 37.86 s  ClockTime = 38 s



Courant Number mean: 0.0040732 max: 2.24028

Interface Courant Number mean: 0.000418459 max: 0.609483

deltaT = 5.54488e-105

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 4.76685e-05, Final residual = 2.01965e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -9.91398e-47  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -3.04471e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.744592, Final residual = 0.0363254, No Iterations 5

time step continuity errors : sum local = 0.000498639, global = 1.99994e-08, cumulative = 0.00744927

DICPCG:  Solving for p_rgh, Initial residual = 0.623883, Final residual = 0.0268522, No Iterations 5

time step continuity errors : sum local = 0.00154667, global = -3.01824e-08, cumulative = 0.00744924

DICPCG:  Solving for p_rgh, Initial residual = 0.613925, Final residual = 9.08759e-08, No Iterations 363

time step continuity errors : sum local = 2.22317e-08, global = -4.95163e-10, cumulative = 0.00744924

ExecutionTime = 38.34 s  ClockTime = 39 s



Courant Number mean: 0.0304221 max: 16.379

Interface Courant Number mean: 0.00306763 max: 4.27059

deltaT = 1.69268e-106

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 4.75175e-05, Final residual = 1.99622e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.24586e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -3.42133e-17  Max(alpha.water) = 1

DICPCG:  Solving for p_rgh, Initial residual = 0.743859, Final residual = 0.0325512, No Iterations 5

time step continuity errors : sum local = 6.12945e-05, global = 4.61408e-09, cumulative = 0.00744925

DICPCG:  Solving for p_rgh, Initial residual = 0.616162, Final residual = 0.0244443, No Iterations 5

time step continuity errors : sum local = 0.000188544, global = -8.24671e-09, cumulative = 0.00744924

DICPCG:  Solving for p_rgh, Initial residual = 0.605392, Final residual = 9.66245e-08, No Iterations 362

time step continuity errors : sum local = 3.08855e-09, global = -3.55315e-11, cumulative = 0.00744924

ExecutionTime = 38.79 s  ClockTime = 39 s



Courant Number mean: 0.0042222 max: 2.28614

Interface Courant Number mean: 0.000412843 max: 0.511151

deltaT = 3.70205e-107

Time = 0.00119069



smoothSolver:  Solving for alpha.water, Initial residual = 4.32103e-05, Final residual = 1.9414e-11, No Iterations 2

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -4.1624e-68  Max(alpha.water) = 1

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.0746171  Min(alpha.water) = -2.01258e-17  Max(alpha.water) = 1

#0  Foam::error::printStack(Foam::Ostream&) at ??:?

#1  Foam::sigFpe::sigHandler(int) at ??:?

#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"

#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

#4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?

#5  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?

#6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/interFoam"

#7  Foam::fvMatrix<double>::solve() in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/interFoam"

#8  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/interFoam"

#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#10  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/interFoam"

Hi Said,

Such small values of deltaT are an indication that the velocity is extremely large. Perhaps it would be a good idea to monitor your minimum and maximum values of velocity in the domain.

That said, I'm also having a similar error but I am using compressibleInterFoam, I was wondering IF and HOW you have managed to solve this error?

I look forward to your reply!

Regards,

JM

Quote:

Originally Posted by mbay101 (Post 443898)

Hi everybody,

I m having as it seems the same Probleme :(. I m using chtMultiRegionSimpleFoam in OP 2.2.0 to simulate a solid Region and a Fluid Region Tempratur Transfer. The Boundary Condition were taking from the OF example HeatTransfert and i did the same with the Files in system. I m using Salome to create my mesh after that i use mappedWall in Boundary to connect the Region.

Code:

#0 Foam::error::printStack(Foam::Ostream&) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#1 Foam::sigSegv::sigHandler(int) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#2 

at sigaction.c:0

#3 Foam::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#4 Foam::polyMesh::tetBasePtIs() const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#5 Foam::mappedPatchBase::facePoints(Foam::polyPatch const&) const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"

#6 Foam::mappedPatchBase::calcMapping() const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"

#7 void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"

#8 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"

#9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"

#10 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"

#11 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"

#12 

in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"

#13 

in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"

#14 __libc_start_main in "/lib64/libc.so.6"

#15 

in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"

./solve: line 5: 10011 Segmentation fault (core dumped) chtMultiRegionSimpleFoam

I attached to this message my system and my 0 directory.
Please!! Please!! take a look at them and if you have any idea what so ever that can help me understand the reason then please shared with me.

In working a Project and i don t have any support from my Team. This is my only option for me to get same Tipps and Help.

Thank you!!
best Regards

Did you have a solution? I have the same question, maybe you can help me somewhere.

Hello everyone!!
I am trying to run a lagrangian simulation using sprayFoam solver and I am getting an error as shown

#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::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5 ? in "/home/riddhideep/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/sprayFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/home/riddhideep/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/sprayFoam"
Floating point exception (core dumped)

Below I am attaching all the files I am using in the form of a google drive link.
Please take a look at the files and help me how to get rid of this error.

https://drive.google.com/drive/folde...-3?usp=sharing

Problem in simulating the liquid jet

Hi,

I made a closed geometry in Solidwords and imported it into openfoam (trisurface).I defined the blockmesh such that it encloses the structure.I am using surfacefeatureextract and snappyhexmesh.Runnning the case just for streamline case now using simplefoam. The error is when i convert to FoamToVTK than it says

size 50000 is not in the range of 1872638.
Can anyone suggest?

Thanks
Monika

Reply to problem simulation liquid jet

Quote:

Originally Posted by Monika Yadav (Post 783803)

Whenever I recieve a message like that it means that I made a change to my mesh after having set initial values. You should check all the files in the 0 directory and make sure that none of them contain a list of cell points with cell points from the old mesh.
My guess is that you ran a simulation, refined the mesh and then not updated your 0 folder.

Hope this helps!

Hello dear foamers,

I have a similar problem with simulating in parallel using rhoCentralFoam. I have tried many changes in my setup, but have not been successful. Here is my error notes and checkMesh notes. I also shared my folders.

I'm looking forward to your suggestions, thanks in advance..

Quote:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1812 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : v1812 OPENFOAM=1812
Arch : "LSB;label=32;scalar=64"
Exec : rhoCentralFoam -parallel -case /okyanus/users/afadil/acf/00B
Date : Mar 26 2022
Time : 12:47:32
Host : s205
PID : 268505
I/O : uncollated
Case : /okyanus/users/afadil/acf/00B
nProcs : 40
Hosts :
(
(s205 40)
)
Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Using LTS
Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
decayControl false;
kInf 0;
omegaInf 0;
}

Creating finite volume options from "system/fvOptions"

Selecting finite volume options type limitTemperature
Source: limitT
- selecting all cells
- selected 1646904 cell(s) with volume 6.10667e-05

Starting time loop

fieldAverage fieldAverage1:
Restarting averaging for fields:
U: starting averaging at time 0
p: starting averaging at time 0

Mean and max Courant Numbers = 0.0014347 0.262613
Flow time scale min/max = 4.28387e-09, 2.03396e-08
Time = 1.125e-09

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 9.91481e-08, No Iterations 159
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 9.53803e-08, No Iterations 131
smoothSolver: Solving for Uz, Initial residual = 0.481511, Final residual = 9.73196e-08, No Iterations 111
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for e, Initial residual = 0.095494, Final residual = 0.000522242, No Iterations 2
GAMG: Solving for omega, Initial residual = 0.017991, Final residual = 8.45975e-05, No Iterations 3
bounding omega, min: -79757.3 max: 3.89408e+09 average: 2.52554e+08
smoothSolver: Solving for k, Initial residual = 1, Final residual = 4.9169e-08, No Iterations 12
ExecutionTime = 2.07 s ClockTime = 3 s

fieldAverage fieldAverage1:
Reading/initialising field UMean
Reading/initialising field pMean
Reading/initialising field UPrime2Mean
Reading/initialising field pPrime2Mean

fieldAverage fieldAverage1 write:
Calculating averages

Mean and max Courant Numbers = 0.00143492 0.261351
Flow time scale min/max = 4.30455e-09, 2.03404e-08
Time = 2.25e-09

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.000994925, Final residual = 2.68726e-08, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.000364341, Final residual = 7.56062e-08, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.0042112, Final residual = 2.08737e-09, No Iterations 2
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for e, Initial residual = 0.000137899, Final residual = 1.74764e-09, No Iterations 1
--> FOAM Warning :
From function virtual void Foam::fv::option::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 125
Source limitT defined for field e but never used
--> FOAM Warning :
From function virtual void Foam::fv::option::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 125
Source limitT defined for field e but never used
GAMG: Solving for omega, Initial residual = 0.000546902, Final residual = 3.751e-07, No Iterations 1
--> FOAM Warning :
From function virtual void Foam::fv::option::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 125
Source limitT defined for field e but never used
--> FOAM Warning :
From function virtual void Foam::fv::option::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 125
Source limitT defined for field e but never used
smoothSolver: Solving for k, Initial residual = 0.330261, Final residual = 3.36906e-08, No Iterations 5
ExecutionTime = 2.54 s ClockTime = 4 s

fieldAverage fieldAverage1 write:
Calculating averages

Mean and max Courant Numbers = 0.00143609 0.365354
Flow time scale min/max = 3.07921e-09, 2.12312e-08
Time = 3.375e-09

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.00117023, Final residual = 4.19306e-08, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.000554007, Final residual = 9.40954e-08, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.00465589, Final residual = 1.67531e-09, No Iterations 2
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for e, Initial residual = 0.000217566, Final residual = 1.90397e-09, No Iterations 1
--------------------------------------------------------------------------
A process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

Local host: [[45705,1],24] (PID 268529)

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
[s205:268529] *** Process received signal ***
[s205:268529] Signal: Floating point exception (8)
[s205:268529] Signal code: (-6)
[s205:268529] Failing at address: 0xb2c000418f1
[s205:268529] [ 0] /lib64/libc.so.6(+0x362f0)[0x2b7f73d292f0]
[s205:268529] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x2b7f73d29277]
[s205:268529] [ 2] /lib64/libc.so.6(+0x362f0)[0x2b7f73d292f0]
[s205:268529] [ 3] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_ EES5_+0x25)[0x2b7f72da8f85]
[s205:268529] [ 4] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam9PBiCGStab5solveERNS_5Fiel dIdEERKS2_h+0xb4b)[0x2b7f72bfedcb]
[s205:268529] [ 5] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarses tLevelERNS_5FieldIdEERKS2_+0x1ee)[0x2b7f72c193be]
[s205:268529] [ 6] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x71b)[0x2b7f72c1ac3b]
[s205:268529] [ 7] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x523)[0x2b7f72c1ce33]
[s205:268529] [ 8] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x120)[0x2b7f6eec1790]
[s205:268529] [ 9] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegr egatedOrCoupledERKNS_10dictionaryE+0x27f)[0x2b7f6e9bf87f]
[s205:268529] [10] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvM atrixIdEERKNS_10dictionaryE+0xf)[0x2b7f6e98719f]
[s205:268529] [11] rhoCentralFoam(_ZN4Foam5solveIdEENS_17SolverPerfor manceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x45)[0x43ebd5]
[s205:268529] [12] /okyanus/progs/OpenFOAM/OpenFOAM.com-v1812/OpenFOAM-v1812/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam13kOmeg aSSTBaseINS_13eddyViscosityINS_8RASModelINS_15Eddy DiffusivityINS_18ThermalDiffusivityINS_27Compressi bleTurbulenceModelINS_11fluidThermoEEEEEEEEEEEE7co rrectEv+0xe10)[0x2b7f71be5890]
[s205:268529] [13] rhoCentralFoam[0x4297dd]
[s205:268529] [14] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b7f73d15445]
[s205:268529] [15] rhoCentralFoam[0x42ccbe]
[s205:268529] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 24 with PID 268529 on node s205 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Quote:

[24] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[24] #1 Foam::sigFpe::sigHandler(int) at ??:?
[24] #2 ? in /lib64/libc.so.6
[24] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[24] #4 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[24] #5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
[24] #6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 ??:?
[24] #7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[24] #8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[24] #9 Foam::fvMatrix<double>::solveSegregatedOrCoupled(F oam::dictionary const&) at ??:?
[24] #10 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[24] #11 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[24] #12 Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASM odel<Foam::EddyDiffusivity<Foam::ThermalDiffusivit y<Foam::CompressibleTurbulenceModel<Foam::fluidThe rmo> > > > > >::correct() at ??:?
[24] #13 ? at ??:?
[24] #14 __libc_start_main in /lib64/libc.so.6
[24] #15 ? at ??:?
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1812 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : v1812 OPENFOAM=1812
Arch : "LSB;label=32;scalar=64"
Exec : reconstructPar -latestTime
Date : Mar 26 2022
Time : 12:47:39
Host : s205
PID : 269361
I/O : uncollated
Case : /okyanus/users/afadil/acf/00B
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM Warning :
From function int main(int, char**)
in file reconstructPar.C at line 242
No times selected

Quote:

Originally Posted by AJAY BHANDARI (Post 601611)

Thanks ehsan for your quick reply. But i was able to solve it today.

Actually in my code the variable orl value was wrongly input by me as 1 which made the denominator term zero. thats why sigFpe error was coming.

But i corrected it and solver ran successfully. But thanks for your reply...

hello Ajay, having the same problem, in serial. this is the error message i get
How did you fix?

Quote:

#0 Foam::error:: printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solveSegregatedOrCoupled(F oam::dictionary const&) at ??:?
#10 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#11 ? in /usr/lib/openfoam/openfoam2106/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in /usr/lib/openfoam/openfoam2106/platforms/linux64GccDPInt32Opt/bin/simpleFoam

Hi gohome,

Thanks for your suggestion. I am solving compressible flow with Peng Robinson equation of state and Sutherland for viscosity. I got a similar error where initial residual for h was close to 0.99 and floating point exception. I used 0.4 relaxation factor for p, rho, h, U, k, e and it worked!
p { line-height: 115%; margin-bottom: 0.25cm; background: transparent }

hi can anyone help me with this error

smoothSolver: Solving for alpha.water, Initial residual = 1.14414e-05, Final residual = 8.06717e-09, No Iterations 1
Phase-1 volume fraction = 0.373702 Min(alpha.water) = -0.0191484 Max(alpha.water) = 1.38326
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.373702 Min(alpha.water) = -0.0191484 Max(alpha.water) = 1.35623
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0158173, No Iterations 50
time step continuity errors : sum local = 4.82275e-07, global = 3.02594e-07, cumulative = 3.8436e+34
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
#10 Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
#11 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
#12 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/interFoam"
Floating point exception (core dumped)

error in openfoam

I am trying to simulate rigid body motion of two cylinder in viv using overset mesh. but i found this error during simulation. Please suggest me some hint to how to handle multi cylinder using overset mesh in viv case.

Thank you in advance...
Mesh Courant Number mean: 0 max: 0
PIMPLE: iteration 1
#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::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:?
#7 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#9 ? at ??:?
#10 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<do uble> >&, Foam::dictionary const&) const at ??:?
#11 ? at ??:?
#12 ? at ??:?
#13 ? in /usr/lib/openfoam/openfoam2206/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam
#14 ? in /lib/x86_64-linux-gnu/libc.so.6
#15 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#16 ? in /usr/lib/openfoam/openfoam2206/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam
Floating point exception (core dumped)

Starting time loop

Courant Number mean: 2.77415e-06 max: 0.00125
deltaT = 0.00120048
Time = 0.00120048

forces forces:
rho: rhoInf
Freestream density (rhoInf) set to 1
Not including porosity effects

forces forces:
rho: rhoInf
Freestream density (rhoInf) set to 1
Not including porosity effects

Rigid-body motion of the cylinder1
Centre of rotation: (16 0 0)
Orientation: (1 0 0 0 1 0 0 0 1)
Linear velocity: (0 0 0)
Angular velocity: (0 0 0)
Rigid-body motion of the cylinder2
Centre of rotation: (5.19615 3 0)
Orientation: (1 0 0 0 1 0 0 0 1)
Linear velocity: (0 0 0)
Angular velocity: (0 0 0)
inverseDistance : detected 2 mesh regions
zone:0 nCells:90000 voxels:(304 304 1) bb:(-30.0002 -40.0002 -0.500197) (150 40.0002 0.500197)
zone:1 nCells:3000 voxels:(304 304 1) bb:(-17.9982 -3.49957 -0.500006) (-14 0.499567 0.500006)
Overset analysis : nCells : 93000
calculated : 92798
interpolated : 170 (interpolated from local:95 mixed local/remote:0 remote:0)
hole : 32

Mesh Courant Number mean: 0 max: 0
PIMPLE: iteration 1
#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::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) at ??:?