Krylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations

Implicit integration factor (IIF) methods were developed for solving time-dependent stiff partial differential equations (PDEs) in literature. In [Jiang and Zhang, Journal of Computational Physics, 253 (2013) 368–388], IIF methods are designed to efficiently solve stiff nonlinear advection–diffusion–reaction (ADR) equations. The methods can be designed for an arbitrary order of accuracy. The stiffness of the system is resolved well, and large-time-step-size computations are achieved. To efficiently calculate large matrix exponentials, a Krylov subspace approximation is directly applied to the IIF methods. In this paper, we develop Krylov IIF methods for solving semilinear fourth-order PDEs. As a result of the stiff fourth-order spatial derivative operators, the fourth-order PDEs have much stricter constraints in time-step sizes than the second-order ADR equations. We analyze the truncation errors of the fully discretized schemes. Numerical examples of both scalar equations and systems in one and higher spatial dimensions are shown to demonstrate the accuracy, efficiency and stability of the methods. Large time-step sizes that are of the same order as the spatial grid sizes have been achieved in the simulations of the fourth-order PDEs.

[1]  Andrew J. Christlieb,et al.  Integral deferred correction methods constructed with high order Runge-Kutta integrators , 2009, Math. Comput..

[2]  Weitao Chen,et al.  Semi-implicit integration factor methods on sparse grids for high-dimensional systems , 2015, J. Comput. Phys..

[3]  Yong-Tao Zhang,et al.  Krylov implicit integration factor methods for spatial discretization on high dimensional unstructured meshes: Application to discontinuous Galerkin methods , 2011, J. Comput. Phys..

[4]  J. Hyman,et al.  THE KURAMOTO-SIV ASIDNSKY EQUATION: A BRIDGE BETWEEN POE'S AND DYNAMICAL SYSTEMS , 1986 .

[5]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[6]  Dong Lu,et al.  Computational Complexity Study on Krylov Integration Factor WENO Method for High Spatial Dimension Convection–Diffusion Problems , 2017, Journal of Scientific Computing.

[7]  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.

[8]  Lloyd N. Trefethen,et al.  Fourth-Order Time-Stepping for Stiff PDEs , 2005, SIAM J. Sci. Comput..

[9]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[10]  Yan Xu,et al.  Local discontinuous Galerkin methods for the Cahn-Hilliard type equations , 2007, J. Comput. Phys..

[11]  Jingfang Huang,et al.  Arbitrary order Krylov deferred correction methods for differential algebraic equations , 2007, J. Comput. Phys..

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

[13]  A. Bourlioux,et al.  High-order multi-implicit spectral deferred correction methods for problems of reactive flow , 2003 .

[14]  Qing Nie,et al.  Numerical methods for stiff reaction-diffusion systems , 2007 .

[15]  Yousef Saad,et al.  Efficient Solution of Parabolic Equations by Krylov Approximation Methods , 1992, SIAM J. Sci. Comput..

[16]  Michael L. Minion,et al.  Conservative multi-implicit spectral deferred correction methods for reacting gas dynamics , 2004 .

[17]  Qing Nie,et al.  Compact integration factor methods for complex domains and adaptive mesh refinement , 2010, J. Comput. Phys..

[18]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[19]  Yan Xu,et al.  Local discontinuous Galerkin methods for the Kuramoto-Sivashinsky equations and the Ito-type coupled KdV equations , 2006 .

[20]  Tian Jiang,et al.  Krylov implicit integration factor WENO methods for semilinear and fully nonlinear advection-diffusion-reaction equations , 2013, J. Comput. Phys..

[21]  Su Zhao,et al.  Operator splitting implicit integration factor methods for stiff reaction-diffusion-advection systems , 2011, J. Comput. Phys..

[22]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[23]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[24]  Lili Ju,et al.  Compact implicit integration factor methods for a family of semilinear fourth-order parabolic equations , 2014 .

[25]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[26]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[27]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

[28]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[29]  Einar M. Rønquist,et al.  An Operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow , 1990 .

[30]  Dong Lu,et al.  Krylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection–Diffusion Equations , 2016, Journal of Scientific Computing.

[31]  Xiaolin Zhong,et al.  Additive Semi-Implicit Runge-Kutta Methods for Computing High-Speed Nonequilibrium Reactive Flows , 1996 .

[32]  Jan S. Hesthaven,et al.  Application of implicit-explicit high order Runge-Kutta methods to discontinuous-Galerkin schemes , 2007, J. Comput. Phys..

[33]  J. M. Keiser,et al.  A New Class of Time Discretization Schemes for the Solution of Nonlinear PDEs , 1998 .

[34]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[35]  M. Minion.  Semi-implicit spectral deferred correction methods for ordinary differential equations , 2003 .

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

[37]  Mark Alber,et al.  Macroscopic dynamics of biological cells interacting via chemotaxis and direct contact. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Tian Jiang,et al.  Krylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations , 2016, J. Comput. Phys..

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

[40]  Philipp Birken,et al.  Numerical Linear Algebra , 2011, Encyclopedia of Parallel Computing.

[41]  Willem Hundsdorfer,et al.  RKC time-stepping for advection-diffusion-reaction problems , 2004 .

[42]  Steven J. Ruuth,et al.  Implicit-Explicit Methods for Time-Dependent PDE''s , 1993 .

[43]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[44]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .