From a01b5bb4a064656d07825359784fa34424c2b894 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 29 Dec 2014 11:06:08 +0000
Subject: [PATCH] Corrected/improved mixed BC formulation to ensure the
 reference temperature is physical Resolves bug-report
 http://www.openfoam.org/mantisbt/view.php?id=1369

---
 ...peratureRadCoupledMixedFvPatchScalarField.C | 18 ++++++++----------
 1 file changed, 8 insertions(+), 10 deletions(-)

diff --git a/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C b/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C
index 1968da38..19f8a971 100644
--- a/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C
+++ b/src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C
@@ -219,7 +219,7 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::updateCoeffs()
     }
     mpp.distribute(KDeltaNbr);
 
-    scalarField KDelta(kappa(*this)*patch().deltaCoeffs());
+    scalarField KDelta(kappa(Tp)*patch().deltaCoeffs());
 
     scalarField Qr(Tp.size(), 0.0);
     if (QrName_ != "none")
@@ -234,17 +234,15 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::updateCoeffs()
         mpp.distribute(QrNbr);
     }
 
-    scalarField alpha(KDeltaNbr - (Qr + QrNbr)/Tp);
-
-    valueFraction() = alpha/(alpha + KDelta);
-
-    refValue() = (KDeltaNbr*TcNbr)/alpha;
+    valueFraction() = KDeltaNbr/(KDeltaNbr + KDelta);
+    refValue() = TcNbr;
+    refGrad() = (Qr + QrNbr)/kappa(Tp);
 
     mixedFvPatchScalarField::updateCoeffs();
 
     if (debug)
     {
-        scalar Q = gSum(kappa(*this)*patch().magSf()*snGrad());
+        scalar Q = gSum(kappa(Tp)*patch().magSf()*snGrad());
 
         Info<< patch().boundaryMesh().mesh().name() << ':'
             << patch().name() << ':'
@@ -254,9 +252,9 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::updateCoeffs()
             << this->dimensionedInternalField().name() << " :"
             << " heat transfer rate:" << Q
             << " walltemperature "
-            << " min:" << gMin(*this)
-            << " max:" << gMax(*this)
-            << " avg:" << gAverage(*this)
+            << " min:" << gMin(Tp)
+            << " max:" << gMax(Tp)
+            << " avg:" << gAverage(Tp)
             << endl;
     }