可压缩两相流求解器动量方程的问题

EricLiu

下面为一个名为hybridCentralFoam的可压缩两相流求解器的动量方程UEqn代码，我仔细阅读了一下，发现在方程等号右边有重力和表面张力项，但没有压力梯度项，即类似 fvc::grad(p) 这样的形式。个人觉得奇怪，这样不就无法描述相界面压差对流动的影响了吗？恳求各位大佬指点一二。
求解器方程这块我是刚开始研究的小白，如果我想要加入压力梯度项的话，具体要怎么加入呢，别的地方还需要修改什么吗？

void Foam::vofTwoPhaseCentralFoam::UEqn()
{
    const auto& rho = mixture_model_.rho();

    surfaceScalarField phiU_own
    (
        vF1face_*phi1_own_ +
        vF2face_*phi2_own_
    );
    surfaceScalarField phiU_nei
    (
        vF1face_*phi1_nei_ +
        vF2face_*phi2_nei_
    );
    phiU_own.rename("phiU_own");
    phiU_nei.rename("phiU_nei");

    E_ = fvc::ddt(rho) + fvc::div(phiU_own) + fvc::div(phiU_nei);
    fvVectorMatrix UEqn
    (
          fvm::ddt(rho,U_)
        + fvm::div(phiU_own,U_) + fvm::div(phiU_nei,U_)
        - fvm::Sp(E_,U_)
        + turbulence_->divDevRhoReff(rho,U_)
    );
    oneByA_  = 1.0/UEqn.A();
    HbyA_ = UEqn.H()*oneByA_;
    HbyA_.boundaryFieldRef() == U_.boundaryField();

    B_ = fvc::reconstruct
    (
        (
           interface_.surfaceTensionForce()
            - ghf_*fvc::snGrad(rho)
        ) * B_.mesh().magSf()
    );
}

李东岳

你要看一下这个文章，再去理解那个压力项。openfoam所有的动量方程都没有压力项 http://dyfluid.com/simple.html

EricLiu

@李东岳好的谢谢东岳老师建议，我去研究一下这个文章

CFD中文网

可压缩两相流求解器动量方程的问题