33 lines
623 B
C
33 lines
623 B
C
fvVectorMatrix UcEqn
|
|
(
|
|
fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
|
|
- fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
|
|
+ continuousPhaseTurbulence->divDevRhoReff(Uc)
|
|
==
|
|
(1.0/rhoc)*cloudSU
|
|
);
|
|
|
|
UcEqn.relax();
|
|
|
|
volScalarField rAUc(1.0/UcEqn.A());
|
|
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));
|
|
|
|
surfaceScalarField phicForces
|
|
(
|
|
(fvc::interpolate(rAUc*cloudVolSUSu/rhoc) & mesh.Sf())
|
|
+ rAUcf*(g & mesh.Sf())
|
|
);
|
|
|
|
if (pimple.momentumPredictor())
|
|
{
|
|
solve
|
|
(
|
|
UcEqn
|
|
==
|
|
fvc::reconstruct
|
|
(
|
|
phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
|
|
)
|
|
);
|
|
}
|