CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Please explain the implementation of species transport Eqn in reactingFoam (http://www.cfd-online.com/Forums/openfoam-solving/58475-please-explain-implementation-species-transport-eqn-reactingfoam.html)

kallipygian October 13, 2008 07:29

I had been working with the bu
 
I had been working with the buoyantFoam solver. I have succeeded in implementing buoyantFoam solver with SIMPLE Algorithm. Now I want to add species transport equation to the buoyant Solver.

On way, I stumbled upon the rectingFoam solver where a generic species transport equation is implemented in $FOAM_SOLVER/combustion/reactingFoam/YEqn.H

Can someone please explain how this implementation takes care of change in enthalpy due to species transport?

My first impression was that the implementation was faulty, but I later realized that the in the createField.H file, the implementation of species also include the enthalpy of the gas.

#################################################
combustionMixture& composition = thermo->composition();
PtrList<volscalarfield>& Y = composition.Y();

word inertSpecie(thermo->lookup("inertSpecie"));
#################################################

Further, in the YEqn.H a new convectionScheme is inherited from the "fv" class called mvConvection, which defines Yi_h. So, somehow the enthalpy is included in the implementation but is very elusive to my inexperienced eyes.

Can someone please help me understand it better?

(**The implementation looks very native to hCombustionThermo.)


All times are GMT -4. The time now is 16:24.