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



  •  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)
    

    我看的动量方程是08ade4b6-514f-4172-bdc2-f803b29b71a3-image.png
    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
     );
    


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

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

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

    这一项就是65293cb6-a527-49eb-bdc4-2355ec82c507-image.png
    我不知道怎么推出的这一项?


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2