# Residual Symmetry, Bäcklund Transformation, and Soliton Solutions of the Higher-Order Broer-Kaup System

Authors: Yarong Xia; Ruoxia Yao; Xiangpeng Xin
Journal: Advances in Mathematical Physics
 (2021)

---

## Abstract

Under investigation in this paper is the higher-order Broer-Kaup(HBK) system, which describes the bidirectional propagation of long waves in shallow water. Via the standard truncated Painlevé expansion method, the residual symmetry of this system is derived. By introducing an appropriate auxiliary-dependent variable, the residual symmetry is successfully localized to Lie point symmetries. Via solving the initial value problems, the finite symmetry transformations are presented. However, the solution which obtained from the residual symmetry is a special group invariant solutions. In order to find more general solution of HBK system, we further generalize the residual symmetry method to the consistent tanh expansion (CTE) method and prove that the HBK system is CTE solvable, then the resonant soliton solutions and interaction solutions among different nonlinear excitations are obtained by the CET method.


---

## Body

## 1. Introduction

Since the Lie group theory was proposed by Sophus Lie to study differential equation, the symmetry theory has been widely developed to study nonlinear equation. As is known, thanks to the classical or nonclassical Lie group method [1, 2], Lie point symmetries of a differential system can be obtained, from which one can reduce the dimensions of partial differential equations and proceed to construct group invariant solutions by similarity reductions [3–8]. Nevertheless, for the integrable system, there may exist nonlocal symmetries which are related to inverse recursion operators [9, 10], Darboux transformation(DT) [11], Bäcklund transformation(BT) [12–15], conformal invariance [16], negative hierarchies [17, 18], and so on. Recently, it is found that Painlevé analysis can also be applied to obtain nonlocal symmetries, and this type of nonlocal symmetries corresponds to the residues with respect to the singular manifold of the truncated Painlevé expansion, which are also called residual symmetries [19–23].One knows that the nonlocal symmetries cannot be used to construct explicit solutions to differential equations directly; naturally, we need to transform the nonlocal symmetries into local ones by introducing suitable prolonged systems. In reference [24, 25], the authors has considered the localizations of the residual symmetries and found related finite transformation. Furthermore, by developing the truncated Painlevé expansion, Lou introduced the definition of consistent Riccati expansion solvable [26], and this method is greatly valid for constructing both possible new integrable systems and interaction solutions between a soliton and other types of nonlinear excitations. A consistent tanh expansion is a special simplified form of the consistent Riccati expansion defined in literature [27]. So, through the relationship of the nonlocal symmetry and the consistent Riccati expansion solvable or the consistent tanh expansion solvable, one can construct the solutions of nonlinear evolution equations. On account of this, there are a lot of papers here to study this problem [28–32].This paper concentrates on investigating the following higher-order Broer-Kaup system(1)ut+4uxx+u3+6uv−3uuxx=0,vt+4vxx+3uvx+3u2v+3v2x=0,which was first derived and studied by Lou and Hu [33]. The system (1) can be seen as an extension of the known Broer-Kaup system which is often used to model the bidirectional propagation of long waves in shallow water. Some special solutions to the system (1) were given in reference [34], huang [35] studied the explicit N-fold Darboux transformation and multisoliton solutions to the system (1), Li et al. [36] investigated the Painlevé analysis and new analytic solutions to the system (1), and the interaction solutions to the system (1) were derived in reference [37]. However, the residual symmetries, the consistent tanh expansion solvable, and interaction solutions of the higher-order Broer-Kaup system have not yet been studied, which is the prime objective of this paper.


## 2. The Residual Symmetry of the Higher-Order Broer-Kaup System and Its Localization

The higher-order Broer-Kaup system is Painlevé integrable, and its truncated Painlevé expansion can be expressed as(2)u=u0+u1ϕ,v=v0+v1ϕ+v2ϕ2,where u0,u1,v0,v1,v2 are functions of x and t, and ϕx,t=0 is the singular manifold equation. Substituting equation (2) into equation (1), then, collecting the power of ϕ and making its coefficients equal to zero, we get two systems of ϕ, and the concert expressions are in the appendix.Vanishing the coefficients of1/ϕ4 in the first one of system (1) and the coefficients of 1/ϕ5 in the second one of system in the appendix, we have
(3)u1=ϕx,v2=−ϕx2.Vanishing the coefficients of1/ϕ3 in the first one of the system (1) and the coefficients of 1/ϕ4 in the second system in the appendix and using (3) leads to
(4)v1=ϕxx.Vanishing the coefficients of1/ϕ2 in the first one of the system (1) and the coefficients of 1/ϕ3 in the second system in the appendix, then, by using (3), (4) leads to
(5)ϕt=−12ϕxu02−12ϕxxu0−24ϕxv0+12ϕxu0,x−4ϕxxx,v0=u0,x,using equations (3), (4), and (5), it is not difficult to find that the coefficient of 1/ϕ2 in the second system of appendix identically becomes zero. The coefficients of ϕ0 in appendix are nothing but the higher-order Broer-Kaup system (1) for u=u0, v=v0, while the coefficient of 1/ϕ in appendix is just the symmetry equation for system (1) with the symmetry
(6)σu0=u1,σv0=v1.That is to say, whenu0,v0 are the solutions of system (1), u1,v1 are the solutions to the corresponding symmetry equation. From above analysis and the standard truncated Painlevé expansion (2), the following Bäcklund transformation theorem and residual symmetry theorem arise.Theorem 1.
If the functionϕ satisfies equation (5), then
(7)u=u0+ϕxϕ=u0+lnϕx,v=v0+ϕxxϕ−ϕx2ϕ2=v0+lnϕxx,are the Bäcklund transformation between solutions u0,v0 and u,v.Theorem 2.
The higher-order Broer-Kaup system (1) has the following residual symmetry
(8)σu=ϕx,σv=ϕxx,where u,v, and ϕ satisfy the Bäcklund transformation (7).Proof.
From the coefficients of1/ϕ in appendix, we find that u1=ϕx,v1=ϕxx is just the solution of the symmetry equation for the system (1), where u1,v1 are just the residual of truncated Painlevé expansion (2) when u=u0,v=v0 are the solution of system (1); so based on the definition of residual symmetry, the nonlocal symmetry (8) is also called residual symmetry.Sinceu1 and v1 are nonlocal symmetries, as we all know, the nonlocal symmetries cannot construct explicit solutions to differential equations directly. One can naturally believe that we need to transform the nonlocal symmetries into local ones. To this end, we introduce new variables f and g to eliminate the space derivatives of ϕ by
(9)f=ϕx,g=fx,then the nonlocal symmetry of the higher-order Broer-Kaup system (1) is localized to the following Lie point symmetries
(10)σu=f,σv=g,σϕ=−ϕ2,σf=−2ϕf,σg=−2f2−2ϕg,for the prolonged systems (1), (7), and (9) with the Lie point symmetry vector
(11)V=f∂u−g∂v−ϕ2∂ϕ−2ϕf∂f−2f2+ϕg∂g.To proceed, we study the finite symmetry transformation of Lie point symmetries (10). According to Lie’s first theorem, by solving the following initial value problems,
(12)du^εdε=f^ε,u^0=u,dv^εdε=g^ε,v^0=v,df^εdε=−2f^εϕ^ε,f^0=f,dg^εdε=−2f∧2ε−2ϕ^εg^ε,g^0=g,dϕ^εdε=−ϕ∧2ε,ϕ^0=ϕ,one can easily obtain the symmetry group transformation theorem as follows:Theorem 3.
Ifu,v,ϕ,f,g is a solution to the prolonged system (1), (7) ,and (9), then so u^,v^,ϕ^,f^,g^ is given by
(13)u^=−fεϕ+1ϕ+uϕ+fϕ,ϕ^=ϕεϕ+1,f^=fεϕ+12,v^=vϕ2+gϕ−f2ε2+2ϕv+gε+v1+εϕ2,g^=gεϕ+12−2f2εεϕ+13.Remark 4.
From Theorem3, we find an interesting phenomenon that the residual symmetry {u1=ϕx,v1=ϕxx} coming from the truncated Painlevé expansion is just the infinitesimal form of the group. Actually, the above group transformation is equivalent to the truncated Painlevé expansion (2) and (3) since the singularity manifold equations (1), (7), and (9) are form invariant under the transformation
(14)1+εϕ→ϕwithεf→ϕx,εg→ϕxx.Remark 5.
For given solutionsu,v of (1), through the finite symmetry transformation (13) will arrive new solutions u^,v^. For example, we take a simple solutions as u=0,v=0 for (1), and then from equation (5), through solving an ordinary differential equation, we get ϕ; further, by equation (9), we get the expression of f, and the concert form of ϕ and f is as follows:
(15)ϕ=C1exp−12b11/3i3+1x+C2exp12b11/3i3−1x+C3b11/3exp3xC4exp−4b1t,f=−12C1b11/3i3+1exp−12b11/3i3+1x+12C2b11/3i3−1xexp12b11/3i3−1x+C3b11/3exp3xC4exp−4b1t,
then, by help of equation (13), we can get a pair of new solution u^ and v^ (the concrete expressions are omitted here since of its prolixity), where Ci,i=1,2,3,4 and b1 are arbitrary constants.However, the solution obtained from (13) is only a special group invariant solutions. Next, in order to obtain more general solutions of equation (1), we generalize the residual symmetry method.In reference [26], Lou start from the truncated Painlevé expansion to propose the definition of consistent Riccati expansion (CRE) solvable. Inspired by this reference, we consider the special case of CRE-consistent Tanh expansion (CTE), which is a more generalized but much simpler method to find interaction solutions between solitons and other nonlinear excitations, such as soliton-resonant solutions, soliton and condial wave, and soliton and sin-cosine wave [38–40].


## 3. Consistent Tanh Expansion Solvability of the Higher-Order Broer-Kaup System

In order to derive the interaction solution of the higher-order Broer-Kaup system (1) by the consistent tanh expansion method, we need to discuss the consistent tanh expansion solvability of the system (1). The first step is to make the leading analysis of the system (1), by balancing the highest order of tanhw, and one can obtain the following truncated tanh function expansions
(16)u¯=u¯0+u¯1tanhw+wx,v¯=v¯0+v¯1tanhw+v¯2tanh2w.For convenience later, in equation (16), we substitute u¯0 for u¯0+wx (see (19), (20) where u¯0 equation has been decoupled). By vanishing all the coefficients of the like powers of tanhw after substituting (16) into (1), we obtain eleven overdetermined equations for only six unknown functions u¯0,u¯1,v¯0,v¯1,v¯2, and w. Vanishing the coefficients of tanh4w,tanh5w leads to
(17)u¯1=wx,v¯2=−wx2,andu¯1=2wx,v¯2=−2wx2,which have to be considered as two possible case separately. To proceed, we consider the first case.Case 6.
Principle branch.u¯1=wx,v¯2=−wx2.In this case, by substituting the expression {u¯1=wx,v¯2=−wx2} into the remained nine equations, one can find that the variables u¯0,v¯0,v¯1, and w are consistent and the concert form reading
(18)v¯1=wxx,v¯0=u¯0,x+wx2+wxx,with u¯0,w satisfying the following two equations
(19)PSTO≔wt+4wxx+3wx2+3u¯0wxx+16wx3+12u¯0wxu¯0+2wx=0,(20)STO≔u¯0,t+4u¯0,xx+3u¯0u¯0,x+u¯03x=0.From (20), we found that wx has been separated out from u¯0 so that the u¯0 equation is independent w, and w equation (19) can also be linearized because it is the potential form of the variable coefficient STO (Sharma-Tasso-Olver) equation—that is PSTO. Then, we have the consistent tanh expansion solvable theorem for the system (1).Theorem 7.
The higher-order Broer-Kaup system (1) is consistent tanh expansion solvable with the consistent tanh expansion as follows:
(21)u=u¯0+wx1+tanhw,v=u¯0,x+wxx1+tanhw+wx21−tanh2w,and the u¯0,w satisfies equations (19) and (20).Case 8.
Auxiliary branch.u¯1=2wx,v¯2=−2wx2.For this case, it is easy to find that the four variablesv¯1,v¯0,u¯0, and w are not consistent, and the concrete forms of u¯,v¯ are
(22)u¯=wx1+2tanhw+wxxwx,v¯=wxx2wx2−wxxxwx+2wxxtanhw+2wx21−tanh2w,while w satisfies the following equations:
(23)wt+44wx3−2wxxx+3wxx2wx=0,4wx2−∂x2wxxxwx−32wxx2wx2+10wxx2=0.It is not difficult to check that this case is not a consistent tanh expansion solvable. However, it is still a useful way to construct the special type of exact solution to the higher-order Broer-Kaup system (1) by solving the overdetermined equations (23).


## 4. Explict Solution to the Higher-Order Broer-Kaup System

The explict solutions of nonlinear PDE play an important role in nonlinear science and engineering, and there are many kinds of effective methods that have been established to construct explict solutions for nonlinear equations, such as symmetry reductions [30, 31], bilinear method [41–43], Darboux transformation [44], and Painlevé analysis [45], but most of them are more difficult to find interaction solutions which are an important and meaningful research topic [46–48]. Fortunately, the above consistent tanh expansion (CTE) method can be easily applied to investigate interaction solutions between a soliton and any other types of nonlinear excitations. Next, we construct the interaction solutions to the HBK system on the basis of CTE solvable Theorem 7.According to Theorem7, we find that the explicit solutions to HBK system can be obtained by solving the w and u0 equation (19), (20). In order to solve the w equation (19), if we choose any fixed solution to the STO equation (20), we will get the solution of w by solving the variable coefficient PSTO equation (19).Now, we first restrict the simplest solutionu0=c to equation (20) and substitute it into equation (19), which will be simplified in the following constant coefficient PSTO equation
(24)wt+4wxx+3wx2+3cwxx+16wx3+12cwxc+2wx=0.From equation (24), we know that w has the following trivial solution,
(25)w=kx+qt+z0,andq=−43c2k+6ck2+4k3,wherez0 is an arbitrary constant, and the corresponding single soliton solution to the HBK system is as follows
(26)u=ktanhkx+qt+z0+k+c,v=k2sech2kx+qt+z0.In order to study the properties of the solution, we plot the structure of the solution. Figures1(a) and 1(b) display this kind of single soliton solution for the fields u and v.Figure 1
Plot (a) refers to the single kink solutions ofu, and (b) refers to bell solitary wave of v in equation (26), with parameters: {c=4/5,k=−4/17,z0=2}.
(a)(b)In what follows, we are interesting in finding another kind of solution in the form(27)w=kx+qt+p2,which reflects the interaction solution between a solition and an STO/PSTO wave. Substituting (27) into equation (24) will lead to
(28)pt+4px3+432px2+3c1px+3c12p+pxxx+12c1px2=0,with c1=c+2k. The corresponding solution to the HBK system (1) is as follows:
(29)u=122k+px1+tanhw+c,v=k+12px2sech2w+12pxx1+tanhw.Now, employing the known solutions to (28), such as resonant soliton solutions and multisoliton solutions from literature [49–52], we will construct the interaction solutions between a soliton and STO/PSTO wave.


### 4.1. The Resonant Multisoliton Solutions

It can be verified that equation (28) possesses the following multiple wave solution
(30)p=ln∑i=1nciexpkix−4ki3+3kic12+3c1ki2t+zi.Substituting (27) and (30) into (29), the n+1 resonant multisoliton solutions to the HBK system can be directly obtained, which displays soliton fission and fusion. For example, if we take n=2 in (30), with parameters: c1=−1/2,g1=1/15,k1=1,z1=1,g2=1/12,k2=−1/13,z2=2 and k>0, then the solution (29) shows the two solitary wave fusion. It shows that two solitons fuse into one soliton after collide, which reflects the inelastic collision between solitons and is displayed in Figure 2.Figure 2
Plot (a) and (b) to the resonant soliton solution ofu and v with k>0. (a) is a graph of u, which shows that two kink solitary waves collide and fuse into a resonant soliton solution at t=0. (b) is a graph of v, which shows that two bell solitary waves collide and fuse into a resonant soliton solution at t=0.
(a)(b)And when thek<0, the solution (29) shows the soliton fission, which also reflects the inelastic collision between solitons and is displayed in Figure 3.Figure 3
Plot (a) and (b) to the resonant soliton solution ofu and v with k<0. (a) is the two kink solitary wave fission solution for u, and (b) is the solitary wave fission for v.
(a)(b)If we choosen=3 in (30), as Figure 4, we derive the three kink solitary wave fusion with parameters: c1=−2/3,g1=2,g2=10,g3=3/5,k1=1.2,k2=5/6,k3=1/3,z1=3.8,z2=3/4,z3=3.Figure 4
Plot (a) and (b) to the resonant soliton solution ofu and v. (a) is a graph of u, which shows that three kink solitary waves collide and fuse into a resonant soliton solution at t=0. (b) is a graph of v, which shows that three bell solitary waves collide and fuse into a resonant soliton solution at t=0.
(a)(b)


### 4.2. Soliton Interactions with Periodic Waves

#### 4.2.1. Soliton Interactions with Sine-Cosine Periodic Waves

It is not difficult to verify that equation (28) possesses the following explicit solutions
(31)p=ln∑i=1naicoski−x+43li2−ki2+3c1+li2texpl1x−4tli3+3lic12−ki2−c1ki2−2li2,whereai,ki,li denote arbitrary constants. Substituting the equation (31) into (29), we can obtain soliton interactions with sine-cosine periodic waves.


#### 4.2.2. Soliton Interactions with Cnoidal Periodic Waves

Interaction solution between soliton and cnoidal periodic waves can display many more interesting physical phenomena, and it can be difficult to obtained by other methods. Here, we takep as the following special Jacobian elliptic function
(32)p=l1x+l2t+μEπsnk1x+k2t,m,n,m,l1,l2,k1,k2,μ, and m are determined later, where snz,m is the usual Jacobian elliptic sine function and
(33)Eπζ,n,m=∫0ζdt1−nt21−t21−m2t2is the third type of incomplete elliptic integral. Substituting (32) into (28) and solving the overdetermined equations with the help of maple will come the following result
(34)k2=−1μ4k13μ3+12c1k12μ2+12k12l1μ2+12c12k1μ+24c1k1l1μ+12k1l12μ+12c12l1+12c1l12+4l13+l2,withn=0, and l1,l2,k1,m,μ are arbitrary constants. Substituting equation (32) and (34) into (29), we can obtain soliton interactions with cnoidal periodic waves. It is concluded that the CTE solutions (29) present the interaction solutions to multiple solitons and multiple periodic waves by substituting (31), (32), and (34). Here, we plot the structure of the solution, and Figures 5(a)–5(d) display the two kinds of interactions solution for the field u.Figure 5
(a) The interaction solution between soliton and sine-cosine wave ofu in equation (29) by substituting (27) and (31), with parameters: {a1=1/3,a2=1/3,c1=−6,k1=1/4,k2=1/4,l1=3/20,l2=3/20}. (b) The interaction solution between soliton-cnoidal wave of u in equation (29) by substituting (27) and (32), with parameters: {μ=50,k1=1/4,l1=−1,l2=0.37,c1=2/3,m=0.7}. (c) and (d) are the overhead view of the wave (a) and (b).This interaction solution reflects the elastic collision.
(a)(b)(c)(d)In a word, according to equations (19) and (20), taking any solution u0 to the STO equation and the variable coefficient PSTO system, one can get the corresponding interaction solutions to the HBK system (1) via the theorem 7.For the non-CTE case, we can still find some special exact solutions for the overdetermined system (23), such as the system (23) possessing the following exact solution
(35)w1=12lns1x+s2−4s1x+s2+4+s3,(36)w2=−12ln2+14lntan2/82x+t+1+12tan2/82x+t+1+1−822.Substituting the expressions (35) and (36) into (22) yields the special rational solutions to a single soliton solution, respectively.


## 4.1. The Resonant Multisoliton Solutions

It can be verified that equation (28) possesses the following multiple wave solution
(30)p=ln∑i=1nciexpkix−4ki3+3kic12+3c1ki2t+zi.Substituting (27) and (30) into (29), the n+1 resonant multisoliton solutions to the HBK system can be directly obtained, which displays soliton fission and fusion. For example, if we take n=2 in (30), with parameters: c1=−1/2,g1=1/15,k1=1,z1=1,g2=1/12,k2=−1/13,z2=2 and k>0, then the solution (29) shows the two solitary wave fusion. It shows that two solitons fuse into one soliton after collide, which reflects the inelastic collision between solitons and is displayed in Figure 2.Figure 2
Plot (a) and (b) to the resonant soliton solution ofu and v with k>0. (a) is a graph of u, which shows that two kink solitary waves collide and fuse into a resonant soliton solution at t=0. (b) is a graph of v, which shows that two bell solitary waves collide and fuse into a resonant soliton solution at t=0.
(a)(b)And when thek<0, the solution (29) shows the soliton fission, which also reflects the inelastic collision between solitons and is displayed in Figure 3.Figure 3
Plot (a) and (b) to the resonant soliton solution ofu and v with k<0. (a) is the two kink solitary wave fission solution for u, and (b) is the solitary wave fission for v.
(a)(b)If we choosen=3 in (30), as Figure 4, we derive the three kink solitary wave fusion with parameters: c1=−2/3,g1=2,g2=10,g3=3/5,k1=1.2,k2=5/6,k3=1/3,z1=3.8,z2=3/4,z3=3.Figure 4
Plot (a) and (b) to the resonant soliton solution ofu and v. (a) is a graph of u, which shows that three kink solitary waves collide and fuse into a resonant soliton solution at t=0. (b) is a graph of v, which shows that three bell solitary waves collide and fuse into a resonant soliton solution at t=0.
(a)(b)


## 4.2. Soliton Interactions with Periodic Waves

### 4.2.1. Soliton Interactions with Sine-Cosine Periodic Waves

It is not difficult to verify that equation (28) possesses the following explicit solutions
(31)p=ln∑i=1naicoski−x+43li2−ki2+3c1+li2texpl1x−4tli3+3lic12−ki2−c1ki2−2li2,whereai,ki,li denote arbitrary constants. Substituting the equation (31) into (29), we can obtain soliton interactions with sine-cosine periodic waves.


### 4.2.2. Soliton Interactions with Cnoidal Periodic Waves

Interaction solution between soliton and cnoidal periodic waves can display many more interesting physical phenomena, and it can be difficult to obtained by other methods. Here, we takep as the following special Jacobian elliptic function
(32)p=l1x+l2t+μEπsnk1x+k2t,m,n,m,l1,l2,k1,k2,μ, and m are determined later, where snz,m is the usual Jacobian elliptic sine function and
(33)Eπζ,n,m=∫0ζdt1−nt21−t21−m2t2is the third type of incomplete elliptic integral. Substituting (32) into (28) and solving the overdetermined equations with the help of maple will come the following result
(34)k2=−1μ4k13μ3+12c1k12μ2+12k12l1μ2+12c12k1μ+24c1k1l1μ+12k1l12μ+12c12l1+12c1l12+4l13+l2,withn=0, and l1,l2,k1,m,μ are arbitrary constants. Substituting equation (32) and (34) into (29), we can obtain soliton interactions with cnoidal periodic waves. It is concluded that the CTE solutions (29) present the interaction solutions to multiple solitons and multiple periodic waves by substituting (31), (32), and (34). Here, we plot the structure of the solution, and Figures 5(a)–5(d) display the two kinds of interactions solution for the field u.Figure 5
(a) The interaction solution between soliton and sine-cosine wave ofu in equation (29) by substituting (27) and (31), with parameters: {a1=1/3,a2=1/3,c1=−6,k1=1/4,k2=1/4,l1=3/20,l2=3/20}. (b) The interaction solution between soliton-cnoidal wave of u in equation (29) by substituting (27) and (32), with parameters: {μ=50,k1=1/4,l1=−1,l2=0.37,c1=2/3,m=0.7}. (c) and (d) are the overhead view of the wave (a) and (b).This interaction solution reflects the elastic collision.
(a)(b)(c)(d)In a word, according to equations (19) and (20), taking any solution u0 to the STO equation and the variable coefficient PSTO system, one can get the corresponding interaction solutions to the HBK system (1) via the theorem 7.For the non-CTE case, we can still find some special exact solutions for the overdetermined system (23), such as the system (23) possessing the following exact solution
(35)w1=12lns1x+s2−4s1x+s2+4+s3,(36)w2=−12ln2+14lntan2/82x+t+1+12tan2/82x+t+1+1−822.Substituting the expressions (35) and (36) into (22) yields the special rational solutions to a single soliton solution, respectively.


## 4.2.1. Soliton Interactions with Sine-Cosine Periodic Waves

It is not difficult to verify that equation (28) possesses the following explicit solutions
(31)p=ln∑i=1naicoski−x+43li2−ki2+3c1+li2texpl1x−4tli3+3lic12−ki2−c1ki2−2li2,whereai,ki,li denote arbitrary constants. Substituting the equation (31) into (29), we can obtain soliton interactions with sine-cosine periodic waves.


## 4.2.2. Soliton Interactions with Cnoidal Periodic Waves

Interaction solution between soliton and cnoidal periodic waves can display many more interesting physical phenomena, and it can be difficult to obtained by other methods. Here, we takep as the following special Jacobian elliptic function
(32)p=l1x+l2t+μEπsnk1x+k2t,m,n,m,l1,l2,k1,k2,μ, and m are determined later, where snz,m is the usual Jacobian elliptic sine function and
(33)Eπζ,n,m=∫0ζdt1−nt21−t21−m2t2is the third type of incomplete elliptic integral. Substituting (32) into (28) and solving the overdetermined equations with the help of maple will come the following result
(34)k2=−1μ4k13μ3+12c1k12μ2+12k12l1μ2+12c12k1μ+24c1k1l1μ+12k1l12μ+12c12l1+12c1l12+4l13+l2,withn=0, and l1,l2,k1,m,μ are arbitrary constants. Substituting equation (32) and (34) into (29), we can obtain soliton interactions with cnoidal periodic waves. It is concluded that the CTE solutions (29) present the interaction solutions to multiple solitons and multiple periodic waves by substituting (31), (32), and (34). Here, we plot the structure of the solution, and Figures 5(a)–5(d) display the two kinds of interactions solution for the field u.Figure 5
(a) The interaction solution between soliton and sine-cosine wave ofu in equation (29) by substituting (27) and (31), with parameters: {a1=1/3,a2=1/3,c1=−6,k1=1/4,k2=1/4,l1=3/20,l2=3/20}. (b) The interaction solution between soliton-cnoidal wave of u in equation (29) by substituting (27) and (32), with parameters: {μ=50,k1=1/4,l1=−1,l2=0.37,c1=2/3,m=0.7}. (c) and (d) are the overhead view of the wave (a) and (b).This interaction solution reflects the elastic collision.
(a)(b)(c)(d)In a word, according to equations (19) and (20), taking any solution u0 to the STO equation and the variable coefficient PSTO system, one can get the corresponding interaction solutions to the HBK system (1) via the theorem 7.For the non-CTE case, we can still find some special exact solutions for the overdetermined system (23), such as the system (23) possessing the following exact solution
(35)w1=12lns1x+s2−4s1x+s2+4+s3,(36)w2=−12ln2+14lntan2/82x+t+1+12tan2/82x+t+1+1−822.Substituting the expressions (35) and (36) into (22) yields the special rational solutions to a single soliton solution, respectively.


## 5. Summary and Discussion

It is summarized that the residual symmetry of the higher-order Broer-Kaup (HBK) system is obtained from the standard truncated Painlevé expansion. For the original nonlinear system, the residual symmetries are nonlocal. However, it can be localized to Lie point symmetries by introducing a properly auxiliary dependent variable, and the finite transformation of the residual symmetry is obtained by solving the standard Lie’s initial value problem, and the group invariant solutions of HBK are derived. Furthermore, in order to study other types of solutions for the HBK system, we did not employ the general symmetry reduction method, but use the CTE method, by discussing the consistent tanh expansion solvable of the HBK system, we derive abundant exact solutions to the HBK system, and these solutions reveal the interaction among solitons and other types of nonlinear excitations especially the STO/PSTO wave such as the multiple resonant solitons and periodic waves. These kinds of solution can be easily applicable to the analysis of physically interesting processes.The method presented in this paper would be applied to other various interesting integrable models. One can also consider the relationship between residual symmetry and other nonlocal symmetries. And more research should be done on the CTE method and more types of the interaction solutions.



---

*Source: 9975303-2021-05-19.xml*