- Research
- Open access
- Published:
On the solutions and conservation laws of the -dimensional higher-order Broer-Kaup system
Boundary Value Problems volume 2013, Article number: 41 (2013)
Abstract
In this paper we obtain exact solutions of the -dimensional higher-order Broer-Kaup system which was obtained from the Kadomtsev-Petviashvili equation by the symmetry constraints. The methods used to determine the exact solutions of the underlying system are the Lie group analysis and the simplest equation method. The solutions obtained are the solitary wave solutions. Moreover, we derive the conservation laws of the -dimensional higher-order Broer-Kaup system by employing the multiplier approach and the new conservation theorem.
1 Introduction
In this paper we study the -dimensional higher-order Broer-Kaup system
which was first introduced by Lou and Hu [1] by considering the symmetry constraints of the Kadomtsev-Petviashvili equation. The system (1.1a) and (1.1b) is in fact an extension of the well-known -dimensional Broer-Kaup system [2–4]
which is used to model the bi-directional propagation of long waves in shallow water. In [5], Fan derived a unified Darboux transformation for the system (1.1a) and (1.1b) with the help of a gauge transformation of the spectral problem and as an application obtained some new explicit soliton-like solutions. Recently, Huang et al. [6] presented a new N-fold Darboux transformations of the -dimensional higher-order Broer-Kaup system with the help of a gauge transformation of the spectral problem and found new explicit multi-soliton solutions of the system (1.1a) and (1.1b).
In the latter half of the nineteenth century, Sophus Lie (1842-1899) developed one of the most powerful methods to determine solutions of differential equations. This method, known as the Lie group analysis method, systematically unifies and extends well-known ad hoc techniques to construct explicit solutions of differential equations. It has proved to be a versatile tool for solving nonlinear problems described by the differential equations arising in mathematics, physics and in other scientific fields of study. For the theory and application of the Lie group analysis methods, see, e.g., the Refs. [7–12].
Conservation laws play a vital role in the solution process of differential equations. Finding conservation laws of the system of differential equations is often the first step towards finding the solution [7]. Also, the conservation laws are useful in the numerical integration of partial differential equations [13], for example, to control numerical errors. The determination of conservation laws of the Korteweg de Vries equation, in fact, initiated the discovery of a number of methods to solve evolutionary equations [14]. Moreover, conservation laws play an important role in the theories of non-classical transformations [15, 16], normal forms and asymptotic integrability [17]. Recently, in [18] the conserved quantity was used to determine the unknown exponent in the similarity solution which cannot be obtained from the homogeneous boundary conditions.
In this paper, we use the Lie group analysis approach along with the simplest equation method to obtain exact solutions of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b). Furthermore, conservation laws will be computed for (1.1a) and (1.1b) using the two approaches: the new conservation theorem due to Ibragimov [19] and the multiplier method [20, 21].
2 Symmetry reductions and exact solutions of (1.1a) and (1.1b)
The symmetry group of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) will be generated by the vector field of the form
Applying the third prolongation [11] to (1.1a) and (1.1b) and solving the resultant overdetermined system of linear partial differential equations one obtains the following three Lie point symmetries:
2.1 One-dimensional optimal system of subalgebras
In this subsection we present an optimal system of one-dimensional subalgebras for the system (1.1a) and (1.1b) to obtain an optimal system of group-invariant solutions. The method which we use here for obtaining a one-dimensional optimal system of subalgebras is that given in [11]. The adjoint transformations are given by
Here is the commutator given by
The commutator table of the Lie point symmetries of the system (1.1a) and (1.1b) and the adjoint representations of the symmetry group of (1.1a) and (1.1b) on its Lie algebra are given in Table 1 and Table 2, respectively. Table 1 and Table 2 are used to construct an optimal system of one-dimensional subalgebras for the system (1.1a) and (1.1b).
From Tables 1 and 2 one can obtain an optimal system of one-dimensional subalgebras given by .
2.2 Symmetry reductions of (1.1a) and (1.1b)
In this subsection we use the optimal system of one-dimensional subalgebras calculated above to obtain symmetry reductions that transform (1.1a) and (1.1b) into a system of ordinary differential equations (ODEs). Later, in the next subsection, we will look for exact solutions of (1.1a) and (1.1b).
Case 1.
The symmetry gives rise to the group-invariant solution
where is an invariant of the symmetry . Substitution of (2.1) into (1.1a) and (1.1b) results in the system of ODEs
Case 2.
The symmetry gives rise to the group-invariant solution of the form
where is an invariant of and the functions F and G satisfy the following system of ODEs:
Case 3.
By solving the corresponding Lagrange system for the symmetry , one obtains an invariant and the group-invariant solution of the form
where the functions F and G satisfy the following system of ODEs:
2.3 Exact solutions using the simplest equation method
In this subsection we use the simplest equation method, which was introduced by Kudryashov [22, 23] and modified by Vitanov [24] (see also [25]), to solve the ODE system (2.2a) and (2.2b), and as a result we will obtain the exact solutions of our -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b). Bernoulli and Riccati equations will be used as the simplest equations.
Let us consider the solutions of the ODE system (2.2a) and (2.2b) in the form
where satisfies the Bernoulli and Riccati equations, M is a positive integer that can be determined by balancing procedure as in [24] and , are parameters to be determined. It is well known that the Bernoulli and Riccati equations are nonlinear ODEs whose solutions can be written in terms of elementary functions.
We consider the Bernoulli equation
where a and b are constants. Its solution is given by
where C is a constant of integration.
For the Riccati equation
where a, b and c are constants, we will use the solutions
and
where and C is a constant of integration.
2.3.1 Solutions of (1.1a) and (1.1b) using the Bernoulli equation as the simplest equation
The balancing procedure [24] yields , so the solutions of (2.2a) and (2.2b) are of the form
Substituting (2.8) into (2.2a) and (2.2b) and making use of (2.6) and then equating all coefficients of the functions to zero, we obtain an algebraic system of equations in terms of , , and , , . Solving the system of algebraic equations with the aid of Mathematica, we obtain the following cases.
Case 1
Thus, a solution of our -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) is
where and C is a constant of integration.
Case 2
and so a solution of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) is
where and C is a constant of integration.
2.3.2 Solutions of (1.1a) and (1.1b) using Riccati equation as the simplest equation
The balancing procedure yields , so the solutions of the ODE system (2.2a) and (2.2b) are of the form
Substituting (2.11) into (2.2a) and (2.2b) and making use of (2.7), we obtain an algebraic system of equations in terms of , , , , , by equating all coefficients of the functions to zero. Solving the algebraic equations, one obtains the following cases.
Case 1
and hence the solutions of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) are
and
where and C is a constant of integration.
Case 2
In this case the solutions of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) are given by
and
where and C is a constant of integration.
Case 3
The solutions in this case are
and
where and C is a constant of integration.
Case 4
and so the solutions are
and
where and C is a constant of integration.
Case 5
The solutions are
and
where and C is a constant of integration.
A profile of the solution (2.21a) and (2.21b) is given in Figure 1.
3 Conservation laws of (1.1a) and (1.1b)
In this section, we derive conservation laws for the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b). Two different approaches will be used. Firstly, we use the new conservation method due to Ibragimov [19] and then employ the multiplier method [20, 21]. We now present some preliminaries that we will need later in this section.
3.1 Preliminaries
In this subsection we briefly present the notation and pertinent results which we utilize below. For details the reader is referred to [8–10, 19–21, 27].
3.1.1 Fundamental operators and their relationship
Consider a k th-order system of PDEs of n independent variables and m dependent variables
where denote the collections of all first, second, …, k th-order partial derivatives, that is, , , …, respectively, with the total derivative operator with respect to given by
where the summation convention is used whenever appropriate.
The Euler-Lagrange operator, for each α, is given by [8–10]
and the Lie-Bäcklund operator is
where is the space of differential functions. The operator (3.4) is an abbreviated form of the infinite formal sum
where the additional coefficients are determined uniquely by the prolongation formulae
in which is the Lie characteristic function given by
The Lie-Bäcklund operator (3.5) can be written in a characteristic form as
The Noether operators associated with the Lie-Bäcklund symmetry operator X are given by
where the Euler-Lagrange operators with respect to derivatives of are obtained from (3.3) by replacing by the corresponding derivatives. For example,
and the Euler-Lagrange , Lie-Bäcklund and Noether operators are connected by the operator identity
The n-tuple vector , , , is a conserved vector of (3.1) if satisfies
The equation (3.12) defines a local conservation law of the system (3.1).
3.1.2 Multiplier method
A multiplier has the property that
hold identically. We consider multipliers of the third-order, that is,
The right-hand side of (3.13) is a divergence expression. The determining equation for the multiplier is
The conserved vectors are calculated via a homotopy formula [20, 21, 26] once the multipliers are obtained.
3.1.3 Variational method for a system and its adjoint
A system of adjoint equations for the system of k th-order differential equations (3.1) is defined by [27]
where
and are new dependent variables.
The following results are given in Ibragimov [19] and recalled here.
Assume that the system of equations (3.1) admits the symmetry generator
Then the system of adjoint equations (3.15) admits the operator
where the operator (3.18) is an extension of (3.17) to the variable and the are obtainable from
Theorem 1 [19]
Every Lie point, Lie-Bäcklund and nonlocal symmetry (3.17) admitted by the system of equations (3.1) gives rise to a conservation law for the system consisting of equation (3.1) and adjoint equation (3.15), where the components of the conserved vector are determined by
with Lagrangian given by
3.2 Construction of conservation laws for (1.1a) and (1.1b)
We now construct conservation laws for the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b) using the two approaches.
3.2.1 Application of the multiplier method
For the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b), after some lengthy calculations, we obtain the third-order multipliers
and
that are given by
where , are arbitrary constants. Corresponding to the above multipliers, we obtain the following seven local conserved vectors of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b):
Remark 1 Higher-order conservation laws of (1.1a) and (1.1b) can be computed by increasing the order of multipliers.
3.2.2 Application of the new conservation theorem
The adjoint equations of (1.1a) and (1.1b), by invoking (3.16), are given by
where and are the new dependent variables. By recalling (3.21), we get the following Lagrangian for the system of equations (1.1a) and (1.1b) and (3.24a) and (3.24b):
Because of the three Lie point symmetries of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b), we have the following three cases to consider:
-
(i)
We first consider the Lie point symmetry generator of the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b). Corresponding to this symmetry, the Lie characteristic function is . Thus, by using (3.20), the components of the conserved vector are given by
-
(ii)
The Lie point symmetry generator has the Lie characteristic function . Hence using (3.20), one can obtain the conserved vector whose components are
-
(iii)
Finally, we consider the symmetry generator . For this case, the Lie characteristic function , and by using (3.20), the components of the conserved vector are given by
Remark 2 The components of the conserved vectors contain the arbitrary solutions ϕ and ψ of adjoint equations (3.24a) and (3.24b), and hence one can obtain an infinite number of conservation laws.
4 Concluding remarks
In this paper we have studied the -dimensional higher-order Broer-Kaup system (1.1a) and (1.1b). Similarity reductions and exact solutions, with the aid of the simplest equation method, were obtained based on optimal systems of one-dimensional subalgebras for the underlying system. We have verified the correctness of the solutions obtained here by substituting them back into the system (1.1a) and (1.1b). Furthermore, conservation laws for the system (1.1a) and (1.1b) were derived by using the multiplier method and the new conservation theorem.
References
Lou SY, Hu XB: Infinitely many Lax pairs and symmetry constraints of the KP equation. J. Math. Phys. 1997, 38: 6401-6427. 10.1063/1.532219
Broer LJF: Approximate equations for long water waves. Appl. Sci. Res. 1975, 31: 377-395. 10.1007/BF00418048
Kaup DJ: A higher-order water wave equation and the method for solving it. Prog. Theor. Phys. 1975, 54: 396-408. 10.1143/PTP.54.396
Kaup DJ: Finding eigenvalue problems for solving nonlinear evolution equations. Prog. Theor. Phys. 1975, 54: 72-78. 10.1143/PTP.54.72
Fan EG: Solving Kadomtsev-Petviashvili equation via a new decomposition and Darboux transformation. Commun. Theor. Phys. 2002, 37: 145-148.
Huang D, Li D, Zhang H: Explicit N -fold Darboux transformation and multi-soliton solutions for the -dimensional higher-order Broer-Kaup system. Chaos Solitons Fractals 2007, 33: 1677-1685. 10.1016/j.chaos.2006.03.015
Bluman GW, Kumei S Applied Mathematical Sciences 81. In Symmetries and Differential Equations. Springer, New York; 1989.
Ibragimov NH 1. In CRC Handbook of Lie Group Analysis of Differential Equations. CRC Press, Boca Raton; 1994.
Ibragimov NH 2. In CRC Handbook of Lie Group Analysis of Differential Equations. CRC Press, Boca Raton; 1995.
Ibragimov NH 3. In CRC Handbook of Lie Group Analysis of Differential Equations. CRC Press, Boca Raton; 1995.
Olver PJ Graduate Texts in Mathematics 107. In Applications of Lie Groups to Differential Equations. 2nd edition. Springer, Berlin; 1993.
Ovsiannikov LV: Group Analysis of Differential Equations. Academic Press, New York; 1982. (English translation by W.F. Ames)
Leveque RJ: Numerical Methods for Conservation Laws. Birkhäuser, Basel; 1992.
Newell AC: The history of the soliton. J. Appl. Mech. 1983, 50: 1127-1137. 10.1115/1.3167195
Mikhailov AV, Shabat AB, Yamilov RI: On an extension of the module of invertible transformations. Sov. Math. Dokl. 1987, 295: 288-291.
Mikhailov AV, Shabat AB, Yamilov RI: Extension of the module of invertible transformations and classification of integrable systems. Commun. Math. Phys. 1988, 115: 1-19. 10.1007/BF01238850
Kodama Y, Mikhailov AV: Obstacles to asymptotic integrability. In Algebraic Aspects of Integrability. Edited by: Gelfand IM, Fokas A. Birkhäuser, Basel; 1996:173-204.
Naz R, Mahomed FM, Mason DP: Comparison of different approaches to conservation laws for some partial differential equations in fluid mechanics. Appl. Math. Comput. 2008, 205: 212-230. 10.1016/j.amc.2008.06.042
Ibragimov NH: A new conservation theorem. J. Math. Anal. Appl. 2007, 333: 311-328. 10.1016/j.jmaa.2006.10.078
Anco SC, Bluman GW: Direct construction method for conservation laws of partial differential equations. Part I: examples of conservation law classifications. Eur. J. Appl. Math. 2002, 13: 545-566.
Hereman W: Symbolic computation of conservation laws of nonlinear partial differential equations in multi-dimensions. Int. J. Quant. Chem. 2006, 106: 278-299. 10.1002/qua.20727
Kudryashov NA: Simplest equation method to look for exact solutions of nonlinear differential equations. Chaos Solitons Fractals 2005, 24: 1217-1231. 10.1016/j.chaos.2004.09.109
Kudryashov NA: Exact solitary waves of the Fisher equation. Phys. Lett. A 2005, 342: 99-106. 10.1016/j.physleta.2005.05.025
Vitanov NK: Application of simplest equations of Bernoulli and Riccati kind for obtaining exact traveling-wave solutions for a class of PDEs with polynomial nonlinearity. Commun. Nonlinear Sci. Numer. Simul. 2010, 15: 2050-2060. 10.1016/j.cnsns.2009.08.011
Vitanov NK, Dimitrova ZI: Application of the method of simplest equation for obtaining exact traveling-wave solutions for two classes of model PDEs from ecology and population dynamics. Commun. Nonlinear Sci. Numer. Simul. 2010, 15: 2836-2845. 10.1016/j.cnsns.2009.11.029
Anthonyrajah M, Mason DP: Conservation laws and invariant solutions in the Fanno model for turbulent compressible flow. Math. Comput. Appl. 2010, 15: 529-542.
Atherton RW, Homsy GM: On the existence and formulation of variational principles for nonlinear differential equations. Stud. Appl. Math. 1975, 54: 31-60.
Acknowledgements
This paper is dedicated to Prof. Ravi P. Agarwal on the occasion of his 65th birthday.
CMK would like to thank the Organizing Committee of ‘International Conference on Applied Analysis and Algebra (ICAAA2012)’ for their kind hospitality during the conference.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author declares that he has no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Khalique, C.M. On the solutions and conservation laws of the -dimensional higher-order Broer-Kaup system. Bound Value Probl 2013, 41 (2013). https://doi.org/10.1186/1687-2770-2013-41
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1687-2770-2013-41