reactingTwoPhaseEulerFoam求解器的相变求解过程

zwl

@东岳今天仔细推导了可压缩两相方程，基本都能与代码都能对的上，pEqnComp代表的就是可压缩项与质量源相的和，但是为什么要返给相方程pEqnComp与p_rgh的积？这跟您之前分析的其它可压缩求解器也不一样？能分析下为什么这么处理吗？公式中没有质量源相，相方程可压缩部分与压力的关系，这样处理的目的和原理是什么？谢谢。

李东岳

为什么要返给相方程pEqnComp与p_rgh的积？

pEqnComp是矩阵系数，与p_rgh的积才是真正的volScalarField

zwl

@东岳那么从物理角度应该怎么理解呢？dgdt*p_rgh什么物理意义？应该不只是体标量场这么简单吧。非常感谢东岳老师耐心的回复。

李东岳

那个是dgdt，不是dgdt*p_rgh哟

zwl

@东岳对，您说的是对的。那我这样表达下我的意思吧：
reactingTwoPhaseEulerFaom的dgdt这里叫dgdt1；compressibleInterFoam里的dgdt这里叫dgdt2；我根据代码解读看到dgdt1=alpha1* alpha2* dgdt2* p_rgh，而且我感觉dgdt1应该定义成alpha1* alpha2* dgdt2才对，膨胀率乘以压力的处理和意义是公式里没有的，所以我才一直有疑惑。不知道我说明白了没有。

zwl

@zwl 在 reactingTwoPhaseEulerFoam求解器的相变求解过程中说：

dgdt1=alpha1alpha2dgdt2p_rgh

乘号漏掉了：dgdt1=alpha1 * alpha2 * dgdt2*p_rgh

李东岳

我根据代码解读看到dgdt1=alpha1* alpha2* dgdt2* p_rgh

是dgdt1=alpha1* alpha2* dgdt2吧，就像我上面说的那个

那个是dgdt，不是dgdt*p_rgh哟

膨胀率乘以压力的处理和意义是公式里没有的

是没有啊，那个是dgdt，不是dgdt乘以p_rgh..

zwl

@东岳代码如下：


    tmp<fvScalarMatrix> pEqnComp1;
    tmp<fvScalarMatrix> pEqnComp2;

    // Construct the compressibility parts of the pressure equation
    if (pimple.transonic())
    {
        if (phase1.compressible())
        {
            surfaceScalarField phid1
            (
                IOobject::groupName("phid", phase1.name()),
                fvc::interpolate(psi1)*phi1
            );

            pEqnComp1 =
                (
                    phase1.continuityError()
                  - fvc::Sp(fvc::ddt(alpha1) + fvc::div(alphaPhi1), rho1)
                )/rho1           //相方程可压缩部分 
              + correction
                (
                    (alpha1/rho1)*
                    (
                        psi1*fvm::ddt(p_rgh)
                      + fvm::div(phid1, p_rgh) - fvm::Sp(fvc::div(phid1), p_rgh)
                    )
                );

            deleteDemandDrivenData(pEqnComp1.ref().faceFluxCorrectionPtr());
            pEqnComp1.ref().relax();
        }

        if (phase1.compressible())
        {
            phase1.divU(-pEqnComp1 & p_rgh);//divU(divU_)函数返回自身参数
        }


    // Construct the dilatation rate source term
    tmp<volScalarField::Internal> tdgdt;

    if (phase1_.divU().valid() && phase2_.divU().valid())
    {
        tdgdt =
        (
            alpha2()
           *phase1_.divU()()()
          - alpha1()
           *phase2_.divU()()()
        );
    }

如此看来是不是dgdt1=alpha1* alpha2* dgdt2* p_rgh？难道我还有遗漏？

李东岳

没理解，就像我之前说的

pEqnComp是矩阵系数，与p_rgh的积才是真正的volScalarField

代码中的dgdt2* p_rgh实际上是公式中的dgdt2，pEqnComp1是矩阵系数并不是值，与p_rgh的积才是dgdt2

zwl

@东岳懂了，谢谢。

CFD中文网

reactingTwoPhaseEulerFoam求解器的相变求解过程