On the risks of using double precision in numerical simulations of spatio-temporal chaos

Abstract

Due to the butterfly-effect, computer-generated chaotic simulations often deviate exponentially from the true solution, so that it is very hard to obtain a reliable simulation of chaos in a long-duration time. In this paper, a new strategy of the so-called Clean Numerical Simulation (CNS) in physical space is proposed for spatio-temporal chaos, which is computationally much more efficient than its predecessor (in spectral space). The strategy of the CNS is to reduce both of the truncation and round-off errors to a specified level by implementing high-order algorithms in multiple-precision arithmetic (with sufficient significant digits for all variables and parameters) so as to guarantee that numerical noise is below such a critical level in a temporal interval $t\in [0,T_c]$ that corresponding numerical simulation remains reliable over the whole interval. Without loss of generality, the complex Ginzburg-Landau equation (CGLE) is used to illustrate its validity. As a result, a reliable long-duration numerical simulation of the CGLE is achieved in the whole spatial domain over a long interval of time $t\in [0,3000]$, which is used as a reliable benchmark solution to investigate the influence of numerical noise by comparing it with the corresponding ones given by the 4th-order Runge-Kutta method in double precision (RKwD). Our results demonstrate that the use of double precision in simulations of chaos might lead to huge errors in the prediction of spatio-temporal trajectories and in statistics, not only quantitatively but also qualitatively, particularly in a long interval of time.

37 pages, 21 figures

References in corpus (5)

• Deep Learning of Vortex Induced Vibrations
• Newton vs the machine: solving the chaotic three-body problem using deep neural networks
• The 1223 new periodic orbits of planar three-body problem with unequal mass and zero angular momentum
• Collisionless periodic orbits in the free-fall three-body problem
• On the origin of intrinsic randomness of Rayleigh-Benard turbulence

Cited by in corpus (15)

• Three-body problem -- from Newton to supercomputer plus machine learning
• Large-scale influence of numerical noises as artificial stochastic disturbances on a sustained turbulence
• On the probability of finding a nonphysical solution through shadowing
• Ultra-chaos: an insurmountable objective obstacle of reproducibility and replicability
• Accurate predictions of chaotic motion of a free fall disk
• A kind of Lagrangian chaotic property of the Arnold-Beltrami-Childress flow
• Is a direct numerical simulation (DNS) of Navier-Stokes equations with small enough grid spacing and time-step definitely reliable/correct?
• One family of 13315 stable periodic orbits of the non-hierarchical unequal-mass triple system
• A Self-Adaptive Algorithm of the Clean Numerical Simulation (CNS) for Chaos
• Noise-expansion cascade: an origin of randomness of turbulence
• Clean numerical simulation (CNS) of three-dimensional turbulent Kolmogorov flow
• Model-free Forecasting of Rogue Waves using Reservoir Computing
• An efficient approach for searching three-body periodic orbits passing through Eulerian configuration
• Physical significance of artificial numerical noise in direct numerical simulation of turbulence
• How Far Can We Trust Chaos? Extending the Horizon of Predictability