On the susceptibility of numerical methods to computational chaos and superstability

Abstract In the present study, the susceptibility of the forward and the backward Euler methods to computational chaos and superstability is investigated via the means of both a theoretical analysis and numerical experiments. A linear stability analysis of the fixed points and the periodic orbits of the maps induced by these methods asserts that, for large enough time-steps Δt , these maps undergo bifurcations and as result the acquired solutions are spurious. More specifically, it is shown that the backward Euler method suppresses chaotic behavior, whereas the forward Euler renders all linearly stable fixed points and periodic orbits of its induced map linearly unstable. Numerical experiments that illustrate the validity of the theoretical analysis are also presented and discussed. For the forward Euler method, in particular, the computation of bifurcation diagrams, the Maximum Lyapunov exponent and the Kolmogorov–Sinai entropy suggest that it can engender computational chaos.

[1]  Y. Pesin.  CHARACTERISTIC LYAPUNOV EXPONENTS AND SMOOTH ERGODIC THEORY , 1977 .

[2]  Robert M Corless,et al.  Numerical methods can suppress chaos , 1991 .

[3]  José Manuel Gutiérrez,et al.  Chaos suppression through changes in the system variables and numerical rounding errors , 1996 .

[4]  Edward N. Lorenz,et al.  Computational chaos-a prelude to computational instability , 1989 .

[5]  P. Grassberger,et al.  Estimation of the Kolmogorov entropy from a chaotic signal , 1983 .

[6]  J. Butcher.  The Numerical Analysis of Ordinary Di erential Equa-tions , 1986 .

[7]  H. C. Yee,et al.  Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I. The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics☆ , 1991 .

[8]  Shigehiro Ushiki,et al.  Central difference scheme and chaos , 1982 .

[9]  de Markus,et al.  Detection of the onset of numerical chaotic instabilities by lyapunov exponents , 2001 .

[10]  Mark J. Ablowitz,et al.  On homoclinic structure and numerically induced chaos for the nonlinear Schro¨dinger equation , 1990 .

[11]  Shigehiro Ushiki,et al.  Chaos in numerical analysis of ordinary differential equations , 1981 .

[12]  F. Vadillo,et al.  On Spurious Fixed Points of Runge-Kutta Methods , 1997 .

[13]  Robert M. Corless,et al.  Defect-controlled numerical methods and shadowing for chaotic differential equations , 1992 .

[14]  H. C. Yee,et al.  On spurious asymptotic numerical solutions of explicit Runge-Kutta methods , 1992 .

[15]  Edward H. Twizell,et al.  Explicit finite-difference methods for non-linear dynamic systems: Froude's pendulum , 1996 .

[16]  W. Tucker.  The Lorenz attractor exists , 1999 .

[17]  Slawomir T. Fryska,et al.  Computer dynamics and shadowing of chaotic orbits , 1992 .

[18]  R. Meyer-Spasche,et al.  Euler-like discrete models of the logistic differential equation , 1998 .

[19]  Arieh Iserles,et al.  A unified approach to spurious solutions introduced by time discretization. Part I: basic theory , 1991 .

[20]  Joseph D. Skufca,et al.  Analysis Still Matters: A Surprising Instance of Failure of Runge-Kutta-Felberg ODE Solvers , 2004, SIAM Rev..

[21]  Lennart Stenflo,et al.  Generalized Lorenz equations for acoustic-gravity waves in the atmosphere , 1996 .

[22]  Robert M. Corless,et al.  What good are numerical simulations of chaotic dynamical systems , 1994 .

[23]  Yoshitaka Saiki,et al.  Numerical detection of unstable periodic orbits in continuous-time dynamical systems with chaotic behaviors , 2007 .

[24]  L. Stenflo.  Equations describing solitary atmospheric waves , 1991 .

[25]  Germund Dahlquist,et al.  Are the numerical methods and software satisfactory for chemical kinetics , 1982 .

[26]  Alicia Serfaty de Markus.  Scaling Laws in numerical Instabilities induced by Fixed-Step Integration of Dynamical Systems , 2001, Int. J. Bifurc. Chaos.

[27]  E. Lorenz.  Deterministic nonperiodic flow , 1963 .

[28]  Ronald E. Mickens,et al.  Suppression of numerically induced chaos with nonstandard finite difference schemes , 1999 .

[29]  Arieh Iserles,et al.  Unified approach to spurious solutions introduced by time discretization Part II: BDF-like methods , 1992 .

[30]  I. Prigogine,et al.  Symmetry Breaking Instabilities in Dissipative Systems. II , 1968 .

[31]  H. Kantz.  A robust method to estimate the maximal Lyapunov exponent of a time series , 1994 .