Related papers: Prolongation-Collocation Variational Integrators
We construct numerical integrators for Hamiltonian problems that may advantageously replace the standard Verlet time-stepper within Hybrid Monte Carlo and related simulations. Past attempts have often aimed at boosting the order of accuracy…
In this paper, we present two Hermite polynomial based approaches to derive one-step numerical integrators for mechanical systems. These methods are based on discretizing the configuration using Hermite polynomials which leads to numerical…
In this paper, we propose two novel fourth-order integrators that exhibit uniformly high accuracy and long-term near conservations for solving the nonlinear Dirac equation (NLDE) in the nonrelativistic regime. In this regime, the solution…
In this paper, we develop a structure-preserving discretization of the Lagrangian framework for electromagnetism, combining techniques from variational integrators and discrete differential forms. This leads to a general family of…
We present a unified framework for the construction of localized exponential integrators that bypasses the traditional trade-off between the accuracy of global spectral methods and the efficiency of sparse finite differences. By evaluating…
Multi-derivative one-step methods based upon Euler-Maclaurin integration formulae are considered for the solution of canonical Hamiltonian dynamical systems. Despite the negative result that simplecticity may not be attained by any…
We consider (3+1)-dimensional second-order evolutionary PDEs where the unknown $u$ enters only in the form of the 2nd-order partial derivatives. For such equations which possess a Lagrangian, we show that all of them have a symplectic…
The Floater--Hormann family of rational interpolants do not have spurious poles or unattainable points, are efficient to calculate, and have arbitrarily high approximation orders. One concern when using them is that the amplification of…
This article revisits the approximation problem of systems of nonlinear delay differential equations (DDEs) by a set of ordinary differential equations (ODEs). We work in Hilbert spaces endowed with a natural inner product including a point…
In the last two decades, significant effort has been put in understanding and designing so-called structure-preserving numerical methods for the simulation of mechanical systems. Geometric integrators attempt to preserve the geometry…
Interconnected systems are an important class of mathematical models, as they allow for the construction of complex, hierarchical, multiphysics, and multiscale models by the interconnection of simpler subsystems. Lagrange--Dirac mechanical…
The discrete gradient methods are integrators designed to preserve invariants of ordinary differential equations. From a formal series expansion of a subclass of these methods, we derive conditions for arbitrarily high order. We derive…
In this paper we show that a variational reduction procedure can be defined for Lagrangian systems subject to scaling symmetries (i.e. Lagrangian systems defined by a homogenous Lagrangian function), in such a way that the trajectories of…
We establish error estimates for semi-Lagrangian schemes for the initial value problem of one-dimensional conservation laws with a dispersive term, including the Korteweg--de Vries equation. The schemes considered in this paper are based on…
We show how to increase the order of one-dimensional discrete gradient numerical integrator without losing its advantages, such as exceptional stability, exact conservation of the energy integral and exact preservation of the trajectories…
Recently, a family of unconventional integrators for ODEs with polynomial vector fields was proposed, based on the polarization of vector fields. The simplest instance is the by now famous Kahan discretization for quadratic vector fields.…
In this paper, discrete analogues of Euler-Poincar\'{e} and Lie-Poisson reduction theory are developed for systems on finite dimensional Lie groups $G$ with Lagrangians $L:TG \to {\mathbb R}$ that are $G$-invariant. These discrete equations…
Using Dirac's approach to constrained dynamics, the Hamiltonian formulation of regular higher order Lagrangians is developed. The conventional description of such systems due to Ostrogradsky is recovered. However, unlike the latter, the…
Conservative symmetric second-order one-step integrators are derived using the Discrete Multiplier Method for a family of vortex-blob models approximating the incompressible Euler's equations on the plane. Conservative properties and second…
We introduce the differential, integral, and variational delta-embeddings. We prove that the integral delta-embedding of the Euler-Lagrange equations and the variational delta-embedding coincide on an arbitrary time scale. In particular, a…