Physics:Joule–Thomson effect

Short description: Phenomenon of non-ideal fluids changing temperature while being forced through small spaces

In thermodynamics, the Joule–Thomson effect (also known as the Joule–Kelvin effect or Kelvin–Joule effect) describes the temperature change of a real gas or liquid (as differentiated from an ideal gas) when it is forced through a valve or porous plug while keeping it insulated so that no heat is exchanged with the environment. This procedure is called a throttling process or Joule–Thomson process. At room temperature, all gases except hydrogen, helium, and neon cool upon expansion by the Joule–Thomson process when being throttled through an orifice; these three gases experience the same effect but only at lower temperatures. Most liquids such as hydraulic oils will be warmed by the Joule–Thomson throttling process.

The gas-cooling throttling process is commonly exploited in refrigeration processes such as liquefiers in air separation industrial process. In hydraulics, the warming effect from Joule–Thomson throttling can be used to find internally leaking valves as these will produce heat which can be detected by thermocouple or thermal-imaging camera. Throttling is a fundamentally irreversible process. The throttling due to the flow resistance in supply lines, heat exchangers, regenerators, and other components of (thermal) machines is a source of losses that limits their performance.

Since it is a constant-enthalpy process, it can be used to experimentally measure the lines of constant enthalpy (isenthalps) on the [math]\displaystyle{ (p, T) }[/math] diagram of a gas. Combined with the specific heat capacity at constant pressure [math]\displaystyle{ c_P = (\partial h / \partial T)_P }[/math] it allows the complete measurement of the thermodynamic potential for the gas.

History

The effect is named after James Prescott Joule and William Thomson, 1st Baron Kelvin, who discovered it in 1852. It followed upon earlier work by Joule on Joule expansion, in which a gas undergoes free expansion in a vacuum and the temperature is unchanged, if the gas is ideal.

Description

The adiabatic (no heat exchanged) expansion of a gas may be carried out in a number of ways. The change in temperature experienced by the gas during expansion depends not only on the initial and final pressure, but also on the manner in which the expansion is carried out.

• If the expansion process is reversible, meaning that the gas is in thermodynamic equilibrium at all times, it is called an isentropic expansion. In this scenario, the gas does positive work during the expansion, and its temperature decreases.
• In a free expansion, on the other hand, the gas does no work and absorbs no heat, so the internal energy is conserved. Expanded in this manner, the temperature of an ideal gas would remain constant, but the temperature of a real gas decreases, except at very high temperature.
• The method of expansion discussed in this article, in which a gas or liquid at pressure P₁ flows into a region of lower pressure P₂ without significant change in kinetic energy, is called the Joule–Thomson expansion. The expansion is inherently irreversible. During this expansion, enthalpy remains unchanged (see proof below). Unlike a free expansion, work is done, causing a change in internal energy. Whether the internal energy increases or decreases is determined by whether work is done on or by the fluid; that is determined by the initial and final states of the expansion and the properties of the fluid.

Sign of the Joule–Thomson coefficient, [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] for N₂. Within the region bounded by the red line, a Joule–Thomson expansion produces cooling ([math]\displaystyle{ \mu_{\mathrm{JT}} \gt 0 }[/math]); outside that region, the expansion produces heating. The gas–liquid coexistence curve is shown by the blue line, terminating at the critical point (the solid blue circle). The dashed lines demarcate the region where N₂ is a supercritical fluid (where properties smoothly transition between liquid-like and gas-like).

The temperature change produced during a Joule–Thomson expansion is quantified by the Joule–Thomson coefficient, [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math]. This coefficient may be either positive (corresponding to cooling) or negative (heating); the regions where each occurs for molecular nitrogen, N₂, are shown in the figure. Note that most conditions in the figure correspond to N₂ being a supercritical fluid, where it has some properties of a gas and some of a liquid, but can not be really described as being either. The coefficient is negative at both very high and very low temperatures; at very high pressure it is negative at all temperatures. The maximum inversion temperature (621 K for N₂) occurs as zero pressure is approached. For N₂ gas at low pressures, [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] is negative at high temperatures and positive at low temperatures. At temperatures below the gas-liquid coexistence curve, N₂ condenses to form a liquid and the coefficient again becomes negative. Thus, for N₂ gas below 621 K, a Joule–Thomson expansion can be used to cool the gas until liquid N₂ forms.

Physical mechanism

There are two factors that can change the temperature of a fluid during an adiabatic expansion: a change in internal energy or the conversion between potential and kinetic internal energy. Temperature is the measure of thermal kinetic energy (energy associated with molecular motion); so a change in temperature indicates a change in thermal kinetic energy. The internal energy is the sum of thermal kinetic energy and thermal potential energy. Thus, even if the internal energy does not change, the temperature can change due to conversion between kinetic and potential energy; this is what happens in a free expansion and typically produces a decrease in temperature as the fluid expands. If work is done on or by the fluid as it expands, then the total internal energy changes. This is what happens in a Joule–Thomson expansion and can produce larger heating or cooling than observed in a free expansion.

In a Joule–Thomson expansion the enthalpy remains constant. The enthalpy, [math]\displaystyle{ H }[/math], is defined as

[math]\displaystyle{ H = U + PV }[/math]

where [math]\displaystyle{ U }[/math] is internal energy, [math]\displaystyle{ P }[/math] is pressure, and [math]\displaystyle{ V }[/math] is volume. Under the conditions of a Joule–Thomson expansion, the change in [math]\displaystyle{ PV }[/math] represents the work done by the fluid (see the proof below). If [math]\displaystyle{ PV }[/math] increases, with [math]\displaystyle{ H }[/math] constant, then [math]\displaystyle{ U }[/math] must decrease as a result of the fluid doing work on its surroundings. This produces a decrease in temperature and results in a positive Joule–Thomson coefficient. Conversely, a decrease in [math]\displaystyle{ PV }[/math] means that work is done on the fluid and the internal energy increases. If the increase in kinetic energy exceeds the increase in potential energy, there will be an increase in the temperature of the fluid and the Joule–Thomson coefficient will be negative.

For an ideal gas, [math]\displaystyle{ PV }[/math] does not change during a Joule–Thomson expansion. As a result, there is no change in internal energy; since there is also no change in thermal potential energy, there can be no change in thermal kinetic energy and, therefore, no change in temperature. In real gases, [math]\displaystyle{ PV }[/math] does change.

The ratio of the value of [math]\displaystyle{ PV }[/math] to that expected for an ideal gas at the same temperature is called the compressibility factor, [math]\displaystyle{ Z }[/math]. For a gas, this is typically less than unity at low temperature and greater than unity at high temperature (see the discussion in compressibility factor). At low pressure, the value of [math]\displaystyle{ Z }[/math] always moves towards unity as a gas expands. Thus at low temperature, [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ PV }[/math] will increase as the gas expands, resulting in a positive Joule–Thomson coefficient. At high temperature, [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ PV }[/math] decrease as the gas expands; if the decrease is large enough, the Joule–Thomson coefficient will be negative.

For liquids, and for supercritical fluids under high pressure, [math]\displaystyle{ PV }[/math] increases as pressure increases. This is due to molecules being forced together, so that the volume can barely decrease due to higher pressure. Under such conditions, the Joule–Thomson coefficient is negative, as seen in the figure above.

The physical mechanism associated with the Joule–Thomson effect is closely related to that of a shock wave, although a shock wave differs in that the change in bulk kinetic energy of the gas flow is not negligible.

The Joule–Thomson (Kelvin) coefficient

Fig. 1 – Joule–Thomson coefficients for various gases at atmospheric pressure

The rate of change of temperature [math]\displaystyle{ T }[/math] with respect to pressure [math]\displaystyle{ P }[/math] in a Joule–Thomson process (that is, at constant enthalpy [math]\displaystyle{ H }[/math]) is the Joule–Thomson (Kelvin) coefficient [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math]. This coefficient can be expressed in terms of the gas's volume [math]\displaystyle{ V }[/math], its heat capacity at constant pressure [math]\displaystyle{ C_{\mathrm{p}} }[/math], and its coefficient of thermal expansion [math]\displaystyle{ \alpha }[/math] as:

[math]\displaystyle{ \mu_{\mathrm{JT}} = \left( {\partial T \over \partial P} \right)_H = \frac V {C_{\mathrm{p}}}(\alpha T - 1)\, }[/math]

See the § Derivation of the Joule–Thomson coefficient below for the proof of this relation. The value of [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] is typically expressed in °C/bar (SI units: K/Pa) and depends on the type of gas and on the temperature and pressure of the gas before expansion. Its pressure dependence is usually only a few percent for pressures up to 100 bar.

All real gases have an inversion point at which the value of [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] changes sign. The temperature of this point, the Joule–Thomson inversion temperature, depends on the pressure of the gas before expansion.

In a gas expansion the pressure decreases, so the sign of [math]\displaystyle{ \partial P }[/math] is negative by definition. With that in mind, the following table explains when the Joule–Thomson effect cools or warms a real gas:

If the gas temperature is	then [math]\displaystyle{ \mu_\text{JT} }[/math] is	since [math]\displaystyle{ \partial P }[/math] is	thus [math]\displaystyle{ \partial T }[/math] must be	so the gas
below the inversion temperature	positive	always negative	negative	cools
above the inversion temperature	negative		positive	warms

Helium and hydrogen are two gases whose Joule–Thomson inversion temperatures at a pressure of one atmosphere are very low (e.g., about 40 K, −233 °C for helium). Thus, helium and hydrogen warm when expanded at constant enthalpy at typical room temperatures. On the other hand, nitrogen and oxygen, the two most abundant gases in air, have inversion temperatures of 621 K (348 °C) and 764 K (491 °C) respectively: these gases can be cooled from room temperature by the Joule–Thomson effect.

For an ideal gas, [math]\displaystyle{ \mu_\text{JT} }[/math] is always equal to zero: ideal gases neither warm nor cool upon being expanded at constant enthalpy.

Theoretical models

For the Dieterici gas, the relation between reduced pressure and reduced inversion temperature is [math]\displaystyle{ \tilde p = (8-\tilde T_I) e^{\frac 52 - \frac 4{8-\tilde T_I}} }[/math]. Reproduced from Fig. 17 of .

For a Van der Waals gas, the coefficient is[math]\displaystyle{ \mu_\text{JT}=-\frac{V_m}{C_{p}} \frac{R T V_m^{2} b-2 a(V_m-b)^{2}}{R T V_m^{3}-2 a(V_m-b)^{2}}. }[/math]with inversion temperature [math]\displaystyle{ \frac{2a}{bR}\left(1 - \frac{b}{V_m}\right)^2 }[/math].

For the Dieterici gas, the reduced inversion temperature is [math]\displaystyle{ \tilde T_I = 8 - 4/\tilde V_m }[/math], and the relation between reduced pressure and reduced inversion temperature is [math]\displaystyle{ \tilde p = (8-\tilde T_I) e^{\frac 52 - \frac 4{8-\tilde T_I}} }[/math]. This is plotted on the right. The critical point falls inside the region where the gas cools on expansion. The outside region is where the gas warms on expansion.

Applications

In practice, the Joule–Thomson effect is achieved by allowing the gas to expand through a throttling device (usually a valve) which must be very well insulated to prevent any heat transfer to or from the gas. No external work is extracted from the gas during the expansion (the gas must not be expanded through a turbine, for example).

The cooling produced in the Joule–Thomson expansion makes it a valuable tool in refrigeration. The effect is applied in the Linde technique as a standard process in the petrochemical industry, where the cooling effect is used to liquefy gases, and in many cryogenic applications (e.g. for the production of liquid oxygen, nitrogen, and argon). A gas must be below its inversion temperature to be liquefied by the Linde cycle. For this reason, simple Linde cycle liquefiers, starting from ambient temperature, cannot be used to liquefy helium, hydrogen, or neon. They must first be cooled to their inversion temperatures, which are -233 C (helium), -71 C (hydrogen), and -42 C (neon).

Proof that the specific enthalpy remains constant

In thermodynamics so-called "specific" quantities are quantities per unit mass (kg) and are denoted by lower-case characters. So h, u, and v are the specific enthalpy, specific internal energy, and specific volume (volume per unit mass, or reciprocal density), respectively. In a Joule–Thomson process the specific enthalpy h remains constant. To prove this, the first step is to compute the net work done when a mass m of the gas moves through the plug. This amount of gas has a volume of V₁ = m v₁ in the region at pressure P₁ (region 1) and a volume V₂ = m v₂ when in the region at pressure P₂ (region 2). Then in region 1, the "flow work" done on the amount of gas by the rest of the gas is: W₁ = m P₁v₁. In region 2, the work done by the amount of gas on the rest of the gas is: W₂ = m P₂v₂. So, the total work done on the mass m of gas is

[math]\displaystyle{ W = mP_1 v_1 - mP_2 v_2. }[/math]

The change in internal energy minus the total work done on the amount of gas is, by the first law of thermodynamics, the total heat supplied to the amount of gas.

[math]\displaystyle{ U - W = Q }[/math]

In the Joule–Thomson process, the gas is insulated, so no heat is absorbed. This means that

[math]\displaystyle{ \begin{align} (mu_2 - mu_1) &- (mP_1 v_1 - mP_2 v_2) = 0 \\ mu_1 + mP_1 v_1 &= mu_2 + mP_2 v_2 \\ u_1 + P_1 v_1 &= u_2 + P_2 v_2 \end{align} }[/math]

where u₁ and u₂ denote the specific internal energies of the gas in regions 1 and 2, respectively. Using the definition of the specific enthalpy h = u + Pv, the above equation implies that

[math]\displaystyle{ h_1 = h_2 }[/math]

where h₁ and h₂ denote the specific enthalpies of the amount of gas in regions 1 and 2, respectively.

Throttling in the T-s diagram

Fig. 2 – T-s diagram of nitrogen. The red dome represents the two-phase region with the low-entropy side (the saturated liquid) and the high-entropy side (the saturated gas). The black curves give the T-s relation along isobars. The pressures are indicated in bar. The blue curves are isenthalps (curves of constant specific enthalpy). The specific enthalpies are indicated in kJ/kg. The specific points a, b, etc., are treated in the main text.

A very convenient way to get a quantitative understanding of the throttling process is by using diagrams such as h-T diagrams, h-P diagrams, and others. Commonly used are the so-called T-s diagrams. Figure 2 shows the T-s diagram of nitrogen as an example. Various points are indicated as follows:

1. T = 300 K, p = 200 bar, s = 5.16 kJ/(kgK), h = 430 kJ/kg;
2. T = 270 K, p = 1 bar, s = 6.79 kJ/(kgK), h = 430 kJ/kg;
3. T = 133 K, p = 200 bar, s = 3.75 kJ/(kgK), h = 150 kJ/kg;
4. T = 77.2 K, p = 1 bar, s = 4.40 kJ/(kgK), h = 150 kJ/kg;
5. T = 77.2 K, p = 1 bar, s = 2.83 kJ/(kgK), h = 28 kJ/kg (saturated liquid at 1 bar);
6. T = 77.2 K, p = 1 bar, s = 5.41 kJ/(kgK), h = 230 kJ/kg (saturated gas at 1 bar).

As shown before, throttling keeps h constant. E.g. throttling from 200 bar and 300 K (point a in fig. 2) follows the isenthalpic (line of constant specific enthalpy) of 430 kJ/kg. At 1 bar it results in point b which has a temperature of 270 K. So throttling from 200 bar to 1 bar gives a cooling from room temperature to below the freezing point of water. Throttling from 200 bar and an initial temperature of 133 K (point c in fig. 2) to 1 bar results in point d, which is in the two-phase region of nitrogen at a temperature of 77.2 K. Since the enthalpy is an extensive parameter the enthalpy in d (h_d) is equal to the enthalpy in e (h_e) multiplied with the mass fraction of the liquid in d (x_d) plus the enthalpy in f (h_f) multiplied with the mass fraction of the gas in d (1 − x_d). So

[math]\displaystyle{ h_d = x_d h_e + (1 - x_d) h_f. }[/math]

With numbers: 150 = x_d 28 + (1 − x_d) 230 so x_d is about 0.40. This means that the mass fraction of the liquid in the liquid–gas mixture leaving the throttling valve is 40%.

Derivation of the Joule–Thomson coefficient

It is difficult to think physically about what the Joule–Thomson coefficient, [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math], represents. Also, modern determinations of [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] do not use the original method used by Joule and Thomson, but instead measure a different, closely related quantity. Thus, it is useful to derive relationships between [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] and other, more conveniently measured quantities, as described below.

The first step in obtaining these results is to note that the Joule–Thomson coefficient involves the three variables T, P, and H. A useful result is immediately obtained by applying the cyclic rule; in terms of these three variables that rule may be written

[math]\displaystyle{ \left(\frac{\partial T}{\partial P}\right)_H\left(\frac{\partial H}{\partial T}\right)_P \left(\frac{\partial P}{\partial H}\right)_T = -1. }[/math]

Each of the three partial derivatives in this expression has a specific meaning. The first is [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math], the second is the constant pressure heat capacity, [math]\displaystyle{ C_{\mathrm{p}} }[/math], defined by

[math]\displaystyle{ C_{\mathrm{p}} = \left(\frac{\partial H}{\partial T}\right)_P }[/math]

and the third is the inverse of the isothermal Joule–Thomson coefficient, [math]\displaystyle{ \mu_{\mathrm{T}} }[/math], defined by

[math]\displaystyle{ \mu_{\mathrm{T}} = \left(\frac{\partial H}{\partial P}\right)_T }[/math].

This last quantity is more easily measured than [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] . Thus, the expression from the cyclic rule becomes

[math]\displaystyle{ \mu_{\mathrm{JT}} = - \frac{\mu_{\mathrm{T}}} {C_p}. }[/math]

This equation can be used to obtain Joule–Thomson coefficients from the more easily measured isothermal Joule–Thomson coefficient. It is used in the following to obtain a mathematical expression for the Joule–Thomson coefficient in terms of the volumetric properties of a fluid.

To proceed further, the starting point is the fundamental equation of thermodynamics in terms of enthalpy; this is

[math]\displaystyle{ \mathrm{d}H = T \mathrm{d}S + V \mathrm{d}P. }[/math]

Now "dividing through" by dP, while holding temperature constant, yields

[math]\displaystyle{ \left(\frac{\partial H}{\partial P}\right)_T = T\left(\frac{\partial S}{\partial P}\right)_T + V }[/math]

The partial derivative on the left is the isothermal Joule–Thomson coefficient, [math]\displaystyle{ \mu_{\mathrm{T}} }[/math], and the one on the right can be expressed in terms of the coefficient of thermal expansion via a Maxwell relation. The appropriate relation is

[math]\displaystyle{ \left(\frac{\partial S}{\partial P}\right)_T= -\left(\frac{\partial V}{\partial T}\right)_P= -V\alpha\, }[/math]

where α is the cubic coefficient of thermal expansion. Replacing these two partial derivatives yields

[math]\displaystyle{ \mu_{\mathrm{T}} = - T V\alpha\ + V. }[/math]

This expression can now replace [math]\displaystyle{ \mu_{\mathrm{T}} }[/math] in the earlier equation for [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] to obtain:

[math]\displaystyle{ \mu_{\mathrm{JT}} \equiv \left( \frac{\partial T}{\partial P} \right)_H = \frac V {C_{\mathrm{p}}} (\alpha T - 1).\, }[/math]

This provides an expression for the Joule–Thomson coefficient in terms of the commonly available properties heat capacity, molar volume, and thermal expansion coefficient. It shows that the Joule–Thomson inversion temperature, at which [math]\displaystyle{ \mu_{\mathrm{JT}} }[/math] is zero, occurs when the coefficient of thermal expansion is equal to the inverse of the temperature. Since this is true at all temperatures for ideal gases (see expansion in gases), the Joule–Thomson coefficient of an ideal gas is zero at all temperatures.

Joule's second law

It is easy to verify that for an ideal gas defined by suitable microscopic postulates that αT = 1, so the temperature change of such an ideal gas at a Joule–Thomson expansion is zero. For such an ideal gas, this theoretical result implies that:

The internal energy of a fixed mass of an ideal gas depends only on its temperature (not pressure or volume).

This rule was originally found by Joule experimentally for real gases and is known as Joule's second law. More refined experiments found important deviations from it.

See also

• Critical point (thermodynamics)
• Enthalpy and Isenthalpic process
• Ideal gas
• Liquefaction of gases
• MIRI (Mid-Infrared Instrument), a J–T loop is used on one of the instruments of the James Webb Space Telescope
• Refrigeration
• Reversible process (thermodynamics)

Bibliography

• M. W. Zemansky (1968). Heat and Thermodynamics; An Intermediate Textbook. McGraw-Hill. pp. 182, 355. https://archive.org/details/heatthermodynami0000zema. 
• D. V. Schroeder (2000). An Introduction to Thermal Physics. Addison Wesley Longman. p. 142. ISBN 978-0-201-38027-9. https://archive.org/details/introductiontoth00schr_817. 
• C. Kittel, H. Kroemer (1980). Thermal Physics. W. H. Freeman. ISBN 978-0-7167-1088-2. 

External links

• Weisstein, Eric Wolfgang, ed. "Joule-Thomson process". http://scienceworld.wolfram.com/physics/Joule-ThomsonProcess.html. 
• Weisstein, Eric Wolfgang, ed. "Joule-Thomson coefficient". http://scienceworld.wolfram.com//physics/Joule-ThomsonCoefficient.html. 
• "Inversion Curve of Joule-Thomson Effect using Peng-Robinson CEOS". Demonstrations Projects of Wolfram Mathematica. http://demonstrations.wolfram.com/InversionCurveOfJouleThomsonUsingThePengRobinsonCubicEquatio/. 
• Joule–Thomson effect module, University of Notre Dame