Thermodynamic Property Relations

Introduction and Context

The first and second laws of thermodynamics provide energy balances and restrictions on energy conversion, but engineering analysis also requires relationships among measurable properties such as pressure, temperature, volume, entropy, and internal energy. Thermodynamic property relations supply these connections. They allow engineers to determine changes in entropy from pressure–volume data, estimate specific heats, calculate phase-equilibrium slopes, predict throttling behavior, and construct accurate real-fluid models.

The central mathematical idea is that thermodynamic properties are state functions. Their differentials are therefore exact, and the equality of mixed partial derivatives produces the Maxwell relations. These relations are especially valuable because they replace difficult-to-measure entropy derivatives with experimentally accessible pressure, volume, and temperature derivatives.

Learning Objectives

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

1. Define the thermodynamic potentials and derive the four Maxwell relations.
2. Apply Maxwell relations to obtain entropy, internal-energy, and enthalpy changes.
3. Derive and use general and special-case specific-heat relations.
4. Derive the Clapeyron and Clausius–Clapeyron equations and identify their limitations.
5. Explain the Joule–Thomson coefficient and distinguish throttling from free expansion.
6. Define residual properties and use them to correct ideal-gas estimates.
7. Interpret these tools in practical systems such as refrigeration and gas processing.

1. Fundamental Thermodynamic Potentials

For a simple compressible, closed system of fixed composition, the fundamental relation is

where U is internal energy (J), T is absolute temperature (K), S is entropy (J K⁻¹), p is pressure (Pa), and V is volume (m³). The relation assumes only p–V work and equilibrium states. For a fixed-mass system, the corresponding specific form is du = T ds − p dv.

Four useful thermodynamic potentials are:

For a multicomponent system, chemical-potential terms must be added. For example, dG = −S dT + V dp + ∑ μ_i dN_i. The single-component, fixed-composition equations are used throughout this chapter.

A differential such as dU = (∂U/∂S)_V dS + (∂U/∂V)_S dV is exact because U is a state function. Consequently,

This mathematical property is the basis of Maxwell relations.

2. Maxwell Relations

2.1 Derivation

From the fundamental equation, the coefficients of dS and dV are

(∂U/∂S)_V = T, (∂U/∂V)_S = −p.

Equality of mixed derivatives gives

Using the other potentials in the same manner gives the complete set:

The third and fourth forms are commonly used because thermal expansion and compressibility data are often available. Their validity requires a simple compressible equilibrium system with differentiable property functions; they are not independent equations of state.

2.2 Physical Interpretation

The relation

states that the entropy response to volume change at constant temperature equals the pressure response to temperature change at constant volume. Although entropy derivatives may be difficult to measure directly, pressure and volume derivatives can often be determined from an equation of state.

For a specific system, the entropy differential can be written as

or, alternatively,

Here c_v and c_p are specific heats at constant volume and pressure (J kg⁻¹ K⁻¹), and v is specific volume (m³ kg⁻¹). These equations are exact and apply to any single-phase simple compressible substance for which the derivatives exist.

3. Property Relations for Internal Energy and Enthalpy

Starting with du = T ds − p dv and substituting the first entropy differential gives

This is a general equation for a simple compressible substance. It shows that internal energy can depend on volume as well as temperature.

Using the second entropy differential and dh = T ds + v dp gives

Thus, enthalpy can depend on pressure as well as temperature.

3.1 Ideal-Gas Specialization

For an ideal gas, pv = RT, where R is the specific gas constant (J kg⁻¹ K⁻¹). Therefore,

(∂p/∂T)_v = R/v = p/T,

and the coefficient of dv in the internal-energy equation is zero:

Likewise, v = RT/p, so (∂v/∂T)_p = R/p = v/T, producing

Hence, for an ideal gas, u and h depend only on temperature. This result is not generally true for real gases.

For constant specific heats over a temperature interval,

For variable specific heats, the correct expressions are

3.2 Thermal Expansion and Isothermal Compressibility

Define the coefficient of volumetric expansion and isothermal compressibility as

where α has units K⁻¹ and κ_T has units Pa⁻¹. Substitution into the enthalpy relation gives

The general specific-heat difference is

Equivalent forms include

and

The signs are consistent because (∂v/∂p)_T < 0 for a mechanically stable substance. For an ideal gas, α = 1/T and κ_T = 1/p, so c_p − c_v = R. For liquids and solids, c_p and c_v are often close, but they are not exactly equal.

Worked Example 1: Specific-Heat Difference of an Ideal Gas

Air is approximated as an ideal gas with R = 0.287 kJ kg⁻¹ K⁻¹ and c_p = 1.005 kJ kg⁻¹ K⁻¹. Determine c_v and the ratio γ = c_p/c_v.

Using c_p − c_v = R,

Then

The units are consistent because both c_p and R are expressed in kJ kg⁻¹ K⁻¹.

4. Phase Equilibrium and the Clapeyron Equation

At equilibrium between two phases, the specific Gibbs free energies are equal:

For a pure substance, dg = −s dT + v dp. Along the coexistence curve, the differentials of the two phase Gibbs energies are equal:

Rearranging produces the Clapeyron equation:

Here h_{βα} is the latent enthalpy, or enthalpy change per unit mass, for the transition from phase α to phase β (J kg⁻¹). The equation applies to any first-order phase boundary in equilibrium, including vaporization, fusion, and sublimation.

For liquid–vapor equilibrium,

where p_sat is saturation pressure, h_fg is latent heat of vaporization, and v_g and v_f are saturated-vapor and saturated-liquid specific volumes.

The sign of the slope depends on both latent heat and the volume change. Vaporization normally has positive slope because h_fg > 0 and v_g > v_f. Ice–water equilibrium is unusual: v_ice is greater than v_liquid, so the fusion curve has a negative slope.

4.1 Clausius–Clapeyron Equation

For vaporization, if the vapor is approximated as an ideal gas and v_g ≫ v_f, then v_g ≈ RT/p_sat. The Clapeyron equation becomes

where R_v is the specific gas constant of the vapor. If h_fg is treated as constant over a limited temperature range, integration between states 1 and 2 gives

or equivalently,

Temperatures must be absolute, pressure units must be identical in the ratio, and h_fg must be expressed in J kg⁻¹ if R_v is in J kg⁻¹ K⁻¹. The constant-latent-heat approximation becomes less reliable over broad temperature ranges, especially near the critical point.

Worked Example 2: Estimating Saturation Pressure

Water vapor has R_v = 0.4615 kJ kg⁻¹ K⁻¹. Let p₁ = 10.0 kPa at T₁ = 318.15 K, and use h_fg = 2,400 kJ kg⁻¹ to estimate p₂ at T₂ = 328.15 K.

The reciprocal-temperature difference is approximately 9.57 × 10⁻⁵ K⁻¹, so

Thus,

This is an estimate because h_fg was treated as constant and the liquid-volume term was neglected.

5. Joule–Thomson Behavior

A steady throttling process through a valve, porous plug, or capillary is commonly modeled as adiabatic, with negligible changes in kinetic and potential energy and no shaft work. The steady-flow energy equation then gives

Although enthalpy remains constant, temperature may change because real-gas enthalpy depends on pressure. The Joule–Thomson coefficient is defined as

Its useful property form is

Using α,

The units are K Pa⁻¹ when pressure is measured in pascals. If μ_JT > 0, a pressure drop at constant enthalpy produces cooling. If μ_JT < 0, the gas warms during throttling. The inversion condition is μ_JT = 0, or αT = 1.

An ideal gas has αT = 1, so μ_JT = 0. Therefore, an ideal gas does not change temperature in an ideal throttling process. Real gases exhibit cooling or heating depending on temperature and pressure. At ordinary temperatures, nitrogen and carbon dioxide generally cool on throttling, whereas hydrogen and helium may warm unless first cooled below their inversion temperatures.

The Joule–Thomson effect must not be confused with free expansion. In free expansion into a vacuum, Q = 0, W = 0, and therefore ΔU = 0 for a closed system. In throttling, the appropriate steady-flow result is ΔH = 0. For an ideal gas both lead to no temperature change, but for real gases the property held constant differs.

6. Residual Properties

An equation of state or property model often gives real-fluid behavior, while ideal-gas data are more readily available. A residual property measures the departure from ideal-gas behavior at the same temperature and pressure:

where M may denote molar volume, enthalpy, entropy, or another property. The superscript R means residual and ig means ideal gas. Residual properties are zero in the ideal-gas limit as pressure approaches zero at fixed temperature.

Important examples are

For a pure fluid with molar volume v̄ and compressibility factor Z = pv̄/(RT), common pressure-integral forms are

Here R is the molar gas constant and the integrals follow an isotherm from zero pressure to the specified pressure. These equations assume a consistent equation of state and a single-phase path. Near phase boundaries, the correct phase-specific relation must be used.

Residual properties are useful in compressors, natural-gas pipelines, refrigeration systems, and chemical-process simulations. They separate temperature-dependent ideal-gas contributions from pressure-dependent real-fluid corrections, improving both calculation transparency and numerical implementation.

7. Integrated Property Evaluation

For a real gas, an enthalpy change can be evaluated as

Similarly,

For an ideal gas with variable specific heat,

This expression assumes an ideal gas and a fixed chemical composition. It is frequently used for compressor and turbine calculations when pressure ratios are known.

8. Applied Case: Throttling in a Natural-Gas Processing Plant

Background

An illustrative natural-gas plant receives a high-pressure methane-rich stream and reduces its pressure through a control valve before entering a separator. The valve is insulated, has no shaft, and is short enough that changes in elevation and velocity are negligible.

Problem

Operators observe that the outlet temperature is lower than the inlet temperature. They must determine whether this is consistent with thermodynamics and whether the temperature change can be predicted using an ideal-gas model.

Analysis

The valve is modeled as a steady throttling device. Therefore,

A pure ideal-gas model would give h = h(T), so equal enthalpy would imply T₂ = T₁. This model cannot explain cooling. A real-gas equation of state, however, predicts a pressure-dependent enthalpy. The observed temperature response is governed locally by

for an isenthalpic path. If μ_JT is positive and dp is negative, then dT is negative.

The engineer should calculate μ_JT at the inlet state using a validated equation of state, evaluate phase stability, and check whether the outlet enters a two-phase region. If flashing occurs, a single-phase Joule–Thomson coefficient is insufficient; an isenthalpic flash calculation is required.

Approach and Outcome

A practical workflow is:

1. Specify inlet pressure, temperature, composition, and flow rate.
2. Select an equation of state suitable for hydrocarbon mixtures.
3. Calculate inlet enthalpy and entropy.
4. Set outlet pressure and solve the valve condition h₂ = h₁.
5. Check phase equilibrium and determine vapor and liquid fractions if necessary.
6. Confirm that predicted temperature and phase conditions are safe for downstream equipment.

The ideal-gas estimate provides a useful screening result but predicts no throttling temperature change. The real-fluid model explains the cooling and identifies possible liquid formation. This distinction is important because liquid carryover can damage compressors and upset flow meters.

Learning Points

The case demonstrates that enthalpy conservation does not imply temperature conservation. It also shows why residual enthalpy, phase-equilibrium calculations, and equipment-specific assumptions must be combined rather than applied independently.

9. Limitations, Best Practices, and Professional Relevance

The relations in this chapter are exact within their thermodynamic assumptions, but engineering results depend on property data and model selection. Important limitations include:

• Maxwell relations do not replace an equation of state or experimental property data.
• Clausius–Clapeyron becomes inaccurate when latent heat varies strongly or the vapor is nonideal.
• The Joule–Thomson coefficient is a local derivative and should not be extrapolated over large pressure changes without integration.
• Residual-property calculations require a thermodynamically consistent equation of state.
• Near the critical point, derivatives can become large and numerical calculations require careful phase and stability checks.
• Mixtures require composition variables, chemical potentials, and appropriate mixing rules.

Best practice is to use absolute temperature, consistent energy and pressure units, clearly stated phase assumptions, and independent checks against trusted property tables or software. In safety-critical applications, property calculations should be traceable, version-controlled, and validated over the expected operating range. Software does not remove the need to inspect phase regions, units, convergence, and physical plausibility.

Modern process simulators use equations of state, automatic differentiation, and flash algorithms to evaluate these relations rapidly. Artificial intelligence may assist with property-model fitting or anomaly detection, but it should not replace conservation-law checks, uncertainty assessment, or engineering responsibility for the selected thermodynamic model.

Chapter Summary

Thermodynamic property relations connect state variables through the mathematical structure of thermodynamic potentials. The four Maxwell relations follow from equality of mixed partial derivatives and convert entropy derivatives into measurable pressure–volume–temperature derivatives. They lead to general expressions for entropy, internal-energy, and enthalpy changes.

The difference between specific heats is c_p − c_v = Tvα²/κ_T, reducing to R for an ideal gas. The Clapeyron equation gives the slope of any equilibrium phase boundary, while the Clausius–Clapeyron equation provides a useful vapor-pressure approximation when the vapor is ideal and the liquid volume is negligible. The Joule–Thomson coefficient describes temperature change during isenthalpic throttling; it is zero for an ideal gas but generally nonzero for real gases. Residual properties quantify real-fluid departures from ideal-gas behavior and are essential in high-pressure equipment and process simulation.

Key Terms

Maxwell relation: An equality between thermodynamic partial derivatives obtained from a thermodynamic potential.

Thermodynamic potential: A state function such as U, H, A, or G used to describe equilibrium systems.

Exact differential: A differential whose integral depends only on end states.

Specific heat: Energy required per unit mass per unit temperature rise under a specified constraint.

Thermal expansion coefficient: α = (1/v)(∂v/∂T)_p.

Isothermal compressibility: κ_T = −(1/v)(∂v/∂p)_T.

Ideal gas: A model satisfying pv = RT with negligible intermolecular effects.

Clapeyron equation: Relation giving the slope of a phase-equilibrium curve.

Clausius–Clapeyron equation: An approximate integrated vapor-pressure relation.

Latent heat: Enthalpy change associated with a phase change at equilibrium.

Saturation pressure: Pressure at which two phases coexist at a specified temperature.

Joule–Thomson coefficient: μ_JT = (∂T/∂p)_h.

Throttling process: A steady-flow, approximately adiabatic process with constant enthalpy.

Inversion condition: State at which the Joule–Thomson coefficient is zero.

Residual property: Difference between a real-fluid property and its ideal-gas value at the same T and p.

Compressibility factor: Z = pv̄/(RT).

Equation of state: Mathematical relation among thermodynamic state variables.

Phase equilibrium: Condition in which phases coexist with equal temperature, pressure, and chemical potential.

Critical point: State at which liquid and vapor become indistinguishable.

Review Questions

Short Answer Questions

1. Why do Maxwell relations exist for thermodynamic state functions?
2. State the four Maxwell relations.
3. Under what assumptions does c_p − c_v = R apply?
4. What physical process is described by the Joule–Thomson coefficient?
5. Define a residual property.

Long Answer Questions

1. Derive the Maxwell relation (∂S/∂V)_T = (∂p/∂T)_V from the Helmholtz free energy.
2. Derive the general differential equation for internal energy in terms of dT and dv.
3. Derive the Clapeyron equation from equality of Gibbs free energy in two phases.
4. Explain the assumptions used to obtain the Clausius–Clapeyron equation from the Clapeyron equation.
5. Compare free expansion and throttling, including the property conserved in each process.

Analytical / Application Questions

1. A gas has α = 1.8 × 10⁻³ K⁻¹, v = 0.020 m³ kg⁻¹, T = 300 K, and c_p = 1.1 kJ kg⁻¹ K⁻¹. Determine μ_JT and state whether a pressure drop causes cooling or heating.
2. Use the Clapeyron equation to explain why the solid–liquid equilibrium curve of water has a negative slope near the melting point.
3. A high-pressure hydrocarbon valve produces a lower outlet temperature than predicted by an ideal-gas model. Identify the missing physical effect and propose a calculation procedure.

Multiple Choice Questions with Answers

Question 1

Which thermodynamic potential directly produces the Maxwell relation (∂S/∂V)_T = (∂p/∂T)_V?

A. Internal energy U

B. Enthalpy H

C. Helmholtz free energy A

D. Gibbs free energy G

Correct Answer: C. Helmholtz free energy A

Question 2

For an ideal gas, which statement is correct?

A. u depends only on pressure

B. h depends only on temperature

C. c_p equals c_v

D. The Joule–Thomson coefficient is always positive

Correct Answer: B. h depends only on temperature

Question 3

The general relation for the specific-heat difference is:

A. c_p − c_v = R for every substance

B. c_p − c_v = Tvα²/κ_T

C. c_p − c_v = κ_T/(Tvα²)

D. c_p − c_v = 0 for every liquid

Correct Answer: B. c_p − c_v = Tvα²/κ_T

Question 4

The Clapeyron equation applies to:

A. Only ideal-gas expansion

B. Any equilibrium first-order phase boundary

C. Only irreversible throttling

D. Only incompressible liquids

Correct Answer: B. Any equilibrium first-order phase boundary

Question 5

In the Clausius–Clapeyron approximation, which assumption is commonly made?

A. The vapor specific volume is negligible

B. The vapor behaves as an ideal gas

C. Latent heat is exactly zero

D. Temperature is measured in Celsius

Correct Answer: B. The vapor behaves as an ideal gas

Question 6

For an ideal gas undergoing an ideal throttling process, the Joule–Thomson temperature change is:

A. Always positive

B. Always negative

C. Zero

D. Undefined because enthalpy changes

Correct Answer: C. Zero

Question 7

If μ_JT is positive, a pressure drop at constant enthalpy generally causes the gas to:

A. Cool

B. Heat

C. Remain at saturation pressure

D. Condense necessarily

Correct Answer: A. Cool

Question 8

A residual enthalpy is defined as:

A. h^ig − h at different temperatures

B. h − h^ig at the same temperature and pressure

C. h divided by pressure

D. The latent heat of vaporization

Correct Answer: B. h − h^ig at the same temperature and pressure

Question 9

For a steady, adiabatic throttling valve with negligible kinetic and potential-energy changes, which quantity is constant?

A. Internal energy

B. Entropy

C. Enthalpy

D. Temperature

Correct Answer: C. Enthalpy

Question 10

For water near its melting point, the solid–liquid coexistence curve has a negative slope mainly because:

A. The latent heat of fusion is negative

B. Liquid water has a larger specific volume than ice

C. Ice has a larger specific volume than liquid water

D. Water vapor is nonideal

Correct Answer: C. Ice has a larger specific volume than liquid water