기체운동론

물리 유체역학 계산물리 마지막 수정: 2026-09-02 04:51:33

1. 개요[편집]

기체운동론 Kinetic Theory of Gases
대상	분자 운동의 통계로 기체의 거시 성질을 유도
핵심 가정	분자 카오스 (Stosszahlansatz)
지배 방정식	볼츠만 방정식
평형해	맥스웰-볼츠만 분포
연속체로 가는 다리	채프먼-엔스코그 전개
유효 범위 판정	크누센수
대표 예측	기체 점성이 압력과 무관하다

분자 $10^{23}$ 개를 세는 대신 확률분포 하나를 세기로 한 순간, 열역학이 역학의 따름정리가 됐다.

기체운동론(kinetic theory of gases)은 기체를 끊임없이 날아다니며 충돌하는 분자들의 집합으로 보고, 그 운동의 통계로부터 압력·온도·점성·열전도 같은 거시 물성을 유도하는 이론 체계다. 열이 무엇인지에 대한 답이자, 나비에-스토크스 방정식이 왜 성립하고 언제 깨지는지에 대한 유일한 원리적 설명이기도 하다.

심위키 관점에서 이 문서는 허브다. 개별 주제는 이미 각자의 자리가 있다.

묻고 싶은 것	가는 곳
지배방정식의 형태와 충돌적분	볼츠만 방정식
평형 속도분포와 대표 속력들	맥스웰-볼츠만 분포
충돌 사이 거리와 그 통계	평균자유행로
연속체 근사가 언제 깨지나	크누센수
볼츠만에서 나비에-스토크스 유도	채프먼-엔스코그 전개
희박기체를 실제로 푸는 코드	직접 시뮬레이션 몬테카를로

그래서 여기서는 계보와 논리 구조 — 어떤 가정에서 출발해 어떤 순서로 무엇이 유도되고, 각 단계에서 무엇을 버렸으며, 그 결과 계산 도구의 사다리가 어떻게 배열되는지 — 를 다룬다.

2. 계보[편집]

이 분야는 두 번 잊혔다가 세 번째에 자리를 잡았다.

1738 — 다니엘 베르누이가 Hydrodynamica 에서 압력을 “입자가 벽을 두들기는 것”으로 설명한다. 아무도 진지하게 받지 않았다. 열은 그때까지 칼로릭이라는 유체였다.
1820~45 — 헤라파스와 워터스턴이 각각 상당히 완성된 형태를 내놓지만 왕립학회에서 게재를 거부당한다. 워터스턴의 원고는 50년 뒤 레일리가 서랍에서 찾아내 발표한다.¹
1857 — 클라우지우스가 「우리가 열이라 부르는 운동의 종류에 관하여」를 내고, 이듬해 평균자유행로 개념을 도입한다. 여기서부터 분야가 굴러가기 시작한다.
1860 — 맥스웰이 대칭성 논증만으로 속도분포를 유도한다. 물리학에 통계적 법칙이 등장한 최초의 순간으로 꼽힌다. 1866년에는 수송이론으로 점성을 계산하고, 그 결과가 반직관적이라 직접 실험까지 한다.
1865 — 로슈미트가 기체 물성에서 분자의 크기와 수밀도를 처음으로 추정한다.
1872 — 볼츠만이 분포함수의 시간 발전 방정식과 H-정리를 세운다. 열역학 제2법칙이 역학에서 유도되는 순간이자, 로슈미트의 가역성 반론과 체르멜로의 재귀 반론이 시작되는 지점이기도 하다.
1909 — 크누센이 희박기체 유동을 체계적으로 측정하고, 미세관 유량의 최솟값(크누센 최소)을 보고한다.
1916~17 — 채프먼과 엔스코그가 독립적으로 볼츠만 방정식에서 수송계수를 계산하는 체계적 전개를 완성한다. 분자에서 나비에-스토크스까지의 다리가 놓인다.
1963 — 버드가 직접 시뮬레이션 몬테카를로를 제안한다. 이론이 계산 도구가 된다.

원자의 실재를 두고 마흐와 오스트발트가 맹렬히 반대하던 시대에 볼츠만이 겪은 고립은 잘 알려져 있다. 원자론이 결정적으로 승리한 것은 1905년 아인슈타인의 브라운 운동 분석과 1908년 페랭의 실험이었고, 볼츠만은 그것을 보지 못하고 1906년에 죽었다.

3. 압력과 온도의 미시적 정의[편집]

운동론이 처음에 판 것은 사실 딱 하나, 압력이 무엇인가였다.

벽에 수직으로 부딪혔다 튀어나오는 분자 하나가 넘겨주는 운동량은 $2mv_x$ 다. 단위 시간에 넓이 $A$ 에 부딪히는 $v_x$ 분자의 수를 세어 전부 더하면, 등방성 $\langle v_x^2\rangle = \tfrac13\langle v^2\rangle$ 을 써서

p = n m \langle v_x^2\rangle = \tfrac13\, n\, m\, \langle v^2\rangle = \tfrac23\, n\, \langle \varepsilon_{\rm tr}\rangle

가 나온다. 압력은 운동에너지 밀도의 2/3이라는 것이다. 여기에 이상기체 법칙 $p = nk_BT$ 를 비교하면 온도의 정체가 곧바로 드러난다.

\langle \varepsilon_{\rm tr}\rangle = \tfrac32 k_B T

온도는 병진 운동에너지의 평균에 비례하는 양이지 별도의 실체가 아니다. 이차 자유도 하나당 $\tfrac12k_BT$ 라는 등분배 정리가 여기서 나오고, 분자동역학 코드가 온도를 정의하는 방식이 정확히 이 식이다. 뒤집어 말하면 유한한 입자 수로 잰 “순간 온도”는 표본 분산이라 요동을 달고 다닌다.

같은 세는 방식이 하나 더 예쁜 결과를 준다. 넓이 $A$ 인 작은 구멍을 단위 시간에 빠져나가는 분자 수는

\Phi = \tfrac14\, n\, \bar v\, A

이고, $\bar v \propto 1/\sqrt m$ 이므로 유출 속도는 분자량의 제곱근에 반비례한다(그레이엄의 법칙). 우라늄 농축의 기체 확산법이 이 $\sqrt{352/349}$ 짜리 미세한 차이를 수천 단 반복해 짜낸 것이고, 진공 장비의 컨덕턴스 계산과 직접 시뮬레이션 몬테카를로의 유입 경계 조건도 같은 식에서 나온다. 여기서 유명한 함정 하나 — 구멍을 통과하는 분자는 빠른 놈일수록 자주 지나가므로, 빔의 속도분포는 챔버 안의 분포가 아니라 플럭스 가중된 분포다. 자세한 것은 맥스웰-볼츠만 분포에 있다.

4. 분자 카오스 — 비가역성이 몰래 들어오는 자리[편집]

지금까지는 평형만 다뤘다. 비평형으로 가려면 분포함수 $f(\mathbf x,\mathbf v,t)$ 의 시간 발전이 필요한데, 여기가 이 분야의 가장 미묘한 지점이다.

정직하게 시작하면 $N$ 입자 리우빌 방정식에서 좌표를 하나씩 적분해 내는 BBGKY 계층이 나온다. 문제는 1체 분포의 방정식에 2체 분포가, 2체의 방정식에 3체가 등장해 결코 닫히지 않는다는 것이다. 볼츠만은 첫 단계에서 이 사슬을 자른다.

f_2(\mathbf v, \mathbf v_*, t) \;\approx\; f(\mathbf v,t)\, f(\mathbf v_*,t)

이것이 분자 카오스 가정(Stosszahlansatz)이며, 충돌 직전 두 분자의 속도가 서로 무상관이라는 요구다.

핵심은 “직전”이다. 충돌 후에는 두 분자가 명백히 상관되어 있고, 그 상관을 가정에서 배제하지 않았다. 시간 방향에 대해 비대칭인 가정을 넣었으므로, 가역적인 뉴턴 역학에서 출발했는데도 결과 방정식이 비가역이 된다. H-정리의 시간의 화살은 여기서 들어온 것이지 어디서 나타난 것이 아니다. 로슈미트의 가역성 반론과 체르멜로의 재귀 반론이 정확히 이 한 줄을 겨눴고, 볼츠만의 답은 “예외적 초기조건이 아닌 한 압도적 다수의 미시상태가 그렇게 행동한다”는 통계적 진술이었다.

이 근사가 언제 정당한지는 수학적으로 부분적인 답만 있다. 볼츠만-그라드 극한( $N\to\infty$ , 분자 지름 $d\to0$ , $Nd^2$ 고정)에서 랜퍼드(1975)가 엄밀 유도에 성공했지만 아주 짧은 시간 동안만이고, 그 시간 제약은 지금도 본질적으로 풀리지 않았다. 방정식의 형태·충돌적분·H-정리의 상세는 볼츠만 방정식이 다룬다.

5. 수송계수 — 어림과 정확값[편집]

운동론이 실무에 팔아먹는 가장 큰 상품이 수송계수를 분자 파라미터에서 계산해 준다는 것이다.

어림 논증은 세 줄이면 된다. 속도 기울기가 있는 기체에서 분자는 $\lambda$ 만큼 떨어진 층의 운동량을 지고 와서 이곳에 내려놓는다. 층당 운동량 차이가 $m\lambda\,du/dy$ , 단위 넓이를 지나는 유량이 $\sim n\bar v$ 이므로

\mu \sim \tfrac13\,\rho\,\bar v\,\lambda, \qquad k \sim \tfrac13\,\rho\,\bar v\,\lambda\, c_v, \qquad D \sim \tfrac13\,\bar v\,\lambda

세 계수가 같은 구조를 공유한다. 여기에 $\lambda \propto 1/(n\sigma)$ 를 넣으면 점성에서 $n$ 이 완전히 약분되어 기체 점성이 압력과 무관해진다는 유명한 결론이 나온다.² 유도와 맥스웰 부부의 검증 실험은 평균자유행로에 있다.

문제는 앞의 $1/3$ 이 근거가 약하다는 것이다. 정확한 계수를 주는 것이 채프먼-엔스코그 전개이고, 강체구 기체의 1차 근사는

\mu = \frac{5}{16 d^2}\sqrt{\frac{m k_B T}{\pi}}

로, 어림값의 계수를 약 $0.499\,\rho\bar v\lambda$ 로 고쳐 준다. 어림이 자릿수를 맞히고 정확한 이론이 계수를 고치는 구도다. 여기서 따라 나오는 관계가 둘 더 있다.

오이켄 관계. 단원자 기체에서 $k = \tfrac52 c_v \mu$ 이고, 따라서 프란틀수가 $\mathrm{Pr} = c_p\mu/k = 2/3$ 로 고정된다. 다원자 기체에서는 내부 자유도가 병진보다 느리게 확산하므로 $k = \tfrac{9\gamma-5}{4}c_v\mu$ 로 보정된다. BGK 모형이 완화율 하나만 쓰는 바람에 $\mathrm{Pr}=1$ 에 묶이는 문제가 정확히 이 값과 어긋나는 것이고, ES-BGK·샤코프 모형이 그걸 고친다.
온도 의존성. 강체구는 $\mu\propto\sqrt T$ 만 낼 수 있지만 실기체는 $\mu\propto T^{\omega}$ , $\omega\approx0.7\sim0.8$ 이다. 분자가 딱딱한 공이 아니라 부드러운 척력을 갖기 때문이며, 그래서 DSMC 는 단면적을 상대속도에 의존시키는 VHS·VSS 모형을 표준으로 쓴다.

이 표가 운동론의 값어치를 요약한다 — 거시 계수를 실험에서 가져오는 대신 분자 모형에서 계산하고, 역으로 측정된 $\mu(T)$ 에서 분자간 힘의 형태를 역추정할 수 있다.

6. 크누센수 사다리와 계산 도구[편집]

운동론이 유체역학과 갈라지는 지점은 하나뿐이다. 분자가 형상을 느끼기 전에 서로 몇 번 부딪히는가. 그 비가 크누센수 $\mathrm{Kn}=\lambda/L$ 이고, 값에 따라 쓸 도구가 갈린다.

$\mathrm{Kn}$	상태	표준 도구
$\lesssim0.01$	연속체	나비에-스토크스 기반 CFD, 격자 볼츠만 방법
$0.01\sim0.1$	미끄럼	CFD + 미끄럼·온도 도약 경계조건
$0.1\sim10$	전이	DSMC, 이산 속도법, 모멘트법
$\gtrsim10$	자유분자	충돌 없는 자유 수송, 시험입자 몬테카를로

여기에 미시 쪽 끝을 하나 더 붙여야 그림이 완성된다. 분자동역학은 궤적을 다 푸는 대신 시공간 규모를 못 키우고, 운동론은 규모를 키우는 대신 분자 카오스를 지불한다. 조밀한 액체에서 분자 카오스가 성립하지 않아 운동론이 못 가는 영역이 정확히 MD 의 영역이고, 그래서 두 방법은 경쟁자가 아니라 이웃이다.

DSMC 로 넘어가는 지점을 실무 기준으로 적으면 이렇다. 전역 $\mathrm{Kn}$ 이 아니라 기울기 기반 국소값 $\mathrm{Kn}_{GL} = \lambda|\nabla Q|/Q$ 가 대략 $0.05$ 를 넘는 영역이 생기면 연속체 해가 신뢰를 잃는다. 충격파 내부, 급팽창부, 벽 근처 크누센층이 대표적이다. 반대로 저속·저Kn 영역에서 DSMC 를 억지로 돌리면 통계 잡음이 신호를 덮어 버리므로, 신호 대 잡음이 나쁜 미세유동에서는 잡음이 없는 이산 속도법이 낫다. 두 영역을 이어 붙이는 하이브리드와 통합 기체운동론 기법(UGKS)이 그 사이를 메운다.

7. 운동론이 아니면 못 보는 것들[편집]

연속체 방정식이 원리적으로 재현할 수 없는 현상들의 목록이 곧 이 분야의 존재 이유다.

속도 미끄럼과 온도 도약. $\mathrm{Kn}$ 이 커지면 벽에서 유체가 미끄러지고, 벽 온도와 접한 기체 온도가 다른 값을 갖는다. 이 보정 계수들은 전부 운동론 계산에서 나온다.
열유동 미끄럼(thermal creep). 압력 차이가 전혀 없어도 온도 기울기만으로 기체가 찬 쪽에서 더운 쪽으로 기어간다. 움직이는 부품 없이 온도차만으로 기체를 미는 크누센 펌프가 이 원리로 돌아간다. 나비에-스토크스에는 이런 항이 존재하지 않는다.
크누센 최소. 긴 미세관의 무차원 질량유량을 $\mathrm{Kn}$ 에 대해 그리면 $\mathrm{Kn}\approx1$ 부근에 최솟값이 나타난다. 희박기체 코드의 국룰 검증 케이스이며, 상세는 크누센수가 다룬다.
충격파 내부 구조. 충격파 두께는 $\lambda$ 규모라 그 안에서 국소 $\mathrm{Kn}$ 이 $O(1)$ 이다. 나비에-스토크스는 고마하수에서 두께를 실측·DSMC 대비 상당히 작게 내놓으며, 파면 내부의 이봉(bimodal) 속도분포는 애초에 표현할 수단이 없다.
열적 비평형. 급팽창이나 강한 충격 뒤에서는 병진·회전·진동 온도가 서로 다른 값으로 갈라진다. “온도”라는 단일 숫자가 의미를 잃는 상태이며, 재진입체 공력가열 예측이 어려운 근본 이유다.

8. 한계와 확장[편집]

기체운동론은 묽고, 단원자이고, 중성이고, 고전적인 기체를 위해 만들어졌다. 각 조건을 풀면 확장판이 하나씩 붙는다.

조밀해지면 분자 지름을 무시할 수 없고 충돌이 상관된다. 엔스코그가 강체구에 부피 배제와 상관인자를 넣은 확장을 만들었지만, 액체에 가까워지면 분자 카오스 자체가 무너져 분자동역학으로 넘어가는 것이 정석이다.
다원자 분자에서는 회전·진동 자유도와 병진 사이의 에너지 교환을 모형화해야 한다. 왕창-울렌벡 방정식이 원리적 틀이고, DSMC 실무에서는 보르그나케-라르센 방식으로 내부 에너지를 확률적으로 재분배한다.
하전입자에서는 쿨롱 힘이 장거리라 이체 충돌만 세는 그림이 깨진다. 충돌항을 무시하면 블라소프 방정식, 소각 산란의 누적을 확산으로 근사하면 란다우·포커-플랑크 형태의 충돌항이 된다.
양자 축퇴가 시작되면 맥스웰-볼츠만 통계 자체가 못 쓴다. 점유수가 1에 근접하는 극저온 원자나 금속 내 전자는 보스-아인슈타인·페르미-디랙 통계를 쓴 운동론(울링-울렌벡 방정식)이 필요하다.

수학적 지위도 아직 미완이다. 힐베르트의 여섯 번째 문제가 원자 모형에서 운동론과 유체역학을 엄밀히 유도하는 것이었는데, 대역 해의 존재성은 디페르나-리옹의 재규격화 해로 얻어졌지만 유일성은 여전히 열려 있다. 인류는 이 방정식의 해가 유일한지 모르는 채로 재진입체와 반도체 진공 공정을 설계하고 있다.³

9. 관련 문서[편집]

워터스턴의 1845년 원고에 붙은 심사평이 물리학사에서 가장 유명한 오심으로 남아 있다 — “이 논문은 난센스일 뿐이며 왕립학회 낭독에도 적합하지 않다.” 레일리가 1892년에 그 원고를 찾아 발표하면서 붙인 서문의 요지는, 학회지가 잘 정립된 분야의 논문에는 관대하고 새로운 방향에는 가혹하다는 것이었다. 130년이 지나도 심사평 게시판의 정서가 크게 다르지 않다는 점이 이 일화의 진정한 교훈. ↩
운동론이 예측한 것 중 당대에 가장 안 믿긴 것이 “기체 점성은 압력과 무관하다”였다. 밀도를 반으로 줄이면 점성도 반으로 줄 것 같은데 그렇지 않다. 맥스웰 본인도 못 믿어서 직접 쟀고, 그게 맞았다. 이론이 상식을 이긴 드문 사례이자, “일단 돌려” 이전에 “일단 재 봐”가 있었다는 V&V 정신의 조상 격 일화다. ↩
이 구도가 나비에-스토크스 방정식의 밀레니엄 문제와 정확히 같고, 사실 한 층 위에서 벌어지는 일이라 더 근본적이다. 위층도 아래층도 존재성·유일성을 모르는데 그 사이의 유도 관계는 짧은 시간만 증명돼 있다. 그래도 비행기는 뜬다. ↩