Commit 0696e9aa authored by Andrew Heather's avatar Andrew Heather
Browse files

corrections to chemical enthalpy transfer

parent 5ccff5ff
......@@ -174,9 +174,12 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::cellValueSourceCorrection
const scalar cpc = td.cpInterp().psi()[cellI];
this->cpc_ = (massCell*cpc + addedMass*cpEff)/massCellNew;
const scalar fCarrier = -1.0/td.constProps().hRetentionCoeff();
const scalar fCarrier =
(1.0 - td.cloud().constProps().hRetentionCoeff())
/td.cloud().constProps().hRetentionCoeff();
const scalar dh =
td.cloud().hsTrans()[cellI] + fCarrier*td.cloud().hcTrans()[cellI];
td.cloud().hsTrans()[cellI] - fCarrier*td.cloud().hcTrans()[cellI];
this->Tc_ += dh/(this->cpc_*massCellNew);
}
......
......@@ -39,7 +39,19 @@ constantProperties
(
dimensionedScalar(this->dict().lookup("hRetentionCoeff")).value()
)
{}
{
if ((hRetentionCoeff_ < 0) || (hRetentionCoeff_ > 1))
{
FatalErrorIn
(
"ReactingMultiphaseParcel<ParcelType>::constantProperties::"
"constantProperties"
) << "hRetentionCoeff must be in the range 0 to 1" << nl
<< exit(FatalError) << endl;
}
hRetentionCoeff_ = max(1e-06, hRetentionCoeff_);
}
template<class ParcelType>
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment