콜-호프 변환

편집 역사 토론
수치해석 전산유체역학 마지막 수정: 2026-08-26 04:14:27

1. 개요[편집]

콜-호프 변환
Cole–Hopf transformation
변환식u = −2ν φx / φ
보내는 곳점성 버거스 방정식 → 열방정식
명명에버하르트 호프(1950) · 줄리언 콜(1951)
성격비선형 PDE의 완전 선형화 (C-적분가능)
부산물ν→0 극한에서 엔트로피 해 자동 선택
실무 함정exp(−1/2ν · …) 언더플로

비선형항을 없애는 방법이 있다. 로그를 씌우는 것이다. 물론 아무 방정식에나 되는 건 아니다.

콜-호프 변환(Cole–Hopf transformation)은 점성 버거스 방정식의 해 uu를 u=−2ν φx/φu=-2\nu\,\varphi_x/\varphi 로 치환해, 비선형 방정식을 근사 없이 정확히 선형 열방정식 φt=νφxx\varphi_t=\nu\varphi_{xx} 로 바꾸는 치환이다.

ut+uux=νuxx→  u=−2ν φx/φ  φt=νφxxu_t + u u_x = \nu u_{xx} \quad\xrightarrow{\;u=-2\nu\,\varphi_x/\varphi\;}\quad \varphi_t = \nu \varphi_{xx}

이게 왜 대단한 일인지는 반대편을 보면 안다. 비선형 편미분방정식은 원칙적으로 중첩 원리가 없고, 해의 존재조차 보장되지 않으며, 3차원 나비에-스토크스 방정식은 매끄러운 해가 있는지조차 아직 아무도 모른다. 그런데 그 나비에-스토크스에서 압력과 횡방향만 지운 버거스는, 치환 한 줄로 학부 2학년이 푸는 열방정식이 되어 버린다. 임의의 초기조건에 대한 명시적 적분해가 존재하고, 충격이 생겼다 사라지는 전 과정을 닫힌 형태로 추적할 수 있다.

이 행운이 수치해석에 주는 것은 두 가지다. 하나는 참값(비선형 충격을 품은 문제인데 정답을 안다)이고, 다른 하나는 엔트로피 해 선택 원리의 구성적 증명이다. ν→0\nu\to0 극한에서 적분해의 라플라스 근사가 물리적으로 옳은 약해 하나를 저절로 집어낸다.1

2. 유도 — 퍼텐셜을 도입하고 로그를 씌운다[편집]

두 단계다. 원래 방정식은 uu에 대해 2계이지만, 우선 보존형으로 쓰고 한 번 적분한다.

ut+(12u2)x=νuxxu_t + \left(\tfrac{1}{2}u^2\right)_x = \nu u_{xx}

전 항이 xx-미분 형태이므로 u=ψxu=\psi_x 로 퍼텐셜을 도입하면 괄호 밖의 ∂x\partial_x 를 통째로 떼어낼 수 있다(적분상수는 ψ\psi의 자유도로 흡수).

ψt+12ψx2=νψxx\psi_t + \tfrac{1}{2}\psi_x^2 = \nu \psi_{xx}

이것이 점성 해밀턴-야코비 방정식이다. 남은 비선형항은 12ψx2\tfrac12\psi_x^2 하나뿐이고, 이 항은 지수함수의 미분이 자기 자신을 곱해서 만들어 내는 그 항과 모양이 같다. 그래서 두 번째 단계:

ψ=−2νln⁡φ⟺φ=exp⁡ ⁣(−ψ2ν)\psi = -2\nu \ln \varphi \quad\Longleftrightarrow\quad \varphi = \exp\!\left(-\frac{\psi}{2\nu}\right)

대입하면 ψt=−2νφt/φ\psi_t=-2\nu\varphi_t/\varphi, ψx=−2νφx/φ\psi_x=-2\nu\varphi_x/\varphi, ψxx=−2ν(φxx/φ−φx2/φ2)\psi_{xx}=-2\nu(\varphi_{xx}/\varphi-\varphi_x^2/\varphi^2) 이고, 12ψx2=2ν2φx2/φ2\tfrac12\psi_x^2=2\nu^2\varphi_x^2/\varphi^2 항이 우변의 +2ν2φx2/φ2+2\nu^2\varphi_x^2/\varphi^2 와 정확히 상쇄된다. 남는 것은

−2νφ(φt−νφxx)=0-\frac{2\nu}{\varphi}\left(\varphi_t - \nu\varphi_{xx}\right) = 0

즉 φt=νφxx\varphi_t=\nu\varphi_{xx}. 그리고 u=ψx=−2ν(ln⁡φ)x=−2νφx/φu=\psi_x=-2\nu(\ln\varphi)_x=-2\nu\varphi_x/\varphi 가 처음에 적어 둔 변환식이다.

상쇄가 일어나는 이유를 한 문장으로 요약하면 이렇다. 로그 미분은 곱을 합으로 바꾸는 대신 제곱항 하나를 만들어 내는데, 버거스의 비선형항이 하필 정확히 그 제곱항이다. 계수 −2ν-2\nu 는 12u2\tfrac12 u^2 의 12\tfrac12 와 확산계수 ν\nu 를 맞추기 위한 것이라 임의로 못 바꾼다. 이류항 계수가 다르거나 확산이 uu에 의존하면 이 마술은 즉시 깨진다.

3. 초기조건과 적분해[편집]

변환은 초기조건도 같이 데리고 간다. ψ0(x)=∫0xu0(s) ds\psi_0(x)=\int_0^x u_0(s)\,ds 이므로

φ(x,0)=φ0(x)=exp⁡ ⁣[−12ν∫0xu0(s) ds]\varphi(x,0) = \varphi_0(x) = \exp\!\left[-\frac{1}{2\nu}\int_0^x u_0(s)\,ds\right]

이고, 열방정식은 가우스 핵으로 명시적으로 풀린다.

φ(x,t)=14πνt∫−∞∞exp⁡ ⁣[−(x−ξ)24νt]φ0(ξ) dξ\varphi(x,t) = \frac{1}{\sqrt{4\pi\nu t}}\int_{-\infty}^{\infty} \exp\!\left[-\frac{(x-\xi)^2}{4\nu t}\right]\varphi_0(\xi)\,d\xi

두 지수를 합치면 −12ν[∫0ξu0+(x−ξ)22t]-\frac{1}{2\nu}\big[\int_0^\xi u_0 + \frac{(x-\xi)^2}{2t}\big] 하나로 정리되므로, u=−2νφx/φu=-2\nu\varphi_x/\varphi 를 계산할 때 ∂x\partial_x 가 지수에서 (x−ξ)/(2νt)(x-\xi)/(2\nu t) 를 끌어내리고 앞의 −2ν-2\nu 와 곱해져 깔끔하게 (x−ξ)/t(x-\xi)/t 가 된다.

u(x,t)=∫−∞∞x−ξt e−G(ξ;x,t)/2ν dξ∫−∞∞e−G(ξ;x,t)/2ν dξ,G(ξ;x,t)=∫0ξu0(s) ds+(x−ξ)22tu(x,t) = \frac{\displaystyle\int_{-\infty}^{\infty}\frac{x-\xi}{t}\,e^{-G(\xi;x,t)/2\nu}\,d\xi} {\displaystyle\int_{-\infty}^{\infty}e^{-G(\xi;x,t)/2\nu}\,d\xi}, \qquad G(\xi;x,t)=\int_0^{\xi}u_0(s)\,ds+\frac{(x-\xi)^2}{2t}

해가 가중평균이라는 점을 주목할 만하다. 분모는 정규화 상수, 분자는 “ξ\xi에서 출발해 시각 tt에 xx에 도착하는 직진 입자의 속도 (x−ξ)/t(x-\xi)/t“이고, 가중치는 e−G/2νe^{-G/2\nu} 다. 즉 콜-호프 해는 모든 가능한 출발점에 대한 볼츠만 가중 평균이며, ν\nu 가 온도 역할을 한다. 통계역학 냄새가 나는 게 우연이 아니라서, 뒤에 나올 KPZ 방정식 연결도 정확히 이 구조에서 온다.

주기 문제에서는 적분 대신 푸리에 급수를 써도 된다. 예컨대 u0=−sin⁡(πx)u_0=-\sin(\pi x) 인 국룰 시험문제에서는 φ0=exp⁡ ⁣[(1−cos⁡πx)/(2πν)]\varphi_0=\exp\!\left[(1-\cos\pi x)/(2\pi\nu)\right] 라 계수가 변형 베셀 함수 In ⁣(12πν)I_n\!\left(\tfrac{1}{2\pi\nu}\right) 로 나오고, 해는 사인·코사인 급수의 비로 적힌다. 예뻐 보이지만 ν\nu 가 작아지면 InI_n 이 지수적으로 커져 곧바로 오버플로한다 — 뒤의 함정 절에서 다시 나온다.

4. ν → 0 — 라플라스 근사가 엔트로피 해를 고른다[편집]

ν\nu 가 작으면 e−G/2νe^{-G/2\nu} 는 GG 의 최소점 근방에만 질량이 몰린 뾰족한 함수가 된다. 이것이 정확히 라플라스 근사(최급강하법)가 다루는 상황이고, 결과는

u(x,t)  ⟶  x−ξ∗(x,t)t,ξ∗(x,t)=arg⁡min⁡ξ[ψ0(ξ)+(x−ξ)22t]u(x,t) \;\longrightarrow\; \frac{x-\xi^\ast(x,t)}{t}, \qquad \xi^\ast(x,t) = \arg\min_{\xi}\left[\psi_0(\xi)+\frac{(x-\xi)^2}{2t}\right]

이다. 이것이 무점성 버거스의 락스-올레이닉(Lax–Oleinik) 공식이며, 우변의 최소화는 ψ0\psi_0 와 포물선의 하한 합성곱(inf-convolution) — 즉 르장드르 변환 세계의 모로 포락선이다. 여기서 세 가지가 공짜로 따라 나온다.

  • 충격의 위치. ξ∗\xi^\ast 는 최소점이 유일한 동안 xx에 대해 연속이지만, 최소점이 둘로 갈리는 xx 에서 점프한다. 그 점프가 바로 충격이고, uu 의 불연속이 거기서 생긴다. 특성선이 교차한 뒤 어느 것을 살릴지 고민할 필요가 없다 — 최소화 문제가 이미 결정했다.
  • 엔트로피 조건. ξ∗\xi^\ast 는 xx의 비감소 함수임을 볼록성으로 보일 수 있다. 그러면 점프에서 ξ∗\xi^\ast 가 위로 뛰므로 u=(x−ξ∗)/tu=(x-\xi^\ast)/t 는 아래로 뛴다. 즉 uL>uRu_L>u_R — 랭킨-위고니오 조건에 붙는 락스 엔트로피 조건이 결론으로 나온다. 팽창 충격은 애초에 생길 수 없다.
  • 올레이닉 단측 립시츠 조건. ξ∗\xi^\ast 가 비감소이므로 ux≤1/tu_x \le 1/t 가 모든 곳에서 성립한다. 초기조건이 아무리 험해도 증가하는 방향의 기울기에는 1/t1/t 라는 보편 상한이 걸린다는 뜻이다(급한 감소, 즉 충격은 허용된다). 초기 데이터와 무관하게 시간만 지나면 해가 저절로 정리된다는 것이고, 이것이 스칼라 보존법칙 해의 정칙화(regularizing effect) 정리의 원형이다.

정리하면 콜-호프는 “점성 소실 극한이 엔트로피 해다”라는 명제를, 존재성 정리가 아니라 계산 가능한 적분 하나로 보여 준다. 스칼라 볼록 보존법칙에서 이 등가성이 그림까지 그려지는 유일한 경로라, 유한체적법 강의의 엔트로피 절이 항상 버거스로 시작하는 이유가 이것이다.

5. 검증용 정확해로 쓰기 — 그리고 언더플로[편집]

콜-호프의 실전 용도 1순위는 검증 및 확인의 참값 공급이다. 충격을 품은 문제의 L1L^1 수렴 차수를 주장하려면 정답이 필요한데, 매끈한 제조해(MMS)로는 충격을 못 만들고 충격을 만들면 정답을 모른다. 버거스는 이 딜레마를 빠져나가는 거의 유일한 시험대다.

콜-호프 변환 u = −2ν φ_x/φ 로 얻은 정확해(굵은 선)와, 같은 초기조건 u₀ = −sin πx 를 4차 중심차분 + RK4 로 직접 전진시킨 수치해(점)를 같은 축에 겹친다 — 정확해 쪽은 시간전진이 없고, 열핵 적분을 x 마다 지수의 최대점을 빼서 안정화한 뒤 심프슨 구적으로 한 번에 계산한 값이다. ν = 0.1 에서는 두 곡선의 전 구간 max 차이가 9.8e−7 로 완전히 겹치지만, 격자를 N = 160 에 고정한 채 ν 를 0.002 까지 내리면 1.5e−1 까지 벌어지고 충격 근처가 2Δx 톱니로 흔들린다(셀 레이놀즈 Δx/ν = 6.25). 아래 패널의 max|u_x| 는 비점성 특성선이 부서지는 t_b = 1/π = 0.3183 에서 치솟다가 점성이 두께 δ ≈ 4.2ν 에서 멈춘다.

문제는 실제로 계산해 보면 시작한다. 흔히 쓰는 시험 조건 ν=0.01/π\nu=0.01/\pi 에서 1/2ν≈1571/2\nu \approx 157 이고, 지수의 인자 GG 가 O(1)O(1) 이면 e−157e^{-157} — 배정밀도 하한(≈10−308\approx 10^{-308})에는 아직 여유가 있다. 그런데 ν\nu 를 한 자릿수만 더 내려 ν=10−4\nu=10^{-4} 로 가면 1/2ν=50001/2\nu=5000 이라 e−5000e^{-5000} — 완전한 언더플로다. 분자와 분모가 동시에 정확히 0이 되고 해는 0/0 = NaN이 된다. 게다가 GG 의 값 자체는 u0u_0 의 크기와 정의역 폭에 비례해 커지므로, ν\nu 가 그리 작지 않아도 초기 진폭이 크면 같은 사고가 난다. “정확해”를 뽑겠다고 짠 코드가 정작 관심 있는 고레이놀즈 영역에서 제일 먼저 죽는 것이다.

처방은 통계·기계학습 쪽에서 익숙한 그것과 같다. 로그 도메인에서 최댓값을 빼는 log-sum-exp 트릭.

u=∑kwkx−ξktexp⁡ ⁣[−Gk−Gmin⁡2ν]∑kwkexp⁡ ⁣[−Gk−Gmin⁡2ν]u = \frac{\sum_k w_k \frac{x-\xi_k}{t} \exp\!\left[-\frac{G_k-G_{\min}}{2\nu}\right]} {\sum_k w_k \exp\!\left[-\frac{G_k-G_{\min}}{2\nu}\right]}

Gmin⁡=min⁡kGkG_{\min}=\min_k G_k 를 분자·분모에서 공통으로 빼면 지수의 최댓값이 정확히 11 이 되고, 나머지는 전부 그보다 작으니 언더플로해도 그냥 0으로 죽어 무해하다. 비율이므로 e−Gmin⁡/2νe^{-G_{\min}/2\nu} 는 어차피 상쇄된다. 실무에서 함께 챙길 것들:

함정처방
지수 언더플로·오버플로Gmin⁡G_{\min} 빼기(log-sum-exp), 로그 도메인 누산
ν\nu 작을 때 피적분함수가 바늘처럼 뾰족ξ∗\xi^\ast 근방으로 적분구간 국소화, 가우스-에르미트 수치적분
t→0t\to0 에서 (x−ξ)2/2t(x-\xi)^2/2t 발산t=0t=0 은 초기조건을 직접 쓰고, 작은 tt 는 ξ\xi 를 xx 기준 스케일 변수로 치환
최소점이 둘인 지점(충격)두 봉우리를 각각 적분해 더한다. 한쪽만 잡으면 충격 위치가 격자 크기만큼 틀어진다
베셀 급수형 사용 시In(z)I_n(z) 대신 지수 스케일된 e−zIn(z)e^{-z}I_n(z) 루틴을 쓴다

이 절차만 지키면 참값을 기계 정밀도 근처까지 뽑을 수 있고, 그때부터 고두노프 도식·MUSCL·WENO 도식의 실제 수렴 차수를 정직하게 잴 수 있다. 반대로 콜-호프를 솔버로 쓰는 것 — 열방정식을 푸리에 변환으로 정확 적분한 뒤 되돌리는 방식 — 은 매력적이지만 같은 언더플로 벽에 부딪히므로, 실제 계산 코드는 여전히 의사스펙트럼법이나 유한체적으로 직접 푼다.2

6. 어디까지 되고 어디부터 안 되는가[편집]

콜-호프가 통하는 확장은 생각보다 넓다.

  • 다차원 퍼텐셜 버거스. u=∇ψ\mathbf{u}=\nabla\psi (비회전)인 벡터 버거스 ut+(u⋅∇)u=ν∇2u\mathbf{u}_t+(\mathbf{u}\cdot\nabla)\mathbf{u}=\nu\nabla^2\mathbf{u} 는 u=−2ν∇φ/φ\mathbf{u}=-2\nu\nabla\varphi/\varphi 로 그대로 φt=ν∇2φ\varphi_t=\nu\nabla^2\varphi 가 된다. 우주론의 접착 모형이 3차원에서도 닫힌 형태로 다뤄지는 이유.
  • 소스·강제항이 있는 경우. ut+uux=νuxx+F′(x)u_t+uu_x=\nu u_{xx}+F'(x) 꼴이면 φ\varphi 방정식에 퍼텐셜 항이 붙어 허수시간 슈뢰딩거 방정식이 된다. 여전히 선형이므로 고유함수 전개가 가능하다.
  • KPZ 방정식. ht=νhxx+λ2hx2+ηh_t=\nu h_{xx}+\tfrac{\lambda}{2}h_x^2+\eta 에 같은 지수 치환을 하면 곱셈 잡음이 붙은 확률 열방정식이 된다. 계면 성장 이론에서 “콜-호프 해”라는 말이 나오는 게 이 맥락이고, 이토 항 때문에 형식적 치환과 엄밀한 해가 어긋나는 문제(호헨버그 상수 이동)가 유명하다.
  • 백룬드 변환의 특수 사례. 두 방정식의 해를 짝지어 주는 미분관계라는 관점에서 보면 콜-호프는 버거스와 열방정식을 잇는 백룬드 변환이다.

안 되는 쪽이 더 중요하다. 콜-호프는 방법이 아니라 사고에 가깝다. 칼로제로(Calogero)의 분류를 빌리면 버거스는 변수변환만으로 선형화되는 C-적분가능계이고, KdV 방정식처럼 역산란 변환이 필요한 S-적분가능계와는 급이 다르다. 그리고 두 부류 모두 적분가능계라는 좁은 섬이다.

  • 계(system)로 못 간다. 오일러 방정식이나 나비에-스토크스에는 대응하는 치환이 없다. 압력이 비국소 구속을 만들고 와도(vorticity) 늘림이 3차원 비선형성의 본질인데, 로그 치환은 그중 무엇도 다루지 못한다.
  • 플럭스가 볼록 2차가 아니면 깨진다. ut+f(u)x=νuxxu_t+f(u)_x=\nu u_{xx} 에서 콜-호프가 성립하는 것은 사실상 f=u2/2f=u^2/2 뿐이다. 그래서 대류-확산 방정식 계열의 다른 비선형 플럭스는 여전히 수치적으로만 다룬다.
  • 경계가 있으면 번거롭다. 무한 구간·주기 구간에서는 열핵이 즉시 쓰이지만, 디리클레·노이만 경계에서는 φ\varphi 의 경계조건이 uu 의 것과 다른 형태로 변환되어 상(image) 전개가 필요하다.

즉 콜-호프의 진짜 교훈은 “비선형도 잘 하면 풀린다”가 아니라, “이렇게까지 잘 맞아떨어지는 사례가 1차원 버거스 하나뿐이라는 사실” 이다. 나머지 전부에 대해 우리는 여전히 격자를 깔고 계산기를 갈아 넣는다.3

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 이름은 에버하르트 호프(E. Hopf, 1950)와 줄리언 콜(J. D. Cole, 1951)에서 왔지만, 두 사람이 최초는 아니다. 같은 치환이 이미 1906년 포사이스(A. R. Forsyth)의 미분방정식 교과서에 실려 있었고, 그 이전에 백룬드 변환 문헌에도 흔적이 있다. 방정식 이름이 뷔르허스가 아니라 베이트먼이었어야 한다는 이야기와 완전히 같은 구조의 사건이 변환 이름에서도 한 번 더 반복된 셈. 참고로 콜의 1951년 논문 제목에 붙은 “occurring in aerodynamics”는 약한 충격파의 천음속 근사 맥락이었다. ↩

  2. “정확해”라는 단어가 주는 안정감이 위험하다. 콜-호프 적분을 부주의하게 짜면 오차가 도식의 오차보다 커져서, 열심히 만든 5차 WENO가 1.8차로 수렴하는 것처럼 보인다. 이 경우 범인은 도식이 아니라 참값이다. 참값 코드에도 격자 수렴 시험을 해야 한다는 뜻인데, 이 자기 참조적 상황을 처음 겪으면 꽤 허탈하다. ↩

  3. 그래서 학계에서는 “새 비선형 방정식을 선형화했다”는 주장에 반사적으로 의심이 붙는다. 대개는 (가) 특수해 몇 개를 찾았거나, (나) 이미 알려진 C-적분가능계를 다시 발견했거나, (다) 어딘가에서 조용히 항을 하나 버렸다. 셋 다 아니면 정말로 큰 발견이지만, 그런 일은 70년에 한 번쯤 일어난다. ↩