Ludwig Gauckler

Ludwig Gauckler

The long-time behaviour of splitting integrators applied to nonlinear Schrödinger equations in a weakly nonlinear setting is studied. It is proven that the energy is nearly conserved on long time intervals. The analysis includes all consistent splitting integrators with real-valued coefficients, in particular splitting integrators of high order. The proof is based on a completely resonant modulated Fourier expansion in time of the numerical solution.

Ludwig Gauckler

An error analysis of trigonometric integrators (or exponential integrators) applied to spatial semi-discretizations of semilinear wave equations with periodic boundary conditions in one space dimension is given. In particular, optimal second-order convergence is shown requiring only that the exact solution is of finite energy. The analysis is uniform in the spatial discretization parameter. It covers the impulse method which coincides with the method of Deuflhard and the mollified impulse method of García-Archilla, Sanz-Serna & Skeel as well as the trigonometric methods proposed by Hairer & Lubich and by Grimm & Hochbruck. The analysis can also be used to explain the convergence behaviour of the Störmer-Verlet/leapfrog discretization in time.

David Cohen, Ludwig Gauckler, Ernst Hairer et al.

For trigonometric and modified trigonometric integrators applied to oscillatory Hamiltonian differential equations with one or several constant high frequencies, near-conservation of the total and oscillatory energies are shown over time scales that cover arbitrary negative powers of the step size. This requires non-resonance conditions between the step size and the frequencies, but in contrast to previous results the results do not require any non-resonance conditions among the frequencies. The proof uses modulated Fourier expansions with appropriately modified frequencies.

Ludwig Gauckler, Ernst Hairer, Christian Lubich

Geometric aspects play an important role in the construction and analysis of structure-preserving numerical methods for a wide variety of ordinary and partial differential equations. Here we review the development and theory of symplectic integrators for Hamiltonian ordinary and partial differential equations, of dynamical low-rank approximation of time-dependent large matrices and tensors, and its use in numerical integrators for Hamiltonian tensor network approximations in quantum dynamics.

Ludwig Gauckler

An error analysis of a splitting method applied to the Zakharov system is given. The numerical method is a Lie-Trotter splitting in time that is combined with a Fourier collocation in space to a fully discrete method. First-order convergence in time and high-order convergence in space depending on the regularity of the exact solution are shown for this method. The main challenge in the analysis is to exclude a loss of spatial regularity in the numerical solution. This is done by transforming the numerical method to new variables and by imposing a natural CFL-type restriction on the discretization parameters.

Ludwig Gauckler, Daniel Weiss

The quadratic nonlinear wave equation on a one-dimensional torus with small initial values located in a single Fourier mode is considered. In this situation, the formation of metastable energy strata has recently been described and their long-time stability has been shown. The topic of the present paper is the correct reproduction of these metastable energy strata by a numerical method. For symplectic trigonometric integrators applied to the equation, it is shown that these energy strata are reproduced even on long time intervals in a qualitatively correct way.

Ludwig Gauckler

Trigonometric integrators for oscillatory linear Hamiltonian differential equations are considered. Under a condition of Hairer & Lubich on the filter functions in the method, a modified energy is derived that is exactly preserved by trigonometric integrators. This implies and extends a known result on all-time near-conservation of energy. The extension can be applied to linear wave equations.

Ludwig Gauckler, Jianfeng Lu, Jeremy L. Marzuola et al.

Trigonometric time integrators are introduced as a class of explicit numerical methods for quasilinear wave equations. Second-order convergence for the semi-discretization in time with these integrators is shown for a sufficiently regular exact solution. The time integrators are also combined with a Fourier spectral method into a fully discrete scheme, for which error bounds are provided without requiring any CFL-type coupling of the discretization parameters. The proofs of the error bounds are based on energy techniques and on the semiclassical Gårding inequality.