Free expansion of an ideal gas


The free expansion of an ideal gas is particularly simple physical system to analyze the thermodynamics of.

Mechanics

Consider an Ideal gas that's contained in a sealed, thermodynamically isolated box by insulating walls. In the real world, this can be approximated by thick walls made of a material with very high heat capacity, such as water (see the figure below).

center

There are two chambers, separated by a movable stopper. Initially, one chamber contains gas at temperature T1T_{1} and volume V1V_{1}, whereas the other is vacuum. The stopper is then removed, letting the gas expand from the first chamber to the second, filling all space. After expansion, the gas has a new temperature T2T_{2} and volume V2V_{2}. Experimentally, it was found that the temperature went unchanged, so that T1=T2T_{1}=T_{2}. Notably, this implies a few things:

  1. Since the temperature was unchanged, no heat was exchanged: ΔQ=0\Delta Q=0.
  2. Since there is no gas to push against during the expansion, no work was performed either: ΔW=0\Delta W=0.
  3. But the first law of thermodynamics therefore gives us ΔU=ΔQΔW=0\Delta U=\Delta Q-\Delta W=0. So the internal energy of the gas also did not change. Thus, we have U1(V1,T)=U2(V2,T)U_{1}(V_{1},T)=U_{2}(V_{2},T).

The last point implies that UU is independent of volume or, in other words, exclusively dependent on temperature: U=U(T)U=U(T). From the heat equations we know

CV=(UT)V=dUdTC_{V}=\left( \frac{ \partial U }{ \partial T } \right)_{V}=\frac{dU}{dT}

where the partial derivative becomes a total derivative since the dependency on VV has been lifted. If CVC_{V} is constant, we can write

U=dU=CV dT=CVT+constU=\int dU=\int C_{V}\ dT=C_{V}T+\text{const}

Note that this is an indefinite integral, not a definite one. The leftover constant vanishes by setting U=0U=0 at T=0T=0, leaving us with

U=CVT\boxed{U=C_{V}T}

Also from the heat equations, we can find the isobaric heat capacity as

CP=(HT)P=((U+PV)T)P=dUdT+(NkBT)T=CV+NkBC_{P}=\left( \frac{ \partial H }{ \partial T } \right)_{P}=\left( \frac{ \partial (U+PV) }{ \partial T } \right)_{P}=\frac{dU}{dT}+ \frac{ \partial (Nk_{B}T) }{ \partial T } =C_{V}+Nk_{B}

where NN is the number of particles in the gas and kBk_{B} is the Boltzmann constant. So the difference between constant-volume and isobaric heat capacities for an ideal gas is

CPCV=NkBC_{P}-C_{V}=Nk_{B}

Transformations

We can work out the equation governing a reversible, adiabatic transformation. Since for an adiabatic transformation dQ=0dQ=0 and work is dW=PdVdW=PdV, the first law of thermodynamics gives us dU=PdVdU=-PdV, but since dU=CVdTdU=C_{V}dT is also true, we get

CVdT+PdV=0C_{V}dT+PdV=0

We can use the ideal gas law PV=NkBTPV=Nk_{B}T to state the infinitesimal form

dT=d(PV)NkB=VdP+PdVNkBdT=\frac{d(PV)}{Nk_{B}}=\frac{VdP+PdV}{Nk_{B}}

so we get

0=CVdT+PdV=CVVdP+PdVNkB+PdV=CV(VdP+PdV)+NkBPdV=CVVdP+(CV+NkB)PdV=CVVdP+CPPdV=dPP+γdVV\begin{align} 0&=C_{V}dT+PdV \\ &=C_{V} \frac{VdP+PdV}{Nk_{B}} + PdV \\ &=C_{V}(VdP+PdV)+Nk_{B}PdV \\ &=C_{V}VdP+(C_{V}+Nk_{B})PdV \\ &=C_{V}VdP+C_{P}PdV \\ &=\frac{dP}{P}+ \gamma\frac{dV}{V} \end{align}

where γCP/CV\gamma\equiv C_{P}/C_{V}. If γ\gamma is constant, we can integrate the former equation to get

0=dPP+γdVV=lnP+γlnV+K0=\int \frac{dP}{P}+ \gamma \int \frac{dV}{V}=\ln P+\gamma \ln V+K

for some constant KK. From this we get

lnP=γlnV+KP=eγlnV+K=eKeln(Vγ)=CVγ\ln P=-\gamma \ln V+K\quad\Rightarrow \quad P=e^{-\gamma \ln V+K}=e^{K}e^{\ln (V^{-\gamma})}=CV^{-\gamma}

where CC is some other constant. Therefore

PVγ=constant\boxed{PV^{\gamma}=\text{constant}}

Using the equation of state PV=NkBTPV=Nk_{B}T, we get the equivalent form

TVγ1=constant\boxed{TV^{\gamma-1}=\text{constant}}

Since PV=constantPV=\text{constant} determines an isotherm (because PVTPV\propto T in the ideal gas law), PVγ=constantPV^{\gamma}=\text{constant} has a similar shape to an isotherm, but since γ>1\gamma>1 as CP=CV+NkB>CVC_{P}=C_{V}+Nk_{B}>C_{V}, it has a steeper slope in a PVPV graph. Thus, for an ideal gas, adiabatic transformations have a steeper PVPV slope.

center