Abstract

The symmetries, dynamics, and control problem of the two-dimensional (2D) Kolmogorov flow are addressed. The 2D Kolmogorov flow is known as the 2D Navier-Stokes (N-S) equations with periodic boundary conditions and with a sinusoidal external force along the -direction. First, using the Fourier Galerkin method on the original 2D Navier-Stokes equations, we obtain a seventh-order system of nonlinear ordinary differential equations (ODEs) which approximates the behavior of the Kolmogorov flow. The dynamics and symmetries of the reduced seventh-order ODE system are analyzed through computer simulations for the Reynolds number range . Extensive numerical simulations show that the obtained system is able to display the different behaviors of the Kolmogorov flow. Then, we design Lyapunov based controllers to control the dynamics of the system of ODEs to different attractors (e.g., a fixed point, a periodic orbit, or a chaotic attractor). Finally, numerical simulations are undertaken to validate the theoretical developments.

1. Introduction

In recent years, a lot of efforts have been devoted to construct dynamical systems that arise from solving the 2D Navier-Stokes equations. In the literature, the dynamics of the Navier-Stokes equations were approximated by using several reduced order models [120]. The concept of approximate inertial manifold (AIM) (see [111] and references therein) and Fourier Galerkin methods [1220] is among the methods used to obtain such reduced order systems with the task of approximating the long-time behavior of the 2D Navier-Stokes equations. The behavior of the 2D Navier-Stokes equations depends on the nature of the forcing. When the force is of a single mode, the 2D N-S equations with periodic boundary conditions are known as the 2D Kolmogorov flow.

During the last three decades, numerous numerical studies of the 2D Kolmogorov flow with different forcing terms have appeared [1129]. In 1981, Franceschini et al. [1417] constructed a system of ODEs approximating the dynamics of the Kolmogorov flow when the force acts on the mode . In [15], a number of steady states and Hopf bifurcation have been observed up to Reynolds number 100. In 1996, Armbruster et al. [25] analyzed the dynamics of the 2D Navier-Stokes equation when the force acts on the mode by exploiting as much as possible the symmetries in the problem and by using the Karhunen-Loeve decomposition (K-L). Later on and in 1997, Smaoui and Armbruster [27] described a computationally effective way to obtain a reduced order equivariant system of the 2D N-S equations using K-L decomposition with symmetries; the reduced order system obtained consists of 12 nonlinear ODEs. They showed that when Re = 16.6, the ODE system exhibits the same dynamics as the original simulation of the 2D N-S equations. However, the obtained Galerkin system does not show signs of a homoclinic behavior as observed in the original PDE simulation [25].

On the other hand, the control problem of the Navier-Stokes equations and especially to the 2D Navier-Stokes equations has not been completely investigated (see [3038] and the references therein). In 2003, using a global pinning coupling strategy, Guan et al. [34] designed an adaptive controller to control flow turbulence governed by the 2D Navier-Stokes equations. In 2009, Gambino et al. [35] designed an adaptive controller to drive the state of the system to the stationary solution. Recently, Smaoui and Zribi [3638] constructed reduced order ODE models using the truncated Fourier expansion method for approximating the dynamics of the 2D N-S equations when the force acts on the mode . They showed that, for , the dynamics of the reduced order models exhibit periodic doubling bifurcation leading to chaotic attractors. In addition, they designed Lyapunov based control laws to drive the states of the reduced order model to the basic state solution and to synchronize two reduced order ODE models having different Reynolds and starting from two different initial conditions.

In this paper, we construct a system of seven ODEs that approximates the dynamics of Kolmogorov flow when the force acts on the mode . It should be noted that the approach used in this paper to construct the system of ODEs is different than the one used in [3638]. In addition, unlike the system obtained earlier by Smaoui and Armbruster [27], this system shows a homoclinic gluing bifurcation similar to the one observed in the original PDE simulations [25]. Furthermore, Lyapunov based control laws are designed to control the dynamics of the system for a given Re.

The paper is organized as follows. In Section 2, the 2D Navier-Stokes equations are presented and a seventh-order nonlinear ODE system is obtained to approximate the behavior of these equations. The reduced order ODE system is also analyzed in Section 2. Section 3 presents the design of a control law which is used to regulate the states of the reduced order ODE system to a desired fixed point. Section 4 deals with the design of a control scheme to synchronize two reduced order systems obtained from the 2D N-S equations having the same or different Reynolds number but they start from different initial conditions. Finally, some concluding remarks are given in Section 5.

2. The 2D Kolmogorov Flow

2.1. The Seven-Mode Reduced Order System of the 2D Kolmogorov Flow

The “basic 2D Kolmogorov flow” was introduced by Kolmogorov in 1958 as an example on which to study transition to turbulence [39]. This basic flow is the solution of the 2D Navier-Stokes equations with periodic boundary conditions in two directions given bywith and force , which is assumed to be stationary and spatially biperiodic. The kinematic viscosity is , where is the Reynolds number and the pressure is .

The perturbed nondimensional vorticity formulation of the Kolmogorov flow is where the scaled time , , , and .

In Smaoui and Zribi [3638], the derivation of the reduced order system of ODEs was based on expanding the stream function in (2) as follows:

In this paper, we use a totally different approach than the one used in [3638]. We derive a system of seven ODEs from the Navier-Stokes equations by expanding in the following form:where is a wave vector with integer components, , and the reality condition must hold.

The equation for iswhere and is the component of with respect to

Next, if the basic flow , then one can consider the set of vectors, , , , , , , and , and their negatives in (5) to obtain the following system of ODEs:

Let , , , , , , and ; then system (6) becomes

After changing the length scale for , the time scale and , and considering the forcing term to act on mode , system (7) becomes

The system given by the equations in (8) can be written in the following form: where the the vector is such that and the diagonal matrix is such that:and the nonlinear vector is such thatwith

Remark 1. It is noted that system (8) is invariant under the following symmetries:where ,  , and are reflection symmetries across the -axis, the -axis, and the origin, respectively. Hence, , , with the identity transformation form an Abelian group: .

2.2. The Dynamics of the Seven-Mode Truncation System

In this subsection, we analyze the dynamics of the seven-mode truncation ODE system presented in (8) for different Reynolds numbers. The DsTool software [40] is used in all numerical simulations presented in this section using the 4th-order Runge-Kutta method as the numerical integrator and with the time step .(i)For , the basic fixed point is the only stable solution and, numerically, it is globally attractive for all . This is a special case of the general results on the theory of Navier-Stokes equations [41].(ii)For there are three fixed points: the old one , which becomes unstable as a result of the crossing of the imaginary axes by one of the eigenvalues of the Lyapunov matrix, and two more, and , that bifurcated from are stable and attracting. Numerical evidence shows that any randomly chosen initial condition is either attracted by or .(iii)For , there are seven fixed points: four asymptotically stable and three unstable. Figure 1 presents the phase portrait of the four asymptotically stable fixed points at . The first asymptotically state fixed point (Figure 1, top left) was obtained using the following initial conditions:The second asymptotically stable fixed point (Figure 1, bottom right) was obtained by applying the reflection symmetry on the initial conditions used for the first asymptotically stable fixed point. That is,The third asymptotically stable fixed point (Figure 1, top right) was generated by applying the reflection symmetry on the initial conditions used for the second asymptotically stable fixed point. That is,The fourth asymptotically stable fixed point (Figure 1, bottom left) was obtained by applying the reflection symmetry on the initial conditions used for the second asymptotically stable fixed point. That is,(iv)For the four asymptotically stable fixed points become unstable because a pair of complex eigenvalues crosses the imaginary axes, and bifurcate via a Hopf bifurcation into a four stable periodic orbits around these four fixed points. The four stable periodic orbits remain stable up to . The other three remain unstable (see Figure 2).(v)For . At each of the four periodic orbits loses stability and bifurcates into a new periodic orbit with double the period.(vi)For , a homoclinic gluing bifurcation similar to the one observed by Armbruster et al. [25] occurs (see Figure 3). It should be noted that the homoclinic gluing bifurcation observed leads to two chaotic or strange attractors each connecting a pair of two stable periodic orbits. Figure 3(b) shows that the dynamics of one of the chaotic attractors follow a quasiperiodic regime; then either bursts to the same attractor or to its symmetric counterpart; then other bursts follow. Numerical simulations indicate that intervals between bursts become shorter as the Reynolds number, , is increased up to .

3. Controlling the Dynamics of the Reduced Order System to a Fixed Point

This section deals with the design of a control scheme to drive the states of the system to a stable or unstable desired fixed point. Hence, control inputs are added to the system in (8). The controlled system of ODEs with control inputs added to the first, the third, the fifth, and the sixth ODE of (8) is as follows.

Let the constant desired fixed point be such that . Since is an equilibrium point of system (18) then it must satisfy the following algebraic equations:Define the errors such that

Using equations (18), (19), and (20), the model for the error system can be written as follows:where the parameters , , , and are such that

Let the gains be positive scalars and let the control gains be positive scalars such that

Theorem 2. The control law,when applied to the error model (21) guarantees the convergence of the errors to zero as tends to infinity. Therefore, system (18) asymptotically converges to a desired fixed point as tends to infinity.

Proof. Consider the Lyapunov function candidate such thatUsing the model of the errors system given in (21), the control law given by (24), and the constraints on the control parameters and given by (23), the derivative of with respect to time is such thatSince the design parameters and are positive constants, then it can be concluded that is negative definite. Therefore, the Lyapunov function in (25) is positive definite, its derivative with respect to time is negative definite, and is radially unbounded. Hence, it can be concluded that the errors in (20) asymptotically converge to zero as tends to infinity. Thus, the states of system (18) asymptotically converge to a desired constant fixed point as tends to infinity.

The simulation results of system (18) when controlled using the proposed control scheme given by (24) are presented below. The simulations are carried out for the Reynolds number . Note that the dynamics at corresponds to four stable periodic solutions of the system of ODEs. These periodic solutions are a result of Hopf bifurcations that took place at . The values of the control gains are such that , , , and . The control gains are taken to be , , , , , , and . The simulation results correspond to an initial condition,and the four desired fixed pointsThese desired fixed points are the asymptotically stable fixed points that became unstable as the result of Hopf bifurcations. At the beginning of the simulations, the controllers are set to zero for the first 200,000. Then the control law given by (24) is switched on and off and the desired fixed point is set to one of the four desired fixed points each time the control law is switched on.

The simulation results are presented in Figures 4 and 5. Figure 4 presents the controlled states versus , and Figure 5 depicts the time evolution of the controlled states . The figures clearly show how the controlled states converge to one of the desired fixed points each time the control law is switched on and to one of the four stable periodic solutions when the control law is switched off. Hence, it can be concluded that the proposed controller given by (24) drives the states of the system in (18) to their desired states.

4. Synchronization of Two ODE Systems Having the Same or Different Reynolds Numbers

In this section, we design a control law to synchronize two ODE systems obtained from the 2D Navier-Stokes equation using the Fourier Expansion truncation method having the same or different Reynolds numbers.

4.1. Model of the Error System

In this subsection, we first derive the error system between the two ODE systems to be synchronized. The model of the first ODE system iswhere is the Reynolds number of the first ODE system. The model of the second ODEs system iswhere are the state variables of the second ODE and is its Reynolds number. Define the errors between the first ODE system in (29) and the second ODE system in (30) as follows:

Using the equations defined in (29)–(31), the model of the error system can be written as follows:where .

4.2. A Lyapunov Based Controller to Synchronize Two Systems Having the Same or Different Reynolds Number

In this subsection, a Lyapunov based controller is designed to drive the states of the system in (30) to asymptotically converge to the states of the system in (29).

Let the gains be positive scalar. In addition, let the control gains be positive scalars such that

The following theorem gives the second control result of the paper.

Theorem 3. The control law,when applied to the error model (32) guarantees the convergence of the errors to zero as tends to infinity. Therefore, the states of system (30) asymptotically converge to the states of system (29) as tends to infinity.

Proof. Consider the Lyapunov function candidate such thatUsing the model of the errors system given in (32), the control law given by (34), and the constraints on the control parameters and given by (33), the derivative of with respect to time is such thatorSince the design parameters and are positive constants, then it is concluded that the Lyapunov function defined by (35) is positive definite and its derivative is negative definite. Therefore, the error functions in (31) asymptotically converge to zero as tends to infinity. Thus, the states of system (30) asymptotically converge to the states of system (29) as tends to infinity.

The performance of the controller in (34) when applied to the ODE system in (18) is simulated. The control gains , and are taken to be 50. Three cases are simulated. For the first case, we choose , and the initial condition    =              for the system in (29) and and the initial condition = for the system in (30). Note that these two states correspond to two stable periodic orbits with different periods. For the second case, we choose , and for the system in (29), and and for the system in (30). The first state in the second case corresponds to a periodic orbit, and the second state corresponds to a chaotic attractor or a strange attractor connecting a pair of two stable periodic orbits. For the third case, we choose and =        for the system in (29), and and for the system in (30). Note that, in this case, the two states correspond to two different chaotic or strange attractors each connecting a pair of two stable periodic orbits. Moreover, at the beginning of simulations in each case, the controllers , and are set to zero for the first 30 seconds. Then the control law presented in (34) is switched on to force system (30) to synchronize with system (29).

The simulation results are presented in Figure 6 for case 1. Figure 6(a) depicts the norm of the error versus time. Also, the states and versus time and and versus time are plotted in Figures 6(b) and 6(c), respectively. Figure 6(d) plots the states and versus and ; the figure shows the effectiveness of the controller to drive the dynamics from one attractor to a different attractor. Figure 7 shows the simulation results for case 2, and Figure 8 shows the results for case 3. In the three cases, it is shown how the error converges to zero. Hence, it can be concluded that the designed control law in (34) is able to synchronize the ODE systems in (29) and (30) when these systems have the same or different Reynolds numbers but they start from two different initial conditions.

5. Conclusion

This paper investigates the symmetries, dynamics, and control of the 2D Kolmogorov flow with a forcing in the mode . First, a seventh-order system of ODEs is constructed using the Fourier series expansion. A computer analysis shows that the constructed system exhibits a chaotic behavior for a certain range of the Reynolds number. In addition, the ODE system shows a homoclinic gluing bifurcation similar to the original PDE simulation. Then, Lyapunov based controllers are designed to control the dynamics of ODEs to a fixed point, a periodic orbit, or a chaotic or strange attractor. Extensive numerical studies are presented to show the effectiveness of the proposed controllers.

Future work will address the design of adaptive control schemes for the reduced order ODE model of the 2D Kolmogorov flow.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported and funded by the Research Sector, Kuwait University, under Grant no. SM07/16.