Search

There exist excellent codes for an efficient numerical treatment of stiff and differential-algebraic problems. Let us mention {\sc Radau5} which is based on the $3$-stage Radau IIA collocation method, and its extension to problems with discrete delays {\sc Radar5}. The aim of the present work is to present a technique that permits a direct application of these codes to problems having a right-hand side with an additional distributed delay term (which is a special case of an integro-differential equation). Models with distributed delays are of increasing importance in pharmacodynamics and pharmacokinetics for the study of the interaction between drugs and the body. The main idea is to approximate the distribution kernel of the integral term by a sum of exponential functions or by a quasi-polynomial expansion, and then to transform the distributed (integral) delay term into a set of ordinary differential equations. This set is typically stiff and, for some distribution kernels (e.g., Pareto distribution), it contains discrete delay terms with constant delay. The original equations augmented by this set of ordinary differential equations can have a very large dimension, and a careful treatment of the solution of the arising linear systems is necessary. The use of the codes {\sc Radau5} and {\sc Radar5} is illustrated at three examples (two test equations and one problem taken from pharmacodynamics). The driver programs for these examples are publicly available from the homepages of the authors., Comment: 20 pages, 3 figures

A basic leapfrog integrator and its energy-preserving and variational / symplectic variants are proposed and studied for the numerical integration of the equations of motion of relativistic charged particles in an electromagnetic field. The methods are based on a four-dimensional formulation of the equations of motion. Structure-preserving properties of the numerical methods are analysed, in particular conservation and long-time near-conservation of energy and mass shell as well as preservation of volume in phase space. In the non-relativistic limit, the considered methods reduce to the Boris algorithm for non-relativistic charged-particle dynamics and its energy-preserving and variational / symplectic variants., Comment: 18 pages, 3 figures

The Boris algorithm, a closely related variational integrator and a newly proposed filtered variational integrator are studied when they are used to numerically integrate the equations of motion of a charged particle in a non-uniform strong magnetic field, taking step sizes that are much larger than the period of the Larmor rotations. For the Boris algorithm and the standard (unfiltered) variational integrator, satisfactory behaviour is only obtained when the component of the initial velocity orthogonal to the magnetic field is filtered out. The particle motion shows varying behaviour over multiple time scales: fast Larmor rotation, guiding centre motion, slow perpendicular drift, near-conservation of the magnetic moment over very long times and conservation of energy for all times. Using modulated Fourier expansions of the exact and numerical solutions, it is analysed to which extent this behaviour is reproduced by the three numerical integrators used with large step sizes.

A modification of the standard Boris algorithm, called filtered Boris algorithm, is proposed for the numerical integration of the equations of motion of charged particles in a strong non-uniform magnetic field in the asymptotic scaling known as maximal ordering. With an appropriate choice of filters, second-order error bounds in the position and in the parallel velocity, and first-order error bounds in the normal velocity are obtained with respect to the scaling parameter. The proof compares the modulated Fourier expansions of the exact and the numerical solutions. Numerical experiments illustrate the error behaviour of the filtered Boris algorithm.

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., Comment: prepared for the Proceedings of ICM 2018 (Christian Lubich's plenary talk)

A semilinear wave equation with slowly varying wave speed is considered in one to three space dimensions on a bounded interval, a rectangle or a box, respectively. It is shown that the action, which is the harmonic energy divided by the wave speed and multiplied with the diameter of the spatial domain, is an adiabatic invariant: it remains nearly conserved over long times, longer than any fixed power of the time scale of changes in the wave speed in the case of one space dimension, and longer than can be attained by standard perturbation arguments in the two- and three-dimensional cases. The long-time near-conservation of the action yields long-time existence of the solution. The proofs use modulated Fourier expansions in time., Comment: 28 pages

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., Comment: 26 pages

We consider multiscale Hamiltonian systems in which harmonic oscillators with several high frequencies are coupled to a slow system. It is shown that the oscillatory energy is nearly preserved over long times eps^{-N} for arbitrary N>1, where eps^{-1} is the size of the smallest high frequency. The result is uniform in the frequencies and does not require non-resonance conditions., Comment: 13 pages

We compare the performances of symplectic and non-symplectic integrators for the computation of normal geodesics and conjugate points in sub-Riemannian geometry at the example of the Martinet case. For this case study we consider first the flat metric, and then a one parameter perturbation leading to non integrable geodesics. From our computations we deduce that near the abnormal directions a symplectic method is much more efficient for this optimal control problem. The explanation relies on the theory of backward error analysis in geometric numerical integration.