# Kunz phase change model and what is done

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 February 24, 2010, 05:37 Kunz phase change model and what is done #1 Senior Member     Roman Thiele Join Date: Aug 2009 Location: London, UK Posts: 368 Rep Power: 13 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

 February 24, 2012, 03:08 mDotp() and mDotAlphal() #2 Member   vahid Join Date: Feb 2012 Location: Mashhad-Iran Posts: 80 Rep Power: 5 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?

 February 28, 2012, 09:03 #3 Senior Member   David Boger Join Date: Mar 2009 Location: Penn State Applied Research Laboratory Posts: 146 Rep Power: 9 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

March 1, 2012, 20:25
myinterPhaseChangeFoam
#4
Member

vahid
Join Date: Feb 2012
Posts: 80
Rep Power: 5
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
 SinghalinterPhaseChangeFoam.zip (22.2 KB, 20 views)

 June 15, 2015, 11:44 #5 Senior Member   Przemek Join Date: Jun 2011 Posts: 133 Rep Power: 7 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

All times are GMT -4. The time now is 08:41.

 Contact Us - CFD Online - Top