From e4cecfcdb5a583b633921b7aea94a92e178e3561 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Sun, 18 May 2014 21:19:26 +0100
Subject: [PATCH] rhoPimpleDyMFoam/correctPhi.H: Add compressibility term in
pcorr equation
---
.../rhoPimpleFoam/rhoPimpleDyMFoam/correctPhi.H | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/correctPhi.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/correctPhi.H
index e916b7a56aa..c704573276b 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/correctPhi.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/correctPhi.H
@@ -46,14 +46,18 @@ if (mesh.changing())
{
fvScalarMatrix pcorrEqn
(
- fvm::laplacian(rAUf, pcorr) == fvc::div(phi) - divrhoU
+ fvm::ddt(psi, pcorr)
+ + fvc::div(phi)
+ - fvm::laplacian(rAUf, pcorr)
+ ==
+ divrhoU
);
pcorrEqn.solve();
if (pimple.finalNonOrthogonalIter())
{
- phi -= pcorrEqn.flux();
+ phi += pcorrEqn.flux();
}
}
}
--
GitLab