chemistry::SpeciesThermo Struct Reference

Species thermodynamics as piecewise polynomials in NASA-9 form, one per temperature range: cp/R = a0/T^2 + a1/T + a2 + a3 T + a4 T^2 + a5 T^3 + a6 T^4, h/(RT) = -a0/T^2 + a1 ln(T)/T + a2 + a3 T/2 + a4 T^2/3 + a5 T^3/4 + a6 T^4/5 + a7/T, s/R = -a0/(2 T^2) - a1/T + a2 ln(T) + a3 T + a4 T^2/2 + a5 T^3/3 + a6 T^4/4 + a8. More...

#include <mechanism.h>

Public Member Functions
size_t	range (double T) const
	Index of the range used at T. As in Cantera, a NASA-7 bound belongs to the range below it, a NASA-9 bound to the range above.

Public Attributes
ThermoModel	model = ThermoModel::NASA7
std::vector< double >	T_bounds
std::vector< std::array< double, 9 > >	coeffs
std::array< double, 4 >	constant_cp = {}

Detailed Description

Species thermodynamics as piecewise polynomials in NASA-9 form, one per temperature range: cp/R = a0/T^2 + a1/T + a2 + a3 T + a4 T^2 + a5 T^3 + a6 T^4, h/(RT) = -a0/T^2 + a1 ln(T)/T + a2 + a3 T/2 + a4 T^2/3 + a5 T^3/4 + a6 T^4/5 + a7/T, s/R = -a0/(2 T^2) - a1/T + a2 ln(T) + a3 T + a4 T^2/2 + a5 T^3/3 + a6 T^4/4 + a8.

NASA-7 coefficients (b0..b6) are stored as (0, 0, b0, ..., b6); constant cp as a2 = cp0/R, a7 = (h0 - cp0 T0)/R, a8 = (s0 - cp0 ln T0)/R. Outside the bounds the outermost polynomial is extrapolated.

Member Function Documentation

range()

size_t chemistry::SpeciesThermo::range

(

double

T

)

const

Index of the range used at T. As in Cantera, a NASA-7 bound belongs to the range below it, a NASA-9 bound to the range above.

Member Data Documentation

coeffs

std::vector<std::array<double, 9> > chemistry::SpeciesThermo::coeffs

constant_cp

std::array<double, 4> chemistry::SpeciesThermo::constant_cp = {}

model

ThermoModel chemistry::SpeciesThermo::model = ThermoModel::NASA7

T_bounds

std::vector<double> chemistry::SpeciesThermo::T_bounds

---

The documentation for this struct was generated from the following files:

• /home/runner/work/mallard/mallard/mallard/src/chemistry/mechanism.h
• /home/runner/work/mallard/mallard/mallard/src/chemistry/mechanism.cpp