43 #include "noPhaseChange.H" 54 int main(
int argc,
char *argv[])
61 #include "initContinuityErrs.H" 63 #include "createFields.H" 64 #include "createFieldRefs.H" 66 #include "initCorrectPhi.H" 73 #include "CourantNo.H" 74 #include "setInitialDeltaT.H" 78 Info<<
"\nStarting time loop\n" <<
endl;
86 #include "setRDeltaT.H" 90 #include "CourantNo.H" 91 #include "alphaCourantNo.H" 92 #include "setDeltaT.H" 107 tmp<volScalarField>
divU;
112 && !isType<twoPhaseChangeModels::noPhaseChange>(
phaseChange)
131 if (
mesh.topoChanging())
143 #include "correctPhi.H" 171 #include "alphaControls.H" 172 #include "alphaEqnSubCycle.H" 192 Info<<
"ExecutionTime = " <<
runTime.elapsedCpuTime() <<
" s" 193 <<
" ClockTime = " <<
runTime.elapsedClockTime() <<
" s"
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.
twoPhaseChangeModel & phaseChange
virtual void preUpdateMesh()
Prepare for mesh update.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const dimensionSet dimTime
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(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
const surfaceScalarField & ghf
const dimensionSet dimMass
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.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & gh
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedVector & g
Execute application functionObjects to post-process existing results.
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
Creates and initialises the velocity field Uf if required.
tmp< surfaceScalarField > talphaPhi1Corr0