关于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是低阶的。