Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 9535294, 8 pages
https://doi.org/10.1155/2019/9535294

Research Article

Investigation of Interaction Solutions for Modified Korteweg-de Vries Equation by Consistent Riccati Expansion Method

Jin-Fu Liang¹ and Xun Wang²

¹School of Physics and Electronic Science, Guizhou Normal University, Guiyang 550025, China
²Key Laboratory of Speech Acoustics and Content Understanding, Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China

Correspondence should be addressed to Jin-Fu Liang; liang.shi2007@163.com

Received 29 March 2019; Accepted 2 July 2019; Published 18 July 2019

Guest Editor: Shaowu Li

Copyright © 2019 Jin-Fu Liang and Xun Wang. 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.

Abstract.

A consistent Riccati expansion (CRE) method is proposed for obtaining interaction solutions to the modified Korteweg-de Vries (mKdV) equation. Using the CRE method, it is shown that interaction solutions such as the soliton-tangent (or soliton-cotangent) wave cannot be constructed for the mKdV equation. More importantly, exact soliton-cnoidal periodic wave interaction solutions are presented. While soliton-cnoidal interaction solutions were found to degenerate to special resonant soliton solutions for the values of modulus (n) closer to one (upper bound of modulus) in the Jacobi elliptic function, a normal kink-shaped soliton was observed for values of n closer to zero (lower bound).

1. Introduction

The direct study of exact solutions to nonlinear evolution equations (NLEEs) has received much attention from many mathematicians and physicists due to the fact that new strides in nonlinear science, which were made possible by a substantial increase in computational platforms such as Mathematica, Maple, and MATLAB, have enabled improvements in the performance of complicated and tedious numerical computational methods. Indeed, several powerful methods such as the Tanh-function method [1–3], F-expansion method [4], Jacobian elliptic function method [5], and variational approach [6, 7] have been proposed for constructing exact solutions to NLEEs. Despite the successful implementation of such methods, it is still challenging to obtain solutions for interactions among different types of nonlinear excitations such as the soliton-soliton interaction.

Recently, some new soliton structure solutions were obtained for nonlinear systems. Chen et al. studied the vortex solitons in Bose-Einstein condensates with spin-orbit coupling and Gaussian optical lattices, based on the analytical and numerical method [8]. Milan et al. found exact fundamental soliton solutions in the spiraling guiding structures by the modified Petviashvilis iteration method [9]. Cheng et al. investigated the formation and propagation of a multipole soliton in a cold atomic gas with a parity-time symmetric potential using the modified square operator method [10]. Liu et al. obtained the three-soliton solutions for high-order nonlinear Schrodinger equation by Hirotas bilinear method [11].

Specially, Lou [12] proposed a consistent Riccati expansion (CRE) method, which is a more generalized yet simpler method to find interaction solutions for various NLEEs [13–17]. The core concept of CRE is the construction of interaction solutions based on the usual Riccati equation method and the consistent equation or the -equation [12]. The CRE method is critical to finding more new solutions to the -equation.

In this study, the CRE method is used to construct several types of interaction solutions for the focusing real modified Korteweg-de Vries (mKdV) equation [18] shown in

where and are arbitrary constants. The mKdV equation plays an important role in describing some physical phenomena, such as optical cycles [19, 20], soliton propagation in plasma [21] and lattices [22], the Schottky barrier transmission lines [23], and fluid mechanics [24].

To provide better insights into these physical phenomena, finding and analyzing exact solutions to the mKdV equation is important. Previously, many powerful methods have been proposed for constructing exact solutions to the mKdV equation. For instance, in 1972, Hirota obtained an exact solution to the mKdV equation for the case of multiple collisions of solitons with different amplitudes [25]. Subsequently, he also derived the exact envelope soliton solution to the mKdV equation [26]. In 1973, Ablowitz et al. obtained exact solutions to the mKdV equation by using the inverse scattering technique [27]. In 1988, Akhmediev et al. used the Darboux transformation scheme to obtain second-order periodic solutions to the mKdV equation [28]. In 2004, Kevrekidis et al. derived some classes of periodic solutions to the mKdV equation by using direct methods [29]. In 2015, Jiao and Lou constructed a new soliton-cnoidal periodic wave interaction solution by using the CRE method [30]. However, they did not investigate how the soliton-cnoidal interaction solutions may be used to derive soliton-soliton or soliton-periodic wave interaction solutions among other types of solutions. Moreover, new interaction solutions to the mKdV equation involving different types of nonlinear waves must be investigated in depth.

The present article is structured as follows. Section 2 introduces the CRE solvability of the mKdV equation. Section 3 describes new explicit interaction solutions such as soliton-soliton, multiple resonant soliton, soliton-cosine wave, and soliton-cnoidal wave solutions to the mKdV equation obtained using the CRE method. Furthermore, it is demonstrated that interaction solutions such as the soliton-tangent wave solution cannot be constructed for the mKdV equation. The last section presents a summary and discussion.

2. CRE Solvability of the mKdV Equation

Consider the following NLEE, shown in (2), with independent variables and a dependent variable

where is a polynomial function of some arguments with the subscripts denoting partial derivatives. We assume that the solution to (2) is the following possible truncated expansion form

where is determined from the leading order analysis of (2). All the expansion coefficient functions () are determined by substituting (3) into (2) and then vanishing all the coefficients for a given power of . Further, and are functions of and satisfies the following simple Riccati equation shown below:

which includes the following five special solutions [31]. For ,

For 0,

For = 0,

Definition. If the equation for and , obtained by vanishing all the coefficients of each power in after the substitution of (3) into (2), is either consistent or not overdetermined, then the expansion in (3) is considered a CRE and the nonlinear system in (2) is said to be CRE solvable [8].

According to the CRE method defined above, one can obtain the following form based on the leading order analysis of the mKdV equation in (1)

where , , and are functions of and satisfies the Riccati equation (see (4) above).

Substituting (8) and (4) into (1) and vanishing all the coefficients of different powers of , one obtains

Based on the definition above, (11) is the consistent equation of the mKdV equation (or the mKdV -equation). If is a solution to the MDWW -equation in (11), the mKdV equation in (1) is CRE solvable. In this study, we set and . Thus, the solutions to the mKdV equation are expressed as follows.

3. Interaction Solutions to the mKdV Equation

Upon the determination of solutions to (11) by using (12), the corresponding solutions to the mKdV equation in (1) can be obtained. In this section, we construct interaction solutions to the mKdV equation by using different types of trivial solutions to (11).

3.1. Soliton-Soliton Interaction Solutions to the mKdV Equation

To obtain soliton-soliton interaction solutions to the mKdV equation, we consider the following form in (13) as the trial solution to (11):

where satisfies the following Riccati equation in

where is an arbitrary constant. This equation has special solutions similar to those in (5a), (5b), (6a), (6b), and (7). By vanishing all the coefficients for each power of after the substitution of (13) and (14) into the mKdV -equation in (11), one can obtain

From (15), it can be seen that both and r are less than 0 when 0. Based on (5a), (5b), (6a), and (6b), Eqs. (4) and (14) have only solitary solutions (viz., (5a) and (5b)) but not tangent or cotangent solutions such as the ones described in (6a) and (6b). This shows that interaction solutions such as soliton-tangent (or soliton-cotangent) wave cannot be constructed for the mKdV equation.

Under condition (15), from (5a), (5b), (12a), and (13), the soliton-soliton interaction solutions of the mKdV equation are expressed as

3.2. Interaction Solution between and Trigonometric Periodic Wave for the mKdV Equation

To investigate the interaction between a soliton and a periodic wave in the interaction solution to the mKdV equation, we consider a solution of the following form (11):

where is determined constant. Substituting (17) into (11) and then vanishing all coefficients of powers , we can obtain

Then, based on (5a), (17), and (12a), one can obtain a soliton-trigonometric periodic wave interaction solution as

where .

3.3. Interaction Solutions between Soliton and Cnoidal Wave for the mKdV Equation

In [26], Jiao and Lou constructed a solution of the following form for (11)

where is the usual Jacobi elliptic sine function and is the third type of incomplete elliptic integral. Jiao and Lou used the following parameters in (21) to obtain a special soliton-cnoidal wave interaction solution to the mKdV equation:

As seen from (21), Jiao and Lou chose the modulus () of Jacobi elliptic function to be 1.5, which is outside the allowed range [13].

In this study, we will further investigate how soliton-cnoidal interaction solutions can be used to derive soliton-soliton and soliton-periodic wave interaction solutions among other types of solutions. To this end, we performed all the substitutions and evaluations by using the Mathematica software.

Consider a trial solution of the following form for solving (11)

where

satisfies the following elliptic equation:

Substituting (22) and (24) into (11), one obtains

From the analysis of (24), we assume the solution of (24) in the following form

Substituting (26) into (24) and setting the coefficient of , cn, equal to zero, one obtains

Based on (25) and (27), one can find a group solution

where , and . Under the substitution of (26), (22), and (5a) into (12a), under the conditions imposed by (28), one obtains

where

To investigate how the soliton-cnoidal interaction solutions could be used to derive soliton-soliton interaction or other types of solutions, we illustrate the following two cases corresponding to the soliton-cnoidal wave interaction solution described in (29) by selecting different sets of parameters. For the first case, the parameters are chosen as

While Figure 1 shows two-dimensional views for interaction solution at and . Figure 2 displays three-dimensional plots for the evolution of soliton-cnoidal wave interaction solution with different values for the modulus in the Jacobian elliptic function, viz., , 0.5, and 0.99999. While exhibited a particular periodic-kink soliton wave interaction, the extreme values of (a value close to the lower modulus limit or 0) showed a normal kink-shaped soliton and (a value close to the upper modulus limit or 1) displayed an interaction between a periodic wave and another periodic wave.

(a)

(b)

Figure 1: (Color online) The evolution of the soliton-cnoidal wave interaction solution shown in (29) obtained by using the parameters listed in (31). The profiles at t = 0 and x = 0 are shown in (a) and (b), respectively. Different color traces represent different values for the modulus n = 0.00001 (black), n = 0.5 (red), and n = 0.99999 (blue).

(a)

(b)

(c)

Figure 2: (Color online) The evolution of the soliton-cnoidal wave interaction solution shown in (29) obtained by using the parameters listed in (31) with different modulus values as a function of both x and t. The modulus values for different panels are (a) n = 0.00001, (b) n = 0.5, and (c) n = 0.99999.

For the second case, the parameters were altered as shown in (32) by changing the angular frequency

Similar to the first case, we illustrate the structures of the soliton-cnoidal wave interaction solution for different values of , 0.5, and 0.99999. Clearly, as shown in Figures 3 and 4, wavenumbers and the amplitudes in the range of x(-40,40) and t(-40,40) are less than that of the first case (cf. Figures 1 and 2 ). While there is still a normal kink soliton in the x-u plot for n = 0.00001, an incomplete kink soliton is observed in the t-u plot in contrast to the first case shown in Figure 1(b). Building on the above two cases, soliton and soliton-soliton wave interaction solutions are derived from the soliton-cnoidal wave interaction solution by making the limit of the modulus approach either 0 or 1.

(a)

(b)

Figure 3: (Color online) The evolution of the soliton-cnoidal wave interaction solution shown in (29) obtained by using the parameters in (32). The profiles at t = 0 and x = 0 are shown in (a) and (b), respectively. Different color traces represent different values for the modulus n = 0.00001 (black), n = 0.5 (red), and n = 0.99999 (blue).

(a)

(b)

(c)

Figure 4: (Color online) The evolution of the soliton-cnoidal wave interaction solution shown in (29) obtained by using the parameters listed in (32) with different modulus values as a function of both x and t. The modulus values for different panels are (a) n = 0.00001, (b) n = 0.5, and (c) n = 0.99999.

4. Summary and Discussion

In this study, we investigated the focusing mKdV equation by using the CRE method. This nonlinear equation was shown to be CRE solvable and interaction solutions; namely, soliton-soliton, soliton-trigonometric periodic waves, and soliton-cnoidal periodic wave for the mKdV equation were explicitly provided by choosing different trial solutions for the mKdV w-equation shown in (11). In addition, analytical solutions for interactions between soliton and cnoidal wave were provided and their properties were discussed graphically. According to the presented analysis, soliton and soliton-soliton wave interaction solutions can be derived from the soliton-cnoidal wave interaction solution by making the limit of the modulus approach either 0 or 1 (i.e., lower or upper bounds for the modulus in the Jacobi elliptical function).

Data Availability

No data were used to support this study.

Conflicts of Interest

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 11864007 and No. 11564006), the Doctor Foundation of Guizhou Normal University (2015), and the Science and Technology Planning Project of Guizhou Province (Grant No. 20185769).

References

H. A. Abdusalam, “On an improved complex tanh-function method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 99–106, 2005.

M. A. Abdou and A. A. Soliman, “Modified extended tanh-function method and its application on nonlinear physical equations,” Physics Letters A, vol. 353, no. 6, pp. 487–492, 2006.

D. L. Sekulic, M. V. Sataric, and M. B. Zivanov, “Symbolic computation of some new nonlinear partial differential equations of nanobiosciences using modified extended tanh-function method,” Applied Mathematics and Computation, vol. 218, no. 7, pp. 3499–3506, 2011.

W. W. Li, Y. Tian, and Z. Zhang, “F method and its application for finding new exact solutions to the sine-Gordon and sinh-Gordon equations,” Applied Mathematics and Computation, vol. 219, no. 3, pp. 1135–1143, 2012.

A. H. Bhrawy, M. A. Abdelkawy, and A. Biswas, “Cnoidal and snoidal wave solutions to coupled nonlinear wave equations by the extended Jacobi's elliptic function method,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 4, pp. 915–925, 2013.

H. Liu, “Variational approach to nonlinear electrochemical system,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 5, no. 1, pp. 95-96, 2004.

J. F. Liang and L. X. Gong, “Complex wave solutions for (2+1)-dimensional modified dispersive water wave system,” Communications in Theoretical Physics, vol. 52, no. 1, pp. 17–22, 2009.

S. F. Chen, Q. Guo, S.-L. Xu et al., “Vortex solitons in bose-einstein condensates with spin-orbit coupling and gaussian optical lattices,” Applied Mathematics Letters, vol. 92, pp. 15–21, 2019.

S. P. Milan, S. Aleksandra, A. Branislav, and R. Milivoj, “Rotating solitons supported by a spiral waveguide,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 98, no. 5, Article ID 063822, 2017.

J. Cheng, S. Xu, M. R. Belić et al., “Multipole solitons in a cold atomic gas with a parity-time symmetric potential,” Nonlinear Dynamics, vol. 95, no. 3, pp. 2325–2332, 2019.

X. Liu, H. Triki, Q. Zhou et al., “Generation and control of multiple solitons under the influence of parameters,” Nonlinear Dynamics, vol. 95, no. 1, pp. 143–150, 2019.

S. Y. Lou, “Consistent Riccati expansion for integrable systems,” Studies in Applied Mathematics, vol. 134, no. 3, pp. 372–402, 2015.

X. Liu, J. Yu, and Z. Lou, “New interaction solutions from residual symmetry reduction and consistent Riccati expansion of the (2+1)-dimensional Boussinesq equation,” Nonlinear Dynamics, vol. 92, no. 4, pp. 1469–1479, 2018.

L. Huang, Z. Qiao, and Y. Chen, “Soliton–cnoidal interactional wave solutions for the reduced Maxwell–Bloch equations,” Chinese Physics B, vol. 27, no. 2, Article ID 020201, 2018.

J. Chen, H. Wu, and Q. Zhu, “Bäcklund transformation and soliton–cnoidal wave interaction solution for the coupled Klein–Gordon equations,” Nonlinear Dynamics, vol. 91, no. 3, pp. 1949–1961, 2018.

Y. Li and H. Hu, “Nonlocal symmetries and interaction solutions of the Benjamin–Ono equation,” Applied Mathematics Letters, vol. 75, pp. 18–23, 2018.

J. Liang and X. Wang, “Consistent Riccati expansion for finding interaction solutions of (2+1)‐dimensional modified dispersive water‐wave system,” Mathematical Methods in the Applied Sciences, 2019.

S. Liu and S. Liu, Nonlinear Equation in Physics, Peking University Press, Beijing, China, 2000.

H. Leblond and D. Mihalache, “Few-optical-cycle solitons: modified korteweg–de vries sine-gordon equation versus other non–slowly-varying-envelope-approximation models,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 79, no. 6, Article ID 063835, 2009.

H. Leblond and D. Mihalache, “Few-optical-cycle dissipative solitons,” Journal of Physics A: Mathematical and Theoretical, vol. 43, no. 37, Article ID 375205, 2010.

N. J. Zabusky and M. D. Kruskal, “Interaction of solitons in acollisionless plasma and the recurrence of initial states,” Physical Review Letters, vol. 15, no. 6, pp. 240–243, 1965.

H. Ono, “Soliton fission in anharmonic lattices with reflectionless inhomogeneity,” Journal of the Physical Society of Japan, vol. 61, no. 12, pp. 4336–4343, 1992.

V. Ziegler, J. Dinkel, C. Setzer, and K. E. Lonngren, “On the propagation of nonlinear solitary waves in a distributed Schottky barrier diode transmission line,” Chaos Solitons and Fractals, vol. 12, no. 9, pp. 1719–1728, 2001.

M. A. Helal, “Soliton solution of some nonlinear partial differential equations and its applications in fluid mechanics,” Chaos Solitons and Fractals, vol. 13, no. 9, pp. 1917–1929, 2002.

R. Hirota, “Exact solution of the korteweg—de vries equation for multiple Collisions of solitons,” Physical Review Letters, vol. 27, no. 18, pp. 1192–1194, 1971.

R. Hirota, “Exact envelope-soliton solutions of a nonlinear wave equation,” Journal of Mathematical Physics, vol. 14, pp. 805–809, 1973.

M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Segur, “Nonlinear-evolution equations of physical significance,” Physical Review Letters, vol. 31, no. 2, pp. 125–127, 1973.

N. Akhmediev, V. I. Korneev, and N. V. Mitskevich, “N-modulation signals in a single-mode optical fiber with allowance for nonlinearity,” Journal of Experimental and Theoretical Physics, vol. 94, pp. 159–170, 1988.

P. G. Kevrekidis, A. Khare, A. Saxena, and G. Herring, “On some classes of mKdV periodic solutions,” Journal of Physics A: Mathematical and General, vol. 37, no. 45, pp. 10959–10965, 2004.

X. Jiao and S. Lou, “CRE method for solving mKdV equation and new interactions between solitons and cnoidal periodic waves,” Communications in Theoretical Physics, vol. 63, no. 1, pp. 7–9, 2015.

S. H. Ma, J. Y. Fang, and C. L. Zheng, “New exact solutions for the (3+1)-dimensional Jimbo-Miwa system,” Chaos Solitons and Fractals, vol. 40, no. 3, pp. 1352–1355, 2009.