Stabilized integrating factor Runge-Kutta method and unconditional preservation of maximum bound principle

by   Jingwei Li, et al.

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.



There are no comments yet.


page 17

page 18

page 19

page 20

page 21


Maximum bound principle preserving integrating factor Runge-Kutta methods for semilinear parabolic equations

A large class of semilinear parabolic equations satisfy the maximum boun...

Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes

The ubiquity of semilinear parabolic equations has been illustrated in t...

A stability property for a mono-dimensional three velocities scheme with relative velocity

In this contribution, we study a stability notion for a fundamental line...

An upwind DG scheme preserving the maximum principle for the convective Cahn-Hilliard model

The design of numerical approximations of the Cahn-Hilliard model preser...

Order reduction and how to avoid it when Lawson methods integrate reaction-diffusion boundary value problems

It is well known that Lawson methods suffer from a severe order reductio...

Optimized integrating factor technique for Schrödinger-like equations

The integrating factor technique is widely used to solve numerically (in...

Quasi-linear analysis of dispersion relation preservation for nonlinear schemes

In numerical simulations of complex flows with discontinuities, it is ne...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.