From 840d43537c9f3c747715552c91f6640d8aa92ef6 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Fri, 5 Apr 2019 09:20:09 -0700
Subject: [PATCH] ENH: Bounding htc for qDot close to zero

---
 .../externalCoupledTemperatureMixedFvPatchScalarField.C         | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C b/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
index 5a06b0dc361..f1d78d5290d 100644
--- a/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
+++ b/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
@@ -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();
 
-- 
GitLab