# Solving periodic semilinear stiff PDEs in 1D, 2D and 3D with exponential integrators

@article{Montanelli2020SolvingPS, title={Solving periodic semilinear stiff PDEs in 1D, 2D and 3D with exponential integrators}, author={Hadrien Montanelli and Niall Bootland}, journal={Math. Comput. Simul.}, year={2020}, volume={178}, pages={307-327} }

Dozens of exponential integration formulas have been proposed for the high-accuracy solution of stiff PDEs such as the Allen-Cahn, Korteweg-de Vries and Ginzburg-Landau equations. We report the results of extensive comparisons in MATLAB and Chebfun of such formulas in 1D, 2D and 3D, focusing on fourth and higher order methods, and periodic semilinear stiff PDEs with constant coefficients. Our conclusion is that it is hard to do much better than one of the simplest of these formulas, the ETDRK4… Expand

#### Figures, Tables, and Topics from this paper

#### 10 Citations

On the convergence of Lawson methods for semilinear stiff problems

- Computer Science, Mathematics
- Numerische Mathematik
- 2020

It is shown that non-stiff order conditions together with appropriate regularity assumptions imply high-order convergence of Lawson methods, which is in some contrast to the fact that they may have a bad convergence behaviour. Expand

On the Stability of Exponential Integrators for Non-Diffusive Equations

- Mathematics, Computer Science
- ArXiv
- 2021

A simple repartitioning approach is proposed that stabilizes the integrators and enables the efficient solution of stiff, non-diffusive equations and, unlike hyperviscosity, it does not require the use of high-order spatial derivatives. Expand

Exponential Polynomial Time Integrators

- Mathematics, Computer Science
- ArXiv
- 2020

A new class of parallel exponential block methods based on the Legendre points is constructed that can be constructed at arbitrary orders of accuracy, have improved stability compared to existing polynomial based exponential linear multistep methods, and offer significant computational savings compared to current state-of-the-art methods. Expand

Numerical analysis and applications of explicit high order maximum principle preserving integrating factor Runge-Kutta schemes for Allen-Cahn equation

- Mathematics
- 2021

Abstract Whether high order temporal integrators can preserve the maximum principle of Allen-Cahn equation has been an open problem in recent years. This work provides a positive answer by designing… Expand

A spectral method for nonlocal diffusion operators on the sphere

- Computer Science, Physics
- J. Comput. Phys.
- 2018

These algorithms are based on the diagonalizability of nonlocal diffusion operators in the basis of spherical harmonics, the computation of their eigenvalues to high relative accuracy using quadrature and asymptotic formulas, and a fast spherical harmonic transform, which lead to an efficient implementation of high-order exponential integrators for time-dependent models. Expand

Fourier-Spectral Method for the Phase-Field Equations

- Physics
- 2020

In this paper, we review the Fourier-spectral method for some phase-field models: Allen–Cahn (AC), Cahn–Hilliard (CH), Swift–Hohenberg (SH), phase-field crystal (PFC), and molecular beam epitaxy… Expand

Mesh adaptation for pseudospectral ultrasound simulations

- Computer Science
- 2018

A bandwidth-based measure of the spatial resolution requirements for a model solution is developed and integrated into a moving mesh method, and a technique is presented for discretising arbitrary acoustic source distributions that greatly improves the accuracy of source discretisations for uniform meshes as well. Expand

Kernel Learning for Robust Dynamic Mode Decomposition: Linear and Nonlinear Disambiguation Optimization (LANDO)

- Physics, Mathematics
- 2021

Research in modern data-driven dynamical systems is typically focused on the three key challenges of high dimensionality, unknown dynamics, and nonlinearity. The dynamic mode decomposition (DMD) has… Expand

Classification of chaotic time series with deep learning

- Computer Science, Engineering
- ArXiv
- 2019

It is observed that a convolutional neural network without batch normalization layers outperforms state-of-the-art neural networks for time series classification and is able to generalise and classify time series as chaotic or not with high accuracy. Expand

#### References

SHOWING 1-10 OF 83 REFERENCES

Fourth-Order Time-Stepping for Stiff PDEs

- Mathematics, Computer Science
- SIAM J. Sci. Comput.
- 2005

A modification of the exponential time-differencing fourth-order Runge--Kutta method for solving stiff nonlinear PDEs is presented that solves the problem of numerical instability in the scheme as… Expand

Evaluating matrix functions for exponential integrators via Carathéodory-Fejér approximation and contour integrals

- Mathematics
- 2007

Among the fastest methods for solving stiff PDE are exponential integrators, which require the evaluation of , where is a negative semidenite matrix and is the exponential function or one of the… Expand

Exponential Time Differencing for Stiff Systems

- Mathematics
- 2002

We develop a class of numerical methods for stiff systems, based on the method of exponential time differencing. We describe schemes with second- and higher-order accuracy, introduce new Runge?Kutta… Expand

Generalized integrating factor methods for stiff PDEs

- Mathematics
- 2005

The integrating factor (IF) method for numerical integration of stiff nonlinear PDEs has the disadvantage of producing large error coefficients when the linear term has large norm. We propose a… Expand

Exponential Integrators for Large Systems of Differential Equations

- Computer Science, Mathematics
- SIAM J. Sci. Comput.
- 1998

This work studies the numerical integration of large stiff systems of differential equations by methods that use matrix--vector products with the exponential or a related function of the Jacobian, and derives methods up to order 4 which are exact for linear constant-coefficient equations. Expand

A new class of exponential propagation iterative methods of Runge-Kutta type (EPIRK)

- Mathematics, Computer Science
- J. Comput. Phys.
- 2011

Introducing the new, more general than previously proposed, ansatz for EPIRK schemes allows for more flexibility in deriving computationally efficient high-order integrators. Expand

Fourth Order Time-Stepping for Kadomtsev-Petviashvili and Davey-Stewartson Equations

- Mathematics, Computer Science
- SIAM J. Sci. Comput.
- 2011

Purely dispersive partial differential equations as the Korteweg-de Vries equation, the nonlinear Schr\"odinger equation and higher dimensional generalizations thereof can have solutions which… Expand

Exponential Rosenbrock-Type Methods

- Computer Science, Mathematics
- SIAM J. Numer. Anal.
- 2008

A new class of exponential integrators for the numerical integration of large-scale systems of stiff differential equations that linearize the flow in each time step and make use of the matrix exponential and related functions of the Jacobian are introduced. Expand

Exponential multistep methods of Adams-type

- Mathematics
- 2011

The paper is concerned with the construction, implementation and numerical analysis of exponential multistep methods. These methods are related to explicit Adams methods but, in contrast to the… Expand

Exponential integrators

- Computer Science
- Acta Numerica
- 2010

The main intention in this article is to present the mathematics behind these methods, and derive error bounds that are independent of stiffness or highest frequencies in the system. Expand