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" 94 Info<<
"\nStarting time loop\n" <<
endl;
96 bool hasChanged =
false;
101 #include "compressibleCourantNo.H" 102 #include "setDeltaT.H" 110 tmagGradP()/
max(tmagGradP())
112 normalisedGradP.writeOpt() = IOobject::AUTO_WRITE;
117 Info<<
"\n\nTime = " << runTime.timeName() <<
endl;
124 PackedBoolList& protectedCell =
125 refCast<dynamicRefineFvMesh>(
mesh).protectedCell();
127 if (protectedCell.empty())
129 protectedCell.setSize(
mesh.nCells());
135 if (
betav[cellI] < 0.99)
137 protectedCell[cellI] = 1;
145 bool meshChanged =
mesh.update();
153 if (runTime.write() && hasChanged)
200 Info<<
"\nExecutionTime = " 201 << runTime.elapsedCpuTime()
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
dimensioned< scalar > mag(const dimensioned< Type > &)
void makeAbsolute(surfaceScalarField &phi, const volVectorField &U)
Make the given flux absolute.
Bound the given scalar field if it has gone unbounded.
const volScalarField & betav
autoPtr< compressible::turbulenceModel > turbulence
const dictionary & pimple
Ostream & endl(Ostream &os)
Add newline and flush stream.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
GeometricField< vector, fvPatchField, volMesh > volVectorField
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Read the control parameters used by setDeltaT.
int main(int argc, char *argv[])
autoPtr< XiModel > flameWrinkling
Create the flame-wrinkling model.
const dimensionedScalar h
Planck constant.