Why RNGkepsilon model gives floating error

shipman · August 31, 2013, 21:19

Dear Foam users;

I am studying cavitating flow inside of the nozzle. As turbulence model i have tried kwsst and kepsilon models. They worked fine, there was no problem. However, when i swtiched to RNGkepslion model, after 2 iterations suddenly velocity reaches crazy value and giving following folating point error;

HTML Code:

in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #12  
[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[0] #10  [3]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[3] #10  [2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[2] #10  
[4]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[4] #9  
[1]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[1] #10  

[7]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[7] #11  __libc_start_main[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[0] #11  __libc_start_main[5]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[5] #10  [6]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[efs-desktop:29586] *** Process received signal ***
[efs-desktop:29586] Signal: Floating point exception (8)
[efs-desktop:29586] Signal code:  (-6)
[efs-desktop:29586] Failing at address: 0x3e800007392
[efs-desktop:29586] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364c0) [0x7fcf5d7894c0]
[efs-desktop:29586] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fcf5d789445]
[efs-desktop:29586] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364c0) [0x7fcf5d7894c0]
[efs-desktop:29586] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditioner15calcReciprocalDERNS_5FieldIdEERKNS_9lduMatrixE+0x97) [0x7fcf5e7bef37]
[efs-desktop:29586] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditionerC2ERKNS_9lduMatrix6solverERKNS_10dictionaryE+0x8e) [0x7fcf5e7bf0be]
[efs-desktop:29586] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam9lduMatrix14preconditioner31addasymMatrixConstructorToTableINS_18DILUPreconditionerEE3NewERKNS0_6solverERKNS_10dictionaryE+0x3c) [0x7fcf5e7bf2fc]
[efs-desktop:29586] [ 6] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam9lduMatrix14preconditioner3NewERKNS0_6solverERKNS_10dictionaryE+0x2ad) [0x7fcf5e7b348d]
[efs-desktop:29586] [ 7] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0x6ba) [0x7fcf5e7b777a]
[efs-desktop:29586] [ 8] cavitatingFoam() [0x44c4ca]
[efs-desktop:29586] [ 9] cavitatingFoam() [0x44cef3]
[efs-desktop:29586] [10] cavitatingFoam() [0x41fec6]
[efs-desktop:29586] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fcf5d77476d]
[efs-desktop:29586] [12] cavitatingFoam() [0x42724d]
[efs-desktop:29586] *** End of error message ***

[3]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[3] #11  __libc_start_main[2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[2] #11  __libc_start_main
 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #12   in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #12   in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #12  
[1]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[1] #11  __libc_start_main--------------------------------------------------------------------------
mpirun noticed that process rank 6 with PID 29586 on node efs-desktop exited on signal 8 (Floating point exception).

to get rid of this trouble i tried following things;

I checked the mesh quality which was ok.
set the initial condition of the k and epsilon properl
2. I decreased courant number to decrease time step.
3. Additionaly decreased time step too.
4. Also tried decreasing of the relaxation factors gradually from 1 to 0.1

But still getting same error. Is there anybody who has idea how to correct and run the RNGkepsilon model?

Thanks in advance.

PS: you can see my initial conditions below.

Epsilon

Code:

internalField   uniform 0.672;
boundaryField
{
    inlet
    {
        type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.001; // 0.1cm - half width of nozzle
        k               k;
        value           uniform 0.672;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    wall
    {
        type            epsilonWallFunction;
        value           uniform 0.672;
    }
}

k

Code:

internalField   uniform 0.0361;

boundaryField
{
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;    // 5% turbulent intensity
        value           uniform 0.0361;
    }
   outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    value           $internalField;
       }
    wall
    {
        type            kqRWallFunction;
        value           uniform 0.0361;
    }

divSchemes

Code:

{
    default         Gauss limitedLinear 1;
    div(phiv,rho)   Gauss limitedLinear 1;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phiv,epsilon) Gauss limitedLinear 1;
    div(phiv,k)     Gauss limitedLinear 1;
}

fvsolution:

Code:

    "(U|k|epsilon)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-06;
        relTol          0.1;
    }

    "(U|k|epsilon)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0;
    }
}

PIMPLE
{
    nCorrectors     3;
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
k 0.1;
epsilon 0.1;
}

wyldckat · September 1, 2013, 09:34

Hi George,

There is some information missing from the crash, because there is no indication of which signal was triggered.

Either way, it does give an indication that something went wrong in the preconditioner, which suggests that either:

The equation/matrix solver itself should be different or reconfigured in "fvSolution".
Or that the k-epsilon or some other field went completely in the wrong direction, for some strange reason.

My suggestion is to try and see the results from the 1st and 2nd iterations, so that you can see the values for the k-epsilon fields and figure out if they have values that are very large or very small.

Best regards,
Bruno

shipman · September 4, 2013, 22:36

Dear Bruno,

First thanks so much for your kind reply. According to your answer, one more time i decreased the courant number and relaxation factors for the k and epsilon. This time i thought that it was working but after 200 iteration again suddenly velocity reached the crazy value and gave floating point errror as follows:

could you tell me your idea how can i run my case with RNGkepsilon model?

HTML Code:

DILUPBiCG:  Solving for rho, Initial residual = 0.000100648, Final residual = 3.11001e-09, No Iterations 1
max-min rho: 146486 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.00280351, Final residual = 1.96098e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.00271059, Final residual = 2.8279e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00281404, Final residual = 1.90802e-11, No Iterations 2
max(U) 3.1788e+10
GAMG:  Solving for p, Initial residual = 5.31511e-07, Final residual = 1.99477e-17, No Iterations 1
Predicted p max-min : 3.20457e+11 -4533.95
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20457e+11 2339
max(U) 3.1788e+10
GAMG:  Solving for p, Initial residual = 9.4014e-10, Final residual = 1.91786e-17, No Iterations 1
Predicted p max-min : 3.20457e+11 2115.91
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20457e+11 2339
max(U) 3.1788e+10
GAMG:  Solving for p, Initial residual = 4.06463e-11, Final residual = 4.06463e-11, No Iterations 0
Predicted p max-min : 3.20457e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20457e+11 2339
max(U) 3.1788e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00494363, Final residual = 4.14931e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00560901, Final residual = 1.95116e-10, No Iterations 2
ExecutionTime = 317905 s  ClockTime = 332312 s

phiv Courant Number mean: 0.000162791 max: 0.125743 acoustic max: 1.66985e-07
deltaT = 3.22642e-16
Time = 0.0021332566986488

DILUPBiCG:  Solving for rho, Initial residual = 0.000100573, Final residual = 3.31574e-09, No Iterations 1
max-min rho: 146534 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.00280365, Final residual = 1.95516e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0027107, Final residual = 2.82825e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00281413, Final residual = 1.90943e-11, No Iterations 2
max(U) 3.19753e+10
GAMG:  Solving for p, Initial residual = 5.31558e-07, Final residual = 2.11003e-17, No Iterations 1
Predicted p max-min : 3.20565e+11 -4502.88
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20565e+11 2339
max(U) 3.19753e+10
GAMG:  Solving for p, Initial residual = 9.32718e-10, Final residual = 2.10365e-17, No Iterations 1
Predicted p max-min : 3.20565e+11 2116.61
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20565e+11 2339
max(U) 3.19753e+10
GAMG:  Solving for p, Initial residual = 4.02187e-11, Final residual = 4.02187e-11, No Iterations 0
Predicted p max-min : 3.20565e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20565e+11 2339
max(U) 3.19753e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00494494, Final residual = 4.14844e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00561015, Final residual = 1.95128e-10, No Iterations 2
ExecutionTime = 317909 s  ClockTime = 332315 s

phiv Courant Number mean: 0.000162769 max: 0.125743 acoustic max: 1.65999e-07
deltaT = 3.20736e-16
Time = 0.0021332566986491

DILUPBiCG:  Solving for rho, Initial residual = 0.000100497, Final residual = 3.55485e-09, No Iterations 1
max-min rho: 146583 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.00280379, Final residual = 1.94943e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0027108, Final residual = 2.8286e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00281422, Final residual = 1.91088e-11, No Iterations 2
max(U) 3.21638e+10
GAMG:  Solving for p, Initial residual = 5.31634e-07, Final residual = 2.03497e-17, No Iterations 1
Predicted p max-min : 3.20672e+11 -4472.83
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20672e+11 2339
max(U) 3.21638e+10
GAMG:  Solving for p, Initial residual = 9.26368e-10, Final residual = 1.88368e-17, No Iterations 1
Predicted p max-min : 3.20672e+11 2117.25
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20672e+11 2339
max(U) 3.21638e+10
GAMG:  Solving for p, Initial residual = 3.98619e-11, Final residual = 3.98619e-11, No Iterations 0
Predicted p max-min : 3.20672e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20672e+11 2339
max(U) 3.21638e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00494624, Final residual = 4.14757e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00561128, Final residual = 1.95141e-10, No Iterations 2
bounding k, min: -7.46097 max: 8.79325e+18 average: 5.77234e+14
ExecutionTime = 317912 s  ClockTime = 332319 s

phiv Courant Number mean: 0.000162747 max: 0.125743 acoustic max: 1.65018e-07
deltaT = 3.18841e-16
Time = 0.0021332566986494

DILUPBiCG:  Solving for rho, Initial residual = 0.000100422, Final residual = 3.83622e-09, No Iterations 1
max-min rho: 146632 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.0028292, Final residual = 5.83431e-10, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.00271901, Final residual = 2.7379e-10, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.00281633, Final residual = 8.69345e-09, No Iterations 4
max(U) 3.23533e+10
GAMG:  Solving for p, Initial residual = 5.33665e-07, Final residual = 1.90162e-17, No Iterations 1
Predicted p max-min : 3.20779e+11 -4443.82
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20779e+11 2339
max(U) 3.23533e+10
GAMG:  Solving for p, Initial residual = 9.22181e-10, Final residual = 1.7469e-17, No Iterations 1
Predicted p max-min : 3.20779e+11 2117.82
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20779e+11 2339
max(U) 3.23533e+10
GAMG:  Solving for p, Initial residual = 3.95845e-11, Final residual = 3.95845e-11, No Iterations 0
Predicted p max-min : 3.20779e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20779e+11 2339
max(U) 3.23533e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00493704, Final residual = 4.13782e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00561454, Final residual = 8.37245e-10, No Iterations 4
bounding k, min: -3.66179e+12 max: 8.89808e+18 average: 5.83999e+14
ExecutionTime = 317916 s  ClockTime = 332323 s

phiv Courant Number mean: 0.000162725 max: 0.125743 acoustic max: 1.64043e-07
deltaT = 3.16958e-16
Time = 0.0021332566986498

DILUPBiCG:  Solving for rho, Initial residual = 0.000100348, Final residual = 4.17185e-09, No Iterations 1
max-min rho: 146680 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.0189912, Final residual = 2.96678e-12, No Iterations 7
DILUPBiCG:  Solving for Uy, Initial residual = 0.00790643, Final residual = 3.04621e-12, No Iterations 7
DILUPBiCG:  Solving for Uz, Initial residual = 0.0064106, Final residual = 2.715e-12, No Iterations 7
max(U) 3.2544e+10
GAMG:  Solving for p, Initial residual = 5.32543e-07, Final residual = 2.03893e-17, No Iterations 1
Predicted p max-min : 3.20886e+11 -4415.93
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20886e+11 2339
max(U) 3.2544e+10
GAMG:  Solving for p, Initial residual = 9.18229e-10, Final residual = 1.91541e-17, No Iterations 1
Predicted p max-min : 3.20886e+11 2118.37
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20886e+11 2339
max(U) 3.2544e+10
GAMG:  Solving for p, Initial residual = 3.93125e-11, Final residual = 3.93125e-11, No Iterations 0
Predicted p max-min : 3.20886e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20886e+11 2339
max(U) 3.2544e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00276504, Final residual = 2.31567e-11, No Iterations 2
bounding epsilon, min: -7.3181e+11 max: 1.0235e+32 average: 4.2102e+27
DILUPBiCG:  Solving for k, Initial residual = 0.00625093, Final residual = 6.81452e-11, No Iterations 7
bounding k, min: -2.6928e+12 max: 9.00417e+18 average: 5.90844e+14
ExecutionTime = 317920 s  ClockTime = 332327 s

phiv Courant Number mean: 0.000162704 max: 0.125743 acoustic max: 1.63074e-07
deltaT = 3.15085e-16
Time = 0.0021332566986501

DILUPBiCG:  Solving for rho, Initial residual = 0.000100274, Final residual = 4.57919e-09, No Iterations 1
max-min rho: 146729 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.00282875, Final residual = 5.21283e-10, No Iterations 7
DILUPBiCG:  Solving for Uy, Initial residual = 0.00273531, Final residual = 2.10038e-09, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.0028631, Final residual = 4.17279e-10, No Iterations 7
max(U) 3.27359e+10
GAMG:  Solving for p, Initial residual = 5.34221e-07, Final residual = 2.29906e-17, No Iterations 1
Predicted p max-min : 3.20993e+11 -4389.24
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20993e+11 2339
max(U) 3.27359e+10
GAMG:  Solving for p, Initial residual = 9.14305e-10, Final residual = 2.16781e-17, No Iterations 1
Predicted p max-min : 3.20993e+11 2118.84
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20993e+11 2339
max(U) 3.27359e+10
GAMG:  Solving for p, Initial residual = 3.90572e-11, Final residual = 3.90572e-11, No Iterations 0
Predicted p max-min : 3.20993e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.20993e+11 2339
max(U) 3.27359e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00494136, Final residual = 4.13749e-11, No Iterations 2
bounding epsilon, min: -8.26292e+24 max: 1.04181e+32 average: 4.28475e+27
DILUPBiCG:  Solving for k, Initial residual = 0.0056153, Final residual = 3.88108e-10, No Iterations 6
bounding k, min: -9.64635e+11 max: 9.11153e+18 average: 5.9777e+14
ExecutionTime = 317924 s  ClockTime = 332331 s

phiv Courant Number mean: 0.000162682 max: 0.125743 acoustic max: 1.62111e-07
deltaT = 3.13223e-16
Time = 0.0021332566986504

DILUPBiCG:  Solving for rho, Initial residual = 0.0001002, Final residual = 5.08403e-09, No Iterations 1
max-min rho: 146777 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.386478, Final residual = 1.81158e-09, No Iterations 18
DILUPBiCG:  Solving for Uy, Initial residual = 0.134064, Final residual = 2.56516e-10, No Iterations 19
DILUPBiCG:  Solving for Uz, Initial residual = 0.197069, Final residual = 4.74468e-09, No Iterations 18
max(U) 3.29289e+10
GAMG:  Solving for p, Initial residual = 5.63346e-07, Final residual = 1.99059e-17, No Iterations 1
Predicted p max-min : 3.211e+11 -4363.93
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.211e+11 2339
max(U) 3.29289e+10
GAMG:  Solving for p, Initial residual = 9.10798e-10, Final residual = 1.63385e-17, No Iterations 1
Predicted p max-min : 3.211e+11 2119.27
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.211e+11 2339
max(U) 3.29289e+10
GAMG:  Solving for p, Initial residual = 3.87974e-11, Final residual = 3.87974e-11, No Iterations 0
Predicted p max-min : 3.211e+11 2339
max-min gamma: 0.999662 0
Phase-change corrected p max-min : 3.211e+11 2339
max(U) 3.29289e+10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000190233, Final residual = 1.5692e-12, No Iterations 2
bounding epsilon, min: -1.08551e+25 max: 1.06046e+32 average: 4.36063e+27
DILUPBiCG:  Solving for k, Initial residual = 0.033393, Final residual = 8.4911e-09, No Iterations 13
bounding k, min: -4.80646e+12 max: 9.22019e+18 average: 6.04779e+14
ExecutionTime = 317929 s  ClockTime = 332337 s

phiv Courant Number mean: 0.000162677 max: 0.125743 acoustic max: 1.61153e-07
deltaT = 3.11372e-16
Time = 0.0021332566986507

DILUPBiCG:  Solving for rho, Initial residual = 0.00010015, Final residual = 5.72728e-09, No Iterations 1
max-min rho: 146826 0.350586
max-min gamma: 0.999662 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.999073, Final residual = 9.81091e-11, No Iterations 22
DILUPBiCG:  Solving for Uy, Initial residual = 0.998753, Final residual = 4.58541e-09, No Iterations 21
DILUPBiCG:  Solving for Uz, Initial residual = 0.999255, Final residual = 1.20042e-10, No Iterations 22
max(U) 2.32822e+14
GAMG:  Solving for p, Initial residual = 0.515752, Final residual = 9.88678e-18, No Iterations 1
Predicted p max-min : 5.14115e+11 -1.06934e+12
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 6.25708e+11 2339
max(U) 2.32822e+14
GAMG:  Solving for p, Initial residual = 0.974, Final residual = 2.84013e-17, No Iterations 1
Predicted p max-min : 8.43301e+13 -1.25483e+14
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 1.12197e+14 2339
max(U) 2.32822e+14
GAMG:  Solving for p, Initial residual = 0.997598, Final residual = 1.19221e-17, No Iterations 1
Predicted p max-min : 1.72386e+16 -7.189e+15
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 1.72386e+16 2339
max(U) 2.32822e+14
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0236629, Final residual = 6.00719e-10, No Iterations 3
bounding epsilon, min: -4.9582e+29 max: 1.07944e+32 average: 4.7102e+27
DILUPBiCG:  Solving for k, Initial residual = 0.995567, Final residual = 7.02795e-11, No Iterations 17
bounding k, min: -8.52859e+24 max: 2.03241e+26 average: 1.72317e+20
ExecutionTime = 317935 s  ClockTime = 332343 s

phiv Courant Number mean: 0.456991 max: 772.891 acoustic max: 1.60201e-07
deltaT = 5.03584e-20
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1010
    Increased the timePrecision from 14 to 16 to distinguish between timeNames at time 0.00213326
Time = 0.002133256698650698

DILUPBiCG:  Solving for rho, Initial residual = 0.00357835, Final residual = 5.35935e-09, No Iterations 1
max-min rho: 7.82735e+09 -365621
max-min gamma: 0.999913 0
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 7.61395, No Iterations 1001
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 6481.94, No Iterations 1001
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 9.76108e+06, No Iterations 1001
max(U) 1.46542e+40
GAMG:  Solving for p, Initial residual = 1, Final residual = 6.0377e-18, No Iterations 1
Predicted p max-min : 1.44697e+40 -1.54529e+40
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 1.44697e+40 2339
max(U) 5.00447e+39
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0207065, No Iterations 1
Predicted p max-min : 1.21434e+64 -1.24781e+64
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 1.21434e+64 2339
max(U) 6.33122e+42
GAMG:  Solving for p, Initial residual = 0.999737, Final residual = 6.40895e-05, No Iterations 1000
Predicted p max-min : 2.65634e+78 -1.85653e+79
max-min gamma: 0.999913 0
Phase-change corrected p max-min : 2.65634e+78 2339
max(U) 8.91897e+52
DILUPBiCG:  Solving for epsilon, Initial residual = 1, Final residual = 4.96062e-09, No Iterations 128
bounding epsilon, min: -5.12582e+96 max: 2.09237e+101 average: 4.43923e+95
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 5.83008e-09, No Iterations 206
bounding k, min: -8.35059e+126 max: 8.35177e+126 average: 2.52563e+117
ExecutionTime = 318097 s  ClockTime = 332511 s

phiv Courant Number mean: 2.86651e+32 max: 3.58169e+36 acoustic max: 2.59093e-11
deltaT = 1.75749e-57
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1010
    Increased the timePrecision from 16 to 17 to distinguish between timeNames at time 0.00213326
Time = 0.0021332566986506977

DILUPBiCG:  Solving for rho, Initial residual = 4.32193e-11, Final residual = 4.32193e-11, No Iterations 0
max-min rho: 1.20598e+72 0.1
max-min gamma: 0.999913 0
[8] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #3  void Foam::multiply<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[8] #4  Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator*<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[8] #5  
[8]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[8] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #7  
[8]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
[efs-desktop:03943] *** Process received signal ***
[efs-desktop:03943] Signal: Floating point exception (8)
[efs-desktop:03943] Signal code:  (-6)
[efs-desktop:03943] Failing at address: 0x3e800000f67
[efs-desktop:03943] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364c0) [0x7fd6e3f134c0]
[efs-desktop:03943] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fd6e3f13445]
[efs-desktop:03943] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364c0) [0x7fd6e3f134c0]
[efs-desktop:03943] [ 3] cavitatingFoam(_ZN4Foam8multiplyINS_6VectorIdEENS_13fvsPatchFieldENS_11surfaceMeshEEEvRNS_14GeometricFieldIT_T0_T1_EERKNS5_IdS7_S8_EERKS9_+0x66) [0x4493e6]
[efs-desktop:03943] [ 4] cavitatingFoam(_ZN4FoammlINS_6VectorIdEENS_13fvsPatchFieldENS_11surfaceMeshEEENS_3tmpINS_14GeometricFieldIT_T0_T1_EEEERKNS6_IdS8_S9_EERKSB_+0x1d0) [0x462c10]
[efs-desktop:03943] [ 5] cavitatingFoam() [0x41fc18]
[efs-desktop:03943] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fd6e3efe76d]
[efs-desktop:03943] [ 7] cavitatingFoam() [0x42724d]
[efs-desktop:03943] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 8 with PID 3943 on node efs-desktop exited on signal 8 (Floating point exception).

wyldckat · September 7, 2013, 08:00

Hi George,

If it has ran for 200 iterations this time, then you better have a look at the results, so visually inspect where the problem is coming from.

Here's an example where visually inspecting the results was the only way to figure out what was happening wrong: http://www.cfd-online.com/Forums/ope...tml#post446350 - post #17

Best regards,
Bruno

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries	NickG	OpenFOAM Installation	3	December 30, 2019 00:21
c++ libraries and solver compiling	vaina74	OpenFOAM Installation	13	February 3, 2012 17:43
How to get the max value of the whole field	waynezw0618	OpenFOAM Running, Solving & CFD	4	June 17, 2008 05:07
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 17:51
error while compiling the USER Sub routine	CFD user	CFX	3	November 25, 2002 15:16

September 1, 2013, 09:34		#2
wyldckat Retired Super Moderator Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 10,975 Blog Entries: 45 Rep Power: 128	Hi George, There is some information missing from the crash, because there is no indication of which signal was triggered. Either way, it does give an indication that something went wrong in the preconditioner, which suggests that either: The equation/matrix solver itself should be different or reconfigured in "fvSolution". Or that the k-epsilon or some other field went completely in the wrong direction, for some strange reason. My suggestion is to try and see the results from the 1st and 2nd iterations, so that you can see the values for the k-epsilon fields and figure out if they have values that are very large or very small. Best regards, Bruno __________________ OpenFOAM: FAQ \| Getting started Forum: How to get help, to post code/output and forum guide Read this before sending me PM

September 7, 2013, 08:00		#4
wyldckat Retired Super Moderator Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 10,975 Blog Entries: 45 Rep Power: 128	Hi George, If it has ran for 200 iterations this time, then you better have a look at the results, so visually inspect where the problem is coming from. Here's an example where visually inspecting the results was the only way to figure out what was happening wrong: http://www.cfd-online.com/Forums/ope...tml#post446350 - post #17 Best regards, Bruno __________________ OpenFOAM: FAQ \| Getting started Forum: How to get help, to post code/output and forum guide Read this before sending me PM