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 .