[논문 리뷰] Stabilized integrating factor Runge-Kutta method and unconditional preservation of maximum bound principle

Jingwei Li, Xiao Li외 2명|arXiv (Cornell University)|2021. 04. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 44인용 수 5

한 줄 요약

이 논문은 선형 안정화 기법과 적분 인자 룬게쿠타(IfRK) 프레임워크를 조합하여, 선형 포텐셜 방정식과 같은 반선형 포아송 방정식에 대해 최대 bound 원칙(MBP)을 조건 없이 유지하는 안정화된 적분 인자 룬게쿠타(sIFRK) 방법을 제안한다. 저자들은 MBP 유지에 필요한 충분조건를 유도하고, 이를 만족하는 최대 3차 정확도의 sIFRK 스킴을 구성하여 이들이 이산 설정에서 조건 없이 MBP를 유지함을 보장한다. 이는 제1차를 초월하는 많은 기존 고차수 스킴, 특히 SSP-sIFRK 변종과는 달리, 조건 없이 MBP를 유지한다.

ABSTRACT

Maximum bound principle (MBP) is an important property for a large class of semilinear parabolic equations, in the sense that the time-dependent solution of the equation with appropriate initial and boundary conditions and nonlinear operator preserves for all time a uniform pointwise bound in absolute value. It has been a challenging problem on how to design unconditionally MBP-preserving high-order accurate time-stepping schemes for these equations. In this paper, we combine the integrating factor Runge-Kutta (IFRK) method with the linear stabilization technique to develop a stabilized IFRK (sIFRK) method, and successfully derive sufficient conditions for the proposed method to preserve MBP unconditionally in the discrete setting. We then elaborate some sIFRK schemes with up to the third-order accuracy, which are proven to be unconditionally MBP-preserving by verifying these conditions. In addition, it is shown that many classic strong stability-preserving sIFRK schemes do not satisfy these conditions except the first-order one. Extensive numerical experiments are also carried out to demonstrate the performance of the proposed method.

연구 동기 및 목표

• 반선형 포아송 방정식에 대해 최대 bound 원칙(MBP)을 조건 없이 유지하는 고차수 시간 스텝 스킴을 개발한다.
• ETD 및 IFRK와 같은 기존 고차수 지수적 통합 방법이 제2차를 초월해도 MBP를 조건 없이 유지하지 못하는 한계를 극복한다.
• 안정화된 적분 인자 룬게쿠타(sIFRK) 스킴에서 조건 없이 MBP를 유지하기 위한 충분조건를 수립한다.
• 이 조건들을 만족하고 MBP를 조건 없이 유지하는 최대 3차 스킴을 식별하고 검증한다.
• 수치 실험을 통해 많은 고전적 강한 안정성 보존(SSP) sIFRK 스킴이 제1차가 아닌 한, MBP를 조건 없이 유지하지 못한다는 것을 입증한다.

제안 방법

• 적분 인자 룬게쿠타(IFRK) 접근법과 선형 안정화 기법을 조합하여 반선형 포아송 방정식용 sIFRK 스킴을 구성한다.
• Butcher 표의 계수와 안정화 매개변수를 기반으로 조건 없이 MBP 유지에 필요한 충분조건를 유도한다.
• 이 조건들은 최대 3차의 sIFRK 스킴에 대해 해석적으로 검증되어, 시간 스텝 크기와 관계없이 이산 MBP 유지가 보장된다.
• 이 방법은 Allen-Cahn 방정식 및 기타 단계 모델에 적용되며, 유한 차분법 또는 유한 요소법을 통해 공간 이산화가 이루어진다.
• 비선형 항에 대해 안정화를 적용하여 강성 제어를 하고, 조건 없이 안정적이고 MBP를 유지하는 시간 적분을 가능하게 한다.
• 정확도, 안정성, MBP 유지 여부를 평가하기 위해 다양한 시간 스텝 크기와 2D 및 3D 영역에서 수치 검증을 수행한다.

Fig. 1: Evolutions of energies (left) and supremum norms (right) calculated by using sIFRK(1,1) (top), sIFRK( $s$ ,2) (middle) and Heun-sIFRK(3,3) (bottom) schemes with $\tau=0.01$ .

실험 결과

연구 질문

• RQ1 반선형 포아송 방정식에 대해 제2차를 초월하는 고차수(초과 제2차) 시간 통합 스킴이 최대 bound 원칙을 조건 없이 유지할 수 있는가?
• RQ2 이산 설정에서 sIFRK 스킴이 조건 없이 MBP를 유지하기 위해 만족해야 할 조건은 무엇인가?
• RQ3 왜 많은 기존의 강한 안정성 보존(SSP) sIFRK 스킴이 제2차일지라도 MBP를 조건 없이 유지하지 못하는가?
• RQ4 제안된 sIFRK 프레임워크는 유한 차분법 및 유한 요소법과 같은 표준 공간 이산화 방법과 호환되는가?
• RQ5 sIFRK 방법은 큰 시간 스텝 크기에서도 정확도와 안정성을 유지하면서 MBP를 유지할 수 있는가?

주요 결과

• 적절한 안정화 매개변수를 갖춘 sIFRK 방법은 시간 스텝 크기와 관계없이 반선형 포아송 방정식에 대해 최대 bound 원칙(MBP)을 조건 없이 유지한다.
• 3차 sIFRK 스킴은 유도된 충분조건를 만족함으로써 조건 없이 MBP를 유지하는 것으로 성공적으로 구성되고 증명되었다.
• 수치 실험 결과, SSP-sIFRK(2,2)는 MBP 유지 조건을 위반하여 t=5.5 이후에 경계를 유지하지 못하는 것으로 나타났지만, 표준 sIFRK(2,2)는 경계를 완벽하게 유지하였다.
• sIFRK(2,2)의 성능은 다양한 시간 스텝 크기(τ=0.1, 0.05, 0.01)에서 일관되게 유지되었으며, 해의 프로파일과 에너지 변화에 있어 유의미한 차이가 없었다.
• Allen-Cahn 방정식의 3D 시뮬레이션에서 sIFRK(2,2)는 모든 테스트 시간 스텝(포함하여 τ=0.1)에서 MBP를 유지하였고, 에너지 감소가 단조로웠다.
• 본 연구는 대부분의 기존 SSP-sIFRK 스킴이 유도된 조건을 만족하지 못하므로 조건 없이 MBP를 유지할 수 없다는 것을 밝혀냈다. 유일한 예외는 제1차 스킴 뿐이었다.

Fig. 2: Initial configuration of the bubble. Left: the surface-project view; right: the cross-section view at $y=0$ .