Poincaré Map Based on Splitting Methods

Firstly, by using the Liouville formula, we prove that the Jacobian matrix determinants of splitting methods are equal to that of the exact flow. However, for the explicit Runge--Kutta methods, there is an error term of order p+1 for the Jacobian matrix determinants. Then, the volume evolution law of a given region in phase space is discussed for different algorithms. It is proved that splitting methods can exactly preserve the sum of Lyapunov exponents invariable. Finally, a Poincaré map and its energy distribution of the Duffing equation are computed by using the second-order splitting method and the Heun method (a second-order Runge--Kutta method). Computation illustrates that the results by splitting methods can properly represent systems' chaotic phenomena.