From 08a9c9da8f9f27e6fc4f0d5b4ad2fc3499c36de8 Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Tue, 7 Sep 2010 09:38:11 +0100
Subject: [PATCH] BUG: ID 0000028 Inconsistency in nuEffa declaration
---
.../solvers/multiphase/twoPhaseEulerFoam/UEqns.H | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
index 047eeac7e2f..c32855904a8 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@@ -4,11 +4,6 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
{
{
volTensorField gradUaT = fvc::grad(Ua)().T();
- volTensorField Rca
- (
- "Rca",
- ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
- );
if (kineticTheory.on())
{
@@ -20,6 +15,13 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
nuEffa = sqr(Ct)*nutb + nua;
}
+ volTensorField Rca
+ (
+ "Rca",
+ ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
+ );
+
+
if (kineticTheory.on())
{
Rca -= ((kineticTheory.lambda()/rhoa)*tr(gradUaT))*tensor(I);
--
GitLab