Dynamics of exact soliton solutions to the coupled nonlinear system using reliable analytical mathematical approaches

Muhammad Bilal,Usman Younas and Jingli Ren

Henan Academy of Big Data/School of Mathematics and Statistics,Zhengzhou University,Zhengzhou 450001,China

Abstract Nonlinear Schr?dinger-type equations are important models that have emerged from a wide variety of fields,such as fluids,nonlinear optics,the theory of deep-water waves,plasma physics,and so on.In this work,we obtain different soliton solutions to coupled nonlinear Schr?dinger-type(CNLST)equations by applying three integration tools known as the-expansion function method,the modified direct algebraic method(MDAM),and the generalized Kudryashov method.The soliton and other solutions obtained by these methods can be categorized as single(dark,singular),complex,and combined soliton solutions,as well as hyperbolic,plane wave,and trigonometric solutions with arbitrary parameters.The spectrum of the solitons is enumerated along with their existence criteria.Moreover,2D,3D,and contour profiles of the reported results are also plotted by choosing suitable values of the parameters involved,which makes it easier for researchers to comprehend the physical phenomena of the governing equation.The solutions acquired demonstrate that the proposed techniques are efficient,valuable,and straightforward when constructing new solutions for various types of nonlinear partial differential equation that have important applications in applied sciences and engineering.All the reported solutions are verified by substitution back into the original equation through the software package Mathematica.

Keywords:soliton solutions,exact solutions,CNLST equations,-expansion function method,MDAM,generalized Kudryashov method

1.Introduction

For many years,investigating the exact solutions of nonlinear partial differential equations(NLPDEs)has turned out to be a charming and challenging area of research for mathematicians and research communities,because they play a broad and significant role in the study of nonlinear physical phenomena in mathematical studies and applied physics,with essential applications in several areas of engineering and natural science including fluid mechanics,chemistry,thermodynamics,physics,electromagnetism,biomathematics,mathematical physics,and so on.Scholars have focused on exact or analytical solutions,due to their essential contribution to the analysis of the real features of nonlinear problems.Due to their wide utilization and applications in the domain of nonlinear sciences,interest in the study of NLPDEs has been increasing[1–9].

Moreover,solitons are also famous as a particular type of solitary wave,which are the solutions to various kinds of NLPDE.This incredible type of solitary waves has various and substantial use in different fields because of its specific characteristic of stability.In short,waves(mainly dispersive in nature)inelastically scatter solitary waves and lose energy due to radiation phenomena;as a result of the solitary wave’s disappearance,the dispersive waves hold their shape and speed after a fully nonlinear connection.Soliton theory has made a significant contribution to the narration and expression of the physical behavior and meaning of nonlinear phenomena.Soliton theory has involved researchers in exploratory investigations due to its use in such diverse fields as media transmission,design,numerical materials science,mathematical physics,and different parts of nonlinear science[10–17].Therefore,it has recently become more appealing for the research community to acquire exact solutions by the use of capable computational packages that alleviate the complex and tedious algebraic computations.In the literature,various powerful computational techniques have been designed to describe the natures of the diverse forms of the solutions[18–21]that have been established for nonlinear physical models.

To the best of our knowledge,it is known that in the existing literature,the solutions of the governing model have not been investigated using the-expansion function method[22],the modified direct algebraic method(MDAM)[23],and the generalized Kudryashov method[24].Motivated by this,we employed these three mathematical techniques to obtain different forms of the solution.The proposed methodologies are powerful,reliable,capable of examining NLPDEs,consistent with computer algebra,and yield more general solutions.The discovered constructed solutions are novel and have potential applications in the nonlinear sciences.

The structure of the rest of this paper is organized as follows.In section 2,the governing equation is presented.In section 3,we discuss the application of the proposed methods.In section 4,the results and a discussion are presented.Finally,the conclusions are revealed in section 5.

2.Governing equation

This paper is concerned with the CNLST system,as presented by Ma and Geng[25],

where ?≠±1 is a real constant,φ and ψ represent real functions of the spatial variable x and the temporal variable t,respectively,while θ is a complex function.

In the following section,the applications of three methods are discussed.

3.Mathematical preliminaries

To solve the above equation(1),we utilize the traveling-wave transformations θ(x,t)=U(η)eiΦ,φ(x,t)=H(η),and ψ(x,t)=Q(η),where Φ=ρx-?t and η=x-c t+ζ0.Here ρ,?,c,and ζ0represent the wave number,frequency,velocity,and phase constant,respectively.Imposing these transformations onto equation(1)yields

By substituting equations(3)–(5)into equation(2),we obtain the following ODE,

Using the homogeneous balance rule that equation(6)yields,n=1.The solution of equation(6)reduces to:

ForSet 1,

Case-1(Trigonometric solutions:

If?Ω>0,

Case-2(Hyperbolic solution:

If?Ω<0,

For the soliton solution,taking Ξ1=Ξ2,we get a singular wave solution as follows:

ForSet 2,

Case-1(Trigonometric solutions:

If?Ω>0,

Case-2(Hyperbolic solution:

If?Ω<0,

For the soliton solution,if we take Ξ1=Ξ2,we get a dark soliton solution,as follows:

ForSet 3,

Case-1(Trigonometric solutions:

If?Ω>0,

For the soliton solution,if we take Ξ1=Ξ2,we get a periodic wave solution,as follows:

Case-2(Hyperbolic solution:

If?Ω<0,

For the soliton solution,if we take Ξ1=Ξ2,we get a combined dark–singular wave solution,as follows:

Case-3(Rational solutions:

If?=0,Ω≠0,

If we take Ξ1=Ξ2,a plane-wave solution is obtained:

ForSet 4,

Case-1(Trigonometric solutions:

If?Ω>0,

For the soliton solution,if we take Ξ1=Ξ2,we get a periodic wave solution,as follows:

Case-2(Hyperbolic solution:

If?Ω<0,

For the soliton solution,if we take Ξ1=Ξ2,we obtain a combined dark–singular wave solution,as follows:

For all the above sets,η=-c t+ζ0+x.

3.2.MDAM

By utilizing the homogeneous balance rule on equation(6),we obtain n=1,which implies the solution of equation(6)is as follows:

where a0,a1,and b1are parameters.In order to find the parameters involved in equation(56),we substitute equation(56)along with(Z′=χ+Z2)into equation(6),and we get a cluster of equations on equating the same power coefficients of of Z.Furthermore,by using Mathematica,we obtain a cluster of solutions,as follows:

For Set 1

Case-1:

When χ<0,we get the following solutions.

Dark soliton structure:

Singular wave structure:

Case-2:

When χ>0,then the following periodic solutions of different forms are obtained:

and

For Set 2

Case-1:

When χ<0,we get the singular and dark soliton structures,respectively:

and

Case-2:

When χ>0,then the periodic solutions are:

and

For Set 3

Case-1:

When χ<0,we get the following mixed hyperbolic solution:

Case-2:

When χ>0,then the periodic solution is expressed as:

For all the above sets,η=-c t+ζ0+x.

3.3.Generalized Kudryashov method

In this section,we apply the above method to obtain the solutions of the governing model.To begin with,we take account of the fact that the homogeneous balance between U3and U″gives the relation T=H+1.In particular,for H=1,we get T=2.

Therefore,equation(6)along with(G′(η)=G2(η)-G(η))take the following form of the solution:

where a0,a1,a2,b0and b1are to be determined.Also,

Now,solving equations(6)and(87),the following system of equations is obtained:

With the assistance of computational software such as Mathematica,we solve system 89,and a variety of solution sets is obtained,as follows:

Set 1

On substituting the above values of parameters into equation(87)and setting S=1 in equation(88),we secure the soliton solution to equation(1)in terms of hyperbolic functions.

Set 2

For S=1,and solving equations(87)and(88)together,we get a soliton solution,as follows:

Set 3

For S=1,and solving equations(87)and(88)together,we get a singular soliton solution as follows:

Set 4

For S=1,and solving equations(87)and(88)together,we get the dark soliton solution as follows:

Set 5

For S=1,and solving equations(87)and(88)together,the plane-wave solution can be expressed as:

Set 6

For S=1,and solving equations(87)and(88)together,we obtain a singular wave solution as follows:

For all the above sets,η=-c t+ζ0+x.

4.Results and discussion

In this section,a comparison is made between our outcomes and some existing results in published works.In this work,we utilized an analytical mathematical method to solve the CNLST equations.A variety of solutions were extracted,such as single(dark,singular),complex,and combined solitons as well as hyperbolic,plane wave,and trigonometric solutions.In published works,we have observed that Elboree[26]only discussed one solution for the CNLST system using He’s semi-inverse variational principle,whereas Ma et al[25]only presented two solutions through the Darboux transformation.Yong et al[27]introduced eight solutions for the CNLST system using the truncated Painleve expansion and the direct quadrature method.Abdelrahman et al[28]discussed exact soliton solutions and employed three methods.If we compare these methods and the three proposed methods in this paper,our methods are more effective at providing many solutions than those other methods.Consequently,these methods are proficient,adequate,and robust for handling other NLPDEs in mathematical physics and applied mathematics.We observe that the results presented in this study could be helpful in explaining the physical meaning of various NLPDEs arising in the different fields of nonlinear sciences.To aid clear and good understanding,the absolute physical behaviors of some of the reported solutions are exhibited through 3D,2D and contour graphs using suitable parameter values.To visualize the dynamics of trigonometric,dark,dark–singular,dark,and complex combined soliton and plane wave solutions which appear in equations(20),(29),(53),(57),(69),(84),(90),and(105),we present figures 1–8,respectively.From the physical description of some solutions and our discussion of the results,we conclude that our modified mathematical methods presented here are fruitful tools for investigating further results of nonlinear wave problems in applied science.

Figure 1.The plots of equation(20)are presented in 3D,2D and contour wave profiles,respectively.

Figure 2.The plots of equation(29)are presented in 3D,2D and contour wave profiles,respectively.

Figure 3.The plots of equation(53)are presented in 3D,2D and contour wave profiles,respectively.

Figure 4.The plots of equation(57)are presented in 3D,2D and contour wave profiles,respectively.

Figure 5.The plots of equation(69)are presented in 3D,2D and contour wave profiles,respectively.

Figure 6.The plots of equation(84)are presented in 3D,2D and contour wave profiles,respectively.

Figure 7.The plots of equation(90)are presented in 3D,2D and contour wave profiles,respectively.

Figure 8.The plots of equation(105)are presented in 3D,2D and contour wave profiles,respectively.

5.Conclusion

In this study,three innovative integration norms have been effectively implemented for the CNLST equations that have emerged in several fields of applied sciences.Different types of solution,such as single(dark,singular),complex,and combined solitons as well as hyperbolic,plane-wave,and trigonometric solutions were successfully obtained.All the solutions obtained were verified by substituting them back into the original equation via the software package Mathematica.The graphical representations of some solutions were also depicted using 2D,3D and contour profiles that displayed the physical behaviors of some of the solutions obtained.These diverse solutions demonstrate the power,capability,consistency,and effectiveness of these methods and can be used to solve many other NLPDEs in mathematical physics.We observe that the results presented in this study could be helpful in explaining the physical meanings of various nonlinear evolution equations that have arisen in the different fields of nonlinear sciences.For instance,hyperbolic functions appears in different areas such as the calculation of special relativity,the Langevin function for magnetic polarization,the gravitational potential of a cylinder,and the calculation of the Roche limit in the profile of a laminar jet[29].The wave results obtained confirm the value of this research and also have significant applications in mathematics and physics.Hence,our techniques,fortified by symbolic computation,provide an active and potent mathematical implement for solving diverse interesting nonlinear wave problems.

Acknowledgments

The authors would like to acknowledge the financial support provided for this research via the National Natural Science Foundation of China(11771407-52071298),ZhongYuan Science and Technology Innovation Leadership Program(214200510010),and the MOST Innovation Methodproject(2019IM050400).They also thank the reviewers for their valuable reviews and kind suggestions.

Conflict of interest

The authors declare that there are no conflicts of interest.