From 382574378fee97e3d79e25211f1dcdf9e6d24db9 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 19 Jun 2014 15:19:13 +0100
Subject: [PATCH] rhoPimpleDyMFoam: Update correctPhi for morphing meshes
---
.../rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H | 4 ---
.../rhoPimpleDyMFoam/rhoPimpleDyMFoam.C | 30 +++++++++++--------
2 files changed, 17 insertions(+), 17 deletions(-)
diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H
index 5270e3d85bb..cb7329026c1 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H
@@ -40,8 +40,6 @@ if (pimple.transonic())
fvOptions(psi, p, rho.name())
);
- fvOptions.constrain(pEqn);
-
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
if (pimple.finalNonOrthogonalIter())
@@ -74,8 +72,6 @@ else
fvOptions(psi, p, rho.name())
);
- fvOptions.constrain(pEqn);
-
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
if (pimple.finalNonOrthogonalIter())
diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
index e633ebee602..fcfc2aa7952 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -24,6 +24,9 @@ License
Application
rhoPimpleFoam
+Group
+ grpCompressibleSolvers grpMovingMeshSolvers
+
Description
Transient solver for laminar or turbulent flow of compressible fluids
for HVAC and similar applications.
@@ -71,17 +74,21 @@ int main(int argc, char *argv[])
#include "setDeltaT.H"
- runTime++;
+ {
+ // Store divrhoU from the previous time-step/mesh for the correctPhi
+ volScalarField divrhoU
+ (
+ "divrhoU",
+ fvc::div(fvc::absolute(phi, rho, U))
+ );
- Info<< "Time = " << runTime.timeName() << nl << endl;
+ runTime++;
+
+ Info<< "Time = " << runTime.timeName() << nl << endl;
- {
// Store momentum to set rhoUf for introduced faces.
volVectorField rhoU("rhoU", rho*U);
- // Store divrhoU from the previous time-step/mesh for the correctPhi
- volScalarField divrhoU(fvc::div(fvc::absolute(phi, rho, U)));
-
// Do any mesh changes
mesh.update();
@@ -102,12 +109,9 @@ int main(int argc, char *argv[])
#include "meshCourantNo.H"
}
- if (pimple.nCorrPIMPLE() <= 1)
- {
- #include "rhoEqn.H"
- Info<< "rhoEqn max/min : " << max(rho).value()
- << " " << min(rho).value() << endl;
- }
+ #include "rhoEqn.H"
+ Info<< "rhoEqn max/min : " << max(rho).value()
+ << " " << min(rho).value() << endl;
// --- Pressure-velocity PIMPLE corrector loop
while (pimple.loop())
--
GitLab