버넷 방정식

유체역학 전산유체역학 수치해석 마지막 수정: 2026-08-15 04:38:52

1. 개요[편집]

버넷 방정식 Burnett Equations
유도	1935~36년, David Burnett
기반	채프먼-엔스코그 전개 Kn² 차수
겨냥한 영역	미끄럼~전이 초입 (Kn ≈ 0.01~0.1)
치명적 결함	보빌레프 불안정성 (짧은 파장 선형 불안정)
대체재	증강 버넷 · BGK-버넷 · 그라드 13 · R13
상위 대안	DSMC

버넷 방정식은 볼츠만 방정식의 채프먼-엔스코그 전개를 크누센수의 2차 항까지 밀어붙여 얻은 연속체 수송방정식이다. 1935~36년 데이비드 버넷이 유도했고,¹ 이름값과 달리 오늘날 실무 코드에서 거의 볼 수 없는 방정식이기도 하다. 이유는 정확도가 부족해서가 아니라 방정식 자체가 수치적으로 불안정하기 때문이다.

계보를 한 줄로 세우면 이렇다.

차수	분포함수	결과 방정식	수송
$\mathrm{Kn}^0$	국소 맥스웰 분포	오일러 방정식	없음 (비점성·단열)
$\mathrm{Kn}^1$	$f^{(0)} + \mathrm{Kn} f^{(1)}$	나비에-스토크스-푸리에	뉴턴 점성 · 푸리에 열전도
$\mathrm{Kn}^2$	$+\,\mathrm{Kn}^2 f^{(2)}$	버넷	2차 비평형 응력·열유속
$\mathrm{Kn}^3$	$+\,\mathrm{Kn}^3 f^{(3)}$	초버넷(super-Burnett)	—

여기서 자연스럽게 드는 기대는 “차수를 올렸으니 더 정확하겠지”다. 버넷 방정식은 그 기대가 왜 순진한지를 가르쳐 주는, 전산유체역학에서 가장 유명한 반례다. 채프먼-엔스코그 급수는 점근 급수이지 수렴 급수가 아니라서, 항을 더한다고 답에 가까워진다는 보장이 없다.

2. 2차 항의 모양과 물리[편집]

버넷 차수에서 응력텐서와 열유속에는 나비에-스토크스에 없던 항들이 무더기로 붙는다. 구조만 보면 이렇다. $\mathbf{S}$ 를 변형률텐서, $\boldsymbol{\Omega}$ 를 회전텐서라 하면, 2차 응력은 $\mu^2/p$ 규모의 계수를 앞에 달고 다음 여섯 종류를 각각 대각합을 제거한 대칭 텐서로 만들어 더한 꼴이다.

\sigma^{(2)}_{ij} \;\sim\; \frac{\mu^{2}}{p}\Big[\, (\nabla\!\cdot\!\mathbf{u})\,S_{ij},\;\; S_{ik}S_{kj},\;\; S_{ik}\Omega_{kj},\;\; \frac{1}{T}\partial_i\partial_j T,\;\; \frac{1}{T^{2}}\,\partial_i T\,\partial_j T,\;\; \frac{1}{\rho T}\,\partial_i p\,\partial_j T \,\Big]

열유속 쪽에는 다섯 종류가 붙는다. 대략 $\partial_i(\nabla\!\cdot\!\mathbf{u})$ , $S_{ij}\partial_j T$ , $S_{ij}\partial_j p$ 같은 조합들이다. 각 항 앞의 계수 $\varpi_1 \dots \varpi_6$ , $\theta_1 \dots \theta_5$ 는 순수한 수이고, 분자 간 상호작용 모형에 따라 값이 정해진다(맥스웰 분자에서는 해석적으로 떨어지고, 강체구는 수치적으로 구한다).

항 목록보다 중요한 건 거기 담긴 물리다.

열응력(thermal stress). $\partial_i\partial_j T$ 와 $\partial_i T\,\partial_j T$ 항은 속도가 0이어도 살아 있다. 즉 정지한 기체 안에 온도장만 불균일하면 응력이 생긴다. 나비에-스토크스에서는 절대 나올 수 없는 결과다. 결과적으로 온도가 불균일한 밀폐 용기 안의 기체는 원리적으로 정지 상태에 머무를 수 없고, 이 열응력이 만드는 미약한 대류를 열응력 유동이라 부른다. 운동론 계산과 DSMC가 실제로 이 흐름을 재현한다. 크룩스 라디오미터의 날개가 도는 원리 계열(열가장자리 유동)이 같은 동네다.
비푸리에 열유속. 열유속이 더 이상 $-k\nabla T$ 가 아니다. 압력 기울기나 변형률만으로도 열이 흐른다. 온도차 없이 열이 흐르는 게 이상해 보이지만, 열역학 제2법칙을 어기는 건 아니고 단지 “온도 기울기가 유일한 구동력”이라는 것이 1차 근사였을 뿐이다.
비등방 수직응력. 강한 압축 영역에서 $\sigma_{xx} \neq \sigma_{yy}$ 가 뚜렷해진다. 충격파 내부 구조가 나비에-스토크스보다 실제에 가까워지는 원인이 이것이다.

그래서 버넷 항은 미끄럼 경계조건의 2차 항, 열전이(thermal transpiration), 크누센 펌프처럼 “나비에-스토크스로는 설명이 안 되는데 실험에서는 보이는” 현상들의 이론적 근거를 제공한다. 이 점에서 버넷 방정식은 계산 도구로 실패했어도 개념 도구로는 살아남았다.

3. 보빌레프 불안정성[편집]

1982년 보빌레프가 선형 안정성 해석으로 못을 박았다. 균일 평형 상태에 파수 $k$ 의 미소 교란을 넣고 버넷 방정식을 선형화하면, $k$ 가 어느 문턱을 넘는 순간 증폭률이 양수가 되고, 그 뒤로는 $k$ 가 커질수록 더 빨리 자란다.² 물리적 소산이 짧은 파장을 죽여야 정상인데 반대로 키우는 것이다.

문턱 파장은 대략 평균자유행로 규모다. 이게 실무에서 무슨 뜻인가 하면,

격자 간격이 평균자유행로보다 커서 짧은 파장 모드를 애초에 표현하지 못하는 동안은 계산이 멀쩡하게 굴러간다.
격자를 촘촘하게 하면 터진다. 격자 수렴 연구를 돌리다가 가장 촘촘한 격자에서만 발산하는 상황을 만나면, 보통은 코드를 의심하지만 여기서는 방정식을 의심해야 한다.

원인은 절단 자체에 있다. 점근 급수를 유한 차수에서 자르면 잘려 나간 고차 항이 담당하던 안정화가 함께 사라지고, 남은 2차 미분 항들이 잘못된 부호로 작동한다. 여기에 더해 특정 조건에서 엔트로피 생성이 음수가 될 수 있다는 것도 알려져 있다. 즉 버넷 방정식은 H-정리와 무조건적으로 양립하지 않는다.

부수적인 실무 장벽도 있다. 방정식에 3계 이상의 공간 미분이 등장하는데, 그 미분들에 대한 추가 경계조건을 운동론이 자연스럽게 제공하지 않는다. 벽에서 무엇을 지정해야 하는지 자체가 미결이라, 경계가 지배하는 미세유동 문제에서는 방정식이 안정해도 문제가 안 닫힌다.

4. 고치려는 시도들[편집]

불안정성이 알려진 뒤로 이 계열을 살리려는 시도가 여러 갈래로 갈렸다.

관례적 버넷 vs 원형 버넷. 채프먼-엔스코그 유도 도중 $f^{(2)}$ 에 나타나는 시간 미분을 낮은 차수(오일러) 방정식으로 치환하면 관례적(conventional) 버넷, 그대로 남기면 원형(original) 버넷이다. 열역학적으로 동등한 차수인데 수치 안정성은 다르게 나온다. 어느 형태를 썼는지 밝히지 않은 논문끼리 결과가 안 맞는 고전적 원인.
증강 버넷(augmented Burnett). 종·매코맥·채프먼이 1993년에 제안한 것으로, 초버넷 차수(3차)의 항 일부를 골라 되돌려 넣어 선형 안정성을 회복시킨다. 차수 논리로는 정당화가 약하지만(“2차 방정식에 왜 3차 항을?”) 극초음속 충격파 구조 계산에서 실제로 돌아간다는 점이 강점이다.
BGK-버넷. 발라크리슈난과 아가르왈이 BGK 충돌 모형에서 출발해 유도한 버전으로, 엔트로피 부등식과의 정합성을 설계 목표에 넣었다. 다만 BGK 모형 자체가 프란틀수를 1로 고정하는 결함을 물려받는다.
모멘트 계열로의 이탈. 근본적인 해법은 전개를 더 미는 대신 미지수를 늘리는 것이다. 1949년 그라드가 제안한 그라드 13 모멘트 방정식은 응력텐서와 열유속을 미분으로 표현하지 않고 독립 변수로 승격시켜 각자의 수송방정식을 준다. 계가 쌍곡형이라 안정적이지만, 마하수가 대략 1.65를 넘으면 충격파 내부에 물리적으로 존재하지 않는 불연속(부수 충격, subshock)이 생기고 크누센층도 못 잡는다. 슈트루흐트루프와 토릴론이 2003년 제시한 **R13 방정식**은 그라드 13에 정칙화 항을 얹어 이 두 문제를 동시에 해결했고, 흥미롭게도 그 정칙화가 만들어 내는 항들이 버넷 항을 포함한다. 즉 R13은 “버넷의 물리는 살리고 불안정성은 버린” 절충안에 가깝다. 뒤이어 R26 등 더 높은 모멘트 계열도 나왔다.

5. 어디까지 쓰고 어디서 갈아타는가[편집]

$\mathrm{Kn} < 0.01$ — 나비에-스토크스로 충분하다. 버넷 항은 계산 비용만 늘린다.
$\mathrm{Kn} \approx 0.01 \sim 0.1$ — 미끄럼 영역. 실무의 첫 선택은 여전히 나비에-스토크스 + 미끄럼/온도 도약 경계조건이다. 버넷 항이 이 구간에서 이론적으로 우월한 건 맞지만, 안정화 장치와 추가 경계조건을 얹는 수고 대비 이득이 크지 않다.
$\mathrm{Kn} \approx 0.1 \sim 1$ — 여기가 버넷 계열(증강 버넷·R13)이 그나마 존재 이유를 갖는 구간이다. 전이 영역 초입에서 DSMC보다 훨씬 싸게 답을 준다.
$\mathrm{Kn} > 1$ — 그냥 DSMC로 간다. 전개를 더 미는 것보다 운동론을 직접 푸는 게 빠르고 정확하다.

충격파 구조가 이 계열의 표준 시험장이다. 나비에-스토크스는 마하 2만 넘어도 충격파를 실제보다 얇게 예측하는데, 버넷 항을 넣으면 눈에 띄게 개선된다. 다만 마하수를 더 올리면 다시 어긋나고, 어느 지점부터는 정상해가 존재하지 않는 문제까지 생긴다. 결국 극초음속 재진입 해석의 정답은 DSMC이거나, 연속체 영역과 희박 영역을 나눠 붙이는 하이브리드 접근이다. 후자에서는 국소 크누센수로 영역을 판정해 경계면에서 두 솔버를 연결한다.

수치 구현에서 통용되는 요령 하나. 버넷 항을 지배 방정식에 완전히 넣어 음함수로 푸는 대신, 나비에-스토크스 솔버의 소스항으로 취급해 직전 스텝 해에서 명시적으로 평가한다. 구현이 훨씬 간단할 뿐 아니라 고차 미분 항이 암묵적으로 지연되면서 실질적인 안정화 효과까지 얻는다. 물론 이건 원리적 해결이 아니라 그냥 잘 듣는 반칙에 가깝다.

6. 왜 여전히 배우는가[편집]

버넷 방정식이 코드에서 밀려난 뒤에도 교과서에서 안 없어지는 이유는 세 가지다.

연속체 서술의 한계선을 정확히 보여 준다. “나비에-스토크스가 왜 안 되는가”에 대한 답은 대개 “고차 항이 있어서”인데, 그 고차 항이 실제로 어떻게 생겼고 무슨 물리를 담고 있는지를 명시적으로 보여 주는 것이 버넷 차수다.
2차 미끄럼 조건과 크누센층 논의의 출발점이다. 미세유동 문헌에서 계수가 문헌마다 다르다고 악명 높은 그 2차 항들이 여기서 나온다.
점근 전개의 위험성에 대한 교보재. 급수의 다음 항을 더하는 것이 항상 개선이 아니라는 사실을, 이만큼 극적으로 보여 주는 사례가 드물다. 같은 교훈은 격자 볼츠만 방법의 고차 격자 설계나 희박기체역학 전반에서 반복해서 등장한다.³

7. 관련 문서[편집]

Burnett, D. (1935~36). Proc. London Math. Soc. 두 편의 논문으로 나뉘어 있다. 정작 본인은 유체 코드를 짤 일이 없었으므로 이 방정식이 90년 뒤 “터지는 방정식”의 대명사가 될 줄은 몰랐을 것이다. 방정식에 자기 이름이 붙는 게 항상 좋은 일만은 아니라는 교훈. ↩
Bobylev, A. V. (1982). 이 결과 때문에 이 불안정성을 보빌레프 불안정성이라 부른다. 재미있는 건 그가 같은 논문에서 그라드 모멘트 방법 쪽은 이런 문제가 없다는 것도 함께 보였다는 점이다. 이후 30년간 희박기체 커뮤니티가 모멘트 계열로 기운 데는 이 한 편의 지분이 크다. ↩
흔한 오해 하나. “버넷 방정식은 크누센수가 커서 틀리는 것”이 아니라 크누센수가 작아도 격자만 촘촘하면 터진다. 불안정성은 물리적 파라미터가 아니라 표현 가능한 최소 파장이 결정한다. 그래서 “우리 문제는 Kn이 작으니 안전하다”는 논리는 안 통하고, 정말로 안전한 유일한 조건은 “격자를 평균자유행로보다 굵게 유지한다”인데 그건 애초에 왜 버넷을 쓰는지를 부정하는 조건이다. ↩