33 namespace laminarFlameSpeedModels
49 Foam::laminarFlameSpeedModels::Gulders::Gulders
62 alpha_(
readScalar(coeffsDict_.lookup(
"alpha"))),
75 inline Foam::scalar Foam::laminarFlameSpeedModels::Gulders::SuRef
82 return W_*
pow(phi, eta_)*
exp(-xi_*
sqr(phi - 1.075));
91 inline Foam::scalar Foam::laminarFlameSpeedModels::Gulders::Su0pTphi
99 static const scalar Tref = 300.0;
100 static const scalar
pRef = 1.013e5;
102 return SuRef(phi)*
pow((Tu/Tref), alpha_)*
pow((p/pRef), beta_)*(1 - f_*Yres);
127 Su0[celli] = Su0pTphi(p[celli], Tu[celli], phi, 0.0);
172 Su0[celli] = Su0pTphi(p[celli], Tu[celli], phi[celli], 0.0);
199 if (psiuReactionThermo_.composition().contains(
"ft"))
201 const volScalarField& ft = psiuReactionThermo_.composition().Y(
"ft");
205 psiuReactionThermo_.p(),
206 psiuReactionThermo_.Tu(),
209 psiuReactionThermo_.lookup(
"stoichiometricAirFuelMassRatio")
210 )*ft/
max(1 - ft, small)
217 psiuReactionThermo_.p(),
218 psiuReactionThermo_.Tu(),
#define forAll(list, i)
Loop across all elements in list.
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
const Boundary & boundaryField() const
Return const-reference to the boundary field.
T & ref() const
Return non-const reference or generate a fatal error.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Macros for easy insertion into run-time selection tables.
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Mesh &, const dimensionSet &, const word &patchFieldType=fvPatchField< scalar >::calculatedType())
Return a temporary field constructed from name, mesh, dimensionSet.
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
dimensionedScalar exp(const dimensionedScalar &ds)
const dictionary & optionalSubDict(const word &) const
Find and return a sub-dictionary if found.
addToRunTimeSelectionTable(laminarFlameSpeed, constant, dictionary)
tmp< volScalarField > operator()() const
Return the laminar flame speed [m/s].
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if successful.
Foam::psiuReactionThermo.
virtual ~Gulders()
Destructor.
const Mesh & mesh() const
Return mesh.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Abstract class for laminar flame speed.
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
A class for managing temporary objects.
defineTypeNameAndDebug(constant, 0)
const dimensionSet dimVelocity