The "seven-gases" option of PHOENICS provides a means of computing the density and temperature of gases involved in single- or two-phase combustion processes; it is useful because it is mainly by way of density and temperature changes that such processes influence velocities, pressure and turbulence.
However, it has a wider significance, as will now be explained.
The seven-gases concept represents an idealization of practically-arising combustion situations, achieved by restricting the number of substances present, and their possible interactions.
Nitrogen (N) is also present, as an inert element; and others, such as silicon (Si) which may be present in ash, are also treated as not engaging in chemical reactions.
The first six may all engage in chemical reactions, i.e. may exchange chemical elements between them.
Nitrogen too may react chemically, as is well-known to be of practical importance. However the extent of the reaction is usually small enough to be neglected in its effects on gas temperature and density.
The seven-gases option can be used in both single-phase (ONEPHS=T) and two-phase (ONEPHS=F) flow simulations. However, more options are available than are implied by these words; for "PHOENICS phases" and "thermodynamicist's phases" are not identical, as will now be explained.
Examples of the choices which could be made are:-
This is indeed the choice which is represented in the coding supplied in the GXDENS.FOR file described below.
It is then possible to deduce the composition of the gas from knowledge of the mass fractions of chemical elements: FO, FC, FH and FN. The following diagram is of assistance.
Figure 1 shows an equilibrium-composition diagram which represents the six reactive gases (nitrogen being therefore excluded), in respect of their elemental oxygen, carbon and hydrogen contents.
It shows three main regions, namely:
This diagram has been drawn on the presumption that oxidation of any of the fuels proceeds in two stages, viz:
Any particular mixture can be located as a point (say M) on the diagram; and its mass fractions of the three elements, FO, FC and FH, are proportional to the perpendicular distances from the point to the side of the triangle opposite to the vertex associated with the element in question.
The sum (FO + FC + FH) would equal unity, were other materials absent;
but in general it equals:
(1 - FN - F_condensed_phase_materials)
the final term being present, of course, only if these materials are being treated as belonging to PHOENICS-phase 1 because they "flow with the gases".
The within-gas chemical reactions do not give rise to sources, because, though they transfer chemical elements from one mixture component to another, they do not change the mass of any element.
Among the interesting consequences are:
When there are n sources of supply of gas-phase material, it is best to solve for n-1 conservation equations, for variables which might be named: MXF1, MXF2, MXF3; and then to deduce the elemental mass fractions of carbon, oxygen and hydrogen from:
FC = sum ( FC1*MXF1 + FC2*MXF2 + ........ + FCn*MXFn )
However they are computed, it is the values of FC, FO, and FH which dictate the relative proportions of O2, CO2, CO, H2O, H2, CxHy (treated as one material) and FN in the gaseous equilibrium mixture.
Species mass fractions begin with:
The diffusion coefficients of the gaseous species are all taken as equal, as are the specific heats; and the reaction rates are diffusion-limited. As a consequence, all species concentrations depend, in piecewise-linear fashion, on the elemental mass fractions.
The values of oxygen fraction FO at which the formulae exbibit discontinuities of slope are called:
T = ( H - HCHX*YCHX - HCOCO2*YCO - HHH2O*YH2 - HCHAR*YCHAR - HVOL*YVOL ) / CPTemperature and concentration are computed as they are needed; they are put into storage arrays if these have been created by instructions in the satellite, as evidenced:
The Fortran coding which computes the gas density according to the seven-gases presumption is in gxdens.for, which can be seen by clicking here.
The numerical values which appear in it are:-
If one or more of the chemical reactions needed to bring about the equilibrium composition represented by Fig. 1 can not reasonably be presumed to proceed with sufficient rapidity, it remains useful to compute the elemental mass fractions; but it is not sufficient.
For example, it might be supposed that the reaction:
2*CO + O2 -> 2*CO2
is the one non-fast reaction.
Then, whereas it might still be presumed that H2 and O2 could not coexist in the red and green regions of Fig.1, CO and O2 could do so.
What values would their concentrations attain? The question could be answered only by solving an additional differential equation, which would have to be provided with a source term representing the finite rate of the CO-oxidation reaction.
It would be immaterial whether the variable solved for were the concentration of CO, of O2, or of CO2; for each can be computed, if one of them is known, from the known values of FO and FC.
Of course, as is well-known, additional gas components such as the "radicles" [O] and [OH], play essential roles in the CO-combustion process; but these are usually such low concentration as not appreciably to affect the density or temperature, which is what the 7-gases option is designed to do.
The subject of how finite chemical reaction rates are to be computed is too extensive to be entered upon here. Attention will be switched instead to the more-frequently encountered reason for departures from equilibrium, namely turbulence.
The time-average mixture in a combustion chamber is often NOT in equilibrium for a reason which is explicable in terms, not of finite-rate chemistry, but of the presence of high-frequency fluctuations of concentration and temperature, i.e. of turbulence.
Precise description of this phenomenon is still beyond the scope of science; but useful idealizations exist, of which one of the most commonly employed is that the gases at any point within a combustion chamber act like a random intermingling of two gases having the same elemental composition, but of which:
It is this concept which underlay the "eddy-break-up (EBU) model", first proposed in 1971. Although EBU has been superseded by the "multi-fluid" (MFM) model, it still allows the major effects of turbulence-chemistry interactions to be reaistically represented.
In MFM terms, EBU represents:
Although the details differ, and the EBU formulation is the more direct, both involve the solution of at least one additional equation transport equation having a source term proportional to the above-mentioned product, from which the mass fractions of the two fluids can be deduced.
From these, and from the densities and temperatures of the individual fluids, the mixture-mean densities and temperatures can be deduced.
In relation to the 7-gases model of density and temperature, it is clear that the composition of the first of these fluids is to be calculated as indicated above under the heading "equilibrium".
And the second? Its composition is deducible for a two-supply-stream situation from the local value of MIXF; and specifically the concentration of each species is equal to:
its concentration in the fuel-supply stream * MIXF
its concentration in the oxidant-supply stream * (1 - MIXF) .
Once FCMX, FHMX, FOMX and FNMX have been computed in this way, the compositions of the two components of the fluid population can be computed.
concentration of species A =
MXF1 * concentration of A in supply stream 1 + MXF2 * concentration of A in supply stream 2 + ..... MXFn * concentration of A in supply stream n
Input-file library cases which exemplify use of the 7-gases model can be seen by using the PHOENICS Commander's library-search facility
Of these, case 477 is a PIL "macro", loaded into other cases by the load(477) command, or simply #477.
Case 477 itself call another macro, namely case 478, which is where, among other things, the non-PIL variables
In case 477, they are assigned to the PIL variables RHO1A, RHO1B and RHO1C for transmission to Earth via EARDAT.