fvVectorMatrix Ueqn

qingdong_wang

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

MPPICInterfoam 中的UEqn == 后面是源项。但是我不太明白的是，为啥这个matrix中间有==呢？因为简单的icoFoam中没有后面的源项，就会感觉比较容易理解，但是加上了等号，我就不能理解了。
还有就是这个cloudSU也是源项吗还是ö粒子力？这个方程在后面不是已经加上particle forces了吗？

UEqn
                 ==
                        fvc::reconstruct
                        (
                                phicForces/rAUcf // is added as particle force
                          +
                                (
                                        fvc::interpolate
                                        (
                                        mixture.sigmaK()
                                        )*fvc::snGrad(alpha1) // can be supplanted by the mixture.surfaceTensionForce()
                                        - ghf*fvc::snGrad(rho)
                                        - fvc::snGrad(p_rgh)
                                ) * mesh.magSf()
            )

李东岳

矩阵求解是Ax=b，等号后面的是b
源项处理一部分隐性，一部分显性，phicForce是显性的部分

qingdong_wang

矩阵是Ax = b, 然后solve（UEqn = fvc：xxx）是求解的Ax=b？那么这样的话A 是前面的UEqn （fvVectorMatrix）？那这个里面的matrix为啥还有等号呢

qingdong_wang

然后还有一点我不太明白的地方，就是这个地方的这个地方的粒子力为啥先乘了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
+ ....

李东岳

pEqn里面用到了rAUc