research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fixed a bug in reporting parameters for Nasa9MultiTempRegion

Browse source 
Changed the Nasa9 reporting parameters routine to be compatible with
the multi temperature region reporting parameters.
This commit is contained in:
Harry Moffat 2007-09-14 16:05:42 +00:00
parent e961587418
commit a8fe60d601
4 changed files with 19 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
Cantera/src/thermo/Nasa9Poly1.cpp
View file

@ -276,8 +276,11 @@ namespace Cantera {
tlow = m_lowT;
thigh = m_highT;
pref = m_Pref;
coeffs[0] = 1;
coeffs[1] = m_lowT;
coeffs[2] = m_highT;
for (int i = 0; i < 9; i++) {
coeffs[i] = m_coeff[i];
coeffs[i+3] = m_coeff[i];
}
}

8
Cantera/src/thermo/Nasa9Poly1.h
View file

@ -196,7 +196,13 @@ namespace Cantera {
* @param thigh output - Maximum temperature
* @param pref output - reference pressure (Pa).
* @param coeffs Vector of coefficients used to set the
* parameters for the standard state.
* parameters for the standard state. There are
* 12 of them, designed to be compatible
* with the multiple temperature formulation.
* coeffs[0] is equal to one.
* coeffs[1] is min temperature
* coeffs[2] is max temperature
* coeffs[3+i] from i =0,9 are the coefficients themselves
*/
virtual void reportParameters(int &n, int &type,
doublereal &tlow, doublereal &thigh,

2
Cantera/src/thermo/Nasa9PolyMultiTempRegion.cpp
View file

@ -343,7 +343,7 @@ namespace Cantera {
doublereal &pref,
doublereal* const coeffs) const {
n = m_index;
type = NASA9;
type = NASA9MULTITEMP;
tlow = m_lowT;
thigh = m_highT;
pref = m_Pref;

7
Cantera/src/thermo/Nasa9PolyMultiTempRegion.h
View file

@ -160,6 +160,13 @@ namespace Cantera {
* @param pref output - reference pressure (Pa).
* @param coeffs Vector of coefficients used to set the
* parameters for the standard state.
* There are 1 + 11*nzones coefficients
* coeffs[0] is equal to nTempZones.
* index = 1
* for each zone:
* coeffs[index] = minTempZone
* coeffs[index+1] = maxTempZone
* coeffs[index+2+i] from i =0,9 are the coefficients themselves
*/
virtual void reportParameters(int &n, int &type,
doublereal &tlow, doublereal &thigh,