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

Kunz phase change model and what is done

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 24, 2010, 05:37
Default Kunz phase change model and what is done
  #1
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Hej,

I have a couple of questions about how the phase change models in interPhaseChangeFoam work.

First:
There are two functions within the solver mDotP() and mDotAlphal()
mDotP()'s description is like this
Code:
//- Return the mass condensation and vaporisation rates as an
//  explicit term for the condensation rate and a coefficient to
//  multiply (p - pSat) for the vaporisation rate
What does this exactly mean?

The same goes for mDotAlphal()
Code:
//- Return the mass condensation and vaporisation rates as a
//  coefficient to multiply (1 - alphal) for the condensation rate
//  and a coefficient to multiply  alphal for the vaporisation rate
What does this exactly mean?

I tried following the steps within the solver and found that for mDotAlphal, the dimensions of the value pair returned are different from anothe.r mDotAlphal[0] (which is mass condensation) is dimensionless and for mDotAlphal[1], the dimension is [kg/s^2/m]. The same is true for the mDotP().

Now, there is a read() function build into the the class as well. As far as I can see, it is never really called, but when following the code if read() would be called, then some of the coefficients are defined differently (namely they include dimensions instead of no dimensions). So, if this read() is called, then the dimensions for all returned values are [kg/m^3/s].

Does anyone know if read() is called while running interPhaseChangeSolver? And what is supposed to be returned by the functions (dimension wise).
__________________
~roman
romant is offline   Reply With Quote

Old   February 24, 2012, 03:08
Default mDotp() and mDotAlphal()
  #2
Member
 
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 4
vahid.najafi is an unknown quantity at this point
Hello OpenFoam users ,
I'm studying on the Kunz model.we know in this model we have two term for mass dest and prod.(m+ and m-)
What mean mDotAlphal() and mDotP() In interphasechange?
What is the difference between these two?
vahid.najafi is offline   Reply With Quote

Old   February 28, 2012, 09:03
Default
  #3
Senior Member
 
David Boger
Join Date: Mar 2009
Location: Penn State Applied Research Laboratory
Posts: 146
Rep Power: 8
boger is on a distinguished road
mDotP and mDotAlpha each have to do with the way that the mass transfer terms are linearized in the p-equation and alpha-equation, respectively. The best way to understand each is to write out the equations in which they are used, looking at pEqn.H and alphaEqn.H. Parts of the terms are treated implicitly, and other parts explicitly. Note that pos and neg are effectively the derivatives of the max and min terms, respectively.
__________________
David A. Boger
boger is offline   Reply With Quote

Old   March 1, 2012, 20:25
Default myinterPhaseChangeFoam
  #4
Member
 
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 4
vahid.najafi is an unknown quantity at this point
Dear David ,
My purpose is adding Zwart and Singhal models to interPhaseChangeFoam solver.
I transform Saur model to Zwart model successfully but transforming Saur model to Singhal model is faild!!!
my error is :

phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C:48: error: no matching function for call to ‘Foam::dimensioned<double>::dimensioned()’

I hope your help to solve this problem.
Thanks a lot.
Attached Files
File Type: zip SinghalinterPhaseChangeFoam.zip (22.2 KB, 13 views)
vahid.najafi is offline   Reply With Quote

Old   June 15, 2015, 11:44
Default
  #5
Member
 
Przemek
Join Date: Jun 2011
Posts: 42
Rep Power: 6
gaza is on a distinguished road
Hi
pos and neg functions are defined in scalar.H file
they are not derivatives! They return 0 or 1 value
Regarding mDotAlphal and mDotP functions they are implemented for cavitation models.
They are inappropriate for heat transfer mass transfer models (condensation, evaporation).
One have to change a bit the code (places where these functions) are used and write own functions as well.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On



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