Parametric excitation of computational mode of the leapfrog scheme applied to the Van der Pol equation (Q1261367)

The leapfrog scheme is applied to the Van der Pol equation. When the amplitude of oscillation of the physical mode exceeds a critical value, the computational mode is parametrically excited by the physical mode. The growth of the computational mode interrupts the integration based on the leapfrog scheme. The critical amplitude of the physical mode is determined by the linear stability analysis and the parametric excitation theory. The Runge-Kutta smoother eliminating the computational mode enables the longtime integration based on the leapfrog scheme.