12.1.16. Thermochemical Properties

MFiX uses the Burcat or NASA 7-Coefficient polynomials to calculate thermochemical properties of species (not the NASA 9-Coefficient polynomials). The directory model/thermochemical contains the database of Burcat and Ruscic (2005) and routines for reading the database. With linkage to this database the users need not manually enter data for molecular weight, specific heat, and heats of reactions. Instead the users only need to enter the names of the species (keyword SPECIES_g and SPECIES_s) in the data file. If such information is already provided in either the data file or a BURCAT.THR file in the run directory then MFiX will not reference the database. That is, MFiX reads the necessary thermochemical data from files in the following order:

  1. mfix.dat or *.mfx

  2. BURCAT.THR file in the run directory

  3. model/thermochemical/BURCAT.THR

The species names are case sensitive and should match the names in BURCAT.THR exactly; alternatively aliases can be defined for common species, such as O2, in read_therm.f. See tests/thermo for a sample case that accesses the database. The format of BURCAT.THR file resembles CHEMKIN format, but with several notable differences. To include thermochemical data in the mfix.dat file then this information must start below a line that starts with THERMO DATA.

Example dataset from BURCAT.THR with notations:

../_images/burcat.png

Each entry in the database starts with a unique CAS identifier (74-82-8) for the species, followed by several lines of comments highlighted in green. The data section starts with the species name in columns 1-18 (CH4 RRHO). Common species names may be followed by strings (RRHO) that identify the method used to determine the coefficients. Additional information follows the species name. The numbers toward the end of the line are the temperature limits (200.000 6000.000) in degrees Kelvin where the property calculation is valid and the molecular weight (16.04246). Unlike CHEMKIN the common temperature for the high and low temperature branches are not recorded; it is always 1000 K. The next three lines give the fourteen coefficients (seven coefficients each for the high and low temperature branches) and the formation enthalpy at 298 K (which is also not included in CHEMKIN format). All the coefficients and the enthalpy of formation are normalized with the gas constant R (cal/mol/K). The low temperature coefficients (\(a_L\)) should be used for temperatures in the range Tlow to 1000K and the high temperature coefficients (\(a_H\)) should be used for temperatures in the range 1000K to Thigh. The coefficients are stored in a fixed format (E15.0) as follows:

Table 12.94 Coefficients

\(a^1_H\)

\(a^2_H\)

\(a^3_H\)

\(a^4_H\)

\(a^5_H\)

\(a^6_H\)

\(a^7_H\)

\(a^1_L\)

\(a^2_L\)

\(a^3_L\)

\(a^4_L\)

\(a^5_L\)

\(a^6_L\)

\(a^7_L\)

\(\frac{\Delta H^{\circ}_f}{R}\)

where \(\Delta H^{\circ}_f\) is the formation enthalpy at 298K.

The normalized specific heat is given by

\[\frac{C_p}{R} = a_1 + a_2 T + a_3 T^2 + a_4 T^3 + a_5 T^4 .\]

The coefficients \(a_6\) and \(a_7\) are not required to calculate \(C_p\). They are used in the \(H_T\), \(S_T\), and \(G_T\) polynomials which are not used by MFiX:

\[ \begin{align}\begin{aligned}\frac{H_T}{RT} = a_1 + \frac{a_2 T}{2} + \frac{a_3 T^2}{3} + \frac{a_4 T^3}{4} + \frac{a_5 T^4}{5} + \frac{a_6}{T}\\\frac{S_T}{R} = a_1 \ln{T} + a_2 T + \frac{a_3 T^2}{2} + \frac{a_4 T^3}{3} + \frac{a_5 T^4}{4} + a_7\\\frac{G_T}{RT} = \frac{H_T}{RT} - \frac{S_T}{R}\end{aligned}\end{align} \]

Instead, MFiX integrates \(C_p\) to determine changes in enthalpy. See Burcat Introduction document for more details, which is also included in the source, model/thermochemical/intro.pdf.

Additional database comments:

  • A number of species in the database have a lower temperature limit of 300K which is 2 degrees above the reference temperature (298 K) used for formation enthalpy calculation. For those species MFiX relaxes the lower limit for \(C_p\) calculations to 298 K to enable heat of reaction calculation (see read_database.f).

  • The database reader is set up such that the database is read only if necessary.

  • If you include thermochemical properties in mfix.dat all keywords defined below the line that starts with THERMO DATA will be ignored!