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

OpenFoam Stopped without error message twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 16, 2020, 01:33
Default OpenFoam Stopped without error message twoPhaseEulerFoam
  #1
New Member
 
Sumit Peh
Join Date: Oct 2018
Location: Beijing
Posts: 20
Rep Power: 8
Jamessmp23 is on a distinguished road
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
Jamessmp23 is offline   Reply With Quote

Old   October 16, 2020, 03:05
Default
  #2
Senior Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 102
Rep Power: 16
ybapat is on a distinguished road
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.
ybapat is offline   Reply With Quote

Old   October 16, 2020, 09:57
Default
  #3
New Member
 
Sumit Peh
Join Date: Oct 2018
Location: Beijing
Posts: 20
Rep Power: 8
Jamessmp23 is on a distinguished road
Quote:
Originally Posted by ybapat View Post
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.
Jamessmp23 is offline   Reply With Quote

Reply

Tags
cfd, crash, openfoam, simulation cfd, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM course for beginners Jibran OpenFOAM Announcements from Other Sources 2 November 4, 2019 09:51
twoPhaseEulerFoam validation OpenFOAM 4.1 jiejie OpenFOAM Running, Solving & CFD 4 April 20, 2019 04:06
OpenFOAM 4.0 Released CFDFoundation OpenFOAM Announcements from OpenFOAM Foundation 2 October 6, 2017 06:40
OpenFOAM Training: Programming CFD Course 12-13 and 19-20 April 2016 cfd.direct OpenFOAM Announcements from Other Sources 0 January 14, 2016 11:19
Modified OpenFOAM Forum Structure and New Mailing-List pete Site News & Announcements 0 June 29, 2009 06:56


All times are GMT -4. The time now is 17:49.