Research Article | Open Access

Ismail Kucuk, "Active Optimal Control of the KdV Equation Using the Variational Iteration Method", *Mathematical Problems in Engineering*, vol. 2010, Article ID 929103, 10 pages, 2010. https://doi.org/10.1155/2010/929103

# Active Optimal Control of the KdV Equation Using the Variational Iteration Method

**Academic Editor:**Jihuan He

#### Abstract

The optimal pointwise control of the KdV equation is investigated with an objective of minimizing a given performance measure. The performance measure is specified as a quadratic functional of the final state and velocity functions along with the energy due to open- and closed-loop controls. The minimization of the performance measure over the controls is subjected to the KdV equation with periodic boundary conditions and appropriate initial condition. In contrast to standard optimal control or variational methods, a direct control parameterization is used in this study which presents a distinct approach toward the solution of optimal control problems. The method is based on finite terms of Fourier series approximation of each time control variable with unknown Fourier coefficients and frequencies. He's variational iteration method for the nonlinear partial differential equations is applied to the problem and thus converting the optimal control of lumped parameter systems into a mathematical programming. A numerical simulation is provided to exemplify the proposed method.

#### 1. Introduction

A modal for planar, unidirectional waves propagating in shallow water was originally introduced by Korteweg and de Vries in 1895 [1]. The modal is expressed by a third-order nonlinear partial differential equation called KdV equation. The KdV equation has been at the center of naval science studies and other physical phenomena such as weakly nonlinear long waves for the last 150 years. Therefore, solving and controlling the behavior of the KdV equation have great implications.

Review of new techniques such as variational approaches, parameter-expanding methods, and parameterized perturbation method for nonlinear problems is presented by He in [2], and a detailed study of He's approaches is given in [3]. In the literature, there are a considerable number of numerical and theoretical aspects of the KdV equation. A survey of results for the KdV equation is given in [4]. Existence and uniqueness of the solution of the KdV equation are given in different forms [5]. Numerical solution of the KdV equation is studied by using modified Bernstein polynomials in [6], Galerkin B-spline finite element method in [7], homotopy perturbation method to find solitary-wave solutions of the combined Korteweg de Vries-Modified Korteweg de Vries Equation in [8], and variational iteration method (VIM) in [9, 10]. Although, the numerical solution of the KdV equation has been studied in depth, to the best knowledge of the author, optimal control aspect of the problem did not attract many researchers. Control applications of nonlinear dispersive wave equations in [11], exact controllability and stability of the KdV equation in [12], theoretical aspect of boundary controllability of KdV on a bounded domain in [13], and stability and numerical aspect of the boundary control of the KdV equation in [14] are discussed.

In this paper, we consider an active control of nonlinear waves expressed by the KdV equation with periodic boundary conditions and initial conditions. To control water waves in a uniform channel, point-wise control actuators in the spatial domain, and linear displacement and displacement-slope feedback controls are implemented in the KdV equation. The dynamic response of the system is measured by performance index functionals that consist of weighted sum of the energy at the terminal time with the total effort of open- and closed-loop controls. The objective of the control problem is to minimize the dynamic response of the system with minimum expenditure of the modified energy. The parameterization of the actuators uses a finite term of orthogonal (or nonorthogonal) functions with unknown coefficients and the solution of state function is expressed as an iterative function with a Lagrange multiplier known as VIM. Thereby, the optimal control problem becomes a mathematical programming for unknown coefficients to be computed optimally while state solution is obtained iteratively. The compact solution of state function is expressed analytically in terms of unknown terms due to applied controls. To compare the effects of different controls, first the open-loop control is applied to the system before both open- and closed-loop controls are applied, or closed-loop control with an optimal actuator is applied.

The computational and graphical results show that the present method has a desired robustness. Moreover, it is observed that closed-loop control with an optimal actuator applied to the system reduces the energy substantially and controls the behavior of the elongation and velocity of waves.

#### 2. Problem Formulation

We consider the KdV equation subject to the following periodic boundary conditions and initial condition: where describes the elongation of the wave at a point , and ; terminal time is ; is the dirac delta function; and are constants that are determined by the nature of a physical application; and are feedback gains; is the number of control actuators placed in the interior of the spatial domain, and is the applied force for that will be obtained optimally.

*Remarks 2.1. *Let the set of Lebesue integrable functions be given as
with the inner product and the norm being
and a Hilbert space is with the inner product.

Our goal in this paper is to reduce the dynamical response of nonlinear waves modelled by the KdV equation by implementing open- and closed-loop controls. The natural attempt is to minimize the energy due to waves that should be achieved with limited expenditure of control energy. The weighted performance of the system based on energy at the terminal time and total control effort is measured by the following performance functional: where in which , for such that and are weighting factors. In (2.5), is the energy at the terminal time, and and are energies for the open-loop and closed-loop controls duration over .

Here we are considering three optimal control problems.(P1)The first is to find optimal for fixed real valued and such that where is the optimum open-loop performance, and subscript indicates open-loop optimization.(P2)Secondly, we seek optimal feedback gains, , for optimal actuators obtained in (P1) such that where is the optimum closed-loop performance, and subscript indicates closed-loop optimization.(P3)Finally, we seek optimal feedback gains, , and optimal actuators, , such that

In problems (2.7), (2.8) and (2.9), is subject to (2.1) and (2.2).

#### 3. Solving Optimal Control Problems

The distributed parameter control problems in (2.7)–(2.9) are transformed into an iterative parameter control problem in which a parameterization of actuators is introduced as a direct method. The transformation is done by VIM. Optimal control problems are then transformed into a mathematical programming problem that consists of unknowns due to parameterization of actuators and the feedback parameters at each search of optimal values the solution of the KdV is obtained iteratively. In the following subsections, the parametrization of actuators is first introduced. It is followed by the solution technique for KdV with a brief introduction to VIM.

##### 3.1. Control Parameterization

A direct control parameterization is introduced using finite terms of a Fourier series approximation. In the approximation, each actuator is expressed as a finite sum of a Fourier series with unknown Fourier coefficients and frequencies. Actuators in (2.1) are given by where and It follows from (3.1) that needs to be calculated optimally in (P1) and (P3).

##### 3.2. Application of He's Variational Iteration Method

It is observed that the VIM converges to the exact solutions, and VIM is used for many nonlinear partial differential equations successfully such as [9, 10, 15]. Although the concern about the convergence rate depends on the accurate calculation of the Lagrange multiplier [16], the nonlinear terms in nonlinear problems are taken as restricted variations in order to determine the Lagrange multipliers. A brief introduction to VIM is given in this section, but readers are referred to a recent review of the method by He in [16, 17] and references therein.

For a nonlinear partial differential equation of the following form where is a linear operator, is a nonlinear operator, and is a nonhomogeneous term, order approximation to is obtained iteratively by a correction function in -direction in which is a general Lagrange multiplier, and is a restricted variation, that is, The general Lagrange multiplier is found optimally via the variational theory [18].

For (2.1), a correctional functional is The dirac delta function in (2.1) is taken as a pseudo dirac delta function for the sake of simplicity in VIM.

If the variation is taken with respect to in (3.5), the following stationary equations are obtained for : from which it follows immediately that the Lagrange multiplier, , is identified as Thus, optimal control problem is reduced to a variational iteration with lumped parameters: where is given by the parametrization (3.1). The initial guess, , is given by a trial function that satisfies boundary conditions.

The closed-form expression for in (3.8) is found with the aid of . The resulting is substituted into (2.7)–(2.9) for the prospective optimization problem. To proceed with problems (P1)–(P3), necessary adjustments are done in (3.8). For (P1), given by (3.2) is calculated optimally for fixed and for which unknown terms have to be determined in the solution in (3.8). For (P3), the solution in (3.8) consists of unknown due to parameterization and feedback gains that are calculated optimally. The iterative terms of are obtained through , and then the terms are placed in appropriate performance functional. Finally, optimization toolbox in is used to calculate the optimal parameters. In calculations, a zero vector for as an initial guess is taken until a convergence is reached in sense, that is, The same procedure is repeated for the calculations of feedback gains. Therefore, the algorithm given in [19] is performed.

#### 4. Numerics

In this section, the proposed technique is examined numerically. In the calculations, the following data is used: , , , , , and a nascent delta function in (2.1). We will conduct the simulations for VIM and Adomian methods.

*Case 1 ( and ). *For the given data above, the nonlinear partial differential equation is solved by using VIM where only one iteration is used to ease the complexity in the calculations without compromising the convergence in the solution. The solution for one iteration is found by using and the performance functions in (2.7)–(2.9) are written in to perform command to find unknown parameters optimally. An optimal result is obtained after three iterations with the updated initial guess obtained at the end of each iteration in the Algorithm [19].

In the problem (P1), optimal actuator (3.2) is obtained as In the problem (P2), optimal feedback gains are found as and for the optimal given by (4.1). For problem (P3), an optimal parameters , , and are found.

The energy for the uncontrolled system is 0.2459 that is reduced by 0.5% when the open-loop control is applied to the system alone and by 49% when the open- and closed-loop controls are applied to the system simultaneously. 78% is the reduction in the energy when the closed-loop is applied to the system along with the optimal actuator obtained in the open-loop control. The elongation of the wave, at a fixed point and at the terminal time are presented for uncontrolled cases in Figures 1 and 2, respectively. The velocity of elongation of the wave, at a fixed point and at the terminal time is presented for uncontrolled cases in Figures 3 and 4, respectively.

#### 5. Discussions and Conclusion

The active control of waves defined by the KdV equation is studied by implementing point-wise actuators in the spatial domain and linear displacement and slope-displacement feedback controls in the system. An energy-based performance measure for control is minimized to get a minimum expenditure in the total control effort. It is observed that open-loop control alone does not reduces the energy, and simultaneously applying the open- closed-loop controls reduce the energy noticeably. On the contrary, closed-loop control with optimal actuators reduces the energy substantially. In Figures 1 and 3, it is observed that when one actuator is applied to the system, the displacement and the velocity show indifferent behaviors with the uncontrolled cases, respectively, that might lead to different results if more than one actuator is placed in the spatial domain. The behavior of the system over time is presented in Figure 2 in which it becomes clear that the open-loop control for one actuator in the spatial domain does not improve the system. The nonlinear involvement of the optimal control problem is solved by using VIM and available toolboxes in . Although applying the VIM is of great achievement in the optimal control problems for a nonlinear phenomena, the computational cost in the calculations increases with the number of iterations taken in VIM. It should also be noted that accurate calculation of the Lagrange multiplier in VIM reduces the number of iterations in the calculations.

#### References

- D. J. Korteweg and F. de Vries, “On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves,”
*Philosophical Magazine*, vol. 39, pp. 422–443, 1895. View at: Google Scholar - J.-H. He, “Some asymptotic methods for strongly nonlinear equations,”
*International Journal of Modern Physics B*, vol. 20, no. 10, pp. 1141–1199, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Some relatively new techniques for nonlinear problems,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 234849, 25 pages, 2009. View at: Google Scholar | Zentralblatt MATH | MathSciNet - R. M. Miura, “The Korteweg-de Vries equation: a survey of results,”
*SIAM Review*, vol. 18, no. 3, pp. 412–459, 1976. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Z. Yin, “On the global weak solutions to an integrable shallow water equation with linear and nonlinear dispersion,”
*Dynamics of Continuous, Discrete & Impulsive Systems. Series A*, vol. 13, no. 3-4, pp. 321–336, 2006. View at: Google Scholar | Zentralblatt MATH | MathSciNet - D. D. Bhatta and M. I. Bhatti, “Numerical solution of KdV equation using modified Bernstein polynomials,”
*Applied Mathematics and Computation*, vol. 174, no. 2, pp. 1255–1268, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - E. N. Aksan and A. Özdeş, “Numerical solution of Korteweg-de Vries equation by Galerkin B-spline finite element method,”
*Applied Mathematics and Computation*, vol. 175, no. 2, pp. 1256–1265, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - H. Mirgolbabaei and D. D. Ganji, “Application of homotopy perturbation method to solve combined Korteweg de Vries-Modified Korteweg de Vries equation,”
*Journal of Applied Sciences*, vol. 9, no. 19, pp. 3587–3592, 2009. View at: Publisher Site | Google Scholar - M. A. Abdou and A. A. Soliman, “New applications of variational iteration method,”
*Physica D*, vol. 211, no. 1-2, pp. 1–8, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Inc, “Numerical simulation of KdV and mKdV equations with initial conditions by the variational iteration method,”
*Chaos, Solitons and Fractals*, vol. 34, no. 4, pp. 1075–1081, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. Zhang,
*Some results for nonlinear dispersive wave equations with applications to control*, Ph.D. thesis, University of Wisconsin-Madison, 1990. - D. L. Russell and B. Y. Zhang, “Exact controllability and stabilizability of the Korteweg-de Vries equation,”
*Transactions of the American Mathematical Society*, vol. 348, no. 9, pp. 3643–3672, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Rosier, “Exact boundary controllability for the Korteweg-de Vries equation on a bounded domain,”
*Control, Optimisation and Calculus of Variations*, vol. 2, pp. 33–55, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. Balogh and M. Krstic, “Boundary control of the Korteweg-de Vries-Burgers equation: further results on stabilization and numerical demonstration,” in
*Proceedings of the 7th Mediterranean Conference on Control and Automation*, pp. 1341–1349, Haifa, Israel, 1999. View at: Google Scholar - M. A. Abdou and A. A. Soliman, “Variational iteration method for solving Burger's and coupled Burger's equations,”
*Journal of Computational and Applied Mathematics*, vol. 181, no. 2, pp. 245–251, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, “Variational iteration method—some recent results and new interpretations,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 3–17, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-H. He, G.-C. Wu, and F. Austin, “The variational iteration method which should be followed,”
*Nonlinear Science Letters A*, vol. 1, no. 1, pp. 1–30, 2010. View at: Google Scholar - J. He, “A new approach to nonlinear partial differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 2, no. 4, pp. 230–235, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH - I. Kucuk and I. Sadek, “An efficient computational method for the optimal control problem for the Burgers equation,”
*Mathematical and Computer Modelling*, vol. 44, no. 11-12, pp. 973–982, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2010 Ismail Kucuk. 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.