|
[Sponsors] |
August 31, 2013, 22:19 |
Why RNGkepsilon model gives floating error
|
#1 |
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 13 |
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).
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; } } 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; } 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; } 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; } |
|
September 1, 2013, 10:34 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
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:
Best regards, Bruno
__________________
|
|
September 4, 2013, 23:36 |
hi
|
#3 |
Senior Member
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 130
Rep Power: 13 |
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). |
|
September 7, 2013, 09:00 |
|
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
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
__________________
|
|
|
|
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 01:21 |
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 06:07 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |
error while compiling the USER Sub routine | CFD user | CFX | 3 | November 25, 2002 16:16 |