[논문 리뷰] Extrapolation-based implicit-explicit general linear methods

Angelamaria Cardone, Z. Jackiewicz외 2명|arXiv (Cornell University)|2013. 04. 08.

Numerical methods for differential equations참고 문헌 13인용 수 6

한 줄 요약

이 논문은 강한 미분방정식의 강성과 비강성 성분으로 자연스럽게 분리되는 경우를 고려하여, 고차수의 암시-명시 일반선형법(IMEX GLMs)을 구성하기 위한 외삽 기반 접근법을 제안한다. 고단계차수 GLMs와 외삽을 활용함으로써 최적의 안정성과 순서 감소를 방지할 수 있으며, 2차원 테스트 문제에서의 수치 결과는 정확성과 효율성을 확인한다.

ABSTRACT

For many systems of differential equations modeling problems in science and engineering, there are natural splittings of the right hand side into two parts, one non-stiff or mildly stiff, and the other one stiff. For such systems implicit-explicit (IMEX) integration combines an explicit scheme for the non-stiff part with an implicit scheme for the stiff part. In a recent series of papers two of the authors (Sandu and Zhang) have developed IMEX GLMs, a family of implicit-explicit schemes based on general linear methods. It has been shown that, due to their high stage order, IMEX GLMs require no additional coupling order conditions, and are not marred by order reduction. This work develops a new extrapolation-based approach to construct practical IMEX GLM pairs of high order. We look for methods with large absolute stability region, assuming that the implicit part of the method is A- or L-stable. We provide examples of IMEX GLMs with optimal stability properties. Their application to a two dimensional test problem confirms the theoretical findings.

연구 동기 및 목표

• 강성 및 비강성 성분을 포함하는 시스템을 위한 고차수 IMEX GLM 쌍 개발
• 기존의 전통적 IMEX 방법에서 흔히 발생하는 순서 감소 문제 해결
• 큰 절대 안정성 영역을 갖는 방법을 구성함으로써 최적의 안정성 특성 확보
• 대규모 과학적 및 공학적 시뮬레이션에 적합한 실용적인 고단계차수 IMEX GLMs 제공
• 2차원 이동-확산-반응 문제에서의 수치 실험을 통해 이론적 결과 검증

제안 방법

• 고단계차수를 갖는 기본 GLMs로부터 고차수 IMEX GLM 쌍을 생성하기 위해 외삽 기반의 구성 기법 사용
• 비강성 성분은 명시적으로, 강성 성분은 암시적으로 처리하는 분할된 GLM 프레임워크 적용
• 외삽 단계의 일致성과 정확성을 보장하기 위해 순서 $ p $ 의 보간 공식 사용
• 결과로 도출된 IMEX GLM에 대해 순서 조건 및 단계차수 조건을 적용하여 $ p $ 차수 정확성과 단계차수 $ q = p $ 또는 $ p-1 $ 보장
• 선형 안정성 분석을 수행하고, 암시 성분에서 $ A $-안정성 또는 $ L $-안정성을 갖는 방법을 구성하여 강성 문제에 대한 강건성 확보
• 두 단계 과정을 통해 방법을 유도함: 먼저 명시적 외삽 계획을 구성한 후, 공유 계수와 절편을 통해 암시적 GLM 성분과 조합

실험 결과

연구 질문

• RQ1 외삽 기법을 효과적으로 활용하여 최적의 안정성 특성을 갖는 고차수 IMEX GLM 쌍을 구성할 수 있는가?
• RQ2 기존의 전통적 IMEX 방법에서 흔히 발생하는 순서 감소 현상과 마찬가지로, 제안된 방법이 강성 문제에서 순서 감소를 피할 수 있는가?
• RQ3 고단계차수가 IMEX GLMs의 순서 조건과 안정성에 어떤 영향을 미치는가?
• RQ4 구성된 IMEX GLMs의 안정성 영역은 기존 방법들과 비교하여 어떻게 다른가?
• RQ5 실제 2차원 이동-확산-반응 문제에서 높은 정확성과 효율성을 유지할 수 있는가?

주요 결과

• 제안된 외삽 기반 IMEX GLM 구성 방법은 최적의 안정성 특성을 확보하며, 큰 절대 안정성 영역을 갖는다.
• 고단계차수 덕분에 추가적인 결합 순서 조건이 필요 없어져 구성 과정이 단순화된다.
• 기존의 다른 IMEX 다단계 방법에서 흔히 발생하는 순서 감소 문제를 피할 수 있으며, 특히 강성 문제가 있는 경우 유리하다.
• 2차원 테스트 문제에서의 수치 실험 결과 이론적 정확도 순서와 방법의 효율성이 확인된다.
• 이론적 분석과 수치 결과는 강성 성분이 있는 경우에도 높은 정확도와 안정성을 유지함을 보여준다.
• 이 방법은 PDE의 공간 반할분할화로부터 유도되는 대규모 시스템에 실용적이고 효과적이다.