Skip Navigation Links
Discontinuity, Nonlinearity, and Complexity

Dimitry Volchenkov (editor), Dumitru Baleanu (editor)

Dimitry Volchenkov(editor)

Mathematics & Statistics, Texas Tech University, 1108 Memorial Circle, Lubbock, TX 79409, USA


Dumitru Baleanu (editor)

Cankaya University, Ankara, Turkey; Institute of Space Sciences, Magurele-Bucharest, Romania


Parallel Computation of Reliable Chaotic Solutions of Saltzman’s Equations by Means of the Clean Numerical Simulation

Discontinuity, Nonlinearity, and Complexity 2(4) (2013) 345--355 | DOI:10.5890/DNC.2013.11.004

Peng Yang$^{2}$, Zhiliang Lin$^{2}$, Shijun Liao$^{1}$,$^{2}$,$^{3}$,$^{4}$

$^{1}$ State Key Laboratory of Ocean Engineering

$^{2}$ School of Naval Architecture, Ocean and Civil Engineering

$^{3}$ Department of Mathematics, Shanghai Jiao Tong University, Shanghai 200240, China

$^{4}$ Nonlinear Analysis and Applied Mathematics Research Group (NAAM), King Abdulaziz University (KAU), Jeddah, Saudi Arabia

Download Full Text PDF



The method of the so-called “Clean Numerical Simulation” (CNS) is applied to gain reliable chaotic solutions of Saltzman’s dynamic system, a simplified model for convection flows of fluid. Based on the high-order Taylor series method with data in multiple precision library and a validation check of global reliability of result, the CNS provides us a practical way to gain reliable, accurate enough solutions of chaotic dynamic systems in a finite but long enough time interval. The parallel computation is used to greatly increase the computational efficiency. The numerical noises of the CNS can be controlled to be so small that even the influence of the micro-level inherent uncertainty of initial conditions can be investigated in details. It is found that the micro-level inherent physical uncertainty (i.e. the unavoidable statistical fluctuation of temperature and velocity of fluid) of initial conditions of chaotic Saltzman’s system transfers into macroscopic randomness. This suggests that chaos might be a bridge between micro-level inherent physical uncertainty and macroscopic randomness. The current work illustrates that the above conclusion holds not only for Lorenz equation with three ODEs but also for Saltzman’s equation with up-to nine ODEs, and thus has general meanings.


Thanks to Dr. P.F. Wang for his kindly assistance on multiple precision library and parallel computation. This work is partly supported by the State Key Lab of Ocean Engineering (Approval No. GKZD010056-6) and the National Natural Science Foundation of China (Approval No. 11272209).


  1. [1]  Lorenz, E.N. (1963), Deterministic nonperiodic flow, Journal of the Atmospheric Sciences, 20, 130-141.
  2. [2]  Lorenz, E.N. (1989), Computational chaos: a prelude to computational instability, Physica D, 35, 299-317.
  3. [3]  Lorenz, E.N. (2006), Computational periodicity as observed in a simple system, Tellus-A, 58, 549-557.
  4. [4]  Yao, L.S. and Hughes, D. (2008), Comment on 'Computational periodicity as observed in a simple system' By Edward N. Lorenz (2006), Tellus-A, 60, 803-805.
  5. [5]  Barrio, R. (2005), Performance of the Taylor series method for ODEs/DAEs, Applied Mathematics and Computation , 163, 525-545.
  6. [6]  Barrio, R., Rodriguez, M., Abad, A., and Blesa, F. (2011), Breaking the limits: the Taylor series method, Applied Mathematics and Computation, 217, 7940-7954.
  7. [7]  Corliss, G. and Chang, Y.F. (1982), Solving ordinary differential equations using Taylor series, Association for Computing Machinery Transactions on Mathematical Software , 8, 114-144.
  8. [8]  Deprit, A. and Zahar, R. (1966), Numerical integration of an orbit and its concomitant variations, Z. Angew. Math. Phys., 17, 425-430.
  9. [9]  Nedialkov, N.S. and Pryce, J.D. (2005), Solving differential-algebraic equations by Taylor series (I): computing Taylot coefficients, BIT, 45, 561-591.
  10. [10]  Liao, S.J. (2009), On the reliability of computed chaotic solution of non-linear differential equations, Tellus-A, 61, 550-564.
  11. [11]  Liao, S.J. (2012), Chaos: A bridge from microscopic uncertainty to macroscopic randomness, Communications in Nonlinear Science and Numerical Simulation, 17, 2564-2569.
  12. [12]  Liao, S.J. (2013), On the numerical simulation of propagation of micro-level uncertainty for chaotic dynamic systems, Chaos, Solitons and Fractals, 47, 1-12.
  13. [13]  Wang, P.F., Li, J.P., and Li, Q. (2012), Computational uncertainty and the application of a high-performance multiple precision scheme to obtaining the correct reference solution of lorenz equations, Numerical Algorithms, 59, 147-159.
  14. [14]  Sarra, S.A. and Meador, C. (2011), On the numerical simulation ofchaotic dynamic systems using extend precision floating point arithmetic and very high order numerical methods, Nonlinear Analysis - Modelling and Control, 16, 340-352.
  15. [15]  Saltzman, B. (1962), Finite amplitude free convection as an initial value problem(1), J. Atmos. Sci, 19, 329-341.
  16. [16]  Nikolaevskaya, E., Khimich, A., and Chistyakova, T. (2012), Programming with Multiple Precision, Springer.