45 Foam::XiGModels::basicSubGrid::basicSubGrid
47 const dictionary& XiGProperties,
48 const psiuReactionThermo&
thermo,
54 k1(XiGModelCoeffs_.
get<scalar>(
"k1")),
69 const objectRegistry& db = Su_.db();
74 tmp<volScalarField> tGtot = XiGModel_->G();
84 Gtot[celli] += k1*
mag(
U[celli])/Lobs[celli];
94 const objectRegistry& db = Su_.db();
100 return XiGModel_->Db()
101 +
rho*Su_*(Xi - 1.0)*mgb*(0.5*Lobs)*Lobs/(mgb*Lobs + 1.0);
109 XiGModelCoeffs_.readEntry(
"k1", k1);
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)
virtual tmp< volScalarField > G() const
Return the flame-wrinkling generation rate.
virtual bool read(const dictionary &XiGProperties)
Update properties from given dictionary.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
GeometricField< vector, fvPatchField, volMesh > volVectorField
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
const dimensionedScalar e
Elementary charge.
virtual tmp< volScalarField > Db() const
Return the flame diffusivity.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
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
const Vector< label > N(dict.get< Vector< label >>("N"))
virtual bool read(const dictionary &XiGProperties)=0
Update properties from given dictionary.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual ~basicSubGrid()
Destructor.
Internal & ref(const bool updateAccessTime=true)
Same as internalFieldRef()
A class for managing temporary objects.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)