From a98fbcd493d1fc4ab000cb4728ff5de047a3afda Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Wed, 14 Oct 2020 13:29:00 -0700
Subject: [PATCH] BUG: Fix to powerLaw model to use nu0 from thermo

---
 .../generalizedNewtonianViscosityModels/powerLaw/powerLaw.C   | 4 +---
 .../generalizedNewtonianViscosityModels/powerLaw/powerLaw.H   | 1 -
 2 files changed, 1 insertion(+), 4 deletions(-)

diff --git a/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.C b/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.C
index 833908189b..bb925aa068 100644
--- a/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.C
+++ b/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.C
@@ -60,7 +60,6 @@ Foam::laminarModels::generalizedNewtonianViscosityModels::powerLaw::powerLaw
 :
     generalizedNewtonianViscosityModel(viscosityProperties),
     powerLawCoeffs_(viscosityProperties.optionalSubDict(typeName + "Coeffs")),
-    k_("k", dimViscosity, powerLawCoeffs_),
     n_("n", dimless, powerLawCoeffs_),
     nuMin_("nuMin", dimViscosity, powerLawCoeffs_),
     nuMax_("nuMax", dimViscosity, powerLawCoeffs_)
@@ -80,7 +79,6 @@ bool Foam::laminarModels::generalizedNewtonianViscosityModels::powerLaw::read
 
     powerLawCoeffs_ = viscosityProperties.optionalSubDict(typeName + "Coeffs");
 
-    powerLawCoeffs_.readEntry("k", k_);
     powerLawCoeffs_.readEntry("n", n_);
     powerLawCoeffs_.readEntry("nuMin", nuMin_);
     powerLawCoeffs_.readEntry("nuMax", nuMax_);
@@ -103,7 +101,7 @@ nu
         min
         (
             nuMax_,
-            k_*pow
+            nu0*pow
             (
                 max
                 (
diff --git a/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.H b/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.H
index 2b441de2d8..62dc0dc29f 100644
--- a/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.H
+++ b/src/TurbulenceModels/turbulenceModels/laminar/generalizedNewtonian/generalizedNewtonianViscosityModels/powerLaw/powerLaw.H
@@ -62,7 +62,6 @@ class powerLaw
 
         dictionary powerLawCoeffs_;
 
-        dimensionedScalar k_;
         dimensionedScalar n_;
         dimensionedScalar nuMin_;
         dimensionedScalar nuMax_;