[논문 리뷰] Method of lines transpose: High order L-stable O(N) schemes for parabolic equations using successive convolution

Matthew F. Causley, Hana Cho외 2명|arXiv (Cornell University)|2015. 08. 13.

Numerical methods for differential equations참고 문헌 42인용 수 5

한 줄 요약

이 논문은 고차수, L-안정성 있는 $Ø(N)$ 해법을 제안하며, 선형 및 비선형 확산 문제(Allen-Cahn 및 Fitzhugh-Nagumo 방정식 포함)를 효과적이고 안정적이며 정확하게 해결할 수 있도록, 선형 미분 방정식의 방법(transpose)과 연속적 복합화(해결 함수 전개)를 결합한다. 이 방법은 공간 및 시간에 대해 임의의 정밀도를 달성할 수 있으며, 강성 감쇠가 입증되고 $Ø(N)$ 계산 비용을 보장한다.

ABSTRACT

We present a new solver for nonlinear parabolic problems that is L-stable and achieves high order accuracy in space and time. The solver is built by first constructing a single-dimensional heat equation solver that uses fast O(N) convolution. This fundamental solver has arbitrary order of accuracy in space, and is based on the use of the Green's function to invert a modified Helmholtz equation. Higher orders of accuracy in time are then constructed through a novel technique known as successive convolution (or resolvent expansions). These resolvent expansions facilitate our proofs of stability and convergence, and permit us to construct schemes that have provable stiff decay. The multi-dimensional solver is built by repeated application of dimensionally split independent fundamental solvers. Finally, we solve nonlinear parabolic problems by using the integrating factor method, where we apply the basic scheme to invert linear terms (that look like a heat equation), and make use of Hermite-Birkhoff interpolants to integrate the remaining nonlinear terms. Our solver is applied to several linear and nonlinear equations including heat, Allen-Cahn, and the Fitzhugh-Nagumo system of equations in one and two dimensions.

연구 동기 및 목표

• 공간 및 시간에 대해 임의의 정밀도를 갖는 고차수, L-안정성 있고 계산적으로 효율적인 포물형 PDE 해법을 개발한다.
• 伝통적인 ADI 및 은밀한 방법의 한계를 해결하기 위해 메모리 사용량과 계산 비용을 줄이면서도 안정성과 정밀도를 유지한다.
• 해결 함수 전개를 사용하여, 선형 미분 방정식의 방법(transpose)을 확장하여 강성 감쇠와 증명 가능한 수렴성을 갖는 포물형 문제에 적용한다.
• 고차수 허미트-버크호프 보간법을 통합 인테그레이팅 팩터 방법과 융합하여 비선형 반응-확산 시스템을 효율적으로 해결한다.
• 1D 및 2D에서 Allen-Cahn 및 Fitzhugh-Nagumo 방정식을 포함한 벤치마크 문제에 대해 이 방법의 강건성과 효율성을 입증한다.

제안 방법

• 이 방법은 선형 미분 방정식의 방법(transpose)을 사용하며, 시간에 대해 먼저 이산화하고, 이를 바탕으로 수정된 헬름홀츠 방정식을 그린 함수와 함께 $Ø(N)$ 빠른 복합화를 통해 해결한다.
• 고차수 공간 정밀도는 경계 적분의 적분 기반 이산화를 통해 달성되며, 스무스 문제에 대해 공간에서 스펙트럼 정밀도를 확보한다.
• 연속적 복합화(해결 함수 전개)를 사용하여 고차수 시간 적분기를 구성하며, 시간에 대해 임의의 정밀도를 확보할 수 있다.
• 해결 함수 전개 프레임워크를 통해 안정성과 수렴성을 증명하며, 모든 정밀도 순서에서 강성 감쇠를 보장한다.
• 다차원 문제의 경우 차원 분할 해법을 순차적으로 적용하여 $Ø(N)$ 복잡도를 유지한다.
• 비선형 항은 통합 인테그레이팅 팩터 방법을 통해 처리되며, 선형 확산 항은 은밀하게 해결되고 비선형 항은 고차수 허미트-버크호프 보간법을 사용해 통합된다.

실험 결과

연구 질문

• RQ1 MOL T와 연속적 복합화를 사용하여 선형 열 방정식에 대해 고차수, L-안정성 있고 $Ø(N)$ 스킴을 구성할 수 있는가?
• RQ2 연속적 복합화 기법이 시간에 대해 임의의 정밀도 순서에서 강성 감쇠와 수렴성을 보장하는가?
• RQ3 이 방법은 Allen-Cahn 및 Fitzhugh-Nagumo와 같은 비선형 포물형 시스템에 대해 고차수 정밀도와 안정성을 갖도록 확장될 수 있는가?
• RQ4 1D 및 2D에서 선형 및 비선형 문제에 대해 수렴 속도와 정밀도 측면에서 이 방법의 성능은 어떠한가?
• RQ5 강성 문제에서 흔히 발생하는 엄격한 시간 간격 제약 조건 하에서도 이 방법은 효율성과 안정성을 유지할 수 있는가?

주요 결과

• 제안된 스킴은 시간 단계당 $\mathcal{O}(N)$ 계산 비용을 보이며, 대규모 문제에 대해 매우 효율적이다.
• 1D 열 방정식의 경우, 시간에 대해 $\mathcal{O}(\Delta t^p)$ 수렴을 보이며, 각각 두 번째, 네 번째, 여섯 번째 차수 스킴에 대해 $p=2,4,6$이다.
• 2D에서는 열 방정식에 대해 시간에 대해 2차 수렴과 공간에 대해 스펙트럼 수렴을 유지한다.
• Allen-Cahn 방정식의 경우, 시간에 대해 2차 수렴을 보이며, $L^\infty$ 오차가 $\Delta t$를 반으로 줄일수록 $5.17 \times 10^{-6}$에서 $1.296 \times 10^{-8}$로 감소한다.
• Fitzhugh-Nagumo 시스템은 $\Delta t = 0.0008$인 두 번째 차수 스킴을 사용하여 나선형 웨이브 형성이 성공적으로 시뮬레이션되었으며, 기존 문헌 결과와 일치한다.
• 모든 테스트 정밀도 순서에서 강성 감쇠가 관찰되어, L-안정성과 강성 있는 포물형 문제에 대한 강건성을 확인한다.