From 1d64955b6955e3f5d69af986205c5401f7282d9d Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 21 Mar 2013 17:05:00 +0000
Subject: [PATCH] rhoReactingFoam: Correct handling of the density correction
 term for the transonic case

---
 .../solvers/combustion/reactingFoam/rhoReactingFoam/pEqn.H      | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/pEqn.H b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/pEqn.H
index 3fb6dfe45e4..7a5c717cff7 100644
--- a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/pEqn.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/pEqn.H
@@ -27,7 +27,7 @@
         fvScalarMatrix pDDtEqn
         (
             fvc::ddt(rho) + fvc::div(phiHbyA)
-          + correction(fvm::ddt(psi, p) + fvm::div(phid, p))
+          + correction(psi*fvm::ddt(p) + fvm::div(phid, p))
         );
 
         while (pimple.correctNonOrthogonal())
-- 
GitLab