CFD Online Discussion Forums - initial condition of kappai at IATE diameter model

Hi,

When IATE diameter model is used for gas phase on phaseProperties, it calls kappai dictionary for boundary and initial condition.
So if I use this model I've to modifie kappai.gas dictionary. There is a description on IATE.H file for kappai and it says:

Code:

 "Interfacial curvature (alpha*interfacial area)"

Unit is m^-1

Also, there is an equation on .C file:

Code:

// Construct the interfacial curvature equation

    fvScalarMatrix kappaiEqn

    (

        fvm::ddt(kappai_) + fvm::div(phase_.phi(), kappai_)

      - fvm::Sp(fvc::div(phase_.phi()), kappai_)

     ==

      - fvm::SuSp(R, kappai_)

    //+ Rph() // Omit the nucleation/condensation term

      + fvOptions(kappai_)

    );

Then, I checked tutorials that uses IATE model to see how they initialize and which boundaries used and I saw some strange numbers on Initial Condition like 13333 and 2000 (Interfacial curvature (alpha*interfacial area) ???).
How they obtain these numbers? I couldn't reached to the original paper that referred on source code (Ishii, M., Kim, S. and Kelly, J., (2005)). Maybe in the paper there is an explanation. But I found some others that might be related to Ishii et al. (2005).

Code:

Ishii, M., Kim, S., & Uhle, J. (2002). Interfacial area transport equation: model development and benchmark experiments. International Journal of Heat and Mass Transfer, 45(15), 3111-3123. https://doi.org/10.1016/S0017-9310(02)00041-8.

--------------------------------------------------------

Kocamustafaogullari, G., & Ishii, M. (1995). Foundation of the interfacial area transport equation and its closure relations. International Journal of Heat and Mass Transfer, 38(3), 481-493. https://doi.org/10.1016/0017-9310(94)00183-V.

However, I still don't have a clue. Any theorical and mathematical explanation would be appretiated.

Best regards,
Saidc.