특이 섭동

1. 개요[편집]

특이 섭동 Singular Perturbation
징후	ε → 0 에서 최고차 미분이 사라진다
결과	경계조건 하나를 못 맞춘다 · 경계층 발생
도구	정합 점근 전개, 반 다이크 정합 원리
친척	WKB, 다중척도법, 재규격화군
대표 응용	프란틀 경계층, QSSA, 강성 ODE
수치적 함의	층을 격자로 못 잡으면 진동한다

작은 수를 0으로 보내면 되지 않나? 그 작은 수가 최고차 미분 앞에 붙어 있다면, 안 된다.

특이 섭동(singular perturbation)은 작은 매개변수 $\varepsilon$ 을 0 으로 보내는 순간 문제의 성격 자체가 바뀌어 버려서, 해를 $\varepsilon$ 의 멱급수로 매끄럽게 전개하는 정칙 섭동이 실패하는 상황을 말한다. 미분방정식에서는 대개 $\varepsilon$ 이 최고차 도함수에 곱해져 있어, $\varepsilon = 0$ 을 넣으면 방정식의 차수가 떨어지고 경계조건 하나를 감당할 수 없게 되는 형태로 나타난다.

물리적으로 이것은 “어디선가 무시했던 항이 좁은 영역에서 되살아난다”는 뜻이다. 그 좁은 영역이 경계층(boundary layer) 혹은 내부층(interior layer)이고, 특이 섭동 이론은 바깥 영역과 층 안쪽에서 각각 다른 근사를 만들어 둘을 이어 붙이는 체계적인 절차를 제공한다. 경계층 이론, 화학의 준정상상태 근사, 강성 방정식의 빠른 모드, 반도체 소자의 공핍층, 얇은 판의 경계 조건 처리가 전부 같은 뼈대를 공유한다.

2. 정칙 섭동이 왜 실패하는가[편집]

먼저 잘 되는 경우를 보자. $x^2 + \varepsilon x - 1 = 0$ 의 근을 $x = x_0 + \varepsilon x_1 + \cdots$ 로 전개하면 $x_0 = \pm 1$ 에서 시작해 항을 차례로 얻는다. 근의 개수가 유지되고 각 근이 $\varepsilon \to 0$ 에서 얌전히 극한을 갖는다. 이것이 정칙 섭동(regular perturbation)이다.

이번엔 $\varepsilon x^2 + x - 1 = 0$. $\varepsilon = 0$ 을 넣으면 $x = 1$ 하나만 남는다. 근이 하나 증발했다. 실제로는 $x \approx -1/\varepsilon$ 이라는 두 번째 근이 $\varepsilon \to 0$ 에서 무한대로 달아난 것이다. 이 근을 붙잡으려면 $x = X/\varepsilon$ 처럼 변수를 다시 재서 균형을 되찾아야 한다. 특이 섭동의 모든 기술은 결국 이 “다시 재기”의 변주다.

미분방정식에서는 사라지는 것이 근이 아니라 자유도다. 2계 방정식에서 $\varepsilon y''$ 를 버리면 1계 방정식이 되고, 1계 방정식은 경계조건을 하나밖에 못 만족시킨다. 그런데 원래 문제는 조건 두 개를 요구한다. 버림받은 경계조건이 있는 자리에 층이 생긴다.¹

3. 표준 예제[편집]

교과서가 늘 같은 문제로 시작하는 데는 이유가 있다.

$$\varepsilon y'' + y' + y = 0, \qquad y(0) = 0, \quad y(1) = 1, \qquad 0 < \varepsilon \ll 1$$

특성근을 보면 무슨 일이 벌어지는지 즉시 알 수 있다. $\varepsilon m^2 + m + 1 = 0$ 의 두 근은

$$m_1 = -1 - \varepsilon + O(\varepsilon^2), \qquad m_2 = -\frac{1}{\varepsilon} + 1 + O(\varepsilon)$$

이다. 하나는 $O(1)$ 의 느린 감쇠, 다른 하나는 $O(1/\varepsilon)$ 의 폭발적으로 빠른 감쇠. 빠른 모드 $e^{-x/\varepsilon}$ 는 $x = 0$ 근처 폭 $O(\varepsilon)$ 안에서만 살아 있고 그 밖에서는 지수적으로 죽는다. 층은 $x = 0$ 에 있다.

3.1. 외부해[편집]

$\varepsilon y''$ 를 버리면 $y' + y = 0$, 즉 $y_{\mathrm{out}} = A e^{-x}$. 이 1계 해가 두 경계조건을 다 만족할 수는 없으니 하나를 골라야 한다. 층이 $x=0$ 에 있으므로 층 바깥의 조건인 $y(1) = 1$ 을 쓴다.

$$y_{\mathrm{out}}(x) = e^{\,1-x}$$

$x = 0$ 에서 이 값은 $e \approx 2.718$ 이라 요구된 0 과 전혀 맞지 않는다. 그 불일치를 메우는 것이 층의 임무다.

3.2. 내부해[편집]

층 안에서는 $y'$ 가 $O(1/\varepsilon)$ 이므로 $x$ 를 그대로 쓰면 안 된다. 늘린 좌표 $X = x/\varepsilon$ 을 도입하면 $d/dx = \varepsilon^{-1} d/dX$ 이고, 방정식은

$$\frac{1}{\varepsilon}\left( \frac{d^2 Y}{dX^2} + \frac{dY}{dX} \right) + Y = 0 \;\;\Longrightarrow\;\; \frac{d^2 Y}{dX^2} + \frac{dY}{dX} + \varepsilon Y = 0$$

가 된다. 최고차 도함수가 살아남았다. 이것이 늘리기의 목적이다. 선도 차수에서 $Y'' + Y' = 0$ 이므로 $Y = B + C e^{-X}$ 이고, 내부 경계조건 $Y(0) = 0$ 에서 $C = -B$,

$$Y(X) = B\left( 1 - e^{-X} \right)$$

늘리는 배율을 $x/\varepsilon$ 로 고른 것이 우연이 아니라는 점은 짚고 넘어갈 만하다. $x = \delta(\varepsilon) X$ 로 두고 세 항의 크기 $\delta^{-2}\varepsilon$, $\delta^{-1}$, $1$ 을 비교하면, 앞의 둘이 같은 크기가 되는 $\delta = \varepsilon$ 에서만 의미 있는 균형이 생긴다. 이런 선택을 구별 극한(distinguished limit)이라 하며, 어떤 문제든 층 두께를 결정하는 것은 이 균형 계산이다.

3.3. 정합과 합성해[편집]

남은 미지수 $B$ 는 정합(matching)으로 정한다. 층에서 멀어질 때의 내부해와, 층에 다가갈 때의 외부해가 겹치는 중간 영역에서 일치해야 한다.

$$\lim_{X \to \infty} Y(X) = B, \qquad \lim_{x \to 0} y_{\mathrm{out}}(x) = e \;\;\Longrightarrow\;\; B = e$$

두 해를 더하고 중복해서 센 공통부분($=e$)을 빼면 전 구간에서 통하는 합성해(composite solution)를 얻는다.

$$y_{\mathrm{comp}}(x) = e^{\,1-x} + e\left(1 - e^{-x/\varepsilon}\right) - e = e^{\,1-x} - e\, e^{-x/\varepsilon}$$

확인해 보면 $y(0) = e - e = 0$, $y(1) = 1 - e\,e^{-1/\varepsilon} \approx 1$ 로 두 경계조건을 다 만족한다. 층 두께가 $O(\varepsilon)$ 이라는 것은 $e^{-x/\varepsilon}$ 이 $x \sim \varepsilon$ 에서 대부분 죽는다는 뜻이고, $\varepsilon = 10^{-3}$ 이면 해가 $x \in [0, 0.003]$ 구간에서 0 에서 2.7 까지 치솟는다는 얘기다.

일반적인 정합 규칙은 반 다이크 정합 원리로 정리된다. 요약하면 “$m$항 내부 전개의 $n$항 외부 전개는 $n$항 외부 전개의 $m$항 내부 전개와 같다”. 실무에서는 이걸 기계적으로 적용해 상수를 결정하며, 고차항으로 갈수록 로그항($\ln \varepsilon$)이 끼어들어 순진한 멱급수로는 정합이 안 되는 경우가 나온다 — 저 레이놀즈수 유동의 오센 문제가 유명한 사례다.

4. 프란틀 경계층이 정확히 같은 구조다[편집]

1904년 프란틀이 한 일이 이 절차 그 자체다. 고 레이놀즈수 유동에서 점성항은 작아 보이므로 버리면 오일러 방정식이 되는데, 오일러 방정식은 2계 항을 잃어 점착 조건을 만족할 수 없다. 그 결과가 달랑베르의 역설(항력이 0)이다.

프란틀은 벽 근처에서 $y$ 를 $\delta$ 로 늘려 점성항을 되살렸다. 나비에-스토크스 방정식에서 대류 $u\,\partial u/\partial x \sim U^2/L$ 과 점성 $\nu\, \partial^2 u/\partial y^2 \sim \nu U/\delta^2$ 이 같은 크기가 되는 구별 극한을 잡으면

$$\delta \sim \sqrt{\frac{\nu L}{U}} = \frac{L}{\sqrt{\mathrm{Re}}}$$

즉 층 두께가 $O(\varepsilon^{1/2})$ 이다($\varepsilon = 1/\mathrm{Re}$). 앞의 ODE 예제는 $O(\varepsilon)$ 이었는데 여기서는 $O(\sqrt{\varepsilon})$ — 층 두께의 지수는 보편 상수가 아니라 구별 극한이 결정한다는 점이 여기서 드러난다. 외부해가 포텐셜 유동, 내부해가 경계층 방정식, 정합이 “경계층 바깥 속도 = 포텐셜 해의 벽면 속도”라는 경계조건, 배제두께를 이용한 반복이 고차 정합에 해당한다. 프란틀은 점근 해석의 언어 없이 물리적 직관만으로 이 구조를 발견했고, 그것이 체계적인 특이 섭동 이론으로 정리된 것은 반 다이크와 콜의 1960년대다.

층류 화염 속도의 ZFK 해석도 같은 골격이다. 거기서는 작은 매개변수가 젤도비치 수의 역수 $1/\beta$ 이고, 예열층이 외부 영역, 얇은 반응층이 내부층이며, 정합으로 결정되는 미지수가 질량유속(=화염 속도)이다.

5. 화학의 QSSA와 강성 ODE[편집]

시간 방향으로 같은 일이 벌어지면 이름이 바뀐다. 라디칼처럼 매우 빠르게 생성·소멸하는 화학종은 지배 방정식이

$$\varepsilon \frac{dR}{dt} = f(R, S), \qquad \frac{dS}{dt} = g(R, S)$$

꼴이 되고, $\varepsilon \to 0$ 이면 미분방정식이 대수방정식 $f(R,S) = 0$ 으로 붕괴한다. 이것이 준정상상태 근사(QSSA)다. 여기서도 초기조건 $R(0)$ 을 만족시킬 수 없는데, 그 자리에 시간 방향의 층 — 초기층(initial layer) — 이 생긴다. 라디칼 농도가 초기값에서 준정상 값으로 급히 달려가는 그 짧은 구간이다. $f_R < 0$ 이라는 안정성 조건 아래 이 붕괴가 정당하다는 것을 보장하는 것이 티호노프 정리이며, 사실상 관성 다양체·느린 다양체 이론의 화학판이다.

이 관점을 자동화한 것이 CSP(computational singular perturbation)와 ILDM 계열이고, 메커니즘 축소의 이론적 뿌리가 여기 있다. “어떤 종을 QSSA 로 놓을까”를 사람이 감으로 고르던 것을, 자코비안의 고윳값·고유벡터를 계산해 빠른 부공간을 식별하는 문제로 바꾼 것이다.

수치 쪽에서 이것을 부르는 이름이 강성이다. 자코비안의 고윳값 스펙트럼이 $-1/\varepsilon$ 과 $O(1)$ 로 갈라져 있으면 강성비가 $1/\varepsilon$ 이고, 명시적 적분기는 안정성 때문에 $\Delta t \lesssim \varepsilon$ 에 묶인다. 관심 있는 현상은 $O(1)$ 시간에 벌어지는데 스텝은 $\varepsilon$ 으로 잘라야 하는 이 부조리가 곧 강성 방정식 문제이고, 암시적 방법(BDF 류)은 초기층을 해상하지 않고 건너뛰어 느린 다양체 위로 곧장 착지하는 방식으로 이를 회피한다. 특이 섭동은 강성의 이론적 얼굴, 강성은 특이 섭동의 계산적 얼굴이라고 보면 대체로 맞다.

6. WKB와 다중척도법[편집]

정합 점근 전개는 층이 국소적일 때 쓰는 도구다. 작은 매개변수의 효과가 정의역 전체에 퍼져 누적되면 다른 무기가 필요하다.

• WKB 근사. $\varepsilon^2 y'' = Q(x)\,y$ 처럼 진동·감쇠가 전 구간에 걸친 문제에서는 해를 $y \sim \exp(S(x)/\varepsilon)$ 로 놓고 $S$ 를 전개한다. 위상이 $1/\varepsilon$ 배로 빨리 도는 것을 지수 안에 미리 넣어 두는 셈이다. $Q(x) = 0$ 인 전환점(turning point)에서 근사가 깨지고, 그 주위에 에어리 함수로 기술되는 국소층이 생겨 양쪽 WKB 해를 정합한다 — 층과 정합이라는 구조가 여기서도 반복된다. 반고전 양자역학의 표준 도구이기도 하다.
• 다중척도법. 약하게 감쇠하는 진동자 $y'' + y + \varepsilon y' = 0$ 를 순진하게 전개하면 $\varepsilon t$ 에 비례하는 세속항(secular term)이 나와 $t \sim 1/\varepsilon$ 에서 근사가 무너진다. $t$ 와 $\tau = \varepsilon t$ 를 독립 변수로 취급하고 세속항이 나타나지 않을 조건을 요구하면, 진폭이 $\tau$ 에 대해 어떻게 천천히 변하는지가 자동으로 결정된다. 비선형 진동·분기 이론의 진폭 방정식이 대개 이 방식으로 유도된다.

셋(정합 전개, WKB, 다중척도)은 별개의 기법처럼 보이지만 “문제에 두 개 이상의 척도가 있고, 순진한 전개가 그 중 하나를 놓친다”는 같은 진단에서 출발한다.

7. 수치적 함의 — 층을 못 잡으면 진동한다[편집]

가장 실용적인 결말이 여기 있다. 앞의 예제에서 $y'$ 를 이류항, $\varepsilon y''$ 를 확산항으로 읽으면 그대로 1차원 **대류-확산 방정식**이다. 균일 격자 $h$ 에 중심차분을 쓰면 셀 페클레 수

$$\mathrm{Pe}_h = \frac{u h}{\varepsilon}$$

가 2 를 넘는 순간 이산해가 격자 단위로 진동한다. 물리적으로는 층 두께 $O(\varepsilon)$ 을 격자가 담지 못했다는 뜻이고, 대수적으로는 이산 방정식의 특성근이 음수가 되어 부호가 교대하기 때문이다. 즉 진동은 버그가 아니라 “층을 해상하지 못했다”는 신호다.

대응은 크게 셋이다.

• 풍상 차분. 1차 상류화는 무조건 단조성을 보장하지만, 실질적으로 확산계수를 $\varepsilon \to \varepsilon + uh/2$ 로 키운 것과 같아 층을 두껍게 번지게 만든다. 진동은 사라지는데 답도 같이 사라지는 셈이라, 수치분산과 함께 CFD 입문자가 처음 배우는 타협이다.
• 지수 적합. 국소적으로 층의 지수 구조를 이산화에 심는다. 1차원 정상 대류-확산이라면 일린-앨런-사우스웰 형태의 스킴이 격자점에서 정확한 해를 준다. 파탄카의 멱법칙 스킴이 이 지수 함수를 값싸게 흉내 낸 것이고, 유한요소 쪽의 SUPG 안정화도 시험함수를 상류로 기울여 같은 효과를 노린다.
• 층에 맞춘 격자. 전이점을 $O(\varepsilon \ln N)$ 에 두는 시시킨 격자 같은 조각별 균일 격자를 쓰면 $\varepsilon$ 에 무관하게 수렴하는(ε-uniform) 오차 상계를 증명할 수 있다. 실무 CFD 의 벽면 프리즘 레이어가 이론적 정당화 없이 하고 있던 일이 사실 이것이다 — 층 두께를 미리 추정해 그 안에 격자를 몇 겹 밀어 넣는 것.

요약하면, 층이 어디에 얼마나 얇게 생기는지를 아는 것이 곧 격자 설계다. 특이 섭동 해석은 계산을 시작하기 전에 그 답을 손으로 알려 준다는 점에서, 순수 이론이라기보다 격자를 깔기 위한 사전 조사에 가깝다.²³

8. 관련 문서[편집]

• 경계층 · 나비에-스토크스 방정식 · 포텐셜 유동
• 강성 방정식 · 관성 다양체 · 메커니즘 축소
• 대류-확산 방정식 · 페클레 수 · 수치분산 · 풍상 차분
• 점근 전개 · 다중척도법 · 분기 이론 · 동역학계
• 층류 화염 속도 · 반응속도론 · 연산자 분리

9. Footnotes[편집]

1. 특이 섭동과 정칙 섭동을 가르는 실전 진단법: $\varepsilon = 0$ 을 대입해 보고 (1) 방정식의 차수가 떨어지거나, (2) 정의역이 무한해지거나, (3) 해가 $t \to \infty$ 에서 균일하지 않으면 특이하다. 셋 중 아무것에도 안 걸리면 마음 놓고 급수를 전개해도 된다. 대개는 걸린다. ↩

2. “격자를 두 배로 촘촘하게 했더니 답이 완전히 달라졌어요”라는 상담의 상당수가 이 문제다. 층을 못 잡던 격자에서 층을 잡기 시작한 격자로 넘어간 것이라, 두 결과가 다른 게 정상이다. 격자 수렴 지수를 계산할 때 이 구간을 걸치고 있으면 수렴 차수가 이론값과 전혀 안 맞는다. ↩

3. 상류화가 층을 번지게 한다는 사실을 두고 “그래서 1차 상류차분은 쓰면 안 된다”고 정리하는 것은 성급하다. 층을 애초에 해상할 예산이 없다면, 진동해서 발산하는 답보다는 번져도 수렴하는 답이 낫다. 문제는 그 답을 정답이라고 보고하는 순간이다. ↩