Semi-implicit Euler method

In mathematics, the semi-implicit Euler method, also called symplectic Euler, semi-explicit Euler, Euler–Cromer, and Newton–Størmer–Verlet (NSV), is a modification of the Euler method for solving Hamilton's equations, a system of ordinary differential equations that arises in classical mechanics. It is a symplectic integrator and hence it yields better results than the standard Euler method.

In mathematics, the semi-implicit Euler method, also called symplectic Euler, semi-explicit Euler, Euler–Cromer, and Newton–Størmer–Verlet (NSV), is a modification of the Euler method for solving Hamilton's equations, a system of ordinary differential equations that arises in classical mechanics. It is a symplectic integrator and hence it yields better results than the standard Euler method. The semi-implicit Euler method can be applied to a pair of differential equations of the form \begin{align} \frac{dx}{dt} &= f(t,v) \\ \frac{dv}{dt} &= g(t,x), \end{align} where f and g are given functions. Here, x and v may be either scalars or vectors. The equations of motion in Hamiltonian mechanics take this form if the Hamiltonian is of the form H = T(t,v) + V(t,x). \, The differential equations are to be solved with the initial condition x(t_0) = x_0, \qquad v(t_0) = v_0. \begin{align} v_{n+1} &= v_n + g(t_n, x_n) \, \Delta t\\[0.3em] x_{n+1} &= x_n + f(t_n, v_{n+1}) \, \Delta t \end{align} where \Delta t is the time step and t_n = t_0 + n\Delta t is the time after n steps. The difference with the standard Euler method is that the semi-implicit Euler method uses v_{n+1} in the equation for x_{n+1} , while the Euler method uses v_n . Applying the method with negative time step to the computation of (x_n, v_n) from (x_{n+1}, v_{n+1}) and rearranging leads to the second variant of the semi-implicit Euler method \begin{align} x_{n+1} &= x_n + f(t_n, v_n) \, \Delta t\\[0.3ex] v_{n+1} &= v_n + g(t_n, x_{n+1}) \, \Delta t \end{align} which has similar properties. The semi-implicit Euler is a first-order integrator, just as the standard Euler method. This means that it commits a global error of the order of Δt. However, the semi-implicit Euler method is a symplectic integrator, unlike the standard method. As a consequence, the semi-implicit Euler method almost conserves the energy (when the Hamiltonian is time-independent). Often, the energy increases steadily when the standard Euler method is applied, making it far less accurate. Alternating between the two variants of the semi-implicit Euler method leads in one simplification to the Störmer-Verlet integration and in a slightly different simplification to the leapfrog integration, increasing both the order of the error and the order of preservation of energy. The stability region of the semi-implicit method was presented by Niiranen although the semi-implicit Euler was misleadingly called symmetric Euler in his paper. The semi-implicit method models the simulated system correctly if the complex roots of the characteristic equation are within the circle defined by \left(x-\frac{2}{t_\Delta}\right)^2 + y^2 which can also be described as a circle centered at (2/t_\Delta,0) with radius 2\sqrt{2}/t_\Delta. This circle is shown below. If the roots of the characteristic polynomial are real numbers s_1, s_2 , then the stability region can extend outside the circle. The criterion for this is given by \begin{cases} s_2 - 2/\Delta t \\ s_2 > \frac{2s_1+\frac{4}{t_\Delta}}{s_1t_\Delta-2} \quad \text{ when } s_1 As an example, if s_1 approaches 0, then the simulation will be stable for s_2 > -2/t_\Delta, which is also the same limit of stability as for the explicit Euler method. As can be seen, the semi-implicit method can simulate correctly both stable systems that have their roots in the left half plane and unstable systems that have their roots in the right half plane. This is clear advantage over forward (standard) Euler and backward Euler. Forward Euler tends to have less damping than the real system when the negative real parts of the roots get near the imaginary axis and backward Euler may show the system be stable even when the roots are in the right half plane. The motion of a spring satisfying Hooke's law is given by \begin{align} \frac{dx}{dt} &= v(t)\\[0.2em] \frac{dv}{dt} &= -\frac{k}{m}\,x=-\omega^2\,x. \end{align} The semi-implicit Euler for this equation is \begin{align} v_{n+1} &= v_n - \omega^2\,x_n\,\Delta t \\[0.2em] x_{n+1} &= x_n + v_{n+1} \,\Delta t. \end{align} Substituting v_{n+1} in the second equation with the expression given by the first equation, the iteration can be expressed in the following matrix form \begin{bmatrix} x_{n+1} \\v_{n+1}\end{bmatrix} = \begin{bmatrix} 1-\omega^2 \Delta t^2 & \Delta t \\ -\omega^2 \Delta t & 1 \end{bmatrix} \begin{bmatrix} x_{n} \\ v_{n} \end{bmatrix}, and since the determinant of the matrix is 1 the transformation is area-preserving. The iteration preserves the modified energy functional E_h(x,v) = \tfrac{1}{2} \left(v^2+\omega^2\,x^2 - \omega^2\Delta t\,vx\right) exactly, leading to stable periodic orbits (for sufficiently small step size) that deviate by O(\Delta t) from the exact orbits. The exact circular frequency \omega increases in the numerical approximation by a factor of 1+\tfrac1{24}\omega^2\Delta t^2+O(\Delta t^4) .

Text derived from https://en.wikipedia.org/w/index.php?oldid=1378314259 under CC BY-SA 4.0; modifications have been made.