Kunz phase change model and what is done
I have a couple of questions about how the phase change models in interPhaseChangeFoam work.
There are two functions within the solver mDotP() and mDotAlphal()
mDotP()'s description is like this
The same goes for mDotAlphal()
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 (which is mass condensation) is dimensionless and for mDotAlphal, 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).
mDotp() and mDotAlphal()
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?
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.
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.
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.
|All times are GMT -4. The time now is 09:36.|