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