#### Abstract

A general coupled KdV equation, which describes the interactions of two long waves with different dispersion relation, is considered. By employing the Hirota’s bilinear method, the bilinear form is obtained, and the one-soliton solution and two-soliton solution are constructed. Moreover, the elasticity of the collision between two solitons is proved by analyzing the asymptotic behavior of the two-soliton solution. Some figures are displayed to illustrate the process of elastic collision.

#### 1. Introduction

In the soliton theory, finding the soliton solutions for the nonlinear partial equations is becoming more and more important since the soliton solutions can describe many complex physical phenomena [1]. Many effective approaches have been proposed, such as the inverse scattering transformation method [2], the Bäcklund transformation method [3], the Darboux transformation method [4–6], the Hirota bilinear method [7–11], and the Riemann-Hilbert method [12]. Among them, the Hirota bilinear method is not only direct but also effective for investigating the soliton solutions.

In the past decades, the coupled Korteweg-de Vries (KdV) equations have been investigated widely and many integrable coupled KdV equations are found. For example, Gurses and Karasu [13] showed that the following coupled KdV equation was integrable and admitted recursion operator and a bi-Hamiltonian structure:

In fact, this equation is Lax integrable, and the Lax pair was firstly given by Drinfeld and Sokolov [14] and then by Bogoyavlenskii [15] and Karasu and Yurduşen [16] independently. Subsequently, this equation was also derived by Satsuma and Hirota [17] as one case of the four-reduction of the KP Hierarchy. Moreover, Karasu and Yurduşen [16] proposed a Bäcklund transformation and some explicit solutions of Equation (1). As far as we know, the soliton solutions and the collision between two solitons have not been investigated. So in this paper, we investigate the following general coupled KdV equation: which is just the coupled KdV Equation (1) for

The rest of this paper is organized as follows. In “The Bilinear Form and Soliton Solutions,”, the bilinear form, the one-soliton, and two-soliton of Equation (2) are obtained based on the Hirota’s direct method. In “Asymptotic Analysis on Two-Soliton Solution,” the asymptotic behaviors are studied to prove that the two-soliton collision is elastic. Finally, conclusions are given in “Conclusion.”

#### 2. The Bilinear Form and Soliton Solutions

We implement the following dependent variable transformation to Equation (2): where and are functions of and Then, the following bilinear equations of Equation (2) are obtained as follows: where the -operators [18] are defined by where both and are integers.

In order to apply the perturbation method to Equation (4) to find the soliton solutions of Equation (2), we expand functions and in power series of a small parameter as where and are functions of and . Substituting Equations (6) and (7) into the bilinear Equation (4) and collecting the coefficients of parameter , we have

##### 2.1. One-Soliton Solution

To obtain the one-soliton solution for the general coupled KdV Equation (2), set where . Substituting it into Equation (8), we have and

Furthermore, from Equation (9), we have

Assuming , it is easy to see that Equation (10) and other equations from the coefficients of parameter are satisfied automatically. So we obtain the following one-soliton solution for the general coupled KdV Equation (2) by setting

Figures 1(a) and 1(b) demonstrate the soliton structures of one-solutions and , respectively, for parameters .

**(a)**

**(b)**

##### 2.2. Two-Soliton Solution

Likewise, to arrive the two-soliton solution, we set where , . Plugging Equation (17) into Equation (8) yields , , and with .

From Equations (9) and (18), we have with

Similarly, from Equations (10), (18), and (19), and can be derived as

with

Moreover, if we plug the above obtained and into Equation (11), we get with .

Assuming , it is observed that the left equations are satisfied automatically. By setting , we obtain the two-soliton solutions for the coupled KdV Equation (2). We i4llustrate the structures of and in Figures 2(a) and 2(b).

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 3. Asymptotic Analysis on Two-Soliton Solution

Now we analyze the two-soliton solution of Equation (2) with long-time asymptotic method. Note that the two-soliton solution can be written as where and are shown in Equations (17)–(23). Without loss of generality, suppose that .

For fixed , note that , then arrive at the following: (i)Solitons-1 before collision (, ) where (ii)Solitons-1 after collision (, ) where .

For fixed , note that , and then we arrive at the following: (i) Solitons-2 before collision (, ) where (ii) Solitons-2 after collision (, ) where .

The above asymptotic analysis can also be seen in Figure 3. Comparing the asymptotic expressions Equation (25) with Equation (26) and Equation (27) with Equation (28), we find that the amplitudes and velocities remain the same, but the phases are changed. To illustrate the collision process exactly, the graphs are presented in Figure 2, which shows that the collisions of the two-soliton waves are exactly elastic.

#### 4. Conclusion

In conclusion, we studied a general coupled KdV Equation (2)via the Hirota’s bilinear method. We first constructed the bilinear form and then the one-soliton solution and two-soliton solution. Furthermore, the asymptotic analysis is given to prove that the collision of the two-soliton solutions is elastic.

#### Data Availability

All data, models, and code generated or used during the study appear in the submitted article.

#### Conflicts of Interest

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

#### Authors’ Contributions

Changhao Zhang is responsible for methodology, data curation, software, and validation. Guiying Chen is responsible for conceptualization, writing-original draft, visualization, writing (review and editing), supervision, and project administration. All authors have read and agreed to the published version of the manuscript.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China under grant No. 62071053 and PHD Research Foundation of Liaocheng University under grant No. 318051431. The authors would like to express their sincere thanks to the editors and reviewers for the kind comments and valuable suggestions.