Please confirm, if this is a bug or a typo
 

Dear Foamers,

Can someone confirm the model of surface reaction using lagrangian particles. I am simulating cola-particle consumption. The consumpotio of coal is as follows:

C + O2 --> CO2

dmC, cMO2 and dmCO2 are being calculated to find the amount of species consumed. Later the particle mass is updated and so is the gas phase species mass.

I am wondering if dMassSolid[CslocalId_] +- dmC; (instead of -dmC)

The part of file, COxidationDiffusionLimitedRate, is as below:
...
...
// Local mass fraction of O2 in the carrier phase
const scalar YO2 = this->owner().mcCarrierThermo().Y(O2GlobalId_)[cellI];
// Change in C mass [kg]
scalar dmC =
4.0*mathematicalConstant::pi*d*D_*YO2*Tc*rhoc/(Sb_*(T + Tc))*dt;
// Limit mass transfer by availability of C
dmC = min(mass*fComb, dmC);
// Change in O2 mass [kg]
const scalar dmO2 = dmC/WC_*Sb_*WO2_;
// Mass of newly created CO2 [kg]
const scalar dmCO2 = dmC + dmO2;
// Update local particle C mass
dMassSolid[CsLocalId_] += dmC;
// Update carrier O2 and CO2 mass
dMassSRCarrier[O2GlobalId_] -= dmO2;
dMassSRCarrier[CO2GlobalId_] += dmCO2;
// Heat of reaction [J]
return -HcCO2_*dmCO2;
}


-------------

Should I change dMassSolid[CsLocalId_] += dmC; to dMassSolid[CsLocalId_] -= dmC; ?

Thanks,
N