BUG: non-orthogonal correction corrected in pEq

2025-11-28 03:28:01 +00:00 · 2011-02-21 16:54:28 +00:00
parent e820c30b60
commit 1ca7097ad6
1 changed files with 10 additions and 5 deletions
--- 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();