61 coeffsDict_((
dict.optionalSubDict(typeName +
"Coeffs"))),
71 for (
const entry& dEntry : functionDicts)
82 totalWaveLength_ += iBands_[nBand][1] - iBands_[nBand][0];
87 for (
const entry& dEntry : specDicts)
89 const word& key = dEntry.keyword();
93 speciesNames_.
insert(key, nSpec);
95 else if (!speciesNames_.
found(key))
98 <<
"specie: " << key <<
" is not in all the bands"
110 coeffsDict_.
found(
"lookUpTableFileName")
111 &&
"none" != coeffsDict_.
get<
word>(
"lookUpTableFileName")
114 lookUpTablePtr_.reset
127 <<
"specie ft is not present to use with "
128 <<
"lookUpTableFileName " <<
nl
139 const word& specieName = iter.key();
140 const label index = iter.val();
146 if (lookUpTablePtr_().found(specieName))
148 const label fieldIndex =
149 lookUpTablePtr_().findFieldIndex(specieName);
151 Info<<
"specie: " << specieName <<
" found on look-up table "
152 <<
" with index: " << fieldIndex <<
endl;
154 specieIndex_[index] = fieldIndex;
159 specieIndex_[index] = 0;
161 Info<<
"specie: " << specieName <<
" is being solved" <<
endl;
166 <<
"specie: " << specieName
167 <<
" is neither in look-up table: "
168 << lookUpTablePtr_().tableName()
169 <<
" nor is being solved" <<
nl
176 specieIndex_[index] = 0;
182 <<
"There is no lookup table and the specie" <<
nl
184 <<
" is not found " <<
nl
232 const label
n = iter();
234 if (specieIndex_[
n] != 0)
239 const List<scalar>& Ynft = lookUpTablePtr_().lookUp(ft[celli]);
242 Xipi = Ynft[specieIndex_[
n]]*
paToAtm(
p[celli]);
252 const label index =
mixture.species().
find(iter.key());
260 scalar Ti =
T[celli];
263 coeffs_[bandi][
n].coeffs(
T[celli]);
265 if (coeffs_[bandi][
n].invTemp())
273 ((((
b[5]*Ti +
b[4])*Ti +
b[3])*Ti +
b[2])*Ti +
b[1])*Ti
318 E.
ref().primitiveFieldRef() =
320 *
Qdot.primitiveField()
321 *(iBands_[bandi][1] - iBands_[bandi][0])
327 E.
ref().primitiveFieldRef() =
329 *
Qdot.primitiveField()
330 *(iBands_[bandi][1] - iBands_[bandi][0])
336 <<
"Incompatible dimensions for Qdot field" <<
endl;
352 for (label j=0; j<nBands_; j++)
354 aLambda[j].primitiveFieldRef() = this->a(j);
357 aLambda[j].primitiveField()
358 *(iBands_[j][1] - iBands_[j][0])
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
bool found(const Key &key) const
Return true if hashed entry is found in table.
bool insert(const Key &key, const T &obj)
Copy insert a new entry, not overwriting existing entries.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
A list of pointers to objects of type <T>, with allocation/deallocation management of the pointers....
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
const word & constant() const
Return constant name.
label find(const T &val, label pos=0) const
Find index of the first occurrence of the value.
Specialization of basicMultiComponentMixture for a mixture consisting of a number for molecular speci...
Abstract base-class for fluid and solid thermodynamic properties.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
T get(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
const dictionary & optionalSubDict(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a sub-dictionary, otherwise return this dictionary.
const dictionary & subDict(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a sub-dictionary.
bool found(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Search for an entry (const access) with the given keyword.
bool readEntry(const word &keyword, T &val, enum keyType::option matchOpt=keyType::REGEX, bool mandatory=true) const
A keyword and a list of tokens is an 'entry'.
A class for handling file names.
Fundamental fluid thermodynamic properties.
Mesh data needed to do the Finite Volume discretisation.
const Time & time() const
Return the top-level database.
A list of lists. Interpolates based on the first dimension. The values must be positive and monotonic...
bool foundObject(const word &name, const bool recursive=false) const
Is the named Type found?
Type * getObjectPtr(const word &name, const bool recursive=false) const
Model to supply absorption and emission coefficients for radiation modelling.
const fvMesh & mesh() const
Reference to the mesh.
const dictionary & dict() const
Reference to the dictionary.
wideBandAbsorptionEmission radiation absorption and emission coefficients for continuous phase.
FixedList< FixedList< absorptionCoeffs, nSpecies_ >, maxBands_ > coeffs_
Absorption coefficients.
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.
virtual ~wideBandAbsorptionEmission()
Destructor.
tmp< volScalarField > eCont(const label bandi=0) const
Emission coefficient for continuous phase.
A class for managing temporary objects.
A class for handling words, derived from Foam::string.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
const word dictName("faMeshDefinition")
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
#define WarningInFunction
Report a warning using Foam::Warning.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionSet dimless
Dimensionless.
const dimensionSet dimEnergy
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVolume(pow3(dimLength))
errorManipArg< error, int > exit(error &err, const int errNo=1)
constexpr scalar paToAtm(const scalar pa) noexcept
Conversion from Pa to atm.
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.
#define forAllConstIters(container, iter)
Iterate across all elements of the container object with const access.
Unit conversion functions.