Numerical methods for Chaotic ODE

TL;DR

This paper investigates backward error analysis methods for numerical solutions of chaotic ODEs, examining how they explain the accuracy of long-term statistical behavior despite numerical challenges.

Contribution

It compares three approaches—residual assessment, modified equations, and shadowing—in analyzing chaotic systems and highlights an open problem in understanding statistical accuracy.

Findings

Methods explain why numerical simulations capture chaotic behavior.

Numerical methods can introduce spurious chaos or suppress chaos.

Open problem: why long-term statistics are often correct despite lack of guarantees.

Abstract

This paper explores backward error analysis for numerical solutions of ordinary differential equations, particularly focusing on chaotic systems. Three approaches are examined: residual assessment, the method of modified equations, and shadowing. We investigate how these methods explain the success of numerical simulations in capturing the behavior of chaotic systems, even when facing issues like spurious chaos introduced by numerical methods or suppression of chaos by numerical methods. Finally, we point out an open problem, namely to explain why the statistics of long orbits are usually correct, even though we do not have a theoretical guarantee why this should be so.