[논문 리뷰] Strong stability preserving explicit linear multistep methods with variable step size

Yiannis Hadjimichael, David I. Ketcheson외 2명|arXiv (Cornell University)|2016. 01. 01.

Numerical methods for differential equations참고 문헌 20인용 수 19

한 줄 요약

이 논문은 비선형 초월 PDE의 시간적 적분을 위한 가변 스텝 크기 변수를 가진 강안정성(Strong Stability Preserving, SSP) 선형 다단계 방법(Explicit Linear Multistep Methods, LMMs)을 제시한다. 이는 비선형 초월 PDE에서 파동 속도가 동적으로 변화하는 상황에서 적응형 시간 적분을 가능하게 한다. 논문은 SSP 계수에 대한 날카운 상한선을 확립하고, 최적의 2차 및 3차 방법을 유도하며, 다양한 파동 속도에도 불구하고 단조성과 안정성 특성을 유지하는 탐욕적 스텝 크기 전략 하에서의 안정성 및 수렴성을 증명한다.

ABSTRACT

Strong stability preserving (SSP) methods are designed primarily for time integration of nonlinear hyperbolic PDEs, for which the permissible SSP step size varies from one step to the next. We develop the first SSP linear multistep methods (of order two and three) with variable step size, and prove their optimality, stability, and convergence. The choice of step size for multistep SSP methods is an interesting problem because the allowable step size depends on the SSP coefficient, which in turn depends on the chosen step sizes. The description of the methods includes an optimal step-size strategy. We prove sharp upper bounds on the allowable step size for explicit SSP linear multistep methods and show the existence of methods with arbitrarily high order of accuracy. The effectiveness of the methods is demonstrated through numerical examples.

연구 동기 및 목표

• 비선형 초월 PDE에서 파동 속도가 동적으로 변화하는 상황에서 시간 적분 문제를 해결하기 위해, 가변 스텝 크기 변수를 가진 최초의 강안정성(Strong Stability Preserving, SSP) 선형 다단계 방법(LMMs)을 개발한다.
• 고정된 스텝 크기 SSP LMMs의 한계를 극복한다. 이는 CFL 제한된 스텝 크기가 통합 중에 변화할 경우 비효율적이거나 실패할 수 있기 때문이다.
• 가변 스텝 크기 LMMs의 SSP 계수에 대한 날카운 이론적 경계를 확립한다. 이 경계는 방법 계수와 스텝 크기 비율에 모두 의존한다.
• SSP 성질이 유지되면서 높은 순서 정확도와 안정성을 확보할 수 있는 최적의 스텝 크기 전략을 설계한다.
• 문제와 스텝 크기 수열에 대한 약한 가정 하에서 제안된 최적 방법의 안정성 및 수렴성을 증명한다.

제안 방법

• 스텝 크기 h_n이 단계마다 변하는 것을 고려해, 시간에 따라 변화하는 계수 α_j,n 및 β_j,n를 사용하여 가변 스텝 크기 LMMs를 수립한다.
• 스텝 크기 비율 ω_j = h_{n−k+j} / h_n 과 누적 합 Ω_j = ∑_{i=1}^j ω_i 를 도입하여, 방법을 상대적 스텝 크기로 표현한다.
• 각 단계에서 SSP 성질이 유지되도록, h_n ≤ min_j (α_j,n / β_j,n) × h_{FE}(u_{n−k+j}) 를 만족하는 탐욕적 스텝 크기 전략을 제안한다.
• 계수 비율의 재귀적 분석을 통해 2차 및 3차 방법에 대한 SSP 계수에 대한 날카운 상한선을 유도하고 증명한다.
• 보조 수열 h±_n 과 그 스케일링된 형태 τ±_n 를 도입하여, 스텝 크기 수열이 안정한 극한으로 수렴함을 증명하는 비교 원리를 적용한다.
• 정리 11(유계 비율을 가진 재귀 수열에 관한)을 적용하여, n→∞ 일 때 τ_n → (k−2)/(k−1) 임을 보이며, 이는 스텝 크기 진동의 유계성과 안정성을 보장한다.

실험 결과

연구 질문

• RQ1 2차 및 3차 가변 스텝 크기 선형 다단계 방법에 대해 SSP 계수에 대한 날카운 이론적 상한선은 무엇인가?
• RQ2 변동하는 스텝 크기 조건 하에서 강안정성과 높은 순서 정확도를 유지하면서 최적의 SSP LMMs를 2차 및 3차로 구성할 수 있는가?
• RQ3 파동 속도 변화로 인해 允許 스텝 크기 h_{FE} 가 변할 경우, SSP 성질을 유지할 수 있는 스텝 크기 전략은 어떻게 설계할 수 있는가?
• RQ4 제안된 가변 스텝 크기 SSP LMMs는 선택된 스텝 크기 전략 하에서 안정성 및 수렴성을 유지하는가?
• RQ5 임의로 높은 순서의 SSP LMMs를 가변 스텝 크기로 구성할 수 있는가? 이 가능성은 어떤 제약 조건에 의해 좌우되는가?

주요 결과

• 2차 및 3차 가변 스텝 크기 LMMs에 대한 SSP 계수에 대한 날카운 상한선이 도출되었고, 이는 최적임을 증명하였으며, 이 경계는 스텝 크기 비율 ω_j 에 따라 달라진다.
• 최적의 2차 및 3차 SSP LMMs가 명시적으로 구성되었으며, 계수는 가변 스텝 크기 조건 하에서 SSP 계수를 최대화하도록 선택되었다.
• 논문은 순서-k 방법에 대해 SSP 계수가 (k−2)/(k−1) 이하로 제한됨을 증명하였고, 이 경계는 제안된 전략 하에서 날카롭고 달성 가능하다.
• 탐욕적 스텝 크기 전략이 각 단계에서 h_n ≤ min_j (α_j,n / β_j,n) × h_{FE}(u_{n−k+j}) 를 만족함으로써 SSP 성질이 유지됨을 증명하였으며, 이 전략은 안정성과 수렴성을 유지한다.
• 약한 가정 하에서 스텝 크기 수열 h_n 이 안정한 극한으로 수렴하며, n→∞ 일 때 τ_n = h_n / µ_n → (k−2)/(k−1) 임을 보였다. 이는 장기적 안정성을 보장한다.
• 제안된 스텝 크기 전략 하에서 방법들이 조건부 안정성 없이 안정성과 수렴성을 보임을 증명하였으며, 수렴 차수는 방법의 순서와 일치한다.