Fast Explicit Integration Factor Methods for Semilinear Parabolic Equations

In this paper, an explicit numerical method and its fast implementation are proposed and discussed for the solution of a wide class of semilinear parabolic equations including the Allen–Cahn equation as a special case. The method combines decompositions of compact spatial difference operators on a regular mesh with stable and accurate exponential time integrators and efficient discrete FFT-based algorithms. It can deal with stiff nonlinearity and both homogeneous and inhomogeneous boundary conditions of different types based on multistep approximations and analytic evaluations of time integrals. Numerical experiments demonstrate effectiveness of the new method for both linear and nonlinear model problems.

[1]  Jie Shen,et al.  Applications of semi-implicit Fourier-spectral method to phase field equations , 1998 .

[2]  Tao Tang,et al.  Stability Analysis of Large Time-Stepping Methods for Epitaxial Growth Models , 2006, SIAM J. Numer. Anal..

[3]  Qing Nie,et al.  Efficient semi-implicit schemes for stiff systems , 2006, J. Comput. Phys..

[4]  S. Krogstad.  Generalized integrating factor methods for stiff PDEs , 2005 .

[5]  C. Loan.  Computational Frameworks for the Fast Fourier Transform , 1992 .

[6]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[7]  L. Evans,et al.  Motion of level sets by mean curvature. II , 1992 .

[8]  Jie Shen,et al.  Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method , 2006, J. Comput. Phys..

[9]  Michelle Schatzman,et al.  Geometrical evolution of developed interfaces , 1995, Emerging applications in free boundary problems.

[10]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[11]  Ricardo H. Nochetto,et al.  A posteriori error control for the Allen–Cahn problem: circumventing Gronwall's inequality , 2004 .

[12]  Lei Zhang,et al.  Array-representation integration factor method for high-dimensional systems , 2014, J. Comput. Phys..

[13]  Andreas Prohl,et al.  Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows , 2003, Numerische Mathematik.

[14]  Qiang Du,et al.  Analysis and Applications of the Exponential Time Differencing Schemes and Their Contour Integration Modifications , 2005 .

[15]  Q. Du,et al.  A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , 2004 .

[16]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[17]  Jian Zhang,et al.  Numerical Studies of Discrete Approximations to the Allen--Cahn Equation in the Sharp Interface Limit , 2009, SIAM J. Sci. Comput..

[18]  Yibao Li,et al.  An unconditionally stable hybrid numerical method for solving the Allen-Cahn equation , 2010, Comput. Math. Appl..

[19]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[20]  Qing Nie,et al.  A compact finite difference method for reaction–diffusion problems using compact integration factor methods in high spatial dimensions , 2008, Advances in Difference Equations.