updating turb coeffs + fix typo

applications/solvers/multiphase/settlingFoam/createFields.H

@@ -56,20 +56,11 @@
     );
 
 
-    dimensionedScalar rhoc
+    dimensionedScalar rhoc(transportProperties.lookup("rhoc"));
-    (
-        transportProperties.lookup("rhoc")
-    );
 
-    dimensionedScalar rhod
+    dimensionedScalar rhod(transportProperties.lookup("rhod"));
-    (
-        transportProperties.lookup("rhod")
-    );
 
-    dimensionedScalar muc
+    dimensionedScalar muc(transportProperties.lookup("muc"));
-    (
-        transportProperties.lookup("muc")
-    );
 
     dimensionedScalar plasticViscosityCoeff
     (
@@ -96,10 +87,7 @@
         transportProperties.lookup("yieldStressOffset")
     );
 
-    Switch BinghamPlastic
+    Switch BinghamPlastic(transportProperties.lookup("BinghamPlastic"));
-    (
-        transportProperties.lookup("BinghamPlastic")
-    );
 
     volScalarField rho
     (
@@ -147,8 +135,8 @@
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
-        muc +
+        muc
-        plasticViscosity
+      + plasticViscosity
         (
             plasticViscosityCoeff,
             plasticViscosityExponent,
@@ -176,10 +164,7 @@
 
     Info<< "Selecting Drift-Flux model " << endl;
 
-    word VdjModel
+    word VdjModel(transportProperties.lookup("VdjModel"));
-    (
-        transportProperties.lookup("VdjModel")
-    );
 
     Info<< tab << VdjModel << " selected\n" << endl;
 
@@ -188,26 +173,13 @@
         transportProperties.subDict(VdjModel + "Coeffs")
     );
 
-    dimensionedVector V0
+    dimensionedVector V0(VdjModelCoeffs.lookup("V0"));
-    (
-        VdjModelCoeffs.lookup("V0")
-    );
 
-    dimensionedScalar a
+    dimensionedScalar a(VdjModelCoeffs.lookup("a"));
-    (
-        VdjModelCoeffs.lookup("a")
-    );
 
-    dimensionedScalar a1
+    dimensionedScalar a1(VdjModelCoeffs.lookup("a1"));
-    (
-        VdjModelCoeffs.lookup("a1")
-    );
-
-    dimensionedScalar alphaMin
-    (
-        VdjModelCoeffs.lookup("alphaMin")
-    );
 
+    dimensionedScalar alphaMin(VdjModelCoeffs.lookup("alphaMin"));
 
 
     IOdictionary RASProperties
@@ -223,15 +195,9 @@
     );
 
 
-    Switch turbulence
+    Switch turbulence(RASProperties.lookup("turbulence"));
-    (
-        RASProperties.lookup("turbulence")
-    );
 
-    dictionary kEpsilonDict
+    dictionary kEpsilonDict(RASProperties.subDictPtr("kEpsilonCoeffs"));
-    (
-        RASProperties.subDictPtr("kEpsilonCoeffs")
-    );
 
     dimensionedScalar Cmu
     (
@@ -273,30 +239,27 @@
         )
     );
 
-    dimensionedScalar alphak
+    dimensionedScalar sigmak
     (
         dimensionedScalar::lookupOrAddToDict
         (
-            "alphaEps",
+            "sigmak",
             kEpsilonDict,
             1.0
         )
     );
 
-    dimensionedScalar alphaEps
+    dimensionedScalar sigmaEps
     (
         dimensionedScalar::lookupOrAddToDict
         (
-            "alphaEps",
+            "sigmaEps",
             kEpsilonDict,
-            0.76923
+            1.3
         )
     );
 
-    dictionary wallFunctionDict
+    dictionary wallFunctionDict(RASProperties.subDictPtr("wallFunctionCoeffs"));
-    (
-        RASProperties.subDictPtr("wallFunctionCoeffs")
-    );
 
     dimensionedScalar kappa
     (

applications/solvers/multiphase/settlingFoam/kEpsilon.H

@@ -15,7 +15,7 @@ if(turbulence)
     tgradU.clear();
 
     volScalarField Gcoef =
-        alphak*Cmu*k*(g & fvc::grad(rho))/(epsilon + epsilon0);
+        Cmu*k/sigmak*(g & fvc::grad(rho))/(epsilon + epsilon0);
 
     #include "wallFunctions.H"
 
@@ -26,7 +26,7 @@ if(turbulence)
       + fvm::div(phi, epsilon)
       - fvm::laplacian
         (
-            alphaEps*mut + mul, epsilon,
+            mut/sigmaEps + mul, epsilon,
             "laplacian(DepsilonEff,epsilon)"
         )
      ==
@@ -51,7 +51,7 @@ if(turbulence)
       + fvm::div(phi, k)
       - fvm::laplacian
         (
-            alphak*mut + mul, k,
+            mut/sigmak + mul, k,
             "laplacian(DkEff,k)"
         )
      ==