From f9c99a109cb9ae333a9688d3e2baf11e9068eafb Mon Sep 17 00:00:00 2001
From: Robert Keser <keser@son>
Date: Wed, 29 Mar 2017 16:49:49 +0200
Subject: [PATCH] Updated cross-coupling terms in coupledKEpsilon turbulence
 model.

---
 .../RAS/coupledKEpsilon/coupledKEpsilon.C           | 13 +++++++------
 1 file changed, 7 insertions(+), 6 deletions(-)

diff --git a/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C b/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C
index 2b7b738ec..6cc5ab1fa 100644
--- a/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C
+++ b/src/turbulenceModels/incompressible/RAS/coupledKEpsilon/coupledKEpsilon.C
@@ -253,9 +253,9 @@ void coupledKEpsilon::correct()
           + fvm::div(phi_, epsilon_)
           + fvm::SuSp(-fvc::div(phi_), epsilon_)
           - fvm::laplacian(DepsilonEff(), epsilon_)
-          + fvm::Sp(2*C2_*epsilon_/k_, epsilon_)
          ==
-            2*C1_*Cmu_*magSqr(symm(fvc::grad(U_)))*k_
+            C1_*G*epsilon_/k_
+          - fvm::Sp(2*C2_*epsilon_/k_, epsilon_)
           + C2_*sqr(epsilon_)/k_
         );
 
@@ -268,8 +268,7 @@ void coupledKEpsilon::correct()
         volScalarField coupling
         (
             "coupling",
-            -2*C1_*Cmu_*magSqr(symm(fvc::grad(U_)))
-           - C2_*sqr(epsilon_/k_)
+            -C2_*sqr(epsilon_/k_)
         );
         scalarField& couplingIn = coupling.internalField();
 
@@ -295,8 +294,10 @@ void coupledKEpsilon::correct()
           + fvm::div(phi_, k_)
           + fvm::SuSp(-fvc::div(phi_), k_)
           - fvm::laplacian(DkEff(), k_)
-          + fvm::Sp(Cmu_*k_/(nut_ + nutSmall), k_)
-          - G
+         ==
+            G
+          + Cmu_*sqr(k_)/(nut_+nutSmall)
+          - fvm::Sp(2*Cmu_*k_/(nut_+nutSmall), k_)
         );
 
         kEqn.relax();