73 int main(
int argc,
char *argv[])
82 #include "readCombustionProperties.H" 83 #include "readGravitationalAcceleration.H" 84 #include "createFields.H" 85 #include "initContinuityErrs.H" 87 #include "compressibleCourantNo.H" 88 #include "setInitialDeltaT.H" 95 Info<<
"\nStarting time loop\n" <<
endl;
97 bool hasChanged =
false;
102 #include "compressibleCourantNo.H" 103 #include "setDeltaT.H" 111 tmagGradP()/
max(tmagGradP())
113 normalisedGradP.writeOpt() = IOobject::AUTO_WRITE;
118 Info<<
"\n\nTime = " << runTime.timeName() <<
endl;
125 PackedBoolList& protectedCell =
126 refCast<dynamicRefineFvMesh>(
mesh).protectedCell();
128 if (protectedCell.empty())
130 protectedCell.setSize(
mesh.nCells());
136 if (
betav[celli] < 0.99)
138 protectedCell[celli] = 1;
146 bool meshChanged =
mesh.update();
154 if (runTime.write() && hasChanged)
201 Info<<
"\nExecutionTime = " 202 << runTime.elapsedCpuTime()
autoPtr< compressible::turbulenceModel > turbulence
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
#define forAll(list, i)
Loop across all elements in list.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
const dictionary & pimple
Ostream & endl(Ostream &os)
Add newline and flush stream.
Read the control parameters used by setDeltaT.
GeometricField< vector, fvPatchField, volMesh > volVectorField
Info<< "Reading strained laminar flame speed field Su\n"<< endl;volScalarField Su(IOobject("Su", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field betav\n"<< endl;volScalarField betav(IOobject("betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Lobs\n"<< endl;volScalarField Lobs(IOobject("Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field CT\n"<< endl;volSymmTensorField CT(IOobject("CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Nv\n"<< endl;volScalarField Nv(IOobject("Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field nsv\n"<< endl;volSymmTensorField nsv(IOobject("nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);IOdictionary PDRProperties(IOobject("PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));autoPtr< PDRDragModel > drag
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const volScalarField & betav
Bound the given scalar field if it has gone unbounded.
const dimensionedScalar h
Planck constant.
void makeAbsolute(surfaceScalarField &phi, const volVectorField &U)
Make the given flux absolute.
dimensioned< scalar > mag(const dimensioned< Type > &)
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
Read the control parameters used by setDeltaT.
autoPtr< XiModel > flameWrinkling
Create the flame-wrinkling model.