Related papers: Volume Preservation by Runge-Kutta Methods
In this paper, we present a new methodology to develop arbitrary high-order structure-preserving methods for solving the quantum Zakharov system. The key ingredients of our method are: (i) the original Hamiltonian energy is reformulated…
We show that symplectic Runge-Kutta methods provide effective symplectic integrators for Hamiltonian systems with index one constraints. These include the Hamiltonian description of variational problems subject to position and velocity…
We consider the numerical simulation of Hamiltonian systems of ordinary differential equations. Two features of Hamiltonian systems are that energy is conserved along trajectories and phase space volume is preserved by the flow. We want to…
Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…
Infinitesimal volumes stretch and contract as they coevolve with classical phase space trajectories according to linearized dynamics. Unless these tangent-space dynamics are modified, chaotic evolution causes the volume spanned by evolving…
Runge-Kutta methods are affine equivariant: applying a method before or after an affine change of variables yields the same numerical trajectory. However, for some applications, one would like to perform numerical integration after a…
In this paper we study the geometric solution of the so called "good" Boussinesq equation. This goal is achieved by using a convenient space semi-discretization, able to preserve the corresponding Hamiltonian structure, then using…
We generalize the idea of relaxation time stepping methods in order to preserve multiple nonlinear conserved quantities of a dynamical system by projecting along directions defined by multiple time stepping algorithms. Similar to the…
We provide a variational description of any Liouville (i.e. volume preserving) autonomous vector fields on a smooth manifold. This is obtained via a ``maximal degree'' variational principle; critical sections for this are integral manifolds…
Symmetric method and symplectic method are classical notions in the theory of Runge-Kutta methods. They can generate numerical flows that respectively preserve the symmetry and symplecticity of the continuous flows in the phase space.…
Many practical problems can be described by second-order system $\ddot{q}=-M\nabla U(q)$, in which people give special emphasis to some invariants with explicit physical meaning, such as energy, momentum, angular momentum, etc. However,…
The framework of inner product norm preserving relaxation Runge-Kutta methods (David I. Ketcheson, \emph{Relaxation Runge-Kutta Methods: Conservation and Stability for Inner-Product Norms}, SIAM Journal on Numerical Analysis, 2019) is…
This work presents a general geometric framework for simulating and learning the dynamics of Hamiltonian systems that are invariant under a Lie group of transformations. This means that a group of symmetries is known to act on the system…
We propose an extended framework for continuous-stage Runge-Kutta methods which enables us to treat more complicated cases especially for the case weighting on infinite intervals. By doing this, various types of weighted orthogonal…
A numerical method is proposed in order to track field lines of three-dimensional divergence free fields. Field lines are computed by a locally valid Hamiltonian mapping, which is computed using a symplectic scheme. The method is…
In this paper, we give precisely the geometric constraint conditions of canonical symplectic form and regular reduced symplectic forms for the dynamical vector fields of a regular controlled Hamiltonian (RCH) system and its regular reduced…
A new approach to the construction of difference schemes of any order for the many-body problem that preserves all its algebraic integrals is proposed. We introduced additional variables, namely, distances and reciprocal distances between…
Recently, the efficient numerical solution of Hamiltonian problems has been tackled by defining the class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs). Their derivation relies on the expansion of…
Given a vector field on a manifold M, we define a globally conserved quantity to be a differential form whose Lie derivative is exact. Integrals of conserved quantities over suitable submanifolds are constant under time evolution, the…
Structure-preserving geometric algorithm for the Vlasov-Maxwell (VM) equations is currently an active research topic. We show that spatially-discretized Hamiltonian systems for the VM equations admit a local energy conservation law in…