compressibleInterFoam pEqn.H

[hongjiewang]

请问有没有前辈了解在compressibleInterFoam求解器中，对可压缩项的这部分处理是怎么弄的，看不太懂~

    p_rghEqnComp1 =
        fvc::ddt(rho1) + psi1*correction(fvm::ddt(p_rgh))
      + fvc::div(phi, rho1) - fvc::Sp(fvc::div(phi), rho1);

    p_rghEqnComp2 =
        fvc::ddt(rho2) + psi2*correction(fvm::ddt(p_rgh))
      + fvc::div(phi, rho2) - fvc::Sp(fvc::div(phi), rho2);

        surfaceScalarField phid1("phid1", fvc::interpolate(psi1)*phi);
        surfaceScalarField phid2("phid2", fvc::interpolate(psi2)*phi);

        p_rghEqnComp1 =
            fvc::ddt(rho1) + fvc::div(phi, rho1) - fvc::Sp(fvc::div(phi), rho1)
          + correction
            (
                psi1*fvm::ddt(p_rgh)
              + fvm::div(phid1, p_rgh) - fvm::Sp(fvc::div(phid1), p_rgh)
            );

[李东岳]

@hongjiewang http://dyfluid.com/rhoPimpleFoam.html 这个你看过没老铁