interPhaseChangeFoam相方程

[飞火流星jyj]

如何理解最后两项？？？
是fvm::Sp(vDotvmcAlphal,alpha1)+vDotcAlphal 不是单一的vDotvmcAlphal，自己看代码好像vDotvmcAlphal已经代表了m*(1/rhol-alpha1*(1/rhol-1/rhov)),这不就是公式变形后的最后一项吗？求老师、大神指导

[CYW]

@飞火流星jyj 你知道mDot 和mDotAlpha 的区别吗

[飞火流星jyj]

@cyw 代码里只有mDotAlpha（）质量源项啊，没有mDot

[CYW]

@飞火流星jyj 啊，不好意思，我说的是另外一个solver，interCondensatingEvaporatingFoam.

[CYW]

@飞火流星jyj 我的理解是
vDotvmcAlphal = vDotvAlpha - vDotcAlpha (凝结与汽化的差)
而相方程的mDot = vDotvAlpha，故又加上了vDotcAlpha。
不知道这样理解对不对

[飞火流星jyj]

@cyw 凝结和气化之差不就已经是体积源了吗，为啥还加vDotcAlpha

[yuan_neu]

这个问题我之前也琢磨过，

加vDotcAlpha是因为后面这个项原型应该为（1-alpha）*vDotcAlpha

其实在湍流模型的k方程里面也有类似的处理，我个人偏向于以下这个结论

尽量对各项做隐性处理，这样可以保证方程的稳定性，

在压力方程里面，mass transfer rate 这一项也是做隐性处理的哦