Source Term calculation for multiphase reactive flows [CoalChemistryFoam]
 

Hello Foam Users,

So I am starting to understand how and where the source terms for coal particle combustion is calculated. I found out that in the solver coalChemistryFoam the UEqn.H and hEqn.H is as follows:

UEqn.H

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
rho.dimensionedInternalField()*g
+ coalParcels.SU()
+ limestoneParcels.SU()
);
UEqn.relax();
if (momentumPredictor)
{
solve(UEqn == -fvc::grad(p));

hEqn.H

fvScalarMatrix hEqn
(
fvm::ddt(rho, h)
+ mvConvection->fvmDiv(phi, h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
DpDt
+ coalParcels.Sh()
+ limestoneParcels.Sh()
+ enthalpySource.Su()
+ radiation->Sh(thermo)
);
hEqn.relax();
hEqn.solve();
thermo.correct();
radiation->correct();

I understand that the source terms from particle momentum is given to CFD cell as coalParcels.SU() and from particle energy equation is given to CFD cell as coalParcels.Sh().

My question is:

1. Where are these source terms calculated.
2. has any one tried to dsitribute these source terms over surrounding cells/nodes

Thanks,
Nir