From e9d239cdef3416beb9555d81f620a13bf4c908bc Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Tue, 16 Sep 2008 15:21:35 +0100
Subject: [PATCH] Correctad ddtPhiCorr for the case where phi is volumetric.

---
 .../CrankNicholsonDdtScheme/CrankNicholsonDdtScheme.C  | 10 +++-------
 1 file changed, 3 insertions(+), 7 deletions(-)

diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicholsonDdtScheme/CrankNicholsonDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicholsonDdtScheme/CrankNicholsonDdtScheme.C
index f591376bf5..10299f6175 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicholsonDdtScheme/CrankNicholsonDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicholsonDdtScheme/CrankNicholsonDdtScheme.C
@@ -1027,8 +1027,7 @@ CrankNicholsonDdtScheme<Type>::fvcDdtPhiCorr
             {
                 dUdt0 = rDtCoef0_(dUdt0)*
                 (
-                    rho.oldTime()*U.oldTime()
-                  - rho.oldTime().oldTime()*U.oldTime().oldTime()
+                    U.oldTime() - U.oldTime().oldTime()
                 ) - offCentre_(dUdt0());
             }
 
@@ -1054,11 +1053,8 @@ CrankNicholsonDdtScheme<Type>::fvcDdtPhiCorr
                       - (
                             fvc::interpolate
                             (
-                                rA*
-                                (
-                                    rDtCoef*rho.oldTime()*U.oldTime()
-                                    + offCentre_(dUdt0())
-                                )
+                                rA*rho.oldTime()
+                               *(rDtCoef*U.oldTime() + offCentre_(dUdt0()))
                             ) & mesh().Sf()
                         )
                     )