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/)
-   -   question about hcTrans calculation in ReactingParcel.C (http://www.cfd-online.com/Forums/openfoam-solving/68405-question-about-hctrans-calculation-reactingparcel-c.html)

wang.zhy September 18, 2009 05:39

question about hcTrans calculation in ReactingParcel.C
 
Hi Foamers
I have got a confused question concerning the calculation of hcTrans in ReactingParcel.C .
// Transfer mass lost from particle to carrier mass source
forAll(dMassPC, i)
{
label gid = td.cloud().composition().localToGlobalCarrierId(0, i);
td.cloud().rhoTrans(gid)[cellI] += np0*dMassPC[i];
td.cloud().hcTrans()[cellI] +=
np0
*dMassPC[i]
*td.cloud().mcCarrierThermo().speciesData()[gid].H(T0);
}
As the comment said ,dMassPC is the mass transfer from particle to carrier.
If dMassPC>0,carrier get mass transfered from particle.
so for mass source: rhoTrans = np0*dMassPC
but ,for hcTrans ,i think it might be different.
in the process that particle give mass to carrier,
carrier shoud give latent heat to particle .
So ,
td.cloud().hcTrans()[cellI] +=
np0
*dMassPC[i]
*td.cloud().mcCarrierThermo().speciesData()[gid].H(T0);

shoud be modified to reflect the Latent heat transport .
this is my modified version
td.cloud().hcTrans()[cellI] -=np0*dMassPC[i]*L;
am i right ?
its important for me !want help!
sorry for bad english .


All times are GMT -4. The time now is 13:36.