Clean numerical simulation (CNS) of three-dimensional turbulent Kolmogorov flow

Abstract: Turbulence holds immense importance across various scientific and engineering disciplines. The direct numerical simulation (DNS) of turbulence proposed by Orszag in 1970 is a milestone in fluid mechanics, which began an era of numerical experiment for turbulence. Many researchers have reported that turbulence should be chaotic, since spatiotemporal trajectories are very sensitive to small disturbance. In 2006 E.D. Lorenz discovered that numerical noise might have a great influence on statistic characteristics of chaotic systems. The above-mentioned facts logically lead to the conclusion that numerical noises of turbulence (as a chaotic system) might have great influence on turbulent flows. This is indeed true for a two-dimensional (2D) Kolmogorov turbulent flow, as currently revealed by a much more accurate algorithm than DNS, namely the ``clean numerical simulation'' (CNS). Different from DNS, CNS can greatly reduce both of truncation error and round-off error to any required small level so that artificial numerical noise can be rigorously negligible throughout a time interval long enough for calculating statistics. However, In physics, 3D turbulent flow is more important than 2D turbulence, as pointed out by Nobel Prize winner T.D. Lee. So, for the first time, we solve a 3D turbulent Kolmogorov flow by means of CNS in this paper, and compare our CNS result with that given by DNS in details. It is found that the spatial-temporal trajectories of the 3D Kolmogorov turbulent flow given by DNS are badly polluted by artificial numerical noise rather quickly, and the DNS result has huge deviations from the CNS benchmark solution not only in the flow field and the energy cascade but also even in statistics.