39 Foam::XiEqModel::XiEqModel
41 const dictionary& XiEqProperties,
42 const psiuReactionThermo&
thermo,
49 XiEqProperties.subDict
51 word(XiEqProperties.
lookup(
"XiEqModel")) +
"Coeffs" 55 turbulence_(turbulence),
91 const fvMesh& mesh =
Su_.
mesh();
99 tmp<volScalarField> tN
106 mesh.time().timeName(),
117 N.primitiveFieldRef() = Nv.primitiveField()*
pow(mesh.V(), 2.0/3.0);
124 mesh.time().timeName(),
137 ns.primitiveFieldRef() = nsv.primitiveField()*
pow(mesh.V(), 2.0/3.0);
153 N.primitiveFieldRef() = upLocal*(
max(scalar(1),
pow(nr, 0.5)) - 1.0);
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
dictionary XiEqModelCoeffs_
Dictionary.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
tmp< volScalarField > calculateSchelkinEffect(const scalar) const
Return the sub-grid Schelkin effect.
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual bool read(const dictionary &XiEqProperties)=0
Update properties from given dictionary.
rhoReactionThermo & thermo
GeometricField< vector, fvPatchField, volMesh > volVectorField
virtual ~XiEqModel()
Destructor.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
stressControl lookup("compactNormalStress") >> compactNormalStress
const dictionary & optionalSubDict(const word &) const
Find and return a sub-dictionary if found.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
defineRunTimeSelectionTable(reactionRateFlameArea, dictionary)
const Mesh & mesh() const
Return mesh.
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)
dimensioned< symmTensor > dimensionedSymmTensor
Dimensioned tensor obtained from generic dimensioned type.
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
fileName::Type type(const fileName &, const bool followLink=true)
Return the file type: DIRECTORY or FILE.
Internal & ref()
Return a reference to the dimensioned internal field.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensioned< scalar > mag(const dimensioned< Type > &)
const volScalarField & Su_
Laminar burning velocity.
const doubleScalar e
Elementary charge.
A class for managing temporary objects.
virtual bool write(const bool valid=true) const
Write using setting from DB.
void writeFields() const
Write fields.