[zqlhzx]

Hi Foamers:
According to the following codes in UEqn.H:
fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
parcels.SU(U)
);
UEqn.relax();
if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
)*mesh.magSf()
)
);
}
The momentum equation that the smallpoolFire2D sloves is
fvm::ddt(rho, U)+ fvm::div(phi, U)+ turbulence->divDevRhoReff(U)
==parcels.SU(U)+fvc::reconstruct((- ghf*fvc::snGrad(rho)- fvc::snGrad(p_rgh))*mesh.magSf()).
I understand all term in the equation except the term of "parcels.SU(U)".
In many example in OF,there is no term of "parcels.SU(U)" in UEqn.H.How can understand the term?Is there any theory about the term?
please some foamers know tell me.I really need help!
Thank you!