chemistryTabulationMethodNew.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2016-2018 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
27 #include "Time.H"
28 
29 // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * //
30 
31 template<class CompType, class ThermoType>
34 (
35  const IOdictionary& dict,
37 )
38 {
39  const dictionary& tabulationDict(dict.subDict("tabulation"));
40 
41  const word methodName(tabulationDict.lookup("method"));
42 
43  Info<< "Selecting chemistry tabulation method " << methodName << endl;
44 
45  const word methodTypeName =
46  methodName + '<' + CompType::typeName + ',' + ThermoType::typeName()
47  + '>';
48 
49  typename dictionaryConstructorTable::iterator cstrIter =
50  dictionaryConstructorTablePtr_->find(methodTypeName);
51 
52  if (cstrIter == dictionaryConstructorTablePtr_->end())
53  {
55  << "Unknown " << typeName_() << " type " << methodName << endl
56  << endl;
57 
58  const wordList names(dictionaryConstructorTablePtr_->toc());
59 
60  wordList thisCmpts;
61  thisCmpts.append(word::null);
62  thisCmpts.append(CompType::typeName);
63  thisCmpts.append
64  (
65  basicThermo::splitThermoName(ThermoType::typeName(), 5)
66  );
67 
68  wordList validNames;
69  forAll(names, i)
70  {
71  const wordList cmpts(basicThermo::splitThermoName(names[i], 7));
72 
73  bool isValid = true;
74  for (label i = 1; i < cmpts.size() && isValid; ++ i)
75  {
76  isValid = isValid && cmpts[i] == thisCmpts[i];
77  }
78 
79  if (isValid)
80  {
81  validNames.append(cmpts[0]);
82  }
83  }
84 
86  << "Valid " << typeName_() << " types for this thermodynamic model "
87  << "are:" << endl << validNames << endl;
88 
89  List<wordList> validCmpts;
90  validCmpts.append(wordList(7, word::null));
91  validCmpts[0][0] = typeName_();
92  validCmpts[0][1] = "reactionThermo";
93  validCmpts[0][2] = "transport";
94  validCmpts[0][3] = "thermo";
95  validCmpts[0][4] = "equationOfState";
96  validCmpts[0][5] = "specie";
97  validCmpts[0][6] = "energy";
98  forAll(names, i)
99  {
100  validCmpts.append(basicThermo::splitThermoName(names[i], 7));
101  }
102 
104  << "All " << validCmpts[0][0] << '/' << validCmpts[0][1]
105  << "/thermoPhysics combinations are:" << endl << endl;
106  printTable(validCmpts, FatalErrorInFunction);
107 
109  }
110 
112  (
113  cstrIter()(dict, chemistry)
114  );
115 
116 }
117 
118 
119 // ************************************************************************* //
dictionary dict
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:428
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
error FatalError
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Definition: error.H:319
void size(const label)
Override size to be inconsistent with allocated storage.
Definition: ListI.H:163
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:256
BasicChemistryModel< rhoReactionThermo > & chemistry
IOdictionary is derived from dictionary and IOobject to give the dictionary automatic IO functionalit...
Definition: IOdictionary.H:52
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
Definition: dictionary.C:692
A class for handling words, derived from string.
Definition: word.H:59
void append(const T &)
Append an element at the end of the list.
Definition: ListI.H:184
void printTable(const List< wordList > &, List< string::size_type > &, Ostream &)
Definition: wordIOList.C:42
List< word > wordList
A List of words.
Definition: fileName.H:54
messageStream Info
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:52
static autoPtr< chemistryTabulationMethod > New(const IOdictionary &dict, TDACChemistryModel< CompType, ThermoType > &chemistry)