The Chapman–Enskog Transport Coefficients for Maxwell Molecules

Contents
  1. The interaction, and why it is the one to take
  2. The linearised collision operator
  3. The shear eigenfunction
  4. The heat-flux eigenfunction
  5. The first-order solution
  6. The transport coefficients

This appendix proves Proposition 61.9 of Experiment: The Kinetic Theory Tested: that for a monatomic gas whose molecules repel as the inverse fifth power of their separation, the Boltzmann transport equation gives

\begin{equation}\tag{A69.1} \frac{\kappa m}{\eta\,c_{V}^{(1)}}=\frac{5}{2} \end{equation}

exactly, with \(c_{V}^{(1)}=\tfrac{3}{2}k_{\mathrm{B}}\) the heat capacity of one molecule. The elementary mean-free-path argument of Equation (61.10) gives \(1\) for this ratio and the relaxation-time model of Proposition 61.7 gives \(\tfrac{5}{3}\).

The proof needs no collision integral. The scattering enters through a single number \(A_{2}\), defined in Equation (A69.6) below, and that number cancels from the ratio — which is why Equation (A69.1) is a pure rational number and not a property of the strength of the force.

The interaction, and why it is the one to take

A Maxwell molecule repels its neighbour with a force \(K/r^{5}\). That is not a choice of convenience: it is the force law Remark 61.3 extracts from Maxwell's own measurement of how the viscosity of a gas varies with temperature, so the model solved here is the model that chapter tests.

Its distinguishing property is a similarity. For a repulsive force proportional to \(r^{-\nu}\) the deflection suffered in a collision depends on the impact parameter \(b\) and the relative speed \(g\) only through the combination \(b\,g^{2/(\nu-1)}\), so the differential cross-section is \(g^{-4/(\nu-1)}\) times a function of the deflection angle \(\chi\) alone. At \(\nu=5\) that exponent is \(-1\), and the product

\begin{equation}\tag{A69.2} g\,\sigma(g,\chi)=F(\chi) \end{equation}

is independent of the relative speed. Equation (A69.2) is the whole of what is used below; \(F\) carries SI dimension \(\mathrm{m}^{3}/\mathrm{s}\), being a speed times an area.

The linearised collision operator

Write the distribution as \(f=f_{0}(1+\Phi)\) with \(f_{0}\) the local Maxwellian built from the local density \(n\), mean velocity \(\vect{u}\) and temperature \(T\), and let \(\vect{C}=\vect{v}-\vect{u}\) be the peculiar velocity. Because a Maxwellian satisfies detailed balance, \(f_{0}'f_{0,1}'=f_{0}f_{0,1}\), the collision integral Equation (59.15) linearises to \(f_{0}\,L[\Phi]\) with

\begin{equation}\tag{A69.3} L[\varphi](\vect{v}) =\int\!\dd^{3}v_{1}\,f_{0}(\vect{v}_{1})\!\int\!\dd\Omega\; g\,\sigma\, \left(\varphi'+\varphi_{1}'-\varphi-\varphi_{1}\right)\ep \end{equation}

Every collision invariant lies in the kernel of \(L\), by Proposition 59.12.

Two pieces of kinematics are used throughout. Writing \(\vect{V}=\tfrac{1}{2}(\vect{v}+\vect{v}_{1})\) for the centre-of-mass velocity and \(\vect{g}=\vect{v}-\vect{v}_{1}\) for the relative velocity,

\begin{equation}\tag{A69.4} \vect{v}=\vect{V}+\tfrac{1}{2}\vect{g}\ec\qquad \vect{v}_{1}=\vect{V}-\tfrac{1}{2}\vect{g}\ec \end{equation}

and an elastic collision leaves \(\vect{V}\) and \(g=\abs{\vect{g}}\) unchanged while turning \(\vect{g}\) into \(\vect{g}'=g\,\vect{n}\) for a unit vector \(\vect{n}\) at angle \(\chi\) to \(\hat{\vect{g}}\). Second, averaging over the azimuth of \(\vect{n}\) about \(\hat{\vect{g}}\) at fixed \(\chi\),

\begin{equation}\tag{A69.5} \avg{\vect{n}\otimes\vect{n}}_{\varepsilon} =\tfrac{1}{3}\identity +P_{2}(\cos\chi)\bigl(\hat{\vect{g}}\otimes \hat{\vect{g}}\bigr)^{0}\ec \end{equation}

with \(P_{2}(x)=\tfrac{1}{2}(3x^{2}-1)\) and a superscript nought denoting the traceless part. Equation (A69.5) says that the only tensors surviving the azimuthal average are \(\identity\) and \(\hat{\vect{g}}\otimes\hat{\vect{g}}\), the coefficient of the second being fixed by contracting both sides with \(\hat{\vect{g}}\otimes\hat{\vect{g}}\).

Define, once and for all,

\begin{equation}\tag{A69.6} A_{2}=\int\!\dd\Omega\;F(\chi)\left[1-P_{2}(\cos\chi)\right]\ec \end{equation}

which by Equation (A69.2) is the same for every collision because \(F\) does not depend on \(g\). It is positive, since \(P_{2}\leq1\) with equality only in the forward direction.

The shear eigenfunction

Lemma A69.1 (The traceless quadratic is an eigenfunction).

For Maxwell molecules,

\begin{equation}\tag{A69.7} L\bigl[(\vect{C}\otimes\vect{C})^{0}\bigr] =-\nu_{\eta}\,(\vect{C}\otimes\vect{C})^{0}\ec\qquad \nu_{\eta}=\tfrac{1}{2}A_{2}n\ep \end{equation}

Rests on Equations (A69.2), (A69.5) and (A69.6).

Proof.

Derives Lemma A69.1. Take \(\varphi=(\vect{C}\otimes\vect{C})^{0}\) and write \(\vect{W}=\vect{V}-\vect{u}\). By Equation (A69.4),

\begin{equation*} \vect{C}\otimes\vect{C}+\vect{C}_{1}\otimes\vect{C}_{1} =2\,\vect{W}\otimes\vect{W} +\tfrac{1}{2}\,\vect{g}\otimes\vect{g}\ec \end{equation*}

the cross terms cancelling between the two molecules. Since \(\vect{W}\) and \(g\) survive the collision, only the last term changes, and

\begin{equation}\tag{A69.8} \varphi'+\varphi_{1}'-\varphi-\varphi_{1} =\tfrac{1}{2}\left[(\vect{g}'\otimes\vect{g}')^{0} -(\vect{g}\otimes\vect{g})^{0}\right]\ep \end{equation}

Integrate over the scattering. With \(\vect{g}'=g\vect{n}\) and Equation (A69.5), the traceless part of \(\vect{g}'\otimes\vect{g}'\) averages over the azimuth to \(g^{2}P_{2}(\cos\chi)(\hat{\vect{g}}\otimes\hat{\vect{g}})^{0} =P_{2}(\cos\chi)(\vect{g}\otimes\vect{g})^{0}\), so by Equation (A69.6)

\begin{equation}\tag{A69.9} \int\!\dd\Omega\;F(\chi)\, \left(\varphi'+\varphi_{1}'-\varphi-\varphi_{1}\right) =-\tfrac{1}{2}A_{2}\,(\vect{g}\otimes\vect{g})^{0}\ep \end{equation}

It remains to average over \(\vect{v}_{1}\). With \(\avg{\vect{C}_{1}}=\vect{0}\) and \(\avg{\vect{C}_{1}\otimes\vect{C}_{1}} =(k_{\mathrm{B}}T/m)\identity\),

\begin{equation*} \avg{\vect{g}\otimes\vect{g}} =\vect{C}\otimes\vect{C} +\frac{k_{\mathrm{B}}T}{m}\identity\ec \end{equation*}

whose traceless part is \((\vect{C}\otimes\vect{C})^{0}\): the isotropic term contributes nothing. Multiplying by the number density \(n\) from \(\int f_{0}\,\dd^{3}v_{1}=n\) gives Equation (A69.7).

The heat-flux eigenfunction

Lemma A69.2 (The Sonine vector is an eigenfunction).

For Maxwell molecules,

\begin{equation}\tag{A69.10} L\left[\left(C^{2} -\frac{5k_{\mathrm{B}}T}{m}\right)\vect{C}\right] =-\nu_{\kappa}\left(C^{2} -\frac{5k_{\mathrm{B}}T}{m}\right)\vect{C}\ec\qquad \nu_{\kappa}=\tfrac{1}{3}A_{2}n\ep \end{equation}

Rests on Equation (A69.2), Equation (A69.5), Equation (A69.6) and Proposition 59.12.

Proof.

Derives Lemma A69.2. The term linear in \(\vect{C}\) is a collision invariant and lies in the kernel of \(L\) by Proposition 59.12, so it may be dropped: it is enough to compute \(L[C^{2}\vect{C}]\).

With \(\vect{W}=\vect{V}-\vect{u}\), expanding \(C^{2}=W^{2}+\vect{W}\cdot\vect{g}+\tfrac{1}{4}g^{2}\) and \(C_{1}^{2}=W^{2}-\vect{W}\cdot\vect{g}+\tfrac{1}{4}g^{2}\) and adding,

\begin{equation*} C^{2}\vect{C}+C_{1}^{2}\vect{C}_{1} =2\vect{W}\left(W^{2}+\tfrac{1}{4}g^{2}\right) +\vect{g}\,(\vect{W}\cdot\vect{g})\ep \end{equation*}

Only the last term changes in the collision, so

\begin{equation}\tag{A69.11} \varphi+\varphi_{1}-\varphi'-\varphi_{1}' =\vect{g}\,(\vect{W}\cdot\vect{g}) -\vect{g}'(\vect{W}\cdot\vect{g}')\ep \end{equation}

Averaging the second term over the azimuth with \(\vect{g}'=g\vect{n}\) and Equation (A69.5),

\begin{equation*} \avg{\vect{g}'(\vect{W}\cdot\vect{g}')}_{\varepsilon} =g^{2}\left[\tfrac{1}{3}\vect{W} +P_{2}(\cos\chi)\bigl(\hat{\vect{g}}\otimes \hat{\vect{g}}\bigr)^{0}\!\cdot\vect{W}\right]\ec \end{equation*}

while the first term is the same expression with \(P_{2}\) replaced by \(1\). Their difference, integrated against \(F\), leaves only the traceless part, the isotropic pieces cancelling:

\begin{equation}\tag{A69.12} \int\!\dd\Omega\;F(\chi) \left(\varphi'+\varphi_{1}'-\varphi-\varphi_{1}\right) =-A_{2}\,(\vect{g}\otimes\vect{g})^{0}\!\cdot\vect{W}\ep \end{equation}

Now average over \(\vect{v}_{1}\). Writing \(a=k_{\mathrm{B}}T/m\) and using \(\avg{\vect{C}_{1}}=\vect{0}\), \(\avg{C_{1,i}C_{1,j}}=a\delta_{ij}\), \(\avg{\vect{C}_{1}C_{1}^{2}}=\vect{0}\) and \(\avg{\vect{C}_{1}(\vect{C}\cdot\vect{C}_{1})}=a\vect{C}\),

\begin{align} \avg{\vect{g}\,(\vect{g}\cdot\vect{W})} &=\tfrac{1}{2}\avg{(\vect{C}-\vect{C}_{1}) \left(C^{2}-C_{1}^{2}\right)} =\tfrac{1}{2}\left(C^{2}\vect{C}-3a\vect{C}\right)\ec \tag{A69.13}\\ \tfrac{1}{3}\avg{g^{2}\vect{W}} &=\tfrac{1}{6}\avg{\abs{\vect{C}-\vect{C}_{1}}^{2} (\vect{C}+\vect{C}_{1})} =\tfrac{1}{6}\left(C^{2}\vect{C}+a\vect{C}\right)\ep \tag{A69.14} \end{align}

Subtracting Equation (A69.14) from Equation (A69.13),

\begin{equation}\tag{A69.15} \avg{(\vect{g}\otimes\vect{g})^{0}\!\cdot\vect{W}} =\tfrac{1}{2}\left(C^{2}-3a\right)\vect{C} -\tfrac{1}{6}\left(C^{2}+a\right)\vect{C} =\tfrac{1}{3}\left(C^{2}-5a\right)\vect{C}\ec \end{equation}

and multiplying by \(A_{2}n\) gives Equation (A69.10).

Remark A69.3 (What makes the ratio a pure number).

Lemmas A69.1 and A69.2 carry the same \(A_{2}\) and differ only in the rational factor \(\tfrac{1}{2}\) against \(\tfrac{1}{3}\), which came from Equation (A69.8) and Equation (A69.15) respectively — that is, from the kinematics and the Gaussian moments alone. Hence

\begin{equation}\tag{A69.16} \frac{\nu_{\kappa}}{\nu_{\eta}}=\frac{2}{3}\ec \end{equation}

whatever the function \(F(\chi)\), and with it whatever the strength \(K\) of the interaction. The scattering law fixes the size of both transport coefficients and has no say in their ratio.

The first-order solution

Lemma A69.4 (The driving term).

To first order in the gradients, with the time derivatives eliminated by the Euler equations,

\begin{equation}\tag{A69.17} \frac{1}{f_{0}}\left(\pp_{t}+\vect{v}\cdot\nabla\right)f_{0} =\left(\mathcal{S}-\tfrac{5}{2}\right) \vect{C}\cdot\nabla\ln T +\frac{m}{k_{\mathrm{B}}T} (\vect{C}\otimes\vect{C})^{0}\!:\!(\nabla\vect{u})^{0}\ec \end{equation}

where \(\mathcal{S}=mC^{2}/(2k_{\mathrm{B}}T)\). Rests on Equation (59.14) and Phenomenon 58.3.

Proof.

Derives Lemma A69.4. From \(\ln f_{0}=\ln n+\tfrac{3}{2}\ln\!\left(m/(2\pi k_{\mathrm{B}}T)\right) -mC^{2}/(2k_{\mathrm{B}}T)\) and \(\dd\vect{C}=-\dd\vect{u}\),

\begin{equation*} \frac{\mathrm{D}f_{0}}{f_{0}} =\mathrm{D}\ln n +\left(\mathcal{S}-\tfrac{3}{2}\right)\mathrm{D}\ln T +\frac{m}{k_{\mathrm{B}}T}\vect{C}\cdot\mathrm{D}\vect{u}\ec \end{equation*}

with \(\mathrm{D}=\pp_{t}+\vect{v}\cdot\nabla =\mathrm{D}_{t}+\vect{C}\cdot\nabla\) and \(\mathrm{D}_{t}\) the material derivative. At zeroth order the Euler equations give \(\mathrm{D}_{t}\ln n=-\nabla\cdot\vect{u}\), \(\mathrm{D}_{t}\ln T=-\tfrac{2}{3}\nabla\cdot\vect{u}\) and \(\mathrm{D}_{t}\vect{u}=-(k_{\mathrm{B}}T/m)\nabla\ln(nT)\), the last from \(\rho\,\mathrm{D}_{t}\vect{u}=-\nabla p\) with \(p=nk_{\mathrm{B}}T\) by Phenomenon 58.3. Substituting, the terms in \(\vect{C}\cdot\nabla\ln n\) cancel identically between the first and third contributions; those in \(\vect{C}\cdot\nabla\ln T\) combine to \(\left(\mathcal{S}-\tfrac{5}{2}\right)\vect{C}\cdot\nabla\ln T\); and the terms in \(\nabla\cdot\vect{u}\) combine with \((m/k_{\mathrm{B}}T)\,\vect{C}\otimes\vect{C}:\nabla\vect{u}\) to leave its traceless part. That is Equation (A69.17).

Because the two driving terms of Equation (A69.17) are precisely the two eigenfunctions of Lemmas A69.1 and A69.2, the equation \(L[\Phi]=\mathrm{D}f_{0}/f_{0}\) is solved by division rather than by expansion:

\begin{equation}\tag{A69.18} \Phi=-\frac{1}{\nu_{\kappa}} \left(\mathcal{S}-\tfrac{5}{2}\right)\vect{C}\cdot\nabla\ln T -\frac{1}{\nu_{\eta}}\frac{m}{k_{\mathrm{B}}T} (\vect{C}\otimes\vect{C})^{0}\!:\!(\nabla\vect{u})^{0}\ep \end{equation}

For Maxwell molecules the first Chapman–Enskog approximation is therefore not an approximation at all: Equation (A69.18) is the exact first-order solution, and no expansion in Sonine polynomials is needed. This is the sense in which \(\nu=5\) is the one interaction for which the transport problem closes.

The transport coefficients

Proposition A69.5 (Viscosity and conductivity).

With \(\Phi\) given by Equation (A69.18),

\begin{equation}\tag{A69.19} \eta=\frac{nk_{\mathrm{B}}T}{\nu_{\eta}}\ec\qquad \kappa=\frac{5nk_{\mathrm{B}}^{2}T}{2m\nu_{\kappa}}\ep \end{equation}

Rests on Equation (A69.18) and Proposition 15.29.

Proof.

Derives Proposition A69.5. Viscosity. The deviatoric pressure tensor is \(\int m(\vect{C}\otimes\vect{C})^{0}f_{0}\Phi\,\dd^{3}v\), and only the second term of Equation (A69.18) contributes, the first being odd in \(\vect{C}\). Using the isotropic fourth moment

\begin{equation}\tag{A69.20} \int (\vect{C}\otimes\vect{C})^{0}_{ij} (\vect{C}\otimes\vect{C})^{0}_{kl}f_{0}\,\dd^{3}v =n\left(\frac{k_{\mathrm{B}}T}{m}\right)^{2} \left(\delta_{ik}\delta_{jl}+\delta_{il}\delta_{jk} -\tfrac{2}{3}\delta_{ij}\delta_{kl}\right)\ec \end{equation}

whose contraction with a symmetric traceless tensor returns twice that tensor, the deviatoric stress is \(-(2nk_{\mathrm{B}}T/\nu_{\eta})(\nabla\vect{u})^{0}\). Comparison with \(-2\eta(\nabla\vect{u})^{0}\) gives the first of Equation (A69.19).

Conductivity. The heat flux is \(\vect{q}=\int\tfrac{1}{2}mC^{2}\vect{C}f_{0}\Phi\,\dd^{3}v\), and now only the first term of Equation (A69.18) contributes. Using \(\int C_{i}C_{j}h(C)f_{0}\,\dd^{3}v =\tfrac{1}{3}\delta_{ij}\int C^{2}h(C)f_{0}\,\dd^{3}v\) and \(C^{2}=(2k_{\mathrm{B}}T/m)\mathcal{S}\),

\begin{equation*} \vect{q}=-\frac{k_{\mathrm{B}}T}{\nu_{\kappa}}\cdot \frac{1}{3}\cdot\frac{2k_{\mathrm{B}}T}{m}\cdot n\, \avg{\mathcal{S}^{2}\left(\mathcal{S} -\tfrac{5}{2}\right)}\,\nabla\ln T\ep \end{equation*}

For a Maxwellian \(\mathcal{S}\) has the density proportional to \(\sqrt{\mathcal{S}}\,\ee^{-\mathcal{S}}\), so \(\avg{\mathcal{S}}=\tfrac{3}{2}\), \(\avg{\mathcal{S}^{2}}=\tfrac{15}{4}\) and \(\avg{\mathcal{S}^{3}}=\tfrac{105}{8}\) by Proposition 15.29, whence \(\avg{\mathcal{S}^{3}}-\tfrac{5}{2}\avg{\mathcal{S}^{2}} =\tfrac{105}{8}-\tfrac{75}{8}=\tfrac{15}{4}\). Substituting and using \(\nabla\ln T=\nabla T/T\) gives \(\vect{q}=-\left(5nk_{\mathrm{B}}^{2}T/(2m\nu_{\kappa})\right)\nabla T\), which against Fourier's law Equation (59.37) is the second of Equation (A69.19).

Corollary A69.6 (The Eucken ratio and the Prandtl number).

For a monatomic Maxwell gas,

\begin{equation}\tag{A69.21} \frac{\kappa}{\eta}=\frac{15k_{\mathrm{B}}}{4m}\ec\qquad \frac{\kappa m}{\eta\,c_{V}^{(1)}}=\frac{5}{2}\ec\qquad \mathrm{Pr}=\frac{c_{p}\eta}{\kappa}=\frac{2}{3}\ep \end{equation}

Rests on Equations (A69.16) and (A69.19).

Proof.

Derives Corollary A69.6. Dividing the two members of Equation (A69.19) and using Equation (A69.16),

\begin{equation*} \frac{\kappa}{\eta} =\frac{5k_{\mathrm{B}}}{2m}\cdot\frac{\nu_{\eta}}{\nu_{\kappa}} =\frac{5k_{\mathrm{B}}}{2m}\cdot\frac{3}{2} =\frac{15k_{\mathrm{B}}}{4m}\ep \end{equation*}

With \(c_{V}^{(1)}=\tfrac{3}{2}k_{\mathrm{B}}\) the second member is \((15/4)/(3/2)=5/2\), and with \(c_{p}=\tfrac{5}{2}k_{\mathrm{B}}/m\) per unit mass the third is \((5/2)/(15/4)=2/3\).

Remark A69.7 (Back to the chapter, and the size of what is left out).

The Chapman–Enskog Transport Coefficients for Maxwell Molecules discharges the derivation owed at Proposition 61.9 of Experiment: The Kinetic Theory Tested, and completes the ladder that chapter sets out. The mean-free-path estimate Equation (61.10) gives \(1\) for the Eucken ratio because it assigns one free path to every transport; the relaxation-time model Proposition 61.7 gives \(\tfrac{5}{3}\) because it assigns one relaxation time to every molecule; and Equation (A69.21) gives \(\tfrac{5}{2}\) because Lemmas A69.1 and A69.2 relax the two fluxes at different rates, in the ratio \(3:2\). The physical content of the factor is exactly that difference: momentum and energy are not carried by the same molecules.

Two limits are worth stating precisely, and both are settled here rather than referred elsewhere.

The interaction. The result is exact for \(\nu=5\) and for no other force law, and the place it fails is identifiable to the line. Everything above rests on Equation (A69.2), that \(g\sigma\) does not depend on the relative speed; that is what let the angular integral in Equation (A69.9) and Equation (A69.12) be taken with \(F\) pulled out, and it is what made the two polynomials eigenfunctions. For \(\nu\neq5\) the product \(g\sigma\) carries a factor \(g^{1-4/(\nu-1)}\), the angular integral no longer factors from the average over \(\vect{v}_{1}\), and neither polynomial is an eigenfunction any more. The first-order equation then has to be solved by expansion instead of by division, \(\tfrac{5}{2}\) becomes the leading term of a series rather than the answer, and the corrections are of the order of the spread in \(g^{1-4/(\nu-1)}\) over a Maxwellian — small, because that exponent is small for every interaction between the two limits \(\nu=5\) and \(\nu\to\infty\), but not zero.

The internal modes. Equation (A69.21) is a monatomic result, while the measured value quoted in Equation (61.8) is \(1.97\) for air, which is diatomic. The correction is short enough to carry out. Suppose the translational energy is transported as it is above, with its coefficient \(\tfrac{5}{2}\), while the internal energy — which the collisions do not redistribute in the same way, being carried by the molecule rather than by its motion — is transported at the rate at which a molecule itself diffuses, whose coefficient is \(\rho D/\eta=1\). Splitting the heat capacity into the two parts, \(c_{V}=c_{\mathrm{tr}}+c_{\mathrm{int}}\) with \(mc_{\mathrm{tr}}=\tfrac{3}{2}k_{\mathrm{B}}\),

\begin{equation}\tag{A69.22} \frac{\kappa m}{\eta} =\frac{5}{2}\left(mc_{\mathrm{tr}}\right) +1\cdot\left(mc_{\mathrm{int}}\right) =\frac{15}{4}k_{\mathrm{B}}+\left(mc_{V} -\frac{3}{2}k_{\mathrm{B}}\right)\ep \end{equation}

Dividing by \(mc_{V}\) and writing \(mc_{V}=k_{\mathrm{B}}/(\gamma-1)\),

\begin{equation}\tag{A69.23} \frac{\kappa m}{\eta\,c_{V}} =(\gamma-1)\left(\frac{15}{4}-\frac{3}{2}\right)+1 =\frac{9\gamma-5}{4}\ec \end{equation}

which is Eucken's relation [Eucken:1913] and which returns \(\tfrac{5}{2}\) at \(\gamma=\tfrac{5}{3}\), as it must. For a diatomic gas at ordinary temperature \(\gamma=\tfrac{7}{5}\) (Equation (59.13)), and Equation (A69.23) gives

\begin{equation}\tag{A69.24} \frac{\kappa m}{\eta\,c_{V}} =\frac{9\times\tfrac{7}{5}-5}{4}=\frac{19}{10}=1.90\ec \end{equation}

against the measured \(1.97\) of Equation (61.8): agreement to \(3.6\,\mathrm{\%}\), from a calculation whose only new ingredient is the assumption that internal energy rides with the molecule. That residual is the crudeness of the assumption — collisions do exchange internal and translational energy, at a rate this argument sets to zero — and it is the honest size of what separates Equation (A69.21) from a measurement on a real gas.