@李东岳 在 SIMPLE(C)解一次動量方程後多次求解壓力方程，壓力方程中計算HbyA的疑問 中说：

pisoFoam

因爲我看到動量預測不是必須的，那我想我能不能某幾步求解一下壓力方程再做一次動量預測，因爲有時候動量方程的增量其實很小，解一次壓力方程對應一次動量方程，是不是有點浪費。

引用自李老師的simple解析。
如果現在每1部做一次動量預測，然後求解20次壓力方程，20次壓力方程的每一步是不是都需要構建新的hybA？ 構建HbyA 需要用到An，同時壓力方程的系數需要用到Ap，是否每一步都需要去更新一次Ap與An？

您可以看看苏州舜云工程软件有限公司哦，主要业务就是做基于粒子法的流体多物理场数值仿真软件研发的，就在苏州市相城区，可以关注公众号：舜云多物理场仿真，了解企业的更多信息，简历投递邮箱是lijuan.cai@shoncloud.com

我也很纠结, 我问过好几个人, 他们也不清楚. 不过我的算例中, 乘和不乘对结果影响不大. 我的是dilute flow

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

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

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

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

 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)

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

@zhe 我能加你一下微信吗？我最近也在做这个问题。因为我直接从interfoam改MPPICInterfoam， 所以我觉得还是有些问题。我的微信是 wqdwqdaaa

谢谢你的指正。
不过我没看湍流模型， 我用的DNS， 不过看湍流的构造函数的时候，他们的确是不同的。

MPPICInterFoam 在左侧加了

fvm : : ddt ( alphacRho ,  U)+ MRF. DDt( alphacRho ,  U）

应该是 When a particle is transformed into a Eulerian structure, its velocity and volume fraction is added to the continuum fluid

右侧是the effect of particle motion on the momentum of surrounding fluid。 也即是我们说的粒子力。

phicForces / rAUcf   //   is   added   as   particle   force 

理论上 ， 当粒子为0的时候， 由于alphac为0， 这两项也应为为0才对。理论上是没有影响的。
不过MPPICInterFoam的alphaEqn改了很多， 说实话，这个我没有看明白。
不过你这两个case是除了solver不一样外， 其他全部一样吗？比如说fvSchemes或者fvSolutions等等

你可以看看下面这个文档， 虽然他是用的openfoam3， 但是感觉会对你有帮助

pg11 http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2016/EbrahimGhahramani/Report_Ghahramani.pdf

然后还有一点我不太明白的地方， 就是这个地方的这个地方的粒子力为啥先乘了rAUc然后又除了rAUc？好像是关于离散吗？

volScalarField rAUc(1.0/UEqn.A());
        surfaceScalarField rAUcf(fvc::interpolate(rAUc));


        surfaceScalarField phicForces
        (
                (fvc::interpolate(rAUc*cloudVolSUSu) & mesh.Sf())
        );

    if (pimple.momentumPredictor())
    {
        solve
        (
                        UEqn
                 ==
                        fvc::reconstruct
                        (
                                phicForces/rAUcf // is added as particle force
+ ....