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

Error while running new solver phaseChangeHeatDyMFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 1, 2018, 03:44
Default Error while running new solver phaseChangeHeatDyMFoam
  #1
New Member
 
Sabrina Carson
Join Date: May 2018
Posts: 15
Rep Power: 8
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   June 10, 2018, 05:31
Default
  #2
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
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.
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Old   June 13, 2018, 07:59
Default
  #3
New Member
 
Sabrina Carson
Join Date: May 2018
Posts: 15
Rep Power: 8
sabrinacarson is on a distinguished road
Quote:
Originally Posted by nimasam View Post
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?
Attached Images
File Type: png Slide1.PNG (191.3 KB, 24 views)

Last edited by sabrinacarson; June 14, 2018 at 00:13.
sabrinacarson is offline   Reply With Quote

Old   June 14, 2018, 13:41
Default
  #4
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
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
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Old   July 10, 2018, 07:28
Default
  #5
New Member
 
Sabrina Carson
Join Date: May 2018
Posts: 15
Rep Power: 8
sabrinacarson is on a distinguished road
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?
sabrinacarson is offline   Reply With Quote

Reply

Tags
error, new solver, 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
New solver not running in loop. silviliril OpenFOAM Running, Solving & CFD 2 May 21, 2018 08:55
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
Error Message when running solver Jiang CFX 4 August 12, 2013 13:37
Statically Compiling OpenFOAM Issues herzfeldd OpenFOAM Installation 21 January 6, 2009 09:38
Kubuntu uses dash breaks All scripts in tutorials platopus OpenFOAM Bugs 8 April 15, 2008 07:52


All times are GMT -4. The time now is 09:18.