From a5aefc54a74d7afbf6ed2976661b07f117f11ba4 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Fri, 12 Jun 2015 09:43:07 +0100
Subject: [PATCH] TurbulenceModels/phaseCompressible: Do not call correctNut in
 constructor if it depends on the availability of the physical properties of
 both phases as the construction of the phase-system may not be complete.

---
 .../phaseCompressible/LES/Niceno/NicenoKEqn.C                 | 4 +++-
 .../phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C | 4 +++-
 .../phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C       | 4 +++-
 3 files changed, 9 insertions(+), 3 deletions(-)

diff --git a/src/TurbulenceModels/phaseCompressible/LES/Niceno/NicenoKEqn.C b/src/TurbulenceModels/phaseCompressible/LES/Niceno/NicenoKEqn.C
index 211e99547..2d3ee8421 100644
--- a/src/TurbulenceModels/phaseCompressible/LES/Niceno/NicenoKEqn.C
+++ b/src/TurbulenceModels/phaseCompressible/LES/Niceno/NicenoKEqn.C
@@ -96,7 +96,9 @@ NicenoKEqn<BasicTurbulenceModel>::NicenoKEqn
 {
     if (type == typeName)
     {
-        correctNut();
+        // Cannot correct nut yet: construction of the phases is not complete
+        // correctNut();
+
         this->printCoeffs(type);
     }
 }
diff --git a/src/TurbulenceModels/phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C b/src/TurbulenceModels/phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C
index 29ba2c732..4eb3187d9 100644
--- a/src/TurbulenceModels/phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C
+++ b/src/TurbulenceModels/phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C
@@ -73,7 +73,9 @@ SmagorinskyZhang<BasicTurbulenceModel>::SmagorinskyZhang
 {
     if (type == typeName)
     {
-        correctNut();
+        // Cannot correct nut yet: construction of the phases is not complete
+        // correctNut();
+
         this->printCoeffs(type);
     }
 }
diff --git a/src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C b/src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
index 7a86dc77e..f523fce5a 100644
--- a/src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
+++ b/src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
@@ -106,7 +106,9 @@ LaheyKEpsilon<BasicTurbulenceModel>::LaheyKEpsilon
 {
     if (type == typeName)
     {
-        correctNut();
+        // Cannot correct nut yet: construction of the phases is not complete
+        // correctNut();
+
         this->printCoeffs(type);
     }
 }