연산자 분리

수치해석 시뮬레이션 마지막 수정: 2026-08-01 04:44:52

1. 개요[편집]

연산자 분리(operator splitting)는 여러 물리 과정이 합쳐진 시간발전 방정식 $\partial_t u = (A + B)u$ 를 $A$ 만 있는 부분문제와 $B$ 만 있는 부분문제로 쪼갠 뒤, 각각에 최적인 적분기를 번갈아 적용해 한 스텝을 진행하는 시간 적분 전략이다. 이류-확산-반응처럼 성격이 전혀 다른 항이 섞인 문제에서, 통짜로 짠 하나의 거대한 솔버 대신 검증된 전용 솔버들을 조립해 쓸 수 있게 해 준다.

범위 정리: 이 문서는 분할 자체의 차수·오차 구조·한계를 다룬다. 단일 방정식에 대한 $\theta$ -도식의 안정성은 크랭크-니콜슨법, 이류·확산 항의 물리와 차분 도식 선택은 대류-확산 방정식, 압력-속도 연성이라는 특수한 분할은 SIMPLE 알고리즘 문서에 있다.

2. 리 트로터 분할과 스트랑 분할[편집]

선형 문제 $u' = (A+B)u$ 의 정확해는 $u(t+\Delta t) = e^{(A+B)\Delta t}u(t)$ 다. 여기서 지수를 그냥 쪼개는 것이 가장 단순한 분할이다.

\Phi^{\text{LT}}_{\Delta t} = e^{A \Delta t} e^{B \Delta t}

리 트로터 분할(Lie–Trotter splitting)로, 국소 절단오차가 $O(\Delta t^2)$ 이고 전역 정확도는 1차다. 절차로 읽으면 “한 스텝 동안 $A$ 만 다 진행시키고, 그 결과를 초기값 삼아 $B$ 만 다 진행시킨다”이다.

대칭으로 감싸면 차수가 하나 올라간다.

\Phi^{\text{S}}_{\Delta t} = e^{A \Delta t/2}\, e^{B \Delta t}\, e^{A \Delta t/2}

이것이 스트랑 분할(Strang splitting, 1968)이고 전역 2차다. 계산비는 사실상 리 트로터와 같다 — 연속된 스텝의 반쪽 $A$ 스텝을 합치면( $e^{A\Delta t/2}e^{A\Delta t/2} = e^{A\Delta t}$ ) 스텝당 $A$ 한 번, $B$ 한 번으로 돌아가기 때문이다. 결국 첫 스텝에서 반 스텝 앞서 출발하고 마지막에 반 스텝 보정하는 것이 전부다. 공짜로 1차가 2차가 되므로, 스트랑을 안 쓸 이유가 거의 없다.

이 대칭 구조는 다른 데서도 익숙하다. 심플렉틱 적분기의 립프로그(=속도 베를레)는 해밀토니안을 $T(p) + V(q)$ 로 나눈 것에 스트랑 분할을 적용한 바로 그 사상이며, “킥-드리프트-킥”이 곧 $e^{V\Delta t/2}e^{T\Delta t}e^{V\Delta t/2}$ 다.

3. 분할오차는 교환자에 비례한다[편집]

왜 쪼개면 틀리는가에 대한 답은 명확하다. 교환자 $[A,B] = AB - BA$ 때문이다. 베이커-캠벨-하우스도르프(BCH) 공식을 쓰면

e^{A\Delta t} e^{B\Delta t} = \exp\!\left( (A+B)\Delta t + \tfrac{1}{2}[A,B]\Delta t^2 + O(\Delta t^3) \right)

이므로, 리 트로터의 선행 오차항이 정확히 $\frac{1}{2}[A,B]\Delta t^2$ 이다. 스트랑은 대칭성 덕에 $\Delta t^2$ 항이 상쇄되고 선행 오차가 $\frac{1}{24}([B,[B,A]] - \frac{1}{2}[A,[A,B]])\Delta t^3$ 꼴의 이중 교환자로 밀린다. 결론은 세 줄로 요약된다.

$[A,B] = 0$ 이면 분할오차는 0이다. 서로 교환하는 연산자는 아무렇게나 쪼개도 정확하다.
오차 크기는 $\Delta t$ 가 아니라 두 물리 과정이 얼마나 서로 간섭하는가가 결정한다. 반응이 국소적이고 확산이 느리면 교환자가 작아 분할이 잘 먹힌다.
강한 반응이 급격한 구배를 만들고 그 구배를 확산이 즉시 되먹이는 상황(화염면, 예리한 반응 전선)에서는 교환자가 커져 $\Delta t$ 를 물리 시간척도까지 줄여야 한다. 흔히 말하는 “분할오차 때문에 화염 두께를 못 맞춘다”가 이 얘기다.¹

또 하나 흔히 놓치는 것 — 분할오차는 정상 상태에서도 사라지지 않는다. 시간 미분이 0이 되어도 분할 순서가 만드는 편향은 남을 수 있어, 정상해를 얻으려는 계산에서 $\Delta t$ 를 키우면 엉뚱한 평형점으로 수렴하기도 한다.

4. 3차 이상의 벽 — 음수 시간 스텝[편집]

“스트랑으로 2차가 됐으니 더 합성하면 4차도 되겠지”는 반쯤만 맞다. 대칭 합성으로 고차 분할을 만들 수 있는 건 사실이지만, 3차 이상의 분할에서는 계수 중에 반드시 음수가 나온다. 셍(Sheng, 1989)과 스즈키(Suzuki, 1991)가 각각 증명한 결과로, 2차를 넘는 분할에서 모든 계수를 양수로 두는 것은 원리적으로 불가능하다.

이게 왜 치명적이냐면, 음수 계수는 시간을 거꾸로 감는 부분 스텝을 뜻하기 때문이다. 해밀토니안 역학처럼 시간 가역인 계에서는 아무 문제가 없어서 요시다 4차 합성이 천체역학 코드의 표준으로 자리잡았다. 반면 **확산은 시간 역전이 심각한 비적정 문제(ill-posed)**다 — 역방향 열방정식은 고주파 성분을 지수적으로 증폭시켜 즉시 폭발한다. 그래서 확산·소산이 들어간 문제에서는 스트랑 2차가 사실상 실용적 상한이며, 더 높은 차수가 필요하면 분할 자체를 포기하고 IMEX 룽게-쿠타 같은 결합 도식으로 넘어간다.²

5. 비선형 문제와 분할 순서[편집]

$e^{A\Delta t}$ 라는 표기는 선형 문제에서만 뜻이 있지만, 분할 자체는 비선형에서도 그대로 쓴다. 지수 대신 각 부분문제의 정확한 흐름(flow map) $\varphi^A_{\Delta t}$ 을 놓고 $\varphi^A_{\Delta t/2}\circ\varphi^B_{\Delta t}\circ\varphi^A_{\Delta t/2}$ 로 합성하면 되며, 차수 결과도 유지된다(리 미분의 교환자로 BCH를 다시 쓰면 된다). 유한체적 쪽에서는 같은 구조를 고두노프 분할이라 부르기도 한다. 실제 코드에서는 부분 흐름을 정확히 못 푸니 각각을 수치 적분기로 대신하는데, 이때 부분 솔버의 차수가 분할 차수보다 낮으면 전체가 거기에 끌려 내려간다 — 스트랑으로 감싸 놓고 안쪽 반응 적분을 오일러법으로 때리면 결국 1차다.

분할에는 이 곱셈형(multiplicative) 말고 덧셈형(additive, 각 부분문제를 같은 초기값에서 병렬로 풀고 증분을 더함)도 있다. 병렬화는 쉽지만 정확도와 안정성이 대체로 나빠 주류는 아니다.

순서 선택도 공짜가 아니다. $A$ - $B$ 와 $B$ - $A$ 는 오차의 부호가 반대인 $\pm\frac{1}{2}[A,B]\Delta t^2$ 를 가지므로, 두 순서를 스텝마다 번갈아 쓰면 오차가 상당 부분 상쇄되어 값싼 2차 효과를 얻는다. 또 가장 강성이 큰 항을 스트랑의 안쪽(전체 스텝)에 두는가 바깥쪽(반 스텝)에 두는가에 따라 상수항이 몇 배씩 달라지는데, 일반적으로 빠른 과정(반응)을 안쪽에 넣는 배치가 안정적이다.

6. 경계조건이 차수를 깎아 먹는다[편집]

이론상 2차인 스트랑 분할이 실제 코드에서 1차로 주저앉는 가장 흔한 원인은 경계조건 처리다. 중간 단계의 값 $u^*$ 는 실제 시각의 해가 아니라 부분 흐름을 거친 중간 상태인데, 여기에 원래 문제의 경계값을 그대로 먹이면 경계 근방에 $O(\Delta t)$ 오차층이 생긴다. 비제차(nonhomogeneous) 디리클레 조건이나 시간 의존 경계값이 있을 때 특히 두드러진다. 처방은 중간 단계용 경계값을 별도로 유도해 넣는 것인데, 유도가 지저분해서 실무에서는 “왜 수렴 차수 그래프가 1.2쯤에서 멈추지”라는 미제 사건으로 남는 경우가 많다.³ 격자 수렴 검증을 할 때 격자 수렴 지수 계산이 이상하다면 여기부터 의심할 만하다.

7. 실전에서의 쓰임새[편집]

강성 반응항의 암시적 처리. 연소·대기화학·연소 후처리에서 화학 소스항의 시간척도는 유동보다 수십 배 짧다. 전체를 암시적으로 풀면 거대한 비선형 연립계가 되지만, 분할하면 이류·확산은 명시적으로 싸게 진행하고 반응만 셀별로 독립적인 작은 상미분방정식계로 떼어 CVODE 같은 강성 방정식 전용 적분기에 던질 수 있다. 셀마다 완전히 독립이라 병렬 컴퓨팅 효율도 좋다. 반응속도론 참고.
ADI 법(교대방향 암시법). 2차원 확산을 $x$ 방향 암시 + $y$ 방향 암시로 쪼개면, 각 반 스텝이 삼중대각 행렬 풀기가 되어 LU 분해 없이 $O(N)$ 에 끝난다. 2차원 문제를 1차원 문제 다발로 낮추는 고전적 분할이며, 사실상 스트랑 분할의 방향 버전이다.
분수 스텝법(fractional step). 비압축성 나비에-스토크스 방정식에서 예측 속도를 구한 뒤 압력 포아송 방정식으로 발산을 걷어내는 투영법이, 이류·점성 항과 압력 제약을 분리한 연산자 분리다. 이때도 압력 경계조건의 차수 손실이 그대로 재현된다.
다물리 연성. 열-구조, 유체-구조처럼 서로 다른 솔버를 번갈아 호출하는 분리형(partitioned) 접근 자체가 연산자 분리다. 유체-구조 연성의 약연성 방식이 리 트로터, 부분 반복을 넣은 강연성이 그 보정에 해당한다.

8. 관련 문서[편집]

그래서 연소 코드에서 $\Delta t$ 를 줄였는데도 화염 속도가 안 맞으면, 반응 적분기의 허용오차부터 조이는 대신 분할 순서(반응 먼저 vs 수송 먼저)를 바꿔 보는 게 진단으로 더 빠를 때가 있다. 두 결과의 차이가 곧 분할오차의 크기다. ↩
역방향 열방정식이 얼마나 나쁜지는 푸리에 모드로 보면 명확하다. 파수 $k$ 성분이 $e^{+\nu k^2 |\Delta t|}$ 로 커지므로, 격자를 촘촘히 할수록(= $k$ 가 커질수록) 더 빨리 터진다. 격자를 늘려서 상황이 나빠지는 몇 안 되는 사례. ↩
수렴 차수 그래프에서 기울기가 애매하게 1.3~1.7쯤 나오면 대개 코드가 반쯤 맞는 것이다. 이때 “격자가 아직 점근 영역에 안 들어왔다”고 결론 내리고 넘어가는 것이 인류의 오랜 전통이지만, 절반은 경계조건이 범인이다. ↩