See also CHEMKIN lecture
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.
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:
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.
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:
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'.
See PHENC entry: Turbulence Models in PHOENICS