29 const dictionary& pimpleDict =
pimple.dict();
35 scalar
maxDeltaT(pimpleDict.lookupOrDefault<scalar>(
"maxDeltaT", GREAT));
40 pimpleDict.lookupOrDefault<scalar>(
"rDeltaTSmoothingCoeff", 0.1)
46 pimpleDict.lookupOrDefault<scalar>(
"rDeltaTDampingCoeff", 1.0)
51 scalar
alphaTemp(pimpleDict.lookupOrDefault(
"alphaTemp", 0.05));
55 scalar
alphaY(pimpleDict.lookupOrDefault(
"alphaY", 1.0));
57 Info<<
"Time scales min/max:" <<
endl;
74 << 1/
gMax(rDeltaT.primitiveField()) <<
", " 75 << 1/
gMin(rDeltaT.primitiveField()) << endl;
81 volScalarField::Internal rDeltaTT
86 Info<<
" Temperature = " 87 << 1/(
gMax(rDeltaTT.field()) + VSMALL) <<
", " 88 << 1/(
gMin(rDeltaTT.field()) + VSMALL) <<
endl;
90 rDeltaT.ref() =
max(rDeltaT(), rDeltaTT);
96 dictionary Yref(pimpleDict.subDict(
"Yref"));
98 volScalarField::Internal rDeltaTY
118 if (Yref.found(Yi.name()))
121 scalar Yrefi =
readScalar(Yref.lookup(Yi.name()));
123 rDeltaTY.field() =
max 138 Info<<
" Composition = " 139 << 1/(
gMax(rDeltaTY.field()) + VSMALL) <<
", " 140 << 1/(
gMin(rDeltaTY.field()) + VSMALL) <<
endl;
142 rDeltaT.ref() =
max(rDeltaT(), rDeltaTY);
147 <<
"Cannot find any active species in Yref " << Yref
153 rDeltaT.correctBoundaryConditions();
167 && runTime.timeIndex() > runTime.startTimeIndex() + 1
178 rDeltaT.correctBoundaryConditions();
181 << 1/
gMax(rDeltaT.primitiveField())
182 <<
", " << 1/
gMin(rDeltaT.primitiveField()) << endl;
scalar maxCo(readScalar(pimpleDict.lookup("maxCo")))
const word & executable() const
Name of executable without the path.
scalar maxDeltaT(pimpleDict.lookupOrDefault< scalar >("maxDeltaT", GREAT))
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Type gMin(const FieldField< Field, Type > &f)
basicMultiComponentMixture & composition
const dictionary & pimple
Ostream & endl(Ostream &os)
Add newline and flush stream.
void smooth(volScalarField &field, const scalar coeff)
volScalarField rDeltaT0("rDeltaT0", rDeltaT)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
#define IOWarningIn(functionName, ios)
Report an IO warning using Foam::Warning.
scalar rDeltaTDampingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTDampingCoeff", 1.0))
Type gMax(const FieldField< Field, Type > &f)
scalar alphaTemp(pimpleDict.lookupOrDefault("alphaTemp", 0.05))
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Info<< "Creating reaction model\"<< endl;autoPtr< combustionModels::psiCombustionModel > reaction(combustionModels::psiCombustionModel::New(mesh))
PtrList< volScalarField > & Y
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< volScalarField > trDeltaT
Foam::argList args(argc, argv)
const label inertIndex(composition.species()[inertSpecie])
scalar rDeltaTSmoothingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTSmoothingCoeff", 0.1))
scalar alphaY(pimpleDict.lookupOrDefault("alphaY", 1.0))