[romant]

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

//- 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()

//- 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).