关于VOF中alphaEqn.H的问题

anubis

         tmp<surfaceScalarField> talphaPhi1Un
         (
             fvc::flux
             (
                 phiCN(),
                 cnCoeff*alpha1 + (1.0 - cnCoeff)*alpha1.oldTime(),
                 alphaScheme
             )
           + fvc::flux
             (
                -fvc::flux(-phir, alpha2, alpharScheme),
                 alpha1,
                 alpharScheme
             )
         );

第二项通量中的-fvc::flux(-phir, alpha2, alpharScheme)为什么要用两个负号？

             MULES::correct
             (
                 geometricOneField(),
                 alpha1,
                 talphaPhi1Un(),
                 talphaPhi1Corr.ref(),
                 Sp,
                 (-Sp*alpha1)(),
                 1,
                 0
             );

调用mules的correct函数的时候，这个talphaPhi1Un()有什么用吗？我看了一下好像correct函数本身没有用到这个参数？

李东岳

flux(-phir, alpha2, alpharScheme)函数插值的时候评估第一个参数的方向，因此-phir和phir是不同的流量方向，因此alpharScheme如果是迎风类或者梯度类格式不同

talphaPhi1Un()是低阶格式的通量，也是用来标记方向

anubis

@东岳
第一个大概明白了，就是用迎风格式的时候
-fvc::flux(-phir, alpha2, alpharScheme) 和 fvc::flux(phir, alpha2, alpharScheme) 会不一样对吧？

第二个的话

tmp<surfaceScalarField> talphaPhi1Corr(talphaPhi1Un() - alphaPhi10);

那alphaPhi10是高阶的吗？我之前以为是talphaPhi1Un是高阶的，alphaPhi10是低阶的。

李东岳

@anubis 哦那我记错了，talphaPhi1Un是高阶的，alphaPhi10是低阶的。