[논문 리뷰] High order schemes based on operator splitting and deferred corrections for stiff time dependent PDEs

Max Duarte, Matthew Emmett|arXiv (Cornell University)|2014. 07. 01.

Numerical methods for differential equations참고 문헌 30인용 수 3

한 줄 요약

이 논문은 연산자 분할과 반복 연기 보정을 조합하여 강한 시간 의존 PDE에 대한 고차 시간 적분 방법을 제시한다. 저자는 체결 노드에서 낮은 차수의 분할 해법을 사용해 초도 근사치를 생성하고, 이를 고차 정확도를 확보하면서도 안정성과 효율성을 유지하는 다항식 기반의 암시적 룬게-쿠타 스킴을 이용해 반복적으로 보정한다.

ABSTRACT

We consider quadrature formulas of high order in time based on Radau-type, L-stable implicit Runge-Kutta schemes to solve time dependent stiff PDEs. Instead of solving a large nonlinear system of equations, we develop a method that performs iterative deferred corrections to compute the solution at the collocation nodes of the quadrature formulas. The numerical stability is guaranteed by a dedicated operator splitting technique that efficiently handles the stiffness of the PDEs and provides initial and intermediate solutions to the iterative scheme. In this way the low order approximations computed by a tailored splitting solver of low algorithmic complexity are iteratively corrected to obtain a high order solution based on a quadrature formula. The mathematical analysis of the numerical errors and local order of the method is carried out in a finite dimensional framework for a general semi-discrete problem, and a time-stepping strategy is conceived to control numerical errors related to the time integration. Numerical evidence confirms the theoretical findings and assesses the performance of the method in the case of a stiff reaction-diffusion equation.

연구 동기 및 목표

• 큰 비선형 시스템을 직접 풀지 않고도 강한 PDE에 대한 고차 시간 적분 스킴을 개발하기 위해.
• L-안정성 있는 라다우 유형의 암시적 룬게-쿠타 스킴의 안정성과 연산자 분할의 효율성을 조합하기 위해.
• 낮은 차수의 분할 근사치에 반복 연기 보정을 적용하여 고차 정확도를 달성하기 위해.
• 강한 준연립 PDE의 시간 적분에서 수치적 안정성과 오차 제어를 확보하기 위해.
• 대규모 PDE 시뮬레이션에 적합한 실용적이고 낮은 복잡도의 고차 시간 적분 프레임워크를 제공하기 위해.

제안 방법

• 고차 정확도를 위한 적분 공식을 정의하기 위해 $s$-단계 $L$-안정성 있는 라다우 유형의 암시적 룬게-쿠타 스킴을 사용한다.
• 각 시간 단계 내에서 $s$개의 체결 노드에서 낮은 차수의 근사를 계산하기 위해 전용 연산자 분할 해법을 활용하여 안정성을 확보하고 강성 제약을 피한다.
• 낮은 차수의 분할 해법을 기반으로 한 해법을 각 체결 노드에서 고차 정확도의 적분 해법으로 향상시키기 위해 반복 연기 보정을 적용한다.
• 반복 보정 과정에서 국소 절단 오차를 관리하기 위해 오차 제어 기반의 시간 단계 전략을 사용한다.
• 스펙트럼 적분 연산자 $S_{t_{i-1}}^{t_i}({oldsymbol{U}})$가 IRK 스킴에서 유도된 체결 노드를 사용하여 적분을 근사하도록 보정 과정을 정의한다.
• 세밀한 전파자(fine propagator)가 분할 기반 수치 적분기로 대체된 SDC(Spectral Deferred Correction)의 변종으로 이 방법을 재해석한다.

실험 결과

연구 질문

• RQ1 큰 비선형 시스템을 직접 풀지 않고도 강한 시간 의존 PDE에 대해 고차 정확도를 달성할 수 있는가?
• RQ2 연산자 분할과 연기 보정을 어떻게 조합하여 안정성을 유지하면서도 고차 수렴을 달성할 수 있는가?
• RQ3 준연립 강한 PDE의 맥락에서 반복 보정 과정의 국소 정확도 차수와 오차 행동은 어떠한가?
• RQ4 이 방법의 시간 단계 전략에서 효과적인 오차 제어를 어떻게 구현할 수 있는가?
• RQ5 정확도와 효율성 측면에서 표준 저차수 분할 또는 완전히 결합된 암시적 룬게-쿠타 스킴에 비해 이 방법이 얼마나 뛰어나게 성능을 발휘하는가?

주요 결과

• 낮은 차수의 분할 해법을 반복 보정하여 고차 정확도를 달성하며, $s$ 단계에 대해 최대 차수 $2s-1$까지 도달한다.
• 강한 반응-확산 방정식에 대한 수치 실험은 이론적 정확도 차수를 확인하고 최적 수렴률을 보여준다.
• 반복 보정 과정은 반복마다 초선형 수렴을 보이며, 이론적 분석과 일치하는 오차 감소를 보인다.
• 오차 제어 기반의 시간 단계 전략은 국소 절단 오차를 성공적으로 관리하고 강한 문제에서 안정성을 유지한다.
• 큰 비선형 시스템을 풀지 않고도 낮은 복잡도의 분할 해법과 반복 보정에 의존함으로써 높은 효율성을 유지한다.
• 이 방법은 세밀한 전파자가 분할 기반 적분기로 대체된 수정된 SDC 스킴과 수학적으로 동치이며, 안정적인 고차 시간 적분을 가능하게 한다.