기하 수치적분

수치해석 계산물리 마지막 수정: 2026-08-03 04:12:41

1. 개요[편집]

기하 수치적분 Geometric Numerical Integration
설계 목표	흐름의 기하 구조를 이산 사상이 정확히 물려받게 하기
보존 대상	심플렉틱 형식 · 1차 적분 · 다양체 제약 · 시간 가역성 · 리 군 구조
이론적 엔진	후진 오차 해석(수정 방정식)
대표 기법	심플렉틱·변분 적분기, RKMK, RATTLE, 이산 그래디언트
대표 한계	구조를 전부 동시에 지킬 수는 없다 (게-마스든)

한 걸음의 오차를 줄이는 게 수치해석이라면, 백만 걸음 뒤에도 무대가 무너지지 않게 하는 게 기하 수치적분이다.

기하 수치적분(geometric numerical integration)은 미분방정식의 정확한 흐름이 가진 기하학적 구조 — 심플렉틱 형식, 위상공간 부피, 에너지·운동량 같은 1차 적분, 다양체 제약, 시간 가역성, 리 군 구조 — 를 이산 사상 $\Phi_h$ 가 오차항 없이 그대로 물려받도록 적분기를 설계하는 분야다. “해를 얼마나 정확히 근사하는가” 대신 “무엇을 정확히 보존하는가”를 설계 변수로 삼는 것이 핵심 전환이다.

이게 왜 필요하냐면, 국소 정확도와 장기 정성적 거동이 서로 다른 물건이기 때문이다. 4차 룽게-쿠타법은 한 스텝 오차가 $O(h^5)$ 로 아름답지만, 분자동역학이나 천체역학처럼 $10^8$ 스텝을 도는 계산에서는 에너지가 한 방향으로 슬금슬금 흘러 계가 혼자 뜨거워지거나 궤도가 나선으로 감긴다. 반면 2차밖에 안 되는 베를레 적분은 같은 시간 동안 에너지를 작은 폭으로 진동만 시킨다.¹ 이 문서는 그런 “구조 보존” 설계 원리 전체의 상위 허브이고, 개별 기법의 세부는 심플렉틱 적분기·변분 적분기·리 군 문서에 위임한다.

2. 보존할 구조의 목록[편집]

기하 적분기 설계는 사실상 “이 계가 가진 구조 중 무엇을 정확히 지킬 것인가”를 고르는 일이다.

구조	대상 계	대표 기법
심플렉틱 형식 $\omega$	해밀토니안 역학계	리프프로그, 가우스-르장드르 RK
변분 원리	라그랑지안계	변분 적분기
1차 적분(에너지 등)	임의의 보존계	투영법, 이산 그래디언트
다양체 제약 $g(q)=0$	구속계, 강체	SHAKE / RATTLE
시간 가역성	$\rho$ -가역계	대칭 적분기
리 군 / 등질공간	자세, 스핀, 게이지장	RKMK, 마그누스, 크라우치-그로스만

왜 이게 통하는가의 답은 후진 오차 해석에 있다. 차수 $p$ 인 적분기의 이산 궤적은 원래 벡터장이 아니라 수정된 벡터장의 정확한 궤적 위에 (거의) 놓이는데, 적분기가 구조를 보존하면 이 수정 벡터장도 같은 구조를 갖는다. 즉 우리가 실제로 푸는 문제가 “원래 문제와 같은 종류의 문제”가 된다. 구조를 안 지키면 수정 벡터장에 소산·발산 항이 붙고, 그게 장기 표류(secular drift)의 정체다.

3. 심플렉틱·변분 계열[편집]

가장 잘 정비된 갈래. 해밀턴 방정식 $\dot z = J\nabla H(z)$ 의 흐름은 야코비안이 $M^\top J M = J$ 를 만족하고, 이 조건을 이산 사상에 그대로 요구한 것이 심플렉틱 적분기다. 가분 해밀토니안 $H=T(p)+V(q)$ 에서는 드리프트-킥 스플리팅으로 공짜에 가깝게 얻어지고, 아니면 음함수 계열(중점법, 가우스-르장드르 RK)을 쓴다.

같은 결과에 라그랑지안 쪽에서 도달하는 것이 변분 적분기다. 작용 적분을 먼저 구적으로 이산화한 뒤 그 이산 작용에 변분 원리를 적용하면, 나오는 사상이 자동으로 심플렉틱이고 게다가 뇌터 정리의 이산판까지 딸려 온다 — 이산 작용이 갖는 대칭에 대응하는 운동량이 기계 정밀도로 보존된다. 회전 대칭이 있으면 각운동량이 스텝 수와 무관하게 지켜지는 식이다.

4. 에너지를 정확히 보존하고 싶다면[편집]

심플렉틱 적분기는 에너지를 정확히 보존하지는 않는다(유계 진동만 한다). 이 둘은 자주 뭉뚱그려지지만 완전히 다른 성질이다.² 에너지 자체를 기계 정밀도로 지켜야 하는 상황도 있는데, 이때 쓰는 것이 두 갈래다.

투영법은 평범한 적분기로 한 스텝 간 뒤 결과를 보존량 준위집합 $\{H(z)=H_0\}$ 위로 라그랑주 승수법을 써서 정사영한다. 구현이 쉽지만 심플렉틱성이 깨지고 위상 오차가 오히려 커지는 부작용이 있다.

**이산 그래디언트 방법**은 좀 더 우아하다. $\bar\nabla H(z,\hat z)$ 를

\bar\nabla H(z,\hat z)^\top (\hat z - z) = H(\hat z) - H(z)

를 만족하도록 정의하고 $(\hat z - z)/h = J\,\bar\nabla H(z,\hat z)$ 로 적분하면, 위 정의식과 $J$ 의 반대칭성만으로 $H(\hat z) = H(z)$ 가 항등적으로 따라 나온다. 대표적인 선택이 평균 벡터장(AVF) 방법 $\bar\nabla H = \int_0^1 \nabla H(z+\xi(\hat z - z))\,\mathrm{d}\xi$ 다. 대신 음함수라 매 스텝 비선형 연립을 풀어야 하고, 심플렉틱하지 않다.

여기서 유명한 벽이 나온다. 게-마스든 정리(Ge Zhong & Marsden, 1988): 고정 스텝 적분기가 심플렉틱이면서 $H$ 까지 정확히 보존하면, 그것은 시간 재매개변수화를 제외하고 정확한 흐름 그 자체다.³ 즉 비적분계에서 “심플렉틱 + 에너지 정확 보존”은 고정 스텝으로 원리적으로 불가능하다. 보존량은 공짜가 아니고, 심지어 골라야 한다.

5. 가역 적분기[편집]

벡터장이 어떤 대합 $\rho$ 에 대해 $\rho f(y) = -f(\rho y)$ 를 만족하면 그 계는 $\rho$ -가역이다(속도 반전 대칭이 있는 역학계 대부분이 여기 해당). 이때 대칭 적분기, 즉 $\Phi_h^{-1} = \Phi_{-h}$ 인 적분기를 쓰면 이산 사상도 가역이 되고, 차수가 자동으로 짝수가 된다는 보너스까지 붙는다. 리프프로그가 심플렉틱이면서 동시에 대칭인 것이 그 인기의 절반쯤을 설명한다.

중요한 점은 가역 KAM 이론이 심플렉틱 KAM 이론과 거의 평행하게 성립한다는 것이다. 그래서 해밀토니안이 아니지만 가역인 계에서는 대칭 적분기가 심플렉틱 적분기와 비슷한 장기 안정성을 준다. KAM 정리와 카오스 이론 쪽 수치실험이 이 성질에 크게 기대고 있다.

6. 다양체 위에서 — 리 군과 제약[편집]

해가 선형공간이 아니라 곡면 위에 살면, 덧셈으로 만든 스텝은 정의상 곡면을 벗어난다. 회전행렬을 $R + h R[\omega]_\times$ 로 갱신하면 $R^\top R$ 이 $I$ 에서 멀어지고, 코드에는 “주기적 재직교화”라는 대증요법이 생긴다.⁴

마그누스 전개. 선형 문제 $Y' = A(t)Y$ 의 해를 $Y(t) = \exp(\Omega(t))Y_0$ 로 두고 $\Omega$ 를 교환자 급수로 전개한다.

\Omega(t) = \int_0^t A(s)\,\mathrm{d}s + \frac{1}{2}\int_0^t\!\!\int_0^{s_1} [A(s_1), A(s_2)]\,\mathrm{d}s_2\,\mathrm{d}s_1 + \cdots

$A$ 가 리 대수 값이면 $\Omega$ 도 리 대수 값이므로 $\exp(\Omega)$ 는 항상 군 원소다. 절단해도 제약 위반이 0이라는 뜻.

RKMK(문테-카스). 비선형 경우로의 일반화. 스텝을 군 위가 아니라 리 대수(선형공간) 위에서 밟는다. $\mathrm{dexp}^{-1}$ 방정식을 세우고 거기에 평범한 룽게-쿠타를 그대로 적용한 뒤, 마지막에 지수사상으로 군에 되돌린다. 고전 RK의 계수표를 그대로 재활용할 수 있다는 게 실무적 장점이다. 자세한 지수/로그 사상 얘기는 리 군 참고.
SHAKE / RATTLE. 홀로노믹 제약 $g(q)=0$ 이 붙은 계에서는 승수 $\lambda$ 를 미지수로 넣어 매 스텝 제약을 정확히 만족시킨다. SHAKE는 위치 제약만 강제해 속도가 제약 접평면을 벗어나는 반면, RATTLE은 숨은 제약 $\nabla g(q)^\top M^{-1}p = 0$ 까지 두 번째 사영으로 강제하며, 제약 다양체 위에서 심플렉틱이다. GROMACS 같은 MD 코드가 결합 길이를 고정해 시간간격을 늘리는 트릭이 정확히 이것이다. 제약 방정식을 그대로 두면 미분대수방정식 문제가 되어 인덱스 축소 얘기로 넘어간다.

7. 한계와 오해[편집]

구조 보존이 정확도를 대체하지 않는다. 심플렉틱 적분기라도 개별 궤도의 위치 오차는 여전히 $O(h^p)$ 로 자라고, 카오스계에서는 지수적으로 갈라진다. 보존되는 건 통계적·정성적 그림이지 궤적 자체가 아니다.
가변 시간간격과 상극이다. 스텝마다 $h$ 가 바뀌면 보존하는 수정 해밀토니안도 바뀌어 장기 이득이 사라진다.
소산계에는 의미가 줄어든다. 서모스탯이나 랑주뱅 동역학을 붙이는 순간 계는 해밀토니안이 아니며, 이때는 목표가 “구조 보존”이 아니라 “올바른 불변측도 표본추출”로 바뀐다. 그래도 해밀토니안 몬테카를로처럼 심플렉틱성과 가역성이 통계적 정확성의 필수 조건이 되는 하이브리드도 있다.
강성(stiff) 진동 문제에서는 구조 보존과 큰 스텝을 동시에 노리는 지수 적분기 계열이 별도 갈래로 발전했다.

8. 관련 문서[편집]

이 분야의 사실상 표준 교과서는 Hairer, Lubich, Wanner, Geometric Numerical Integration (2nd ed., 2006). 표지에 케플러 궤도 그림이 있고, 1장부터 “RK4로 태양계를 돌리면 이렇게 됩니다” 하고 망한 그림을 보여주며 시작한다. 설득력이 상당하다. ↩
그래서 “이 코드는 에너지 보존형입니다”라는 홍보 문구를 보면 반사적으로 되물어야 한다. 심플렉틱인가, 에너지 보존인가, 아니면 그냥 잔차가 작다는 뜻인가. 셋은 전부 다른 말이다. ↩
증명의 취지를 한 줄로 옮기면 “그 정도로 완벽한 적분기가 있다면 그건 미분방정식을 이미 풀어놓은 것”이다. 공짜 점심은 없고, 심지어 메뉴판에서 하나만 골라야 한다. ↩
재직교화는 잘못이 아니라 증상이다. 매 스텝 그람-슈미트를 돌리고 있다면, 그건 적분기가 다양체를 모르고 있다는 신호로 읽는 게 맞다. ↩