버거스 방정식

수치해석 전산유체역학 마지막 수정: 2026-08-24 04:12:41

1. 개요[편집]

버거스 방정식 Burgers' equation
형태	u_t + u u_x = ν u_xx
이름	요하네스 마르티뉘스 뷔르허스 (1939·1948)
구성 요소	비선형 이류 + 선형 확산
ν → 0 극한	스칼라 쌍곡형 보존법칙, 충격파 형성
정확해	콜-호프 변환으로 열방정식과 동치
주 용도	수치기법 검증의 국룰 시험대

나비에-스토크스에서 압력과 3차원을 빼고 남은 뼈다귀. 그런데 이 뼈다귀는 정확해가 있다.

버거스 방정식(Burgers’ equation)은 비선형 이류항과 선형 확산항만을 남긴 1차원 편미분방정식으로,

\frac{\partial u}{\partial t} + u\,\frac{\partial u}{\partial x} = \nu\,\frac{\partial^2 u}{\partial x^2}

의 꼴을 갖는다. $\nu=0$ 이면 무점성 버거스 방정식, $\nu>0$ 이면 점성 버거스 방정식이라 부른다.

나비에-스토크스 방정식에서 압력항·비압축 구속·횡방향 성분을 전부 지우면 정확히 이 식이 남는다. 즉 유체 방정식의 두 핵심 성질 — 자기 자신을 실어 나르는 비선형성과 평활화하는 확산 — 만 골라 담은 최소 모형이다. 그런데 놀랍게도 이 조합은 콜-호프 변환으로 열방정식이 되어 버려서, 임의의 초기조건에 대한 정확해가 존재한다.

비선형 충격을 만들면서도 정답을 알고 있다는 이 조합은 수치해석자에게 축복이다. 새 이산화 도식을 만들면 일단 버거스에 던져 보는 것이 이 바닥의 국룰이고, 검증 및 확인에서 “충격이 있는 문제의 수렴 차수”를 재는 표준 시험문제가 대체로 버거스다.¹

2. 무점성 극한 — 특성선이 서로를 덮친다[편집]

$\nu=0$ 이면 버거스는 $u$ 를 스스로의 전파 속도로 갖는 스칼라 이류 방정식이다. 특성곡선법을 적용하면 $dx/dt = u$ 위에서 $du/dt = 0$ , 즉 각 유체 입자가 자기 초기값을 그대로 들고 등속 직진한다.

u(x,t) = u_0(\xi), \qquad x = \xi + u_0(\xi)\,t

빠른 입자가 느린 입자를 뒤에서 따라잡으므로 특성선이 언젠가 교차한다. 교차 시각은 위 사상 $\xi \mapsto x$ 의 야코비안 $1 + u_0'(\xi)t$ 가 처음 0이 되는 시각이며

t_b = -\frac{1}{\min_\xi u_0'(\xi)}

이다. $u_0'$ 가 음수인 지점, 즉 속도가 감소하는 구간(압축 구간)이 하나라도 있으면 유한 시간에 $u_x \to -\infty$ 가 되고 고전해가 죽는다. 매끈한 사인파 $u_0=\sin x$ 면 $\min u_0' = -1$ 이므로 $t_b=1$ — 학부 과제에서 ” $t=1$ 이후로는 왜 코드가 이상해지나요”의 정체가 이것이다.

$t>t_b$ 에서는 다가 함수를 잘라내고 약해(weak solution)를 받아들여야 한다. 그러려면 방정식을 반드시 보존형으로 써야 한다.

u_t + \left(\frac{u^2}{2}\right)_x = 0

플럭스 $f(u)=u^2/2$ 가 볼록하므로 버거스는 볼록 스칼라 보존법칙의 교과서적 대표 선수다.

3. 랭킨-위고니오와 엔트로피[편집]

불연속을 가로지르는 적분형 보존을 요구하면 랭킨-위고니오 조건이 나온다. 충격 속도 $s$ 는

s = \frac{f(u_L)-f(u_R)}{u_L-u_R} = \frac{u_L+u_R}{2}

즉 좌우 값의 산술평균이다. 버거스 충격이 유난히 다루기 편한 이유가 이 한 줄이다.

문제는 랭킨-위고니오만으로는 답이 하나로 정해지지 않는다는 것. $u_L < u_R$ 인 초기 불연속에 대해 “그 자리에 그대로 서 있는 불연속”도, “부채꼴로 벌어지는 희박팬”도 둘 다 약해다. 물리적으로 옳은 쪽을 고르는 것이 엔트로피 조건이며, 볼록 플럭스에서는 특성선이 충격으로 들어가야 한다는 락스(Lax) 조건

u_L > s > u_R

으로 간단히 표현된다. 압축(좌 > 우)만 충격이고 팽창(좌 < 우)은 희박팬이다. 이 판정을 어기는 해를 엔트로피 위반 해라 부르고, 고두노프 도식의 소박한 변형이나 로(Roe) 근사 리만 솔버가 음속점(sonic point, $u_L<0<u_R$ )에서 만들어 내는 “충격처럼 서 있는 팽창 불연속”이 그 대표 사례다. 엔트로피 픽스가 필요한 이유를 설명할 때 예제로 등장하는 것이 거의 항상 버거스인 것도 그 때문.²

물리적으로 옳은 해는 점성 소실 극한(vanishing viscosity)으로도 정의된다. $\nu>0$ 인 해의 $\nu\to0^+$ 극한이 유일한 엔트로피 해라는 정리가 있고, 이 사실이 버거스를 “이론과 계산이 정확히 맞물리는 몇 안 되는 문제”로 만든다.

4. 콜-호프 변환 — 비선형을 선형으로 접기[편집]

점성 버거스가 특별한 이유는 다음 치환 때문이다.

u = -2\nu\,\frac{\varphi_x}{\varphi}

이걸 대입하면(중간 단계에서 $u=\psi_x$ 로 퍼텐셜을 도입해 한 번 적분한다) $\varphi$ 가 만족하는 식은

\varphi_t = \nu\,\varphi_{xx}

— 그냥 열방정식이다. 초기조건도 같은 변환으로 옮기면 되고, 열방정식은 가우스 핵으로 명시적으로 풀리므로

u(x,t) = \frac{\displaystyle\int_{-\infty}^{\infty} \frac{x-\xi}{t}\, \exp\!\left[-\frac{1}{2\nu}G(\xi;x,t)\right] d\xi} {\displaystyle\int_{-\infty}^{\infty} \exp\!\left[-\frac{1}{2\nu}G(\xi;x,t)\right] d\xi}, \qquad G = \int_0^\xi u_0(s)\,ds + \frac{(x-\xi)^2}{2t}

라는 적분 표현이 나온다. 남는 계산은 수치적분 하나뿐이라, 참값을 기계 정밀도 수준으로 뽑아 놓고 도식의 오차를 재는 것이 가능하다. 새 도식의 $L^1$ 수렴 차수를 주장할 때 버거스가 동원되는 실질적 이유가 여기 있다.

$\nu\to0$ 에서는 위 적분이 최급강하법으로 평가되어 $G$ 를 최소화하는 $\xi$ 가 지배하는데, 그 최소점이 둘로 갈리는 지점이 정확히 충격 위치다. 콜-호프는 엔트로피 해 선택 원리를 적분 하나에 내장하고 있는 셈이다.³

5. 정상 충격 프로파일과 N파[편집]

테일러 충격. 점성이 있으면 불연속 대신 유한한 두께의 진행파가 선다. $u\to u_1$ ( $x\to-\infty$ ), $u\to u_2$ ( $x\to+\infty$ ), $u_1>u_2$ 인 진행파 해는

u(x,t) = \frac{u_1+u_2}{2} - \frac{\Delta u}{2}\tanh\!\left[\frac{\Delta u\,(x-st)}{4\nu}\right], \qquad \Delta u = u_1-u_2,\; s=\frac{u_1+u_2}{2}

이다. 속도는 랭킨-위고니오와 정확히 같고, 두께는 $\delta \sim 4\nu/\Delta u$ — 점성에 비례하고 충격 강도에 반비례한다. 이 스케일이 실무에서 중요한 이유는 격자 요구조건을 직접 규정하기 때문이다. 충격을 “해상”하려면 $\Delta x \lesssim \nu/\Delta u$ 가 필요하고, 그게 불가능해서 충격을 포착(capturing)하기로 타협하는 것이 유한체적법 기반 압축성 코드의 세계관 전체다.

N파와 자기상사. 유한한 구간에만 실려 있던 초기 교란은 시간이 지나면 초기 형상을 잊고 보편적인 꼴로 수렴한다.

한 봉우리(면적 $A=\int u\,dx>0$ ) — 무점성 극한에서 삼각파 $u = x/t$ ( $0<x<\sqrt{2At}$ )로 간다. 진폭 $\sqrt{2A/t} \propto t^{-1/2}$ , 폭 $\propto t^{1/2}$ . 면적 $A$ 가 보존되므로 국소 레이놀즈수 $R=A/2\nu$ 가 상수로 유지되고, 이 자기상사는 영원히 깨지지 않는다.
N파(양·음 두 로브, 총면적 0) — 로브 면적이 보존되지 않고 서서히 새어나간다. 비선형 단계에서는 마찬가지로 진폭 $t^{-1/2}$ ·폭 $t^{1/2}$ 이지만, 국소 레이놀즈수가 1 아래로 떨어지는 순간 비선형이 죽고 선형 확산 지배로 넘어가 진폭이 $t^{-1}$ 로 붕괴한다. 초음속 비행체의 소닉붐이 지면까지 오는 동안 왜 저 모양으로 늘어지는지가 이 해의 이야기다.

6. 버거스 난류[편집]

$u_t + uu_x = \nu u_{xx} + f(x,t)$ 처럼 백색잡음에 가까운 강제를 넣으면 버거스 난류(burgulence)라는 별종 난류가 나온다. 실제 3차원 난류와 결정적으로 다른 점은 소용돌이 늘림이 없다는 것이지만, 그 대신 충격이라는 강한 특이구조가 통계를 지배하는 세계가 열린다.

에너지 스펙트럼 $E(k)\propto k^{-2}$ . 점프 불연속의 푸리에 계수가 $1/k$ 로 감쇠하므로 스펙트럼은 자동으로 $k^{-2}$ 가 된다. 콜모고로프의 $k^{-5/3}$ 과 달리 유도에 차원해석이 필요 없다.
극단적 간헐성. 구조함수 $S_p(\ell)=\langle|\delta u|^p\rangle$ 가 $p\ge1$ 에서 $\ell^1$ 로, $p<1$ 에서 $\ell^p$ 로 스케일한다. 지수가 두 직선으로 꺾이는 이중 프랙탈(bifractal)이며, 실제 난류의 완만하게 휘는 다중 프랙탈과 대비된다. 매끈한 구간과 충격이라는 두 종류의 구조밖에 없으니 당연한 결과.
KPZ 방정식과의 동치. $u = -h_x$ 로 두면 강제 버거스는 정확히 1차원 KPZ 방정식이 된다. 계면 성장 문제와 1차원 압축성 난류가 같은 방정식이라는 사실은 통계물리에서 꽤 유명한 다리다.

우주론에도 같은 방정식이 나온다. 대규모 구조 형성의 젤도비치 근사는 입자가 초기 속도로 자유 비행하는 무점성 버거스 그 자체이고, 여기에 인공 점성을 넣어 교차한 입자들이 서로 달라붙게 만든 것이 접착 모형(adhesion model)이다. 그때 생기는 충격 그물망이 은하 필라멘트 구조에 대응한다.

7. 수치적으로 버거스가 가르치는 것[편집]

보존형이 아니면 충격 속도가 틀린다. 이 교훈을 가장 값싸게 배우는 무대가 버거스다. 비보존형 $u_t + u\,u_x = 0$ 을 그대로 이산화하면

u_j^{n+1} = u_j^n - \frac{\Delta t}{\Delta x}\,u_j^n\,(u_j^n - u_{j-1}^n)

같은 식이 나오는데, 매끈한 영역에서는 멀쩡하게 1차 수렴하다가 충격이 생기는 순간 충격이 엉뚱한 속도로 간다. 격자를 아무리 줄여도 안 고쳐진다 — 수렴하기는 하는데 틀린 답으로 수렴한다. 반면 보존형

u_j^{n+1} = u_j^n - \frac{\Delta t}{\Delta x}\left(F_{j+1/2} - F_{j-1/2}\right)

으로 짜면 이산 보존 덕분에 충격 속도가 자동으로 랭킨-위고니오를 만족한다(락스-벤드로프 정리). 초심자가 “ $uu_x$ 랑 $(u^2/2)_x$ 는 미분하면 같은데요?”라고 항변하는 지점이고, 답은 “매끈한 해에서만 같다”다.⁴

나머지 교훈도 대부분 버거스에서 처음 만난다.

현상	버거스에서 보는 방법
수치 소산	1차 풍상차분의 충격 두께가 $\nu$ 대신 $\Delta x$ 로 정해진다
수치 분산·진동	중심차분·락스-벤드로프가 충격 뒤에 링잉을 만든다
단조성 한계	고두노프 정리 — 선형 2차 이상 도식은 단조를 못 지킨다
해법	전변분 감소 제한자, MUSCL, WENO 도식
시간 간격 제약	이류는 CFL 조건, 확산은 $\nu\Delta t/\Delta x^2 \le 1/2$ — 둘 중 빡센 쪽이 이긴다

$\nu$ 가 작아질수록 확산 안정조건이 잔인해지므로, 실무 코드는 이류를 명시적으로·확산을 암시적으로 처리하는 IMEX나 연산자 분리를 쓴다. 주기 영역이라면 의사스펙트럼법으로 확산을 정확 적분하고 비선형항만 명시적으로 밟는 것이 표준인데, 이때 $uu_x$ 의 앨리어싱을 3/2 규칙으로 잘라내지 않으면 계산이 조용히 폭발한다.

버거스와 형제인 방정식들도 같은 자리에서 배운다. 확산 대신 분산 $u_{xxx}$ 을 넣으면 KdV 방정식이 되어 충격 대신 솔리톤이 나오고, 음의 확산 $-u_{xx}$ 에 4차 감쇠 $-u_{xxxx}$ 를 얹으면 쿠라모토-시바신스키 방정식이 되어 시공간 카오스로 간다. 부호와 차수 하나 차이로 완전히 다른 세계가 열리는 셈. 확산항 대신 압력·에너지를 얹어 계로 확장하면 그것이 오일러 방정식이고, 그때부터는 리만 솔버 이야기가 시작된다.

8. 관련 문서[편집]

이름은 네덜란드 유체역학자 요하네스 마르티뉘스 뷔르허스(J. M. Burgers)에서 왔다. 다만 이 방정식이 뷔르허스보다 먼저 해리 베이트먼(1915)의 논문에 등장했다는 점은 문헌에 잘 정리되어 있어서, 요즘은 예의상 베이트먼-버거스 방정식이라고 부르는 논문도 종종 보인다. 마태 효과의 표본. ↩
음속점 엔트로피 위반은 정확히 “격자 한 칸짜리 계단이 이론적으로 정상 상태가 되어 버리는” 사고다. 도식이 잘못 계산한 게 아니라, 그 도식이 허용하는 정상해가 하나 더 있는 것이다. 그래서 격자를 줄여도 계단은 사라지지 않고 그냥 얇아진다. 잔차 그래프는 아주 예쁘게 수렴한다 — 틀린 답으로. ↩
그렇다고 콜-호프가 만능은 아니다. $\nu$ 가 아주 작으면 지수 안의 $1/2\nu$ 가 폭발해 배정밀도로는 적분이 언더플로한다. “정확해”를 뽑겠다고 짠 코드가 정작 관심 있는 고레이놀즈 영역에서 NaN을 뱉는 광경은 이 바닥의 흔한 통과의례다. 로그 영역에서 최댓값을 빼고 계산하는 log-sum-exp 트릭이 여기서도 답이다. ↩
이 사실을 정리로 못 박은 것이 후·르플로슈(Hou & LeFloch, 1994)의 결과다. 비보존형 도식의 극한은 원래 보존법칙의 해가 아니라 소스항이 붙은 측도값 해가 되고, 그 소스가 충격 위치에 델타로 얹혀 속도를 틀어 놓는다. “틀린 답으로 정확히 수렴한다”는 표현이 문학적 과장이 아니라는 뜻. ↩