DeardorffDiffStress: updated to correspond more closely to the original paper

Added references to the generalized gradient diffusion model

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

@@ -93,7 +93,7 @@ DeardorffDiffStress<BasicTurbulenceModel>::DeardorffDiffStress
         (
             "Ce",
             this->coeffDict_,
-            1.048
+            1.05
         )
     ),
     Cs_
@@ -102,7 +102,7 @@ DeardorffDiffStress<BasicTurbulenceModel>::DeardorffDiffStress
         (
             "Cs",
             this->coeffDict_,
-            0.22
+            0.25
         )
     )
 {
@@ -189,12 +189,11 @@ void DeardorffDiffStress<BasicTurbulenceModel>::correct()
     (
         fvm::ddt(alpha, rho, R)
       + fvm::div(alphaRhoPhi, R)
-      - fvm::laplacian(Cs_*(k/this->epsilon())*R, R)
+      - fvm::laplacian(I*this->nu() + Cs_*(k/this->epsilon())*R, R)
       + fvm::Sp(Cm_*alpha*rho*sqrt(k)/this->delta(), R)
      ==
         alpha*rho*P
-      + (4.0/5.0)*alpha*rho*k*D // Deardorff
+      + (4.0/5.0)*alpha*rho*k*D
-    //- 0.6*alpha*rho*dev(P) // LRR
       - ((2.0/3.0)*(1.0 - Cm_/this->Ce_)*I)*(alpha*rho*this->epsilon())
     );
 
@@ -203,7 +202,6 @@ void DeardorffDiffStress<BasicTurbulenceModel>::correct()
 
     this->boundNormalStress(this->R_);
     correctNut();
-    this->correctWallShearStress(this->R_);
 }
 
 

src/TurbulenceModels/turbulenceModels/LES/DeardorffDiffStress/DeardorffDiffStress.H

@@ -42,6 +42,21 @@ Description
     This SGS model uses a full balance equation for the SGS stress tensor to
     simulate the behaviour of B.
 
+    This implementation is as described in the above paper except that the
+    triple correlation model of Donaldson is replaced with the generalized
+    gradient diffusion model of Daly and Harlow:
+    \verbatim
+        Daly, B. J., & Harlow, F. H. (1970).
+        Transport equations in turbulence.
+        Physics of Fluids (1958-1988), 13(11), 2634-2649.
+    \endverbatim
+    with the default value for the coefficient Cs of 0.25 from
+    \verbatim
+        Launder, B. E., Reece, G. J., & Rodi, W. (1975).
+        Progress in the development of a Reynolds-stress turbulence closure.
+        Journal of fluid mechanics, 68(03), 537-566.
+    \endverbatim
+
 SourceFiles
     DeardorffDiffStress.C
 
@@ -134,15 +149,6 @@ public:
         //- Return the turbulence kinetic energy dissipation rate
         virtual tmp<volScalarField> epsilon() const;
 
-        //- Return the effective diffusivity for B
-        tmp<volScalarField> DBEff() const
-        {
-            return tmp<volScalarField>
-            (
-                new volScalarField("DBEff", this->nut_ + this->nu())
-            );
-        }
-
         //- Correct sub-grid stress, eddy-Viscosity and related properties
         virtual void correct();
 };

src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C

@@ -165,7 +165,17 @@ Foam::ReynoldsStress<BasicTurbulenceModel>::ReynoldsStress
         ),
         this->mesh_
     )
-{}
+{
+    if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
+    {
+        FatalErrorIn
+        (
+            "ReynoldsStress::ReynoldsStress"
+        )   << "couplingFactor = " << couplingFactor_
+            << " is not in range 0 - 1" << nl
+            << exit(FatalError);
+    }
+}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //