35 namespace radiationModels
37 namespace absorptionEmissionModels
62 coeffsDict_(
dict.subDict(modelName +
"Coeffs")),
64 specieIndex_(
label(0)),
73 if (!iter().isDict())
continue;
82 if (!iter().isDict())
continue;
84 const word& key = iter().keyword();
94 <<
"specie: " << key <<
" is not in all the bands"
123 <<
"specie ft is not present to use with "
124 <<
"lookUpTableFileName " <<
nl
143 Info<<
"specie: " << iter.key() <<
" found on look-up table "
144 <<
" with index: " << index <<
endl;
153 Info<<
"specie: " << iter.key() <<
" is being solved" <<
endl;
158 <<
"specie: " << iter.key()
159 <<
" is neither in look-up table: "
161 <<
" nor is being solved" <<
nl
174 <<
" there is no lookup table and the specie" <<
nl
176 <<
" is not found " <<
nl
222 if (specieIndex_[
n] != 0)
227 const List<scalar>& Ynft = lookUpTablePtr_().lookUp(ft[celli]);
230 Xipi = Ynft[specieIndex_[
n]]*
paToAtm(
p[celli]);
248 scalar Ti =
T[celli];
251 coeffs_[bandi][
n].coeffs(
T[celli]);
253 if (coeffs_[bandi][
n].invTemp())
261 ((((
b[5]*Ti +
b[4])*Ti +
b[3])*Ti +
b[2])*Ti +
b[1])*Ti
299 for (
label j=0; j<nBands_; j++)
301 aLambda[j].primitiveFieldRef() = this->a(j);
304 aLambda[j].primitiveField()
305 *(iBands_[j][1] - iBands_[j][0])
#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.
Internal::FieldType & primitiveFieldRef()
Return a reference to the internal field.
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.
bool found(const Key &) const
Return true if hashedEntry is found in table.
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
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.
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
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.
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.
void correct(volScalarField &a, PtrList< volScalarField > &aLambda) const
Correct rays.
UPtrList< volScalarField > Yj_
Pointer list of species being solved involved in the absorption.
FixedList< FixedList< absorptionCoeffs, nSpecies_ >, maxBands_ > coeffs_
Absorption coefficients.
tmp< volScalarField > aCont(const label bandi=0) const
Absorption coefficient for continuous phase.
HashTable< label > speciesNames_
Hash table with species names.
wideBand(const dictionary &dict, const fvMesh &mesh, const word &modelName=typeName)
Construct from components.
autoPtr< interpolationLookUpTable > lookUpTablePtr_
Look-up table of species related to ft.
FixedList< label, nSpecies_ > specieIndex_
Indices of species in the look-up table.
tmp< volScalarField > ECont(const label bandi=0) const
Emission contribution for continuous phase.
tmp< volScalarField > eCont(const label bandi=0) const
Emission coefficient for continuous phase.
virtual ~wideBand()
Destructor.
FixedList< Vector2D< scalar >, maxBands_ > iBands_
Bands.
dictionary coeffsDict_
Absorption model dictionary.
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.
const dimensionSet dimless
const dimensionSet dimLength
scalar paToAtm(const scalar pa)
Conversion from atm to Pa.
word name(const complex &)
Return a string representation of a complex.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
Unit conversion functions.