Chemical & Process Engineering / Gas Solubility Henry's Law
FREE TO EXPLORE
LABORATORY 39 · CHEMICAL & PROCESS ENGINEERING

Gas Solubility Henry's Law theory

Understand gas solubility henry's law: governing transport equations, heat and mass balances, reaction kinetics, and unit operations.

Chemical engineering transport & equilibrium foundations of Gas Solubility Henry's Law.

Gas solubility in liquids is a core physical phenomenon across chemical engineering, environmental science, beverage carbonation, bioprocessing, and carbon capture and sequestration (CCS). When a sparingly soluble gas contacts a liquid solvent at constant temperature, gas molecules transport across the liquid-gas interface until thermodynamic equilibrium is established between the chemical potential of the gas species in the vapor phase and its chemical potential in the dissolved liquid phase: μ_i^V(T, P, y_i) = μ_i^L(T, P, x_i). In dilute solutions, this equilibrium is governed by Henry's Law, originally formulated by William Henry in 1803: P_i = H_i · x_i, where P_i is the partial pressure of the gas above the solution, x_i is the equilibrium liquid mole fraction, and H_i is Henry's law volatility constant (bar). Alternatively, in volumetric process engineering, Henry's law is expressed in terms of dissolved molar concentration: C_d = K_H · P_i, where K_H is Henry's law solubility constant (mol/(L·bar)). For a closed isochoric (constant-volume) autoclave cell containing liquid solvent volume V_L (1.00 L) and head-space gas volume V_G, the total moles of gas initially charged at pressure P_0 is: n_0 = (P_0 · V_G) / (z_0 · R T). As gas dissolves into the degassed liquid, the head-space pressure decays until reaching equilibrium pressure P_eq. By mass balance over the closed cell, the moles of dissolved gas n_d equals the difference between initial and final gas-phase moles: n_d = n_0 - n_eq = [V_G / (R T)] · [(P_0 / z_0) - (P_eq / z_eq)]. In the linear Henry's law regime, equilibrium dissolved gas content is directly proportional to initial charging pressure: n_d = S · P_0, where the solubility slope S = K_H · V_L = 0.0340 mol/bar. For CO_2 dissolving into 1.00 L of water at 25.00 °C, K_H = 0.0340 mol/(L·bar) (equivalent to Henry's volatility constant H = 1630 bar). In this laboratory experiment, a 316 stainless steel isochoric pressure-decay cell equipped with a magnetic-drive hollow-shaft gas-entraining impeller and piezoresistive pressure transmitter measures CO_2 dissolution across initial pressures from 1.0 to 6.0 bar to evaluate n_d and determine S = 0.0340 mol/bar.

n_d = S · P_0 = 0.0340 · P_0

Phase Equilibrium Thermodynamics of Dilute Gas-Liquid Solutions

Phase equilibrium between a gas phase ($V$) and a liquid solution ($L$) containing dissolved gas solute $i$ requires equality of chemical potentials:

\[ \mu_i^V(T, P, y_i) = \mu_i^L(T, P, x_i) \implies \hat{f}_i^V(T, P, y_i) = \hat{f}_i^L(T, P, x_i) \]

In the gas phase, fugacity is related to total pressure and vapor composition by the fugacity coefficient $\hat{\Phi}_i^V$:

\[ \hat{f}_i^V = y_i \, P \, \hat{\Phi}_i^V \]

In the liquid phase, for a sparingly soluble gas solute ($x_i \to 0$), fugacity is represented using Henry's law reference state at infinite dilution:

\[ \hat{f}_i^L = x_i \, \gamma_i^* \, H_{i,s}(T, P) \]

where $H_{i,s}$ is Henry's law constant of solute $i$ in solvent $s$, and $\gamma_i^*$ is the asymmetric activity coefficient ($\lim_{x_i \to 0} \gamma_i^* = 1.0$). At moderate pressures ($P \le 10\text{ bar}$), $\hat{\Phi}_i^V \approx 1.00$, $\gamma_i^* \approx 1.00$, and the Poynting correction is unity, producing Henry's law:

\[ P_i = H_{i,s} \, x_i \]

In terms of volumetric liquid molar concentration $C_d = n_d / V_L$:

\[ C_d = K_H \, P_i \implies n_d = K_H \, V_L \, P_i \]

Isochoric Pressure-Decay Mass Conservation

In an isochoric autoclave of total volume $V_{\text{cell}}$ containing liquid volume $V_L$ and head-space gas volume $V_G = V_{\text{cell}} - V_L$, pure gas is initially charged to pressure $P_0$. Initial gas moles charged:

\[ n_0 = \frac{P_0 \, V_G}{z_0 \, R T} \]

As gas dissolves, head-space pressure decays to equilibrium pressure $P_{\text{eq}}$. The final equilibrium gas moles remaining in the head-space is:

\[ n_{\text{eq}} = \frac{P_{\text{eq}} \, V_G}{z_{\text{eq}} \, R T} \]

By total mass conservation within the closed vessel, the equilibrium dissolved gas content $n_d$ is:

\[ n_d = n_0 - n_{\text{eq}} = \frac{V_G}{R T} \left( \frac{P_0}{z_0} - \frac{P_{\text{eq}}}{z_{\text{eq}}} \right) \]

Linear Solubility Regression & Henry's Parameter

Because the equilibrium partial pressure $P_{\text{eq}}$ is related to initial pressure $P_0$ by the system volume ratio and solubility constant ($P_{\text{eq}} = P_0 / [1 + K_H R T (V_L / V_G)]$), the dissolved gas moles $n_d$ scales linearly with initial charging pressure $P_0$:

\[ n_d = S \cdot P_0 = (K_H \cdot V_L) \cdot P_0 = 0.0340 \cdot P_0 \]

Plotting dissolved moles $n_d$ versus initial pressure $P_0$ yields a straight line through the origin with slope $S = 0.0340\text{ mol/bar}$, confirming Henry's law solubility constant $K_H = 0.0340\text{ mol/(L}\cdot\text{bar)}$ for $V_L = 1.00\text{ L}$.

Viva voce questions & analytical theory

How does temperature affect the solubility of gases in liquids?

For most common gases (CO_2, O_2, N_2) dissolving in water, dissolution is an exothermic process (ΔH_sol < 0); according to Le Chatelier's principle and the van 't Hoff equation d(ln K_H)/dT = -ΔH_sol/(R T^2), increasing temperature decreases gas solubility.

What is the difference between Henry's volatility constant H_i and solubility constant K_H?

Henry's volatility constant H_i (bar) relates partial pressure to mole fraction via P_i = H_i · x_i (higher H means lower solubility); Henry's solubility constant K_H (mol/(L·bar)) relates concentration to partial pressure via C_i = K_H · P_i (higher K_H means higher solubility); they are reciprocally related by K_H = 1 / (H_i · V_m,solvent).

Why does a hollow-shaft gas-entraining impeller accelerate equilibrium in pressure-decay cells?

The high-speed rotation of the hollow impeller generates a localized low-pressure zone at the impeller blades, sucking gas from the head-space through the hollow shaft and dispersing it as millions of microbubbles, creating massive interfacial contact area a (m^2/m^3).

Under what conditions does Henry's Law fail?

Henry's Law fails at very high pressures (> 20–50 bar) where gas fugacity deviates strongly from pressure and liquid non-ideality occurs; for highly soluble gases that chemically react and ionize extensively (e.g. HCl or NH_3); and when solute liquid concentrations exceed ~2–3 mol%.

How does the chemical hydration of dissolved CO_2 affect its apparent solubility?

Dissolved molecular CO_2(aq) undergoes reversible hydration to carbonic acid: CO_2(aq) + H_2O <-> H_2CO_3 <-> H^+ + HCO_3^-; however, at 25 °C and neutral pH, over 99.8% remains as physical dissolved molecular CO_2(aq), making Henry's physical law strictly valid.

How does water salinity ('salting-out' effect) alter Henry's law constant?

Dissolved ionic salts (e.g. NaCl) attract polar water molecules into tight hydration shells, reducing the concentration of free water molecules available to dissolve nonpolar or weakly polar gas molecules; this increases the gas activity coefficient and decreases solubility (described by the Setchenow equation).

Experimental precautions and process safety

  • Ensure the autoclave closure bolts are tightened in a diametrically opposed cross pattern with a calibrated torque wrench before pressurization.
  • Verify all dissolved atmospheric air (O_2 and N_2) is removed by vacuum boiling prior to introducing CO_2 to prevent mixed gas partial pressure errors.
  • Operate within the pressure vessel maximum allowable working pressure (MAWP = 10.0 MPa), and verify pressure relief rupture disc integrity.
  • Maintain stirring at constant 1000 RPM until pressure remains strictly flat for at least 15 minutes to guarantee complete phase equilibrium.

Theoretical assumptions & boundary condition limits

  • The transport and thermodynamic models apply within the specified calibrated operating window (1.0 to 6.0 bar).
  • Assumes ideal continuum flow and homogeneous phase mixing without unmodeled parasitic wall losses.
  • Continuous sensor readings reflect ordinary least-squares regression across successive observations.
  • High-temperature or extreme concentration excursions require non-ideal equation-of-state corrections.
↳
Connect transport theory to empirical process plant data.

Evaluated using a 2.0 L 316SS Buchi high-pressure autoclave reactor with hollow-shaft gas-entrainment Rushton impeller (1000 RPM) and Rosemount 3051S precision pressure transmitter (±0.005 bar). Pure CO_2 (99.99%) dissolving into 1.000 L degassed deionized water at 25.000 ± 0.01 °C.