멜니코프 방법

1. 개요[편집]

멜니코프 방법 Melnikov method
제시	V. K. 멜니코프 (1963)
대상	주기 섭동을 받는 평면 해밀턴 계
재는 것	Wu와 Ws 사이의 1차 거리
판정	단순 영점 존재 → 횡단 교차 → 말굽
대표 결과	감쇠·강제 더핑의 카오스 문턱 공식
한계	1차 섭동 — ε가 작아야 하고, 문턱은 하한

멜니코프 방법은 호모클리닉 고리를 가진 평면계에 작은 주기 섭동을 얹었을 때, 깨진 안정·불안정 다양체가 서로 횡단적으로 교차하는지를 1차 섭동으로 판정하는 기법이다. 판정은 하나의 적분값 — 멜니코프 함수 $M(t_0)$ — 의 부호 변화 여부로 끝나고, 단순 영점이 있으면 스메일 말굽의 존재, 따라서 카오스적 불변집합의 존재가 따라 나온다.

이 방법의 가치는 카오스 문턱을 손으로 계산할 수 있게 해 준다는 데 있다. 카오스는 보통 수치 실험으로 “돌려 보고 알게 되는” 대상인데, 멜니코프는 감쇠·강제 진폭·진동수만 넣으면 닫힌 형태의 부등식을 뱉는다. 비선형 진동, 해양 구조물 전복, 플러터의 극한주기진동 같은 곳에서 여전히 쓰이는 이유가 이것이다.

호모클리닉 궤도 자체의 정의와 얽힘·말굽의 기하는 호모클리닉 궤도에 있으니 여기서는 적분이 어디서 나오고, 실제로 어떻게 계산하며, 어디까지 믿어도 되는지에 집중한다.

2. 무대 설정[편집]

비섭동 계는 평면 해밀턴 계다.

$$\dot{\mathbf x} = \mathbf f(\mathbf x), \qquad \mathbf x\in\mathbb R^2,\quad \nabla\cdot\mathbf f = 0$$

여기에 쌍곡 안장 $\mathbf p_0$ 와 그것에 붙은 호모클리닉 궤도 $\mathbf q_0(t)$ 가 있다고 하자. 해밀턴 계에서는 보존량이 조건 하나를 자동으로 채워 주므로 이런 고리는 매개변수 조율 없이 공짜로 존재한다(진자의 분리선이 표본).

여기에 주기 $T = 2\pi/\omega$ 의 작은 섭동을 얹는다.

$$\dot{\mathbf x} = \mathbf f(\mathbf x) + \varepsilon\,\mathbf g(\mathbf x, t), \qquad 0<\varepsilon\ll1$$

시간을 상태로 승격시켜 $\mathbb R^2\times S^1$ 위의 자율계로 보면, 안장 $\mathbf p_0$ 는 쌍곡 주기궤도가 되고 그 안정·불안정 다양체는 2차원 면이 된다. 시간 $t_0$ 마다 한 번씩 찍는 스트로보스코프 사상(주기 $T$ 의 푸앵카레 단면) $P_{t_0}^{\varepsilon}$ 을 보면, 안장은 쌍곡 고정점 $\mathbf p_\varepsilon = \mathbf p_0 + O(\varepsilon)$ 이 되고 그 다양체는 원래 고리에 $O(\varepsilon)$ 가까이 놓인다. 법선 쌍곡성 덕분에 이 지속성이 보장된다는 점은 불변 다양체 문서의 페니첼 이야기 그대로다.

문제는 이것이다. 깨진 두 다양체가 서로 만나는가, 아니면 스쳐 지나가는가. 비섭동 상태에서는 완벽히 포개져 있던 것이 섭동으로 갈라지므로, 갈라진 부호 있는 간격을 1차까지 계산하면 답이 나온다.

3. 멜니코프 함수는 어디서 나오나[편집]

기준점 $\mathbf q_0(0)$ 에서 궤도에 수직한 선을 긋고, 그 선 위에서 $W^u$ 위의 점과 $W^s$ 위의 점의 부호 있는 거리를 재면

$$d(t_0) = \varepsilon\,\frac{M(t_0)}{\lVert \mathbf f(\mathbf q_0(0))\rVert} + O(\varepsilon^2)$$

가 되고, 여기 등장하는 것이 멜니코프 함수다.

$$M(t_0) = \int_{-\infty}^{\infty} \mathbf f\big(\mathbf q_0(t)\big) \wedge \mathbf g\big(\mathbf q_0(t),\, t+t_0\big)\, dt, \qquad \mathbf a\wedge\mathbf b = a_1b_2 - a_2b_1$$

유도의 골자는 의외로 짧다. 섭동된 다양체 위의 점을 $\mathbf q_0(t) + \varepsilon\,\mathbf q_1(t) + O(\varepsilon^2)$ 로 전개하고, 갈라짐을 재는 양

$$\Delta(t) = \mathbf f\big(\mathbf q_0(t)\big)\wedge \big(\mathbf q_1^{u}(t) - \mathbf q_1^{s}(t)\big)$$

를 만들어 미분하면, 1차 변분방정식으로부터

$$\frac{d\Delta}{dt} = \operatorname{tr} D\mathbf f\big(\mathbf q_0(t)\big)\,\Delta + \mathbf f\wedge\mathbf g$$

가 나온다. 비섭동 계가 해밀턴이라 $\operatorname{tr}D\mathbf f = 0$ 이므로 첫 항이 통째로 사라지고, 남는 것은 단순한 적분뿐이다. 멜니코프 공식이 유난히 깔끔한 이유가 이 한 줄이다. 비섭동 계가 해밀턴이 아니면 피적분 함수에 $\exp\!\big(-\int_0^t \operatorname{tr}D\mathbf f\big)$ 라는 가중치가 붙는다.

적분의 수렴은 공짜로 얻는다. $t\to\pm\infty$ 에서 $\mathbf q_0(t)$ 는 안장으로 지수적으로 접근하고 거기서 $\mathbf f\to 0$ 이므로, 피적분 함수가 지수적으로 죽는다. 수치적으로는 $\pm T_{\max}$ 에서 잘라도 오차가 지수적으로 작다는 뜻이라, 실무 계산이 아주 편하다.

4. 판정 규칙[편집]

$M$ 이 단순 영점을 가지면 — 즉 $M(t_0^*)=0$ 이고 $M'(t_0^*)\ne 0$ 인 $t_0^*$ 가 있으면 — 충분히 작은 $\varepsilon$ 에 대해 $W^u$ 와 $W^s$ 가 횡단적으로 교차한다. 음함수 정리를 $d(t_0)$ 에 적용한 결과이며, 미분이 0이 아니라는 조건이 “스치는 접촉이 아니라 가로지름”을 보장한다. 그다음은 호모클리닉 궤도에서 본 사슬이 자동으로 굴러간다.

$$\text{단순 영점} \;\Rightarrow\; \text{횡단 교차} \;\Rightarrow\; \text{얽힘} \;\Rightarrow\; \text{스메일-버코프} \;\Rightarrow\; \text{말굽}$$

반대로 $M$ 이 0에서 멀리 떨어져 부호를 유지하면 두 다양체는 갈라진 채 만나지 않고, 그 방향으로는 얽힘이 없다. $M$ 이 0을 갖되 미분도 함께 0인 경계 상황이 접촉 분기(homoclinic tangency)에 해당하며, 이 값이 실제 파라미터 평면에서의 문턱 곡선을 준다.

반드시 붙여야 할 단서 하나. 말굽은 카오스적 불변집합이지 끌개가 아니다. 말굽 위의 집합은 대개 측도 0의 안장형 집합이라, 문턱을 막 넘긴 계에서 관측되는 것은 지속적 카오스가 아니라 과도 카오스와 프랙탈 유역 경계인 경우가 많다. 즉 멜니코프 문턱은 “카오스가 보이기 시작하는 값”이 아니라 **“흡인 유역 경계가 얽히기 시작하는 값”**으로 읽는 것이 정확하다.

5. 고전 예제 — 감쇠·강제 더핑 방정식[편집]

가장 많이 인용되는 계산이 쌍우물 더핑 방정식이다.

$$\ddot x - x + x^3 = \varepsilon\,(\gamma\cos\omega t - \delta\,\dot x)$$

비섭동 계 $\ddot x = x - x^3$ 은 원점이 안장, $x=\pm1$ 이 중심인 쌍우물 해밀턴 계이고, 원점에 붙은 호모클리닉 해가 닫힌 형태로 나온다.

$$\mathbf q_0(t) = \big(\sqrt2\,\operatorname{sech} t,\ -\sqrt2\,\operatorname{sech} t\,\tanh t\big)$$

$\mathbf f = (y,\ x-x^3)$, $\mathbf g = (0,\ \gamma\cos\omega t - \delta y)$ 이므로 $\mathbf f\wedge\mathbf g = y\,(\gamma\cos\omega t - \delta y)$ 이고,

$$M(t_0) = \gamma\!\int_{-\infty}^{\infty}\! y_0(t)\cos\omega(t+t_0)\,dt \;-\; \delta\!\int_{-\infty}^{\infty}\! y_0(t)^2\,dt$$

두 적분 모두 초등적으로 계산된다. $\int y_0^2\,dt = 2\int \operatorname{sech}^2 t\,\tanh^2 t\,dt = 4/3$ 이고, $y_0$ 가 홀함수라 $\cos$ 성분은 죽고 $\sin$ 성분만 남아 $\int \operatorname{sech} t\,\tanh t\,\sin\omega t\,dt = \pi\omega\operatorname{sech}(\pi\omega/2)$ 를 쓴다. 결과는

$$M(t_0) = \sqrt2\,\pi\,\gamma\,\omega\,\operatorname{sech}\!\Big(\frac{\pi\omega}{2}\Big)\sin(\omega t_0) \;-\; \frac{4\delta}{3}$$

감쇠는 $t_0$ 에 무관한 음의 상수, 강제는 진폭 $\sqrt2\pi\gamma\omega\operatorname{sech}(\pi\omega/2)$ 의 사인파라는 구조가 그대로 드러난다. 단순 영점 조건은 사인파의 진폭이 상수항을 넘는 것이므로

$$\frac{\gamma}{\delta} \;>\; \frac{2\sqrt2}{3\pi\omega}\cosh\!\Big(\frac{\pi\omega}{2}\Big)$$

가 된다. $\omega=1$ 이면 우변이 약 $0.753$ 이다. 이 문턱 곡선은 $\omega\to0$ 과 $\omega\to\infty$ 에서 모두 발산하고 중간에 최솟값을 갖는데, $(\pi\omega/2)\tanh(\pi\omega/2)=1$ 을 풀면 $\omega\approx0.76$ 에서 $\gamma/\delta\approx0.711$ 로 최소가 된다. 즉 이 계를 가장 싸게 얽히게 만드는 강제 진동수가 존재하며, 그것은 우물 바닥의 선형 고유진동수와 같지 않다. 공진이 아니라 호모클리닉 궤도의 시간 폭이 정하는 값이라는 점이 재미있는 대목이다.¹

물리적 해석도 깔끔하다. $M$ 은 사실상 한 주기 동안 강제가 넣어 준 에너지와 감쇠가 빼 간 에너지의 차다. 강제가 감쇠를 이기는 위상 $t_0$ 가 존재하면 궤도가 고리 밖으로 밀려나며 다양체가 교차한다. 감쇠가 없으면($\delta=0$) $M$ 은 항상 부호를 바꾸므로 아무리 약한 주기 강제라도 얽힘이 생긴다 — 보존계가 섭동에 취약하다는 KAM 이론의 정서와 같은 방향이다.

6. 변종들[편집]

• 준조화 멜니코프 함수(subharmonic Melnikov). 호모클리닉 고리 안쪽의 주기궤도족을 대상으로 같은 적분을 유한 구간에서 계산하면, 주기 $m T/n$ 의 준조화 궤도가 어디서 태어나는지를 준다. 얽힘 문턱 아래에서 나타나는 공진 띠와 주기배가의 순서를 예측하는 데 쓰인다.
• 헤테로클리닉 버전. 두 안장을 잇는 궤도에 대해서도 같은 형태의 적분이 정의된다. 다만 이때는 두 다양체가 서로 다른 안장에 속해 있어 조건이 하나 더 붙는다.
• 2차 멜니코프. $M\equiv0$ 이 되는 퇴화 상황(대칭 때문에 1차 항이 통째로 상쇄되는 경우가 흔하다)에서는 $\varepsilon^2$ 항까지 가야 한다. 계산량이 급증하고, 이 지점에서 손계산을 포기하고 수치로 넘어가는 것이 보통이다.
• 고차원 확장. 다차원 계나 법선 쌍곡 불변 다양체에 붙은 호모클리닉 구조에도 멜니코프류 적분이 정의되며, 이것이 아널드 확산 논의의 기본 도구다. 3차원 안장-초점에 붙은 호모클리닉이 만드는 카오스는 성격이 달라 실니코프 분기 쪽 이론이 필요하다.

7. 수치적으로 다루기[편집]

적분 계산. 비섭동 호모클리닉 해가 닫힌 형태로 안 나오는 경우가 대부분이다. 그럴 때는 절단+사영 경계조건 BVP로 $\mathbf q_0(t)$ 를 수치적으로 구한 뒤(호모클리닉 궤도의 계산 절), 그 궤도를 따라 $M(t_0)$ 를 수치적분한다. 피적분 함수가 양 끝에서 지수적으로 죽으므로 절단 오차는 걱정거리가 아니고, 대신 궤도의 정확도가 그대로 $M$ 의 정확도가 된다.

얽힘 눈으로 보기. 문턱을 넘겼는지 확인하는 가장 직관적인 방법은 스트로보스코프 사상을 만들어 다양체를 대역적으로 그리는 것이다. 절차는 이렇다. 주기 $T$ 마다 상태를 기록하는 사상 $P$ 를 정의하고, 쌍곡 고정점 $\mathbf p_\varepsilon$ 을 뉴턴법으로 찾고, 그 고유벡터 방향으로 미세한 선분을 잡아 $P$ 를 반복 적용해 $W^u$ 를 성장시킨다. $W^s$ 는 $P^{-1}$ 로 같은 일을 한다(즉 시간을 거꾸로 적분). 문턱 아래에서는 두 곡선이 나란히 갈라져 있고, 문턱을 넘기면 로브가 무한히 늘어나며 서로를 관통하는 그림이 나온다.

유역으로 검증. 다양체를 그리기 귀찮으면 흡인 유역을 그려도 된다. 얽힘이 시작되는 순간부터 유역 경계가 프랙탈 손가락 모양으로 침식되기 시작하므로, 파라미터를 훑으며 불확실성 지수가 1에서 떨어지는 지점을 찾으면 그것이 실측 문턱이다. 구조·해양 분야에서 안전 유역 침식 연구가 멜니코프 문턱과 나란히 인용되는 이유가 이것이다.

파라미터 평면 지도. 실무 산출물은 대개 $(\omega,\gamma)$ 평면 위에 해석적 문턱 곡선을 그리고 그 위에 수치 실험 결과(최대 랴푸노프 지수 부호, 유역 프랙탈성)를 겹쳐 놓은 그림이다. 해석 곡선은 하한이므로 관측된 카오스 영역보다 안쪽에 놓이는 것이 정상이고, 둘이 정확히 겹치기를 기대하면 안 된다.²

8. 오용 주의[편집]

• $\varepsilon$ 은 정말 작아야 한다. 멜니코프는 1차 섭동 이론이다. 감쇠비 0.3짜리 계에 문턱 공식을 그대로 대입한 결과는 예측이 아니라 희망이다. 감쇠와 강제가 둘 다 $O(\varepsilon)$ 이어야 한다는 전제도 자주 잊힌다.
• 문턱은 필요조건이 아니라 충분조건 쪽이다. $M$ 에 단순 영점이 있으면 말굽이 있다. 그러나 $M$ 에 영점이 없다고 카오스가 없다는 보장은 못 한다 — 2차 이상의 효과나 다른 기구로 얼마든지 생길 수 있다.
• 말굽 ≠ 관측되는 카오스. 앞에서 말한 대로다. “멜니코프 문턱을 넘었으니 우리 계는 카오스”라는 문장은 절반만 맞다. 정확히는 “카오스적 안장이 생겼고, 그 결과 최종상태 예측이 어려워졌다”이다.
• 적분 안에 들어가는 것은 비섭동 궤도다. 섭동된 계의 수치 궤도를 넣고 계산하는 실수가 종종 나온다. 그렇게 얻은 값은 1차 이론의 항이 아니라 정체불명의 수다.³

9. 관련 문서[편집]

• 호모클리닉 궤도 · 스메일 말굽 · 카오스 이론
• 불변 다양체 · 흡인 유역 · 끌개 · 안장점
• 푸앵카레 단면 · 랴푸노프 지수 · 분기 이론
• 더핑 방정식 · 실니코프 분기 · 극한 순환
• KAM 정리 · 아널드 확산 · 해밀토니안 역학
• 동역학계 · 수치 연속법 · 플러터

10. Footnotes[편집]

1. 여기서 자주 나오는 오해가 “카오스는 공진에서 시작된다”는 직관이다. 쌍우물 더핑의 얽힘 문턱을 최소로 만드는 진동수는 우물 바닥의 소진폭 고유진동수($\sqrt2$)와도, 우물을 넘나드는 큰 진동의 진동수와도 다르다. 문턱을 정하는 것은 $\operatorname{sech}(\pi\omega/2)$ 라는 인자, 즉 호모클리닉 궤도가 안장을 떠났다 돌아오는 데 걸리는 시간 폭의 푸리에 변환이다. 강제가 그 시간 폭과 맞아떨어질 때 에너지가 제일 잘 들어간다는 뜻이라, 공진은 공진인데 선형 공진이 아니다. ↩

2. 그래서 이런 그림을 볼 때는 “곡선이 실험을 얼마나 잘 맞히나”가 아니라 “곡선 아래쪽에서는 정말 아무 일도 안 일어나는가”를 봐야 한다. 문턱 아래에서 카오스가 관측되면 그건 이론의 실패이고, 문턱 위에서 한동안 얌전한 것은 이론의 실패가 아니다. 이 비대칭을 모르고 곡선을 회귀선처럼 다루는 보고서가 은근히 많다. ↩

3. 비슷한 부류로 “닫힌 형태 호모클리닉 해가 없으니 대충 비슷한 함수로 근사했다”도 있다. $\operatorname{sech}$ 계열 함수로 피팅해서 넣으면 적분값이 몇 배씩 틀어진다. 지수적으로 죽는 꼬리의 감쇠율이 그대로 $\operatorname{sech}(\pi\omega/2)$ 인자의 지수로 들어가기 때문에, 꼬리를 대충 맞춘 근사는 문턱을 자릿수 단위로 망친다. ↩