A method that includes the use of modified equations of state has been developed to enable noniterative, efficient computation of the thermodynamic behaviors of gases and liquids at temperatures >100 K and pressures from 1 to 100 MPa. The method is intended to be particularly useful for calculating such quantities as partial molar volumes, enthalpies, entropies, compressibilities, and thermal expansivities of oxidant/fuel mixtures at given pressures, temperatures, and mass fractions in gas turbine and rocket engines.

The method is based partly on the concept of a reference state of a real gas, which state is similar to that of a perfect gas at the same temperature and pressure and is defined with respect to the condition of the real gas at a relatively low reference pressure. The method is also based partly on the use of departure functions to represent the deviations of thermodynamic functions of the real gas at high pressure from reference-state values of those functions.

The form of the departure functions is obtained from the Peng-Robinson equation of state:

p=RT/(v-bm)-am/(v²+2bmv-b²m),

where p is the absolute pressure, R is the molar ideal-gas constant,T is the absolute temperature, and v is the molar volume. The parameters am and bm are semiempirical terms calculated from critical-state properties; these parameters follow conventional mixing rules; namely,

am=ΣΣxixjaij and bmxibi,

where xi is the mole fraction of the ith molecular species.

One first applies the departure-function formalism to each of the pure constituents of a mixture, then reuses this formalism for the mixture as a whole. In the application to each pure constituent, one expresses the reference-state enthalpy and entropy as finite series of temperature- and pressure-dependent terms, the coefficients of which are established by least-squares fits to the best available empirical data on enthalpy and entropy (see table). This method follows the standard practice in calculating the properties of nonideal mixtures by use of accurately known properties of the pure constituents along with excess Gibbs energy and/or fugacity coefficients, which are defined by use of conventional mixing rules. 