웰밸런스드 도식

1. 개요[편집]

웰밸런스드 도식 Well-balanced scheme
대상	소스항이 있는 쌍곡계 $q_t + f(q)_x = s(q,x)$
요구	정지 평형해를 이산 수준에서 정확히 보존
표준 시험	lake at rest — $h+b=\text{const}$, $u=0$
성질 이름	C-property (Bermúdez & Vázquez-Cendón, 1994)
대표 처방	정수압 재구성 · 표면 경사법 · 확장 리만 솔버
동반 요구	양의 수심 보존, 마름-젖음 처리

아무 일도 일어나지 않는 상황을 정확히 재현하는 것이, 무언가 일어나는 상황을 재현하는 것보다 어렵다.

웰밸런스드 도식(well-balanced scheme)은 소스항이 있는 쌍곡형 보존계에서, 플럭스 기울기와 소스항이 정확히 상쇄되어 성립하는 정지 평형해를 이산 수준에서도 오차 없이 보존하도록 설계한 수치 도식이다. 원래의 보존계

$$\frac{\partial q}{\partial t} + \frac{\partial f(q)}{\partial x} = s(q,x)$$

에서 좌변 둘째 항과 우변이 서로 상쇄하는 상태가 정상해인데, 문제는 이 둘을 서로 다른 방식으로 이산화한다는 데 있다. 플럭스는 셀 경계에서 리만 솔버로, 소스는 셀 중심에서 중앙차분으로 — 각각은 얌전히 수렴하지만 둘의 차가 0이 되어야 할 자리에서 0이 되지 않는다.

가장 유명한 사례가 얕은 물 방정식의 정지한 호수(lake at rest)다. 지형 $b(x)$가 아무리 울퉁불퉁해도 수면이 수평이고($h+b=\text{const}$) 물이 멈춰 있으면($u=0$) 그것은 정확한 정상해여야 한다. 실제로 운동량 방정식의 정상 조건은

$$\frac{\partial}{\partial x}\!\left(\tfrac{1}{2}gh^2\right) = -gh\frac{\partial b}{\partial x} \;\Longleftrightarrow\; gh\frac{\partial (h+b)}{\partial x} = 0$$

으로 깔끔하게 성립한다. 그런데 소박하게 이산화한 도식에 이 초기조건을 넣고 돌리면, 가만히 있어야 할 물이 출렁이기 시작한다.

2. 가짜 유속이 왜 치명적인가[편집]

불균형이 만드는 가짜 유속의 크기는 도식의 절단오차 크기다. 1차 도식이면 $O(\Delta x)$, 2차면 $O(\Delta x^2)$ — 격자를 줄이면 줄어들기는 한다. 문제는 그 크기가 우리가 잡으려는 신호보다 크다는 것이다.

심해를 지나는 해일의 진폭은 수십 cm인데 수심은 4 km다. 상대 크기로 $10^{-4}$ 수준이다. 지형 기복이 수백 m 단위로 널려 있는 해저에서 불균형 오차가 지형 경사에 비례해 생기면, 그 가짜 신호가 잡으려는 파보다 커진다. 하천 범람 계산도 마찬가지여서, 실제 수심 변화 몇 cm를 보려는데 하상 기복이 만든 가짜 유속이 m/s 단위로 나오면 계산 결과 전체가 무의미해진다.

게다가 이 오차는 누적된다. 한 스텝의 불균형은 작지만, 정상해 주위에서 그것이 계속 같은 방향으로 운동량을 주입하면 유속이 자라고, 자란 유속이 자유표면을 흔들고, 흔들린 표면이 다시 불균형을 키운다. 지형이 험한 영역에서 며칠짜리 계산을 돌리면 “초기조건에서 아무것도 안 했는데 도메인 전체가 끓고 있는” 결과를 만나게 된다.¹

3. C-property — 무엇을 정확히 보존한다는 말인가[편집]

베르무데스와 바스케스-센돈이 1994년에 이 요구를 형식화하면서 C-property(conservation property)라는 이름을 붙였다. 정의는 두 단계다.

• 정확한 C-property(exact): 임의의 지형 $b(x)$에 대해, 정지 평형 초기조건을 넣으면 도식이 그 상태를 반올림 오차 수준까지 정확히 유지한다. 즉 이산 잔차가 대수적으로 0이다.
• 근사 C-property(approximate): 잔차가 0은 아니지만 도식의 절단오차보다 높은 차수로 작다. 예컨대 1차 도식인데 정지 상태 잔차가 $O(\Delta x^2)$면 실용상 문제가 없다.

“웰밸런스드”라고 부를 때는 보통 exact 쪽을 뜻한다. 중요한 것은 이 성질이 정확도 차수와 독립이라는 점이다. 2차 정확도 도식이 웰밸런스드가 아닐 수 있고, 1차 도식이 웰밸런스드일 수 있다. 격자 수렴 시험은 둘을 구별해 주지 못하므로, 별도의 시험을 따로 돌려야 한다.

또 하나 짚어야 할 것은 어떤 평형을 보존하느냐다. 정지수(still water) 평형만 보존하는 것이 표준이고, 물이 흐르는 채로 정상인 상태(moving equilibrium)까지 보존하는 것은 훨씬 어려운 별개의 문제다(아래 참조).

4. 정수압 재구성[편집]

가장 널리 쓰이는 처방이 오뒤스·부슈·브리스토·클라인·페르탐(2004)의 정수압 재구성(hydrostatic reconstruction)이다. 발상은 “소스항을 잘 이산화하자”가 아니라 **“소스항을 플럭스 안으로 밀어 넣자”**에 가깝다.

셀 $i$와 $i+1$ 사이의 경계에서 다음을 정의한다.

$$b_{i+1/2} = \max(b_i,\, b_{i+1}), \qquad h^{*}_{L} = \max\!\left(0,\; h_i + b_i - b_{i+1/2}\right), \qquad h^{*}_{R} = \max\!\left(0,\; h_{i+1} + b_{i+1} - b_{i+1/2}\right)$$

즉 경계에서 더 높은 쪽 바닥을 기준으로 삼아 양쪽 수심을 다시 잰다. 이 재구성된 상태 $U^{*}_{L}, U^{*}_{R}$로 리만 솔버를 호출하고, 원래 셀 값과의 차이를 정수압 보정항으로 되돌려 준다.

$$\mathcal{F}^{L}_{i+1/2} = F(U^{*}_{L}, U^{*}_{R}) + \begin{pmatrix} 0 \\[2pt] \tfrac{1}{2}g h_i^2 - \tfrac{1}{2}g (h^{*}_{L})^2 \end{pmatrix}$$

오른쪽 셀에는 $h_{i+1}, h^{*}_{R}$로 만든 대응 보정항을 쓴다. 왜 이게 통하는지는 정지 상태를 넣어 보면 한 줄에 보인다. $h+b=\eta$가 일정하면 $h^{*}_{L} = h^{*}_{R} = \max(0,\eta - b_{i+1/2})$로 좌우가 같아지므로 리만 문제가 자명해지고, 리만 솔버는 정확한 플럭스 $(0,\ \tfrac{1}{2}g(h^{*})^2)$를 돌려준다. 보정항을 더하면 $\mathcal{F}^{L}_{i+1/2} = (0,\ \tfrac{1}{2}gh_i^2)$이고, 셀 $i$의 왼쪽 경계에서도 같은 계산이 $(0,\ \tfrac{1}{2}gh_i^2)$를 준다. 두 플럭스의 차가 정확히 0 — 소스항을 따로 이산화할 필요조차 없이 균형이 대수적으로 성립한다.

이 재구성의 또 다른 미덕은 $\max(0,\cdot)$ 덕분에 수심이 음수가 되지 않는다는 것이다. 웰밸런스드와 양의 수심 보존을 한 장치로 동시에 얻는 셈이고, 그래서 범람 계산의 사실상 기본값이 됐다.

2차 이상으로 올릴 때는 셀 안에서 $h$가 아니라 수면고도 $\eta = h+b$를 재구성한 뒤 위 절차를 경계값에 적용하고, 셀 내부 지형 경사에 대한 중앙형 소스항을 하나 더 붙인다. 이 조합이 exact C-property를 유지한다는 것이 원논문의 결과다.

5. 표면 경사법[편집]

저우·코슨·밍엄·잉그램(2001)의 표면 경사법(surface gradient method, SGM)은 훨씬 소박하지만 유한체적법 코드에 얹기 쉬워 널리 쓰인다. 핵심은 재구성 변수를 바꾸는 것 하나다.

• 기존: 셀 평균 $h_i$를 MUSCL 방식으로 선형 재구성 → 경계값 $h_{i+1/2}^{\pm}$.
• SGM: 셀 평균 $\eta_i = h_i + b_i$를 선형 재구성 → 경계값 $\eta^{\pm}_{i+1/2}$ → $h^{\pm} = \eta^{\pm} - b_{i+1/2}$.

정지 상태에서는 $\eta$가 상수라 재구성 기울기가 정확히 0이 되고, 경계 수심이 지형과 자동으로 정합한다. 여기에 소스항을 $-g\,h_i\,(b_{i+1/2}-b_{i-1/2})/\Delta x$처럼 같은 경계 지형값으로 이산화하면 플럭스 차와 상쇄된다. “물리적으로 상수인 양을 재구성하라”는 한 문장으로 요약되는 이 원리는 유한체적법에서 반복적으로 나타난다 — 상수여야 하는 것을 상수로 다루면 대부분의 균형 문제가 사라진다.

SGM의 약점은 마른 셀이다. 셀의 일부만 잠긴 상태에서 $\eta$를 그대로 재구성하면 $h^{\pm}<0$이 나오므로, 부분 침수 셀에서 재구성을 국소적으로 손보는 보정이 반드시 따라붙는다.

6. 리만 문제에 지형을 넣기[편집]

셋째 노선은 아예 지형을 방정식의 일부로 승격시키는 것이다. $b(x)$를 시간에 대해 불변인 추가 변수로 놓으면 $b_t = 0$이라는 방정식이 하나 늘고, 계는 $(h, hu, b)$의 확장 시스템이 된다. 이 계에는 속도 0인 정지파(stationary wave)가 하나 생기는데, 지형 점프가 그 정지파를 가로지르는 불연속으로 표현되고 소스항이 그 점프 조건 안으로 흡수된다. 정지 평형은 “정지파만 있고 나머지 파의 진폭이 0인 상태”가 되어, 리만 솔버가 그것을 정확히 재현하면 웰밸런스드가 저절로 따라온다.

이 계열의 출발점은 르베크(1998)의 준정상 파동전파 방법으로, 셀 중앙에 지형 점프를 갖는 가상의 리만 문제를 하나 더 놓아 소스를 파 분해 안에 넣는다. 이후 확장(증강) 로 솔버, 경로 보존(path-conservative) 형식 등으로 발전했고, GeoClaw 같은 지구물리 범람 코드가 이 노선을 쓴다.

대가는 이론적 함정이다. 확장 계의 고유값은 $u-\sqrt{gh},\ 0,\ u+\sqrt{gh}$인데, $|u| \to \sqrt{gh}$ 즉 프루드 수 1 근처에서 정지파가 비선형파와 겹친다. 계가 엄밀 쌍곡형이 아니게 되는 이 공명(resonance) 상황에서는 리만 문제의 해가 유일하지 않을 수도 있다. 임계 흐름 근처의 수문·둑 문제에서 도식마다 답이 갈리는 이유가 여기 있고, 여전히 완전히 정리되지 않은 영역이다.²

7. 마름-젖음과 양수성의 상호작용[편집]

범람·처오름 계산에서는 $h \to 0$인 셀이 계속 생기고, 여기서 웰밸런스드와 양수성 보존이 서로 발목을 잡는다.

• 웰밸런스드를 위해 $\eta$를 재구성하면 부분 침수 셀에서 $h<0$이 나올 수 있다. 그대로 두면 $\sqrt{gh}$에서 NaN.
• 반대로 양수성만 지키려고 마른 셀의 재구성 기울기를 0으로 눌러 버리면, 그 셀 주변에서 균형이 깨져 물가 근처에서 가짜 유속이 튄다.
• 속도 복원 $u = (hu)/h$는 $h$가 작을수록 폭발한다. 임계 수심 이하에서 속도를 죽이는 비특이화(desingularization)를 쓰는데, 이 조작 자체가 운동량을 미세하게 바꿔 균형을 흔든다.

그래서 실무 도식은 (1) $\max(0,\cdot)$가 들어간 정수압 재구성, (2) 부분 침수 셀 전용 재구성 보정, (3) 양수성을 보장하는 $\Delta t$ 제한, (4) 비특이화 임계값을 한 세트로 설계한다. 하나만 빼도 물가에서 터지거나 조용히 틀린다. 마찰항이 마른 셀에서 강성이 되는 문제까지 겹치므로 반음해 처리도 대개 같이 들어간다. 자세한 물리적 배경은 개수로 유동과 얕은 물 방정식 참고.

8. 움직이는 평형, 그리고 다른 물리로의 일반화[편집]

정지수만이 평형은 아니다. 1차원 정상 유동에서는 유량과 에너지가 보존된다.

$$hu = \text{const}, \qquad \frac{u^2}{2} + g(h+b) = \text{const}$$

이 움직이는 평형(moving equilibrium)을 정확히 보존하는 도식은 훨씬 어렵다. 평형 조건이 비선형 대수식이라 각 셀에서 그것을 풀어야 하고, 아임계·초임계 두 해 중 어느 쪽인지 골라야 하며, 임계점 근처에서 그 선택이 불안정해진다. 뇔레·싱·슈(2007) 이후 고차 웰밸런스드 도식들이 나와 있지만 정지수만큼 표준화되지는 않았다. 마찰이 있는 수로에서 바닥 경사와 마찰이 균형을 이룬 정상수심 상태를 보존하는 문제도 같은 계열이다.

발상 자체는 얕은 물에 국한되지 않는다. 소스항이 있는 쌍곡계에서 물리적으로 중요한 평형이 존재한다면 같은 프로그램이 반복된다.

• 중력 하의 오일러 방정식. 정수압 평형 $\partial p/\partial z = -\rho g$를 이산 수준에서 보존하지 못하면, 별 내부나 대기 모델에서 정지 대기가 스스로 대류를 시작한다. 실제 대류 섭동이 배경 압력의 $10^{-6}$ 수준인 항성 대류 계산에서는 이게 결정적이라, 국소 정수압 재구성 계열의 웰밸런스드 도식이 표준이 됐다.
• 반응·이완 소스. 화학 평형이나 완화 항이 있는 계에서 평형 상태를 정확히 보존하지 못하면, 아무 일도 없어야 할 영역에서 가짜 반응이 진행된다.
• 회전·자기장 평형 등도 같은 틀에서 다뤄진다.

9. 검증하는 법[편집]

웰밸런스드 여부는 격자 수렴 시험으로 드러나지 않으므로, 전용 시험을 돌린다.

1. 정지수 시험. 불규칙한 지형(젖은 셀만 있는 경우와 마른 셀이 섞인 경우 각각) 위에 수평 수면을 놓고 오래 돌린다. $\max|u|$와 $\|h+b-\eta_0\|_\infty$가 배정밀도 기준 $10^{-14}$ 근처에 머물면 exact C-property를 만족하는 것이고, $10^{-3}$쯤에서 정체하면 아니다. 지형 경사를 키워 가며 반복하는 것이 핵심으로, 평평한 바닥에서는 안 웰밸런스드 도식도 통과한다.
2. 미소 섭동 시험. 위 정지 상태에 수면 섭동 $10^{-3}$을 얹어 지형(예: 매끈한 융기) 위를 지나가게 한다. 웰밸런스드 도식은 섭동이 깨끗한 좌우 파로 갈라져 지나가고, 아닌 도식에서는 지형이 만든 가짜 신호에 섭동이 파묻힌다. 두 결과를 나란히 그려 놓은 그림이 이 분야 논문의 국룰 삽화다.
3. 장시간 보존. 수천~수만 스텝을 돌려 오차가 선형으로 자라는지 본다. 근사 C-property 도식은 여기서 걸린다.

정리하면 웰밸런스드는 정확도의 문제가 아니라 구조 보존(structure-preserving)의 문제다. 고들리노프 도식 계열이 보존성을 구조적으로 보장하듯, 웰밸런스드는 평형을 구조적으로 보장한다. 격자를 아무리 조밀하게 해도 살 수 없는 성질을, 도식을 고쳐서 공짜로 얻는 몇 안 되는 사례다.³

10. 관련 문서[편집]

• 얕은 물 방정식 · 개수로 유동 · 프루드 수 · 수충격
• 유한체적법 · 리만 솔버 · 차분 도식
• 수치 소산과 분산 · CFL 조건 · 적응 격자 세분화
• 특성곡선법 · 충격파 · 전산유체역학
• 검증 및 확인 · 원시방정식
• 양수성 보존 · 고들리노프 도식 · MUSCL

11. Footnotes[편집]

1. 이 증상을 처음 보면 대개 초기조건이나 경계조건을 의심하며 며칠을 태운다. “물을 가만히 놔뒀는데 왜 움직이죠”라는 질문에 “도식이 그렇게 생겼습니다”라는 답이 돌아오는 순간, 수치해석이 물리와 별개의 학문이라는 사실이 몸으로 이해된다. ↩

2. 공명 조건이 얄미운 이유는 하필 실무에서 중요한 지점에서 발생한다는 것이다. 수문 아래, 둑 위, 여울목 — 임계 흐름은 하천 공학에서 제일 관심 있는 단면들이다. 이론이 제일 흔들리는 곳과 답이 제일 궁금한 곳이 정확히 일치하는, 흔한 종류의 저주. ↩

3. 그래서 논문의 셀 수 대신 정지수 시험의 $\max|u|$가 $10^{-15}$인지 $10^{-4}$인지를 먼저 확인하는 습관이 생긴다. 격자 수렴 그래프는 웬만하면 예쁘게 나오고, 예쁘게 나오는 지표는 정보량이 적다. ↩