关于VOF中alphaEqn.H的问题



  •          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()是低阶格式的通量,也是用来标记方向



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

    第二个的话

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

    那alphaPhi10是高阶的吗?我之前以为是talphaPhi1Un是高阶的,alphaPhi10是低阶的。



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2