분기 이론

1. 개요[편집]

분기 이론(bifurcation theory)은 동역학계의 매개변수를 매끄럽게 바꿀 때 해의 개수·안정성·위상적 구조가 질적으로 바뀌는 지점을 분류하고 추적하는 이론이다. 여기서 “질적”이 핵심이다. 매개변수를 아무리 바꿔도 그냥 해가 조금씩 움직일 뿐이라면 분기가 아니다. 고정점이 새로 생기거나 사라지거나 안정성을 바꾸거나, 없던 주기궤도가 튀어나올 때 비로소 분기라고 부른다.

$\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x};\mu)$ 에서 분기가 일어나려면 필요조건이 있다. 고정점 $\mathbf{x}^*$ 에서의 야코비안 $D\mathbf{f}$ 가 허수축 위에 고유값을 가져야 한다. 실고유값 하나가 0을 지나면 정상(steady) 분기 계열, 복소켤레 쌍이 $\pm i\omega$ 로 축을 가로지르면 호프 분기다. 이 사실은 공학적으로 무겁다 — 좌굴에서 접선강성행렬이 특이해지는 순간, 플러터에서 감쇠가 0을 통과하는 순간이 전부 이 조건의 다른 얼굴이기 때문이다.

2. 국소 분기 사대장과 정규형[편집]

중심다양체 정리에 따르면, 분기점 근처에서 계의 본질적 거동은 허수축 고유값에 대응하는 저차원 다양체 위로 축약된다. 그 위에서 좌표변환으로 불필요한 항을 다 죽인 최소 형태가 정규형(normal form)이다. 1-매개변수 계에서 여차원 1로 만나는 것은 사실상 네 개뿐이다.

안장-마디 분기(saddle-node, fold): $\dot{x} = \mu + x^{2}$. $\mu<0$ 이면 고정점이 $x = \pm\sqrt{-\mu}$ 두 개(안정 하나, 불안정 하나), $\mu=0$ 에서 충돌해 소멸하고 $\mu>0$ 이면 아무것도 없다. 해가 존재하다가 그냥 없어진다는 점에서 가장 폭력적인 분기다. 극한점(limit point)·턴백 점이라고도 부르며, 구조해석의 스냅스루가 바로 이것.

가로지름 분기(transcritical): $\dot{x} = \mu x - x^{2}$. 고정점 $x=0$ 과 $x=\mu$ 가 $\mu=0$ 에서 만나 안정성을 교환한다. 개체수 모형처럼 $x=0$ 해가 물리적으로 항상 존재해야 하는 계에서 나온다.

갈래 분기(pitchfork): 대칭성 $x \to -x$ 가 있는 계 전용이다.

$$\dot{x} = \mu x - x^{3} \quad (\text{초임계}), \qquad \dot{x} = \mu x + x^{3} \quad (\text{아임계})$$

초임계는 $\mu>0$ 에서 $x = \pm\sqrt{\mu}$ 두 안정해가 부드럽게 자란다. 아임계는 반대로 불안정한 가지가 $\mu<0$ 쪽으로 뻗어 있어, 임계점에서 계가 멀리 튄다. 완전한 오일러 기둥의 좌굴이 초임계 갈래 분기이며, 여기에 초기 결함을 넣으면 대칭이 깨져 갈래가 풀리고(imperfect bifurcation) 매끄러운 곡선 + 떨어져 나온 가지가 된다. 이것이 좌굴 문서의 결함 민감도 이야기와 같은 그림이다.

호프 분기(Hopf/Andronov-Hopf): 복소켤레 고유값 쌍이 허수축을 가로지르며 극한주기궤도가 태어난다. 극좌표 정규형은

$$\dot{r} = \mu r + a r^{3}, \qquad \dot{\theta} = \omega + b r^{2}$$

이고, $a<0$ 이면 초임계(진폭이 $\sqrt{\mu}$ 로 부드럽게 자람), $a>0$ 이면 아임계(작은 불안정 궤도가 있고, 임계점을 넘으면 대진폭 진동으로 점프 + 이력). 공학에서 진동이 “스르륵 커지는” 계와 “어느 날 갑자기 터지는” 계를 가르는 것이 바로 이 3차 계수의 부호다.¹

3. 분기 다이어그램과 수치 연속법[편집]

분기 다이어그램은 가로축에 매개변수 $\mu$, 세로축에 해의 대표값(진폭·노름 등)을 놓고 안정 가지는 실선, 불안정 가지는 점선으로 그린 지도다. 문제는 이 지도를 어떻게 그리느냐다.

가장 순진한 방법은 $\mu$ 를 조금씩 키우며 뉴턴-랩슨법으로 $\mathbf{f}(\mathbf{x};\mu)=0$ 을 푸는 자연 매개변수 연속법이다. 이전 해를 초기추정으로 쓰니 수렴도 빠르다. 그런데 안장-마디 점에서 이 방법은 반드시 죽는다. 그 점에서 $D_x\mathbf{f}$ 가 특이해져 뉴턴 보정이 발산하고, 애초에 그 너머에 해가 없으므로 $\mu$ 를 밀어붙일 대상 자체가 없다.

해법은 $\mu$ 를 미지수로 승격시키고, 대신 해곡선을 따라 잰 호의 길이를 새 매개변수로 쓰는 의사 호길이 연속법(pseudo-arclength continuation, Keller 1977)이다. 방정식 하나를 추가한다.

$$\mathbf{f}(\mathbf{x},\mu) = \mathbf{0}, \qquad \dot{\mathbf{x}}_0^{\mathsf{T}}(\mathbf{x}-\mathbf{x}_0) + \dot{\mu}_0(\mu-\mu_0) - \Delta s = 0$$

확장된 야코비안은 극한점에서도 정칙이므로 곡선을 뒤로 접힌 구간까지 그대로 따라간다. 구조해석에서 후좌굴 경로를 쫓는 호장법이 바로 이 아이디어의 유한요소판이며, 동역학계 쪽에서는 AUTO·MatCont·COCO 같은 패키지가 표준이다. 가지 위에서 야코비안의 고유값을 계속 감시하다가 부호가 바뀌는 지점을 시험함수(test function)의 영점으로 잡아내면 분기점이 자동 검출된다.

4. 주기배가 캐스케이드와 파이겐바움[편집]

이산 사상에서는 고유값이 $-1$ 을 통과하는 주기배가 분기(period doubling)가 추가로 등장한다. 교과서 표본이 로지스틱 사상 $x_{n+1} = r\,x_n(1-x_n)$ 이다.

$r<1$ 이면 0으로 사멸, $1<r<3$ 이면 $x^*=1-1/r$ 로 수렴, $r=3$ 에서 주기 2로 갈라지고, $r_2 \approx 3.4495$ 에서 주기 4, $r_3 \approx 3.5441$ 에서 주기 8 — 이 배가가 기하급수적으로 빨라지며 $r_\infty \approx 3.569946$ 에 쌓인다. 그 너머는 카오스 밴드이고, 중간중간 주기 3·5·6 창(window)이 박혀 있다.

미첼 파이겐바움이 1975년에 발견한 것은 이 간격의 비가 보편상수로 수렴한다는 사실이다.

$$\delta = \lim_{n\to\infty} \frac{r_n - r_{n-1}}{r_{n+1} - r_n} = 4.669201609\ldots$$

로지스틱이든 사인 사상이든, 극값이 이차인 단봉 사상이면 전부 같은 $\delta$ 가 나온다. 상전이의 보편성 부류와 정확히 같은 냄새가 나는데, 실제로도 재규격화 논증으로 설명된다.

5. 랴푸노프 지수 — 카오스의 계량[편집]

분기 다이어그램이 “무엇이 보이나”라면 랴푸노프 지수는 “얼마나 민감한가”를 숫자로 준다. 1차원 사상에서는

$$\lambda(r) = \lim_{N\to\infty} \frac{1}{N}\sum_{n=0}^{N-1} \ln\big|f'(x_n)\big| = \big\langle \ln|r(1-2x)| \big\rangle$$

이며, 이웃한 두 초기조건의 간격이 $e^{\lambda n}$ 으로 벌어진다는 뜻이다. $\lambda<0$ 이면 주기 끌개로 수렴(간격이 줄어듦), $\lambda>0$ 이면 결정론적 카오스, $\lambda=0$ 은 분기점이다. 그래서 분기 다이어그램 위에 $\lambda(r)$ 을 겹쳐 그리면 배가점마다 곡선이 0에 닿았다 내려가고, $r_\infty$ 이후에 0 위로 올라온 뒤 주기 창마다 다시 음수로 처박히는 그림이 나온다.

연속계에서는 지수가 자유도 수만큼 존재하고(랴푸노프 스펙트럼), 최대 지수의 부호로 카오스를 판정한다. 유의할 점: 유한한 적분 구간에서 얻은 $\lambda$ 는 항상 추정치이며, 과도응답을 충분히 버리지 않으면 부호가 뒤집혀 나온다. 부동소수점 연산 오차가 지수적으로 증폭되므로 개별 궤적은 어차피 참값이 아니지만, 그림자 정리(shadowing) 덕에 통계적 성질은 살아남는다는 것이 카오스 수치해석의 위태로운 위안이다.²

6. 공학에서 분기점은 어디에나 있다[편집]

• 좌굴: 완전 기둥은 초임계 갈래 분기, 얇은 셸은 아임계라 스냅스루로 무너진다. 후좌굴 경로 추적이 호장법의 존재 이유. 자세한 건 좌굴 문서.
• 플러터: 구조 감쇠가 공력 감쇠에 잡아먹혀 0을 지나는 호프 분기다. 아임계 호프면 임계속도 아래에서도 큰 교란에 의해 극한주기진동(LCO) 이 촉발될 수 있어서, 선형 플러터 속도만 믿으면 안 된다. 플러터·항공탄성학 참고.
• 연소 불안정: 열방출과 음향 모드의 되먹임(레일리 판정)이 임계를 넘으면 호프 분기로 열음향 진동이 자란다. 로켓 엔진과 가스터빈에서 하드웨어를 실제로 부순다. 연소 시뮬레이션 참고.
• 유동 천이: 원형 실린더 뒤 카르만 와열의 개시($Re \approx 47$)는 정상해가 잃는 호프 분기다. 그 뒤로 분기가 겹겹이 쌓여 결국 난류 천이에 이른다(뤼엘-타켄스 시나리오).
• 다중해: 같은 경계 조건에서 CFD 해가 두 개 이상 존재하고 초기값에 따라 다른 답이 나오는 상황도 분기의 결과다. “왜 케이스마다 답이 다르냐”는 질문의 절반은 수렴 문제가 아니라 지역 최적해처럼 실재하는 다중 가지 문제다.

즉 “어느 값을 넘기면 갑자기 못 쓰게 되는” 모든 공학 문제는 분기점 찾기 문제다. 그 지점을 매개변수 스윕으로 무식하게 훑을 것인지, 연속법으로 곡선을 따라가며 시험함수 영점을 잡을 것인지가 갈릴 뿐이다.³

7. 관련 문서[편집]

• 좌굴 · 호장법 · 비선형 구조해석
• 플러터 · 항공탄성학 · 난류 천이
• 상전이 · 자기조직화 · 셀룰러 오토마타
• 카오스 이론 · 랴푸노프 지수 · 고유값 문제
• 뉴턴-랩슨법 · 카르만 와열 · 계산물리

8. Footnotes[편집]

1. 그래서 실험 데이터에서 “진동 진폭이 속도의 제곱근으로 자란다”가 보이면 초임계 호프를 의심하고, “임계속도를 넘자마자 대진폭이 튀어나오고 속도를 낮춰도 안 꺼진다”면 아임계 호프 + 이력을 의심하면 된다. 물론 현장에서는 이 관찰을 할 여유가 있기 전에 시편이 먼저 부러진다. ↩

2. 그림자 정리는 “수치 궤적 근처에 정확한 궤적이 하나 존재한다”는 것이지 “네 궤적이 맞다”가 아니다. 로렌츠 계에서 며칠짜리 적분 결과를 궤적 단위로 신뢰하는 건 여전히 도박이며, 이래서 수치기상예보는 단일 예보 대신 앙상블 예보로 간다. ↩

3. 실무 CAE에서 하중을 100 스텝으로 쪼개 돌리다가 87 스텝에서 “수렴 실패”를 보고 스텝을 더 잘게 쪼개는 사람을 자주 본다. 그게 극한점이면 스텝을 아무리 쪼개도 영원히 안 넘어간다. 하중 제어를 버리고 호장법으로 갈아타는 게 답이다. ↩