47 const dictionary& XiEqProperties,
48 const psiuReactionThermo&
thermo,
50 const volScalarField&
Su
54 XiEqCoef_(XiEqModelCoeffs_.
get<scalar>(
"XiEqCoef")),
55 XiEqExp_(XiEqModelCoeffs_.
get<scalar>(
"XiEqExp")),
56 lCoef_(XiEqModelCoeffs_.
get<scalar>(
"lCoef")),
58 uPrimeCoef_(XiEqModelCoeffs_.
get<scalar>(
"uPrimeCoef")),
59 subGridSchelkin_(XiEqModelCoeffs_.
get<
bool>(
"subGridSchelkin")),
62 Su.
mesh().lookupObject<IOdictionary>(
"combustionProperties"),
84 up.primitiveFieldRef() += calculateSchelkinEffect(uPrimeCoef_);
94 tmp<volScalarField> tXiEq
114 if (Ma[celli] > 0.01)
117 XiEqCoef_*
pow(
K[celli]*Ma[celli], -XiEqExp_)*upBySu[celli];
123 forAll(xieq.boundaryField(), patchi)
127 const scalarField& Map = Ma.boundaryField()[patchi];
128 const scalarField& upBySup = upBySu.boundaryField()[patchi];
132 if (Ma[facei] > 0.01)
135 XiEqCoef_*
pow(Kp[facei]*Map[facei], -XiEqExp_)
149 XiEqModelCoeffs_.readEntry(
"XiEqCoef", XiEqCoef_);
150 XiEqModelCoeffs_.readEntry(
"XiEqExp", XiEqExp_);
151 XiEqModelCoeffs_.readEntry(
"lCoef", lCoef_);
152 XiEqModelCoeffs_.readEntry(
"uPrimeCoef", uPrimeCoef_);
153 XiEqModelCoeffs_.readEntry(
"subGridSchelkin", subGridSchelkin_);
CGAL::Exact_predicates_exact_constructions_kernel K
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Internal & ref(const bool updateAccessTime=true)
Return a reference to the dimensioned internal field.
GeometricBoundaryField< scalar, fvPatchField, volMesh > Boundary
Type of boundary fields.
virtual bool read()
Re-read model coefficients if they have changed.
Simple SCOPEXiEq model for XiEq based on SCOPEXiEqs correlation with a linear correction function to ...
virtual ~SCOPEXiEq()
Destructor.
virtual tmp< volScalarField > XiEq() const
Return the flame-wrinkling XiEq.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
A class for managing temporary objects.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
compressible::turbulenceModel & turbulence
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
tmp< GeometricField< Type, faPatchField, areaMesh > > average(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Area-weighted average a edgeField creating a areaField.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionSet dimless
Dimensionless.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
static constexpr const zero Zero
Global zero (0)
#define forAll(list, i)
Loop across all elements in list.