CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    3
    5
    907
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • Q
      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)
      

      我看的动量方程是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)
      
      
      

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

      1 条回复 最后回复 回复 引用
      • 星
        星星星星晴 最后由 编辑

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

        m.sui20@foxmail.com

        1 条回复 最后回复 回复 引用
        • 星
          星星星星晴 最后由 编辑

          对不起 看错了

          m.sui20@foxmail.com

          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 编辑

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

            CFD高性能服务器 http://dyfluid.com/servers.html
            2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

            1 条回复 最后回复 回复 引用
            • Q
              qingdong_wang 最后由 编辑

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

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

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

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

              1 条回复 最后回复 回复 引用
              • First post
                Last post