[Jamessmp23]

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.

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

[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.

[Jamessmp23]

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.