Poisson integrators

An overview of Hamiltonian systems with noncanonical Poisson structures is given. Examples of bi-Hamiltonian ode's, pde's and lattice equations are presented. Numerical integrators using generating functions, Hamiltonian splitting, symplectic Runge-Kutta methods are discussed for Lie-Poisson systems and Hamiltonian systems with a general Poisson structure. Nambu-Poisson systems and the discrete gradient methods are also presented.Comment: 30 page

Stability in the Numerical Treatment of Volterra Integral and Integro-Differential Equations with emphasis on Finite Recurrence Relations.

In the last two decades the theory of Volterra integral equations and of integro-differential equations has developed extensively. New classes of methods for the numerical solution of such equations have been developed and at the same time there have been advances in the qualitative theory of these equations. More frequent use is being made of Volterra equations to model various physical and biological phenomenon as new insight has occurred into the asymptotic behaviour of solutions. In consequence, there has emerged a need for reliable and efficient methods for the numerical treatment of such equations. This thesis is concerned with an aspect of numerical solution of Volterra integral and integro-differential equations. In Chapters 1 and 2 we are concerned with background material. We provide results on the classical theory of Volterra equations in Chapter 1 and on numerical methods in Chapter 2. The original material is contained in Chapters 3, 4 and 5. Here, stability results which involve the construction and analysis of finite-term recurrence relations are presented. The techniques relate to the treatment of Volterra integral and integro-differential equations. They permit the analysis of classical and 7-modified numerical methods. The results presented should be viewed as a contribution towards an understanding of numerical stability for the methods considered. The area is one in which further work (subsequent to the present investigation and involving advanced techniques) has been performed and where open questions still remain. The techniques which are employed in this thesis are applicable in other areas of numerical analysis and therefore have intrinsic interest

Variational Partitioned Runge–Kutta Methods for Lagrangians Linear in Velocities

In this paper, we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the “Hamiltonian” equations of motion can be formulated as an index-1 differential-algebraic system. We also construct variational Runge–Kutta methods and analyze their properties. The general properties of Runge–Kutta methods depend on the “velocity” part of the Lagrangian. If the “velocity” part is also linear in the position coordinate, then we show that non-partitioned variational Runge–Kutta methods are equivalent to integration of the corresponding first-order Euler–Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge–Kutta method are retained. If the “velocity” part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We verified our results through numerical experiments for various dynamical systems

Parallel algorithm with spectral convergence for nonlinear integro-differential equations

We discuss a numerical algorithm for solving nonlinear integro-differential equations, and illustrate our findings for the particular case of Volterra type equations. The algorithm combines a perturbation approach meant to render a linearized version of the problem and a spectral method where unknown functions are expanded in terms of Chebyshev polynomials (El-gendi's method). This approach is shown to be suitable for the calculation of two-point Green functions required in next to leading order studies of time-dependent quantum field theory.Comment: 15 pages, 9 figure