45 Foam::XiEqModels::Gulder::Gulder
47 const dictionary& XiEqProperties,
48 const psiuReactionThermo&
thermo,
54 XiEqCoef_(XiEqModelCoeffs_.
get<scalar>(
"XiEqCoef")),
56 uPrimeCoef_(XiEqModelCoeffs_.
get<scalar>(
"uPrimeCoef")),
57 subGridSchelkin_(XiEqModelCoeffs_.
get<bool>(
"subGridSchelkin"))
76 up.primitiveFieldRef() += calculateSchelkinEffect(uPrimeCoef_);
90 return (1.0 + XiEqCoef_*
sqrt(up/(Su_ + SuMin_))*Reta);
98 XiEqModelCoeffs_.readEntry(
"XiEqCoef", XiEqCoef_);
99 XiEqModelCoeffs_.readEntry(
"uPrimeCoef", uPrimeCoef_);
100 XiEqModelCoeffs_.readEntry(
"subGridSchelkin", subGridSchelkin_);
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
tmp< GeometricField< Type, faPatchField, areaMesh > > average(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Area-weighted average a edgeField creating a areaField.
virtual tmp< volScalarField > XiEq() const
Return the flame-wrinkling XiEq.
virtual bool read(const dictionary &XiEqProperties)=0
Update properties from given dictionary.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
virtual bool read(const dictionary &XiEqProperties)
Update properties from given dictionary.
const dimensionedScalar e
Elementary charge.
virtual ~Gulder()
Destructor.
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
defineTypeNameAndDebug(combustionModel, 0)
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
A class for managing temporary objects.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)