[sabrinacarson]

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-24, cumulative = 8.21085e-24
GAMG:  Solving for pcorr, Initial residual = 0.434647, Final residual = 5.36898e-05, No Iterations 7
time step continuity errors : sum local = 8.03698e-22, global = -2.9247e-23, cumulative = -2.10362e-23
Courant Number mean: 0.00060848 max: 0.000819

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0.00060848 max: 0.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-22, cumulative = 8.81193e-22
GAMG:  Solving for pcorr, Initial residual = 0.983295, Final residual = 9.04202e-05, No Iterations 7
time step continuity errors : sum local = 2.74338e-20, global = 7.31173e-22, cumulative = 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.sh: line 37: 27615 Floating point exception(core dumped) phaseChangeHeatDyMFoam 


The boundary conditions, fvschemes and fvsolution are given below in folder:
Click here : https://drive.google.com/open?id=1IT...FO7H-EDNlHlXWj

[regard]

Originally Posted by sabrinacarson

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-24, cumulative = 8.21085e-24
GAMG:  Solving for pcorr, Initial residual = 0.434647, Final residual = 5.36898e-05, No Iterations 7
time step continuity errors : sum local = 8.03698e-22, global = -2.9247e-23, cumulative = -2.10362e-23
Courant Number mean: 0.00060848 max: 0.000819

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0.00060848 max: 0.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-22, cumulative = 8.81193e-22
GAMG:  Solving for pcorr, Initial residual = 0.983295, Final residual = 9.04202e-05, No Iterations 7
time step continuity errors : sum local = 2.74338e-20, global = 7.31173e-22, cumulative = 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.sh: line 37: 27615 Floating point exception(core dumped) phaseChangeHeatDyMFoam 


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?