Related papers: Time-Symmetrized Kustaanheimo-Stiefel Regularizati…
We present the Kustaanheimo-Stiefel (KS) regularization of the elliptic restricted three-body problem (ER3BP) at the secondary body $P_2$, and discuss its use to study a category of transits through its Hill's sphere (fast close…
We discuss the implementation of a new regular algorithm for simulation of the gravitational few-body problem. The algorithm uses components from earlier methods, including the chain structure, the logarithmic Hamiltonian, and the…
The Kustaanheimo-Stiefel (KS) transformation maps the non-linear and singular equations of motion of the three-dimensional Kepler problem to the linear and regular equations of a four-dimensional harmonic oscillator. It is used extensively…
The article introduces a new algorithm for solving a class ofequilibrium problems involving strongly pseudomonotone bifunctions with Lipschitz-type condition. We describe how to incorporate the proximal-like regularized technique with…
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…
We introduce a numerical method to integrate the stochastic Landau-Lifshitz-Gilbert equation in spherical coordinates for generic discretization schemes. This method conserves the magnetization modulus and ensures the approach to…
We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for…
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…
Numerical solutions to Newton's equations of motion for chaotic self gravitating systems of more than 2 bodies are often regarded to be irreversible. This is due to the exponential growth of errors introduced by the integration scheme and…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms…
The time-symmetric block time--step (TSBTS) algorithm is a newly developed efficient scheme for $N$--body integrations. It is constructed on an era-based iteration. In this work, we re-designed the TSBTS integration scheme with dynamically…
In this paper, we study a second-order accurate and linear numerical scheme for the nonlocal Cahn-Hilliard equation. The scheme is established by combining a modified Crank-Nicolson approximation and the Adams-Bashforth extrapolation for…
In this note we approach the classical, Newtonian, gravitational $N$-body problem by mean of a new, original numerical integration method. After a short summary of the fundamental characteristics of the problem, including a sketch of some…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
This work considers the {\em gravitational} $N$-body problem and introduces global time-renormalization {\em functions} that allow the efficient numerical integration with fixed time-steps. First, a lower bound of the radius of convergence…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
This is the first in a series of papers devoted to fully general-relativistic $N$-body simulations applied to late-time cosmology. The purpose of this paper is to present the combination of a numerical relativity scheme, discretization…
We first discuss two-body and chain regularization methods for direct N-body simulations on HARP-2 and GRAPE-6. The former is used for accurate integration of perturbed binaries and hierarchies, whereas the latter deals with strong…
The importance and usefulness of renormalization are emphasized in nonrelativistic quantum mechanics. The momentum space treatment of both two-body bound state and scattering problems involving some potentials singular at the origin…