Joule-Thomson Expansion

Turn the valve on a canister of compressed CO₂ and the nozzle frosts over instantly. Turn the valve on a hydrogen cylinder under similar conditions. The nozzle gets warmer instead. Both processes are the same kind of event. A gas is forced through a throttle from high pressure to low pressure with no heat exchange. Why do they cool or heat depending on the gas and the starting temperature? This throttling process is called Joule-Thomson expansion. It is the working principle behind gas liquefaction, refrigerators, and the frost you see on a compressed-air nozzle. Unlike the free (Joule) expansion of an ideal gas, throttling conserves enthalpy rather than internal energy. It also exposes real intermolecular forces that an ideal gas hides completely. This article derives the Joule-Thomson coefficient μJT\mu_{JT} from the First Law. It connects this coefficient to measurable quantities like CpC_p and the equation of state. The article also works two full numerical examples. It sets up the idea of an inversion temperature. The next article develops this into the practical theory of gas liquefaction.

Key takeaways
  • Joule-Thomson expansion is adiabatic throttling through a porous plug or valve at constant enthalpy, not constant internal energy.
  • The Joule-Thomson coefficient is defined as μJT=(∂T/∂p)H\mu_{JT} = (\partial T/\partial p)_H, measured directly by experiment.
  • For an ideal gas μJT=0\mu_{JT} = 0 exactly; any nonzero value is a real-gas effect from intermolecular forces.
  • μJT\mu_{JT} can be positive (cooling on expansion) or negative (heating on expansion), depending on temperature and the gas.
  • The sign flips at the inversion temperature, the boundary that governs whether throttling can liquefy a gas.
  • Key idea: Throttling conserves enthalpy, so temperature changes with pressure only because real molecules interact — measured by μJT=(∂T/∂p)H\mu_{JT} = (\partial T/\partial p)_H.
The core insight

Throttling a gas through a valve keeps enthalpy constant because the flow work done pushing gas in exactly equals the flow work done pushing gas out; any temperature change on the low-pressure side is then a direct, quantitative fingerprint of intermolecular attractions and repulsions, vanishing identically for an ideal gas and reversing sign at a characteristic inversion temperature.

The physical picture

Picture a rigid, thermally insulated tube divided by a porous plug or a narrow throttle valve. Gas at pressure p1p_1 and temperature T1T_1 sits upstream; gas emerges downstream at a lower pressure p2p_2 and some new temperature T2T_2. No heat crosses the insulated walls, and no shaft work is extracted — the only work done is the flow work of pushing gas through the plug. This is fundamentally different from the free Joule expansion into a vacuum, where the prerequisite result is that an ideal gas undergoes no temperature change because no external work is done and no attractive forces exist to overcome. In throttling, the gas does do net work on itself, and that work is done against intermolecular forces if any exist. Consider a fixed amount of gas, one mole, pushed steadily through the plug. Upstream, a piston-like action does work p1V1p_1V_1 pushing the gas in; downstream, the gas pushes with work p2V2p_2V_2 expanding into the low-pressure region. Since the tube is insulated and there is no other heat or work exchange, the First Law applied to this steady-flow segment gives conservation of the sum U+pVU + pV, which is enthalpy. The process is isenthalpic: H1=H2H_1 = H_2, even though pressure and volume both change. Because HH is conserved rather than TT, any temperature change reflects how enthalpy is partitioned between UU and pVpV as pressure drops. For an ideal gas, HH depends only on TT, so constant HH forces constant TT: no temperature change at all. For a real gas, HH depends on both TT and pp because intermolecular potential energy contributes to UU. Dropping the pressure changes the average intermolecular separation, which changes the potential energy contribution to UU, and the temperature must shift to keep HH fixed. Physically: if attractive forces dominate (gas well below its Boyle-like regime), pulling molecules apart during expansion costs energy, which is drawn from kinetic energy — the gas cools. If repulsive forces dominate (gas already very compressed, or at high temperature where kinetic energies swamp attractive wells), expansion releases stored repulsive potential energy into kinetic energy — the gas warms.

Figure
Throttling apparatus: adiabatic, constant-enthalpy expansion across a porous plug
upstream: p1, T1, V1 downstream: p2, T2, V2 porous plug H1 = H2 (isenthalpic), p1 > p2, adiabatic, no shaft work
Figure
Joule-Thomson coefficient versus temperature at fixed pressure, showing the inversion temperature
Temperature T (K) mu JT (K per atm) T_inv ≈ 621 K (N2) mu JT > 0: cooling mu JT < 0: heating

The core relations

H1=H2(throttling is isenthalpic)H_1 = H_2 \quad \text{(throttling is isenthalpic)}

where: H1,H2H_1, H_2 are the molar (or total) enthalpies upstream and downstream of the plug.

μJT≡(∂T∂p)H\mu_{JT} \equiv \left(\frac{\partial T}{\partial p}\right)_H

where: μJT\mu_{JT} is the Joule-Thomson coefficient, the slope of an isenthalp on a TT–pp diagram, typically in K atm⁻¹ or K bar⁻¹.

μJT=−1Cp(∂H∂p)T\mu_{JT} = -\frac{1}{C_p}\left(\frac{\partial H}{\partial p}\right)_T

where: CpC_p is the constant-pressure heat capacity, always positive, so the sign of μJT\mu_{JT} is opposite the sign of (∂H/∂p)T(\partial H/\partial p)_T.

μJT=1Cp[T(∂V∂T)p−V]\mu_{JT} = \frac{1}{C_p}\left[T\left(\frac{\partial V}{\partial T}\right)_p - V\right]

where: VV is molar volume; this form uses only CpC_p and the equation of state, both experimentally accessible.

μJT≈1Cp(2aRT−b)(van der Waals, low-density limit)\mu_{JT} \approx \frac{1}{C_p}\left(\frac{2a}{RT} - b\right) \quad \text{(van der Waals, low-density limit)}

where: a,ba, b are van der Waals constants; this predicts μJT=0\mu_{JT} = 0 at the inversion temperature Tinv=2a/(Rb)T_{inv} = 2a/(Rb).

The core relations

Derivation AFrom the First Law to μ_JT and the van der Waals estimate

Goal. Show that throttling conserves enthalpy, then derive μJT\mu_{JT} in terms of measurable quantities, and finally evaluate it for a van der Waals gas.

Step 1 — Throttling is isenthalpic. Consider one mole of gas pushed steadily through the plug. Take the gas between two moving boundaries: one just upstream, one just downstream, each enclosing a fixed amount of substance as it crosses. Upstream, an imaginary piston pushes the mole of gas from volume V1V_1 down to zero, doing work w1=p1V1w_1 = p_1V_1 on the system (work done on the gas by whatever is behind it). Downstream, the gas expands from zero to V2V_2 against the plug's other side, doing work w2=−p2V2w_2 = -p_2V_2 (work done by the gas). No heat crosses the insulated walls, so q=0q = 0. By the First Law,

ΔU=U2−U1=w1+w2=p1V1−p2V2.\Delta U = U_2 - U_1 = w_1 + w_2 = p_1V_1 - p_2V_2 .

Rearranging,

U2+p2V2=U1+p1V1⟹H2=H1.U_2 + p_2V_2 = U_1 + p_1V_1 \quad \Longrightarrow \quad H_2 = H_1 .

Enthalpy is conserved even though neither UU nor pVpV individually is conserved.

Step 2 — Express μJT\mu_{JT} using CpC_p. Treat H=H(T,p)H = H(T,p) as a state function and write its total differential,

dH=(∂H∂T)pdT+(∂H∂p)Tdp=Cp dT+(∂H∂p)Tdp.dH = \left(\frac{\partial H}{\partial T}\right)_p dT + \left(\frac{\partial H}{\partial p}\right)_T dp = C_p\, dT + \left(\frac{\partial H}{\partial p}\right)_T dp .

Along an isenthalp, dH=0dH = 0, so

0=Cp dT+(∂H∂p)Tdp⟹(∂T∂p)H=−1Cp(∂H∂p)T.0 = C_p\, dT + \left(\frac{\partial H}{\partial p}\right)_T dp \quad \Longrightarrow \quad \left(\frac{\partial T}{\partial p}\right)_H = -\frac{1}{C_p}\left(\frac{\partial H}{\partial p}\right)_T .

This is the general definition of μJT\mu_{JT} in terms of the isothermal pressure dependence of enthalpy — a quantity that vanishes identically for an ideal gas because HidealH_{ideal} depends on TT alone.

Step 3 — Convert (∂H/∂p)T(\partial H/\partial p)_T to equation-of-state variables. Starting from H=U+pVH = U + pV and using a Maxwell relation derived from dG=−SdT+VdpdG = -SdT + Vdp (specifically (∂S/∂p)T=−(∂V/∂T)p(\partial S/\partial p)_T = -(\partial V/\partial T)_p), one obtains the standard thermodynamic identity

(∂H∂p)T=V−T(∂V∂T)p.\left(\frac{\partial H}{\partial p}\right)_T = V - T\left(\frac{\partial V}{\partial T}\right)_p .

Substituting into Step 2,

μJT=−1Cp[V−T(∂V∂T)p]=1Cp[T(∂V∂T)p−V].\mu_{JT} = -\frac{1}{C_p}\left[V - T\left(\frac{\partial V}{\partial T}\right)_p\right] = \frac{1}{C_p}\left[T\left(\frac{\partial V}{\partial T}\right)_p - V\right].

For an ideal gas, V=RT/pV = RT/p, so T(∂V/∂T)p=RT/p=VT(\partial V/\partial T)_p = RT/p = V, and the bracket vanishes identically: μJT=0\mu_{JT} = 0 for all T,pT, p. This confirms that any measured μJT≠0\mu_{JT} \neq 0 is entirely a real-gas effect.

Step 4 — Evaluate for a van der Waals gas. The van der Waals equation of state is

(p+aVm2)(Vm−b)=RT.\left(p + \frac{a}{V_m^2}\right)(V_m - b) = RT .

To leading order in the small quantities a/(VmRT)a/(V_mRT) and b/Vmb/V_m (moderate pressure, not too close to condensation), one can solve for VmV_m and differentiate to show

T(∂Vm∂T)p−Vm≈2aRT−b.T\left(\frac{\partial V_m}{\partial T}\right)_p - V_m \approx \frac{2a}{RT} - b .

Substituting into Step 3's result,

μJT≈1Cp(2aRT−b)\boxed{\mu_{JT} \approx \frac{1}{C_p}\left(\frac{2a}{RT} - b\right)}

Setting this to zero gives the (Boyle-type) inversion temperature

Tinv=2aRb.T_{inv} = \frac{2a}{Rb} .

Below TinvT_{inv}, attractive forces (aa) dominate and μJT>0\mu_{JT} > 0: the gas cools on expansion. Above TinvT_{inv}, the size term (bb, repulsive core) dominates and μJT<0\mu_{JT} < 0: the gas warms on expansion.

Worked examples

Example 1 — Foundational: computing μ_JT from tabulated data for nitrogen. At 300 K and 1 atm, nitrogen gas has Cp,m=29.1 J K−1mol−1C_{p,m} = 29.1\ \text{J K}^{-1}\text{mol}^{-1}. Experimental data show that near this state, throttling nitrogen from 20 atm to 1 atm at an average temperature of 300 K produces a temperature drop of ΔT=−21.5\Delta T = -21.5 K per 19 atm drop (a standard tabulated result consistent with μJT(N2,300 K)≈0.25\mu_{JT}(N_2, 300\text{ K}) \approx 0.25 K atm⁻¹). Compute the expected temperature change for a smaller, controlled throttle from p1=5.00p_1 = 5.00 atm to p2=4.00p_2 = 4.00 atm at 300 K, assuming μJT\mu_{JT} is roughly constant over this narrow interval. Step 1. Use the defining relation in finite-difference form, valid since Δp\Delta p is small compared to the pressure range over which μJT\mu_{JT} varies appreciably:

ΔT≈μJT Δp.\Delta T \approx \mu_{JT}\,\Delta p .

Step 2. Insert numbers, with Δp=p2−p1=4.00−5.00=−1.00\Delta p = p_2 - p_1 = 4.00 - 5.00 = -1.00 atm:

ΔT≈(0.25 K atm−1)(−1.00 atm)=−0.25 K.\Delta T \approx (0.25\ \text{K atm}^{-1})(-1.00\ \text{atm}) = -0.25\ \text{K}.

Step 3. Sanity check. The pressure drop is small (1 atm out of 5), so linear extrapolation of μJT\mu_{JT} is reasonable, and the sign is correct: nitrogen at room temperature is well below its inversion temperature (621 K, from Table data with a=1.408a = 1.408 atm L² mol⁻², b=0.03913b = 0.03913 L mol⁻¹, giving Tinv=2a/(Rb)=2(1.408)/(0.08206×0.03913)≈876T_{inv} = 2a/(Rb) = 2(1.408)/(0.08206 \times 0.03913) \approx 876 K by the crude van der Waals estimate — real N2N_2's measured inversion temperature is closer to 621 K because the low-density van der Waals approximation overestimates TinvT_{inv}). Either way, T≪TinvT \ll T_{inv}, so cooling on expansion, μJT>0\mu_{JT} > 0, is expected. The predicted −0.25-0.25 K cooling for a 1 atm drop is small and physically sensible — throttling refrigerators use multi-atmosphere pressure drops precisely because the per-atmosphere effect is modest.

Example 2 — Applied: van der Waals estimate and liquefaction feasibility for CO₂ versus H₂. Use the low-density van der Waals formula to estimate μJT\mu_{JT} and TinvT_{inv} for carbon dioxide and hydrogen, and decide whether room-temperature throttling can cool each gas. Data: for CO₂, a=3.640a = 3.640 atm L² mol⁻², b=0.04267b = 0.04267 L mol⁻¹, Cp,m≈37.1C_{p,m} \approx 37.1 J K⁻¹mol⁻¹ near 300 K. For H₂, a=0.2476a = 0.2476 atm L² mol⁻², b=0.02661b = 0.02661 L mol⁻¹, Cp,m≈28.8C_{p,m} \approx 28.8 J K⁻¹mol⁻¹ near 300 K. Use R=0.08206R = 0.08206 L atm K⁻¹ mol⁻¹ for the van der Waals terms and R=8.314R = 8.314 J K⁻¹mol⁻¹ for energy units. Step 1. Inversion temperatures:

Tinv(CO2)=2aRb=2(3.640)(0.08206)(0.04267)=7.2800.003502≈2079 K.T_{inv}(\text{CO}_2) = \frac{2a}{Rb} = \frac{2(3.640)}{(0.08206)(0.04267)} = \frac{7.280}{0.003502} \approx 2079\ \text{K}.
Tinv(H2)=2aRb=2(0.2476)(0.08206)(0.02661)=0.49520.002184≈227 K.T_{inv}(\text{H}_2) = \frac{2a}{Rb} = \frac{2(0.2476)}{(0.08206)(0.02661)} = \frac{0.4952}{0.002184} \approx 227\ \text{K}.

Step 2. Interpretation at T=300T = 300 K. For CO₂, 300 K≪2079 K300\ \text{K} \ll 2079\ \text{K}, so CO₂ sits deep in its cooling regime; μJT>0\mu_{JT} > 0 strongly. For H₂, 300 K>227 K300\ \text{K} > 227\ \text{K}, so room-temperature hydrogen sits above its inversion temperature; μJT<0\mu_{JT} < 0 and hydrogen warms on throttling at 300 K. Step 3. Quantify μJT\mu_{JT} at 300 K for each gas. Convert aa from atm L² mol⁻² to SI (J m³ mol⁻²) is unnecessary if we keep consistent units: work in L·atm for the bracket, then convert to J using 1 L atm=101.3251\ \text{L atm} = 101.325 J. For CO₂:

2aRT−b=2(3.640)(0.08206)(300)−0.04267=7.28024.62−0.04267=0.2957−0.04267=0.2530 L mol−1.\frac{2a}{RT} - b = \frac{2(3.640)}{(0.08206)(300)} - 0.04267 = \frac{7.280}{24.62} - 0.04267 = 0.2957 - 0.04267 = 0.2530\ \text{L mol}^{-1}.

Convert to J mol⁻¹ atm⁻¹... actually the bracket has units of volume (L mol⁻¹); dividing by CpC_p (J K⁻¹mol⁻¹) requires converting the bracket to J atm⁻¹ mol⁻¹ using 1 L atm=101.3251\ \text{L atm} = 101.325 J:

0.2530 L mol−1×101.325 J L−1atm−1=25.64 J atm−1mol−1.0.2530\ \text{L mol}^{-1} \times 101.325\ \text{J L}^{-1}\text{atm}^{-1} = 25.64\ \text{J atm}^{-1}\text{mol}^{-1}.
μJT(CO2)=25.64 J atm−1mol−137.1 J K−1mol−1≈0.69 K atm−1.\mu_{JT}(\text{CO}_2) = \frac{25.64\ \text{J atm}^{-1}\text{mol}^{-1}}{37.1\ \text{J K}^{-1}\text{mol}^{-1}} \approx 0.69\ \text{K atm}^{-1}.

This is close to the accepted experimental value near 1.1 K atm⁻¹ at 300 K (the low-density van der Waals estimate is approximate, but the sign and order of magnitude are correct). For H₂:

2aRT−b=2(0.2476)(0.08206)(300)−0.02661=0.495224.62−0.02661=0.02012−0.02661=−0.00649 L mol−1.\frac{2a}{RT} - b = \frac{2(0.2476)}{(0.08206)(300)} - 0.02661 = \frac{0.4952}{24.62} - 0.02661 = 0.02012 - 0.02661 = -0.00649\ \text{L mol}^{-1}.

Convert: −0.00649×101.325=−0.658 J atm−1mol−1-0.00649 \times 101.325 = -0.658\ \text{J atm}^{-1}\text{mol}^{-1}.

μJT(H2)=−0.65828.8≈−0.023 K atm−1.\mu_{JT}(\text{H}_2) = \frac{-0.658}{28.8} \approx -0.023\ \text{K atm}^{-1}.

Step 4. Sanity check and conclusion. Both signs match qualitative prediction: CO₂ cools substantially (about 0.7-1 K per atmosphere dropped, consistent with why CO₂ fire extinguishers frost heavily on discharge), while hydrogen very slightly warms per atmosphere at 300 K, consistent with the historical fact that hydrogen and helium must first be pre-cooled below their (low) inversion temperatures before Joule-Thomson throttling can liquefy them — the technological puzzle the next article resolves in detail.

Limiting cases & common traps

SituationBehavior
Ideal gas, any T,pT, pμJT=0\mu_{JT} = 0 exactly; no temperature change on throttling
T<TinvT < T_{inv}μJT>0\mu_{JT} > 0; gas cools on expansion (attractive forces dominate)
T>TinvT > T_{inv}μJT<0\mu_{JT} < 0; gas warms on expansion (size/repulsion dominates)
Free (Joule) expansionConstant UU, not constant HH; ideal gas shows no ΔT\Delta T either way, but the conserved quantity differs
Very high pressure, dense gasSimple van der Waals low-density formula for μJT\mu_{JT} breaks down; use full equation of state

The chemist behind it

James Prescott Joule had already established, with his free-expansion experiment, that an ideal gas shows no cooling when expanded without doing external work — a crucial but frustratingly inconclusive test for real gases, since the heat capacity of the surrounding water bath was far larger than that of the gas itself, masking any small real-gas effect. Joule then collaborated with William Thomson (Lord Kelvin), who had been developing an absolute temperature scale and needed a more sensitive way to detect deviations from ideal-gas behavior. Together they redesigned the experiment: instead of free expansion into vacuum, they forced gas continuously through a porous plug under steady pressures on each side, measuring the small but reproducible temperature change directly in the flowing gas rather than in a bath. This porous-plug experiment, reported in the 1850s, gave the first clean, quantitative evidence of intermolecular forces in real gases and directly enabled the practical liquefaction technology that followed within decades.

Origins · 1 minWilliam Thomson, 1st Baron Kelvin and how this was discovered
WT
William Thomson, 1st Baron Kelvin
1824-1907 · Co-developed the porous-plug throttling experiment with Joule and used it to refine the absolute temperature scale and detect real-gas deviations

Ask the tutor: How does the inversion curve of a real gas look on a full T-p diagram, and why does it eventually close at high pressure?

Why is enthalpy conserved in throttling but not internal energy or volume?

Throttling involves flow work: gas upstream does work p1V1p_1V_1 pushing the sample through, and the sample does work p2V2p_2V_2 expanding out the other side. These two work terms do not cancel individually, so UU changes by exactly p1V1−p2V2p_1V_1 - p_2V_2. Rearranged, that statement is precisely U2+p2V2=U1+p1V1U_2 + p_2V_2 = U_1 + p_1V_1, which defines conservation of H=U+pVH = U + pV. Volume is not conserved either, since the gas is expanding. Only the specific combination U+pVU + pV stays fixed, which is why enthalpy — not energy — is the natural variable for steady-flow processes at constant pressure segments.

How does μ_JT relate to Cp and Cv from the earlier relation between them?

The Cp-Cv relation shows that Cp−Cv=RC_p - C_v = R for an ideal gas and involves (∂p/∂T)V(\partial p/\partial T)_V and (∂V/∂T)p(\partial V/\partial T)_p for real gases through the same equation-of-state derivatives that appear in μJT\mu_{JT}. Both quantities trace back to how VV depends on TT at fixed pp, so a gas with strongly nonideal (∂V/∂T)p(\partial V/\partial T)_p typically shows both a large Cp−CvC_p - C_v difference and a large μJT\mu_{JT}. They are not numerically equal, but they share the same underlying thermodynamic derivative and both vanish in the ideal-gas limit for the same reason: V=RT/pV = RT/p is perfectly linear in TT.

Can μ_JT ever change sign twice for the same gas?

Yes. A full equation of state (not just the low-density van der Waals approximation) generally gives an upper and a lower inversion temperature at each pressure, and the inversion curve itself is a closed loop in the TT-pp plane. Below the lower inversion temperature and above the upper inversion temperature, μJT\mu_{JT} can be negative, while it is positive in between. At very high pressures the entire inversion curve can disappear, meaning throttling produces heating at every temperature. This richer structure is explored fully when discussing inversion curves and liquefaction cycles.

Is the Joule-Thomson coefficient the same thing as the cooling used in a household refrigerator?

Not exactly. Vapor-compression refrigerators (the common household type) rely mainly on latent heat absorbed during evaporation of a liquid refrigerant, with throttling used only to drop the pressure of the liquid before it evaporates. The temperature drop from that throttling step is a genuine Joule-Thomson effect, but the dominant cooling in the cycle comes from the enthalpy of vaporization, not from μJT\mu_{JT} alone. Pure gas liquefaction plants for air, nitrogen, or helium, by contrast, rely on repeated Joule-Thomson throttling of the gas itself as the primary cooling mechanism.

Before this

Where this leads