Related papers: Towards exact symplectic integrators from Liouvill…
The theory of integrable systems of Hamiltonian PDEs and their near-integrable deformations is used to study evolution equations resulting from vertical-averages of the Euler system for two-layer stratified flows in an infinite 2D channel.…
Riemannian manifold Hamiltonian Monte Carlo is traditionally carried out using the generalized leapfrog integrator. However, this integrator is not the only choice and other integrators yielding valid Markov chain transition operators may…
Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
We follow up on our previous works which presented a possible approach for deriving symplectic schemes for a certain class of highly oscillatory Hamiltonian systems. The approach considers the Hamilton-Jacobi form of the equations of…
We develop a reduced-order framework for optimizing mixing in two-dimensional incompressible flows. Instead of optimizing the full transport PDE, the method maximizes the length of advected material interfaces, leading to a…
We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These…
In this paper we present novel integrable symplectic maps, associated with ordinary difference equations, and show how they determine, in a remarkably diverse manner, the integrability, including Lax pairs and the explicit solutions, for…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
Integrals of the Liouville $1$-form, known as the first Poincar\'e integral invariant, provide a computable figure of merit for monitoring the conservation of symplecticity in the numerical integration of Hamiltonian systems. These…
By exploiting the error functions of explicit symplectic integrators for solving separable Hamiltonians, I show that it is possible to develop explicit, time-reversible symplectic integrators for solving non-separable Hamiltonians of the…
In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby…
We present a method for explicit leapfrog integration of inseparable Hamiltonian systems by means of an extended phase space. A suitably defined new Hamiltonian on the extended phase space leads to equations of motion that can be…
The guiding center approximation represents a very powerful tool for analyzing and modeling a charged particle motion in strong magnetic fields. This approximation is based on conservation of the adiabatic invariant, magnetic moment.…
This paper focuses on the numerical approximation of the linearized shallow water equations using hybridizable discontinuous Galerkin (HDG) methods, leveraging the Hamiltonian structure of the evolution system. First, we propose an…
We present a set of new, efficient high-order symplectic methods designed for Hamiltonian systems with cubic or quartic potentials. By demonstrating that polynomial potentials require fewer order conditions, we develop schemes that…
We construct a symplectic integrator for non-separable Hamiltonian systems combining an extended phase space approach of Pihajoki and the symmetric projection method. The resulting method is semiexplicit in the sense that the main time…
In this paper an approach is outlined. With this approach some explicit algorithms can be applied to solve the initial value problem of $n-$dimensional damped oscillators. This approach is based upon following structure: for any…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
Contact Hamiltonian systems extend symplectic Hamiltonian mechanics to dissipative settings while retaining geometric structure. We develop a structure-preserving splitting framework for contact Hamiltonian systems on $J^1(\mathbb{R}^n)$…