最近在看DPMFoam的代码，但是发现动量方程好像和其数学表达式对不上
fvVectorMatrix UcEqn
(
fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
- fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
+ continuousPhaseTurbulence->divDevTau(Uc)
==
(1.0/rhoc)*cloudSU
);
UcEqn.relax();
fvOptions.constrain(UcEqn);
volScalarField rAUc(1.0/UcEqn.A());
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));
surfaceScalarField phicForces
(
fvc::flux(rAUc*cloudVolSUSu/rhoc) + rAUcf*(g & mesh.Sf())
);
if (pimple.momentumPredictor())
{
solve
(
UcEqn
==
fvc::reconstruct
(
phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
)
);
fvOptions.correct(Uc);
}
此处通过fvc::reconstrct返回的
fvc::snGrad(p)*mesh.magSf()
是否需要乘上alpha
