Symplectic integrators in the shearing sheet
نویسندگان
چکیده
منابع مشابه
The accuracy of symplectic integrators
We judge symplectic integrators by the accuracy with which they represent the Hamil-tonian function. This accuracy is computed, compared and tested for several diierent methods. We develop new, highly accurate explicit fourth-and fth-order methods valid when the Hamiltonian is separable with quadratic kinetic energy. For the near-integrable case, we connrm several of their properties expected f...
متن کاملSymplectic Integrators in Numerical Relativity
The purpose of this note is to point out that a naive application of symplectic integration schemes such as SHAKE or RATTLE which preserve holonomic constraints encounters difficulties when applied to the numerical treatment of the equations of general relativity. It is well known that the equations of General Relativity (GR) can be derived from a variational principle and that they can be cast...
متن کاملNonaxisymmetric stability in the shearing sheet approximation
Aims. To quantify the transient growth of nonaxisymmetric perturbations in unstratified magnetized and stratified non-magnetized rotating linear shear flows in the shearing sheet approximation of accretion disc flows. Methods. The Rayleigh quotient in modal approaches for the linearized equations (with time-dependent wavenumber) and the amplitudes from direct shearing sheet simulations using a ...
متن کاملSymplectic integrators for spin systems.
We present a symplectic integrator, based on the implicit midpoint method, for classical spin systems where each spin is a unit vector in R{3}. Unlike splitting methods, it is defined for all Hamiltonians and is O(3)-equivariant, i.e., coordinate-independent. It is a rare example of a generating function for symplectic maps of a noncanonical phase space. It yields a new integrable discretizatio...
متن کاملSymplectic integrators from composite operator
I derive fourth order symplectic integrators by factorizing the exponential of two operators in terms of an additional higher order composite operator with positive coefficients. One algorithm requires only one evaluation of the force and one evaluation of the force and its gradient. When applied to Kepler’s problem, the energy error function associated with these algorithms are approximately 1...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Monthly Notices of the Royal Astronomical Society
سال: 2011
ISSN: 0035-8711
DOI: 10.1111/j.1365-2966.2011.18939.x