Related papers: A Symplectic Integrator for Hill's Equations
Reliable studies of the long-term dynamics of planetary systems require numerical integrators that are accurate and fast. The challenge is often formidable because the chaotic nature of many systems requires relative numerical error bounds…
We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body…
We introduce Spiral, a third-order integration algorithm for the rotational motion of extended bodies. It requires only one force calculation per time step, does not require quaternion normalization at each time step, and can be formulated…
The mathematical theory underlying Hamiltonian mechanics is called symplectic geometry. So symplectic geometry arose from the roots of mechanics and is seen as one of the most valuable links between physics and mathematics today. Symplectic…
This article is a survey of classical and quantum completely integrable systems from the viewpoint of local ``phase space'' analysis. It advocates the use of normal forms and shows how to get global information from glueing local pieces.…
We study the spatial isosceles three-body problem from the perspective of Symplectic Dynamics. For certain choices of mass ratio, angular momentum, and energy, the dynamics on the energy surface is equivalent to a Reeb flow on the tight…
We present a brief tutorial on the nuts and bolts computation of a multisymplectic particle-in-cell algorithm using the discretized Lagrangian approach. This approach, originated by Marsden, Shadwick, and others, brings the benefits of…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
The paper intends to lay out the first steps towards constructing a unified framework to understand the symplectic and spectral theory of finite dimensional integrable Hamiltonian systems. While it is difficult to know what the best…
In this paper, we propose a numerical method to solve isotropic elliptic equations on point cloud by generalizing the point integral method. The idea of the point integral method is to approximate the differential operators by integral…
Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in…
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…
In this work, we introduce a spectral-infinite element method for solving Einstein's constraint equations in hyperbolic form. As an application of this, we use this method for computing asymptotically flat perturbations of a Kerr black hole…
The integration of the equations of motion in gravitational dynamical systems -- either in our Solar System or for extra-solar planetary system -- being non integrable in the global case, is usually performed by means of numerical…
We report a few sumerical tests comparing some newly defined energy-preserving integrators and symplectic methods, using either constant and variable stepsize.
We develop a unified approach for construction of symplectic forms for 1D integrable equations with the periodic and rapidly decaying initial data. As an example we consider the cubic nonlinear Schr\"{o}dinger equation.
In this paper an approach is outlined. With this approach some explicit algorithms can be applied to solve the initial value problem of $n-$dimensional damped oscillators. This approach is based upon following structure: for any…
This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show…
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…
We compute the normal forms for the Hamiltonian leading to the epicyclic approximations of the (perturbed) Kepler problem in the plane. The Hamiltonian setting corresponds to the dynamics in the Hill synodic system where, by means of the…