Encyclopaedia Index


See also CHEMKIN lecture


  1. Introduction
  2. The transport model
  3. Physical properties
  4. Chemical reaction models
  5. Further problem-specification facilities
  6. Additional print-out
  7. Model set-up
  8. Exemplification
  9. Developments and modifications
  10. Sources of further information

March 1996

6. Additional Print-out

It is useful, indeed it is often the reason to perform a simulation, to be able to print-out and to plot reaction rates, rates of production of species, the heat-release rate, and the elemental composition. The print-out is activated as follows:

STORE(QDOT) - activates storage and print-out of heat- release rate in erg/s/cm^3.
STORE(Cn+) - activates storage and print-out of the rate of production of concentration Cn, eg. STORE(C12+) for rate of production of C12
STORE(n&) - activates storage and print-out of the rate of reaction n, eg. STORE(4&) for rate of reaction 4
STORE(ELxx) - activates storage and print-out of the mass composition of element xx, eg. STORE(ELN) for nitrogen

In addition, the elemental composition is printed at the monitor point whenever the monitor data is printed.

The elemental composition is calculated as follows:

Ej = sum (Ajk*Wj*Yk/Wk)

where the summation is over all species k, Ej is the mass composition of element j, Aij is the number of atoms of the element j in the kth species, Yk is the mass fraction of species k, Wj is the molecular mass of element k, and Wk is the molecular mass of species k.

Mole-fractions are computed and printed at completion, and whenever the NPRINT settings etc. indicate that the STOREd dependent variables are to be printed. The mole-fraction print-out is always generated, and does not require storage of extra variables.

However, if desired, the user may STORE the mole fractions for printing in the RESULT file and also for plotting purposes with PHOTON and AUTOPLOT. This facility is activated by setting:


which activates storage and print-out of the mole fraction of the species xx, eg. STORE(MH2) for the SOLVEd hydrogen species H2.

The mole fractions Xi are calculated from:

Xi = Yi*W/Wi

where W is the mixture molecular mass, and Wi is the molecular mass of species i.

7 Model Setup

In order to ease the setup of models employing the CHEMKIN interface, a SATLIT code, CHEMST, and two PIL fragments have been supplied. Prior to calling CHEMST the user should, where appropriate:

When called, the routine CHEMST does the following:

9 Developments and modifications


9.1 Developments
9.2 Modifications Made to CHEMKIN Routines

9.1 Developments

There are a number of areas in which work remains to be done, and these are:-

Much more can be done in many areas, but what should be done depends upon the demand for further deveopments.

9.2 Modifications Made to CHEMKIN Routines

No modifications have been made to any subroutines in the CHEMKIN or TRANLIB libraries. The only software associated with CHEMKIN to be have been modified is the TWOPNT solver. The modifications made are:

10. Sources of Further Information

Further information can be found in the references cited below, and a lecture is provided under the POLIS entries: 'Lectures on PHOENICS', 'Turbulence and Combustion'. Additional information on running CHEMKIN can be found under the POLIS entries 'About PHOENICS 2.2', 'How to run PHOENICS modules'.

  1. R.J.Kee, F.M.Rupley, J.A.Miller, 'CHEMKIN-II: A FORTRAN Chemical Kinetics Package for the Analysis of Gas-Phase Chemical Kinetics', SAND89-8009B, UC-706, Sandia National Laboratories, Albuquerque, New Mexico, (1993a).
  2. R.J.Kee, G.Dixon-Lewis, J.Warnatz, M.E.Coltrin, and J.A.Miller, 'A FORTRAN Computer Code for the Evaluation of Gas-Phase Multicomponent Transport Properties", Sandia Report SAND86-8246, UC-401, (1993b).
  3. R.J.Kee, F.M.Rupley and J.A.Miller, 'The CHEMKIN Thermodynamic Data Base', Sandia Report SAND87-8215B, UC-4, (1993c).
  4. P.Glarborg, R.J.Kee, J.F.Grear, J.A.Miller,'PSR: A FORTRAN Program for Modelling Well-Stirred Reactors', SAND86-8209, Sandia National Laboratories, Albuquerque, New Mexico, (1992).
  5. J.F.Grear, 'The Twopnt Program for Boundary Value Problems', Sandia Report SAND91-8230, (1992)