Commit ca540320 authored by sergio's avatar sergio
Browse files

BUG: Correcting htc definition to avoid negative values

parent cbcc8754
......@@ -283,9 +283,18 @@ void Foam::externalCoupledTemperatureMixedFvPatchScalarField::writeData
const scalarField Tfluid(tfluid);
// Heat transfer coefficient [W/m2/K]
const scalarField htc(qDot/(max(Twall - Tfluid), 1e-3));
// This htc needs to be always larger or equal to zero
//const scalarField htc(qDot/max(Twall - Tfluid, 1e-3));
scalarField htc(qDot.size(), 0);
forAll (qDot, i)
{
scalar deltaT = mag(Twall[i] - Tfluid[i]);
if (deltaT > 1e-3)
{
htc[i] = sign(qDot[i])*qDot[i]/deltaT;
}
}
const Field<scalar>& magSf = this->patch().magSf();
......
Markdown is supported
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