특이 섭동과 경계층의 정합

1. 작은 계수가 경계조건 하나를 책임질 수 있다

양 끝을 고정한 곡선을 생각해 봅시다. 대부분의 구간에서는 완만하지만, 한쪽 끝의 고정값에 맞추려면 아주 좁은 곳에서 급하게 휘어야 할 수 있습니다. 이때 작은 매끄럽게 하는 항을 지워 버리면 급한 변화가 사라집니다. 계수가 작다는 사실과 그 계수에 미분값을 곱한 항이 작다는 사실은 다릅니다.

이번 장에서는 정확히 풀리는 경계값 문제로 이 현상을 확인합니다. 한 점에서의 정확도, 전체 구간의 정확도, 미분값의 정확도, 미분방정식 잔차의 크기를 구분하는 것이 목표입니다. 같은 그림이 그럴듯해 보여도 네 질문의 답은 다를 수 있습니다.

앞선 두 장의 비선형 평형과 유한 적분에서는 차수별 근사를 비교했습니다. 여기서는 공간의 일부에서만 전개가 실패하는 특이 극한을 다룹니다.

2. 방정식·부호·경계값부터 고정한다

무차원화한 문제는 다음과 같습니다.

$$\varepsilon u''(x)+a u'(x)=f,\qquad 0<x<1,\quad u(0)=A,\quad u(1)=B,\qquad \varepsilon>0.$$

a와 f는 상수이며 프라임은 x에 대한 미분입니다. 실험 입력은 유리수로 받습니다. 0<ε≤16, |a|·|f|·|A|·|B|≤16, 관측 위치는 [0,1], 등간격 격자 패널은 4–256개입니다. a=0은 별도 식으로 계산합니다. ε=0은 극한 방정식이며 실제 계산에 넣는 값이 아닙니다.

수송 방정식의 부호를 주의해야 합니다. 양변에 −1을 곱하면 흔히 쓰는 −εu″+v u′=s 꼴이고 v=−a, s=−f입니다. 이 장의 규약에서는 a>0일 때 감쇠하는 경계층이 왼쪽에 생깁니다. 다른 자료에서 기억한 흐름 방향을 부호 확인 없이 가져오지 마세요.

3. ε를 0으로 두면 무엇을 잃는가

a≠0에서 최고차 미분항을 지우면 a u′=f라는 일차 방정식이 됩니다. 적분상수는 하나인데 원래 문제는 양 끝값 두 개를 요구합니다. 두 경계값이 우연히 그 기울기와 맞지 않으면 모두 만족시킬 수 없습니다.

a>0일 때 먼쪽 경계값을 만족시키는 바깥해는

$$u_{\mathrm{out}}(x)=B+\frac{f}{a}(x-1)$$

입니다. 왼쪽 값은 B−f/a입니다. D=B−A−f/a라고 쓰면 왼쪽 경계 오차는 D입니다. D가 ε와 무관한 0이 아닌 수라면, ε를 아무리 줄여도 이 끝점 오차는 줄지 않습니다. “ε가 매우 작다”는 말만으로 경계조건을 무시할 수 없는 이유입니다.

4. 흐름 부호가 달라도 같은 좌표로 비교한다

b=|a|로 둡니다. a>0이면 y=x, 가까운 경계값 N=A, 먼 경계값 F=B로 정합니다. a<0이면 y=1−x, N=B, F=A로 반사합니다. 일차 미분의 부호가 바뀌므로 변환한 문제의 계수는 항상 양수입니다.

$$\varepsilon v_{yy}+b v_y=f,\qquad v(0)=N,\quad v(1)=F.$$

k=f/b, D=F−N−k, δ=ε/b라고 쓰겠습니다. y는 경계층이 있는 경계에서부터의 거리이지 언제나 x=0으로부터의 거리는 아닙니다. 반사하면 물리적 일차 미분에는 −1이 붙고 이차 미분은 그대로입니다. 관측점의 기울기를 비교할 때 이 부호도 유지합니다.

5. 먼저 정확해를 유도한다

w=v_y라고 두면 εw_y+bw=f라는 일차 방정식입니다. 해는 상수 k와 감쇠 지수함수의 합입니다. 한 번 적분하고 양 끝값을 적용하면

$$v(y)=N+ky+D\frac{1-e^{-y/\delta}}{1-e^{-1/\delta}}$$

를 얻습니다. 유한한 δ>0에서 분모는 양수입니다. 이 식은 δ가 작지 않아도 정확합니다. y=0에서는 비율이 정확히 0이고 y=1에서는 정확히 1입니다. 구현도 이 두 점에서는 가까운 부동소수점 수를 빼는 대신 정확한 유리수 경계값을 반환합니다.

기본 예제 ε=1/10, a=1, f=0, A=0, B=1에서는 (1−e^(−10x))/(1−e^(−10))입니다. x=1/20이면 늘린 좌표는 1/2이고 값은 약 0.39349입니다. 바깥해는 상수 1이지만 왼쪽 근처에서는 큰 변화가 남습니다.

다음 표에서 바깥해는 세 점 모두 1입니다. 작은 x에서는 바깥해가 크게 틀리지만 합성해는 처음부터 좁은 오차 안에 있습니다. 반대로 합성해의 최대 오차는 층의 가운데가 아니라 먼 끝점에서 얻어진다는 것도 확인하세요. 표의 소수는 반올림한 설명값이며 정확한 포함 구간은 계산 표에서 읽습니다.

x	정확해	합성해	정확해−합성해
0.05	0.3934872046	0.3934693403	0.0000178643
0.1	0.6321492584	0.6321205588	0.0000286995
0.5	0.9933071491	0.9932620530	0.0000450961

6. 지배적 균형으로 경계층 폭을 찾는다

값이 크기 1 정도 바뀌는 공간 폭을 ℓ이라고 가정합시다. 일차 미분은 1/ℓ, 이차 미분은 1/ℓ² 규모입니다. 두 미분항을 모두 남기려면 ε/ℓ²와 b/ℓ가 같은 규모여야 합니다.

$$\ell\sim\delta=\frac{\varepsilon}{b},\qquad X=\frac{y}{\delta},\qquad v_{XX}+v_X=\delta k.$$

b>0과 나머지 자료를 고정하고 ε→0으로 보내면 δ→0이고 δk도 작아집니다. 선도 안쪽 방정식은 v_XX+v_X=0입니다. 반대로 ε가 작아도 b가 더 작으면 얇은 층이라는 해석이 실패할 수 있습니다. ε=0.1, b=0.001이면 δ=100입니다. 길이 1인 구간에 비해 좁은 경계층이 아닙니다.

7. 안쪽해와 바깥해를 정합한다

정준 좌표의 바깥해는 F+k(y−1)입니다. 가까운 경계로 접근한 값은 F−k입니다. 선도 안쪽해는 가까운 경계값을 만족하고, X가 커질 때 이 공통값으로 가도록 정합니다.

$$v_{\mathrm{in}}(X)=F-k-D e^{-X}.$$

X=0에서 N이고 X→∞에서 F−k입니다. 정합은 선도 안쪽해에 반대편 경계 X=1/δ의 조건까지 억지로 부과하는 일이 아닙니다. 두 근사는 서로 다른 공간 영역을 설명합니다.

겹치는 부분은 상수 F−k입니다. δ≪y≪1인 영역이 존재한다는 설명은 δ→0의 점근적 상황에 해당합니다. δ가 크더라도 식의 계산은 가능하지만 이러한 겹침 영역이 있다고 해석해서는 안 됩니다.

8. 공통 부분을 한 번 빼서 합성해를 만든다

두 근사를 더한 뒤 중복으로 더해진 공통값을 뺍니다.

$$v_{\mathrm{comp}}(y)=v_{\mathrm{out}}(y)+v_{\mathrm{in}}(y/\delta)-(F-k) =F+k(y-1)-D e^{-y/\delta}.$$

가까운 경계에서는 v_comp(0)=N입니다. q=e^(−1/δ)라고 쓰면 먼 경계의 잔차는 v_comp(1)−F=−Dq입니다. 합성해는 전체 구간의 좋은 근사이지만, 정확히 정규화한 해와 같지는 않습니다.

기본 예제의 합성해는 1−e^(−10x)입니다. 정확해와 거의 겹치는 그림만으로 오차를 판단하지 말고 다음 절에서 차이를 직접 계산해 봅시다.

9. 전체 구간 오차를 증명한다

정확해에서 합성해를 빼면

$$v(y)-v_{\mathrm{comp}}(y) =Dq\frac{1-e^{-y/\delta}}{1-q},\qquad q=e^{-1/\delta}$$

입니다. 비율은 [0,1]에서 0부터 1까지 단조롭게 증가합니다. 따라서 절댓값 최대는 먼 끝점에서 얻어집니다.

$$\|v-v_{\mathrm{comp}}\|_\infty=|D|q.$$

이 유한 모형에서는 상계만이 아니라 정확한 등식입니다. 그림에 찍힌 몇 개 점 중 최대 오차를 고르는 것보다 강한 보증입니다. 다만 상수 계수 방정식과 주어진 경계조건으로 증명한 결과입니다. 모든 정합 점근 전개가 언제나 지수적으로 작은 오차를 갖는다는 일반 명제로 확대하지 않습니다.

10. 점별 수렴은 크기 1의 경계 오차를 숨길 수 있다

바깥해와 정확해의 차이는

$$v(y)-v_{\mathrm{out}}(y) =-D\frac{e^{-y/\delta}-q}{1-q},\qquad \|v-v_{\mathrm{out}}\|_\infty=|D|$$

입니다. 고정한 y>0에서는 δ→0일 때 0으로 가지만 y=0에서는 −D로 남습니다. 그러므로 경계층이 있는 끝점을 제외한 점별 수렴은 성립해도, D≠0이면 닫힌 구간 전체에서 균일하게 수렴하지 않습니다. 끝점의 값과 안쪽에서 끝점으로 접근한 극한의 순서를 바꾸면 결과가 달라집니다.

평균 오차도 이 문제를 숨길 수 있습니다. 위 식을 제곱하고 분자의 절댓값을 e^(−y/δ)로 제한하면 L² 노름은 |D|√(δ/2)/(1−q) 이하입니다. 평균 손실이 작아져도 최악의 경계 오차가 작다는 보증은 아닙니다.

11. 작은 값 오차와 큰 기울기를 구분한다

정확해를 미분하면

$$v_y=k+\frac{D}{\delta}\frac{e^{-y/\delta}}{1-q},\qquad v_{yy}=-\frac{D}{\delta^2}\frac{e^{-y/\delta}}{1-q}$$

입니다. 고정한 D≠0에서 경계층의 기울기와 곡률은 각각 1/δ, 1/δ² 규모입니다. 작은 ε가 큰 곡률에 곱해지므로, 이 항을 지우면 경계 근처의 균형 자체가 사라집니다.

합성해의 일차 미분 오차 최대는 |D|q/[δ(1−q)]이며 가까운 끝점에서 얻어집니다. 값의 오차 |D|q와는 다른 양입니다. 반사한 경우 u_x=−v_y, u_xx=v_yy라는 점도 기억하세요. 기울기를 절댓값만으로 표시하면 변화 방향을 놓칠 수 있습니다.

12. 내부 방정식 잔차가 0이어도 경계값 문제를 푼 것은 아니다

바깥해와 합성해는 모두 내부 잔차 εv_yy+bv_y−f가 정확히 0입니다. 그런데 바깥해는 가까운 경계를 D만큼, 합성해는 먼 경계를 −Dq만큼 틀립니다. 미분방정식 잔차만으로 경계값 문제의 해를 결정할 수 없습니다.

AI 대리모형이나 PINN을 평가한다면 경계값, 내부 잔차, 값 오차, 미분 오차를 따로 확인해야 합니다. 경계층 밖에서 주로 뽑은 표본은 급한 변화를 놓치고도 좋아 보일 수 있습니다. 이 예제는 이를 확인할 수 있는 통제된 기준 문제이며, 특정 신경망이나 최적화기가 반드시 실패한다는 실험 결과는 아닙니다.

13. 작은 매개변수가 언제나 눈에 보이는 층을 만들지는 않는다

D=0이면 정확해는 N+ky이고 바깥해·합성해와 전 구간에서 같습니다. 방정식의 특이 구조는 경계층을 허용하지만, 이 경계 자료는 지수 모드를 여기하지 않습니다.

a=1, f=1, A=0, B=1이면 D=0이고 모든 ε>0에서 u=x입니다. a=1, f=1, A=B=0으로 바꾸면 D=−1이라 경계층이 남습니다. D의 부호는 보정 방향을 바꾸지만, 감쇠 모드가 위치하는 경계는 a의 부호가 결정합니다.

14. 흐름 계수 0은 다른 퇴화 문제다

a=0에서는 a나 b로 나누지 않습니다. 두 번 직접 적분하면

$$u(x)=A+(B-A)x+\frac{f}{2\varepsilon}x(x-1)$$

입니다. f≠0인 채 ε→0이면 내부 값은 보통 1/ε 규모로 커집니다. f=0이면 양 끝을 잇는 직선입니다. ε=0인 축약 방정식은 f≠0일 때 모순인 0=f이고, f=0일 때에는 해를 정하지 못하는 0=0입니다. 비영 흐름 계수의 바깥 문제와 같은 것으로 취급하지 않습니다.

계산기는 아주 작은 양의 a를 0으로 반올림하지 않습니다. 유리수 입력과 방향을 보존한 지수함수 포함 구간으로 이 구분을 유지합니다. 단순한 부동소수점 식에서 비슷한 수를 빼며 정확도를 잃을 수 있다는 점과 물리적으로 a=0이라는 주장은 다릅니다.

15. 격자 해상도와 점근 근사의 적용 조건을 분리한다

m개 등간격 패널이면 h=1/m이고 h/δ=b/(mε)입니다. 한 물리 격자 칸이 늘린 좌표에서 얼마나 먼 거리를 건너뛰는지 나타냅니다. 첫 칸 안에 들어가는 지수 모드 변화의 비율은 1−e^(−h/δ)입니다.

화면은 mδ≥4인지 표시합니다. e배 감쇠 거리당 패널 4개 이상이라는 명시적인 시각화 기준이며, 유한차분법 수렴 정리가 아닙니다. δ≤0.1을 비교적 얇은 층으로 표시하는 것도 화면상의 분류 규약입니다.

물리 좌표 65점 그림은 아주 얇은 층을 놓칠 수 있습니다. 별도의 안쪽 좌표 65점 그림은 0≤X≤min(8,1/δ)를 표본화합니다. 두 관점을 함께 봐야 합니다. X=8에서 그림을 끝냈다고 지수함수가 정확히 0이 된 것은 아닙니다. 전체 구간의 증명된 오차 보증도 이 두 그림의 격자에 의존하지 않습니다.

16. 수치 포함 구간과 정확한 항등식을 읽는 법

매개변수는 유리수로 유지하고, 음의 지수함수는 192비트 출력 격자에서 방향을 보존한 포함 구간으로 계산합니다. 이후 유리수 구간 연산으로 반올림 오차를 전파합니다. 이 연산 격자가 최종 해의 정확한 192비트를 약속하지는 않습니다. 작은 정규화 분모로 나누거나 큰 계수를 곱하면 최종 구간은 넓어질 수 있습니다.

지수함수 구간의 아래 끝점이 0인 것은 표현 해상도보다 작다는 뜻일 수 있으며, 지수함수가 정확히 사라졌다는 뜻은 아닙니다. 양의 ε와 비영 흐름 계수라는 분류는 그대로 유지합니다. 정확한 경계값·흐름 계수 0의 다항식·D=0인 항등식은 대수적으로 처리합니다.

독립 참조 계산은 원래 부호가 있는 좌표의 적분인자 해와 고정밀 expm1을 사용합니다. 생산 계산은 반사와 감쇠 지수함수를 사용합니다. 수치 일치는 검산이고, 보증의 이유는 앞에서 유도한 항등식과 방향을 보존한 연산입니다.

17. 손계산 예제와 해설 문제

1. 기본 왼쪽 층. ε=0.1, a=1, f=0, A=0, B=1에서 δ, D, 합성해 최대 오차를 구하세요. 각각 0.1, 1, e^(−10)≈0.00004540입니다. 바깥해 최대 오차는 여전히 1입니다.
2. 흐름 반전. a만 −1로 바꾸세요. 가까운 경계는 x=1이고 N=1, F=0, D=−1입니다. 층이 오른쪽으로 이동합니다. 정준 기울기는 음수지만 반사 때문에 물리적 기울기는 양수입니다.
3. 맞는 경계조건. a=f=B=1, A=0이면 D=0이어서 u=x입니다. 명목상의 δ가 작다는 사실만으로 비영 경계층 진폭을 결론 내릴 수 없습니다.
4. 잔차의 함정. 첫 문제에서 u=1을 방정식에 대입하세요. 내부 잔차는 0이지만 u(0)=1이 지정한 0을 위반합니다. 경계 검사가 빠지면 해가 아닌 함수를 통과시킵니다.
5. 더 얇은 층과 격자. δ=0.01, m=4이면 h/δ=25라 변화 대부분을 한 칸이 건너뜁니다. m=256이면 0.390625이며 e배 거리당 2.56패널이어서 화면 목표 4에는 아직 못 미칩니다. 어느 표시도 수치 해법의 오차 증명은 아닙니다.
6. 약한 흐름. ε=0.1, a=0.001이면 δ=100입니다. ε<1만을 근거로 얇은 층의 선도 근사라고 해석하면 안 됩니다.
7. 흐름 없이 힘만 있는 경우. a=0, f=1, ε=0.1, A=0, B=1이면 중앙값은 1/2−1/(8ε)=−0.75입니다. 음수로 내려가는 것은 이 방정식 부호에서의 정확한 해이며 계산 불안정의 증거가 아닙니다.
8. 미분항의 크기. 왜 작은 εu″를 남기나요? 진폭이 1 정도이고 δ=ε/b인 층에서는 εu″와 bu′가 모두 b²/ε 규모입니다. 두 항의 상쇄가 축약 방정식에서 사라진 물리적 균형입니다.

18. 연구 연결과 남은 범위

실험에서는 그림을 보기 전에 층의 위치와 폭부터 예측하세요. 물리 좌표와 늘린 좌표, 바깥해 오차와 합성해 보증, 내부 잔차와 양 끝 경계 오차를 나란히 비교합니다. 입력을 하나씩 바꾸고 어떤 가정이 달라졌는지 기록하세요.

이 장은 상수 계수 선형 경계값 문제를 다룹니다. 가변 계수, 전환점, 여러 경계층이나 내부층, 비선형 정합, WKB, 일반 PDE 수치해법의 안정성 보증은 추가 학습이 필요합니다. 다음 섭동 장에서는 또 다른 비균일 극한인 세속항과 장시간 위상 오차를 다룹니다.

경계층의 척도와 정합에 대한 배경은 John K. Hunter §4.3, 균일 점근성의 정의는 NIST DLMF §2.1(iv)를 참고하세요. 이 장의 방정식과 전체 구간 보증은 본문에서 직접 유도했습니다. Hunter의 예시 방정식에는 이 장과 다른 영차 항이 포함되어 있습니다.