안장-마디 분기

수치해석 계산물리 물리 마지막 수정: 2026-08-09 04:37:05

1. 개요[편집]

안장-마디 분기 Saddle-node bifurcation
다른 이름	접힘(fold) · 극한점(limit point) · 턴백점
여차원	1 (매개변수 하나로 일반적으로 만남)
정규형	ẋ = μ − x²
조건	비퇴화 f_xx ≠ 0 · 횡단성 f_μ ≠ 0
스케일링	가지 x ∝ √μ · 병목 통과시간 ∝ 1/√μ
공학적 얼굴	스냅스루 · 하중 극한점 · 착화/소화

안장-마디 분기(saddle-node bifurcation)는 매개변수를 바꿀 때 안정 평형점과 불안정 평형점이 서로 다가와 충돌한 뒤 함께 사라지는 분기다. 접힘(fold), 극한점(limit point), 턴백점 등 분야마다 이름이 다르지만 전부 같은 사건을 가리킨다. 분기 이론의 여차원 1 분기 네 가지 중 해 자체가 소멸한다는 점에서 가장 폭력적이며, 그래서 공학 사고 보고서에 제일 자주 등장한다.

다른 분기들은 해가 남아 있다. 가로지름 분기는 안정성을 교환할 뿐이고, 갈래 분기는 가지가 늘어나며, 호프 분기는 진동으로 갈아탄다. 안장-마디만이 갈아탈 곳을 남기지 않는다. 그 순간 계는 근처에 아무 평형점도 없는 상태가 되어 어딘가 먼 다른 끌개로 튄다. 셸이 뒤집히고, 반응기가 착화하고, 하중 제어 FEM이 “수렴 실패”를 뱉는 것이 전부 이 장면이다.

2. 정규형과 성립 조건[편집]

1차원 표준형은 다음과 같다.

\dot x = \mu - x^2

$\mu>0$ 이면 평형점이 $x = \pm\sqrt{\mu}$ 두 개다. $f'(x) = -2x$ 이므로 $x=+\sqrt\mu$ 는 안정(마디), $x=-\sqrt\mu$ 는 불안정. $\mu=0$ 에서 둘이 원점에서 충돌하고 $\mu<0$ 이면 실 평형점이 없다. 부호 관습이 책마다 달라 $\dot x = \mu + x^2$ 로 쓰기도 하는데, $(x,\mu) \to (-x,-\mu)$ 치환으로 서로 옮겨 가는 같은 정규형이다.¹

일반적인 $\dot x = f(x,\mu)$ 가 $(x_0,\mu_0)$ 에서 안장-마디 분기를 하려면 다음이 필요하다.

f = 0, \quad f_x = 0 \quad (\text{분기 조건}); \qquad f_{xx} \neq 0 \quad (\text{비퇴화}), \qquad f_\mu \neq 0 \quad (\text{횡단성})

이 조건들 아래 음함수 정리를 $\mu$ 에 대해 적용하면 해 가지가 $\mu$ 의 그래프가 아니라 $x$ 의 그래프로 매끄럽게 풀리고, 국소적으로

\mu - \mu_0 \simeq -\frac{f_{xx}}{2 f_\mu}\,(x-x_0)^2 \quad \Longleftrightarrow \quad x - x_0 \simeq \pm\sqrt{-\frac{2f_\mu}{f_{xx}}\,(\mu-\mu_0)}

즉 제곱근 스케일링이 나온다. 분기점에 다가갈수록 두 해가 $\sqrt{|\mu-\mu_0|}$ 로 가까워지고, 가지의 접선이 수직이 된다는 뜻이다. 실험 데이터에서 두 상태가 제곱근으로 붙으며 사라지면 안장-마디를 의심하면 된다.

$n$ 차원에서는 야코비안 $A$ 가 단순 영고유값 하나를 갖고 나머지는 허수축 밖에 있어야 한다. $Aq = 0$ , $p^{\mathsf T}A = 0$ , $p^{\mathsf T}q = 1$ 로 좌·우 영벡터를 잡으면 조건이 이렇게 번역된다.

p^{\mathsf T} B(q,q) \neq 0 \ (\text{비퇴화}), \qquad p^{\mathsf T} f_\mu \neq 0 \ (\text{횡단성})

여기서 $B$ 는 2계 미분의 이중선형 형식이다. 중심 다양체 정리로 1차원으로 내린 뒤 정규형 이론으로 계수를 뽑으면 앞의 스칼라 식이 그대로 재생된다. 나머지 방향은 하트만-그로브만 정리가 표준 안장으로 규격화해 주므로 신경 쓸 필요가 없다.

3. 유령, 병목, 그리고 SNIC[편집]

평형점이 사라진 직후가 오히려 흥미롭다. $\dot x = \mu + x^2$ 에서 $\mu$ 가 작은 양수면 평형점은 없지만 원점 근처에서 $\dot x$ 가 아주 작아 궤도가 오래 머문다. 통과 시간은 정확히 계산된다.

T = \int_{-\infty}^{\infty} \frac{dx}{\mu + x^2} = \frac{\pi}{\sqrt{\mu}}

이 $1/\sqrt\mu$ 발산이 병목(bottleneck) 또는 유령(ghost)이라 불리는 현상이다. 사라진 평형점이 유령처럼 남아 궤도를 붙잡는다.

이 유령이 원형 궤도 위에 놓이면 불변원 위의 안장-마디 분기(saddle-node on an invariant circle, SNIC)가 된다. 안정 평형점과 안장이 불변원 위에서 만나 소멸하면, 그 즉시 원 전체가 하나의 주기궤도가 된다. 태어난 극한 순환의 주기는 병목이 지배하므로

T \propto \frac{1}{\sqrt{\mu - \mu_c}}

로, 임계점에서 주기가 무한대로 발산한다. 진폭은 처음부터 유한한데 진동수만 0에서 연속적으로 자라는 것이 호프 분기와 정반대다. 신경 모형에서 이 성질이 I형 흥분성(임의로 낮은 발화율이 가능)의 정체이고, 적분-발화 모형의 세타 뉴런 표준형이 정확히 SNIC 정규형이다. 흥분성 매질과 피츠휴-나구모 모형에서 문턱값 거동을 만드는 것도 근처의 안장과 그 불변 다양체다.

4. 이산 사상에서의 접힘[편집]

사상 $x_{n+1} = g(x_n,\mu)$ 에서도 같은 분기가 있다. 조건은 고정점의 승수가 $+1$ 을 가로지르는 것이고, 이 경우 접선 분기(tangent bifurcation)라고도 부른다. 이름 그대로 $g$ 의 그래프가 대각선 $y=x$ 에 접했다가 떨어지는 그림이기 때문이다.

가장 유명한 사례가 로지스틱 사상의 주기 3 창이다. $f^3$ 의 그래프가 대각선에 접하는 지점

r = 1 + 2\sqrt{2} \approx 3.828427

에서 안정 주기 3 궤도와 불안정 주기 3 궤도가 한꺼번에 태어나고, 카오스 밴드 한복판에 주기 창이 열린다. $r$ 을 이 값보다 아주 조금 작게 두면 궤도가 거의 주기 3처럼 보이다가 이따금 카오스로 터지는 I형 간헐성(type-I intermittency)이 나타나는데, 이때 규칙적 구간의 평균 길이가 $(r_c-r)^{-1/2}$ 로 발산한다. 앞 절에서 본 병목 통과 시간 $\pi/\sqrt\mu$ 와 정확히 같은 지수이며, 실제로 같은 정규형에서 나온 결과다. 간헐성·카오스 이론 참고.

5. 연속법이 죽는 자리[편집]

수치적으로 안장-마디가 악명 높은 이유는 거기서 뉴턴법이 원리적으로 실패하기 때문이다. $f_x$ 가 특이해지므로 뉴턴-랩슨법의 보정량이 폭발하고, 매개변수를 조금만 더 밀면 애초에 해가 없어서 어떤 초기추정으로도 수렴하지 않는다. 하중을 100단계로 쪼개 돌리다 87단계에서 죽었을 때 단계를 더 잘게 쪼개는 것이 헛수고인 이유가 이것이다. 분할을 아무리 늘려도 없는 해로는 수렴하지 못한다.

해법은 $\mu$ 를 미지수로 올리고 방정식 하나를 보태는 것, 즉 수치 연속법의 의사 호길이 방식이다. 핵심은 확장(테두리) 야코비안이

\begin{pmatrix} f_x & f_\mu \\ \dot{x}_0^{\mathsf T} & \dot{\mu}_0 \end{pmatrix}

접힘점에서도 정칙이라는 사실이다. $f_x$ 가 특이해도 횡단성 $f_\mu \neq 0$ 덕분에 마지막 열이 결손 랭크를 메워 준다. 그래서 확장된 계에 대한 뉴턴법은 접힘점을 아무 일도 없다는 듯 통과하고, 곡선은 뒤로 접힌 채 계속 이어진다. 구조해석의 호장법(Riks–Crisfield)이 같은 아이디어의 유한요소판이다.

부수적으로 얻는 진단 기준도 있다. 접힘점과 갈래(분기)점은 둘 다 $f_x$ 를 특이하게 만들지만, 접힘점에서는 $f_\mu$ 가 $f_x$ 의 상공간 밖에 있고 갈래점에서는 안에 있다. 연속법 패키지가 LP(limit point)와 BP(branch point)를 구분해 찍어 주는 근거가 이 한 줄이다.² 검출 자체는 가지를 따라가며 시험함수(보통 테두리 계로 정의한 매끄러운 스칼라)의 부호 변화를 보는 방식으로 한다.

6. 공학에 박혀 있는 접힘점들[편집]

스냅스루 좌굴. 얕은 아치나 얕은 셸에 하중을 걸면 하중-변위 곡선이 최대점을 지나 아래로 꺾인다. 그 최대점이 안장-마디이고, 하중 제어 실험에서는 그 순간 구조가 뒤집힌 형상으로 튄다(스냅스루). 변위 제어나 호장법으로는 아래로 꺾인 구간까지 정적으로 추적할 수 있다. 좌굴·비선형 구조해석 참고.
하중 극한점. 접선강성행렬 $K_T$ 의 행렬식이 0이 되는 점이다. 완전 기둥의 갈래 분기와 달리 극한점은 결함이 있는 실제 구조에서 훨씬 흔하게 나타난다. 초기 결함이 갈래 분기를 풀어 놓으면 결국 남는 것이 접힘점이기 때문.
반응기의 착화와 소화. 연속 교반 탱크 반응기(연속 교반 탱크 반응기)의 정상상태 온도를 체류시간에 대해 그리면 S자 곡선이 나오고, 위아래 두 접힘점 사이에서 저온해·중간해·고온해 세 개가 공존한다. 체류시간을 늘리다 아래쪽 접힘을 지나면 착화(고온 가지로 점프), 줄이다 위쪽 접힘을 지나면 소화가 일어나며, 두 사건이 다른 값에서 일어나므로 이력(hysteresis)이 생긴다. 열발생이 아레니우스 방정식의 지수 함수라서 생기는 구조이고, 무차원 정리에서는 담쾰러 수가 분기 매개변수 역할을 한다. 반응속도론·연소 시뮬레이션 참고.
극한주기궤도의 접힘. 평형점만의 이야기가 아니다. 안정 극한궤도와 불안정 극한궤도가 충돌해 소멸하는 주기궤도의 접힘(LPC, fold of cycles)도 같은 정규형을 따른다. 아임계 호프 분기 뒤에 이 접힘이 붙어 있어서 플러터의 극한주기진동이 임계속도 아래에서도 유지되는 이력 고리가 만들어진다.
다중 정상해. 같은 경계 조건에서 CFD 정상해가 두 개 이상 나오는 상황(확대관 유동의 대칭 파괴 후 재대칭화, 자연대류 방향 반전 등)도 접힘점 사이 구간에 들어와 있다는 뜻이다. 초기장에 따라 다른 답이 나오는 것은 버그가 아니라 계의 성질이다.³

7. 관련 문서[편집]

그래서 논문마다 부호가 달라도 싸울 필요가 없다. 정말 싸울 일은 “이게 접힘이냐 갈래냐”이고, 그건 부호가 아니라 $f_\mu$ 가 0인지 아닌지로 갈린다. ↩
실무에서는 이걸 몰라서 대칭성 있는 모형의 갈래점을 접힘점으로 오인하고 “왜 호장법으로 넘어갔는데 다른 가지가 안 나오죠”라고 묻는 일이 흔하다. 접힘은 넘어가면 되지만 갈래는 가지 전환(branch switching)을 따로 해 줘야 한다. ↩
이 답을 상사에게 설명하는 것이 제일 어렵다. “물리적으로 답이 두 개입니다”는 대체로 “그래서 어느 게 맞냐”는 재질문으로 돌아오는데, 정답은 “둘 다 맞고, 어디로 갈지는 초기조건과 교란이 정한다”이다. 이 대답이 통하는 회사는 드물다. ↩