57 int main(
int argc,
char *argv[])
64 #include "createControl.H" 66 #include "../interFoam/interDyMFoam/createDyMControls.H" 67 #include "initContinuityErrs.H" 68 #include "createFields.H" 78 IOobject::READ_IF_PRESENT,
86 #include "CourantNo.H" 87 #include "setInitialDeltaT.H" 93 Info<<
"\nStarting time loop\n" <<
endl;
97 #include "../interFoam/interDyMFoam/readControls.H" 104 #include "CourantNo.H" 105 #include "setDeltaT.H" 109 Info<<
"Time = " << runTime.timeName() <<
nl <<
endl;
116 scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
122 Info<<
"Execution time for mesh.update() = " 123 << runTime.elapsedCpuTime() - timeBeforeMeshUpdate
135 #include "correctPhi.H" 147 #include "alphaControls.H" 163 #include "alphaEqnSubCycle.H" 182 Info<<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s" 183 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s"
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
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Creates and initialises the velocity velocity field Uf.
const dimensionedVector & g
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating phaseChangeTwoPhaseMixture\"<< 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.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & gh
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Execute application functionObjects to post-process existing results.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
Read the control parameters used by setDeltaT.
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
volScalarField rAU(1.0/UEqn.A())