From 93587afefc2ec98bb8f25c505e0d86f589fee9d2 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 15 Nov 2017 12:58:58 -0800
Subject: [PATCH] ENH: Avoidding 0/0 conflict in
 externalCoupledTemperatureMixed increaging the zero order of denomitaror

---
 .../externalCoupledTemperatureMixedFvPatchScalarField.C       | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C b/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
index 70128b8800..8961922b4b 100644
--- a/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
+++ b/src/functionObjects/field/externalCoupled/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C
@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2013-2016 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2017 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -196,7 +196,7 @@ void Foam::externalCoupledTemperatureMixedFvPatchScalarField::writeData
     const scalarField Tc(patchInternalField());
 
     // Heat transfer coefficient [W/m2/K]
-    const scalarField htc(qDot/(Tp - Tc + ROOTVSMALL));
+    const scalarField htc(qDot/(Tp - Tc + 1e-3));
 
     const Field<scalar>& magSf(this->patch().magSf());