It is well-known that chaotic dynamic systems, e.g., three-body system and turbulent flow, have sensitive dependence on the initial conditions(SDIC). Unfortunately,numerical noises, i.e., truncation error and round-of...It is well-known that chaotic dynamic systems, e.g., three-body system and turbulent flow, have sensitive dependence on the initial conditions(SDIC). Unfortunately,numerical noises, i.e., truncation error and round-off error, always exist in practice. Thus,due to the SDIC, the long-term accurate prediction of chaotic dynamic systems is practically impossible. In this paper, a new strategy for chaotic dynamic systems, i.e., the clean numerical simulation(CNS), is briefly described, and applied to a few Hamiltonian chaotic systems. With negligible numerical noises, the CNS can provide convergent(reliable) chaotic trajectories in a long enough interval of time. This is very important for Hamiltonian systems, and thus should have many applications in various fields. It is found that the traditional numerical methods in double precision cannot give not only reliable trajectories but also reliable Fourier power spectra and autocorrelation functions(ACFs). In addition, even the statistic properties of chaotic systems cannot be correctly obtained by means of traditional numerical algorithms in double precision, as long as these statistics are time-dependent. The CNS results strongly suggest that one had better be very careful on the direct numerical simulation(DNS) results of statistically unsteady turbulent flows, although DNS results often agree well with experimental data when the turbulent flow is in a statistical stationary state.展开更多
Turbulence is strongly associated with the vast majority of fluid flows in nature and industry.Traditionally,results given by the direct numerical simulation(DNS)of Navier-Stokes(NS)equations that relate to a famous m...Turbulence is strongly associated with the vast majority of fluid flows in nature and industry.Traditionally,results given by the direct numerical simulation(DNS)of Navier-Stokes(NS)equations that relate to a famous millennium problem are widely regarded as‘reliable’benchmark solutions of turbulence,as long as grid spacing is fine enough(i.e.less than the minimum Kolmogorov scale)and time-step is small enough,say,satisfying the Courant-Friedrichs-Lewy condition(Courant number<1).Is this really true?In this paper a two-dimensional sustained turbulent Kolmogorov flow driven by an external body force governed by the NS equations under an initial condition with a spatial symmetry is investigated numerically by the two numerical methods with detailed comparisons:one is the traditional DNS,the other is the‘clean numerical simulation’(CNS).In theory,the exact solution must have a kind of spatial symmetry since its initial condition is spatially symmetric.However,it is found that numerical noises of the DNS are quickly enlarged to the same level as the‘true’physical solution,which finally destroy the spatial symmetry of the flow field.In other words,the DNS results of the turbulent Kolmogorov flow governed by the NS equations are badly polluted mostly.On the contrary,the numerical noise of the CNS is much smaller than the‘true’physical solution of turbulence in a long enough interval of time so that the CNS result is very close to the‘true’physical solution and thus can remain symmetric,which can be used as a benchmark solution for comparison.Besides,it is found that numerical noise as a kind of artificial tiny disturbances can lead to huge deviations at large scale on the two-dimensional Kolmogorov turbulence governed by the NS equations,not only quantitatively(even in statistics)but also qualitatively(such as spatial symmetry of flow).This highly suggests that fine enough spatial grid spacing with small enough time-step alone could not guarantee the validity of the DNS of the NS equations:it is only a necessary condition but not sufficient.All of these findings might challenge some of our general beliefs in turbulence:for example,it might be wrong in physics to neglect the influences of small disturbances to NS equations.Our results suggest that,from physical point of view,it should be better to use the Landau-Lifshitz-Navier-Stokes(LLNS)equations,which consider the influence of unavoidable thermal fluctuations,instead of the NS equations,to model turbulent flows.展开更多
基金Project supported by the National Natural Science Foundation of China(No.91752104)
文摘It is well-known that chaotic dynamic systems, e.g., three-body system and turbulent flow, have sensitive dependence on the initial conditions(SDIC). Unfortunately,numerical noises, i.e., truncation error and round-off error, always exist in practice. Thus,due to the SDIC, the long-term accurate prediction of chaotic dynamic systems is practically impossible. In this paper, a new strategy for chaotic dynamic systems, i.e., the clean numerical simulation(CNS), is briefly described, and applied to a few Hamiltonian chaotic systems. With negligible numerical noises, the CNS can provide convergent(reliable) chaotic trajectories in a long enough interval of time. This is very important for Hamiltonian systems, and thus should have many applications in various fields. It is found that the traditional numerical methods in double precision cannot give not only reliable trajectories but also reliable Fourier power spectra and autocorrelation functions(ACFs). In addition, even the statistic properties of chaotic systems cannot be correctly obtained by means of traditional numerical algorithms in double precision, as long as these statistics are time-dependent. The CNS results strongly suggest that one had better be very careful on the direct numerical simulation(DNS) results of statistically unsteady turbulent flows, although DNS results often agree well with experimental data when the turbulent flow is in a statistical stationary state.
文摘Turbulence is strongly associated with the vast majority of fluid flows in nature and industry.Traditionally,results given by the direct numerical simulation(DNS)of Navier-Stokes(NS)equations that relate to a famous millennium problem are widely regarded as‘reliable’benchmark solutions of turbulence,as long as grid spacing is fine enough(i.e.less than the minimum Kolmogorov scale)and time-step is small enough,say,satisfying the Courant-Friedrichs-Lewy condition(Courant number<1).Is this really true?In this paper a two-dimensional sustained turbulent Kolmogorov flow driven by an external body force governed by the NS equations under an initial condition with a spatial symmetry is investigated numerically by the two numerical methods with detailed comparisons:one is the traditional DNS,the other is the‘clean numerical simulation’(CNS).In theory,the exact solution must have a kind of spatial symmetry since its initial condition is spatially symmetric.However,it is found that numerical noises of the DNS are quickly enlarged to the same level as the‘true’physical solution,which finally destroy the spatial symmetry of the flow field.In other words,the DNS results of the turbulent Kolmogorov flow governed by the NS equations are badly polluted mostly.On the contrary,the numerical noise of the CNS is much smaller than the‘true’physical solution of turbulence in a long enough interval of time so that the CNS result is very close to the‘true’physical solution and thus can remain symmetric,which can be used as a benchmark solution for comparison.Besides,it is found that numerical noise as a kind of artificial tiny disturbances can lead to huge deviations at large scale on the two-dimensional Kolmogorov turbulence governed by the NS equations,not only quantitatively(even in statistics)but also qualitatively(such as spatial symmetry of flow).This highly suggests that fine enough spatial grid spacing with small enough time-step alone could not guarantee the validity of the DNS of the NS equations:it is only a necessary condition but not sufficient.All of these findings might challenge some of our general beliefs in turbulence:for example,it might be wrong in physics to neglect the influences of small disturbances to NS equations.Our results suggest that,from physical point of view,it should be better to use the Landau-Lifshitz-Navier-Stokes(LLNS)equations,which consider the influence of unavoidable thermal fluctuations,instead of the NS equations,to model turbulent flows.
