Related papers: A novel energy-conserving scheme for eight-dimensi…
This paper develops the theory of abelian Routh reduction for discrete mechanical systems and applies it to the variational integration of mechanical systems with abelian symmetry. The reduction of variational Runge-Kutta discretizations is…
A novel symplectic integrator for Hamiltonian equations on $S_2^n \times T^{\ast} \RR^m$ is developed and studied. Partitioned Runge--Kutta methods for Hamiltonian systems on products of Hamiltionian manifolds are studied, specifically,…
In this paper, we develop high-order splitting methods for linear port-Hamiltonian systems, focusing on preserving their intrinsic structure, particularly the dissipation inequality. Port-Hamiltonian systems are characterized by their…
The reduction of Hamiltonian systems aims to build smaller reduced models, valid over a certain range of time and parameters, in order to reduce computing time. By maintaining the Hamiltonian structure in the reduced model, certain…
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…
In previous papers, explicit symplectic integrators were designed for nonrotating black holes, such as a Schwarzschild black hole. However, they fail to work in the Kerr spacetime because not all variables can be separable, or not all…
In this paper we propose a novel approach to discretize linear port-Hamiltonian systems while preserving the underlying structure. We present a finite element exterior calculus formulation that is able to mimetically represent conservation…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
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…
This paper aims to develop a linearly implicit structure-preserving numerical scheme for the space fractional sine-Gordon equation, which is based on the newly developed invariant energy quadratization method. First, we reformulate the…
This work concerns the design and analysis of a limiting technique that allows the preservation of invariant domains for high-order numerical approximations of nonlinear hyperbolic systems of conservation laws. The method can be applied to…
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…
In this article, we propose a novel discretization method based on numerical integration for discretizing continuous systems, termed the $\alpha\beta$-approximation or Scalable Bilinear Transformation (SBT). In contrast to existing methods,…
We introduce a class of general purpose linear multisymplectic integrators for Hamiltonian wave equations based on a diamond-shaped mesh. On each diamond, the PDE is discretized by a symplectic Runge--Kutta method. The scheme advances in…
Fully implicit timestepping methods have several potential advantages for atmosphere/ocean simulation. First, being unconditionally stable, they degrade more gracefully as the Courant number increases, typically requiring more solver…
In this manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid…
In this work, a kernel-based surrogate for integrating Hamiltonian dynamics that is symplectic by construction and tailored to large prediction horizons is proposed. The method learns a scalar potential whose gradient enters a…
The rotation-two-component Camassa--Holm system, which possesses strongly nonlinear coupled terms and high-order differential terms, tends to have continuous nonsmooth solitary wave solutions, such as peakons, stumpons, composite waves and…
This work discusses the model reduction problem for large-scale multi-symplectic PDEs with cubic invariants. For this, we present a linearly implicit global energy-preserving method to construct reduced-order models. This allows to…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…