Fugacity measurements

The function / is called the fugacity of the gas. The fugacity measures the Gibbs energy of a real gas in the same way as the pressure measures the Gibbs energy of an ideal gas. 

Consider a binary mixture of two liquids that mix in all proportions. We assume that only component A is appreciably volatile. By measuring the fugacity of A in a gas phase equilibrated with the binary mixture, we can evaluate its activity coefficient based on a pure-liquid reference state ya = /a/( a/a) (Table 9.4). We wish to use the same fugacity measurements to determine the activity coefficient of the nonvolatile component, B. 

In Chapter 2 we discuss briefly the thermodynamic functions whereby the abstract fugacities are related to the measurable, real quantities temperature, pressure, and composition. This formulation is then given more completely in Chapters 3 and 4, which present detailed material on vapor-phase and liquid-phase fugacities, respectively. 

A component in a vapor mixture exhibits nonideal behavior as a result of molecular interactions only when these interactions are very wea)c or very infrequent is ideal behavior approached. The fugacity coefficient (fi is a measure of nonideality and a departure of < ) from unity is a measure of the extent to which a molecule i interacts with its neighbors. The fugacity coefficient depends on pressure, temperature, and vapor composition this dependence, in the moderate pressure region covered by the truncated virial equation, is usually as follows ... 

We have repeatedly observed that the slowly converging variables in liquid-liquid calculations following the isothermal flash procedure are the mole fractions of the two solvent components in the conjugate liquid phases. In addition, we have found that the mole fractions of these components, as well as those of the other components, follow roughly linear relationships with certain measures of deviation from equilibrium, such as the differences in component activities (or fugacities) in the extract and the raffinate. 

The determination of equilibria is done theoretically via the calculation of free energies. In practice, the concept of fugacity is used for which the unit of measurement is the bar. The equation linking the fugacity to the free energy is written as follows >... 

The properties of the solids most commonly encountered are tabulated. An important problem arises for petroleum fractions because data for the freezing point and enthalpy of fusion are very scarce. The MEK (methyl ethyl ketone) process utilizes the solvent s property that increases the partial fugacity of the paraffins in the liquid phase and thus favors their crystallization. The calculations for crystallization are sensitive and it is usually necessary to revert to experimental measurement. 

Only those components which are gases contribute to powers of RT. More fundamentally, the equiUbrium constant should be defined only after standard states are specified, the factors in the equiUbrium constant should be ratios of concentrations or pressures to those of the standard states, the equiUbrium constant should be dimensionless, and all references to pressures or concentrations should really be references to fugacities or activities. Eor reactions involving moderately concentrated ionic species (>1 mM) or moderately large molecules at high pressures (- 1—10 MPa), the activity and fugacity corrections become important in those instances, kineticists do use the proper relations. In some other situations, eg, reactions on a surface, measures of chemical activity must be introduced. Such cases may often be treated by straightforward modifications of the basic approach covered herein. 

The residual Gibbs energy and the fugacity coefficient are useful where experimental PVT data can be adequately correlated by equations of state. Indeed, if convenient treatment or all fluids by means of equations of state were possible, the thermodynamic-property relations already presented would suffice. However, liquid solutions are often more easily dealt with through properties that measure their deviations from ideal solution behavior, not from ideal gas behavior. Thus, the mathematical formahsm of excess properties is analogous to that of the residual properties. 

Solute/Solvent Systems The gamma/phi approach to X T.E calculations presumes knowledge of the vapor pressure of each species at the temperature of interest. For certain binary systems species I, designated the solute, is either unstable at the system temperature or is supercritical (T > L). Its vapor pressure cannot be measured, and its fugacity as a pure liquid at the system temperature/i cannot be calculated by Eq. (4-281). 

In Chapter 5, we considered systems in which composition becomes a variable, and defined and described the chemical potential. We showed that the chemical potential provides the condition for spontaneity or equilibrium. It is the potential that drives the flow of mass in a chemical process, A useful quantity related to the chemical potential is the fugacity. It can also be thought of as a measure of the flow of mass in a chemical process, and can be used to determine the point of equilibrium. It is often known as the escaping tendency since it can be used to describe the ease with which mass flows from one phase to another, particularly the flow from a solid or liquid phase to a gas phase. 

Binary (vapor + liquid) equilibria studies involve the determination of / as a function of composition. the mole fraction in the liquid phase. Of special interest is the dependence of/ on composition in the limit of infinite dilution. In the examples which follow, equilibrium vapor pressures, p,. are measured and described. These vapor pressures can be corrected to vapor fugacities using the techniques described in the previous section. As stated earlier, at the low pressures involved in most experiments, the difference between p, and / is very small, and we will ignore it unless a specific application requires a differentiation between the two. 

The vapor pressures (fugacities) shown in Figure 6.14 were reported by J. J. Fritz and C. R. Fuget, Vapor Pressure of Aqueous Hydrogen Chloride Solutions", Chem. Eng. Data Ser., 1, 10-12 (1956). The vapor pressures are too small to measure directly. The values reported were calculated from the results of emf measurements made on an electrochemical cell. In Chapter 9, we will describe this and other cells in detail. 

The vapor pressure of a polymer is, of course, far too small to measure We may, nevertheless, insist that such a vapor pressure exists, however small it may be. Or we may resort to the use of the escaping tendency, or fugacity, in place of the partial vapor pressure in the development given above, in accordance with usual thermodynamic procedures applied to the treatment of solutions. The treatment given here is in no way restricted to volatile solutes. 

The other state variables are the fugacity of dissolved methane in the bulk of the liquid water phase (fb) and the zero, first and second moment of the particle size distribution (p0, Pi, l )- The initial value for the fugacity, fb° is equal to the three phase equilibrium fugacity feq. The initial number of particles, p , or nuclei initially formed was calculated from a mass balance of the amount of gas consumed at the turbidity point. The explanation of the other variables and parameters as well as the initial conditions are described in detail in the reference. The equations are given to illustrate the nature of this parameter estimation problem with five ODEs, one kinetic parameter (K ) and only one measured state variable. 

Some economies are possible if equilibrium is assumed between selected compartments, an equal fugacity being assignable. This is possible if the time for equilibration is short compared to the time constant for the dominant processes of reaction or advection. For example, the rate of chemical uptake by fish from water can often be ignored (and thus need not be measured or known within limits) if the chemical has a life time of hundreds of days since the uptake time is usually only a few days. This is equivalent to the frequently used "steady state" assumption in chemical kinetics in which the differential equation for a short lived intermediate species is set to zero, thus reducing the equation to algebraic form. When the compartment contains a small amount of chemical or adjusts quickly to its environment, it can be treated algebraically. 

An attractive feature of K<)A is that it can replace the liquid or supercooled liquid vapor pressure in a correlation. K,-ja is an experimentally measurable or accessible quantity, whereas the supercooled liquid vapor pressure must be estimated from the solid vapor pressure, the melting point and the entropy of fusion. The use of KOA thus avoids the potentially erroneous estimation of the fugacity ratio, i.e., the ratio of solid and liquid vapor pressures. This is especially important for solutes with high melting points and, thus, low fugacity ratios. 

In the multimedia models used in this series of volumes, an air-water partition coefficient KAW or Henry s law constant (H) is required and is calculated from the ratio of the pure substance vapor pressure and aqueous solubility. This method is widely used for hydrophobic chemicals but is inappropriate for water-miscible chemicals for which no solubility can be measured. Examples are the lower alcohols, acids, amines and ketones. There are reported calculated or pseudo-solubilities that have been derived from QSPR correlations with molecular descriptors for alcohols, aldehydes and amines (by Leahy 1986 Kamlet et al. 1987, 1988 and Nirmalakhandan and Speece 1988a,b). The obvious option is to input the H or KAW directly. If the chemical s activity coefficient y in water is known, then H can be estimated as vwyP[>where vw is the molar volume of water and Pf is the liquid vapor pressure. Since H can be regarded as P[IC[, where Cjs is the solubility, it is apparent that (l/vwy) is a pseudo-solubility. Correlations and measurements of y are available in the physical-chemical literature. For example, if y is 5.0, the pseudo-solubility is 11100 mol/m3 since the molar volume of water vw is 18 x 10-6 m3/mol or 18 cm3/mol. Chemicals with y less than about 20 are usually miscible in water. If the liquid vapor pressure in this case is 1000 Pa, H will be 1000/11100 or 0.090 Pa m3/mol and KAW will be H/RT or 3.6 x 10 5 at 25°C. Alternatively, if H or KAW is known, C[ can be calculated. It is possible to apply existing models to hydrophilic chemicals if this pseudo-solubility is calculated from the activity coefficient or from a known H (i.e., Cjs, P[/H or P[ or KAW RT). This approach is used here. In the fugacity model illustrations all pseudo-solubilities are so designated and should not be regarded as real, experimentally accessible quantities. 

Saturation properties such as solubility in water and vapor pressure can be measured directly for solids and liquids. For certain purposes it is useful to estimate the solubility that a solid substance would have if it were liquid at a temperature below the melting point. For example, naphthalene melts at 80°C and at 25°C the solid has a solubility in water of 33 g/m3 and a vapor pressure of 10.9 Pa. If naphthalene was a liquid at 25°C it is estimated that its solubility would be 115 g/m3 and its vapor pressure 38.1 Pa, both a factor of 3.5 greater. This ratio of solid to liquid solubilities or vapor pressures is referred to as the fugacity ratio. It is 1.0 at the melting point and falls, in this case at lower temperatures to 0.286 at 25°C. 

Measurements of gas chromatographic retention time are often used as a fast and easy method of estimating vapor pressure. These estimated pressures are related to the gas/substrate partition coefficient, which can be regarded as a ratio of solubility of the substance in the gas to that in the substrate, both solubilities being of the substance in the liquid state. As a result the estimated vapor pressures are of the liquid state. To obtain the solid vapor pressure requires multiplication by the fugacity ratio. It is important to establish if the estimated and reported property is of the vapor or liquid. 

As was discussed earlier in Section 1.2.8 a complication arises in that two of these properties (solubility and vapor pressure) are dependent on whether the solute is in the liquid or solid state. Solid solutes have lower solubilities and vapor pressures than they would have if they had been liquids. The ratio of the (actual) solid to the (hypothetical supercooled) liquid solubility or vapor pressure is termed the fugacity ratio F and can be estimated from the melting point and the entropy of fusion. This correction eliminates the effect of melting point, which depends on the stability of the solid crystalline phase, which in turn is a function of molecular symmetry and other factors. For solid solutes, the correct property to plot is the calculated or extrapolated supercooled liquid solubility. This is calculated in this handbook using where possible a measured entropy of fusion, or in the absence of such data the Walden s Rule relationship suggested by Yalkowsky (1979) which implies an entropy of fusion of 56 J/mol-K or 13.5 cal/mol-K (e.u.)... 

Fig. 6.1. pH of surface seawater from the western Pacific Ocean (Skirrow, 1965), as measured in situ during oceanographic cruises (various symbols). Line shows pH predicted by the model for seawater in equilibrium with atmospheric CO2 at a fugacity of 10-3-5. Dashed lines show pH values that result from assuming larger fugacities of 10-33 and... 

The log K for this reaction increases from —2.12 at 250°C to —1.45 at 25 °C. The final C02 fugacity is about 15, corresponding to a partial pressure considerably in excess of atmospheric pressure. We would certainly need to take extraordinary measures to prevent the fluid from effervescing, if we were actually performing this experiment instead of simulating it. 

Artola-Garicano et al. [27] compared their measured removals of AHTN and HHCB [24] to the predicted removal of these compounds by the wastewater treatment plant model Simple Treat 3.0. Simple Treat is a fugacity-based, nine-box model that breaks the treatment plant process into influent, primary settler, primary sludge, aeration tank, solid/liquid separator, effluent, and waste sludge and is a steady-state, nonequilibrium model [27]. The model inputs include information on the emission scenario of the FM, FM physical-chemical properties, and FM biodegradation rate in activated sludge. 

---