강성 행렬 미분 방정식의 저차원 스트랭 분할법 수렴성 분석

초록

본 논문은 Sylvester 형태를 갖는 강성 행렬 미분 방정식에 대해, 선형 부분을 정확히 행렬 지수로 처리하고 비선형 부분을 2차 동적 저차원(BUG2) 적분기로 풀어내는 스트랭 분할 스키마를 제안한다. 가정 하에 전체 스키마가 2차 정확도를 보이며, 수치 실험을 통해 이론적 수렴 결과와 효율성을 검증한다.

상세 분석

이 연구는 시간 독립 행렬 A와 비강성 비선형 연산자 G(t,X) 로 구성된 Sylvester 형태의 방정식 Ẋ = AX + XA* + G(t,X)를 대상으로 한다. 저차원 근사에 적합한 구조적 특성을 활용해, 해를 고정된 순위 r 의 매니폴드 M_r 위에 제한한다. 선형 부분은 정확히 Φ_A^{τ}(·)=e^{τA}(·)e^{τA*} 로 풀어내어 저차원 형태가 보존됨을 보이고, 비선형 부분은 동적 저차원 프레임워크를 이용해 접공간 T_Y M_r 로 투사한다. 투사된 방정식 Ẏ = P(Y)G(t,Y)를 2차 중점 BUG2 스키마로 적분함으로써, 시간 정확도 2차와 저차원 유지라는 두 목표를 동시에 달성한다.

수렴 분석은 크게 세 부분으로 나뉜다. 첫째, 전체(풀랭크) 스트랭 분할 S_τ = Φ_A^{τ/2}∘Φ_G^{τ}∘Φ_A^{τ/2} 가 강성에 무관하게 전역 2차 정확도 ‖X(t₀+nτ)−S_τⁿ(X₀)‖ ≤ C τ²(1+|log τ|) 를 만족함을 기존 이론(예: