콜-호프 변환

오류 제보
수치해석 전산유체역학 마지막 수정: 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년에 한 번쯤 일어난다. ↩

이 문서가 문제를 푸는 데 도움이 됐나요? 한 번의 클릭이 다음 검수 순서를 정합니다.

의견을 반영할게요. 감사합니다.

다음으로 읽기 검증 및 확인 약 5분 더 탐색 →