2. Properties, Equations of State, and Phase Behavior

Introduction and Context

Thermodynamic analysis depends on the ability to determine the state of a substance and to evaluate its properties accurately. Pressure, temperature, specific volume, internal energy, enthalpy, and entropy are not merely tabulated quantities; they are the variables through which energy conversion, heat-transfer performance, fluid-flow behavior, and equipment design are predicted. In practical engineering, the quality of an analysis is often limited less by the governing energy equation than by the accuracy with which properties are obtained.

A pure substance may exist as a compressed liquid, saturated liquid, two-phase liquid–vapor mixture, saturated vapor, or superheated vapor. Its properties change sharply near phase boundaries and may depart substantially from ideal-gas behavior at high pressure or low temperature. Engineers therefore use several complementary tools: phase diagrams, property tables, equations of state, generalized compressibility correlations, corresponding-states methods, and computational property packages.

This chapter develops a systematic understanding of these tools. Particular attention is given to the meaning of thermodynamic properties, the geometry of phase behavior, real-gas models, vapor–liquid equilibrium, and the selection of an appropriate property method. The discussion provides the foundation for later analysis of cycles, heat exchangers, turbines, compressors, refrigeration systems, and reacting or multiphase equipment.

Learning Objectives

After studying this chapter, the reader should be able to:

1. Distinguish intensive, extensive, specific, and molar properties and explain the state postulate for a simple compressible pure substance.

2. Identify compressed-liquid, saturated, two-phase, and superheated states using phase diagrams and property data.

3. Determine quality and mixture properties in the liquid–vapor region.

4. Use pressure, temperature, and volume tables appropriately, including interpolation and reference-state conventions.

5. Explain the physical meaning of the critical point, triple point, metastable states, and phase equilibrium.

6. Apply cubic equations of state and compressibility-factor methods to real gases.

7. Interpret fugacity, activity, and vapor–liquid equilibrium at an introductory graduate level.

8. Select suitable property software and assess the reliability of calculated results.

9. Recognize common property-estimation errors and make technically defensible engineering decisions.

The Thermodynamic State and Its Properties

A thermodynamic system is described by measurable macroscopic variables. A property is a quantity whose value depends only on the state of the system, not on the path used to reach that state. Pressure, temperature, density, internal energy, enthalpy, and entropy are properties. Heat and work are not properties because their values depend on the process.

Properties are classified as intensive or extensive. Intensive properties, such as temperature, pressure, density, and specific volume, do not depend on the amount of matter. Extensive properties, such as mass, total volume, total internal energy, and total entropy, are proportional to system size. Dividing an extensive property by mass produces a specific property; dividing it by amount of substance produces a molar property. Thus,

v = V/m, u = U/m, h = H/m

where lower-case symbols denote specific quantities. For a fixed composition, specific and molar properties are intensive.

For a simple compressible pure substance at equilibrium, the state is fixed by two independent intensive properties. This is the practical form of the state postulate. Pressure and temperature are usually independent in a single-phase region, but they are not independent along a saturation curve. In the two-phase region, pressure and temperature are related by saturation conditions; one additional independent property, such as quality or specific volume, is required.

The specific enthalpy is defined by

h = u + pv.

Although enthalpy is not a directly measured energy stored in a fluid in the same sense as internal energy, it is especially convenient in steady-flow devices because the combination u + pv appears naturally in the energy equation. Entropy measures the dispersal of energy and provides the basis for the second-law analysis of real processes. Absolute values of u, h, and s are not usually required; differences are. Consequently, property tables adopt reference states, and different data sources can use different zero points without changing physically meaningful energy balances, provided the same source is used consistently.

Phase Behavior of a Pure Substance

A pure substance has uniform chemical composition throughout the system, although it may contain more than one phase. Water and refrigerant R134a are pure substances when each consists of only the stated chemical species, even when liquid and vapor coexist.

When a liquid is heated at a fixed pressure, its temperature rises until it reaches the saturation temperature. Further heat addition produces vaporization at nearly constant temperature and pressure. The reverse process occurs during condensation. The saturated liquid state is the state at which the first infinitesimal amount of vapor is about to form. The saturated vapor state is the state at which the last infinitesimal amount of liquid has disappeared.

Between these two states lies the saturated liquid–vapor mixture region. The quality, or dryness fraction, is defined as

x = m_v/(m_l + m_v),

where m_v is vapor mass and m_l is liquid mass. Quality is meaningful only in the two-phase region. It is not a general measure of vapor content in superheated or compressed states.

For any extensive or specific property that is approximately additive between phases,

Y = (1 − x)Y_f + xY_g,

or, equivalently,

Y = Y_f + xY_fg,

where Y may represent v, u, h, or s; Y_f is the saturated-liquid value, Y_g is the saturated-vapor value, and Y_fg = Y_g − Y_f. Therefore,

x = (Y − Y_f)/Y_fg.

For example, if a water mixture at a specified pressure has h_f = 640 kJ/kg, h_g = 2740 kJ/kg, and h = 1690 kJ/kg, its quality is approximately (1690 − 640)/(2740 − 640) = 0.50. This means that 50 percent of the mass is vapor, not that 50 percent of the volume is vapor. Because vapor has a much larger specific volume than liquid, its volume fraction may be close to unity even when its mass fraction is modest.

Phase Diagrams and Critical Behavior

A pressure–temperature diagram shows the regions in which solid, liquid, and vapor are stable. The sublimation curve separates solid and vapor, the fusion curve separates solid and liquid, and the vaporization curve separates liquid and vapor. These curves meet at the triple point, where all three phases coexist in equilibrium. The liquid–vapor curve terminates at the critical point, defined by the critical temperature and critical pressure. Above the critical temperature, a gas cannot be liquefied by pressure alone.

A pressure–specific-volume or temperature–specific-volume diagram provides additional information about the amount of each phase. The dome-shaped two-phase region is bounded by the saturated-liquid line and saturated-vapor line. The peak of the dome corresponds to the critical point. Within the dome, pressure and temperature are not independent for a pure substance, and quality varies from zero on the saturated-liquid boundary to one on the saturated-vapor boundary.

Near the critical point, liquid and vapor densities become equal, latent heat approaches zero, compressibility becomes highly sensitive to state, and ordinary property correlations may lose accuracy. Equipment operating near critical conditions requires validated property models and careful numerical treatment.

Property Tables and Diagrams as Engineering Tools

Property tables are often the most reliable and convenient source for water, steam, refrigerants, and selected industrial fluids. Common tables include saturation tables arranged by temperature or pressure, superheated-vapor tables, and compressed-liquid tables. A disciplined procedure is essential:

1. Identify the known independent properties and the likely phase region.

2. Compare the specified temperature and pressure with saturation data.

3. Select the appropriate table rather than assuming the substance is an ideal gas.

4. Interpolate only between nearby values and retain appropriate significant figures.

5. Check whether the resulting state is physically consistent.

At a specified pressure, if T is below T_sat, the state is generally a compressed liquid; if T equals T_sat, it may be saturated or two-phase; if T exceeds T_sat, it is superheated vapor. At a specified temperature, analogous comparisons can be made using saturation pressure.

Compressed-liquid tables are less commonly available because liquid properties vary relatively weakly with pressure. A frequently used approximation is

v(T,p) ≈ v_f(T),

and, for enthalpy,

h(T,p) ≈ h_f(T) + v_f(T)[p − p_sat(T)],

when pressure effects are important enough to include but detailed data are unavailable. This approximation should not be used blindly near the critical region or for strongly pressure-sensitive liquids.

Interpolation is not a substitute for a physical model. Linear interpolation is usually acceptable over small table intervals, but errors may become significant where properties change rapidly, especially near saturation and the critical point. A calculated quality outside the interval 0 ≤ x ≤ 1 signals either an incorrect phase assumption, inconsistent data, or an arithmetic error.

Equations of State for Real Fluids

An equation of state relates pressure, temperature, volume, and composition. The ideal-gas equation,

pv = RT,

is useful at low density and sufficiently high temperature, but it neglects molecular attraction and finite molecular size. Real-gas effects become important at high pressure, near saturation, and near the critical point.

The compressibility factor is defined as

Z = pv/(RT).

For an ideal gas, Z = 1. Values below one generally indicate that attractive forces reduce the pressure relative to ideal behavior at the same T and v. At sufficiently high density, repulsive effects dominate and Z may exceed one. The real-gas relation can therefore be written as

pv = ZRT.

For mixtures, the molar volume and the universal gas constant are used, and the molecular weight is required when converting between mass-specific and molar quantities.

Cubic Equations of State

Cubic equations of state provide a practical compromise between simplicity and physical realism. The van der Waals equation is

[p + a/v_m²](v_m − b) = RT,

where a represents attractive interactions and b represents excluded molecular volume. It explains qualitatively why real gases depart from ideality and predicts a liquid–vapor loop below the critical temperature. The mechanically unstable portion of that loop is replaced in practice by a constant-pressure phase-equilibrium line using the Maxwell equal-area construction.

More accurate cubic models include the Redlich–Kwong, Soave–Redlich–Kwong, and Peng–Robinson equations. Their attraction and repulsion parameters are fitted using critical properties and, often, the acentric factor. The Peng–Robinson equation is widely used for hydrocarbon vapor–liquid calculations because it offers useful predictions for both vapor and liquid densities over engineering ranges.

A cubic equation may yield three real volume roots in a two-phase-like region. The smallest root is commonly associated with a liquid-like state and the largest with a vapor-like state, but roots alone do not determine the equilibrium phase split. Phase stability and equal chemical potential, or equivalent fugacity criteria, must be considered. This distinction is important in process simulation: selecting a mathematically available root without checking phase stability can produce a physically incorrect result.

Corresponding States and Generalized Correlations

The principle of corresponding states states that fluids at the same reduced temperature and reduced pressure exhibit approximately similar behavior. Reduced properties are defined as

T_r = T/T_c, p_r = p/p_c.

Generalized compressibility charts provide Z as a function of reduced pressure and temperature. A third parameter, the acentric factor, improves the correlation for fluids whose molecular shape and polarity differ from the reference fluids used in the original two-parameter theory.

Corresponding-states methods are valuable for preliminary design, hand calculations, and checking software output. They are less reliable for strongly polar, associating, or chemically reacting fluids, and near the critical region. Their uncertainty should be considered when safety margins or equipment sizing depend strongly on density or enthalpy.

Vapor–Liquid Equilibrium and Chemical Potential

At equilibrium between liquid and vapor, temperature and pressure are uniform and the chemical potential of each component is equal in both phases. For a pure substance, this condition generates the saturation relationship p_sat(T). During vaporization at fixed pressure, energy is required to overcome intermolecular organization; the latent heat is

h_fg = h_g − h_f.

For mixtures, equilibrium is described by component fugacities. Fugacity is an effective pressure that accounts for nonideal molecular behavior. For component i,

f_i^L = f_i^V

at vapor–liquid equilibrium. In an ideal vapor and ideal liquid solution, this may reduce to Raoult's law,

y_i p = x_i p_i^{sat}(T),

where x_i and y_i are liquid- and vapor-phase mole fractions. Real mixtures require activity coefficients in the liquid phase and fugacity coefficients in the vapor phase. At high pressure, both phases may be significantly nonideal.

This framework explains why boiling temperature depends on composition and pressure, why distillation is possible, and why a pressure reduction can cause flashing. In a throttling valve, enthalpy is approximately conserved, but the pressure drop can move the fluid into the two-phase region. The resulting vapor fraction must be found from an isenthalpic flash calculation, not assumed from temperature alone.

Property Estimation in Engineering Practice

Property determination is best treated as a model-selection problem. The engineer should first establish the substance, composition, expected phase, pressure and temperature range, and required accuracy. A practical hierarchy is:

• Use validated tables for well-characterized pure substances in common operating ranges.
• Use an appropriate equation of state for high-pressure gases, hydrocarbons, and phase-equilibrium calculations.
• Use activity-coefficient models for strongly nonideal liquid mixtures at moderate pressure.
• Use generalized correlations for preliminary estimates or independent checks.
• Use specialized multiparameter formulations when near-critical accuracy, cryogenic performance, or reference-quality data are required.

Software packages such as REFPROP, CoolProp, and process simulators implement equations of state, transport-property correlations, and flash algorithms. Their use does not remove the need for engineering judgment. A software result should be accompanied by the selected fluid model, unit system, reference convention where relevant, phase specification, and an assessment of plausibility. A solver can converge to a mathematically valid but physically inappropriate phase root.

Useful validation checks include the ideal-gas limit at low density, monotonicity of density with pressure in a stable single phase, energy consistency in a cycle, and comparison with an independent data source. Units deserve particular attention: pressure may be entered as absolute or gauge pressure, and gas constants may be mass-specific or universal. Errors in these conversions can produce apparently reasonable but incorrect answers.

Applied Case: Flashing of a High-Pressure Water Stream

Background

Consider an illustrative process in which high-pressure liquid water is throttled through a valve to a lower pressure before entering a separator. The separator must be sized for the vapor generated by flashing.

Situation and Problem

Water enters the valve as a compressed liquid with negligible kinetic and potential-energy changes. Heat transfer is negligible, and the valve performs no shaft work. The inlet enthalpy is approximately 850 kJ/kg. Downstream pressure is 10 bar, for which representative saturation data are h_f = 763 kJ/kg and h_fg = 2015 kJ/kg. Determine the approximate vapor quality after throttling.

Analysis and Approach

For a steady, adiabatic throttling process,

h_1 ≈ h_2.

Because the downstream state lies in the two-phase region, write

h_2 = h_f + x h_fg.

Therefore,

x = (h_2 − h_f)/h_fg = (850 − 763)/2015 ≈ 0.043.

Approximately 4.3 percent of the outlet mass is vapor. Although this mass fraction is small, the vapor occupies a disproportionately large volume and can control separator dimensions and downstream pressure losses.

Outcome

The valve produces a two-phase mixture rather than a uniformly cooled liquid. A separator or suitable two-phase-flow design is therefore required. If a designer incorrectly treated the outlet as a compressed liquid, the predicted density, velocity, and pressure drop could be seriously wrong.

Learning Points

The case illustrates four principles: throttling is approximately isenthalpic; pressure reduction can create vapor without external heat input; quality is determined from downstream saturation properties; and vapor mass fraction does not indicate vapor volume fraction. In an actual design, the property data, valve coefficient, metastable effects, pressure losses, and separator performance would require further evaluation.

Limitations, Risks, and Common Mistakes

Several errors recur in property calculations. The most serious is identifying the phase region from temperature or pressure alone without comparing the state with saturation conditions. Another is using quality outside the two-phase region. Quality is not defined for a single-phase superheated vapor or compressed liquid.

Treating all gases as ideal can lead to large errors in compressor work, storage capacity, pipeline pressure drop, and relief-system sizing. Conversely, using an elaborate equation of state when only a rough estimate is needed can obscure the analysis and introduce unnecessary model uncertainty.

Reference-state inconsistency is a frequent source of apparent disagreement between data sources. Absolute enthalpy values from two sources may differ while enthalpy differences remain identical. Mixing values from unrelated references in one energy balance is unsafe.

Near phase boundaries, small measurement or rounding errors can cause a large change in calculated quality or density. Numerical solvers may also fail near the critical point because derivatives become ill-conditioned. Engineers should report the property model and uncertainty whenever the result affects equipment safety or regulatory compliance.

Professional and Regulatory Considerations

Property data influence pressure-vessel design, relief-valve sizing, refrigeration safety, emissions estimation, and process hazard analysis. Conservative design requires recognized data sources and documented methods. For refrigerants, toxic or flammable fluids, and high-pressure systems, the selected property formulation should be compatible with applicable standards and regulatory requirements.

Professional practice also requires traceability. A calculation should identify fluid composition, input conditions, units, software or table source, equation of state, phase assumption, and validation checks. When proprietary software is used, reproducible input files and version information should be retained. Engineers must not present uncertain estimates as precise measurements, particularly when decisions affect public safety or environmental performance.

Current Developments and Best Practices

Modern property packages increasingly combine multiparameter equations of state, molecular simulation, automated phase-stability analysis, and data-driven corrections. Digital twins and process simulators can update fluid properties as operating conditions change, while machine-learning models are being investigated for rapid surrogate calculations. Such methods are useful only when trained and validated within the relevant thermodynamic domain; extrapolation beyond the data range can be dangerous.

Best practice is to begin with a physical phase assessment, select the simplest model that meets the accuracy requirement, validate the result independently, and document assumptions. Near critical points, in highly nonideal mixtures, and in multiphase equipment, model selection should be reviewed by a specialist rather than delegated entirely to default software settings.

Chapter Summary

Thermodynamic properties define the state of a pure substance and provide the quantitative basis for energy and entropy balances. A simple compressible pure substance requires two independent intensive properties in a single-phase region, while the two-phase region requires an additional specification such as quality. Saturation tables, phase diagrams, and mixture relations allow properties to be determined efficiently when reliable data are available.

Real fluids depart from ideal-gas behavior because of molecular attraction and finite molecular volume. The compressibility factor, cubic equations of state, and corresponding-states correlations provide progressively more capable tools for estimating real-fluid behavior. Vapor–liquid equilibrium is governed by equality of chemical potential or fugacity, and flash calculations determine phase amounts after pressure or temperature changes.

Property software is powerful but must be used with physical judgment. Correct phase identification, consistent units and reference states, model validation, and awareness of uncertainty are essential. These practices connect textbook property determination with the safe and efficient design of thermal and fluid systems.

Key Terms

Intensive property: A property independent of the amount of substance.

Extensive property: A property proportional to the amount of substance.

Specific property: An extensive property expressed per unit mass.

State postulate: The principle that two independent intensive properties fix the state of a simple compressible pure substance.

Saturation state: An equilibrium state at which phase change is about to begin or end.

Quality: The mass fraction of vapor in a saturated liquid–vapor mixture.

Compressed liquid: A liquid at a pressure higher than its saturation pressure at the specified temperature, or at a temperature below saturation at the specified pressure.

Superheated vapor: Vapor at a temperature above saturation for its pressure.

Critical point: The termination of the liquid–vapor equilibrium curve.

Triple point: The state at which solid, liquid, and vapor coexist in equilibrium.

Latent heat: Energy required for phase change at equilibrium without a change in temperature for a pure substance at fixed pressure.

Compressibility factor: The ratio pv/RT, measuring departure from ideal-gas behavior.

Equation of state: A relation connecting pressure, temperature, volume, and composition.

Reduced property: A property divided by its critical value.

Acentric factor: A molecular-shape parameter used to improve corresponding-states correlations.

Fugacity: An effective pressure representing a species' escaping tendency in a nonideal phase.

Vapor–liquid equilibrium: Coexistence of liquid and vapor with equal chemical potential for each species.

Flash calculation: Determination of equilibrium phase amounts and compositions at specified conditions.

Reference state: The arbitrarily selected zero level for tabulated thermodynamic properties.

Review Questions

Short Answer Questions

1. Why are heat and work not classified as thermodynamic properties?

2. Under what conditions is quality defined, and what does it represent?

3. What is the physical significance of the compressibility factor Z?

4. Why are pressure and temperature not independent in the two-phase region of a pure substance?

5. State two reasons why reference-state consistency matters in energy calculations.

Long Answer Questions

1. Explain the phase behavior of a pure substance using a pressure–temperature diagram and a temperature–specific-volume diagram.

2. Derive the mixture-property relation Y = Y_f + xY_fg and discuss its application to enthalpy and specific volume.

3. Compare ideal-gas, corresponding-states, and cubic-equation-of-state methods for real-fluid property estimation.

4. Explain the thermodynamic basis of vapor–liquid equilibrium using chemical potential and fugacity.

5. Describe a rigorous procedure for selecting and validating a property model in process-simulation software.

Analytical / Application Questions

1. A saturated water mixture has h_f = 500 kJ/kg, h_g = 2700 kJ/kg, and h = 1600 kJ/kg. Determine its quality and explain why the result cannot be interpreted as a volume fraction.

2. A gas occupies 0.08 m³ at 12 MPa and 350 K. Explain how you would decide whether the ideal-gas equation is adequate and identify the data or correlations required for a better estimate.

3. A throttling valve reduces the pressure of a high-pressure liquid stream and produces a downstream two-phase mixture. Develop the calculation procedure for outlet quality and discuss the design consequences for a downstream separator.

22. Multiple Choice Questions with Answers

Question 1

A pure substance in the liquid–vapor region is specified by pressure and quality. Which statement is correct?

A. Temperature and pressure are independent variables in this region.

B. Quality is the vapor volume fraction.

C. Pressure and temperature are related by saturation conditions.

D. The substance must be superheated.

Correct Answer: C. Pressure and temperature are related by saturation conditions.

Question 2

A two-phase mixture has h_f = 600 kJ/kg, h_g = 2600 kJ/kg, and h = 1100 kJ/kg. Its quality is:

A. 0.20

B. 0.25

C. 0.40

D. 0.55

Correct Answer: B. 0.25

Question 3

The compressibility factor of a gas is less than one primarily when:

A. Repulsive forces dominate molecular behavior.

B. Attractive forces significantly reduce the pressure.

C. The gas has zero molecular volume.

D. The gas is necessarily at its critical point.

Correct Answer: B. Attractive forces significantly reduce the pressure.

Question 4

Which property method is generally most appropriate for a preliminary estimate of a nonpolar gas at moderate reduced pressure when limited data are available?

A. A generalized compressibility correlation

B. A saturated-liquid table only

C. A quality relation

D. A solid–liquid phase diagram

Correct Answer: A. A generalized compressibility correlation

Question 5

During an idealized adiabatic throttling process, which property is approximately conserved?

A. Temperature

B. Entropy in all cases

C. Enthalpy

D. Pressure

Correct Answer: C. Enthalpy

Question 6

A calculated quality of −0.12 for a supposed saturated mixture most directly indicates that:

A. The vapor mass fraction is negative physically.

B. The state assumption or property data are inconsistent.

C. The mixture is exactly at the critical point.

D. The fluid must be an ideal gas.

Correct Answer: B. The state assumption or property data are inconsistent.

Question 7

At vapor–liquid equilibrium for a component in a nonideal mixture, the fundamental equilibrium condition is equality of:

A. Density in the two phases

B. Specific volume in the two phases

C. Fugacity in the two phases

D. Quality in the two phases

Correct Answer: C. Fugacity in the two phases

Question 8

Why can a small vapor mass fraction create a large vapor volume fraction in a liquid–vapor mixture?

A. Vapor usually has a much larger specific volume than liquid.

B. Liquid has zero density at saturation.

C. Quality is defined as volume fraction.

D. Latent heat increases vapor density.

Correct Answer: A. Vapor usually has a much larger specific volume than liquid.

Question 9

Above the critical temperature, a pure substance:

A. Can always be liquefied by increasing pressure.

B. Cannot be liquefied by pressure alone.

C. Must exist as saturated liquid.

D. Has a compressibility factor of exactly zero.

Correct Answer: B. Cannot be liquefied by pressure alone.

Question 10

Which practice best improves the reliability of property-software results?

A. Accepting the default phase without checking it

B. Mixing enthalpy values from different reference conventions

C. Documenting the model, units, phase assumption, and validation checks

D. Using the most complex equation of state regardless of the fluid

Correct Answer: C. Documenting the model, units, phase assumption, and validation checks