[cyw]

Dear all,

Does anyone use this solver? I wonder if anyone could help me to understand this solver.

in this solver mDot equals
coeffC_*mixture_.rho2()*limitedAlpha2*max(TSat - T.oldTime(), T0),
coeffE_*mixture_.rho1()*limitedAlpha1*max(T.oldTim e() - TSat, T0)

mDotAlpha equals
coeffC_*mixture_.rho2()*max(TSat - T.oldTime(), T0),
-coeffE_*mixture_.rho1()*max(T.oldTime() - TSat, T0)

That is to say, mDotAlpha*Alpha=mDot? And mDotAlpha is used in alphaEqn.H. Why not mDot used in alphaEqn.H?

If i implemented Schrage phase change model, how to calculate mDotAlpha?