Related papers: Numerical solution of perturbed Kepler problem usi…
Nonlinear elliptic problems arise in many fields, including plasma physics, astrophysics, and optimal transport. In this article, we propose a novel operator-splitting/finite element method for solving such problems. We begin by introducing…
In a previous work, we developed the idea to solve Kepler's equation with a CORDIC-like algorithm, which does not require any division, but still multiplications in each iteration. Here we overcome this major shortcoming and solve Kepler's…
The existence of elliptic periodic solutions of a perturbed Kepler problem is proved. The equations are in the plane and the perturbation depends periodically on time. The proof is based on a local description of the symplectic group in two…
We use the Melnikov integral method to prove that the Hamiltonian flow on the zero-energy manifold for the Kepler problem perturbed by a quadrupole moment is chaotic, irrespective of the perturbation being of prolate or oblate type. This…
We consider the numerical integration of the Gross-Pitaevskii equation with a potential trap given by a time-dependent harmonic potential or a small perturbation thereof. Splitting methods are frequently used with Fourier techniques since…
This overview is devoted to splitting methods, a class of numerical integrators intended for differential equations that can be subdivided into different problems easier to solve than the original system. Closely connected with this class…
It is common in classical mechanics to encounter systems whose Hamiltonian $H$ is the sum of an often exactly integrable Hamiltonian $H_0$ and a small perturbation $\epsilon H_1$ with $\epsilon\ll1$. Such near-integrability can be exploited…
We present a new symplectic integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine…
We investigate the Kepler problem using a symplectic structure consistent with the commutation rules of the noncommutative quantum mechanics. We show that a noncommutative parameter of the order of $10^{-58} \text m^2$ gives observable…
In view of the existing limitations of sequential computing, parallelization has emerged as an alternative in order to improve the speedup of numerical simulations. In the framework of evolutionary problems, space-time parallel methods…
In the Schrodinger picture of the Dirac quantum mechanics, defined in charts with spatially flat Robertson-Walker metrics and Cartesian coordinates the perturbation theory is applied to the interacting part of the Hamiltonian operator…
The interval approach to computation of dynamics of celestial bodies in the planetary problem has been considered. It is based on the refusal from idealization of infinitely high resolving capacity of measuring tools, and forms an…
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 introduce a novel numerical method to integrate partial differential equations representing the Hamiltonian dynamics of field theories. It is a multi-symplectic integrator that locally conserves the stress-energy tensor with an excellent…
The explicit split-operator algorithm is often used for solving the linear and nonlinear time-dependent Schr\"{o}dinger equations. However, when applied to certain nonlinear time-dependent Schr\"{o}dinger equations, this algorithm 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…
To solve the spinor-spinor Bethe-Salpeter equation in Euclidean space we propose a novel method related to the use of hyperspherical harmonics. We suggest an appropriate extension to form a new basis of spin-angular harmonics that is…
We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps…
Splitting methods have emerged as powerful tools to address complex problems by decomposing them into smaller solvable components. In this work, we develop a general approach to forward-backward splitting methods for solving monotone…
There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we…