multiphaseEulerFoam: Use pressureControl

pMin and pMax settings are now available in multiphaseEulerFoam in the
PIMPLE section of the system/fvOptions file. This is consistent with
other compressible solvers. The pMin setting in system/phaseProperties
is no longer read, and it's presence will result in a warning.
This commit is contained in:
Will Bainbridge
2021-03-16 17:21:23 +00:00
parent 96f07b0137
commit 762fb48ddf
74 changed files with 155 additions and 184 deletions

View File

@ -11,13 +11,6 @@ autoPtr<phaseSystem> fluidPtr
phaseSystem& fluid = fluidPtr();
phaseSystem::phaseModelList& phases = fluid.phases();
dimensionedScalar pMin
(
"pMin",
dimPressure,
fluid
);
#include "gh.H"
volScalarField& p = phases[0].thermoRef().p();
@ -36,32 +29,41 @@ volScalarField p_rgh
mesh
);
label pRefCell = 0;
scalar pRefValue = 0.0;
if (fluid.found("pMin"))
{
IOWarningInFunction(fluid)
<< "Pressure limits, pMin and pMax, are now read from "
<< pimple.dict().name() << endl;
}
pressureControl pressureControl
(
p,
p_rgh,
fluid.rho(),
pimple.dict(),
fluid.incompressible()
);
if (fluid.incompressible())
{
p = max(p_rgh + fluid.rho()*gh, pMin);
if (p_rgh.needReference())
{
setRefCell
(
p,
p_rgh,
pimple.dict(),
pRefCell,
pRefValue
);
p += dimensionedScalar
(
"p",
p.dimensions(),
pRefValue - getRefCellValue(p, pRefCell)
);
p_rgh = p - fluid.rho()*gh;
}
p = p_rgh + fluid.rho()*gh;
pressureControl.limit(p);
}
if (p_rgh.needReference() && fluid.incompressible())
{
p += dimensionedScalar
(
"p",
p.dimensions(),
pressureControl.refValue()
- getRefCellValue(p, pressureControl.refCell())
);
}
p_rgh = p - fluid.rho()*gh;
mesh.setFluxRequired(p_rgh.name());
PtrList<volScalarField> rAUs;

View File

@ -38,6 +38,7 @@ Description
#include "phaseSystem.H"
#include "phaseCompressibleMomentumTransportModel.H"
#include "pimpleControl.H"
#include "pressureControl.H"
#include "localEulerDdtScheme.H"
#include "fvcSmooth.H"

View File

@ -257,7 +257,11 @@ while (pimple.correct())
if (fluid.incompressible())
{
pEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));
pEqn.setReference
(
pressureControl.refCell(),
pressureControl.refValue()
);
}
pEqn.solve();
@ -327,7 +331,8 @@ while (pimple.correct())
}
// Update and limit the static pressure
p = max(p_rgh + rho*gh, pMin);
p = p_rgh + rho*gh;
pressureControl.limit(p);
// Account for static pressure reference
if (p_rgh.needReference() && fluid.incompressible())
@ -336,7 +341,8 @@ while (pimple.correct())
(
"p",
p.dimensions(),
pRefValue - getRefCellValue(p, pRefCell)
pressureControl.refValue()
- getRefCellValue(p, pressureControl.refCell())
);
}

View File

@ -248,7 +248,11 @@ while (pimple.correct())
if (fluid.incompressible())
{
pEqn.setReference(pRefCell, getRefCellValue(p_rgh, pRefCell));
pEqn.setReference
(
pressureControl.refCell(),
pressureControl.refValue()
);
}
pEqn.solve();
@ -305,7 +309,8 @@ while (pimple.correct())
}
// Update and limit the static pressure
p = max(p_rgh + rho*gh, pMin);
p = p_rgh + rho*gh;
pressureControl.limit(p);
// Account for static pressure reference
if (p_rgh.needReference() && fluid.incompressible())
@ -314,7 +319,8 @@ while (pimple.correct())
(
"p",
p.dimensions(),
pRefValue - getRefCellValue(p, pRefCell)
pressureControl.refValue()
- getRefCellValue(p, pressureControl.refCell())
);
}

View File

@ -187,7 +187,4 @@ turbulentDispersion
(
);
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -53,6 +53,9 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}

View File

@ -144,7 +144,4 @@ turbulentDispersion
(
);
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -53,6 +53,9 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}

View File

@ -183,8 +183,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}
relaxationFactors

View File

@ -184,8 +184,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}
relaxationFactors

View File

@ -226,8 +226,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 1;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}
relaxationFactors

View File

@ -151,8 +151,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
frozenFlow yes;
}

View File

@ -198,6 +198,4 @@ virtualMass
wallLubrication
();
pMin 10000;
// ************************************************************************* //

View File

@ -57,7 +57,10 @@ PIMPLE
{
nOuterCorrectors 3;
nCorrectors 2;
faceMomentum true;
pMin 1e4;
}
relaxationFactors

View File

@ -145,8 +145,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
frozenFlow yes;
}

View File

@ -189,8 +189,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -56,10 +56,12 @@ solvers
PIMPLE
{
nOuterCorrectors 1;
nCorrectors 1;
nCorrectors 1;
pRefCell 0;
pRefValue 1e5;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}
relaxationFactors

View File

@ -141,7 +141,4 @@ turbulentDispersion
(
);
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -53,6 +53,9 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
pRefCell 0;
pRefValue 1e5;
}

View File

@ -168,8 +168,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -70,6 +70,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -141,8 +141,4 @@ phaseTransfer
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -74,7 +74,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -115,8 +115,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -89,7 +89,10 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -168,8 +168,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -70,6 +70,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -213,8 +213,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -80,6 +80,8 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -278,8 +278,4 @@ turbulentDispersion
}
);
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -70,7 +70,10 @@ PIMPLE
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
faceMomentum true;
pMin 1e4;
}
relaxationFactors

View File

@ -114,8 +114,4 @@ turbulentDispersion
(
);
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -88,7 +88,10 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
faceMomentum no;
pMin 1e4;
}
relaxationFactors

View File

@ -176,8 +176,4 @@ phaseTransfer
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -65,7 +65,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -216,8 +216,4 @@ phaseTransfer
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -73,7 +73,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -259,8 +259,4 @@ phaseTransfer
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -82,7 +82,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -373,8 +373,4 @@ surfaceTension
phaseTransfer
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -82,7 +82,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -206,6 +206,4 @@ turbulentDispersion
interfaceCompression
();
pMin 10000;
// ************************************************************************* //

View File

@ -61,6 +61,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -173,8 +173,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -61,6 +61,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -224,8 +224,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -71,6 +71,8 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -232,7 +232,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -68,6 +68,8 @@ PIMPLE
nOuterCorrectors 5;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -191,8 +191,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -60,6 +60,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -197,8 +197,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -77,6 +77,8 @@ PIMPLE
{
nCorrectors 3;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -111,8 +111,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -88,7 +88,10 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
faceMomentum no;
pMin 1e4;
}
relaxationFactors

View File

@ -169,8 +169,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -61,6 +61,8 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pMin 1e4;
}
relaxationFactors

View File

@ -468,8 +468,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -80,6 +80,8 @@ PIMPLE
nCorrectors 3;
nNonOrthogonalCorrectors 0;
pMin 1e4;
pRefCell 0;
pRefValue 0;
}

View File

@ -468,8 +468,4 @@ wallLubrication
turbulentDispersion
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -79,6 +79,7 @@ PIMPLE
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
correctPhi yes;
pRefCell 0;

View File

@ -138,8 +138,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -82,7 +82,10 @@ PIMPLE
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 2;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -213,8 +213,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -84,7 +84,10 @@ PIMPLE
nCorrectors 1;
nEnergyCorrectors 1;
nNonOrthogonalCorrectors 0;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -225,8 +225,4 @@ turbulentDispersion
interfaceCompression
();
// Minimum allowable pressure
pMin 10000;
// ************************************************************************* //

View File

@ -84,7 +84,10 @@ PIMPLE
nCorrectors 1;
nEnergyCorrectors 1;
nNonOrthogonalCorrectors 0;
faceMomentum yes;
pMin 1e4;
}
relaxationFactors

View File

@ -191,6 +191,4 @@ turbulentDispersion
interfaceCompression
();
pMin 10000;
// ************************************************************************* //

View File

@ -61,7 +61,10 @@ PIMPLE
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
partialElimination true;
pMin 1e4;
}
relaxationFactors