[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

[nimasam]

Dear Sabrina

it seems the error is related to the library smoothInterfaceProperties which is a modified version of interfaceProperties. Therefore, i recommend to use the original library of OpenFOAM, and go forward.
smoothInterfaceProperties as i mentioned in my published paper is to eliminate parasite current. it is not vital for your implementation.

[sabrinacarson]

Originally Posted by nimasam

Dear Sabrina

it seems the error is related to the library smoothInterfaceProperties which is a modified version of interfaceProperties. Therefore, i recommend to use the original library of OpenFOAM, and go forward.
smoothInterfaceProperties as i mentioned in my published paper is to eliminate parasite current. it is not vital for your implementation.

Thanks. It actually worked. I used interfaceProperties.H instead of smoothInterfaceProperties.H

But I am getting spurious currents around the interface. Interface is not smooth, (smoothinterfaceproperties also gave spurious currents kind of result in phaseChangeHeatDyMFoam)
Please find the attachment to see the results.

What can be the solution?

[nimasam]

i guess you may need to update phi and rhoPhi after mesh refinement, it may helps you to get better result, it is far more than parasite current problem

[sabrinacarson]

Hi,

Thanks for your reply. I was trying few other ways to make the solver correct, but didn't work.

How the values of phi and rhoPhi will be updated? In which particular file these values are being updated?