# OpenFoam Stopped without error message twoPhaseEulerFoam

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

 October 16, 2020, 00:33 OpenFoam Stopped without error message twoPhaseEulerFoam #1 New Member   Sumit Peh Join Date: Oct 2018 Location: Beijing Posts: 20 Rep Power: 7 Hi Foamer ! I searched and solved many errors on my case in this community, but I didn't found the solution on this case. I'm running modified bubble column case with my own geometry. I used the setup case from bubble column tutorial in openfoam, mesh by the cfmesh with .fms geometry file and modified the meshDict, setFieldsDict, g(correct z axis),d0(bubble size), controlDict(Endtime 60s). Everything seems to run smoothly on the cartesianMesh, SetFields and even the twoPhaseEulerFoam. But it randomly stopped at "Constructing momentum equations" here is the end part of the terminal. Code: ```Courant Number mean: 0.192367 max: 54.0446 Max Ur Courant Number = 34.4552 Time = 6.825 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.251251 Min(alpha.air) = -1.23474e-20 Max(alpha.air) = 1.06091 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.993034, Final residual = 5.72382e-05, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.900522, Final residual = 0.00115984, No Iterations 1000 min T.air 296.75 min T.water 299.824 GAMG: Solving for p_rgh, Initial residual = 0.00734846, Final residual = 3.66422e-09, No Iterations 9 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.251258 Min(alpha.air) = -2.51025e-21 Max(alpha.air) = 1.02043 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.29879, Final residual = 0.0027419, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.665088, Final residual = 0.00319244, No Iterations 1000 min T.air 296.691 min T.water 299.824 GAMG: Solving for p_rgh, Initial residual = 0.00326417, Final residual = 4.35379e-09, No Iterations 8 PIMPLE: iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.251249 Min(alpha.air) = -7.92949e-21 Max(alpha.air) = 1.20173 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.21708, Final residual = 0.00331834, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.178944, Final residual = 0.00354534, No Iterations 1000 min T.air 295.475 min T.water 301.303 GAMG: Solving for p_rgh, Initial residual = 0.00583616, Final residual = 8.10231e-09, No Iterations 9 bounding km, min: 9.94883e-16 max: 19599.2 average: 50.575 bounding epsilonm, min: 2.0515e-21 max: 5.27789e+07 average: 30345.7 smoothSolver: Solving for epsilonm, Initial residual = 0.000395396, Final residual = 8.97474e-08, No Iterations 22 bounding epsilonm, min: -9904.44 max: 5.27789e+07 average: 39637.7 smoothSolver: Solving for km, Initial residual = 2.44555e-09, Final residual = 3.26389e-10, No Iterations 1 bounding km, min: -0.562915 max: 16106.6 average: 153.064 ExecutionTime = 115.67 s ClockTime = 121 s fieldAverage fieldAverage1 write: Calculating averages Courant Number mean: 0.447823 max: 54.8887 Max Ur Courant Number = 369.284 Time = 6.83 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.251188 Min(alpha.air) = -3.96738e-21 Max(alpha.air) = 1.36108 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.194015, Final residual = 0.00185957, No Iterations 1000 smoothSolver: Solving for e.water, Initial residual = 0.00960397, Final residual = 9.90323e-08, No Iterations 328 min T.air 191.056 min T.water 301.337 GAMG: Solving for p_rgh, Initial residual = 0.0732878, Final residual = 4.43483e-09, No Iterations 11 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.251121 Min(alpha.air) = -7.1303e-22 Max(alpha.air) = 4.08928 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 0.981585, Final residual = 9.4538e-08, No Iterations 126 smoothSolver: Solving for e.water, Initial residual = 0.0150729, Final residual = 6.45716e-07, No Iterations 1000 min T.air -2387.77 min T.water 286.602 GAMG: Solving for p_rgh, Initial residual = 0.60359, Final residual = 5.60784e-09, No Iterations 17 PIMPLE: iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = -1.31932 Min(alpha.air) = -3583.82 Max(alpha.air) = 11737.8 Constructing momentum equations jameslinux@DESKTOP-XXXXXX:~/Practice_files/bbccaselessmesh\$``` I don't know the reason know why it stopped without any error message, please guide me the reason. ps. my drive space still has 30 GB left, if that the reason. Thank you Foamer

 October 16, 2020, 02:05 #2 Senior Member   Yogesh Bapat Join Date: Oct 2010 Posts: 102 Rep Power: 15 Hello, If you look at log and volume fraction values at the end min value is negative and max value is much greater than 1. The case has diverged and may have got a floating point exception.

October 16, 2020, 08:57
#3
New Member

Sumit Peh
Join Date: Oct 2018
Location: Beijing
Posts: 20
Rep Power: 7
Quote:
 Originally Posted by ybapat Hello, If you look at log and volume fraction values at the end min value is negative and max value is much greater than 1. The case has diverged and may have got a floating point exception.
Thanks for the replied. I solved this problem by changed the magnitude of U.air to Z-axis. the bubble column tutorial geometry are in the Y-axis length.
but yet, I still got another problem to fixed on Negative initial temperature T0.
I'll try to fix it and post the solution on my next thread if i possible to manage it.

remarks for others : check the g, U.air in the correct axis // p p_rho with 1e5.

 Tags cfd, crash, openfoam, simulation cfd, twophaseeulerfoam