CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Why RNGkepsilon model gives floating error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 31, 2013, 21:19
Default Why RNGkepsilon model gives floating error
  #1
Senior Member
 
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 6
shipman is on a distinguished road
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;
  1. I checked the mesh quality which was ok.
  2. set the initial condition of the k and epsilon properl
  3. 2. I decreased courant number to decrease time step.
  4. 3. Additionaly decreased time step too.
  5. 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;
}
shipman is offline   Reply With Quote

Old   September 1, 2013, 09:34
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,816
Blog Entries: 39
Rep Power: 108
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
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
__________________
wyldckat is offline   Reply With Quote

Old   September 4, 2013, 22:36
Default hi
  #3
Senior Member
 
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 6
shipman is on a distinguished road
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).
shipman is offline   Reply With Quote

Old   September 7, 2013, 08:00
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,816
Blog Entries: 39
Rep Power: 108
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
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
__________________
wyldckat is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
compile error in dynamicMesh and thermophysicalModels libraries NickG OpenFOAM Installation 2 August 30, 2013 07:42
c++ libraries and solver compiling vaina74 OpenFOAM Installation 13 February 3, 2012 18: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 16:16


All times are GMT -4. The time now is 08:22.