From 448cd8b97b5989525a0e4f659fc8cfbe99bde00d Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Sat, 12 Mar 2011 15:49:38 +0000
Subject: [PATCH] MULES: Corrected the handling of density variation for the
 static mesh case

---
 .../fvMatrices/solvers/MULES/MULESTemplates.C | 19 ++++++++++++++-----
 1 file changed, 14 insertions(+), 5 deletions(-)

diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
index d7a92be388..ea7de07931 100644
--- a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
+++ b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
@@ -457,12 +457,20 @@ void Foam::MULES::limiter
         tmp<volScalarField::DimensionedInternalField> V0 = mesh.Vsc0();
 
         psiMaxn =
-            V*((rho.field()/deltaT - Sp.field())*psiMaxn - Su.field())
+            V
+           *(
+               (rho.field()/deltaT - Sp.field())*psiMaxn
+             - Su.field()
+            )
           - (V0().field()/deltaT)*rho.oldTime().field()*psi0
           + sumPhiBD;
 
         psiMinn =
-            V*(Su.field() - (rho.field()/deltaT - Sp.field())*psiMinn)
+            V
+           *(
+               Su.field()
+             - (rho.field()/deltaT - Sp.field())*psiMinn
+            )
           + (V0().field()/deltaT)*rho.oldTime().field()*psi0
           - sumPhiBD;
     }
@@ -472,16 +480,17 @@ void Foam::MULES::limiter
             V
            *(
                (rho.field()/deltaT - Sp.field())*psiMaxn
-             - (rho.oldTime().field()/deltaT)*psi0
              - Su.field()
+             - (rho.oldTime().field()/deltaT)*psi0
             )
           + sumPhiBD;
 
         psiMinn =
             V
            *(
-               (rho.field()/deltaT)*psi0
-             - (rho.oldTime().field()/deltaT - Sp.field())*psiMinn + Su.field()
+               Su.field()
+             - (rho.field()/deltaT - Sp.field())*psiMinn
+             + (rho.oldTime().field()/deltaT)*psi0
             )
           - sumPhiBD;
     }