MPPIC interfoam 依据的动量方程是那一个呢？

qingdong_wang

 fvVectorMatrix UEqn
 (
     fvm::ddt(alphacRho, U)
   + MRF.DDt(alphacRho, U)
   - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
   + fvm::div(rhoPhi, U)
   + turbulence->divDevRhoReff(U)
   ==
     fvOptions(rho, U)
   + cloudSU
 );

这是MPPIC Intefoam 1912里面的动量方程matrix，但是这一行代表什么意思呢

 - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

我看的动量方程是
fluent的theory也是这个，所以我觉得方程是不是

fvVectorMatrix UEqn
    (
      //  fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
      //+ MRF.DDt(rho, U)
     fvm::div(alphaRhoPhic, U)
     + fvm::ddt(alphacRho, U)
     + MRF.DDt(alphacRho, U)
      + turbulence->divDevRhoReff(U)
     ==
        fvOptions(rho, U)
        + cloudSU
    );
....

 -alphac * ghf*fvc::snGrad(rho)
- alphac * fvc::snGrad(p_rgh)

我没有看懂原先方程的意义，或者说他是怎么推倒出来的呢？如果根据文献的这个方程的话。

星星星星晴

代表着source term 和fluent里面的cloudSU是一样的

星星星星晴

对不起看错了

李东岳

fvVectorMatrix UEqn
 (
     fvm::ddt(alphacRho, U)//
   + MRF.DDt(alphacRho, U)
   - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
   + fvm::div(rhoPhi, U)
   + turbulence->divDevRhoReff(U)
   ==
     fvOptions(rho, U)
   + cloudSU
 );

这个方程时间步里面考虑了相分数，但是对流项等并没有考虑，你试过改成下面这种么？

fvVectorMatrix UEqn
 (
     fvm::ddt(rho, U)//
   + MRF.DDt(alphacRho, U)
   - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
   + fvm::div(rhoPhi, U)
   + turbulence->divDevRhoReff(U)
   ==
     fvOptions(rho, U)
   + cloudSU
 );

qingdong_wang

我找着我第一次发的方程改得，就是第二个block的代码，算出来的结果和MPPICInterfoam的基本一样。当然有可能我的粒子直径比较小，对计算结果影响不大？

我没明白您下面写的这个是什么意思。如果不看mrf的话，你这个和interfoam相比就是减上了一项

fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

这一项就是
我不知道怎么推出的这一项？