Linearly implicit methods for nonlinear PDEs with linear dispersion and dissipation

The linear stability of IMEX (IMplicit-EXplicit) methods and exponential integrators for stiff systems of ODEs arising in the discrete solution of PDEs is examined for nonlinear PDEs with both linear dispersion and dissipation, and a clear method of visualization of the linear stability regions is proposed. Predictions are made based on these visualizations and are supported by a series of experiments on five PDEs including quasigeostrophic equations and stratified Boussinesq equations. The experiments, involving 24 IMEX and exponential methods of third and fourth order, confirm the predictions of the linear stability analysis, that the methods are typically limited by small eigenvalues of the linear term and by eigenvalues on or near the imaginary axis rather than by large eigenvalues near the negative real axis. The experiments also demonstrate that IMEX methods achieve comparable stability to exponential methods, and that exponential methods are significantly more accurate only when the problem is nearly linear. Novel IMEX predictor-corrector methods are also derived.

[1]  H. A. Ashi,et al.  Comparison of methods for evaluating functions of a matrix exponential , 2009 .

[2]  I. Beresnev,et al.  A model for nonlinear seismic waves in a medium with instability , 1993 .

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

[4]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[5]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[6]  Marlis Hochbruck,et al.  Explicit Exponential Runge-Kutta Methods for Semilinear Parabolic Problems , 2005, SIAM J. Numer. Anal..

[7]  Willem Hundsdorfer,et al.  Stability of implicit-explicit linear multistep methods , 1997 .

[8]  T. Driscoll,et al.  Regular Article: A Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion , 1999 .

[9]  A. Staniforth,et al.  Semi-Lagrangian integration schemes for atmospheric models - A review , 1991 .

[10]  Willem Hundsdorfer,et al.  IMEX extensions of linear multistep methods with general monotonicity and boundedness properties , 2007, J. Comput. Phys..

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

[12]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

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

[14]  G. Russo,et al.  Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation , 2005 .

[15]  Jean Côté,et al.  A Two-Time-Level Semi-Lagrangian Semi-implicit Scheme for Spectral Models , 1988 .

[16]  Elena Celledoni,et al.  Semi-Lagrangian Runge-Kutta Exponential Integrators for Convection Dominated Problems , 2009, J. Sci. Comput..

[17]  Willem Hundsdorfer,et al.  High-order linear multistep methods with general monotonicity and boundedness properties , 2005 .

[18]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[19]  D. A. Voss,et al.  A linearly implicit predictor–corrector scheme for pricing American options using a penalty method approach , 2006 .

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

[21]  A. Ostermann,et al.  A Class of Explicit Exponential General Linear Methods , 2006 .

[22]  Takuji Kawahara,et al.  Approximate Equations for Long Nonlinear Waves on a Viscous Fluid , 1978 .

[23]  S. Koikari.  Rooted tree analysis of Runge-Kutta methods with exact treatment of linear terms , 2005 .

[24]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[25]  M. Calvo,et al.  Linearly implicit Runge—Kutta methods for advection—reaction—diffusion equations , 2001 .

[26]  E. T. WHITTAKER,et al.  Partial Differential Equations of Mathematical Physics , 1932, Nature.

[27]  F. Dias,et al.  One-dimensional wave turbulence , 2004 .

[28]  Andrew J. Majda,et al.  A one-dimensional model for dispersive wave turbulence , 1997 .

[29]  S. Cox,et al.  Nikolaevskiy equation with dispersion. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  H. Bateman.  Partial Differential Equations of Mathematical Physics , 1932 .

[31]  Sebastiano Boscarino,et al.  On an accurate third order implicit-explicit Runge--Kutta method for stiff problems , 2009 .

[32]  B. Minchev,et al.  A review of exponential integrators for first order semi-linear problems , 2005 .

[33]  Brynjulf Owren,et al.  B-series and Order Conditions for Exponential Integrators , 2005, SIAM J. Numer. Anal..

[34]  D. A. Voss,et al.  A linearly implicit predictor-corrector method for reaction-diffusion equations , 1999 .

[35]  Toshiyuki Koto,et al.  Stability of implicit-explicit linear multistep methods for ordinary and delay differential equations , 2009 .

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

[37]  Yu Ito,et al.  Stability of Localized Pulse Trains with a Long Tail in the Generalized Kuramoto-Sivashinsky Equation(Electromagnetism, optics, acoustics, heat transfer, classical mechanics, and fluid mechanics) , 2008 .

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

[39]  J. M. Sanz-Serna,et al.  Symplectic Methods Based on Decompositions , 1997 .

[40]  T. Driscoll.  A composite Runge-Kutta method for the spectral solution of semilinear PDEs , 2002 .

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

[42]  Thor Gjesdal.  Implicit--explicit methods based on strong stability preserving multistep time discretizations , 2003 .

[43]  Willem Hundsdorfer,et al.  Partially Implicit BDF2 Blends for Convection Dominated Flows , 2000, SIAM J. Numer. Anal..

[44]  Self-similarity in decaying two-dimensional stably stratified adjustment , 2006, physics/0607045.

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

[46]  Håvard Berland,et al.  NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET , 2005 .

[47]  Mari Paz Calvo,et al.  A class of explicit multistep exponential integrators for semilinear problems , 2006, Numerische Mathematik.