36 namespace radiationModels
38 namespace absorptionEmissionModels
63 coeffsDict_(
dict.subDict(modelName +
"Coeffs")),
65 specieIndex_(
label(0)),
70 if (!isA<basicSpecieMixture>(
thermo_))
73 <<
"Model requires a multi-component thermo package"
85 const word& key = iter().keyword();
110 <<
"specie ft is not present to use with "
111 <<
"lookUpTableFileName " <<
nl
129 Info<<
"specie: " << iter.key() <<
" found on look-up table "
130 <<
" with index: " << index <<
endl;
139 Info<<
"specie: " << iter.key() <<
" is being solved" <<
endl;
144 <<
"specie: " << iter.key()
145 <<
" is neither in look-up table: "
147 <<
" nor is being solved" <<
nl
160 <<
" there is not lookup table and the specie" <<
nl
162 <<
" is not found " <<
nl
194 "aCont" +
name(bandI),
197 extrapolatedCalculatedFvPatchVectorField::typeName
209 if (specieIndex_[
n] != 0)
215 const List<scalar>& Ynft = lookUpTablePtr_().lookUp(ft[celli]);
217 Xipi = Ynft[specieIndex_[
n]]*
paToAtm(
p[celli]);
235 scalar Ti =
T[celli];
237 if (coeffs_[
n].invTemp())
244 ((((
b[5]*Ti +
b[4])*Ti +
b[3])*Ti +
b[2])*Ti +
b[1])*Ti
249 ta.
ref().correctBoundaryConditions();
#define forAll(list, i)
Loop across all elements in list.
#define forAllConstIter(Container, container, iter)
Iterate across all elements in the container object of type.
Macros for easy insertion into run-time selection tables.
Generic GeometricField class.
static tmp< GeometricField< Type, PatchField, GeoMesh > > New(const word &name, const Internal &, const PtrList< PatchField< Type >> &)
Return a temporary field constructed from name,.
bool insert(const Key &, const T &newElmt)
Insert a new hashedEntry.
static const word & constant()
Return constant name.
Specialisation of basicMixture for a mixture consisting of a number for molecular species.
A list of keyword definitions, which are a keyword followed by any number of values (e....
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
bool found(const word &, bool recursive=false, bool patternMatch=true) const
Search dictionary for given keyword.
A class for handling file names.
Base-class for fluid thermodynamic properties.
Mesh data needed to do the Finite Volume discretisation.
const Time & time() const
Return the top-level database.
Type & lookupObjectRef(const word &name) const
Lookup and return the object reference of the given Type.
bool foundObject(const word &name) const
Is the named Type in registry.
An abstract base class for physical properties.
void initialise(const dictionary &)
Model to supply absorption and emission coefficients for radiation modelling.
const dictionary & dict() const
Reference to the dictionary.
const fvMesh & mesh() const
Reference to the mesh.
virtual tmp< volScalarField > ECont(const label bandI=0) const
Emission contribution for continuous phase.
const fluidThermo & thermo_
Thermo package.
UPtrList< volScalarField > Yj_
Pointer list of species in the registry involved in the absorption.
HashTable< label > speciesNames_
Hash table of species names.
autoPtr< interpolationLookUpTable > lookUpTablePtr_
Look-up table of species related to ft.
virtual ~greyMean()
Destructor.
FixedList< label, nSpecies_ > specieIndex_
Indices of species in the look-up table.
tmp< volScalarField > eCont(const label bandI=0) const
Emission coefficient for continuous phase.
absorptionCoeffs coeffs_[nSpecies_]
tmp< volScalarField > aCont(const label bandI=0) const
Absorption coefficient for continuous phase.
tmp< volScalarField > ECont(const label bandI=0) const
Emission contribution for continuous phase.
dictionary coeffsDict_
Absorption model dictionary.
greyMean(const dictionary &dict, const fvMesh &mesh, const word &modelName=typeName)
Construct from components.
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
A class for handling words, derived from string.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.name(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
defineTypeNameAndDebug(greyMeanCombustion, 0)
addToRunTimeSelectionTable(absorptionEmissionModel, greyMeanCombustion, dictionary)
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Ostream & endl(Ostream &os)
Add newline and flush stream.
word name(const bool)
Return a word representation of a bool.
errorManip< error > abort(error &err)
const dimensionSet dimless
const dimensionSet dimLength
scalar paToAtm(const scalar pa)
Conversion from atm to Pa.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
Unit conversion functions.