在pEqn.H中对U进行更新时有
const surfaceScalarField phig(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf());
U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf);
我想请问其中的p_rghEqn.flux()代表什么意思,以及为什么phig中要先乘rhorAUf再在U中除以rhorAUf;
以及之后的根据密度差判断收敛性中,我发现compressibleContinuityErrs.H代码如下,只计算了误差值,但没有判断是否符合要求,想请问各位大佬是在哪里判断的呢
dimensionedScalar totalMass = fvc::domainIntegrate(rho);
scalar sumLocalContErr =
(fvc::domainIntegrate(mag(rho - thermo.rho()))/totalMass).value();
scalar globalContErr =
(fvc::domainIntegrate(rho - thermo.rho())/totalMass).value();
cumulativeContErr += globalContErr;
Info<< "time step continuity errors : sum local = " << sumLocalContErr
<< ", global = " << globalContErr
<< ", cumulative = " << cumulativeContErr
<< endl;
}
谢谢大家!