sabrinacarson |
May 27, 2018 02:01 |
Error while running annular boiling case
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
|