Related papers: Symplectic Runge-Kutta Methods for Hamiltonian Sys…
We propose an implementation of symplectic implicit Runge-Kutta schemes for highly accurate numerical integration of non-stiff Hamiltonian systems based on fixed point iteration. Provided that the computations are done in a given floating…
We prove that a class of A-stable symplectic Runge--Kutta time semidiscretizations (including the Gauss--Legendre methods) applied to a class of semilinear Hamiltonian PDEs which are well-posed on spaces of analytic functions with analytic…
We propose a family of integrators, Flow-Composed Implicit Runge-Kutta (FCIRK) methods, for perturbations of nonlinear ordinary differential equations, consisting of the composition of flows of the unperturbed part alternated with one step…
Numerical integrators could be used to form interpolation conditions when training neural networks to approximate the vector field of an ordinary differential equation (ODE) from data. When numerical one-step schemes such as the Runge-Kutta…
We introduce a new class of parametricization structure-preserving partitioned Runge-Kutta ($\alpha$-PRK) methods for Hamiltonian systems with holonomic constraints. When the scalar parameter $\alpha=0$, the methods are reduced to the usual…
In recent years, the efficient numerical solution of Hamiltonian problems has led to the definition of a class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs). Such methods admit an interesting…
This work develops a symplectic framework for quantum computing to be applied to classical Hamiltonian systems, exploiting the intrinsic geometric compatibility between unitary quantum evolution and symplectic phase-space dynamics in a…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
We study the construction and convergence of semi-explicit and iterative decoupling schemes for an elliptic-parabolic problem using higher-order Runge-Kutta methods. For the semi-explicit schemes, which are constructed using a nearby delay…
We reconsider the variational derivation of symplectic partitioned Runge-Kutta schemes. Such type of variational integrators are of great importance since they integrate mechanical systems with high order accuracy while preserving the…
We are concerned with the efficient implementation of symplectic implicit Runge-Kutta (IRK) methods applied to systems of (non-necessarily Hamiltonian) ordinary differential equations by means of Newton-like iterations. We pay particular…
New time integration methods are proposed for simulating incompressible multiphase flow in pipelines described by the one-dimensional two-fluid model. The methodology is based on 'half-explicit' Runge-Kutta methods, being explicit for the…
Most physical processes posses structural properties such as constant energies, volumes, and other invariants over time. When learning models of such dynamical systems, it is critical to respect these invariants to ensure accurate…
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic…
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…
We consider a Lagrangian system $L(q,\dot q) = \sum_{l=1}^{N}L^{\{l\}}(q,\dot q)$, where the $q$-variable is treated by a Generalized Additive Runge--Kutta (GARK) method. Applying the technique of discrete variations, we show how to…
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.…
We develop Chebyshev symplectic methods based on Chebyshev orthogonal polynomials of the first and second kind separately in this paper. Such type of symplectic methods can be conveniently constructed with the newly-built theory of weighted…
An error analysis of Runge-Kutta convolution quadrature based on Gauss methods applied to hyperbolic operators is given. The order of convergence relies heavily on the parity of the number of stages, a more favourable situation arising for…
Hamiltonian systems are differential equations which describe systems in classical mechanics, plasma physics, and sampling problems. They exhibit many structural properties, such as a lack of attractors and the presence of conservation…