ENH: Bounding htc for qDot close to zero

This commit is contained in:
sergio 2019-04-05 09:20:09 -07:00 committed by Andrew Heather
parent 49ba6a8e94
commit 968b95d0cf

View File

@ -281,7 +281,7 @@ void Foam::externalCoupledTemperatureMixedFvPatchScalarField::writeData
// Heat transfer coefficient [W/m2/K]
const scalarField htc(qDot/(Twall - Tfluid + 1e-3));
const scalarField htc(qDot/(max(Twall - Tfluid), 1e-3));
const Field<scalar>& magSf = this->patch().magSf();