Leapfrog integration

In mathematics leapfrog integration is a method for numerically integrating differential equations of the form

{\ddot {x}}=d^{2}x/dt^{2}=F(x)

or equivalently of the form

{\dot {v}}=dv/dt=F(x),\;{\dot {x}}=dx/dt=v

particularly in the case of a dynamical system of classical mechanics.

The method is known by different names in different disciplines. In particular, it is similar to the velocity Verlet method, which is a variant of Verlet integration. Leapfrog integration is equivalent to updating positions $x(t)$ and velocities $v(t)={\dot {x}}(t)$ at interleaved time points, staggered in such a way that they "leapfrog" over each other.

Leapfrog integration is a second-order method, in contrast to Euler integration, which is only first-order, yet requires the same number of function evaluations per step. Unlike Euler integration, it is stable for oscillatory motion, as long as the time-step $\Delta t$ is constant, and $\Delta t\leq 2/\omega$ .^[1]

In leapfrog integration, the equations for updating position and velocity are

{\begin{aligned}x_{i}&=x_{i-1}+v_{i-1/2}\,\Delta t,\\a_{i}&=F(x_{i}),\\v_{i+1/2}&=v_{i-1/2}+a_{i}\,\Delta t,\end{aligned}}

where $x_{i}$ is position at step $i$ , $v_{i+1/2\,}$ is the velocity, or first derivative of $x$ , at step $i+1/2\,$ , $a_{i}=F(x_{i})$ is the acceleration, or second derivative of $x$ , at step $i$ , and $\Delta t$ is the size of each time step. These equations can be expressed in a form that gives velocity at integer steps as well:^[2]

{\begin{aligned}x_{i+1}&=x_{i}+v_{i}\,\Delta t+{\tfrac {1}{2}}\,a_{i}\,\Delta t^{\,2},\\v_{i+1}&=v_{i}+{\tfrac {1}{2}}(a_{i}+a_{i+1})\,\Delta t.\end{aligned}}

However, even in this synchronized form, the time-step $\Delta t$ must be constant to maintain stability.^[3]

One use of this equation is in gravity simulations, since in that case the acceleration depends only on the positions of the gravitating masses (and not on their velocities), although higher-order integrators (such as Runge–Kutta methods) are more frequently used.

There are two primary strengths to leapfrog integration when applied to mechanics problems. The first is the time-reversibility of the Leapfrog method. One can integrate forward n steps, and then reverse the direction of integration and integrate backwards n steps to arrive at the same starting position. The second strength is its symplectic nature, which implies that it conserves the (slightly modified) energy of dynamical systems. This is especially useful when computing orbital dynamics, as many other integration schemes, such as the (order-4) Runge-Kutta method, do not conserve energy and allow the system to drift substantially over time.

Because of its time-reversibility, and because it is a symplectic integrator, leapfrog integration is also used in Hamiltonian Monte Carlo, a method for drawing random samples from a probability distribution whose overall normalization is unknown.^[4]

References

↑ C. K. Birdsall and A. B. Langdon, Plasma Physics via Computer Simulations, McGraw-Hill Book Company, 1985, p. 56.
↑ 4.1 Two Ways to Write the Leapfrog
↑ [Skeel, R. D., "Variable Step Size Destabilizes the Stömer/Leapfrog/Verlet Method", BIT Numerical Mathematics, Vol. 33, 1993, p. 172–175.]
↑ Bishop, Christopher (2006). Pattern Recognition and Machine Learning. New York: Springer-Verlag. pp. 548–554. ISBN 978-0-387-31073-2.

External links

The Leapfrog Integrator, Drexel University Physics

Numerical methods for integration

First-order methods	Euler method Backward Euler Semi-implicit Euler Exponential Euler

Second-order methods	Verlet integration Velocity Verlet Trapezoidal rule Beeman's algorithm Midpoint method Heun's method Newmark-beta method Leapfrog integration

Higher-order methods	Exponential integrators Runge–Kutta methods List of Runge–Kutta methods Linear multistep method General linear methods Backward differentiation formula

This article is issued from Wikipedia - version of the 6/3/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.

Leapfrog integration

See also

References

External links