We aim to construct exact and explicit solutions to a generalized Bogoyavlensky-Konopelchenko equation through the Maple computer algebra system. The considered nonlinear equation is transformed into a Hirota bilinear form, and symbolic computations are made for solving both the nonlinear equation and the corresponding bilinear equation. A few classes of exact and explicit solutions are generated from different ansätze on solution forms, including traveling wave solutions, two-wave solutions, and polynomial solutions.

1. Introduction

One of the fundamental problems in the theory of differential equations is to determine a solution of a differential equation satisfying what are known as initial values. There are two systematical approaches in a linear world: Laplace’s method for solving linear ordinary differential equations and the Fourier transform method, for linear partial differential equations [1, 2]. In the modern theory of integrable systems, the isomonodromic transform method and the inverse scattering transform method have been created for attempting initial value problems for nonlinear ordinary and partial differential equations, respectively [3, 4].

However, only the simplest differential equations, often linear, are solvable precisely. It is definitely not an easy task for us to find exact solutions to nonlinear differential equations, either ordinary or partial. Nevertheless, the Lie group method and the Hirota bilinear method are among effective approaches for finding exact solutions to nonlinear differential equations. The Lie group method is to determine Lie symmetries, which are used to solve ordinary differential equations immediately or used to reduce partial differential equations and to solve simpler reduced ones [5]. The Hirota bilinear method is to transform differential equations into bilinear counterparts and then solve the resulting bilinear ones [6, 7].

Based on Hirota bilinear forms, one can find solitons—a kind of analytic solutions exponentially localized [3, 4]. Some recent studies have also been made on another kind of interesting explicit solutions called lumps, originated from solving integrable equations [815]. Lump solutions are a class of analytical rational function solutions localized in all directions in space [9]. A bilinear framework for getting soliton solutions in the (2+1)-dimensional case is as follows. Suppose that a determines a Hirota bilinear form where , , and are Hirota’s bilinear derivatives, for a given (2+1)-dimensional partial differential equation:Through the Hirota bilinear technique, soliton solutions can be often formulated as follows: where stands for the sum over all possibilities for taking either 0 or 1 and the wave variables and the phase shifts are determined byandwith , , and satisfying the corresponding dispersion relation and being arbitrary translation shifts. Solitons contain various kinds of exact solutions to integrable equations, and taking long wave limits of -soliton solutions can generate special lumps [16].

In this paper, we would like to look for exact solutions to a (2+1)-dimensional generalized Bogoyavlensky-Konopelchenko equationwhere , and , , , , and are constant coefficients, through Maple symbolic computations. Based on its Hirota bilinear form, a few solution ansätze will be analyzed to compute exact solutions to the nonlinear equation and its bilinear counterpart. Moreover, starting from the nonlinear equation itself, we will do a thorough symbolic computation by Maple within our capacity to generate a few classes of exact and explicit solutions, including traveling wave solutions, two-wave solutions, and polynomial solutions. Conclusions and remarks will be given in the last section.

2. One-Wave Type and Two-Wave and Polynomial Solutions

Let us consider the (2+1)-dimensional generalized Bogoyavlensky-Konopelchenko (gBK) equation (6). Equivalently, the (2+1)-dimensional gBK equation (6) can be written as follows: This is a generalization of the (2+1)-dimensional Bogoyavlensky-Konopelchenko (BK) equation a special case of which was introduced as a (2+1)-dimensional version of the KdV equation in [17] and described as the interaction of a long wave propagating along the -axis and a Riemann wave propagating along the -axis [18]. For the (2+1)-dimensional BK equation (8), a Darboux transformation has been given, together with some traveling wave solutions [19], and a few particular properties have been explored (see, e.g., [20, 21]).

By a direct computation, we can show that the (2+1)-dimensional gBK equation (6) can be written as a Hirota bilinear form [22]: under the logarithmic transformations Such transformations play a prominent role in Bell polynomial theories for soliton equations and their generalized counterparts (see, e.g., [23, 24]). Precisely, we can haveAnd, thus, when solves the bilinear gBK equation (9), and will solve the nonlinear gBK equation (6).

Beginning with the gBK equation (6) and its Hirota bilinear form (9), we can compute various exact solutions to the gBK equation (6), through carrying out our searches via symbolic computations. For example, lower-order lump solutions have been presented to the gBK equation (6) in [22].

In what follows, we start with some special ansätze for and to construct new exact and explicit solutions to the gBK equation (6) (or equivalently (7)) and the bilinear gBK equation (9). Plugging each ansatz into the nonlinear or bilinear gBK equation leads to a system of algebraic equations on the parameters and the coefficients. Then, conduct symbolic computations with Maple to obtain solutions to the algebraic system and further exact solutions to the (2+1)-dimensional gBK equation (6).

2.1. One-Wave Type Solutions

First, trying one-wave type solutions, we can easily determine the following five exact and explicit solutions to the bilinear gBK equation (9):where , are arbitrary functions, = sin or cos, and the constant parameters , are arbitrary provided that every term in the solutions makes sense. Though all classes of solutions above are interesting solutions to the bilinear gBK equation (9), only the first class of solutions can lead to nontrivial exact solutions to the nonlinear gBK equation (6). Amazingly, the first class of solutions presents traveling wave solutions involving an arbitrary function for the (2+1)-dimensional gBK equation (6), which is a special characteristic of the the gBK equation (6) [22]. Also from the above first class of exact solutions, we can easily formulate various lump-type solutions (i.e., rational and analytical function solutions that are localized in almost all directions in space, under the Lebesgue measure) to the gBK equation (6), by taking to be positive polynomial functions.

2.2. Two-Wave Solutions

To search for two-wave solutions, let us now set where , are constant parameters to be determined. Try an ansatz of two-wave solutions to the bilinear gBK equation (9): where is another constant parameter to be determined. and we can show that the resulting system of algebraic equations has two classes of explicit solutions:andwhere the constant satisfies and the other parameters could be arbitrary provided that the solutions of and presented by (10) will be well defined.

If we try another ansatz for two-wave solutions: and then we can havewhere the constant needs to satisfy a quadratic equation This is different from (17). When , these are two quadratic equations and so each of them has two roots for ; but when , they become linear and so only one solution for is possible.

Then, through the transformations in (10), we can obtain diverse classes of exact and explicit two-wave solutions to the (2+1)-dimensional gBK equation (6).

2.3. Polynomial Solutions

Let us thirdly try an ansatz on polynomial solutions for :where and , are constants coefficients to be determined. By symbolic computations for the case of , we can get the following two classes of exact polynomial solutions to the gBK equation (7):andwhere the constant is given byand the other involved parameters are arbitrary.

The first class of polynomial solutions does not depend on any coefficient of the nonlinear terms in the gBK equation (7), but the second class depends critically on the coefficient of the second group containing nonlinear terms in the gBK equation (7) and involves two more monomials: and . It is also surprising that the cases of are just too complicated for us to work out any nontrivial explicit polynomial solutions.

Taking in the first solution (22) andin the second solution (23), we obtain two quadratic function solutions immediately: andwhere is defined by (24).

2.4. An Illustrative Example

Let us now takeThen, from (7), we obtain the following specific gBK equation Based on our presented solutions, this nonlinear equation has a traveling wave solution a two-wave solution and a polynomial solution

All the solutions computed above via symbolic computations supplement the solution theories available on soliton solutions and dromion-type solutions, generated through powerful existing approaches such as the Hirota perturbation method and symmetry constraints including symmetry reductions (see, e.g., [2530]).

3. Concluding Remarks

We have constructed a few classes of exact and explicit solutions, including two-wave solutions and polynomial solutions, to a (2+1)-dimensional generalized Bogoyavlensky-Konopelchenko (gBK) equation. Symbolic computations with Maple are the adopted technique and the Hirota bilinear form is a basis for getting one-wave type and two-wave solutions.

The obtained result enriches the existing studies on the (2+1)-dimensional Bogoyavlensky-Konopelchenko equation [1922]. We remark that there are plenty of interaction solutions to integrable equations (see, e.g., [31]), particularly between lumps and other kinds of exact solutions to (2+1)-dimensional nonlinear integrable equations (see, e.g., [3235] for lump-kink interaction solutions and [3639] for lump-soliton interaction solutions).

However, we couldn’t find any interaction solutions between lump (or lump-type) solutions and kink (or soliton) solutions, and any nontrivial -wave solutions, including -soliton solutions, where , for the (2+1)-dimensional gBK equation (6). We believe that the existence of interaction solutions and three-wave solutions to partial differential equations should be a characteristic property of complete integrability for the differential equations under consideration. It is definitely interesting to search for exact solutions and interaction solutions to partial differential equations in both (1+1)- and (2+1)-dimensions. Another important problem is to identify nonlinear partial differential equations that possess soliton solutions and interaction solutions, to understand integrable properties of nonlinear equations.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


The work was supported in part by NSFC under Grants 11301454, 11301331, 11371086, 11571079, and 51771083, NSF under Grant DMS-1664561, the Jiangsu Qing Lan Project for Excellent Young Teachers in University (2014), Six Talent Peaks Project in Jiangsu Province (2016-JY-081), the Natural Science Foundation for Colleges and Universities in Jiangsu Province (17KJB110020), Emphasis Foundation of Special Science Research on Subject Frontiers of CUMT under Grant No. 2017XKZD11, and the Distinguished Professorships by Shanghai University of Electric Power, China, and North-West University, South Africa.