TurbulenceModels: Correct nut during construction for single-phase solvers only.

For multiphase solvers the phase construction is not complete at this point.

src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LamBremhorstKE/LamBremhorstKE.C

@@ -180,8 +180,15 @@ LamBremhorstKE::LamBremhorstKE
 
     if (type == typeName)
     {
-        correctNut();
         printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LienCubicKE/LienCubicKE.C

@@ -330,8 +330,15 @@ LienCubicKE::LienCubicKE
 
     if (type == typeName)
     {
-        correctNut();
         printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LienLeschziner/LienLeschziner.C

@@ -224,8 +224,15 @@ LienLeschziner::LienLeschziner
 
     if (type == typeName)
     {
-        correctNut();
         printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/ShihQuadraticKE/ShihQuadraticKE.C

@@ -222,8 +222,15 @@ ShihQuadraticKE::ShihQuadraticKE
 
     if (type == typeName)
     {
-        correctNut();
         printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C

@@ -207,8 +207,15 @@ qZeta::qZeta
 
     if (type == typeName)
     {
-        correctNut();
         printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C

@@ -106,10 +106,10 @@ LaheyKEpsilon<BasicTurbulenceModel>::LaheyKEpsilon
 {
     if (type == typeName)
     {
+        this->printCoeffs(type);
+
         // Cannot correct nut yet: construction of the phases is not complete
         // correctNut();
-
-        this->printCoeffs(type);
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/DeardorffDiffStress/DeardorffDiffStress.C

@@ -110,9 +110,16 @@ DeardorffDiffStress<BasicTurbulenceModel>::DeardorffDiffStress
 {
     if (type == typeName)
     {
-        this->boundNormalStress(this->R_);
-        correctNut();
         this->printCoeffs(type);
+        this->boundNormalStress(this->R_);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/Smagorinsky/Smagorinsky.C

@@ -113,8 +113,15 @@ Smagorinsky<BasicTurbulenceModel>::Smagorinsky
 {
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/SpalartAllmarasDES/SpalartAllmarasDES.C

@@ -318,8 +318,15 @@ SpalartAllmarasDES<BasicTurbulenceModel>::SpalartAllmarasDES
 {
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/WALE/WALE.C

@@ -143,8 +143,15 @@ WALE<BasicTurbulenceModel>::WALE
 {
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/dynamicKEqn/dynamicKEqn.C

@@ -185,8 +185,15 @@ dynamicKEqn<BasicTurbulenceModel>::dynamicKEqn
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C

@@ -122,8 +122,15 @@ dynamicLagrangian<BasicTurbulenceModel>::dynamicLagrangian
 {
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/LES/kEqn/kEqn.C

@@ -113,8 +113,15 @@ kEqn<BasicTurbulenceModel>::kEqn
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/LRR/LRR.C

@@ -200,11 +200,19 @@ LRR<BasicTurbulenceModel>::LRR
 {
     if (type == typeName)
     {
+        this->printCoeffs(type);
+
         this->boundNormalStress(this->R_);
         bound(epsilon_, this->epsilonMin_);
         k_ = 0.5*tr(this->R_);
-        correctNut();
+
-        this->printCoeffs(type);
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/LaunderSharmaKE/LaunderSharmaKE.C

@@ -204,8 +204,15 @@ LaunderSharmaKE<BasicTurbulenceModel>::LaunderSharmaKE
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/RNGkEpsilon/RNGkEpsilon.C

@@ -205,8 +205,15 @@ RNGkEpsilon<BasicTurbulenceModel>::RNGkEpsilon
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/SSG/SSG.C

@@ -209,11 +209,19 @@ SSG<BasicTurbulenceModel>::SSG
 {
     if (type == typeName)
     {
+        this->printCoeffs(type);
+
         this->boundNormalStress(this->R_);
         bound(epsilon_, this->epsilonMin_);
         k_ = 0.5*tr(this->R_);
-        correctNut();
+
-        this->printCoeffs(type);
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/SpalartAllmaras/SpalartAllmaras.C

@@ -255,8 +255,15 @@ SpalartAllmaras<BasicTurbulenceModel>::SpalartAllmaras
 {
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C

@@ -188,7 +188,14 @@ kEpsilon<BasicTurbulenceModel>::kEpsilon
     if (type == typeName)
     {
         this->printCoeffs(type);
-        correctNut();
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/kOmega/kOmega.C

@@ -147,8 +147,15 @@ kOmega<BasicTurbulenceModel>::kOmega
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/kOmegaSST/kOmegaSST.C

@@ -350,8 +350,15 @@ kOmegaSST<BasicTurbulenceModel>::kOmegaSST
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/realizableKE/realizableKE.C

@@ -226,8 +226,15 @@ realizableKE<BasicTurbulenceModel>::realizableKE
 
     if (type == typeName)
     {
-        correctNut();
         this->printCoeffs(type);
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }
 

src/TurbulenceModels/turbulenceModels/RAS/v2f/v2f.C

@@ -239,7 +239,14 @@ v2f<BasicTurbulenceModel>::v2f
     if (type == typeName)
     {
         this->printCoeffs(type);
-        correctNut();
+
+        // Correct nut for single-phase solvers only.
+        // For multiphase solvers the phase construction is not complete
+        // at this point.
+        if (isType<geometricOneField>(alpha))
+        {
+            correctNut();
+        }
     }
 }