53 int main(
int argc,
char *argv[])
60 #include "initContinuityErrs.H" 61 #include "createControl.H" 62 #include "createFields.H" 64 #include "createControls.H" 65 #include "CourantNo.H" 66 #include "setInitialDeltaT.H" 76 Info<<
"\nStarting time loop\n" <<
endl;
80 #include "readControls.H" 88 #include "CourantNo.H" 89 #include "setDeltaT.H" 93 Info<<
"Time = " << runTime.timeName() <<
nl <<
endl;
95 scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
102 Info<<
"Execution time for mesh.update() = " 103 << runTime.elapsedCpuTime() - timeBeforeMeshUpdate
106 gh = (
g &
mesh.C()) - ghRef;
115 #include "correctPhi.H" 132 #include "alphaEqnsSubCycle.H" 156 p_rgh.correctBoundaryConditions();
160 Info<<
"ExecutionTime = " 161 << runTime.elapsedCpuTime()
autoPtr< compressible::turbulenceModel > turbulence
interfaceProperties interface(alpha1, U, mixture())
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
const dictionary & pimple
Ostream & endl(Ostream &os)
Add newline and flush stream.
const surfaceScalarField & ghf
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Creates and initialises the velocity velocity field Uf.
const dimensionedVector & g
Info<< "Reading field p_rgh\n"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating phaseChangeTwoPhaseMixture\n"<< endl;autoPtr< phaseChangeTwoPhaseMixture > mixture
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.
const volScalarField & gh
volScalarField divU(fvc::div(fvc::absolute(phi, U)))
MULES: Multidimensional universal limiter for explicit solution.
Execute application functionObjects to post-process existing results.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.