From ef49aaf03d09d4d88693ae827371c396598cb2de Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 28 Feb 2012 15:54:02 +0000
Subject: [PATCH] twoPhaseEulerFoam: Remove k from the Reynolds stress because
 there is only a single pressure in the system and to avoid the generation of
 flow from the gradient of k the pressure gradient must balance the k term in
 each U equation

---
 applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
index 0a590e6bf9..a4d428e176 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@@ -18,7 +18,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
         volTensorField Rca
         (
             "Rca",
-            ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
+            (((2.0/3.0)*I)*nuEffa)*tr(gradUaT) - nuEffa*gradUaT
         );
 
         if (kineticTheory.on())
@@ -62,7 +62,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
         volTensorField Rcb
         (
             "Rcb",
-            ((2.0/3.0)*I)*(k + nuEffb*tr(gradUbT)) - nuEffb*gradUbT
+            (((2.0/3.0)*I)*nuEffb)*tr(gradUbT) - nuEffb*gradUbT
         );
 
         surfaceScalarField phiRb