Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 2850651, 13 pages
http://dx.doi.org/10.1155/2016/2850651

Research Article

Low Model Analysis and Synchronous Simulation of the Wave Mechanics

Wenyuan Duan, Heyuan Wang, and Meng Kan

College of Science, Liaoning University of Technology, Jinzhou 121001, China

Received 11 May 2016; Accepted 13 July 2016

Academic Editor: Bernard Bonello

Copyright © 2016 Wenyuan Duan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The dynamic behavior of a chaotic system in the internal wave dynamics and the problem of the tracing and synchronization are investigated, and the numerical simulation is carried out in this paper. The globally exponentially attractive set and positive invariant set of the chaotic system are studied via constructing the positive definite and radial unbounded Lyapunov function. There are no equilibrium positions, periodic solutions, quasi-period motions, wandering recovering motions, and other chaotic attractors of the system out of the globally exponentially attractive set. Strange attractors can only locate in the globally exponentially attractive set. A feedback controller is designed for the chaotic system to realize the control of the unstable point. The second method of Lyapunov is used to discuss theoretically the rationality of the design of the controller. The driving-response synchronization method is used to realize the globally exponential synchronization. The numerical simulation is carried out by MATLAB software, and the simulation results show that the method is effective.

1. Introduction

Since 1960s, in the meteorological numerical research conducted by American meteorologist Lorenz, the three-dimensional truncation is carried out for the Rayleigh-Bernard heat convection problem of the infinite dimensional dynamical system to get the famous Lorenz system and the chaos phenomenon is found accidently [1]. A lot of researches on chaos are conducted by many scholars since the discovery of chaos [1–10]. Lorenz system reveals the essential phenomenon of the nonlinear complexity as the first chaotic system [8, 9]. The chaotic attractor was formally proposed by Ruelle and Takens in 1971 and the nonperiodic flow appearing in the dissipative system was called the strange attractor by them [1]. The capture area of the strange attractor and the global stability of the chaotic system are hot spots of common concern to the people. The low mode analysis is usually used to explain and predict the chaotic phenomenon of infinite dimensional dynamical system [10]. The inertial manifold and approximate inertial manifold theories are the theoretical foundation and basis of the low mode analysis (they are considered to be a low dimensional smooth manifold with exponent attracting all tracks containing the global attractor). The complex dynamic behavior of the infinite dimensional dynamical system usually originated from a simple origin and can be distinguished by the simple equation. The three functions of the internal gravity wave in the atmosphere are the following [11]: transmission and storage of energy and momentum, starting and organization of the convective activity, and occurrence and modulation of the turbulence. The thermal convection and turbulence caused by the internal wave are studied in essence in Lorenz’s paper [10]. The chaos caused by the match of the driving factor with the dissipation factor of the turbulence appears in the internal gravity wave dynamics. The dissipation factor is too weak causing the absence of chaos when the Prandtl number σ is equal to 1. The concept of chaos can help explain the occurrence of atmosphere turbulence. During the day, unstable temperature stratification is the driving factor of turbulence, and the molecular viscosity and certain velocity shear are dissipation factor. At night, effective surface long wave radiation is the driving factor of turbulence. The chaos can still appear under the condition that the real atmosphere Prandtl number is equal to 0.7 after introducing the velocity shear damping. The smaller the Richardson number is, the more likely the chaos appears at night [11]. Chaos is too sensitive to the initial value and the unpredictability of the trend of the long time of chaos exists, so the chaos system is considered to be uncontrollable and two chaos systems are considered to be synchronous impossibly by people for a long time. The OGY method was proposed since 1990 and the chaos synchronization [12–15] was used by the United States Navy laboratories to get the secure communication firstly in 1990 that changed people’s original idea. The control and synchronization of chaos has become a hot spot in the field of nonlinear science in recent years.

A new chaos system is found by Liu through numerical calculation in the study of the internal wave dynamics [11]. But he did not have a comprehensive analysis of the appearance and disappearance of chaotic attractor. Based on the research of Liu, we have carried out a comprehensive and detailed discussion and numerical simulation in this system; the corresponding attractors of the system are given; at the same time, we obtain the bifurcation diagram, power spectrum, Poincare section, and return mapping of the system. The general characteristics of the behavior of the chaotic system are detected via simulation, and we discuss the global stability of the system. The globally exponentially attractive set and positive invariant set of the system are given by constructing the positive and radial unbounded Lyapunov function. According to the concept of globally exponential tracking, for any periodic solution or equilibrium, we design some simple feedback controllers to control all the trajectories of the chaotic systems tracking the special periodic solution globally. The adaptive control method is used to realize the generalized index synchronization and the effectiveness of the method is validated by the numerical simulation.

2. Brief Introduction of the Behavior of the Dynamics and Three-Mode System of the Internal Wave Dynamics

Three-dimensional nonlinear system is achieved by the truncated spectral method in [11] as follows:

where , , and are constants, is the variational parameter, and , , are the spectral expansion coefficient.

The linear stability analysis of system (1) is carried out by [11] and the main conclusions are as follows:

(1)When , the equilibrium state is stable.(2)When , the equilibrium state is unstable and the equilibrium states and are stable, where , the real part of the equation is negative, and the pitch fork bifurcation occurs when .(3)When , the equilibrium states and lose the stability, where . But the real part of the equation is positive. The subcritical Hopf bifurcation occurs when .

For system (1), the change rate of the volume in phase space is

From formula (2), system (1) is the dissipative system of volume shrinkage on the whole, so movement will eventually tend to a certain attractor.

On the basis, in [11], the parameters can be selected as , and the initial values can be selected as , step size . The simulation calculation and analysis on the internal wave dynamics triple modular system (1) are carried out and the projection of a part of the trajectory in the phase space in the two-dimensional plane is simulated. The space orbit of the three-dimensional space is changed into the one-dimensional map by the Poincare section method. The chaotic phenomenon exists in system (1) when that is illustrated by analysis and calculation.

3. Numerical Simulation of the Dynamic Behavior of the System

More detailed numerical simulation and theoretical analysis of the dynamic behavior of system (1) are carried out in this section. The parameters can be selected as ; with the increase of the , we obtain the attractor diagram, bifurcation diagram of system (1) by simulation, and we obtain the Poincare section, power spectrum, and return mapping of system (1) by simulation, further validating the chaotic behavior of the system. The conclusions combining the simulation results with theoretical analysis are as follows:

(1)When , there is only one fixed point , and it is stable.(2)When , point becomes the unstable saddle-node point with one unstable direction and two stable directions. The eigenvalues of and all are negative real roots, so and are the asymptotically stable nodes.(3)When , the equation has a negative real root and conjugate complex roots with negative real part that represents the point in one direction of and being asymptotically stable, while there is stable focus in the plane, which is perpendicular to this direction. Numerical calculation shows that they are global attractors, as shown in Figure 1.(4)When , there is still a negative real root and others are conjugate complex roots with the positive real part that represents the point in the direction of and being stable and the focus in the plane vertical to the stable direction is unstable focus. Then Hopf bifurcation [16–18] appears when and strange attractors generated to cause chaos. Eventually, the process of chaos is shown in Figures 2, 3, 4, and 5.(5)When , Figures 6 and 7 show that the attractor gradually starts to shrink into the limit cycle which is an inverted bifurcation process, and the numerical results show that the bifurcation point satisfies the Feigenbaum constant.(6)When , chaos appears firstly and then gradually disappears to shrink into the torus which is still an inverted bifurcation process; this process is shown in Figures 8, 9, and 10, and the bifurcation point satisfies the Feigenbaum constant.(7)The bifurcation diagram of system (1) can be seen in Figure 11. From the simulation results of the image, the whole process from the starting to the end of the chaotic phenomenon can be observed. The Poincare section, power spectrum, and return mapping in Figures 12, 13, and 14 when all show the chaotic character of the system.

Figure 1: .

Figure 2: .

Figure 3: .

Figure 4: .

Figure 5: .

Figure 6: .

Figure 7: .

Figure 8: .

Figure 9: .

Figure 10: .

Figure 11: The bifurcation diagram of system (1).

Figure 12: The Poincare section of system (1).

Figure 13: The power spectrum of system (1).

Figure 14: The return mapping of system (1).

4. The Globally Exponentially Attractive Set and Positive Invariant Set

In the following we discuss the global stability of the chaotic system (1) in this section to get the estimation of globally exponentially attractive set and the positive invariant set of system (1).

First of all, we give some related definitions [19]. Denote by the state vector of system (1), and as a compact set (bounded and closed set). Let be the initial time. By we denote the solution of system (1) satisfying the initial value . If there is no confusion, we denote this solution by . Define the distance between the solution vector and the set as . Define as any set that includes the set ; that is, .

Definition 1. If there exists a compact set in the space such that as , , then one calls the set a globally attractive set of system (1). In other words, there exists , such that for . A system with the globally attractive set is called a globally asymptotically Lagrangian stable system or a dissipative system in ultimately bounded sense. If, , holds for all , the set is called a positively invariant set of system (1).

From the definition, we conclude that if is a globally attractive set then all are the globally attractive sets too.

Theorem 2. Let . For any constant , ellipsoidal

is a globally attractive set and positively invariant set of system (1).

Proof. It is easy to show that if the unique equilibrium of system (1) is globally exponentially stable. If , the equilibrium is stable but not asymptotically stable. Now we construct a family of generalized radically infinite and positive definite Lyapunov functions as

when is arbitrary constant. Since holds for , , is not positive with respect to the state vector. The point is not the equilibrium of system (1). Therefore, the family of Lyapunov functions cannot be employed to study stability of any equilibrium of system (1). So, it is called the family of generalized positive definite Lyapunov functions.
Computing the derivative of along the positive half-trajectory of system (1), we have

where , and define the function

computing the extreme value of Lagrange of function on , because is the quadratic function and the local maxima are the global maximum, and letting

we have .
Then computing the two-order derivative of , we obtain

so .
When ,  , so we have

now we prove (9) by contradiction. Otherwise, suppose that the trajectory of system (1) stays outside the set . Since is monotonic decreasing outside the set , the limit

exists. Let

where are constants. So we have

Thus, we have for . This contradiction implies that expression (9) is true.
From Definition 1, we confirm that is a globally attractive set and positively invariant set of system (1).

Remark 3. (1) Taking , we get a new cylindrate formula for system (1) as

(2) Taking , we get a new ellipsoidal formula for system (1) as

(3) Taking , we get another new ellipsoidal formula for system (1) as

(4) Taking , we get another ellipsoidal formula for system (1) as

Similarly, we can get more estimation formulas by taking different value of .
When , our estimation (14) is more accurate than estimation (15) with respect to and . In fact, we can observe the result from the following inequalities:

On the other hand, when , from the following inequality, we can see that our estimation (16) is better than estimation (14), with respect to :

When , for the variable , the estimation (16) is equal to estimation (14).
From the intersection of set theory, we construct different globally attractive sets to get effective results.

Theorem 4. Suppose , , , . The cylinder given by

is a globally attractive set and positively invariant set of system (1), where is any constant.

Proof. At first, we consider the case of . For the second and third equations of system (1), we construct a generalized positive definite and radical unbounded Lyapunov function

We compute the derivate of along the trajectory of system (1) as

for . Thus the cylinder is a globally attractive set and positively invariant set of system (1) with respect to and . Then we have .
Substituting the ultimate boundedness of into the first equation of system (1) and employing constant variation formula, we have

Therefore,

for , ; that is, is greater than and tends to exponentially.
At the same time, we have

Thus, one can see that  ,  for . That is, is less than and tends to exponentially. Consequently, the set is the globally attractive set and positively invariant set of system (1).
Consider the case of . Employing Lyapunov function (20) again, we compute the derivate of along the trajectory of system (1) as

where , and define the function

Computing the extreme value of Lagrange of function on , because is the quadratic function and the local maxima are the global maximum, and letting

we have , .
Computing the second partial derivative of ,

therefore,

When , . Now we prove by contradiction similar to Theorem 2 that , so the set is the globally attractive set and positively invariant set of system (1).
From

we have . Similar to (23) and (24), we can get , so the set is the globally attractive set with respect to that shows that is globally attractive set and positively invariant set of system (1).

5. Applications to Global and Exponential Tracking (or Stabilizing) of Any -Period Solution

There are two aims of the chaotic control. One is to stabilize an unstable (or stable but not asymptotically stable) equilibrium so that it can become a globally asymptotically stable or even a globally exponentially stable one. The other is to track a special unstable ω-period solution, that is, to control all the solutions and let them converge to this ω-period solution. Due to existence of the globally attractive set, all the solutions will enter this globally attractive compact set. In the compact set, the right side of system (1) must satisfy Lipschitz condition. This fact ensures the existence of some linear feedback controllers so that the controlled system can track some ω-period solution or can be stabilized.

Suppose that is any periodic solution or any equilibrium of system (1). Let . Then satisfies

System (31) with controllers can be written as

Now we design some simplest linear feedback controllers so that the zero solution of the controlled system (32) is globally exponentially stable. When is a periodic solution, we will confirm that all the solutions of system (1) track globally and exponentially, and when is an equilibrium , we will prove that can be globally exponentially stabilized.

Theorem 5. The control parameters always can be selected as and suitably large so that the zero solution of system (32) is globally exponentially stable. In other words, all the solutions of system (1) track the periodic solution globally and exponentially, or the equilibrium can be globally and exponentially stabilized.

Proof. (1) Suppose that is a known periodic solution. Select the parameters as

It implies that the following matrix of periodic functions is negative definite [20–23]:

Denote by the maximal eigenvalues of . Then, we have .
Construct a positive definite and radical unbounded Lyapunov function for system (32) as

We compute the derivative of it and have

According to Gronwall inequation, formula (36) is transformed as

Formula (37) implies the conclusion.
(2) Suppose that the expression of is unknown and the amplitude of it is known or can be estimated.
Suppose , , are amplitudes of the 3 components of the periodic solution. The parameters can be selected as

It implies that the following constant matrix is negative definite:

By , we denote the maximal eigenvalues of matrix ; then . Employing the Lyapunov function (35), we have

Hence, we have

This inequality implies the conclusion.
(3) Suppose the amplitude of is unknown. The periodic solution is assumed to be in the globally attractive set. We can use the boundedness of globally attractive set as the boundedness of .
If , we select and

This guarantees that the following constant matrix is negative definite:

Denote by the maximal eigenvalues of matrix . Employing the Lyapunov function (35), we have

Hence, we obtain

This ensures the conclusion.
When , we take

It ensures that the following constant matrix is negative definite:

Denote by the matrix eigenvalues of matrix . Employing the Lyapunov function (35), we have

Hence, we obtain

This ensures the conclusion. The proof of Theorem 5 is complete.

6. Globally Exponential Synchronization of Two Chaotic Systems

Definition 6. Two nonlinear dynamic systems [24] are given as follows:

where and the function is an -dimensional nonlinear function. Function is an -dimensional control input function. We call system (50) the drive system and system (51) the response system. If , then system (50) and system (51) are exponentially synchronized.

Consider two similar systems of system (1); the variables of the drive system are marked with the subscript 1 and the variables of the response system are marked with the subscript 2.

Suppose that the drive system is

The response system is

where ,  , and are the controllers designed satisfying .

Let ,  ,  ,  ; then we have the error dynamic system as follows:

The aim is to design the effective controllers ,  , and to ensure that the zero solution of the error dynamic system (54) is globally exponentially stable. It means that the drive system (52) and the response system (53) are globally exponentially synchronized. That is to say that .

We suppose that ,  ,   because the chaotic system is bounded.

Definition 7. If there exists constant satisfying for arbitrary , then the origin of the system is said to be exponentially stable.

Theorem 8. For the error dynamic system (54), the controller is designed as follows:

We can select the parameter to make the matrix

positive definite. Then the zero solution of the error dynamic system (54) is globally exponentially synchronized. It means that the drive system (52) and the response system (53) are globally exponentially synchronized.

Proof. We construct a positive definite and radical unbounded Lyapunov function as

Computing the time derivative of along the positive half-trajectory of system (1), we have

where

Obviously, the zero solution of the error dynamic system (54) is globally exponentially stable when the matrix is positive definite under the condition that the following inequalities are established:

we can have .
Thus, when , the matrix is positive definite and is negative definite. It can be observed from formula (58) and the knowledge of advanced algebra that

Therefore,

When , . It ensures that the zero solution of the error dynamic system (54) is globally exponentially stable. It means that the drive system (52) and the response system (53) are globally exponentially synchronized.

7. Numerical Simulation

In this section, we will do some numerical simulations to verify the effectiveness of the proposed methods by using the Runge-Kutta method. In the numerical simulations, the time step size is , and the parameters are chosen as , , to ensure the existence of system (1).

For the numerical simulation, we consider the systems given as (52) and (53). The initial values of the drive system and response system are taken as , . Hence, the initial error is .

Considering controller in Theorem 8, we select the feedback gain as the controller law of system (53) for the linear feedback synchronization. The simulation results are illustrated in Figures 15–17. In Figure 15, we can notice that the state variants , , of system (1) are unstable; this is when there is no controller. Figure 16 displays the time evolution of the drive system and the response system when we add the controller law. From Figure 16, we can see that the two systems starting from different conditions are synchronized with each other. The dynamics of synchronization error is defined as , whose time evolution is shown in Figure 17. It is clear that the error signals converge to zero finally.

Figure 15: State of variants , , of system (1) without control law.

Figure 16: Linear feedback synchronization .

Figure 17: Error with the change of time .

8. Conclusions

In this paper, a low mode system in the internal wave dynamics is studied. Dynamical behaviors of the system are analyzed including some basic dynamical properties, such as the attractor diagram, bifurcation diagram, Poincare section, power spectrum, and return mapping, and we obtain the estimation of the globally exponentially attractive set and the positive invariant set of the system via selecting appropriate positive definite and radical unbounded Lyapunov function. On this basis, for any periodic solution, we design some feedback controllers to control all the trajectories of the chaotic system tracking the special periodic solution globally. The globally exponential synchronization of the chaotic system is realized by the method of drive-response synchronization, and the synchronization conditions are given. The effectiveness of the method is shown by theoretical analysis and numerical simulation.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China under Project nos. 11572146 and 11526105, the Research Fund of Liaoning Provincial Education Department under Project no. L2013248, and the Jinzhou Science and Technology Special Fund under Project no. 13A1D32.

References

E. N. Lorenz, “Deterministic non-periods flows,” Atoms Science, vol. 20, pp. 130–141, 1963.

Y. Shu, Y. Zhang, and H. Xu, “Control and synchronization of general Lorenz chaotic system,” Journal of Chongqing Institute of Technology (Natural Sciences), vol. 22, no. 8, pp. 54–61, 2008.

D. G. Thomas, B. Khomami, and R. Sureshkumar, “Nonlinear dynamics of viscoelastic Taylor-Couette flow: effect of elasticity on pattern selection, molecular conformation and drag,” Journal of Fluid Mechanics, vol. 620, pp. 353–382, 2009.

W.-S. Kim, “Application of Taylor vortex to crystallization,” Journal of Chemical Engineering of Japan, vol. 47, no. 2, pp. 115–123, 2014.

H. L. Swinney and J. P. Gollub, Hydrodynamic Instabilities and the Transition to Turbulence, Springer, Berlin, Germany, 1981.

G. I. Taylor, “Stability of a viscons liquid contained between two rotating cylinders,” Philosophical Transactions of the Royal Society of London A, vol. 223, pp. 289–343, 1923.

H. Liu and B. Jiang, “Chaos control of gear system with elostomeric web based on multi-parameter multistep method,” Journal of Harbin Institute of Technology, vol. 19, no. 5, pp. 23–30, 2012.

G. R. Chen and X. Dong, From Chaos to Order: Methodologies, Perspectives and Applications, World Scientific, Singapore, 1998.

H. L. Swinney and J. P. Gollub, Hydrodynamic, Instablilities and the Transition to Turbulence, Springer, Berlin, Germany, 1981.

H. Y. Wang, “The Chaos behavior and simulation of three-model systems of Couette-Taylor flow,” Acta Mathematica Scientia. Series A, vol. 35, no. 4, pp. 769–779, 2015.

S. Liu, “The wave mechanics in chaos and atmospheric turbulence,” Science in China (Series B), vol. 5, 1986.

L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990.

L. M. Pecora and T. L. Carroll, “Driving systems with chaotic signals,” Physical Review A, vol. 44, no. 4, pp. 2374–2383, 1991.

D.-Y. Chen, C. Wu, C.-F. Liu, X.-Y. Ma, Y.-J. You, and R.-F. Zhang, “Synchronization and circuit simulation of a new double-wing chaos,” Nonlinear Dynamics, vol. 67, no. 2, pp. 1481–1504, 2012.

D. Chen, W. Zhao, X. Liu, and X. Ma, “Synchronization and antisynchronization of a class of chaotic systems with nonidentical orders and uncertain parameters,” Journal of Computational and Nonlinear Dynamics, vol. 10, no. 1, Article ID 011003, 2015.

J. Guckenheimer, Nonlinear Oscillation, Dynamical Systems, and Bifurcation of Vector Fields, Springer, Berlin, Germany, 1983.

J. H. Ge and J. Xu, “Hopf bifurcation and chaos in an inertial neuron system with coupled delay,” Science China Technological Sciences, vol. 56, no. 9, pp. 2299–2309, 2013.

H. A. El-Saka, E. Ahmed, M. I. Shehata, and A. M. A. El-Sayed, “On stability, persistence, and Hopf bifurcation in fractional order dynamical systems,” Nonlinear Dynamics, vol. 56, no. 1-2, pp. 121–126, 2009.

X. Liao, “On the new results of global attractive set and positive invariant set of the Lorenz chaotic system and the applications to chaos control and synchronization,” Science in China, Series E: Information Sciences, vol. 34, no. 12, pp. 1404–1419, 2004.

L. Huang, Theory and Application of the Stability and Robustness, National Defend Industy Press, Beijing, China, 2000.

A. A. Martynyuk, “Elements of the theory of stability of hybrid systems (review),” International Applied Mechanics, vol. 51, no. 3, pp. 243–302, 2015.

X. Liao, The Stability Theory and Application of Dynamical System, National Defend Industy Press, Beijing, China, 2000.

V. I. Vorotnikov and Yu. G. Martyshenko, “On partial stability theory of nonlinear dynamic systems,” Journal of Computer and Systems Sciences International, vol. 49, no. 5, pp. 702–709, 2010.

J. Fang, “Chaos control method, synchronization principle and its application prospect in nonlinear system (two),” Progress in Physics, vol. 16, no. 2, pp. 137–159, 1996.