Solitons and breathers for nonisospectral mKdV equation with Darboux transformation

Solitons and breathers for nonisospectral mKdV equation with Darboux transformation

Abstract

Under investigation in this paper is the nonisospectral and variable coefficients modified Kortweg-de Vries (vc-mKdV) equation, which manifests in diverse areas of physics such as fluid dynamics, ion acoustic solitons and plasma mechanics. With the degrees of restriction reduced, a simplified constraint is introduced, under which the vc-mKdV equation is an integrable system and the spectral flow is time-varying. The Darboux transformation for such equation is constructed, which gives rise to the generation of variable kinds of solutions including the double-breather coherent structure, periodical soliton-breather and localized solitons and breathers. In addition, the effect of variable coefficients and initial phases is discussed in terms of the soliton amplitude, polarity, velocity and width, which might provide feasible soliton management with certain conditions taken into account.

PACS numbers: 05.45.Yv, 02.30.Ik, 47.35.Fg

Keywords: Nonisospectral modified Korteweg-de Vries equation; Darboux transformation; Breather; Soliton management

I. Introduction

The modified Kortweg-de Vries (mKdV) type equations arise in diverse areas of physics, including fluid dynamics, ion acoustic solitons and plasma mechanics [[1], [3], [5], [4], [2]]. For instance, the dynamics of the interfacial waves in a two-layer fluid of slowly varying depth is studied by the formulated standard mKdV equation [[3]], the ion-acoustic solitary waves in certain unmagnetized plasma is described by the cylindrical and spherical mKdV equation [[4]] and the propagation of circularly polarized few-cycle pulses with wave polarization is investigated by a non-integrable complex mKdV equation [[5]]. With the inhomogeneous environmental density and boundary conditions taken into account [[6]], the constant coefficients mKdV equation can be extended as followings:

(1)

where a(t), b(t), and d(t) are functions of variable t.

With different selections of coefficients, Eq. (1) has been investigated for physical interest [[7], [8], [9], [10]]. Thereinto, periodic wave solutions are constructed in bilinear forms based on a multidimensional Riemann theta function [[7]], group classification is carried out and all the classes under consideration are normalized [[8]], the pulse waves in thin walled prestressed elastic blood vessels is described [[9]] and a modified perturbation technique is applied to the problem of the structural instability of algebraic solitons [[10]].

For constructing soliton solutions from the soliton equations like Eq. (1), there exist many methods such as Hirota bilinear method [[11], [12], [13]], associated scattering problem [[14]] and Darboux transformation (DT) [[15], [16], [17], [18], [19]], in which the DT has proven itself a purely algebraic iterative tool [[15], [16]]. The main difficulty in finding a proper spectral problem of constructing DT has been investigated and many spectral problems have been searched for, including different forms of DT of Eq. (1) under certain conditions [[17], [16], [18], [19]].

Except for soliton solutions, Eq. (1) also has solutions in the form of oscillating packets (breathers) [[10]], which along with the solitons, determine the asymptotics of the wave field [[20]]. It is believed that some temporal and spatial variability that has been observed in oceanic internal soliton fields may be due to the breather interaction [[21]] and several investigations have been carried out [[21], [22], [23], [20]]. The circumstances supporting the formation of breathers are determined with various piecewise-constant initial conditions [[22]], some detailed examinations of the breather-soliton interaction process are analyzed by the Hirota bilinear method [[21]] and the breather in numerical simulations using the full nonlinear Euler equations for stratified fluid is presented under several forms of the initial disturbance [[23]].

Meanwhile, Eq. (1) becomes integrable when the variable coefficients satisfy certain constraint conditions, in which and , leading to analytical discussion related to different branches of physics [[24], [25], [11], [26], [27], [28], [29]]. Different from above investigation, hereby we introduce a constraint reducing the degrees of restriction

(2)

under which Eq. (1) is an integrable system and there can be an integrable constant which is normalized to be 1 in this case.

However, to our best knowledge, the Lax pair and DT for Eq. (1) under constraint (2) have not been obtained. Therefore, in this paper, we will construct the DT with spectral problem for Eq. (1) and generate multi-soliton and breather solutions from the seed ones.

The outline of the present paper is as follows. In Section II and III, the Lax pair and DT for Eq. (1) under constraint (2) are constructed respectively. In Section IV, one-soliton solution is generated by a DT and the characteristic line with soliton amplitude is given. In Section V, two-soliton and three-soliton solutions are investigated analytically by appropriate selection of variable coefficients. In Section VI, different kinds of breathers as well as the interaction between breather and soliton are thrived, including the double-breather coherent structure, periodical soliton-breather and localized breathers. In Section VII, the conclusions are given.

II. Lax pair

Based on AKNS procedure [[18]], the Lax pair of Eq. (1) under the general constraint (2) is constructed with the undetermined coefficients corresponding to time and space

(3)
(4)

where with T representing the transpose of the vector. The should comply with the compatibility condition, which leads to a zero curvature equations [[15]]

(5)

Substituting Eqs. (3) and (4) into the zero curvature equations and comparing the coefficients of the power of , we have

(6)
(7)
(8)
(9)

Meanwhile, spectral parameter should satisfy the following relation

(10)

III. Construction of DT

Now we introduce a gauge transformation for the spectral problem

(11)

where is defined by

(12)
(13)

The gauge transformation is called DT [[17]], if the spectral problem (3) and (4) can be transformed into

(14)
(15)

where and has the same form as U and V and the old potential is mapped into a new potential .

Suppose

(16)

where .

Substituting Expression (16) into Eqs. (14) and (15), by a direct calculation, we have

(17)
(18)
(19)

where .

The transformation between old potentials and new ones is given as below

(20)

So far, we have obtained the DT (16)-(20), which transforms the matrix spectral problem into another spectral problem of the same type and can generate new solutions from seed ones by purely algebraic iteration.

It is worth noting with certain modification, the DT can be extended as following

(21)
(22)
(23)

where and are linearly independent solutions of spectral problems.

IV. One soliton solution

Employing as the original solution for Eq. (1), we can obtain the solutions of spectral problems (3) and (4) as below

(24)
(25)
(26)

where is an integration constant representing soliton initial phase, in every iteration which can change to a certain value.

For a given spectral parameters and initial phase , the explicit solutions can be obtained in terms of Expressions (24)-(26), which leads to a new one soliton solution with the aid of DT (20)

(27)

where

(28)
(29)
(30)

Meanwhile, the soliton amplitude and characteristic line can be respectively derived as

(31)
(32)

which indicates one soliton amplitude is only affected by when is given. It is worth noting that the polarity of soliton also refers to the sign of A. Meanwhile, the soliton velocity can be obtained by the derivation of Expression (32).

V. Multi-soliton solutions

For general application in various fields related to soliton dynamics, we need to consider two soliton interaction, which can be generated by a second DT based on the one-soliton solution. For soliton characters describing their physical features, soliton polarity (or phase) should be taken into account. It is worth mentioning that two soliton polarities are allowed in the mKdV framework due to the isotropic nonliearity [[20]] and the bipolar soliton interaction plays a role in obtaining soliton management. Therefore, with a positive coefficient of the cubic nonlinear term in Eq. (1), the effect of initial phase on the soliton polarities will be discussed in the following section.

()                                                                                    ()
Fig.  1: Solitonic propagation and interaction with coefficients , , , , and initial phase: (a) , (b)

Compared with Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation(a), an initial phase shift value of takes place along with the soliton inverse polarity in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation(b). The depression, which has a negative amplitude, inverses its polarity while the elevation having the positive amplitude remains unchanged. As a result, the group changes from the bipolar solitons into unipolar solitons. The soliton cross disappears after the inverse polarity of the depression.

()                                                                                      ()
Fig.  2: Solitonic propagation and interaction with variable coefficients and spectral parameters , , , , , and line-damping coefficient: (a) , (b)

Under the analytical discussion before, soliton amplitude can be interpreted by the effect of line-damping coefficient when is given. As demonstrated in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation, a periodic value independent to time is applied for and correspondingly the soliton dynamics presents amplitude periodicity. Coefficient is capable of describing nonuniformity of media and should be evaluated in the nonisospectral problems. For positive and negative implication of , the solitons propagation dynamics differs in terms of its width. With being 1 and -1 in this case, the soliton compress and swell respectively along the propagation.

()                                                                                    ()

()                                                                                    ()

Fig.  3: Interaction among elevation and depression with variable coefficients and spectral parameters , , , , , , and , (a) , (b) . (c) and (d) profile (a) and (b) respectively when t=-1(solid line) and t=1(dash line)

Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation depicts three soliton solution, including the interaction among depression and elevation. Comparing Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation(a) with (b), in addition to change the value of initial phase, we could also manage to replace the elevation by depression via a sign reversal of spectral parameter from positive to negative. The three soliton propagation is exhibited by profiles in Figs. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation(c) and (d) at certain time, with high-amplitude elevation (depression) firstly leaving behind the low-amplitude elevation (depression) and then surpassing when time is approaching to 1. Such phenomena are also indicated in ref [[11]] with Eq. (1) under an extra ().

VI. Breathers and breather-soliton solutions

Breather solutions can be constructed by DT under a pair of complex conjugate spectral parameters (). Similar to the definition that breather can be regarded as a central valley with two small hills of elevation adjacently, which should reverse its polarities every half a cycle later [[21]], the breather generated by DT can be also regarded as valley-hill feature in symmetric form at certain time.

()                                                                                  ()
Fig.  4: (a) Breather propagation , , , , and ; (b) Profiles (a) when t=-0.566(solid line), t=-0.267(bold dashed line), t=0.155(bold solid line) and t=0.454(dashed line)

As illustrated in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation, the propagation of breathers demonstrates periodically pulsating and isolated wave forms by profiles at certain time. Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation(b) illustrates a downward displacement at center and the horizontal symmetry of breather structure. With line-damping coefficient taken into consideration, breather propagation presents a character of time-space locality.

()                                                                                      ()
Fig.  5: (a) Breather propagation with spectral parameters , , initial phases and coefficients , and ; (b) Contour plot of (a)

Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation shows a kind of periodic soliton solution for the existence of two periodic coefficient and supporting the periodically oscillating breather. From the analytical discussion, the character line is determined only by when the other coefficients and parameters are given, and the time dependent gives rise to the periodical soliton amplitude. In a conclusion, there exist three kinds of periodism with respect to breather, character line and the soliton amplitude.

()                                                                                      ()
Fig.  6: (a) Interaction between a breather and a soliton , , , , , and , (b) Contour plot of (a)

The contour plot in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation (b) exhibits the effects of the dispersive term and line-damping term on the interaction between a breather and a elevation. Being interacted by breather propagation, the elevation undergoes a similar periodicity corresponding to the breather and propagate without interrupting the pecks of breather.

()                                                                                      ()
Fig.  7: (a) Interaction between a breather and a soliton , , , , , , and ; (b) Contour plot of (a)

Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation depicts the collision between a soliton and breather with approximately x=0 being the intersection point, which witnesses the mutual interaction. The character line is overlapped and swell phenomena occur corresponding to the soliton and breather width. Compared with Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation, without periodical time-varying function , Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation demonstrates another kind of character line.

()                                                                                      ()
Fig.  8: (a) Interaction between two breather with spectral parameters , , , and and coefficients , and ; (b) Contour plot of (a)

With b(t) being periodical, we can observe two-breather periodical oscillation in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation. The continuous interaction between two breathers in one period is demonstrated in the contour plot, which indicates the two-breather structure in this case can be regarded as a coherent structure.

()                                                                                      ()
Fig.  9: (a) Interaction between a breather and two solitons with spectral parameters , , , and and coefficients , and ; (b) Contour plot of (a)

Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation illustrates the interaction between a breather and two solitons. It is worth noting that with the breather-soliton convergence at time approaching to zero, soliton amplitude also reaches its peak value due to the line-damping term . Compared with the structure including a soliton and a breather in Figs. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation and Solitons and breathers for nonisospectral mKdV equation with Darboux transformation, Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation presents three character lines, which individually support two solitons and a breather in its propagation.

VII. Conclusions

Based on the introduced constraint (2), under which Eq. (1) is integrable and becomes independent from , the work of our paper can be concluded as following:

(1) The Lax pair is generated by means of AKNS process with a nonisospectral flow under constraint (2). Meanwhile, DT is constructed based on the Lax pair, by which multi-soliton and breather solutions could be iterated from the seed ones.

(2) Analytical discussion corresponding to the characteristic line is applied, in which the one soliton amplitude (polarity), width (wave number) and velocity can be obtained. Meanwhile, with the aid of one soliton solution, multi-soliton solutions are generated with their graphical illustration. Initial phases are discussed by employing a value shift, which gives rise to the soliton inverse polarity. With appropriate selection of coefficients and spectral parameters, such reverse polarity might prevent soliton cross. In addition, variable coefficients influencing soliton amplitude, velocity and width are investigated and the interaction among three solitons is demonstrated.

(3) N-th iteration by DT can generate N-soliton solutions, including the breathers in the periodically pulsating and isolated wave forms, which are generated by employing a pair of complex conjugate spectral parameters () during iteration. For example, the four-soliton solution can generate two breathers with periodical oscillation as shown in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation or a breather and two solitons with three character lines as shown in Fig. Solitons and breathers for nonisospectral mKdV equation with Darboux transformation. The inherent periodism of breathers and the variable coefficients of Eq. (1) have coupling effects, which can yield abundant structures of breather and soliton solutions, such as a cyclical breather of three kinds of periodism, a breather-soliton interaction with time-space locality and two-breathers in a coherent structure. Such results can be extended to multi-soliton solutions by DT and similar phenomena could be observed.

Acknowledgements

We express our sincere thanks to all the members of our discussion group for their valuable comments. This work has been supported by the National Natural Science Foundation of China under Grant No. 11302014, and by the Fundamental Research Funds for the Central Universities under Grant Nos. 50100002013105026 and 50100002015105032 (Beijing University of Aeronautics and Astronautics).

References

  1. Z. T. Fu, S. D. Liu and S. K. Liu, Phys. Lett. A 326, 364, 2004.
  2. Y. L. Wang, Z. X. Zhou, X. Q. Jiang, X. D. Ni, Y. Zhang, J. Shen and P. Qian, Phys. Lett. A 373, 2944, 2009.
  3. K. Helfrich, W. K. Melville and J. W. Miles, J. Fluid. Mech. 149, 305, 1984.
  4. D. K. Ghosh, G. Mandal, P. Chatterjee and U. N. Ghosh, IEEE Trans. Plas. Sci. 41, 5, 2013.
  5. H. Leblond, H. Triki, F. Sanchez and D. Mihalache, Opt. Commun. 285, 356, 2012.
  6. R. Grimshaw, D. Pelinovsky, E. Pelinovsky and T. Talipova, Phys. D 159, 35, 2001.
  7. Y. Zhang, Z. L. Cheng and X. H. Hao, Chin. Phys. B 21, 12, 2012.
  8. O. Vaneeva, Commun. Nonlinear Sci. Numer. Simulat. 17, 611, 2014.
  9. X. L. Gai, Y. T. Gao, L. Wang, D. X. Meng, X. L, Z. Y. Sun and X. Yu, Commun. Nonlinear Sci. Numer. Simulat. 16, 1776, 2011.
  10. D. E. Pelinovsky and R. H. J. Grimshaw, Phys. Lett. A 229, 165, 1997.
  11. Z. Y. Sun, Y. T. Gao, Y. Liu and X. Yu, Phys. Rev. E 84, 026606, 2011.
  12. C. J. Wang, Z. D. Dai, S. Q. Lin and G. Mu, Appl. Math. Comput. 216, 341, 2010.
  13. Y. Zhang, J. B. Li and Y. N. Lv, Ann. Phys. 323, 3059, 2008.
  14. H. C. Hu and Q. P. Liu, Chaos Solitons Fractals 17, 921, 2013.
  15. Q. L. Zha, Commun. Nonlinear Sci. Numer. Simulat 57, 083506, 2016.
  16. Q. L. Zha and Z. B. Li, Chin. Phys. Lett. 25, 8, 2008.
  17. L. Wang, Y. T. Gao and F. H. Qi, Ann. Phys. 327, 1974, 2012.
  18. X. L. Gai, Y. T. Gao, D. X. Meng, L. Wang, Z. Y. Sun, X. L, Q. Feng, M. Z. Wang, X. Yu and S. H. Zhu, Commun. Theor. Phys. 53, 673, 2010.
  19. J. L. Ji and Z. N. Zhu, Commun. Nonlinear Sci. Numer. Simulat. 42, 699, 2017.
  20. A. V. Slyunyaev, J. Exp. Theor. Phys 92, 529, 2001.
  21. K. W. Chow, R. H. J. Grimshaw and E. Ding, Wave Motion 43, 158, 2005.
  22. S. Clarke, R. Grimshaw, P. Miller, E. Pelinovsky and T. Talipova, Chaos 10, 383, 2000.
  23. K. G. Lamb, O. Polukhina, T. Talipova, E. Pelinovsky, W. T. Xiao and A. Kurkin, Phys. Rev. E 75, 046306, 2007.
  24. W. L. Chan and K. S. Li, J. Phys. A 27, 883, 1994.
  25. W. L. Chan and X. Zhang, J. Phys. A 28, 407, 1995.
  26. C. Q. Dai, J. M. Zhu and J. F. Zhang, Chaos Solitons Fractals 27, 881, 2006.
  27. Q. Feng, Y. T. Gao, X. H. Meng, X. Yu, Z. Y. Sun, T. Xu and B. Tian, Int. J. Mod. Phys. B 25, 723, 2011.
  28. S. Zhang and T. C. Xia, Phys. Lett. A 372, 1741, 2008.
  29. Z. Y. Yan, Commun. Nonlinear Sci. Numer. Simulat. 4, 284, 1999.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minumum 40 characters
   
Add comment
Cancel
Loading ...
104562
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description