From 3bffa0e552fc998d322117985228d28dde206b24 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 28 Mar 2011 10:31:24 +0100
Subject: [PATCH] ENH: (re-)added under-relaxation to pEqn in
LTSReactingParcelFoam solver
---
.../lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C | 5 +++++
applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H | 4 +++-
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
index dc1c09e7440..9cb7908a147 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C
@@ -93,6 +93,11 @@ int main(int argc, char *argv[])
bool finalIter = oCorr == nOuterCorr-1;
#include "addFinalIter.H"
+ if (nOuterCorr != 1)
+ {
+ p.storePrevIter();
+ }
+
turbulence->correct();
#include "UEqn.H"
diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
index c25c63df0b3..97ba0c5c781 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
@@ -59,7 +59,7 @@
}
}
- Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
+ p.relax();
// Second part of thermodynamic density update
thermo.rho() += psi*p;
@@ -75,4 +75,6 @@
rho = min(rho, rhoMax);
#include "setPressureWork.H"
+
+ Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
}
--
GitLab