Related papers: On conservative difference schemes for the many-bo…
It is well-known that a numerical method which is at the same time geometric structure-preserving and physical property-preserving cannot exist in general for Hamiltonian partial differential equations. In this paper, we present a novel…
We employ a port-Hamiltonian approach to model nonlinear rigid multibody systems subject to both position and velocity constraints. Our formulation accommodates Cartesian and redundant coordinates, respectively, and captures kinematic as…
It is well known that symplectic Runge-Kutta and Partitioned Runge-Kutta methods exactly preserve {\em quadratic} first integrals (invariants of motion) of the system being integrated. While this property is often seen as a mere curiosity…
Conservative symmetric second-order one-step schemes are derived for dynamical systems describing various many-body systems using the Discrete Multiplier Method. This includes conservative schemes for the $n$-species Lotka-Volterra system,…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
This paper introduces a new difference scheme to the difference equations for N-body type problems. To find the non-collision periodic solutions and generalized periodic solutions in multi-radial symmetric constraint for the N-body type…
We show a novel systematic way to construct conservative finite difference schemes for quasilinear first-order system of ordinary differential equations with conserved quantities. In particular, this includes both autonomous and…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
In this work, we present a modification of explicit Runge-Kutta temporal integration schemes that guarantees the preservation of any locally-defined quasiconvex set of bounds for the solution. These schemes operate on the basis of a…
Stochastic Hamiltonian partial differential equations, which possess the multi-symplectic conservation law, are an important and fairly large class of systems. The multi-symplectic methods inheriting the geometric features of stochastic…
In a recent series of papers, the class of energy-conserving Runge-Kutta methods named Hamiltonian BVMs (HBVMs) has been defined and studied. Such methods have been further generalized for the efficient solution of general conservative…
Recently, the class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs), has been proposed for the efficient solution of Hamiltonian problems, as well as for other types of conservative problems. In…
In this paper, we present a novel class of high-order energy-preserving schemes for solving the Zakharov-Rubenchik equations. The main idea of the scheme is first to introduce an quadratic auxiliary variable to transform the Hamiltonian…
This work is devoted to a systematic exposition of the dynamics of a rigid body, considered as a system with kinematic constraints. Having accepted the variational problem in accordance with this, we no longer need any additional postulates…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
There are periodic solutions to the equal-mass three-body (and N-body) problem in Newtonian gravity. The figure-eight solution is one of them. In this paper, we discuss its solution in the first and second post-Newtonian approximations to…
In this paper, we propose linearly implicit and arbitrary high-order conservative numerical schemes for ordinary differential equations with a quadratic invariant. Many differential equations have invariants, and numerical schemes for…
We use Lagrangian formalism and jet spaces to derive a computational model to simulate multibody dynamics with holonomic constraints. Our approach avoids the traditional problems of drift-off and spurious oscillations. Hence even long…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. This study explores…