logo

기체운동론으로 유도되는 이상기체방정식 📂열물리학

기체운동론으로 유도되는 이상기체방정식

정의1

면적 MM에 작용하는 유체의 압력 pp를 유체가 면적 MM에 수직으로 작용하는 힘 FF와 면적 MM의 비율로 정의한다.

p:=FM[N/m2] p:=\frac{F}{M} \left[ \mathrm{N/m^{2}} \right]

공식

기체의 부피를 VV, 온도TT, 분자 수를 NN이라고 하자. 그러면 기체의 압력 pp는 다음과 같은 식을 만족한다.

pV=NkBT p V = N k_{B}T

이때 kB=1.3807×1023J/Kk_{B} = 1.3807 \times 10^{-23} J / K볼츠만 상수Boltzmann constant이다.

설명

위 식은 이상기체 방정식라고 불린다. 처음에는 실험법칙으로부터 유도되었지만 후에 아래와 같이 수식적으로 유도되었다. 이는 기체의 압력을 수식적으로 정리하는 과정으로부터 나온 결과이다.

유도 과정을 이해하기 위해서는 확률과 전체에 대한 비율이 같은 개념이라는 점을 유념해야한다. 가령 두 주사위를 던져서 나온 눈의 합이 7인 경우는 다음의 표에서 알 수 있듯이 총 3636가지의 경우 중에서 66가지이다. 이 비율은 곧 두 주사위를 던져서 나온 눈의 합이 77일 확률 736\dfrac{7}{36}이 된다.

눈의 합23456789101112
경우의 수1234565432136

유도

단위부피당 분자의 수를 n=N/Vn = N / V이라고 표기하자. 기체 분자의 속력 분포의 확률밀도함수는 다음과 같이 주어진다.

맥스웰-볼츠만 분포

f(v)=4π(m2kBT)3/2v2emv2/2kBT f(v) = \dfrac{4}{\sqrt{ \pi}} \left( \dfrac{m}{2 k_{B} T} \right)^{3/2} v^{2} e^{-mv^2 / 2k_{B}T }

f(v)f(v)는 기체 분자의 속력이 vv일 확률, 즉 전체 분자의 수 중에서 속력이 vv인 분자의 비율을 의미한다. 따라서 다음의 식은 단위 부피당 속력이 vv인 분자의 수 를 의미한다.

nf(v) nf (v)

이제 아래와 같이 특정한 방향을 고정했다고 하자. 이를 zz축이라 하면 θ\theta는 구면좌표계의 변수 (r,θ,ϕ)(r,\theta, \phi)θ\theta와 같다.

1.jpg

입체각Ω\Omega인 영역 안에 기체 분자들이 있을 비율은 다음과 같다.

Ω4π=12sinθ \dfrac{\Omega}{4\pi} = \dfrac{1}{2}\sin \theta

따라서 다음의 식은 단위 부피당, 입체각 0Ω0\sim\Omega 사이에서 운동하고 있는, 속력이 vv인 분자들의 수 를 의미한다.

(nf(v))(12sinθ) \left( n f(v) \right) \left( \dfrac{1}{2} \sin\theta\right)

이제 다음 그림과 같이 θ\theta의 각도로 면적이 AA인 벽을 때리는 속력이 vv인 분자를 고려하자.

3.pgn

속력이 vv인 분자가 tt초 후에 벽을 때렸다고 하자. 그러면 tt시간 분자들이 휩쓸고간 영역(진한 부분)을 분자의 수라고 둘 수 있다. 평행사변형의 넓이는 밑변×\times높이이므로, 시간 tt동안 면적 AA를 때리는 분자의 수는 다음과 같다.

Avtcosθ Avt \cos\theta

따라서 단위면적을 단위시간동안 θ\theta의 각도로 때리는 속력이 vv인 분자들의 수 는 다음과 같다.

(nf(v))(12sinθ)(vcosθ) \left( n f(v) \right) \left( \dfrac{1}{2} \sin\theta\right) \left( v \cos \theta \right)

힘은 운동량의 변화량이므로 벽과 충돌한 분자들의 운동량 변화량을 구해보자. 운동량 변화는 벽에 수직한 방향으로만 일어난다. 벽과 가까워 지는 쪽을 ++라고 두면 다음과 같다.

mvcosθ(mvcosθ)=2mvcosθ mv\cos\theta - (-mv\cos\theta) = 2mv \cos\theta

따라서 다음의 식은 단위시간당, θ\theta방향으로 운동하는 속력이 vv인 분자들이 수직한 벽에 가하는 힘 이 된다.

(nf(v))(12sinθ)(vcosθ)(2mvcosθ) \left( n f(v) \right) \left( \dfrac{1}{2} \sin\theta\right) \left( v \cos \theta \right) \left( 2mv \cos\theta \right)

따라서 이를 모든 속력 v=0v = 0 \sim \infty, 모든 θ=0π/2\theta = 0 \sim \pi/2에 대해서 적분하면 기체가 면적에 가하는 압력이 된다.

p=v=0θ=0π/2(nf(v))(12sinθ)(vcosθ)(2mvcosθ)dvdθ=nmv=0f(v)v2dvθ=0π/2cos2θsinθdθ \begin{align*} p &= \int_{v = 0} ^{\infty} \int_{\theta=0}^{\pi/2} \left( n f(v) \right) \left( \dfrac{1}{2} \sin\theta\right) \left( v \cos \theta \right) \left( 2mv \cos\theta \right) dv d\theta \\ &= nm \int_{v = 0} ^{\infty} f(v) v^{2} dv \int_{\theta=0}^{\pi/2} \cos^{2} \theta \sin\theta d\theta \end{align*}

vv에 대한 적분은 v2v^{2}기댓값이다. v=0f(v)v2dv=v2=3kBTm \int_{v = 0} ^{\infty} f(v) v^{2} dv = \left\langle v^{2} \right\rangle = \dfrac{3 k_{B} T}{m}

θ\theta에 대한 적분은 cosθ=x\cos \theta = x로 치환하면 sinθdθ=dx-\sin\theta d\theta = dx이므로 다음과 같다.

10x2dx=13x310=13 \int_{1}^{0} - x^{2} dx = \left. -\dfrac{1}{3}x^{3} \right|_{1}^{0} = \dfrac{1}{3}

따라서 압력 pp는 다음과 같다.

p=nm3kBTm13=nkBT p = n m \dfrac{3 k_{B} T}{m} \dfrac{1}{3} = n k_{B} T

이때 n=N/Vn=N/V은 단위부피당 분자 수였으므로 다음의 결과를 얻는다.

p=NkBTV    pV=NkBT p=\dfrac{Nk_{B} T}{V} \implies pV = Nk_{B}T


  1. Stephen J. Blundell and Katherine M. Blundell, 열 물리학(Concepts in Thermal Physics, 이재우 역) (2nd Edition, 2014), p71-75 ↩︎