Efficient low-storage Runge-Kutta schemes with optimized stability regions

A variety of numerical calculations, especially when considering wave propagation, are based on the method-of-lines, where time-dependent partial differential equations (PDEs) are first discretized in space. For the remaining time-integration, low-storage Runge-Kutta schemes are particularly popular due to their efficiency and their reduced memory requirements. In this work, we present a numerical approach to generate new low-storage Runge-Kutta (LSRK) schemes with optimized stability regions for advection-dominated problems. Adapted to the spectral shape of a given physical problem, those methods are found to yield significant performance improvements over previously known LSRK schemes. As a concrete example, we present time-domain calculations of Maxwell's equations in fully three-dimensional systems, discretized by a discontinuous Galerkin approach.

[1]  Kurt Busch,et al.  Spatio-spectral characterization of photonic meta-atoms with electron energy-loss spectroscopy , 2011 .

[2]  J. Williamson.  Low-storage Runge-Kutta schemes , 1980 .

[3]  P. Houwen.  Explicit Runge-Kutta formulas with increased stability boundaries , 1972 .

[4]  M. Nallasamy,et al.  High-Accuracy Large-Step Explicit Runge-Kutta (HALE-RK) Schemes for Computational Aeroacoustics , 2006 .

[5]  W. Habashi,et al.  2N-Storage Low Dissipation and Dispersion Runge-Kutta Schemes for Computational Acoustics , 1998 .

[6]  M. Ali,et al.  Some Variants of the Controlled Random Search Algorithm for Global Optimization , 2006 .

[7]  R. Lewis,et al.  Low-storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations , 2000 .

[8]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

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

[10]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[11]  K. Busch,et al.  Computing electron energy loss spectra with the Discontinuous Galerkin Time-Domain method , 2011 .

[12]  C. Balanis.  Advanced Engineering Electromagnetics , 1989 .

[13]  M. E. Galassi,et al.  GNU SCIENTI C LIBRARY REFERENCE MANUAL , 2005 .

[14]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[15]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[16]  K. Busch,et al.  Analysis of light propagation in slotted resonator based systems via coupled-mode theory. , 2011, Optics express.

[17]  P. Houwen,et al.  On the Internal Stability of Explicit, m‐Stage Runge‐Kutta Methods for Large m‐Values , 1979 .

[18]  Kurt Busch,et al.  Higher-order time-domain methods for the analysis of nano-photonic systems , 2009 .

[19]  Manuel Torrilhon,et al.  Essentially optimal explicit Runge–Kutta methods with application to hyperbolic–parabolic equations , 2007, Numerische Mathematik.

[20]  Robert Vichnevetsky,et al.  New stability theorems concerning one-step numerical methods for ordinary differential equations , 1983 .

[21]  William G. Gray,et al.  One step integration methods of third-fourth order accuracy with large hyperbolic stability limits , 1984 .

[22]  Timothy C. Warburton,et al.  Nodal discontinuous Galerkin methods on graphics processors , 2009, J. Comput. Phys..

[23]  I. P. E. Kinnmark.  A principle for construction of one-step integration methods with maximum imaginary stability limits , 1987 .

[24]  Kurt Busch,et al.  Comparison of Low-Storage Runge-Kutta Schemes for Discontinuous Galerkin Time-Domain Simulations of Maxwell's Equations , 2010 .

[25]  David I. Ketcheson,et al.  Highly Efficient Strong Stability-Preserving Runge-Kutta Methods with Low-Storage Implementations , 2008, SIAM J. Sci. Comput..