46 int main(
int argc,
char *argv[])
54 #include "createFields.H" 55 #include "createFieldRefs.H" 59 #include "CourantNo.H" 60 #include "setInitialDeltaT.H" 65 pimple.dict().lookupOrDefault<Switch>(
"faceMomentum",
false)
67 Switch partialElimination
69 pimple.dict().lookupOrDefault<Switch>(
"partialElimination",
false)
76 Info<<
"\nStarting time loop\n" <<
endl;
78 while (
pimple.run(runTime))
84 pimple.dict().lookupOrDefault<
int>(
"nEnergyCorrectors", 1)
89 #include "setRDeltaT.H" 97 #include "CourantNo.H" 98 #include "setDeltaT.H" 106 tmp<volScalarField>
divU;
125 Info<<
"Time = " << runTime.userTimeName() <<
nl <<
endl;
137 if (
pimple.thermophysics())
141 fluid.correctContinuityError();
207 fluid.correctContinuityError();
209 if (
pimple.thermophysics())
218 if (
pimple.thermophysics())
229 if (
pimple.thermophysics())
237 fluid.correctKinematics();
241 fluid.correctTurbulence();
248 Info<<
"ExecutionTime = " 249 << runTime.elapsedCpuTime()
pressureReference & pressureReference
#define forAll(list, i)
Loop across all elements in list.
pimpleNoLoopControl & pimple
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
virtual void correct()
Correct the fvModels.
Ostream & endl(Ostream &os)
Add newline and flush stream.
virtual void preUpdateMesh()
Prepare for mesh update.
PtrList< fvScalarMatrix > pEqnComps(phases.size())
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);if(fluid.found("pMin")){ IOWarningInFunction(fluid)<< "Pressure limits, pMin and pMax, are now read from "<< pimple.dict().name()<< endl;}pressureReference pressureReference(p, p_rgh, pimple.dict(), fluid.incompressible());if(fluid.incompressible()){ p=p_rgh+fluid.rho() *gh;}if(p_rgh.needReference() &&fluid.incompressible()){ p+=dimensionedScalar("p", p.dimensions(), pressureReference.refValue() - getRefCellValue(p, pressureReference.refCell()));}p_rgh=p - fluid.rho() *gh;mesh.schemes().setFluxRequired(p_rgh.name());PtrList< volScalarField > rAUs
const surfaceScalarField & ghf
PtrList< surfaceScalarField > rAUfs
tmp< volScalarField > divU
Foam::fvModels & fvModels
Calculates and outputs the mean and maximum Courant Numbers.
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
phaseSystem::phaseModelList & phases
const volScalarField & gh
const dimensionedVector & g
Execute application functionObjects to post-process existing results.
Provides functions smooth spread and sweep which use the FvFaceCellWave algorithm to smooth and redis...