수정 방정식

수치해석 전산유체역학 마지막 수정: 2026-07-20 04:21:04

1. 개요[편집]

수정 방정식 Modified Equation
다른 이름	등가 미분방정식(equivalent differential equation)
분석 대상	차분 스킴이 실제로 푸는 PDE
주요 도구	테일러 급수 전개
알려 주는 것	수치 소산 · 수치 분산 · 정밀도 차수
제안	Warming & Hyett (1974)

수정 방정식(modified equation)은 어떤 차분 스킴이 겉으로 이산화하려던 원래 편미분방정식이 아니라 실제로는 어떤 PDE를 정확히 풀고 있는가를 밝혀내는 분석 기법이다. 원 방정식에 절단 오차(truncation error) 항들이 미분항의 형태로 더 붙은 등가 PDE를 유도하여, 스킴이 만드는 오차의 물리적 성격을 대수적으로 읽어낸다. 1974년 워밍(R. F. Warming)과 하이엇(B. J. Hyett)이 체계화했다.¹

발상은 반전 매력이 있다. 우리는 예컨대 대류-확산 방정식을 풀려고 차분식을 세웠는데, 막상 컴퓨터가 밟고 있는 진짜 방정식은 원 방정식 + 여분의 미분항들이다. 이 여분 항이 짝수 차 미분(2계, 4계…)이면 확산처럼 해를 뭉개고(수치 소산), 홀수 차 미분(3계, 5계…)이면 파를 어긋나게 흘려보낸다(수치 분산). 즉 수정 방정식은 수치 소산과 분산이 “어디서 오는가”에 대한 가장 직접적인 답을 준다.

2. 유도 절차[편집]

수정 방정식을 얻는 표준 절차는 테일러 급수 전개 한 판이다. 1차 상류 차분(upwind)으로 이산화한 선형 대류 방정식 $u_t + a u_x = 0$ ( $a>0$ )을 예로 들자. 차분식은

\frac{u_i^{n+1} - u_i^{n}}{\Delta t} + a\,\frac{u_i^{n} - u_{i-1}^{n}}{\Delta x} = 0

이다. 각 격자점 값을 $u(x_i, t^n)$ 주위로 테일러 전개해 대입하고 정리하면, 이산 관계가 아니라 연속 함수 $u$ 가 만족하는 방정식으로 다음이 나온다.

u_t + a u_x = \frac{a\,\Delta x}{2}\left(1 - \nu\right) u_{xx} + O(\Delta x^2), \qquad \nu = \frac{a\,\Delta t}{\Delta x}

우변 첫 항이 바로 수치 점성(numerical viscosity)이다. 여기 등장하는 무차원수 $\nu$ 가 CFL 조건의 그 쿠랑수(Courant number)라는 점이 결정적이다.

한 가지 실무 팁: 유도 과정에서 시간 미분( $u_{tt}$ 등)은 원 방정식 관계 $u_t = -a u_x$ 를 반복 사용해 공간 미분으로 치환해야 한다. 이 치환을 빼먹으면 계수가 틀려 엉뚱한 결론에 도달한다.²

3. 소산·분산과 정밀도 차수 읽기[편집]

수정 방정식의 우변 항을 미분 차수로 분류하면 스킴의 성격이 한눈에 보인다.

짝수 차 미분 항( $u_{xx}, u_{xxxx}, \dots$ ): 소산(dissipation). 계수가 양이면 확산처럼 진폭을 감쇠시켜 수치적으로 안정화하지만 해를 뭉갠다. 계수가 음이면 반확산(anti-diffusion)이라 폭발한다.
홀수 차 미분 항( $u_{xxx}, u_{xxxxx}, \dots$ ): 분산(dispersion). 진폭은 살리되 파수(wavenumber)마다 위상 속도를 다르게 만들어, 급경사 뒤에 **비물리적 진동(wiggle)**을 남긴다.

우변에서 살아남는 최저 차 항의 $\Delta x$ 멱수가 곧 그 스킴의 공간 정밀도 차수다. 위 상류 차분은 우변 최저항이 $\Delta x^1$ 이므로 1차 정밀도이고, 그 항이 하필 소산항이라 “안정하지만 뭉갠다”는 상류법의 악명이 대수적으로 설명된다. 차분 도식 문서의 QUICK·TVD 논의가 노리는 것도 결국 이 소산·분산 균형을 손보는 일이다.

4. 관련 안정성 분석과의 관계[편집]

수정 방정식은 폰 노이만 안정성 해석, CFL 조건과 형제 관계지만 각자 보는 각도가 다르다. 셋을 헷갈리지 않는 것이 중요하다.

기법	도구	답하는 질문
폰 노이만 안정성 해석	푸리에 모드 증폭계수	오차가 시간에 따라 터지는가
CFL 조건	의존 영역 비교	안정에 필요한 $\Delta t$ 상한은
수정 방정식	테일러 급수 전개	실제로 푸는 PDE는 무엇이며 오차의 성격은

세 가지는 서로를 뒷받침한다. 폰 노이만 해석은 증폭계수의 크기로 안정성을, 위상으로 분산을 정량화하는데, 그 결과가 수정 방정식이 예측하는 소산·분산 항과 정확히 일치한다.³ 예컨대 수정 방정식의 소산 계수 $\frac{a\Delta x}{2}(1-\nu)$ 가 음이 아니려면 $\nu \le 1$ 이어야 하는데, 이 조건이 바로 상류 스킴의 CFL 안정 한계와 같다. 세 렌즈가 같은 진실의 다른 면인 셈이다. 요컨대 폰 노이만이 “터지느냐”를, CFL이 “얼마나 작게 밟느냐”를, 수정 방정식이 “무엇을 밟고 있느냐”를 답한다.

5. 활용과 한계[편집]

수정 방정식이 실무에서 빛나는 지점은 다음과 같다.

스킴 진단: 전산유체역학 결과에서 충격파가 번지면 소산 항을, 뒤에 진동이 끼면 분산 항을 지목해 원인을 대수적으로 짚는다.
인공 점성 설계: 필요한 만큼만 소산을 의도적으로 주입할 때, 수정 방정식이 그 양의 가이드가 된다.
고차 스킴 설계: 저차 항을 상쇄하도록 계수를 맞춰 정밀도를 끌어올리는 Lax–Wendroff류 설계의 이론적 뼈대.

다만 한계도 분명하다. 유도 자체가 테일러 급수에 기대므로 매끄러운 해를 가정한다. 충격파처럼 불연속이 있으면 전개가 무너져 정량적 예측력이 떨어진다. 또 비선형 방정식에서는 항이 급격히 지저분해져 손계산이 사실상 자동 미분/기호 연산의 영역으로 넘어간다.⁴ 그럼에도 “내 스킴이 진짜로 뭘 풀고 있나”를 묻는 순간, 수정 방정식만큼 정직하게 답해 주는 도구는 없다.

6. 관련 문서[편집]

Warming, R. F. & Hyett, B. J. (1974). “The modified equation approach to the stability and accuracy analysis of finite-difference methods.” J. Comput. Phys. 이들은 수정 방정식의 우변 부호만으로 안정성까지 판정하는 규칙을 제시했다. 폰 노이만 해석과 놀랍도록 잘 맞는다. ↩
이 시간 미분 → 공간 미분 치환을 빼먹는 것이 대학원생 수정 방정식 유도의 국룰 실수다. 원 PDE를 시간으로 한 번 더 미분해 $u_{tt} = a^2 u_{xx}$ 를 얻는 식으로 소거해야 한다. 안 하면 계수에 엉뚱한 부호가 붙어 “안정한 스킴이 불안정해 보이는” 괴담이 탄생한다. ↩
다르게 말하면 폰 노이만 해석의 증폭계수 $G(k)$ 를 작은 파수에서 테일러 전개한 것이 수정 방정식이다. $|G|$ 의 저차항이 소산, $\arg G$ 의 저차항이 분산에 대응한다. 두 기법이 같은 정보를 주파수 공간과 물리 공간에서 각각 보는 것. ↩
손으로 3계 항까지 전개하다 보면 부호 하나 놓쳐 밤을 새우게 된다. 요즘은 SymPy 같은 기호 연산으로 돌리는 게 인지상정. 그래도 개념을 모르면 나온 항이 맞는지 검산할 수가 없다. ↩