Next Article in Journal
Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses
Previous Article in Journal
Accelerated Subgradient Extragradient Algorithm for Solving Bilevel System of Equilibrium Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique

1
Department of Mathematics and Statistics, The University of Lahore, Lahore 54000, Pakistan
2
Department of Mathematics, The University of The Gambia, Serrekunda P.O. Box 3530, The Gambia
3
Department of Computer Science and Mathematics, Lebanese American University, Beirut P.O. Box 13-5053, Lebanon
4
Department of Mathematics, Art and Science Faculty, Siirt University, 56100 Siirt, Turkey
5
Department of Mathematics, Mathematics Research Center, Near East University, Near East Boulevard, Mersin 10, 99138 Nicosia, Turkey
6
Multidisciplinary Institute of Mathematics, Universitat Politènica de València, 46022 València, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(9), 1682; https://doi.org/10.3390/sym15091682
Submission received: 12 August 2023 / Revised: 28 August 2023 / Accepted: 31 August 2023 / Published: 1 September 2023
(This article belongs to the Section Mathematics)

Abstract

:
This paper investigates the ion-acoustic wave structures in fluid ions for the Benjamin–Bona–Mahony–Peregrine–Burgers (BBMPB) equation. The various types of wave structures are extracted including the three-wave hypothesis, breather wave, lump periodic, mixed-type wave, periodic cross-kink, cross-kink rational wave, M-shaped rational wave, M-shaped rational wave solution with one kink wave, and M-shaped rational wave with two kink wave solutions. The Hirota bilinear transformation is a powerful tool that allows us to accurately find solutions and predict the behaviour of these wave structures. Through our analysis, we gain a better understanding of the complex dynamics of ion-acoustic waves and their potential applications in various fields. Moreover, our findings contribute to the ongoing research in plasma physics that utilize ion-acoustic wave phenomena. To show the physical behaviour of the solutions, some 3D plots and their respective contour level are shown, choosing different values of the parameters.

1. Introduction

The study of wave phenomena is an important aspect of modern physics and has been a subject of interest for many years. Waves can be observed in various physical systems such as fluids, plasmas, and solids, and their dynamics can be described by a variety of mathematical models. In particular, the behaviour of waves in plasmas has been an active area of research due to its relevance for a wide range of applications, including space physics, fusion research, and plasma processing.
One of the most important wave phenomena in plasmas is the ion-acoustic wave [1], which is a compressional wave that is primarily driven by the motion of the ions in the plasma. The ion-acoustic wave is a fundamental mode of oscillation in plasmas and is characterized by its dispersion relation, which relates the wave frequency to the wave number. The study of ion-acoustic waves is of great importance in plasma physics because it provides insights into the basic plasma processes, such as energy transport, wave–particle interactions, and turbulence [2,3,4,5].
The BBMPB model is a nonlinear partial differential equation that can be used to explore nonlinear wave patterns. The equation is given by:
v t v x x t α v x x + λ v x + θ v v x + β v x x x = 0 ,
where v ( x , t ) is the unknown function, and α , β , λ , and θ are constants.
The BBMPB equation has been studied extensively by researchers in the field of nonlinear wave theory. In particular, there have been many studies on the existence and stability of solitary wave solutions for this equation. One of the earliest studies on the BBMPB equation was conducted by Benjamin and Peregrine [6], where they derived the BBMPB equation as a model for long waves in shallow water. They showed that the equation has solitary wave solutions that are stable under certain conditions. As example of some of the recent studies of this equation, Yang et al. constructed the travelling wave solutions for the Zakhrov–Kuznetsov–Benjamin–Bona–Mahony equation [7], while Akcagil et al. found the exact travelling wave solutions of nonlinear pseudoparabolic equations by using the G G expansion method [8]. Overall, the BBMPB equation has been studied extensively in the literature, and there have been many important results regarding its properties and solutions.
There are many techniques to obtain a different soliton solution for physical systems. Zou used the Riemann–Hilbert approach and gained the soliton interactions and position shift for the higher-order Gerdjikov–Ivanov equation [9,10], and Shen et al. considered the nonlocal nonlinear Schrödinger equation to obtain the Gaussian-like, nearly flat-topped, multipeak, and four-peak form solitons [11] and complex-valued astigmatic cosine–Gaussian soliton solutions [12]. Li and Guo explored optical solitons in the form of breathers, rogue waves, and semirational solutions on periodic backgrounds for the coupled Lakshmanan–Porsezian–Daniel equations [13], and Song et al. studied Laguerre–Gaussian and Hermite–Gaussian solitons in the nonlocal nonlinear Schrödinger equation [14]. Zhang and Xu worked on the localized symmetric and asymmetric solitary wave solutions using the Darboux transformation [15]. There are many mathematical techniques to explore the soliton solution such as the direct algebraic method [16], the sine-Gordon expansion method [17], the new MEDA method [18,19], the Riccati equation mapping (REM) method [20], the Sardar subequation method [21], the Jacobi elliptic functions method [22]. However, in this study, we apply Hirota’s direct strategy [23] that provides us with the purely soliton solutions. The Hirota bilinear approach is one of them and is a crucial and straightforward method. It can solve both integrable and nonintegrable equations. The Hirota bilinear approach has the advantage of being an algebraic method as opposed to an analytical one, and it has been used to solve a lot of soliton problems, including the nonlinear Schrödinger equation, the KdV equation, the mKdV equation, the sine-Gordon equation, etc.

Discussion of the Model and Wave Structures

The BBMPB model is a mathematical model that describes the dynamics of ion-acoustic wave structures in a plasma consisting of fluid ions. This equation is a modified version of the classical Burgers equation [8], which is a nonlinear partial differential equation describing the propagation of shocks in a fluid. The BBMPB equation incorporates the effects of dispersion, diffusion, and nonlinearity, which are all important factors that affect the behaviour of ion-acoustic waves in a plasma. The BBMPB equation is given by:
v t v x x t α v x x + λ v x + θ v v x + β v x x x = 0 ,
where v ( x , t ) is the ion-acoustic wave amplitude, t is time, x is the position, and the constants α , λ , θ , and β are parameters that depend on the properties of the plasma. The first term on the left-hand side of the equation represents the dispersive effects of the plasma, while the second term represents the diffusive effects. The third and fourth terms represent the linear damping and driving forces, respectively, while the fifth term represents the nonlinear effects that arise due to the interaction between the waves.
The BBMPB equation can be used to study a wide range of ion-acoustic wave phenomena in plasmas, including the formation of solitons, shock waves, and other nonlinear structures. In particular, the equation can be used to study the propagation of ion-acoustic waves in a plasma with a spatially varying ion density profile, which is a common feature of many plasma systems. The equation can also be used to study the effects of external forces, such as electric fields or magnetic fields, on the dynamics of ion-acoustic waves.
One of the most important features of ion-acoustic waves in plasmas is their ability to form coherent structures, such as solitons and shock waves, that can propagate over long distances without dissipating [1,24]. These structures arise due to the interplay between the dispersive, diffusive, and nonlinear effects in the plasma, and their properties can be studied using the BBMPB equation. Solitons are stable, localized wave packets that maintain their shape as they propagate through the plasma, while shock waves are characterized by a rapid increase in wave amplitude and are often associated with energy dissipation [24,25].
The BBMPB equation can be solved analytically in some special cases, such as when the plasma is homogeneous or when the nonlinearity is weak [26]. The BBMPB equation has been extensively studied in the literature due to its rich and diverse wave structures. In this context, analytical methods have been used to investigate different types of wave structures. In the following, we briefly introduce some of the wave structures that are studied in this research for the BBMPB equation.
Three-wave hypothesis: The three-wave hypothesis is a well-known phenomenon in nonlinear science, which describes the interaction of three waves that satisfy certain resonance conditions [1,27,28].
  • Breather wave: A breather wave is a localized and oscillatory solution that maintains its shape over time [1,9,29].
  • Lump periodic waves: lump periodic waves are periodic solutions that consist of a sequence of wave packets [1,9,29].
  • Mixed–type wave solutions: mixed-type wave solutions are complex and diverse wave structures that have both soliton-like and oscillatory components [1,28,30].
  • Periodic cross kink: a periodic cross kink is a localized wave structure that exhibits a crossing behaviour [1,9,30].
  • Cross-kink rational wave solution: the cross-kink rational wave solution is a type of nonlinear wave solution that is characterized by the presence of two perpendicular kinks in the wave profile, which cross each other at a single point [31,32].
  • M-shaped rational wave solution: the M-shaped rational wave solution is another type of nonlinear wave solution that is characterized by an “M”-shaped profile [31,32].
  • M-shaped rational wave solution with one kink wave: The M-shaped rational wave solution with one kink wave is a variant of the M-shaped solution that includes a single kink in the wave profile [31,32].
  • M-shaped rational wave solution with two kink waves: The M-shaped rational wave solution with two kink waves is a variant of the M-shaped solution that includes two kinks in the wave profile [31,32].

2. Glimpse of the Method

In this section, we use the method to gain the solutions of the ion-acoustic wave structures in fluid ions, described by the nonlinear partial differential equation (NLPDE) in two variables, x and t, given by
Ψ ( v , v t , v x , v x t , v t t , v x x , ) = 0 ,
where v ( x , t ) is the ion-acoustic wave amplitude, t is time, x is the position , Ψ is a polynomial in v = v ( x , t ) with different partial derivatives, in which the highest-order derivatives and nonlinear terms are involved.
We use the transformation v ( x , t ) = V ( η ) , where η = x c t , and V = f ( η ) η f ( η ) in (2), and obtain the following ordinary differential equation:
Ψ ( V , c V , V , c V , c 2 V , V , ) = 0 ,
where Ψ is a polynomial function of V, and we integrate (3) with respect to η and let the integral constants 0 for our convenience; simplifying the resulting equation, we obtain
( c + β ) V α V + ( λ c ) V θ 2 V 2 = 0 .
To find the different forms of solutions to (1), we use the following transformation [1]
V = f η f ,
4 ( β 1 ) f ( η ) 3 + f ( η ) f ( η ) ( 2 α + θ ) f ( η ) ( β 1 ) f ( η ) + 2 f ( η ) 2 ( β 1 ) f ( 3 ) ( η ) α f ( η ) + ( λ + 1 ) f ( η ) = 0 .
In bilinear form, we have
2 f 2 α 2 f η 2 + ( β 1 ) 3 f η 3 + ( λ + 1 ) f η + f f η ( 2 α + θ ) f η 6 ( β 1 ) 2 f η 2 + 4 ( β 1 ) f η 3 = 0 .
Our focus now is on Equation (7), which we use to find the different solutions of the wave structures discussed for Equation (1).

3. Finding the Solutions of the Wave Structures

1.
Multiwave solutions: With the help of the following transformation [1], we are able to use the three wave hypothesis to generate different types of solutions:
f = n 2 cos η ϵ 3 + ϵ 4 + n 1 cosh η ϵ 1 + ϵ 2 + n 3 cosh η ϵ 5 + ϵ 6 .
Substituting Equation (8) in Equation (7), simplifying and collecting similar terms with trigonometric and hyperbolic functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 3 = 0 , ϵ 1 = λ + 1 2 β 1 and ϵ 3 = λ 1 2 β 1 , and putting them in Equation (8) and then in Equation (5) gives
V 1 , 1 ( η ) = λ + 1 n 1 sinh λ + 1 η 2 β 1 ϵ 2 2 β 1 λ 1 n 2 sin λ 1 η 2 β 1 + ϵ 4 2 β 1 n 2 cos λ 1 η 2 β 1 + ϵ 4 + n 1 cosh λ + 1 η 2 β 1 ϵ 2 .
The multiwave solution of Equation (1) is extracted as
v 1 , 1 ( x , t ) = λ + 1 n 1 sinh λ + 1 ( x c t ) 2 β 1 ϵ 2 λ 1 n 2 sin λ 1 ( x c t ) 2 β 1 + ϵ 4 2 β 1 n 2 cos λ 1 ( x c t ) 2 β 1 + ϵ 4 + n 1 cosh λ + 1 ( x c t ) 2 β 1 ϵ 2 .
Set 2: Setting n 1 = 0 , ϵ 3 = λ 1 2 β 1 , ϵ 5 = λ + 1 2 β 1 , and putting them in Equation (8) and then in Equation (5) gives
V 1 , 2 ( η ) = λ + 1 n 3 sinh λ + 1 η 2 β 1 ϵ 6 2 β 1 λ 1 n 2 sin λ 1 η 2 β 1 + ϵ 4 2 β 1 n 2 cos λ 1 η 2 β 1 + ϵ 4 + n 3 cosh λ + 1 η 2 β 1 ϵ 6 .
The multiwave solution of Equation (1) is extracted as
v 1 , 2 ( x , t ) = λ + 1 n 3 sinh λ + 1 ( x c t ) 2 β 1 ϵ 6 λ 1 n 2 sin λ 1 ( x c t ) 2 β 1 + ϵ 4 2 β 1 n 2 cos λ 1 ( x c t ) 2 β 1 + ϵ 4 + n 3 cosh λ + 1 ( x c t ) 2 β 1 ϵ 6 .
2.
Interaction via double exponential form: With the help of the following transformation [1], we generate different types of solutions:
f = n 1 exp η ϵ 1 + ϵ 2 + n 2 exp η ϵ 3 + ϵ 4 .
Substituting Equation (13) in Equation (5), simplifying and collecting similar terms with exponential functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting ϵ 1 = 2 ( λ + 1 ) θ , ϵ 3 = 0 , α = θ 2 , β = 1 , and putting them in Equation (13) and in Equation (5), we obtain
V 2 , 1 ( η ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) η θ θ n 1 e ϵ 2 2 ( λ + 1 ) η θ + n 2 e ϵ 4 .
Thus, the solution of Equation (1) is extracted as
v 2 , 1 ( x , t ) = 2 ( λ + 1 ) n 1 e ϵ 2 θ n 2 e 2 ( λ + 1 ) ( x c t ) θ + ϵ 4 + n 1 e ϵ 2 .
Set 2: Setting ϵ 1 = 0 , ϵ 3 = 2 ( λ + 1 ) θ , α = θ 2 , β = 1 , and putting them in (13) and Equation (5), we obtain
V 2 , 2 ( η ) = 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) η θ θ n 2 e ϵ 4 2 ( λ + 1 ) η θ + n 1 e ϵ 2 .
Thus, the solution of Equation (1) is extracted as
v 2 , 2 ( x , t ) = 2 ( λ + 1 ) n 2 e ϵ 4 θ n 1 e 2 ( λ + 1 ) ( x c t ) θ + ϵ 2 + n 2 e ϵ 4 .
Set 3: Setting ϵ 1 = 2 ( λ + 1 ) θ , ϵ 3 = 2 ( λ + 1 ) θ , and putting them in (13) in Equation (5), we obtain
V 2 , 3 ( η ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) η θ θ 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) η θ θ n 1 e ϵ 2 2 ( λ + 1 ) η θ + n 2 e ϵ 4 2 ( λ + 1 ) η θ .
Thus, the solution of Equation (1) is extracted as
v 2 , 3 ( x , t ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) ( x c t ) θ θ 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) ( x c t ) θ θ n 1 e ϵ 2 2 ( λ + 1 ) ( x c t ) θ + n 2 e ϵ 4 2 ( λ + 1 ) ( x c t ) θ .
3.
Homoclinic breather approach: With the help of the following transformation [1], we generate different types of solutions:
f = n 1 exp r η ϵ 3 + ϵ 4 + exp r η ϵ 1 + ϵ 2 + n 2 cos r η ϵ 5 + ϵ 6 .
Substituting Equation (20) in Equation (7), simplifying and collecting similar terms with exponential, trigonometric, and exponential–trigonometric functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 1 = 0 , ϵ 1 = ( λ + 1 ) 2 ϵ 5 λ + 1 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , r = 2 λ 2 2 λ 1 θ ϵ 5 , and putting them in Equation (20) and then in Equation (5), we obtain
V 3 , 1 ( η ) = 2 ( λ + 1 ) n 1 e r ϵ 4 2 ( λ + 1 ) η θ r θ 2 ( λ + 1 ) e r 2 ( λ + 1 ) η θ r + ϵ 2 θ n 1 e r ϵ 4 2 ( λ + 1 ) η θ r + n 2 cos r ϵ 6 + e r 2 ( λ + 1 ) η θ r + ϵ 2 .
Thus, the homoclinic breather solution of Equation (1) is extracted as
v 3 , 1 ( x , t ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) 2 ϵ 2 θ ϵ 5 n 2 e 2 ( λ + 1 ) ( x c t ) θ sinh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 θ n 2 e 2 ( λ + 1 ) ( x c t ) θ cosh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 + e 2 ( λ + 1 ) 2 ϵ 2 θ ϵ 5 .
Set 2: Setting ϵ 1 = ( λ + 1 ) 2 ϵ 5 λ + 1 , ϵ 3 = ( λ + 1 ) 2 ϵ 5 λ 1 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , r = 2 λ 2 2 λ 1 θ ϵ 5 , and putting them in Equation (20) and then in Equation (5), we obtain
V 3 , 2 ( η ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 5 n 2 e 2 λ η + ( λ + 1 ) 2 ϵ 4 ϵ 5 + η θ sinh G + n 1 θ e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 5 + n 2 e 2 λ η + ( λ + 1 ) 2 ϵ 4 ϵ 5 + η θ cosh G + n 1 ,
where G = 2 ( λ + 1 ) η ϵ 5 + ϵ 6 θ ϵ 5 . Thus, the homoclinic breather solution of Equation (1) is extracted as
v 3 , 2 ( x , t ) = 2 ( λ + 1 ) n 2 exp R sinh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 + e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 5 + n 1 θ n 2 exp R cosh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 + e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 5 + n 1 ,
where R = 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 4 ϵ 5 θ .
Set 3: Setting n 1 = 0 , ϵ 1 = ( λ + 1 ) 2 ϵ 5 λ 1 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( γ + 1 ) , r = 2 λ 2 2 λ 1 θ ϵ 5 and putting them in Equation (20) and then in Equation (5), we obtain
V 3 , 3 ( η ) = 2 ( λ + 1 ) n 2 e 2 λ η + ( λ + 1 ) 2 ϵ 2 ϵ 5 + η θ sinh 2 ( λ + 1 ) η ϵ 5 + ϵ 6 θ ϵ 5 1 θ + θ n 2 e 2 λ η + ( λ + 1 ) 2 ϵ 2 ϵ 5 + η θ cosh 2 ( λ + 1 ) η ϵ 5 + ϵ 6 θ ϵ 5 .
Thus, the homoclinic breather solution of Equation (1) is extracted as
v 3 , 3 ( x , t ) = 2 ( λ + 1 ) n 2 exp 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 2 ϵ 5 θ sinh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 1 θ n 2 exp 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 2 ϵ 5 θ cosh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 + θ .
4.
Mixed-type solutions: With the help of the following transformation [1], we generate different types of solutions:
f = n 1 exp r η ϵ 1 + ϵ 2 + n 2 exp r η ϵ 1 + ϵ 2 + n 3 sin r η ϵ 3 + ϵ 4 + n 4 sinh r η ϵ 5 + ϵ 6 .
Substituting in Equation (27) and then in Equation (7), simplifying and collecting similar terms with exponential, trigonometric, and exponential–trigonometric functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 1 = 0 , n 3 = 0 , ϵ 1 = ( λ + 1 ) 2 ϵ 5 λ 1 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , r = 2 λ 2 + 2 λ + 1 θ ϵ 5 , and putting them in Equation (27) and then in Equation (5) yields
V 4 , 1 ( η ) = 2 ( λ + 1 ) 2 λ 2 + 2 λ + 1 n 2 exp Q ( λ 1 ) θ 2 λ 2 + 2 λ + 1 n 4 cosh 2 λ 2 + 2 λ + 1 η ϵ 5 + ϵ 6 θ ϵ 5 θ n 2 exp Q n 4 sinh 2 λ 2 + 2 λ + 1 η ϵ 5 + ϵ 6 θ ϵ 5 ,
where Q = 2 λ 2 + 2 λ + 1 ( λ + 1 ) 2 η ϵ 5 λ 1 + ϵ 2 θ ϵ 5 .
The mixed-type solution of Equation (1) is extracted as
v 4 , 1 ( x , t ) = 2 ( λ + 1 ) n 2 exp W 2 ( λ + 1 ) 2 n 4 cosh 2 ( λ + 1 ) ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 θ n 2 exp W n 4 sinh 2 ( λ + 1 ) 2 ϵ 5 ( x c t ) + ϵ 6 θ ϵ 5 ,
where W = 2 ( λ + 1 ) ( c t x ) + ( λ + 1 ) 2 ϵ 2 ϵ 5 θ .
Set 2: Setting n 1 = 0 , n 4 = 0 , ϵ 1 = 2 ( λ + 1 ) θ r , ϵ 3 = 2 λ 2 2 λ 1 θ r , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , and putting them in Equation (27) and then in Equation (5) yields
V 4 , 2 ( η ) = 2 λ 2 2 λ 1 n 3 cos r 2 λ 2 2 λ 1 η θ r + ϵ 4 θ 2 ( λ + 1 ) n 2 e r 2 ( λ + 1 ) η θ r + ϵ 2 θ n 3 sin r 2 λ 2 2 λ 1 η θ r + ϵ 4 + n 2 e r 2 ( λ + 1 ) η θ r + ϵ 2 .
The mixed-type solution of Equation (1) is extracted as
v 4 , 2 ( x , t ) = 2 ( λ + 1 ) 2 n 3 e 2 ( λ + 1 ) ( x c t ) θ + r ϵ 2 cos 2 ( λ + 1 ) 2 ( x c t ) θ + r ϵ 4 2 ( λ + 1 ) n 2 θ n 3 e 2 ( λ + 1 ) ( x c t ) θ + r ϵ 2 sin 2 ( λ + 1 ) 2 ( x c t ) θ + r ϵ 4 + n 2 .
5.
Periodic cross kink: With the help of the following transformation [1], we generate different types of solutions:
f = n 1 exp r η ϵ 3 + ϵ 4 + exp r η ϵ 1 + ϵ 2 + n 2 cos r η ϵ 5 + ϵ 6 + n 3 cosh r η ϵ 7 + ϵ 8 + ϵ 9 .
Substituting Equation (32) in Equation (7), simplifying and collecting similar terms with exponential, trigonometric, and exponential–trigonometric functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 2 = 0 , ϵ 1 = ( λ + 1 ) 2 ϵ 7 λ 1 , ϵ 3 = ( λ + 1 ) 2 ϵ 7 λ + 1 , ϵ 9 = 0 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , r = 2 λ 2 + 2 λ + 1 θ ϵ 7 , and putting them in Equation (32) and then in Equation (5) yields
V 5 , 1 ( η ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 7 n 3 e 2 λ η + ( λ + 1 ) 2 ϵ 4 ϵ 7 + η θ sinh 2 ( λ + 1 ) η ϵ 7 + ϵ 8 θ ϵ 7 + n 1 θ e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 7 + n 3 e 2 λ η + ( λ + 1 ) 2 ϵ 4 ϵ 7 + η θ cosh 2 ( λ + 1 ) η ϵ 7 + ϵ 8 θ ϵ 7 + n 1 .
Thus, the periodic cross-kink solution of Equation (1) is extracted as
v 5 , 1 ( x , t ) = 2 ( λ + 1 ) n 3 exp L sinh 2 ( λ + 1 ) ϵ 7 ( x c t ) + ϵ 8 θ ϵ 7 + e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 7 + n 1 θ n 3 exp L cosh 2 ( λ + 1 ) ϵ 7 ( x c t ) + ϵ 8 θ ϵ 7 + e 2 ( λ + 1 ) 2 ϵ 2 + ϵ 4 θ ϵ 7 + n 1 ,
where L = 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 4 ϵ 7 θ .
Set 2: Setting n 1 = 0 , n 2 = 0 , ϵ 1 = ( λ + 1 ) 2 ϵ 7 λ + 1 , ϵ 9 = 0 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , r = 2 λ 2 + 2 λ + 1 θ ϵ 7 , and putting them in Equation (32) and then in Equation (5) yields
V 5 , 2 ( η ) = 2 ( λ + 1 ) n 3 e 2 λ η + ( λ + 1 ) 2 ϵ 2 ϵ 7 + η θ sinh 2 ( λ + 1 ) η ϵ 7 + ϵ 8 θ ϵ 7 1 θ + θ n 3 e 2 λ η + ( λ + 1 ) 2 ϵ 2 ϵ 7 + η θ cosh 2 ( λ + 1 ) η ϵ 7 + ϵ 8 θ ϵ 7 .
The periodic cross-kink solution of Equation (1) is extracted as
v 5 , 2 ( x , t ) = 2 ( λ + 1 ) n 3 exp 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 2 ϵ 7 θ sinh 2 ( λ + 1 ) ϵ 7 ( x c t ) + ϵ 8 θ ϵ 7 1 θ n 3 exp 2 ( λ + 1 ) ( x c t ) + ( λ + 1 ) 2 ϵ 2 ϵ 7 θ cosh 2 ( λ + 1 ) ϵ 7 ( x c t ) + ϵ 8 θ ϵ 7 + θ .
Set 3: Setting n 1 = 0 , n 2 = 0 , ϵ 1 = 2 ( λ + 1 ) θ r , ϵ 7 = 2 λ 2 + 2 λ + 1 θ r , ϵ 9 = 0 , α = θ 8 , β = 32 λ + θ 2 + 32 32 ( λ + 1 ) , and putting them in Equation (32) and then in Equation (5) yields
V 5 , 3 ( η ) = 2 λ 2 + 2 λ + 1 n 3   sinh   r 2 λ 2 + 2 λ + 1 η θ r + ϵ 8 θ 2 ( λ + 1 ) e r 2 ( λ + 1 ) η θ r + ϵ 2 θ n 3 cosh r 2 λ 2 + 2 λ + 1 η θ r + ϵ 8 + e r 2 ( λ + 1 ) η θ r + ϵ 2 .
Thus, the periodic cross-kink solution of Equation (5) is extracted as
v 5 , 3 ( x , t ) = 2 ( λ + 1 ) 2 n 3 e 2 ( λ + 1 ) ( x c t ) θ + r ϵ 2 sinh 2 ( λ + 1 ) 2 ( x c t ) θ + r ϵ 8 + λ + 1 θ n 3 e 2 ( λ + 1 ) ( x c t ) θ + r ϵ 2 cosh 2 ( λ + 1 ) 2 ( x c t ) θ + r ϵ 8 + θ .
6.
Cross-kink rational wave solution: With the help of the following transformation [31], we generate different types of solutions:
f = n 1 exp η ϵ 1 + ϵ 2 + exp η ϵ 1 + ϵ 2 + η r 1 + r 2 2 + η r 3 + r 4 2 + r 5 .
Substituting Equation (39) in Equation (7), simplifying and collecting similar terms with exponential functions and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 1 = 0 , ϵ 1 = 2 ( λ + 1 ) θ , r 2 = r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) , r 3 = 0 , α = θ 2 , β = 1 , and putting them in Equation (39) and in Equation (5), we obtain
V 6 , 1 ( η ) = 4 ( λ + 1 ) 3 2 ( λ + 1 ) 2 θ + θ 3 r 1 2 e 2 ( λ + 1 ) η θ + ϵ 2 .
Thus, the cross-kink rational wave solution of Equation (1) is extracted as
v 6 , 1 ( x , t ) = 4 ( λ + 1 ) 3 2 ( λ + 1 ) 2 θ + θ 3 r 1 2 e 2 ( λ + 1 ) ( x c t ) θ + ϵ 2 .
Set 2: Setting ϵ 1 = 0 , r 2 = r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) , r 3 = 0 , α = θ 2 , β = 1 , and putting them in Equation (39) and in Equation (7), we obtain
V 6 , 2 ( η ) = 2 r 1 ( 2 λ + 2 ) r 1 2 ( λ + 1 ) η r 1 r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) n 1 e ϵ 2 + η r 1 r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) 2 + r 4 2 + r 5 + e ϵ 2 .
Thus, the cross-kink rational wave solution of Equation (1) is extracted as
v 6 , 2 ( x , t ) = 2 r 1 ( 2 λ + 2 ) r 1 2 ( λ + 1 ) r 1 ( x ct ) r 1 ( 2 λ ( x ct ) + 2 ( x ct ) + θ ) 2 ( λ + 1 ) r 1 ( x ct ) r 1 ( 2 λ ( x ct ) + 2 ( x ct ) + θ ) 2 ( λ + 1 ) 2 + n 1 e ϵ 2 + r 4 2 + r 5 + e ϵ 2 .
Set 3: Setting n 1 = 0 , ϵ 1 = 2 ( λ + 1 ) θ , r 1 = i r 3 , r 4 = i r 2 , r 5 = 0 , and putting them in Equation (39) and in Equation (7), we obtain
V 6 , 3 ( η ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) η θ ϵ 2 θ + 2 r 3 η r 3 + i r 2 2 i r 3 r 2 i η r 3 e 2 ( λ + 1 ) η θ ϵ 2 + η r 3 + i r 2 2 + r 2 i η r 3 2 .
Thus, the cross-kink rational wave solution of Equation (1) is extracted as
v 6 , 3 ( x , t ) = 2 r 3 r 3 ( x c t ) + i r 2 2 i r 3 r 2 i r 3 ( x c t ) 2 ( λ + 1 ) e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 θ r 3 ( x c t ) + i r 2 2 + r 2 i r 3 ( x c t ) 2 + e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 .
7.
M-shaped rational wave solution: With the help of the following transformation [31], we generate different types of solutions:
f = η r 1 + r 2 2 + η r 3 + r 4 2 + r 5
Substituting Equation (46) in Equation (7), simplifying and collecting similar terms, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting r 2 = η r 1 2 η r 3 2 r 3 r 4 r 1 , r 5 = r 1 2 + r 3 2 η r 3 + r 4 2 r 1 2 , and putting them in Equation (46) and in Equation (5), we obtain
V 7 , 1 ( η ) = 2 r 1 + r 1 2 r 3 2 r 1 η r 1 + η r 1 2 η r 3 2 r 3 r 4 r 1 + 2 r 3 η r 3 + r 4 2 r 3 r 1 2 + r 3 2 η r 3 + r 4 r 1 2 η r 1 + η r 1 2 η r 3 2 r 3 r 4 r 1 2 + η r 3 + r 4 2 r 1 2 + r 3 2 η r 3 + r 4 2 r 1 2 .
Thus, the M-shaped rational wave solution of Equation (1) is extracted as
v 7 , 1 ( x , t ) = 2 r 1 + r 1 2 r 3 2 r 1 H + 2 r 3 r 3 ( x c t ) + r 4 2 r 3 r 1 2 + r 3 2 r 3 ( x c t ) + r 4 r 1 2 H 2 + r 3 ( x c t ) + r 4 2 r 1 2 + r 3 2 r 3 ( x c t ) + r 4 2 r 1 2 ,
where H = r 1 2 ( ( x c t ) ) η r 3 2 r 3 r 4 r 1 + r 1 ( x c t ) .
Set 2: Setting r 2 = i r 4 , r 3 = i r 1 , and putting them in Equation (46) and in Equation (5), we obtain
V 7 , 2 ( η ) = 2 r 1 η r 1 + i r 4 2 i r 1 r 4 i η r 1 r 4 i η r 1 2 + η r 1 + i r 4 2 + r 5 .
Thus, the M-shaped rational wave solution of Equation (1) is extracted as
v 7 , 2 ( x , t ) = 2 r 1 r 1 ( x ct ) + i r 4 2 i r 1 r 4 i r 1 ( x ct ) r 4 i r 1 ( x ct ) 2 + r 1 ( x ct ) + i r 4 2 + r 5 .
Set 3: Setting r 2 = η r 1 2 η r 3 2 r 3 r 4 r 1 , α = 0 , and putting them in Equation (46) and in Equation (7), we obtain
Ψ 7 , 3 ( η ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) η θ ϵ 2 θ + 2 r 3 η r 3 + i r 2 2 i r 3 r 2 i η r 3 e 2 ( λ + 1 ) η θ ϵ 2 + η r 3 + i r 2 2 + r 2 i η r 3 2 .
Thus, the M-shaped rational wave solution of Equation (1) is extracted as
v 7 , 3 ( x , t ) = 2 r 3 r 3 ( x c t ) + i r 2 2 i r 3 r 2 i r 3 ( x c t ) 2 ( λ + 1 ) e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 θ r 3 ( x c t ) + i r 2 2 + r 2 i r 3 ( x c t ) 2 + e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 .
8.
M-shaped rational wave solution with one kink wave: With the help of the following transformation [31], we generate different types of solutions:
f = n 1 exp η ϵ 1 + ϵ 2 + η r 1 + r 2 2 + η r 3 + r 4 2 + r 5 .
Substituting Equation (53) in Equation (7), simplifying and collecting similar terms with exponential functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting n 1 = 0 , ϵ 1 = 2 ( λ + 1 ) θ , r 2 = r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) , r 3 = 0 , α = θ 2 , β = 1 , and putting them in Equation (53) and in Equation (5), we obtain
V 8 , 1 ( η ) = 4 ( λ + 1 ) 3 2 ( λ + 1 ) 2 θ + θ 3 r 1 2 e 2 ( λ + 1 ) η θ + ϵ 2 .
Thus, the solution of Equation (1) is extracted as
v 8 , 1 ( x , t ) = 4 ( λ + 1 ) 3 2 ( λ + 1 ) 2 θ + θ 3 r 1 2 e 2 ( λ + 1 ) ( x c t ) θ + ϵ 2 .
Set 2:  ϵ 1 = 0 , r 2 = r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) , r 3 = 0 , α = θ 2 , β = 1 , and putting them in Equation (53) in Equation (7), we obtain
V 8 , 2 ( η ) = 2 r 1 ( 2 λ + 2 ) r 1 2 ( λ + 1 ) η r 1 r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) n 1 e ϵ 2 + η r 1 r 1 ( 2 λ η + 2 η + θ ) 2 ( λ + 1 ) 2 + r 4 2 + r 5 + e ϵ 2 .
Thus, the solution of Equation (1) is extracted as
v 8 , 2 ( x , t ) = 2 r 1 ( 2 λ + 2 ) r 1 2 ( λ + 1 ) r 1 ( x ct ) r 1 ( 2 λ ( x ct ) + 2 ( x ct ) + θ ) 2 ( λ + 1 ) r 1 ( x ct ) r 1 ( 2 λ ( x ct ) + 2 ( x ct ) + θ ) 2 ( λ + 1 ) 2 + n 1 e ϵ 2 + r 4 2 + r 5 + e ϵ 2 .
Set 3: Setting n 1 = 0 , ϵ 1 = 2 ( λ + 1 ) θ , r 1 = i r 3 , r 4 = i r 2 , r 5 = 0 , and putting them in Equation (53) and in Equation (7), we obtain
V 8 , 3 ( η ) = 2 ( λ + 1 ) e 2 ( λ + 1 ) η θ ϵ 2 θ + 2 r 3 η r 3 + i r 2 2 i r 3 r 2 i η r 3 e 2 ( λ + 1 ) η θ ϵ 2 + η r 3 + i r 2 2 + r 2 i η r 3 2 .
Thus, the solution of Equation (1) is extracted as
v 8 , 3 ( x , t ) = 2 r 3 r 3 ( x c t ) + i r 2 2 i r 3 r 2 i r 3 ( x c t ) 2 ( λ + 1 ) e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 θ r 3 ( x c t ) + i r 2 2 + r 2 i r 3 ( x c t ) 2 + e 2 ( λ + 1 ) ( x c t ) θ ϵ 2 .
9.
M-shaped rational wave solution with two kink waves: With the help of the following transformation [31], we generate different types of solutions:
f = n 1 exp η ϵ 1 + ϵ 2 + n 2 exp η ϵ 3 + ϵ 4 + η r 1 + r 2 2 + η r 3 + r 4 2 + r 5 .
Substituting Equation (60) in Equation (7), simplifying and collecting similar terms with exponential functions, and equating the coefficients of each obtained expression to zero, we obtained a system of equations and simplified it with the help of Mathematica to gain the following different sets of unknown constants:
Set 1: Setting ϵ 1 = 2 ( λ + 1 ) θ , ϵ 3 = 2 ( λ + 1 ) θ , r 1 = 0 , r 3 = 0 , r 5 = r 2 2 r 4 2 , and putting them in Equation (60) and in Equation (5), we obtain
V 9 , 1 ( η ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) η θ θ 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) η θ θ n 1 e ϵ 2 2 ( λ + 1 ) η θ + n 2 e ϵ 4 2 ( λ + 1 ) η θ .
Thus, the solution of Equation (1) is extracted as
v 9 , 1 ( x , t ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) ( x c t ) θ θ 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) ( x c t ) θ θ n 1 e ϵ 2 2 ( λ + 1 ) ( x c t ) θ + n 2 e ϵ 4 2 ( λ + 1 ) ( x c t ) θ .
Set 2: Setting ϵ 1 = 2 ( λ + 1 ) θ , ϵ 3 = 2 ( λ + 1 ) θ , r 2 = i r 4 , r 3 = i r 1 , r 5 = 0 , and putting them in Equation (60) and in Equation (5), we obtain
V 9 , 2 ( η ) = 2 ( λ + 1 ) n 1 e ϵ 2 2 ( λ + 1 ) η θ θ 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) η θ θ + 2 i r 1 r 4 + i η r 1 + 2 r 1 η r 1 i r 4 n 1 e ϵ 2 2 ( λ + 1 ) η θ + n 2 e ϵ 4 2 ( λ + 1 ) η θ + r 4 + i η r 1 2 + η r 1 i r 4 2 .
Thus, the solution of Equation (1) is extracted as
v 9 , 2 ( x , t ) = 2 ( λ + 1 ) n 1 e ϵ 2 A θ 2 ( λ + 1 ) n 2 e ϵ 4 A θ + 2 i r 1 r 4 + i r 1 ( x ct ) + 2 r 1 r 1 ( x ct ) i r 4 n 1 e ϵ 2 A + n 2 e ϵ 4 A + r 4 + i r 1 ( x ct ) 2 + r 1 ( x ct ) i r 4 2 ,
where A = 2 ( λ + 1 ) ( x ct ) θ .
Set 3: Setting n 1 = 0 , ϵ 3 = 2 ( λ + 1 ) θ , r 1 = i r 3 , r 4 = i r 2 , r 5 = 0 , and putting them in Equation (60) and in Equation (5), we obtain
V 9 , 3 ( η ) = 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) η θ θ + 2 r 3 η r 3 + i r 2 2 i r 3 r 2 i η r 3 n 2 e ϵ 4 2 ( λ + 1 ) η θ + η r 3 + i r 2 2 + r 2 i η r 3 2 .
Thus, the solution of Equation (1) is extracted as
v 9 , 3 ( x , t ) = 2 ( λ + 1 ) n 2 e ϵ 4 2 ( λ + 1 ) ( x c t ) θ θ + 2 r 3 r 3 ( x c t ) + i r 2 2 i r 3 r 2 i r 3 ( x c t ) r 3 ( x c t ) + i r 2 2 + r 2 i r 3 ( x c t ) 2 + e ϵ 4 2 ( λ + 1 ) ( x c t ) θ .

4. Graphical Presentations

Finding soliton solutions was the main goal of this research. Optics is one of the sciences that extensively studies the intriguing physics phenomenon known as a soliton. A self-reinforcing solitary wave is referred to as a soliton if it can keep its shape and speed while moving across a medium without separating or dissipating. Solitons have an unusual behaviour that makes them remarkably stable and enables them to maintain their shape across great distances. These solitons can arise in a variety of waveguide designs, including planar waveguides or ion-acoustic waveguides. Solitons’ capacity to keep their shape and propagate unaltered even in the presence of nonlinear influences and dispersion is their defining feature. Solitons have a physical meaning when they can balance the opposing effects of dispersion and nonlinearity. Over time, dispersion tends to spread out a pulse, causing it to enlarge and change shape. The pulse can be compressed by a self-focusing action brought on by nonlinearity, which, on the other hand, can combat dispersion. For solitons to emerge, nonlinearity and dispersion must coexist in harmony. Ionic-acoustic solitons are created when the nonlinear effects properly balance out the dispersion, creating a localised waveform that is stable and unaltered at vast distances. Solitons are advantageous for numerous applications, including optical signal processing, ultrafast laser pulse propagation, and high-capacity optical communications. Both theoretical and experimental physics have greatly benefited from the discovery and understanding of solitons. They have aided in the growth of uses of ion-acoustic waves and their prospective applications in a variety of fields by providing insights into the interaction between linear and nonlinear effects in wave propagation. Solitons can behave very differently depending on the physical impacts of various dispersion parameters on their dynamics or characteristics.
Dispersion refers to the dependence of the wave’s phase or group velocity on its frequency, and different types of dispersion can arise in different waveguide structures or propagation media. We explored the graphical representations of the aforesaid solutions. The various wave structures below demonstrate graphical depictions of some of the solutions obtained above. Graphs have been extensively utilized to explain the dynamics and distinctive appearance of solutions generated from the BBMPB equation. The results obtained here show that the system has a very diversified wave shape for use in fluid ion applications. Figure 1 and Figure 2 were drawn for the multiwave solutions while Figure 3 shows the kink-type solution for the double rational form solution. Figure 4 and Figure 5 show the breather waves, and Figure 6 corresponds to the mixed-type solution. Figure 7, Figure 8 and Figure 9 represent the physical behaviour of the periodic cross-kink solutions. These solutions are effective in the further study of dynamical systems.

5. Conclusions

The BBMPB model is a useful model for understanding the behaviour of ion-acoustic waves in fluid ions. The ion-acoustic wave is a fundamental mode of oscillation in plasmas and is characterized by its dispersion relation, which relates the wave frequency to the wave number. The study of ion-acoustic waves is of great importance in plasma physics because it provides insights into the basic plasma processes, such as energy transport, wave–particle interactions, and turbulence. The Hirota bilinear transformation technique was used to construct different types of the solutions, and this method gave us the exact ionic wave structures. Different types of wave structures were constructed successfully in the forms of breather waves, lump periodic waves, mixed-type wave solutions, cross-kink rational wave solutions, M-shaped rational wave solutions, and M-shaped rational wave solutions with one or two kink waves. Overall, the study of ion-acoustic wave structures using the BBMPB equation provides valuable insights into the behaviour of plasma waves and the underlying physics of ion-acoustic wave phenomena. The different types of wave solutions studied through this equation gives us a clear understanding of the complex dynamics of these wave structures and may have important implications for the design of plasma-based technologies and experiments.

Author Contributions

The individual contributions of the authors are as follows: methodology, B.C.; writing—original draft preparation, A.A.; software, B.C.; formal analysis, M.Z.B.; investigation, M.Z.B.; resources, N.A. and A.A.; data curation, N.A.; writing—review and editing, A.C. and J.R.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the anonymous reviewers for their comments and suggestions, which have improved the final version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Younas, U.; Ren, J.; Baber, M.Z.; Yasin, M.W.; Shahzad, T. Ion-acoustic wave structures in the fluid ions modeled by higher dimensional generalized Korteweg-de Vries–Zakharov–Kuznetsov equation. J. Ocean Eng. Sci. 2022, 7, 1–13. [Google Scholar] [CrossRef]
  2. Goswami, A.; Singh, J.; Kumar, D.; Gupta, S. An efficient analytical technique for fractional partial differential equations occurring in ion acoustic waves in plasma. J. Ocean Eng. Sci. 2019, 4, 85–99. [Google Scholar] [CrossRef]
  3. Dubinov, A.E. Gas-dynamic approach to the theory of non-linear ion-acoustic waves in plasma with Kaniadakis’ distributed species. Adv. Space Res. 2023, 71, 1108–1115. [Google Scholar] [CrossRef]
  4. Usman, M.; Hussain, A.; Zaman, F.D.; Khan, I.; Eldin, S.M. Reciprocal Bäcklund transformations and travelling wave structures of some nonlinear pseudo-parabolic equations. Partial Differ. Equ. Appl. Math. 2023, 7, 100490. [Google Scholar] [CrossRef]
  5. Lipatov, A.S. The Hybrid Multiscale Simulation Technology: An Introduction with Application to Astrophysical and Laboratory Plasmas; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  6. Xiang, C.; Wang, H. New Exact Solutions for Benjamin-Bona-Mahony-Burgers Equation. Open J. Appl. Sci. 2020, 10, 543–550. [Google Scholar] [CrossRef]
  7. Yang, L. Application of classification of traveling wave solutions to the Zakhrov-Kuznetsov-Benjamin-Bona-Mahony equation. Appl. Math. 2014, 5, 1432. [Google Scholar] [CrossRef]
  8. Akcagil, S.; Aydemir, T.; Gozukizil, O.F. Exact travelling wave solutions of nonlinear pseudoparabolic equations by using the G G Expansion Method. New Trends Math. Sci. 2016, 4, 51–66. [Google Scholar] [CrossRef]
  9. Liu, J.G.; Du, J.Q.; Zeng, Z.F.; Ai, G.P. Exact periodic cross-kink wave solutions for the new (2+1)-dimensional KdV equation in fluid flows and plasma physics. Chaos Interdiscip. J. Nonlinear Sci. 2016, 26, 103114. [Google Scholar] [CrossRef]
  10. Song, L.M.; Yang, Z.J.; Li, X.L.; Zhang, S.M. Coherent superposition propagation of Laguerre–Gaussian and Hermite–Gaussian solitons. Appl. Math. Lett. 2020, 102, 106114. [Google Scholar] [CrossRef]
  11. Shen, S.; Yang, Z.J.; Pang, Z.G.; Ge, Y.R. The complex-valued astigmatic cosine-Gaussian soliton solution of the nonlocal nonlinear Schrödinger equation and its transmission characteristics. Appl. Math. Lett. 2022, 125, 107755. [Google Scholar] [CrossRef]
  12. Shen, S.; Yang, Z.; Li, X.; Zhang, S. Periodic propagation of complex-valued hyperbolic-cosine-Gaussian solitons and breathers with complicated light field structure in strongly nonlocal nonlinear media. Commun. Nonlinear Sci. Numer. Simul. 2021, 103, 106005. [Google Scholar] [CrossRef]
  13. Li, X.L.; Guo, R. Interactions of Localized Wave Structures on Periodic Backgrounds for the Coupled Lakshmanan–Porsezian–Daniel Equations in Birefringent Optical Fibers. Ann. Phys. 2023, 535, 2200472. [Google Scholar] [CrossRef]
  14. Zou, Z.; Guo, R. The Riemann–Hilbert approach for the higher-order Gerdjikov–Ivanov equation, soliton interactions and position shift. Commun. Nonlinear Sci. Numer. Simul. 2023, 124, 107316. [Google Scholar] [CrossRef]
  15. Zhang, S.; Zhu, F.; Xu, B. Localized Symmetric and Asymmetric Solitary Wave Solutions of Fractional Coupled Nonlinear Schrödinger Equations. Symmetry 2023, 15, 1211. [Google Scholar] [CrossRef]
  16. Taghizadeh, N.; Mirzazadeh, M. The direct algebraic method to complex nonlinear partial differential equations. Int. J. Appl. Math. Comput. 2013, 5, 12–16. [Google Scholar]
  17. Sulaiman, T.A.; Younas, U.; Younis, M.; Ahmad, J.; Rehman, S.U.; Bilal, M.; Yusuf, A. Modulation instability analysis, optical solitons and other solutions to the (2+ 1)-dimensional hyperbolic nonlinear Schrodinger’s equation. Comput. Methods Differ. Equ. 2022, 10, 179–190. [Google Scholar]
  18. Seadawy, A.R.; Younis, M.; Baber, M.Z.; Rizvi, S.T.; Iqbal, M.S. Diverse acoustic wave propagation to confirmable time–space fractional KP equation arising in dusty plasma. Commun. Theor. Phys. 2021, 73, 115004. [Google Scholar] [CrossRef]
  19. Shahzad, T.; Baber, M.Z.; Ahmad, M.O.; Ahmed, N.; Akgül, A.; Ali, S.M.; Ali, M.; El Din, S.M. On the analytical study of predator–prey model with Holling-II by using the new modified extended direct algebraic technique and its stability analysis. Results Phys. 2023, 51, 106677. [Google Scholar] [CrossRef]
  20. Baber, M.Z.; Seadway, A.R.; Iqbal, M.S.; Ahmed, N.; Yasin, M.W.; Ahmed, M.O. Comparative analysis of numerical and newly constructed soliton solutions of stochastic Fisher-type equations in a sufficiently long habitat. Int. J. Mod. Phys. 2023, 37, 2350155. [Google Scholar] [CrossRef]
  21. Younas, U.; Sulaiman, T.A.; Ren, J. On the optical soliton structures in the magneto electro-elastic circular rod modeled by nonlinear dynamical longitudinal wave equation. Opt. Quantum Electron. 2022, 54, 688. [Google Scholar] [CrossRef]
  22. Feng, Q. A new approach for seeking coefficient function solutions of conformable fractional partial differential equations based on the Jacobi elliptic equation. Chin. J. Phys. 2018, 56, 2817–2828. [Google Scholar] [CrossRef]
  23. Hirota, R. The Direct Method in Soliton Theory (No. 155); Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  24. Ghosh, S.; Bharuthram, R. Ion acoustic solitons and double layers in electron–positron–ion plasmas with dust particulates. Astrophys. Space Sci. 2008, 314, 121–127. [Google Scholar] [CrossRef]
  25. Petviashvili, V.I.; Pokhotelov, O.A. Solitary Waves in Plasmas and in the Atmosphere; Taylor & Francis: London, UK, 1992. [Google Scholar]
  26. Yuan, W.; Xiao, B.; Wu, Y.; Qi, J. The general traveling wave solutions of the Fisher type equations and some related problems. J. Inequal. Appl. 2014, 2014, 1–15. [Google Scholar] [CrossRef]
  27. Wang, X.; Cao, J.; Chen, Y. Higher-order rogue wave solutions of the three-wave resonant interaction equation via the generalized Darboux transformation. Phys. Scr. 2015, 90, 105. [Google Scholar] [CrossRef]
  28. Liu, J.G.; Du, J.Q.; Zeng, Z.F.; Nie, B. New three-wave solutions for the (3+1)-dimensional Boiti–Leon–Manna–Pempinelli equation. Nonlinear Dyn. 2017, 88, 655–661. [Google Scholar] [CrossRef]
  29. Mihalache, D. Localized structures in optical and matter-wave media: A selection of recent studies. Rom. Rep. Phys. 2021, 73, 403. [Google Scholar]
  30. Ma, H.; Zhang, C.; Deng, A. New periodic wave, cross-kink wave, breather, and the interaction phenomenon for the (2+ 1)-dimensional Sharmo–Tasso–Olver equation. Complexity 2020, 2020, 4270906. [Google Scholar] [CrossRef]
  31. Alsallami, S.A.; Rizvi, S.T.; Seadawy, A.R. Study of stochastic–fractional Drinfel’d–Sokolov–Wilson equation for M-shaped rational, homoclinic breather, periodic and Kink-Cross rational solutions. Mathematics 2023, 11, 1504. [Google Scholar] [CrossRef]
  32. Seadawy, A.R.; Rizvi, S.T.; Younis, M.; Ashraf, M.A. Breather, multi-wave, periodic-cross kink, M-shaped and interactions solutions for perturbed NLSE with quadratic cubic nonlinearity. Opt. Quantum Electron. 2021, 53, 1–14. [Google Scholar] [CrossRef]
Figure 1. The graph depicts the solution of v 1 , 1 ( x , t ) by choosing β = −0.5, c = 1.1, λ = 1.5, n 1 = 1.5, n 2 = 1.5, ϵ 2 = 0.5, and ϵ 4 = 1.5.
Figure 1. The graph depicts the solution of v 1 , 1 ( x , t ) by choosing β = −0.5, c = 1.1, λ = 1.5, n 1 = 1.5, n 2 = 1.5, ϵ 2 = 0.5, and ϵ 4 = 1.5.
Symmetry 15 01682 g001
Figure 2. The graph depicts the solution of v 1 , 2 ( x , t ) by choosing β = 0.5, c = 2.1, λ = 1.5, n 2 = 0.5, n 3 = 0.5, ϵ 4 = 1.5, and ϵ 6 = 2.
Figure 2. The graph depicts the solution of v 1 , 2 ( x , t ) by choosing β = 0.5, c = 2.1, λ = 1.5, n 2 = 0.5, n 3 = 0.5, ϵ 4 = 1.5, and ϵ 6 = 2.
Symmetry 15 01682 g002
Figure 3. The graph depicts the solution of v 2 , 1 ( x , t ) by choosing c = 1.01, λ = 1.5, θ = 1.2, n 1 = 0.5, n 2 = 0.5, ϵ 2 = 0.5, and ϵ 4 = 0.5.
Figure 3. The graph depicts the solution of v 2 , 1 ( x , t ) by choosing c = 1.01, λ = 1.5, θ = 1.2, n 1 = 0.5, n 2 = 0.5, ϵ 2 = 0.5, and ϵ 4 = 0.5.
Symmetry 15 01682 g003
Figure 4. The graph depicts the solution of v 3 , 1 ( x , t ) by choosing c = 1.1, λ = 0.2, θ = 2.5, n 2 = −1.0, ϵ 2 = 1.5, ϵ 5 = 0.5, and ϵ 6 = 0.1.
Figure 4. The graph depicts the solution of v 3 , 1 ( x , t ) by choosing c = 1.1, λ = 0.2, θ = 2.5, n 2 = −1.0, ϵ 2 = 1.5, ϵ 5 = 0.5, and ϵ 6 = 0.1.
Symmetry 15 01682 g004
Figure 5. The graph depicts the solution of v 3 , 2 ( x , t ) by choosing c = 0.985, λ = 1.5, θ = 1.2, n 1 = 1.5, n 2 = 0.67, ϵ 2 = 0.5, ϵ 4 = 1.5, ϵ 5 = 0.9, and ϵ 6 = 0.1.
Figure 5. The graph depicts the solution of v 3 , 2 ( x , t ) by choosing c = 0.985, λ = 1.5, θ = 1.2, n 1 = 1.5, n 2 = 0.67, ϵ 2 = 0.5, ϵ 4 = 1.5, ϵ 5 = 0.9, and ϵ 6 = 0.1.
Symmetry 15 01682 g005
Figure 6. The graph depicts the solution of v 4 , 1 ( x , t ) by choosing c = 1.1, λ = 0.5, θ = 1.2, n 2 = 0.5, n 3 = 0.05, r = 2.1, ϵ 2 = −0.5, and ϵ 4 = −1.5.
Figure 6. The graph depicts the solution of v 4 , 1 ( x , t ) by choosing c = 1.1, λ = 0.5, θ = 1.2, n 2 = 0.5, n 3 = 0.05, r = 2.1, ϵ 2 = −0.5, and ϵ 4 = −1.5.
Symmetry 15 01682 g006
Figure 7. The graph depicts the solution of v 5 , 1 ( x , t ) by choosing c = 1.1, λ = 0.5, θ = 1.2, n 1 = 1.5, n 3 = 3.5, ϵ 2 = 0.5, ϵ 4 = 1.5, ϵ 7 = 1.5, and ϵ 8 = 1.0.
Figure 7. The graph depicts the solution of v 5 , 1 ( x , t ) by choosing c = 1.1, λ = 0.5, θ = 1.2, n 1 = 1.5, n 3 = 3.5, ϵ 2 = 0.5, ϵ 4 = 1.5, ϵ 7 = 1.5, and ϵ 8 = 1.0.
Symmetry 15 01682 g007
Figure 8. The graph depicts the solution of v 5 , 2 ( x , t ) by choosing c = −1.1, λ = 0.5, θ = 2.2, n 3 = 1.5, ϵ 2 = 0.5, ϵ 7 = 3.5, and ϵ 8 =1.0.
Figure 8. The graph depicts the solution of v 5 , 2 ( x , t ) by choosing c = −1.1, λ = 0.5, θ = 2.2, n 3 = 1.5, ϵ 2 = 0.5, ϵ 7 = 3.5, and ϵ 8 =1.0.
Symmetry 15 01682 g008
Figure 9. The graph depicts the solution of v 5 , 3 ( x , t ) by choosing c = −1.5, λ = 1.5, θ = 2.2, n 3 = 1.5, r = 2.1, ϵ 2 = 0.5, and ϵ 8 =1.0.
Figure 9. The graph depicts the solution of v 5 , 3 ( x , t ) by choosing c = −1.5, λ = 1.5, θ = 2.2, n 3 = 1.5, r = 2.1, ϵ 2 = 0.5, and ϵ 8 =1.0.
Symmetry 15 01682 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ceesay, B.; Baber, M.Z.; Ahmed, N.; Akgül, A.; Cordero, A.; Torregrosa, J.R. Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique. Symmetry 2023, 15, 1682. https://doi.org/10.3390/sym15091682

AMA Style

Ceesay B, Baber MZ, Ahmed N, Akgül A, Cordero A, Torregrosa JR. Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique. Symmetry. 2023; 15(9):1682. https://doi.org/10.3390/sym15091682

Chicago/Turabian Style

Ceesay, Baboucarr, Muhammad Zafarullah Baber, Nauman Ahmed, Ali Akgül, Alicia Cordero, and Juan R. Torregrosa. 2023. "Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique" Symmetry 15, no. 9: 1682. https://doi.org/10.3390/sym15091682

APA Style

Ceesay, B., Baber, M. Z., Ahmed, N., Akgül, A., Cordero, A., & Torregrosa, J. R. (2023). Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique. Symmetry, 15(9), 1682. https://doi.org/10.3390/sym15091682

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop