A question about UEqn.H in smallpoolfire2D
 

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!

I'm not an expert in the fire simulation, but from the source which you provided (fireFoam) i understand that term parcels.SU(U) is an explicit contribution to the moment conservation of the fluid from parcels.

This parcels (i guess!!!) are moving in the main fluid and term parcels.SU(U) used to include their drag and other forces.

To obtain more information, you have to make a more deep look at the fireFoam solver source code

Thank you for your reply,now I make a more deep look at the fireFoam solver first!