void Foam::ThermoParcel<ParcelType>::cellValueSourceCorrection ( TrackData& td, const scalar dt, const label cellI ) { this->Uc_ += td.cloud().UTrans()[cellI]/this->massCell(cellI); const scalar CpMean = td.CpInterp().psi()[cellI]; Tc_ += td.cloud().hsTrans()[cellI]/(CpMean*this->massCell(cellI)); if (Tc_ < td.cloud().constProps().TMin()) { if (debug) { WarningIn ( "void Foam::ThermoParcel<ParcelType>::cellValueSourceCorrection" "(" "TrackData&, " "const scalar, " "const label" ")" ) << "Limiting observed temperature in cell " << cellI << " to " << td.cloud().constProps().TMin() << nl << endl; } Tc_ = td.cloud().constProps().TMin(); } }
void Foam::ThermoParcel<ParcelType>::setCellValues ( TrackData& td, const scalar dt, const label cellI ) { ParcelType::setCellValues(td, dt, cellI); tetIndices tetIs = this->currentTetIndices(); Cpc_ = td.CpInterp().interpolate(this->position(), tetIs); Tc_ = td.TInterp().interpolate(this->position(), tetIs); if (Tc_ < td.cloud().constProps().TMin()) { if (debug) { WarningIn ( "void Foam::ThermoParcel<ParcelType>::setCellValues" "(" "TrackData&, " "const scalar, " "const label" ")" ) << "Limiting observed temperature in cell " << cellI << " to " << td.cloud().constProps().TMin() << nl << endl; } Tc_ = td.cloud().constProps().TMin(); } }