CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Error while running annular boiling case

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 27, 2018, 02:01
Default Error while running annular boiling case
  #1
New Member
 
Sabrina Carson
Join Date: May 2018
Posts: 15
Rep Power: 5
sabrinacarson is on a distinguished road
Hello,

I have modified custom solver phaseChangeHeatFoam with dyanmic mesh capabilities. New solver is called phaseChangeHeatDyMFoam.
The solver is working fine when I am working with boiling of static liquid, as the phase is changing the mesh is adapting according to alpha value at the interface.

But, If I am working with annular two phase flow of gas and liuid in cylinder having degree of superheat applied on the wall (fluid is having velocity now, its not static as earlier). as soon as it starts adapting the mesh the solver stops. It shows following error:

PHP Code:
Create time

Create mesh 
for time 0

Selecting dynamicFvMesh dynamicRefineFvMesh
Reading field p_rgh

Reading field T

Reading field alpha1

Reading field U

Reading
/calculating face flux field phi

Creating phaseChangeTwoPhaseMixture

Selecting phaseChange model Hardt
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar

Reading g

Reading hRef
Calculating field g
.h

No finite volume options present

Creating field kinetic energy K


PIMPLE
Operating solver in PISO mode

time step continuity errors 
sum local 1.36275e-21, global = 8.21085e-24cumulative 8.21085e-24
GAMG
:  Solving for pcorrInitial residual 0.434647, Final residual 5.36898e-05No Iterations 7
time step continuity errors 
sum local 8.03698e-22, global = -2.9247e-23cumulative = -2.10362e-23
Courant Number mean
0.00060848 max0.000819

Starting time loop

Interface Courant Number mean0 max0
Courant Number mean
0.00060848 max0.000819
deltaT 
1e-06
Time 
1e-06

Selected 6760 cells 
for refinement out of 27885.
Refined from 27885 to 75205 cells
.
Selected 0 split points out of a possible 6760.
Execution time 
for mesh.update() = 0.35 s
time step continuity errors 
sum local 1.04334e-19, global = 9.02229e-22cumulative 8.81193e-22
GAMG
:  Solving for pcorrInitial residual 0.983295, Final residual 9.04202e-05No Iterations 7
time step continuity errors 
sum local 2.74338e-20, global = 7.31173e-22cumulative 1.61237e-21
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = 0  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -5.31557e-37  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -9.65276e-37  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -5.9622e-37  Max(alpha1) = 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#4  Foam::smoothInterfaceProperties::calculateK() at ??:?
#5  
 
at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 
at ??:?
RUN.shline 3727615 Floating point exception(core dumpedphaseChangeHeatDyMFoam 
The boundary conditions, fvschemes and fvsolution are given below in folder:
Click here : https://drive.google.com/open?id=1IT...FO7H-EDNlHlXWj
sabrinacarson is offline   Reply With Quote

Old   April 19, 2019, 14:39
Default
  #2
Member
 
...
Join Date: May 2018
Posts: 37
Rep Power: 5
regard is on a distinguished road
Quote:
Originally Posted by sabrinacarson View Post
Hello,

I have modified custom solver phaseChangeHeatFoam with dyanmic mesh capabilities. New solver is called phaseChangeHeatDyMFoam.
The solver is working fine when I am working with boiling of static liquid, as the phase is changing the mesh is adapting according to alpha value at the interface.

But, If I am working with annular two phase flow of gas and liuid in cylinder having degree of superheat applied on the wall (fluid is having velocity now, its not static as earlier). as soon as it starts adapting the mesh the solver stops. It shows following error:

PHP Code:
Create time

Create mesh 
for time 0

Selecting dynamicFvMesh dynamicRefineFvMesh
Reading field p_rgh

Reading field T

Reading field alpha1

Reading field U

Reading
/calculating face flux field phi

Creating phaseChangeTwoPhaseMixture

Selecting phaseChange model Hardt
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar

Reading g

Reading hRef
Calculating field g
.h

No finite volume options present

Creating field kinetic energy K


PIMPLE
Operating solver in PISO mode

time step continuity errors 
sum local 1.36275e-21, global = 8.21085e-24cumulative 8.21085e-24
GAMG
:  Solving for pcorrInitial residual 0.434647, Final residual 5.36898e-05No Iterations 7
time step continuity errors 
sum local 8.03698e-22, global = -2.9247e-23cumulative = -2.10362e-23
Courant Number mean
0.00060848 max0.000819

Starting time loop

Interface Courant Number mean0 max0
Courant Number mean
0.00060848 max0.000819
deltaT 
1e-06
Time 
1e-06

Selected 6760 cells 
for refinement out of 27885.
Refined from 27885 to 75205 cells
.
Selected 0 split points out of a possible 6760.
Execution time 
for mesh.update() = 0.35 s
time step continuity errors 
sum local 1.04334e-19, global = 9.02229e-22cumulative 8.81193e-22
GAMG
:  Solving for pcorrInitial residual 0.983295, Final residual 9.04202e-05No Iterations 7
time step continuity errors 
sum local 2.74338e-20, global = 7.31173e-22cumulative 1.61237e-21
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = 0  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -5.31557e-37  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -9.65276e-37  Max(alpha1) = 1
MULES
Solving for alpha1
Liquid phase volume fraction 
0.305544  Min(alpha1) = -5.9622e-37  Max(alpha1) = 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#4  Foam::smoothInterfaceProperties::calculateK() at ??:?
#5  
 
at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 
at ??:?
RUN.shline 3727615 Floating point exception(core dumpedphaseChangeHeatDyMFoam 
The boundary conditions, fvschemes and fvsolution are given below in folder:
Click here : https://drive.google.com/open?id=1IT...FO7H-EDNlHlXWj
hi

is this solver(phasechangeheatfoam) can simulation of subcooled flow boling?
regard is offline   Reply With Quote

Reply

Tags
adaptive meshing, developing new solver, dynamic meshing, openfoam

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
interFoam running blowing up sandy13 OpenFOAM Running, Solving & CFD 2 May 5, 2015 07:16
Running bubbleFoam for turbulent case... vishal3 OpenFOAM 7 March 12, 2013 06:13
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
How to save a case running in background us FLUENT 0 July 6, 2005 10:43
Annular Diffusor Test Case Jonas Larsson Main CFD Forum 0 July 2, 2003 09:34


All times are GMT -4. The time now is 01:12.