Related papers: Comment on `conservative discretizations of the Ke…
We present a remarkable discretization of the classical Kepler problem which preserves its trajectories and all integrals of motion. The points of any discrete orbit belong to an appropriate continuous trajectory.
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
We study four particular 3-dimensional natural Hamiltonian systems defined in conformally Euclidean spaces. We prove their superintegrability and we obtain, in the four cases, the maximal number of functionally independent integrals of…
We discuss the exact discretization of the classical harmonic oscillator equation (including the inhomogeneous case and multidimensional generalizations) with a special stress on the energy integral. We present and suggest some numerical…
We show that any second order linear ordinary diffrential equation with constant coefficients (including the damped and undumped harmonic oscillator equation) admits an exact discretization, i.e., there exists a difference equation whose…
Recently, a family of unconventional integrators for ODEs with polynomial vector fields was proposed, based on the polarization of vector fields. The simplest instance is the by now famous Kahan discretization for quadratic vector fields.…
Traditional numerical discretizations of conservative systems generically yield an artificial secular drift of any nonlinear invariants. In this work we present an explicit nontraditional algorithm that exactly conserves these invariants.…
Discretizations of the Euler top sharing the integrals of motion with the continuous time system are studied. Those of them which are also Poisson with respect to the invariant Poisson bracket of the Euler top are characterized. For all…
We present a method to obtain arbitrarily accurate solutions for conservative classical oscillators. The method that we propose here works both for small and large nonlinearities and provides simple analytical approximations. A comparison…
The accelerated Kepler problem is obtained by adding a constant acceleration to the classical two-body Kepler problem. This setting models the dynamics of a jet-sustaining accretion disk and its content of forming planets as the disk loses…
In this paper we show, in a systematic way, how to relate the Kepler problem to the isotropic harmonic oscillator. Unlike previous approaches, our constructions are carried over in the Lagrangian formalism dealing with second order vector…
A fundamental relation in celestial mechanics is Kepler's equation, linking an orbit's mean anomaly to its eccentric anomaly and eccentricity. Being transcendental, the equation cannot be directly solved for eccentric anomaly by…
We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms…
We develop a numerical scheme for the Kepler problem that preserves exactly all first integrals: angular momentum, total energy, and the Laplace-Runge-Lenz vector. This property ensures that orbital trajectories retain their precise shape…
We consider secular perturbations of nearly Keplerian two-body motion under a perturbing potential that can be approximated to sufficient accuracy by expanding it to second order in the coordinates. After averaging over time to obtain the…
Continuous Galerkin Petrov time discretization scheme is tested on some Hamiltonian systems including simple harmonic oscillator, Kepler's problem with different eccentricities and molecular dynamics problem. In particular, we implement the…
This paper investigates the equations of motion for a relativistic charged particle in a general magnetic field. By reformulating the dynamics in four-dimensional spacetime and separating the linear and nonlinear parts, we construct an…
The problem of Kepler dynamics on a conformable Poisson manifold is addressed. The Hamiltonian function is defined and the related Hamiltonian vector field governing the dynamics is derived, which leads to a modified Newton second law.…
A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the…
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…