CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Source Term calculation for multiphase reactive flows [CoalChemistryFoam] (http://www.cfd-online.com/Forums/openfoam/78018-source-term-calculation-multiphase-reactive-flows-coalchemistryfoam.html)

 N. A. July 9, 2010 13:37

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)
{

hEqn.H

fvScalarMatrix hEqn
(
fvm::ddt(rho, h)
+ mvConvection->fvmDiv(phi, h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
DpDt
+ coalParcels.Sh()
+ limestoneParcels.Sh()
+ enthalpySource.Su()
);
hEqn.relax();
hEqn.solve();
thermo.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

 All times are GMT -4. The time now is 14:28.