Skip to content
Snippets Groups Projects
Commit 38257437 authored by Henry's avatar Henry Committed by Andrew Heather
Browse files

rhoPimpleDyMFoam: Update correctPhi for morphing meshes

parent 61c85c97
No related branches found
No related tags found
No related merge requests found
...@@ -40,8 +40,6 @@ if (pimple.transonic()) ...@@ -40,8 +40,6 @@ if (pimple.transonic())
fvOptions(psi, p, rho.name()) fvOptions(psi, p, rho.name())
); );
fvOptions.constrain(pEqn);
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter()))); pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
if (pimple.finalNonOrthogonalIter()) if (pimple.finalNonOrthogonalIter())
...@@ -74,8 +72,6 @@ else ...@@ -74,8 +72,6 @@ else
fvOptions(psi, p, rho.name()) fvOptions(psi, p, rho.name())
); );
fvOptions.constrain(pEqn);
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter()))); pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
if (pimple.finalNonOrthogonalIter()) if (pimple.finalNonOrthogonalIter())
......
...@@ -2,7 +2,7 @@ ...@@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
...@@ -24,6 +24,9 @@ License ...@@ -24,6 +24,9 @@ License
Application Application
rhoPimpleFoam rhoPimpleFoam
Group
grpCompressibleSolvers grpMovingMeshSolvers
Description Description
Transient solver for laminar or turbulent flow of compressible fluids Transient solver for laminar or turbulent flow of compressible fluids
for HVAC and similar applications. for HVAC and similar applications.
...@@ -71,17 +74,21 @@ int main(int argc, char *argv[]) ...@@ -71,17 +74,21 @@ int main(int argc, char *argv[])
#include "setDeltaT.H" #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. // Store momentum to set rhoUf for introduced faces.
volVectorField rhoU("rhoU", rho*U); 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 // Do any mesh changes
mesh.update(); mesh.update();
...@@ -102,12 +109,9 @@ int main(int argc, char *argv[]) ...@@ -102,12 +109,9 @@ int main(int argc, char *argv[])
#include "meshCourantNo.H" #include "meshCourantNo.H"
} }
if (pimple.nCorrPIMPLE() <= 1) #include "rhoEqn.H"
{ Info<< "rhoEqn max/min : " << max(rho).value()
#include "rhoEqn.H" << " " << min(rho).value() << endl;
Info<< "rhoEqn max/min : " << max(rho).value()
<< " " << min(rho).value() << endl;
}
// --- Pressure-velocity PIMPLE corrector loop // --- Pressure-velocity PIMPLE corrector loop
while (pimple.loop()) while (pimple.loop())
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment