songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bugfix: added missing multiply by nParticle

Browse Source
This commit is contained in:
andy 2009-04-16 19:37:52 +01:00
parent ef39baacfc
commit 3f4dae4695
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
View File

@ -453,7 +453,7 @@ Foam::scalar Foam::ReactingMultiphaseParcel<ParcelType>::calcDevolatilisation
dMassDV[i] = YVolatile[i]*dMassTot;
}
td.cloud().addToMassDevolatilisation(dMassTot);
td.cloud().addToMassDevolatilisation(this->nParticle_*dMassTot);
return -dMassTot*td.constProps().LDevol();
}
@ -515,7 +515,8 @@ Foam::scalar Foam::ReactingMultiphaseParcel<ParcelType>::calcSurfaceReactions
td.cloud().addToMassSurfaceReaction
(
sum(dMassSRGas) + sum(dMassSRLiquid) + sum(dMassSRSolid)
this->nParticle_
*(sum(dMassSRGas) + sum(dMassSRLiquid) + sum(dMassSRSolid))
);
// Add enthalpy of consumed components to the carrier phase enthalpy

2
src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
View File

@ -274,7 +274,7 @@ Foam::scalar Foam::ReactingParcel<ParcelType>::calcPhaseChange
scalar dMassTot = sum(dMassPC);
// Add to cumulative phase change mass
td.cloud().addToMassPhaseChange(dMassTot);
td.cloud().addToMassPhaseChange(this->nParticle_*dMassTot);
// Effective latent heat of vaporisation
scalar LEff = td.cloud().composition().L(idPhase, YComponents, pc_, T);