희박기체역학

유체역학 계산물리 전산유체역학 마지막 수정: 2026-08-22 04:39:07

1. 개요[편집]

희박기체역학 Rarefied Gas Dynamics
판정 척도	크누센수 Kn = λ / L
지배 방정식	볼츠만 방정식 (6차원 위상공간 + 시간)
주력 수치기법	DSMC · 이산속도법 · 모멘트법 · UGKS
용어 정립	H. S. Tsien, Superaerodynamics (1946)
주 무대	진공 공정 · 고고도 재진입 · MEMS · 우주 추진

유체역학 첫 시간에 배운 “유체는 연속체다”라는 문장에는 유효기간이 있었다. 그 기한이 지난 뒤의 세계가 여기다.

희박기체역학(rarefied gas dynamics)은 기체 분자의 평균자유행로가 문제의 대표 길이와 견줄 만큼 커져 연속체 가정이 무너지는 영역의 기체 유동을 다루는 분야다. 판정 척도는 크누센수 $\mathrm{Kn} = \lambda/L$ 이고, 지배 방정식은 나비에-스토크스 방정식이 아니라 그 상위 이론인 볼츠만 방정식이다.

“희박”이라는 말이 오해를 부르는데, 기체가 묽어야만 여기 오는 게 아니다. 상온 1기압 공기의 $\lambda$ 는 약 68 nm이므로, 압력을 낮추지 않아도 형상을 100 nm까지 줄이면 그대로 희박 영역이다. 진공 챔버와 고고도 대기는 $\lambda$ 를 키워서, MEMS와 나노 기공은 $L$ 을 줄여서 같은 지점에 도착한다. 그래서 이 분야의 응용 목록에는 재진입 캡슐과 마이크로 가속도계가 나란히 앉아 있다.

이 문서는 분야 전체의 지형도를 다룬다. 무차원수 자체와 영역 구분표·미끄럼 조건의 상세는 크누센수가, 대표 수치기법 하나의 내부 구현은 직접 시뮬레이션 몬테카를로가, 나비에-스토크스가 볼츠만의 1차 근사임을 보이는 전개는 채프먼-엔스코그 전개가 각각 소스 오브 트루스다.

2. 무엇을 미지수로 삼는가[편집]

연속체 역학과 희박기체역학의 진짜 차이는 방정식이 아니라 미지수다. 연속체에서는 각 점마다 $\rho$ , $\mathbf{u}$ , $T$ 라는 5개의 장(field)이면 충분하다. 이 5개로 충분한 이유는 분자 속도 분포가 국소적으로 맥스웰-볼츠만 분포에 아주 가깝고, 그 분포가 딱 5개의 모수로 결정되기 때문이다.

희박해지면 그 전제가 깨진다. 분자가 벽이나 충격파를 지나면서 겪은 이력이 충돌로 지워지기 전에 관측점에 도달하므로, 속도 분포가 두 봉우리를 갖거나 심하게 비틀린다. 그래서 미지수는 속도분포함수 $f(\mathbf{x}, \mathbf{v}, t)$ 자체가 된다 — 6차원 위상공간 위의 함수다. 거시량은 전부 이것의 모멘트로 정의된다.

\rho = m\!\int\! f\,d\mathbf{v}, \qquad \rho\mathbf{u} = m\!\int\! \mathbf{v} f\,d\mathbf{v}, \qquad \frac{3}{2}\rho R T = \frac{m}{2}\!\int\! |\mathbf{v}-\mathbf{u}|^2 f\,d\mathbf{v}

$f$ 의 시간 전개를 지배하는 것이 볼츠만 방정식이다.

\frac{\partial f}{\partial t} + \mathbf{v}\cdot\nabla_{\mathbf{x}} f + \frac{\mathbf{F}}{m}\cdot\nabla_{\mathbf{v}} f = Q(f, f)

좌변은 그냥 자유 비행이라 쉽고, 모든 어려움은 우변 충돌 적분에 있다.

Q(f,f) = \int_{\mathbb{R}^3}\!\!\int_{S^2} \left(f' f_1' - f f_1\right) g\,\sigma(g,\Omega)\,d\Omega\,d\mathbf{v}_1

한 점에서 이 값을 구하려면 5중 적분을 해야 하고, 그것을 6차원 격자의 모든 점에서, 매 시간 스텝마다 해야 한다. 게다가 이 항은 $f$ 에 대해 이차(quadratic)이고 비국소적이다. 희박기체 수치해석의 역사는 사실상 “이 적분을 어떻게 피할 것인가”의 역사다.¹

3. 충돌항을 대신하는 모형들[편집]

가장 널리 쓰인 우회로는 충돌 적분을 평형으로의 완화로 바꿔치기하는 것이다. 바트나가르-그로스-크룩(1954)의 BGK 모형이 원조다.

Q(f,f) \;\longrightarrow\; -\frac{f - f^{\mathrm{M}}}{\tau}, \qquad \tau = \frac{\mu}{p}

$f^{\mathrm{M}}$ 은 현재의 $\rho, \mathbf{u}, T$ 로 만든 국소 맥스웰 분포다. 이 한 줄로 5중 적분이 사라지고, 질량·운동량·에너지 보존과 H-정리를 모두 지킨다. 값싸고 정직한 대체재다.

문제는 프란틀수다. BGK는 점성과 열전도를 같은 완화시간 $\tau$ 하나로 만들기 때문에 $\mathrm{Pr} = 1$ 이 강제된다. 단원자 기체의 실제 값은 $2/3$ 이므로, 열유속이 계통적으로 50% 과대평가된다. 공력가열을 계산하겠다고 이 모형을 쓰면 곤란해진다. 두 갈래 수리법이 있다.

ES-BGK(Holway, 1966): 완화 대상을 등방 맥스웰이 아니라 응력 텐서를 섞은 비등방 가우스 분포로 바꾼다. 자유 모수 하나로 $\mathrm{Pr}$ 을 맞출 수 있고, 2000년에 앙드리에·르 탈레크 등이 H-정리를 증명해 열역학적 정당성까지 확보했다.
샤호프 모형(Shakhov, 1968): 맥스웰 분포에 열유속 보정항을 곱해 $\mathrm{Pr}$ 을 맞춘다. 구현이 더 간단해 널리 쓰이지만 H-정리 증명은 없다.

다원자 기체는 여기에 회전·진동 자유도가 붙는다. 러셀-플라워스, 롱-무스 같은 확장 모형이나, 병진과 내부 자유도에 서로 다른 완화시간을 주는 다중 완화 모형을 쓴다.

격자 볼츠만 방법이 여기 어디쯤인지를 정확히 아는 게 중요하다. LBM은 BGK 모형 방정식을 속도 공간에서 저차 가우스-에르미트 구적으로 이산화한 것이다. D2Q9는 3차 에르미트 정확도밖에 없어서 응력까지만 정확히 담고 열유속부터 틀리며, 결과적으로 복원되는 것은 나비에-스토크스다. 유한 $\mathrm{Kn}$ 을 노리려면 고차 격자(D3Q39, D3Q121 등)와 운동론적 벽 반사 모형을 함께 갖춰야 하고, 그래도 전이 영역 상부에서는 신뢰도가 떨어진다.

4. 방법론 지형도[편집]

계열	대표 기법	잡음	강한 영역	약점
확률적 입자	DSMC	통계 잡음 있음	고속·전이·자유분자	저속 미세유동에서 잡음 폭발
결정론적 속도격자	이산속도법, 고속 스펙트럼법	없음	저속 미세유동, 정상해	속도격자 6차원 비용
모멘트	그라드 13, R13	없음	미끄럼~전이 하부	닫힘 붕괴, 강한 충격에 취약
점근 보존	UGKS, DUGKS	없음	전 영역 이음매 없이	구현 복잡, 비용 중간
하이브리드	CFD + DSMC	접합면에 잡음	재진입 궤적 전체	접합 기준·통계 결합

**직접 시뮬레이션 몬테카를로**는 전이 영역의 사실상 표준이다. 입자를 이동과 충돌로 분리해 전진시키고 충돌 쌍을 확률적으로 뽑는데, 1992년 바그너가 입자 수 무한 극한에서 볼츠만 방정식의 해로 수렴함을 증명해 “몬테카를로 땜질”이 아님이 확정됐다. 아킬레스건은 통계 잡음이 $1/(\mathrm{Ma}\sqrt{M})$ 로 스케일한다는 것 — 마하 20 재진입은 며칠이면 끝나는데 초속 1 mm짜리 마이크로채널은 몇 주를 돌려도 잡음에 묻힌다.

**이산속도법**은 정반대의 성격이다. 속도 공간을 격자로 잘라 $f$ 를 직접 저장하고 시간 전진하므로 잡음이 원리적으로 0이다. 저속 미세유동, 크누센 펌프, 미세관 유량 곡선처럼 잡음이 신호를 덮는 문제에서 DSMC를 대체한다. 대가는 메모리다 — 3차원이면 $N_x^3 \times N_v^3$ 개의 값을 들고 있어야 하고, 고속 유동에서는 속도격자가 넓은 범위를 덮어야 해서 급격히 비싸진다. 완전한 볼츠만 충돌 적분을 유지하고 싶다면 무오·파레스키(2006)의 고속 스펙트럼법이 있는데, 푸리에-갈러킨 표현으로 5중 적분을 $\mathcal{O}(N\log N)$ 에 계산한다.

모멘트법은 $f$ 대신 유한 개의 모멘트만 시간 전진시킨다. 그라드의 13모멘트가 원형이고, 이를 정칙화한 R13 방정식(슈트루흐트루프·토릴론, 2003)이 미끄럼 영역과 전이 영역 하부에서 좋은 성적을 낸다. 다만 모멘트 계층은 반드시 어딘가에서 닫아야 하고, 그 닫힘이 깨지는 지점이 존재한다. 무한히 밀어붙이면 되지 않는 이유는 버넷 방정식의 사례가 보여 준다 — 짧은 파장 교란에 대해 선형 불안정이라 격자를 촘촘히 할수록 발산한다.

UGKS(통합 기체운동론 기법, 쉬쿤·황 2010)와 그 유한체적 사촌 DUGKS는 접근이 우아하다. 셀 경계의 플럭스를 운동론 방정식의 적분해로부터 직접 구성하는데, 이 적분해가 $\mathrm{Kn}$ 이 작으면 자동으로 나비에-스토크스 플럭스로, 크면 자유 수송 플럭스로 수렴한다. 그래서 격자와 시간간격을 평균자유행로에 맞추지 않아도 연속체 극한에서 옳은 답을 준다(점근 보존, asymptotic preserving). 재진입 궤적처럼 한 계산 안에서 $\mathrm{Kn}$ 이 몇 자릿수 변하는 문제에 특히 어울린다.

5. 기체-표면 상호작용[편집]

희박기체 계산에서 결과를 가장 많이 흔드는 것은 지배 방정식이 아니라 벽 모형이다. 분자 하나가 벽에 부딪혀 어떤 속도로 튀어나오는지에 대한 완전한 제1원리 모형이 없기 때문이다.

관례적 표준은 맥스웰(1879)의 두 갈래 혼합이다. 확률 $\sigma$ 로는 확산 반사(벽 온도의 반쪽 맥스웰 분포에서 새 속도를 뽑음, 입사 이력을 완전히 잊음), 확률 $1-\sigma$ 로는 경면 반사(법선 속도만 뒤집음, 접선 운동량 보존). $\sigma$ 를 수용계수(accommodation coefficient)라 하고, 실제 공학 표면은 흡착층 때문에 $\sigma \approx 0.8 \sim 1.0$ 으로 거의 완전 확산에 가깝다. 확산 반사에서 나오는 각도 분포가 벽 법선을 축으로 하는 크누센 코사인 법칙이며, 진공 장비의 몬테카를로 컨덕턴스 계산이 이걸 그대로 쓴다.

수용계수는 거시 경계조건으로도 번역된다. 맥스웰의 1차 속도 미끄럼 조건과 짝을 이루는 것이 스몰루호프스키의 온도 도약 조건이다.

T_g - T_w = \frac{2-\sigma_T}{\sigma_T}\,\frac{2\gamma}{\gamma+1}\,\frac{\lambda}{\mathrm{Pr}}\left.\frac{\partial T}{\partial n}\right|_w

벽에 닿은 기체의 온도가 벽 온도와 다르다는 것이 요점이다. 연속체 해석에서 당연시하던 $T_g = T_w$ 가 깨지므로 열유속 예측이 통째로 달라진다. 속도 미끄럼 쪽 표현과 열유동 미끄럼(thermal creep) 항은 크누센수 문서를 참고. 더 정교한 모형으로는 법선·접선 방향의 에너지 수용을 따로 주는 CLL(체르치냐니-람피스-로드) 모형이 있고, 원자 수준에서 분자동역학으로 산란 커널을 직접 뽑아 쓰는 시도도 있다.²

6. 자유분자 영역의 단순함[편집]

$\mathrm{Kn} > 10$ 이 되면 오히려 문제가 쉬워진다. 분자끼리 안 부딪히니 충돌항을 통째로 버릴 수 있고, 남는 것은 “벽에서 출발한 분자가 다음에 어느 벽에 닿는가”뿐이다. 즉 유체 문제가 가시성(visibility) 기하 문제로 바뀐다. 컴퓨터 그래픽스의 형태계수(view factor) 계산과 수학적으로 같은 구조이며, 실제로 진공 장비 설계자는 광선 추적과 다를 바 없는 코드를 돌린다.

이 영역의 대표 결과 몇 가지.

분출(effusion) 플럭스: 작은 구멍을 통해 빠져나가는 단위면적당 분자 수는 $\tfrac{1}{4} n \bar{c}$ 이고, $\bar{c} = \sqrt{8k_BT/\pi m}$ 이다. 압력계와 질량분석기 이온원의 기본 공식.
컨덕턴스: 긴 원형관의 자유분자 컨덕턴스는 $d^3/L$ 에 비례하고 압력과 무관하다. 연속체 영역에서 유량이 압력에 비례하던 것과 완전히 다르다. 진공 배관을 굵고 짧게 하라는 국룰이 여기서 나온다.³
자유분자 공력·가열: 경계층이라는 개념이 없다. 입사하는 분자의 운동량·에너지 플럭스와 반사 분자의 그것을 벽면 위에서 직접 더하면 항력과 열유속이 나온다. 인공위성의 궤도 감쇠(공력 항력) 계산이 이 방식이다.

문제는 그 사이의 전이 영역이다. 자유분자 해와 연속체 해를 $\mathrm{Kn}$ 의 함수로 매끄럽게 잇는 브리징 함수(오차함수 꼴 보간이 흔하다)가 실무에서 널리 쓰이는데, 이론적 근거가 약해서 결국 DSMC로 검증해야 한다.

7. 응용[편집]

진공 공정: 스퍼터링·CVD·ALD 챔버 내 기체 분포, 터보분자펌프 날개 설계, 배기 배관 컨덕턴스. 반도체 장비 설계에서 이 계산 없이 챔버 균일도를 논할 수 없다. 관련해 반도체 소자 시뮬레이션의 공정 쪽 대응물이라고 봐도 된다.
고고도 재진입 공력: 궤적 하나에서 자유분자 → 전이 → 연속체를 순서대로 다 겪는다. 고고도 구간에서는 공력계수가 연속체 값과 두 배 넘게 차이 나기도 하고, 자세 제어 스러스터의 플룸이 태양전지판을 때리는 문제도 여기 속한다. 물리는 극초음속 유동 참고.
MEMS 미세유동: 마이크로채널 유량, 가속도계·자이로의 스퀴즈 필름 감쇠, 하드디스크 헤드-디스크 간극(수 nm에서 $\mathrm{Kn} \approx 20$ ). 여기서는 속도가 느려 DSMC 잡음이 문제이므로 결정론적 기법이나 미끄럼 보정 연속체 해석이 선호된다.
우주 환경·천체: 혜성 코마, 행성 외기권 탈출, 이오·엔셀라두스 플룸. 밀도가 낮아 자유분자에 가깝지만 중력장이 있어 좌변 외력항이 살아난다.
셰일 가스·나노 기공: 나노미터 기공에서 미끄럼 때문에 겉보기 투과도가 다르시 값보다 커진다(다공성 매질 유동).

검증 케이스는 오래전에 굳어졌다. 미세관 질량유량의 크누센 최소( $\mathrm{Kn}\approx 1$ 부근에서 유량이 극소가 되는 현상), 두 평판 사이의 쿠에트·푸아죄유 유동, 그리고 충격파 내부 구조(두께가 $\lambda$ 규모라 나비에-스토크스가 두께를 계통적으로 과소평가한다)가 3대장이다.

8. 관련 문서[편집]

이 분야의 이름을 정착시킨 건 첸쉐썬(H. S. Tsien)의 1946년 논문 Superaerodynamics, Mechanics of Rarefied Gases 다. 연속체·미끄럼·전이·자유분자라는 네 구간 구분도 여기서 정리됐다. 제트추진연구소를 세운 그 사람이 맞고, 나중에 매카시즘에 밀려 중국으로 건너가 중국 우주개발의 아버지가 된 것도 그 사람이 맞다. ↩
수용계수 $\sigma$ 는 이 분야의 대표적인 “논문 뒤에 숨은 손잡이”다. 결과가 실험과 안 맞으면 $\sigma$ 를 0.9에서 0.85로 내리면 얼추 맞는데, 그 값에 물리적 근거가 있느냐고 물으면 대개 “문헌에서 그렇게 쓴다”는 답이 돌아온다. 표면을 아르곤으로 세척했는지 대기에 노출했는지에 따라 실제로 달라지는 값이라 틀린 관행이라 하기도 어렵다는 게 더 곤란하다. ↩
그래서 진공 실험실의 오래된 훈계가 “펌프를 키우지 말고 배관을 키워라”다. 직렬 컨덕턴스는 전기 회로의 병렬 저항처럼 역수 합으로 붙기 때문에, 배관 하나가 좁으면 뒤에 아무리 큰 펌프를 달아도 실효 배기속도가 그 좁은 관의 값에 눌린다. 압력과 무관하다는 성질 덕에 이 계산은 다행히 한 번만 하면 된다. ↩