From 692242f574427183cd2384add7657d09ad050ffc Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Fri, 10 Mar 2017 11:08:10 +0000
Subject: [PATCH] CrankNicolsonDdtScheme: Corrected fvcDdtUfCorr for
 compressible flow

Resolves bug-report https://bugs.openfoam.org/view.php?id=2491
---
 .../CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C           | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
index 0c7fe053b..dd12b2552 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
@@ -1255,7 +1255,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtUfCorr
             ddt0_<GeometricField<Type, fvPatchField, volMesh>>
             (
                 "ddtCorrDdt0(" + rho.name() + ',' + U.name() + ')',
-                U.dimensions()
+                rho.dimensions()*U.dimensions()
             );
 
         DDt0Field<GeometricField<Type, fvsPatchField, surfaceMesh>>& dUfdt0 =
@@ -1350,7 +1350,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtPhiCorr
             ddt0_<GeometricField<Type, fvPatchField, volMesh>>
             (
                 "ddtCorrDdt0(" + rho.name() + ',' + U.name() + ')',
-                U.dimensions()
+                rho.dimensions()*U.dimensions()
             );
 
         DDt0Field<fluxFieldType>& dphidt0 =