From 3dc6bb4ecd91acfaea86f72d82e57b293ad73787 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 21 Feb 2011 10:45:33 +0000
Subject: [PATCH] BUG: Corrected non-ortho correction in pEqn
---
.../lagrangian/steadyReactingParcelFoam/pEqn.H | 15 ++++++++++-----
1 file changed, 10 insertions(+), 5 deletions(-)
diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H
index 322289684dd..fb6c4f76d9b 100644
--- a/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H
@@ -23,16 +23,21 @@
);
}
+ fvScalarMatrix DrhoDtEqn
+ (
+ fvc::ddt(rho) + psi*correction(fvm::ddt(p))
+ + fvc::div(phi)
+ ==
+ parcels.Srho()
+ + massSource.SuTot()
+ );
+
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pEqn
(
- fvc::ddt(rho) + psi*correction(fvm::ddt(p))
- + fvc::div(phi)
+ DrhoDtEqn
- fvm::laplacian(rho*rAU, p)
- ==
- parcels.Srho()
- + massSource.SuTot()
);
pEqn.solve();
--
GitLab