# 3.8: Thermodynamics Applied to Fluids

- Page ID
- 18787

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In the rest of this chapter we apply the general results of section 3.6, “The Thermodynamic Dance”, to particular concrete situations. We begin with fluids, that is, systems whose thermodynamic states are adequately described by giving the variables temperature *T*, volume *V*, and particle number *N*. Furthermore, in this section we will not allow particles to enter or leave our system, so only two variables, *T* and *V*, are needed.

*3.8.1 Heat capacities *

Is there any relation between

\[ C_{p}(T, p)=T \frac{\partial S}{\partial T} )_{p},\]

the heat capacity at constant pressure, and

\[ C_{V}(T, V)=T \frac{\partial S}{\partial T} )_{V},\]

the heat capacity at constant volume? Remembering that entirely different experiments are used to measure the two quantities, you might guess that there is not. But mathematically, the difference between *C _{p}* and

*C*is related to a change in variable from (

_{V}*T, V*) to (

*T, p*), so you might begin to suspect a relation. [You might want to prepare for this section by working problem 3.12.]

Begin with the mathematical definition of the total differential of entropy,

\[ d S=\frac{\partial S}{\partial T} )_{V} d T+\frac{\partial S}{\partial V} )_{T} d V.\]

The above holds for any infinitesimal change. We restrict ourselves to infinitesimal changes at constant pressure, and divide by *dT*, to find

\[ \frac{\partial S}{\partial T} )_{p}=\frac{\partial S}{\partial T} )_{V}+\frac{\partial S}{\partial V} )_{T} \frac{\partial V}{\partial T} )_{p}.\]

Multiplying both sides by *T* gives a heat capacity relation

\[ C_{p}=C_{V}+T \frac{\partial S}{\partial V} )_{T} \frac{\partial V}{\partial T} )_{p}.\]

This relationship is correct but not yet written in its most convenient form. For example, reference books provide tabulations of *C _{p}* and

*C*, but not of ∂

_{V}*S*/∂

*V*)

_{T}or ∂

*V*/∂

*T*)

_{p}.

The first step in changing the heat capacity relation into “standard form” is an easy one. Recall from problem 1.2 that the expansion coefficient (an intensive tabulated quantity^{6}) is defined by

\[ \beta(T, p)=\frac{1}{V} \frac{\partial V}{\partial T} )_{p}\]

whence

\[ C_{p}=C_{V}+T V \beta \frac{\partial S}{\partial V} )_{T}.\]

The second step is less obvious. At equation (3.71) I discussed the Maxwell relation

\[ \frac{\partial S}{\partial V} )_{T}=\frac{\partial p}{\partial T} )_{V},\]

and mentioned that the right hand side was much easier to measure. Thus we have the preferable expression

\[C_{p}=C_{V}+T V \beta \frac{\partial p}{\partial T} )_{V}.\]

This expression is still not very convenient, however, because the quantity ∂*p*/∂*T*)_{V} is neither named nor tabulated. Although easier to measure than ∂*S*/∂*V*)_{T} is, its measurement still requires a constant-volume strong box. Measurements at constant pressure are easier to perform and hence more frequently tabulated. We now write ∂*p*/∂*T*)_{V} in terms of such quantities.

The total differential of *V*(*T, p*) is

\[ d V=\frac{\partial V}{\partial T} )_{p} d T+\frac{\partial V}{\partial p} )_{T} d p\]

or, in terms of named quantities,

\[ d V=V \beta d T-V \kappa_{T} d p.\]

(The isothermal compressibility

\[ \kappa_{T}(T, V) \equiv-\frac{1}{V} \frac{\partial V}{\partial p} )_{T}\]

was defined in problem 1.2.) Restricting the above total differential to changes at constant V gives

\[ V \beta dT = V \kappa_T dp\]

or

\[ \frac{\beta}{\kappa_{T}}=\frac{\partial p}{\partial T} )_{V}.\]

This immediately gives us the final form for our relationship between heat capacities,

\[ C_{p}=C_{V}+T V \frac{\beta^{2}}{\kappa_{T}}.\]

This is a far from obvious result that is true for all fluids and, I emphasize, was derived assuming only that entropy exists and without the use of any explicit formula for the entropy.

Note that *T*, *V*, and *κ _{T}* are all positive quantities. The expansion coefficient

*β*may be either positive or negative, but it enters the relationship as

*β*. Thus

^{2}\[ C_{p} \geq C_{V}.\]

Can we understand these results physically? Heat capacity is the amount of heat required to raise the temperature of a sample by one kelvin. If the sample is heated at constant volume, all of the heat absorbed goes into increasing the temperature. But if the sample is heated at constant pressure, then generally the substance will expand as well as increase its temperature, so some of the heat absorbed goes into increasing the temperature and some is converted into expansion work. Thus we expect that more heat will be needed for a given change of temperature at constant pressure than at constant volume, i.e. that *C _{p}* ≥

*C*. We also expect that the relation between

_{V}*C*and

_{p}*C*will depend upon the expansion coefficient

_{V}*β*. It is hard, however, to turn these qualitative observations into the precise formula (3.117) without recourse to the formal mathematics that we have used. It is surprising, for example, that

*β*enters into the formula as

*β*

^{2}, so that materials which contract with increasing temperature still have

*C*≥

_{p}*C*.

_{V}For solids and liquids, there is little expansion and thus little expansion work, whence *C*_{p }≈ *C _{V}*. But for gases, \( C_{p} \gg C_{V}\). What happens at two-phase coexistence?

So far we have discussed only the difference between the two heat capacities. It is conventional to also define the “heat capacity ratio”

\[ \gamma \equiv \frac{C_{p}}{C_{V}},\]

which can be regarded as a function of *T* and *V* or as a function of *T* and *p*. But the result above shows that for all values of *T* and *V*,

\[ \gamma(T, V) \geq 1.\]

We can apply these results to the ideal gas with equation of state

\[ pV = Nk_B T.\]

Problem 1.2 used this equation to show that for an ideal gas,

\[ \beta=\frac{1}{T} \quad \text { and } \quad \kappa_{T}=\frac{1}{p}.\]

Thus the heat capacity relation (3.117) becomes

\[ C_{p}=C_{V}+T V \frac{p}{T^{2}}=C_{V}+\frac{p V}{T}=C_{V}+k_{B} N,\]

and we have

\[ \gamma=\frac{C_{p}}{C_{V}}=1+\frac{k_{B} N}{C_{V}}.\]

## 3.8.2 Energy as a function of temperature and volume

You know that *E*(*S, V*) is a master function with the famous master equation

\[dE = T dS − p dV.\]

For the variables *T* and *V* the master function is the Helmholtz potential F(*T, V*) with master equation

\[dF = −S dT − p dV.\]

Thus in terms of the variables *T* and *V* the energy in no longer a master function, but that doesn’t mean that it’s not a function at all. It is possible to consider the energy as a function of *T* and *V*, and in this subsection we will find the total differential of *E*(*T, V*). We begin by finding the total differential of *S*(*E, V*) and finish off by substituting that expression for *dS* into the master equation (3.125). The mathematical expression for that total differential is

\[ d S=\frac{\partial S}{\partial T} )_{V} d T+\frac{\partial S}{\partial V} )_{T} d V,\]

but we showed in the previous subsection that this is more conveniently written as

\[ d S=\frac{C_{V}}{T} d T+\frac{\beta}{\kappa_{T}} d V.\]

Thus

\[ d E=T d S-p d V=T\left[\frac{C_{V}}{T} d T+\frac{\beta}{\kappa_{T}} d V\right]-p d V\]

and, finally,

\[ d E=C_{V} d T+\left[T \frac{\beta}{\kappa_{T}}-p\right] d V.\]

The most interesting consequence of this exercise is that the heat capacity *C _{V}*, which was defined in terms of an entropy derivative, is also equal to an energy derivative:

\[ C_{V}(T, V) \equiv T \frac{\partial S}{\partial T} )_{V}=\frac{\partial E}{\partial T} )_{V}.\]

For an ideal gas

\[ \left[T \frac{\beta}{\kappa_{T}}-p\right]=0,\]

whence

\[ d E=C_{V} d T+0 d V.\]

Thus the energy, which for most fluids depends upon both temperature and volume, depends in the ideal gas only upon temperature:

\[ E(T, V ) = E(T).\]

and the heat capacity ratio

\[ \gamma(T)=1+\frac{k_{B} N}{C_{V}(T)}.\]

Thermodynamics proves that for an ideal gas *C _{V}* depends only upon

*T*and not

*V*, but experiment shows that in fact

*C*is almost independent of

_{V}*T*as well. From the point of view of thermodynamics this is just an unexpected experimental result and nothing deeper. From the point of view of statistical mechanics this is not just a coincidence but is something that can be proved. . . we will do so in chapter 5.

## 3.8.3 Quasistatic adiabatic processes

This topic, titled by a string of polysyllabic words, sounds like an arcane and abstract one. In fact it has direct applications to the speed of sound, the manufacture of refrigerators, to the inflation of soccer balls, and to the temperature of mountain tops. [You might want to prepare for this section by working problem 3.13.]

If a fluid is changed in a quasistatic, adiabatic manner, its entropy will remain constant during the change. So, for example, if the state of the system is described by the two variables temperature and volume, then during a quasistatic adiabatic process the system will move along contour lines^{7} of constant entropy in the (*T, V*) plane, such as those below.

Our aim in this section is to find an equation V(T) for such constant-entropy contours.

We begin with the now-familiar total differential for the entropy as a function of temperature and volume,

\[ d S=\frac{C_{V}}{T} d T+\frac{\beta}{\kappa_{T}} d V.\]

This equation holds for any infinitesimal change, but we will apply it now to changes along an entropy contour, i.e. changes for which *dS* = 0. For such changes

\[ \frac{\beta}{\kappa_{T}} d V=-\frac{C_{V}}{T} d T\]

whence the differential equation for an entropy contour is

\[ \frac{d V(T)}{d T}=-\frac{C_{V}(T, V) \kappa_{T}(T, V)}{T \beta(T, V)}.\]

The quantities on the right are usually (not always) positive, so the entropy contours usually (not always) slope downward. This is as much as we can say for a general fluid.

But for an ideal gas we can fill in the functions on the right to find

\[ \frac{d V(T)}{d T}=-\frac{C_{V}(T)}{p(T, V)}=-\frac{C_{V}(T) V}{N k_{B} T}.\]

Using equation (3.136), this result is often written in the form

\[ \frac{d V(T)}{d T}=\frac{1}{1-\gamma(T)} \frac{V}{T}.\]

This is as far as we can go for an arbitrary ideal gas.

But for an ideal gas with constant heat capacity, and hence with constant γ, this differential equation is readily solved using separation of variables. We write

\[ (1-\gamma) \frac{d V}{V}=\frac{d T}{T}\]

with the immediate solution

\[ (1-\gamma) \ln V=\ln T+\text { constant, }\]

giving

\[ V^{1-\gamma}=\text { constant } T \quad \text { or } \quad T V^{\gamma-1}=\text { constant. }\]

Remembering that γ ≥ 1, this confirms the downward slope of the entropy contours on the (*V, T*) plane. This result is most frequently used with the variables *p* and *T*, rather than *V* and *T*:

\[ p V^{\gamma}=\text { constant }.\]

The above equation holds only for quasistatic, adiabatic processes in ideal gases that have heat capacities independent of temperature. You might think that it is such a restricted result that it holds little interest. This equation is used frequently in applications, in the graduate record examination, and in physics oral examinations! You should memorize it and know how to work with it.

## Problems

*3.29 Intensive vs. extensive variables*

Equation (3.134), *E*(*T, V*) = *E*(*T*), states that for an ideal gas the energy is a function of temperature alone. How is it possible for *E*, an extensive quantity, to be a function of only *T*, an intensive quantity?

*3.30 Heat capacity at constant pressure *

Equation (3.131) shows that the heat capacity at constant volume, which is defined in terms of an entropy derivative, is also equal to an energy derivative. You might suspect a similar relation between *C _{p}*(

*T, p, N*) and

\( \frac{\partial E(T, p, N)}{\partial T} )_{p, N}.\)

Show that such a suspicion is not correct, and instead find an expression for *C _{p}* in terms of a derivative of enthalpy.

*3.31 Qualitative heat capacity relations *

(This problem is stolen from a GRE Physics test.)

For an ideal gas, *C _{p}* is greater than

*C*because:

_{V}a. The gas does work on its surroundings when its pressure remains constant while its temperature is increased.

b. The heat input per unit increase in temperature is the same for processes at constant pressure and at constant volume.

c. The pressure of the gas remains constant when its temperature remains constant.

d. The increase in the gas’s internal energy is greater when the pressure remains constant than when the volume remains constant.

e. The heat needed is greater when the volume remains constant than when the pressure remains constant.

*3.32 Heat capacities in a magnetic system *

For a magnetic system (see equation (3.100)), show that

\[ C_{H}=T \frac{\partial S}{\partial T} )_{H}, \quad C_{M}=T \frac{\partial S}{\partial T} )_{M}, \quad \beta=\frac{\partial M}{\partial T} )_{H}, \quad \text { and } \quad \chi_{T}=\frac{\partial M}{\partial H} )_{T}\]

are related through

\[ C_{M}=C_{H}-T \beta^{2} / \chi_{T}\]

*3.33 Isothermal compressibility *

a. Show that the isothermal compressibility, defined in problem 1.2 as

\[ \kappa_{T}=-\frac{1}{V} \frac{\partial V}{\partial p} )_{T, N}\]

is also given by

\[ \kappa_{T}=\frac{1}{\rho} \frac{\partial \rho}{\partial p} )_{T}=\frac{1}{\rho^{2}} \frac{\partial \rho}{\partial \mu} )_{T}=\frac{1}{\rho^{2}} \frac{\partial^{2} p}{\partial \mu^{2}} )_{T},\]

where *ρ* is the number density *N/V*. Clue:

\[ \frac{\partial \rho}{\partial p} )_{T}=\frac{\partial \rho}{\partial \mu} )_{T} \frac{\partial \mu}{\partial p} )_{T}.\]

b. What does this result tell you about the relation between density and chemical potential?

c. In part (a.) we began with a description in terms of the three variables *T, p, N* and then reduced it to an intensive-only description, which requires just two variables, such as *µ* and *T*. Reverse this process to show that

\[ \kappa_{T}=\frac{V}{N^{2}} \frac{\partial N}{\partial \mu} )_{T, V}=-\frac{1}{N} \frac{\partial V}{\partial \mu} )_{T, N}.\]

*3.34 Pressure differential *

By regarding pressure as a function of temperature *T* and number density *ρ*, show that

\[ d p=\frac{\beta}{\kappa_{T}} d T+\frac{1}{\rho \kappa_{T}} d \rho.\]

*3.35 Isothermal vs. adiabatic compressibility *

In class we derived a remarkable relation between the heat capacities *C _{p}* and

*C*. This problem uncovers a similar relation between the isothermal and adiabatic compressibilities,

_{V}\[ \kappa_{T}=-\frac{1}{V} \frac{\partial V}{\partial p} )_{T} \quad \text { and } \quad \kappa_{S}=-\frac{1}{V} \frac{\partial V}{\partial p} )_{S}.\]

The adiabatic compressibility *κ _{S}* is the compressibility measured when the fluid is thermally insulated. (It is related to the speed of sound: see problem 3.44.)

a. Use

\[ d S=\frac{\partial S}{\partial T} )_{p} d T+\frac{\partial S}{\partial p} )_{T} d p\]

to show that

\[ \frac{\partial T}{\partial p} )_{S}=\frac{\beta T}{C_{p} / V}.\]

Sketch an experiment to measure this quantity directly.

b. From the mathematical relation

\[ d V=\frac{\partial V}{\partial p} )_{T} d p+\frac{\partial V}{\partial T} )_{p} d T\]

derive the multivariate chain rule

\[ \frac{\partial V}{\partial p} )_{S}=\frac{\partial V}{\partial p} )_{T}+\frac{\partial V}{\partial T} )_{p} \frac{\partial T}{\partial p} )_{S},\]

whence

\[\kappa_{S}=\kappa_{T}-\frac{\beta^{2} T}{C_{p} / V}.\]

c. Finally, show that

\[ \gamma \equiv \frac{C_{p}}{C_{V}}=\frac{\kappa_{T}}{\kappa_{S}}.\]

*3.36 Change of chemical potential with temperature *

Prove that

\[ \frac{\partial \mu}{\partial T} )_{p, N}=-\frac{S}{N}\]

and that

\[ \frac{\partial \mu}{\partial T} )_{V, N}=-\frac{\partial S}{\partial N} )_{T, V}=-\frac{S}{N}+\frac{\beta}{\rho \kappa_{T}}.\]

How's that for weird?

^{6}It is easy to see why ∂*V*/∂*T*)_{p} itself is not tabulated. Instead of requiring a tabulation for iron, you would need a tabulation for a five gram sample of iron, for a six gram sample of iron, for a seven gram sample of iron, etc.

^{7}Such contours are sometimes called isoentropic curves or, even worse, adiabats. The latter name is not just a linguistic horror, it also omits the essential qualifier that only quasistatic adiabatic processes move along curves of constant entropy