On traveling wave solutions with stability and phase plane analysis for the modified Benjamin-Bona-Mahony equation

Abstract

The modified Benjamin-Bona-Mahony (mBBM) model is utilized in the optical illusion field to describe the propagation of long waves in a nonlinear dispersive medium during a visual illusion (Khater 2021). This article investigates the mBBM equation through the utilization of the rational -expansion technique to derive new analytical wave solutions. The analytical solutions we have obtained comprise hyperbolic, trigonometric, and rational functions. Some of these exact solutions closely align with previously published results in specific cases, affirming the validity of our other solutions. To provide insights into diverse wave propagation characteristics, we have conducted an in-depth analysis of these solutions using 2D, 3D, and density plots. We also investigated the effects of various parameters on the characteristics of the obtained wave solutions of the model. Moreover, we employed the techniques of linear stability to perform stability analysis of the considered model. Additionally, we have explored the stability of the associated dynamical system through the application of phase plane theory. This study also demonstrates the efficacy and capabilities of the rational -expansion approach in analyzing and extracting soliton solutions from nonlinear partial differential equations.

1. Introduction

The study of nonlinear partial differential equations (NLPDEs) is now regarded as one of the most important and fundamental areas of science, particularly for uncovering previously unknown characteristics of complex phenomena in a variety of disciplines, including chemical physics, biology, plasma physics, hydrodynamics, condensed matter physics, fluid mechanics, ocean physics, nuclear physics, optical physics, and more. In order to reflect many processes in nature, physics, and engineering, numerous nonlinear models have been developed. Numerous researchers have developed different methods and approaches to examine analytical solutions considering the significance of these models and their varied applications. The F-expansion approach [1], Lie symmetry technique [2], Kudryashov process [3], -expansion technique [4], rational sine–cosine scheme [5], modified Jacobi elliptic expansion method [6], exponential-expansion scheme [7], generalized exponential method [8], the hyperbolic function method [9], positive quadratic function method [10], the inverse scattering approach [11], generalized Khater technique [12], test function method [13], improved tanh method [14], different form of the -expansion methods [15–18], and other such methods are examples of these approaches.

The study on solutions, structures, interactions, and further properties of solitary wave and soliton attained much concentration, and different meaningful results have been successfully derived [19–26]. In recent years, significant advancements have been made in the study of soliton solutions through the utilization of the Hirota bilinear method [27]. Noteworthy contributions have emerged from several recent studies, shedding new light on the method’s applicability. For example, the analysis of multi-soliton solutions and soliton molecules in the context of the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli equation underscores the usefulness of the Hirota bilinear method [28]. Furthermore, the exploration of novel soliton solutions to the extended (3+1)-dimensional Jimbo-Miwa equation using the Hirota bilinear method demonstrated its effectiveness in uncovering complex nonlinear dynamics [29].

The regularized long-wave equation, commonly referred to as the Benjamin–Bona–Mahony (BBM) equation, is given by(1)It was originally explored by Benjamin et al. [30] as a novel model for describing long waves, providing an alternative to the Korteweg–de Vries (KdV) equation. To discover new properties of soliton waves, Eq (1) is derived as an improved model of the well-known KdV model. Significant importance is captured by this equation in the modeling of nonlinear dispersive systems. The BBM equation is applied to study hydromagnetic waves in cold plasma, long-wavelength surface waves in liquids, acoustic waves in anharmonic crystals, and acoustic-gravity waves in compressible fluids [31]. This model is utilized to study the dynamical properties of extended surface gravity waves with small amplitudes as they propagate in (1+1)-dimensions. Additionally, the BBM equation has also been used in the study of undular bores. The generalized Benjamin–Bona–Mahony equation is given by(2)where β is an arbitrary constant. For β = 1, the Eq (2) is called modified Benjamin–Bona–Mahony (mBBM) equation. In this paper, we have investigated the mBBM equation, which can be expressed as [32, 33]:(3)where α, β and γ denote any real constants. Similarly, the Eq (3) can be viewed as an alternative formulation of the modified Korteweg-de Vries equation. The mBBM equation also finds application in the realm of optical illusions, where light, color, and patterns are skillfully combined to craft images capable of deceiving the human brain, as demonstrated in prior studies [34–36]. An optical illusion is defined as an illusion created by the visual system and characterized by a perception of the visual world that may appear to deviate from reality.

So far, numerous scholars have made significant contributions to solve the mBBM equation. The -function method is utilized by Aslan [32] to generate exact and explicit solutions. Abbasbandy and Shirzadi [31] employed the first integral method to solve mBBM model, and the exp(−ϕ(ξ))-expansion process was employed by Khan and Akbar [37] for its study. Some new solitonary solutions for the mBBM model were obtained by Noor et al. [38] using exp-function method. Song [39] employed the bifurcation method of dynamical system to investigate the problem, while Tariq and Seadawy [40] adopted the auxiliary equation method to ascertain exact solutions. In two different studies, Wang [41, 42] investigated the mBBM model using different methodologies to obtain various wave structures, including bright waves, bright-dark waves, bright-like waves, kinky-bright waves, and periodic waves. Furthermore, Khater et al. [34] utilized the modified Khater method, and Naeem et al. [43] introduced the Sardar-subequation method to uncover abundant exact solutions for mBBM equation. Additionally, Khater and Salama [44] applied the modified Khater method to determine its solutions. To date, the mBBM Eq (1) has not been explored by the rational -expansion method [16].

The goal of this research is to identify new analytical wave solutions, such as kink, anti-kink, anti-peakon, singular periodic, anti-bell and V-shaped soliton solutions of the mBBM model (3) through the rational -expansion method. Varying parameters can cause changes in the qualitative structure of the flow. This can result in the creation or destruction of fixed points, or changes in their stability. These changes in dynamics are known as bifurcations, and the specific parameter values at which they occur are called bifurcation points. Bifurcations are significant as they offer a model for understanding transitions and instabilities that occur when a control parameter is changed. Studying the movement of orbits in a phase plane can offer valuable insights into the dynamics of the system and aid in forecasting its future behavior. By examining both the behavior of the orbit in the phase plane and the corresponding solutions of the components with respect to the wave variable, one can gain a comprehensive understanding of the system’s dynamics and make predictions about its future behavior. Furthermore, we thoroughly investigate the stability of the model, analyze the parametric effects on the wave solutions, and conduct a comprehensive phase-plane analysis to gain profound insights into complex systems across diverse scientific domains.

This paper is constructed into the subsequent sections. We have introduced the rational -expansion method in Section 2 to find exact solitary wave solutions from NLPDE. In Section 3, we applied this method to solve the mBBM equation to obtain soliton solutions. Section 4 offers a graphical interpretation of the uncovered solutions to discuss their characteristics. Section 5 delves into an examination of the stability of the mBBM equation and compares our solutions with those obtained using other methods. Section 6 demonstrates the effects of different parameters on wave profiles. Finally, before concluding the discussion, Section 7 explores phase plane analysis.

2. Brief outline of the rational -expansion method

Consider the NLPDE featuring independent variables x, y, z, ⋯, t and the dependent variable u as(4)where P represents a polynomial in its variables, encompassing the highest-order nonlinear and derivative terms. The transformation u(x, t) = V(ζ), ζ = lx + ωt converts Eq (4) to the ordinary differential equation (ODE)(5)where prime denotes the differentiation with respect to ζ. According to the rational -expansion method, the solution of Eq (3) is assumed as follows:(6)where α_j, β_j (j = 1, 2, ⋯, n) are all unknown constants and G = G(ζ) is the solution of the ensuing ODE:(7)where ν and μ are any constants. The general solution of Eq (7) is(8)To estimate the unknown wave function u(x, t) = V(ζ), we follow the succeeding steps.

1. Step 1. Establishing homogeneous balance between the highest-order derivative of V and the highest-order nonlinear term in Eq (5) allows for the determination of the integer n.

2. Step 2. After determining n, the Eqs (6) and (7) are replaced into Eq (5) and then the Eq (5) becomes a polynomial of . Equating the coefficients of the same powers of , this polynomial yields a set of nonlinear algebraic equations involving α_j, β_j, ν, μ, k, and ω.

3. Step 3. Solving the derived system for α_j, β_j, k, and ω, and subsequently substituting the determined values of α_j, β_j, along with the general solution of (7) into Eq (6), yield the sought-after solutions for Eq (4).

Although numerous attempts have been made to discover different techniques for solving both integrable and non-integrable NLPDEs, developing a unified method for solving these equations is still difficult to achieve. The rational -expansion method excels in generating exact traveling wave solutions for NLPDEs, delivering solutions of high precision and accuracy. These precise solutions play a crucial role in elucidating the fundamental principles underlying physical phenomena. The method we are discussing offers advantages compared to other approaches. The primary advantage of the rational -expansion method compared to other methods lie in its provision of more general solutions featuring several free parameters. Through judicious parameter selection, these solutions can coincide with established solutions obtained through existing methods. Again, the proposed technique handles NLPDEs in a direct manner, requiring no initial/boundary conditions or initial trial function at the outset, unlike numerical methods. The capacity of the suggested approach to incorporate a larger number of arbitrary constants compared to other existing techniques renders it especially valuable in the realm of mathematical physics for the analysis of nonlinear wave phenomena.

Despite its numerous benefits, it is important to consider the limitations of the method. This method occasionally generates solutions that resemble well-known solutions obtainable through alternative techniques. Furthermore, depending on the efficiency of the computer, the rational -expansion approach may not guarantee the existence of solutions to the resulting algebraic equations when the balance number of the reduced ordinary differential equation exceeds two.

3. Analytical wave solutions of the mBBM equation

In this section the rational -expansion method is applied on the mBBM Eq (3) to find its analytical wave solutions. The wave variable transformations u(x, t) = V(ζ), ζ = lx + ωt converts the stated model (3) into the following ODE:(9)where l and ω represents the wave number and speed of the wave respectively. Balancing the term V³ with V′′ which are present in Eq (9), we get n = 1. Then the solution for this equation can be written as(10)where, the constants a₀, a₁, b₀, and b₁ are unspecified and will be determined subsequently.

Now using step 2 of the technique and replacing the values of Eqs (10) along with (7) in (9), we obtain the subsequent system of equations:(11a) (11b) (11c) (11d)By solving the above system using Maple, we obtain the following sets of solutions:

Set 1:(12)Set 2:(13)Set 3:(14)Set 4:(15)Substituting the provided estimates from (12)–(15) into the solution (10), we get(16)where .(17)where .(18)where .(19)where .

The traveling wave solutions for the mBBM model (3) are derived by employing the outcomes presented in (8) based on the solution (16) in the following manner:

When ν² − 4μ > 0, the hyperbolic trigonometric solution is obtained as(20)where A, B are integral constants and .

When ν² − 4μ < 0, the trigonometric solution is achieved as(21)where A, B are arbitrary constants and .

When ν² − 4μ = 0, the rational solution is gained as(22)where A, B are free constants and ζ = l(x − αt).

Incorporating the solutions organized in Eq (8) into Eq (17), we obtain subsequent soliton evaluations for the mBBM Eq (3):

For ν² − 4μ > 0, the hyperbolic trigonometric solution is obtained as(23)and for ν² − 4μ < 0, the trigonometric solution is achieved as(24)where A, B are random constants and .

By employing the solutions provided in (8) within the context of (18), we derive the subsequent progressive wave solutions for the mBBM Eq (3):

If ν² − 4μ > 0, the hyperbolic trigonometric solution is simplified as(25)If ν² − 4μ < 0, the trigonometric solution is calculated as(26)If ν² − 4μ = 0, the rational solution is simplified as(27)where A, B are subjective constants and .

Lastly, by utilizing the solutions compiled in (8) based on (18), we ascertain the subsequent exact wave solutions for the mBBM Eq (3):

Because of ν² − 4μ > 0, the hyperbolic function solution is ascertained as(28)Because of ν² − 4μ < 0, the trigonometric function solution is attained as(29)where A, B are any constants and .

4. Characteristics of the wave solutions

In this section, we delve into the attributes of certain solutions identified from the mBBM model. The rational -expansion method was utilized for the mBBM equation, yielding wave solutions characterized by various free parameters. This section explores the graphical representations and explanations of these solutions in the context of the mBBM equation. Various 3D, density and 2D plots illustrating the obtained solutions are generated using Matlab and Maple software, allowing for a comprehensive analysis of their nature. To emphasize this point, by assigning different values to the free parameters we have generated plots depicting the solutions (20), (21), (22), (23), (25), (27), (28) and (29). This exploration showcases various types of solitons, including kink, anti-kink, anti-peakon, singular periodic, anti-bell (dark) and V-shaped solitons. Different combinations of the free parameters utilized in the study may result in the emergence of diverse soliton solutions.

The solution (20) describes a soliton with anti-kink shape, where parameters α = 1, β = −1, l = 1, γ = −1, ν = 1, μ = −1, A = 0 and B = 1 are employed. This soliton is illustrated in Fig 1 for the limit −10 ≤ x, t ≤ 10. Fig 1(b) illustrates the contour plot of the solution (20), while Fig 1(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

The expression (21) represents a singular periodic soliton observed for parameters α = 1, β = 10, l = 1, γ = 1, ν = 4, μ = 5, A = 0 and B = 1 within −5 ≤ x, t ≤ 5, as shown in Fig 2. Fig 2(b) illustrates the contour plot of the solution (20), while Fig 2(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

The expression (23) describes the features of the kink soliton when α = 1, β = −0.5, l = 1, γ = −1, ν = 1, μ = −1, A = 1 and B = 0 over the interval −10 ≤ x, t ≤ 10, as exhibited in Fig 3. Fig 3(b) illustrates the contour plot of the solution (20), while Fig 3(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

The dark soliton, manifested in Fig 4, is expressed by the absolute value of the solution (25) with parameters α = 1, β = 1, l = 1, γ = 1, ν = 1, μ = −1, b₀ = 1, b₁ = 1, A = 1 and B = 0, confined within the boundaries −5 ≤ x, t ≤ 5. Fig 4(b) illustrates the contour plot of the solution (20), while Fig 4(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

In Fig 5, it is evident that the modulus of the solution (28) produces an anti-peakon soliton for the given parameters α = 1, β = 1, γ = 5, ν = 2, μ = −1, A = 1 and B = 0, within the specified range of −10 ≤ x, t ≤ 10. Fig 5(b) illustrates the contour plot of the solution (20), while Fig 5(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

The V-shaped soliton, as depicted in Fig 6, is obtained from the absolute value of the solution (29) under the conditions α = −0.01, β = −0.01, γ = 6, ν = 2, μ = 2, A = 1 and B = 0. This soliton is observed within the domain −2 ≤ x, t ≤ 2. Fig 6(b) illustrates the contour plot of the solution (20), while Fig 6(c) presents the 2D plots at different time levels.

[Figure omitted. See PDF.]

The analysis presented above demonstrates the effectiveness of the described method in generating accurate wave solutions for the mBBM model, encompassing diverse profiles. Among various wave structures, kink-type solitons are traveling wave solutions that transition from one phase to another by rising or descending. As they propagate to infinity, kink solitons eventually stabilize at a constant value. Additionally, kink-type solitons play a vital role in various fields such as nonlinear optics and condensed matter physics. They represent stable, localized transitions between distinct states within materials, offering valuable insights into defects in crystalline structures, phase transitions, and optical communication systems.

A periodic soliton is a robustly formed solitary wave capable of traveling long distances with consistent velocity and intensity. It maintains a clearly defined amplitude and phase, representing a highly stable solution. Such solitons are employed for transmitting signals across long distances.

Singular solitons are significant for understanding the mathematical modeling of physical problems, but they have limited applicability in practical physical scenarios. These are a form of solitary wave that has a singularity, which is typically an infinite discontinuity. Singularities can appear practically everywhere, and they are surprisingly common in the mathematical models that physicists employ to study the universe. Furthermore, singular solitons denote sudden, concentrated disturbances within physical systems, playing a crucial role in diverse fields such as nonlinear optics and fluid dynamics. They embody localized regions of intense energy or complex structures, providing insights into intricate phenomena, shockwaves, and compact information carriers within nonlinear dynamic wave equations [45].

Dark solitons are dark pulses that exist within the interior of a continuous wave, maintaining a distinct structure while preserving their propagation characteristics. They are pulses that propagate through various media such as liquids, gases, plasmas, and solids.

An anti-peakon soliton is characterized by a discontinuous first derivative at its lowest peak point. Similar to the kink soliton, an anti-peakon soliton maintains a constant value as time approaches infinity. It is also known as dark soliton in optical fiber studies.

V-shaped dark solitons offer soliton control in optical fiber studies. They suddenly change direction at the zero point, maintaining the original pulse duration and velocity. This behavior may have applications in the design of all optical switches. Therefore, the importance and essential role of different solitons become undeniable in elucidating nonlinear and intricate physical phenomena.

5. Comparison and stability analysis

5.1 Comparison

In this su-bsection, we have compared the solutions derived in this article with those documented in the existing literature for the mBBM equation. Many investigations in existing literature have explored the analysis of the mBBM equation employing various methodologies. Specifically, we contrasted the solutions obtained through the technique [32] and the Sardar-subequation method [43] with the outcomes presented in this work. The method and the Sardar-subequation method are used to analyze the NLPDE in mathematical physics, each of which possesses distinct characteristics and advantages.

1. i. Aslan [32] employed the -expansion method to discover exact and explicit solutions for (3). He had successfully derived one set of solutions for the mBBM model. In our approach, when b₀ = 1 and b₁ = 0, it reduces to the -expansion method. Notably, we have identified four distinct set of solutions using this modified method. Among these sets, the solutions in Set 1 perfectly align with all the solutions previously derived in [32]. Consequently, it can be asserted that the rational method exhibits greater efficacy compared to the conventional -expansion method. Table 1 presents a comparison between the solutions obtained by Aslan [32] and our acquired solutions for the mBBM equation. From the table, it is evident that the obtained solutions V₁₁(ζ), V₁₂(ζ), and V₁₃(ζ) are same as the solutions presented by Aslan [32]. Therefore, the additional solutions for the mBBM equation, derived using the aforementioned method, are all new and significant.

2. ii. The mBBM model (3), with α = 1, β = μ, and γ = −1, was investigated by Naeem et al. [43] using the Sardar-sub equation method to obtain diverse wave solutions. They identified exponential function solutions, which subsequently transformed into hyperbolic and trigonometric solutions. In contrast, our method yields hyperbolic, trigonometric, and rational function solutions. Moreover, our approach involves more parameters than theirs. Describing wave structures becomes more comfortable with solutions containing a greater number of parameters. Therefore, we assert that our method is more successful than theirs in analyzing the mBBM model. Table 2 presents a correlation between the solutions derived by Naeem et al. [43] and our obtained solutions for the mBBM model. It becomes evident that some of the acquired solutions align with those from previous studies under suitable configurations of arbitrary parameters, while others are novel and standard.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

5.2 Stability analysis

Numerous nonlinear phenomena exhibit instability, leading to the modulation of steady states due to the interplay between nonlinear and dispersive effects [46]. In this section, we explore the modulation instability of the mBBM Eq (3) through standard linear stability analysis [46, 47]. The linear stability begins with a given steady state solution along with its associated small perturbation, formulated as follows [47, 48]:(30)where δ << 1 denotes the magnitude of the perturbation, and it is evident that any constant value H₀ serves as a steady-state solution to the Eq (3). By introducing the perturbation δw(x, t), we aim to analyze the stability of the steady state under small disturbances. Substituting the solution (30) into (3) and then linearizing in δ, we have(31) Eq (31) is a linear constant coefficient partial differential equation with unknown w(x, t). Contemplate a solution to this equation in the following form [47]:(32)where k is the wave number of normalized form. Putting (32) in (31) and then solving for c, we obtain(33)It is evident from Eq (33) that any superposition of the solutions will exhibit a decaying behavior for negative values of c, indicating stable dispersion. The dispersion would be unstable otherwise.

6. Parametric effects on the wave profiles

The importance of the model parameters in the mBBM equation is pivotal in determining the behaviour and temporal evolution of solutions. By delving into the examination of the influence of parameters, we aim to offer valuable insights into the dynamics of waves. The model parameters, α, β, γ, substantially impact various aspects of the system, including stability, bifurcation phenomena (Figs 14–17), wave propagation, amplitude modulation (Figs 7–13), and the overall dynamics. A profound comprehension of the effect of the parameters are imperative for both theoretical analysis and practical applications. To visually depict the impact of the parameters, Figs 7–13 present illustrations of amplitude modulation on the solution curves, namely v₁₁(x, t). These figures serve as visual aids in understanding how parameters shapes the behavior of the solutions in the given context. The divergence of the nonlinear coefficient β towards either positive or negative infinity indicates a corresponding decrease in wave amplitude v to zero. Conversely, as the linear coefficients (α, γ) approach positive or negative infinity, they lead to an augmentation in wave amplitude. Specifically, a positive value for β results in the formation of anti-kink waves, while a negative value for β yields in the generation of kink waves. Analogously, the influence of model parameters on alternative solutions can be demonstrated in a similar manner.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

7. Phase plane analysis

Varying parameters can cause changes in the qualitative structure of the flow. This can result in the creation or destruction of fixed points, or changes in their stability. Examining the movement of orbits in a phase plane can be highly informative in understanding the dynamics of the system and can assist in predicting its future behaviour. By analysing both the behaviour of the orbits in the phase plane and the corresponding solutions of the components concerning the wave variable, one can gain a comprehensive understanding of the dynamics of the system and make predictions about its future behaviour. We initiate the phase plane analysis [49] of the mBBM model, following the assumption that W = X and . Then we can express Eq (9) in the form of a first-order dynamical system as follows:(34)

The Hamiltonian associated with the planar system (34) is(35)

If ω = −lα, the system (34) possesses a single equilibrium point at (X*, Y*) = (0, 0). However, if ω ≠ −lα, the system exhibits three equilibrium points: (0, 0), and . It is to be noted that lβ ≠ 0.

The Jacobian matrix is defined asThe corresponding charactersitic equation of the matrix J iswhere trace(J) = 0 and . To examine the stability [49, 50] of the equilibrium points, we consider the following notations:

Case 1: Stability of (0, 0)

For this case we have p = 0, and . If , i.e., q < 0, then the fixed point (0, 0) is a saddle point that is unstable in nature. If , i.e., Δ < 0 and q > 0, then the equilibrium is a stable center. Hence, the variations in the parameter values have the potential to shift the stability of the equilibrium point (0, 0), transforming it from an unstable saddle to a stable center (Figs 14–17).

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Case 2: Stability of

In this case we have p = 0, and . If , then we have q < 0, indicating that the stationary points are unstable saddles. If then we have Δ < 0 and q > 0 which means that the equilibria are stable centers. Therefore, the change in the parameter values have the potential to shift the stability of the equilibrium points , transforming them from unstable saddle to stable center (Figs 14–15).

8. Conclusion

In this study, we delved into the investigation of the mBBM equation utilizing the rational -expansion method. Our examination led to the discovery of a rich array of soliton solutions, including kink, anti-kink, anti-peakon, singular kink, singular periodic, singular bell, dark, and V-shaped solitons. These diverse solutions serve as windows into the intricate physical behaviors inherent in the governing model, offering valuable insights for understanding its dynamics. To ensure the credibility of our findings, we rigorously cross-verified the obtained solutions by substituting them into the original differential equation, leveraging mathematical software, such as Maple. This meticulous validation process underscores the robustness of our results. Moreover, we systematically explored the influence of various parameters on the wave profile of the governing equation through the generation of 2D and 3D plots (please see Figs 7–13). Notably, our simulations revealed intriguing patterns of the wave profile in the behavior of the nonlinear wave coefficient β and the linear coefficients (α, γ) as they approach positive or negative infinity. Specifically, we observed a consequential reduction in wave amplitude v to zero with increasing asymptotic values of β, while the linear coefficients exhibited an amplification in wave amplitude under similar conditions (please see Figs 7–13). These insights deepen our understanding of the underlying mechanisms governing the mBBM model and its response to parameter variations.

A rigorous validation process using mathematical software Maple, coupled with the comprehensive analysis of parameter influences on the wave profile, enhances the credibility and novelty of our results. Additionally, our exploration of stability considerations and phase plane analysis adds another layer of original contribution to the study. By elucidating the influence of parameter variations on the dynamics of soliton solutions, we uncover novel insights into the behavior of the mBBM model, further distinguishing our work from existing research.

In essence, our study contributes to the comprehensive understanding of the mBBM equation and its soliton solutions, elucidating the intricate interplay between parameters and wave dynamics. These findings hold implications for diverse fields where nonlinear wave phenomena are encountered, offering avenues for further exploration and application.

Citation: Sagib M, Hossain MA, Saha BK, Khan K (2024) On traveling wave solutions with stability and phase plane analysis for the modified Benjamin-Bona-Mahony equation. PLoS ONE 19(7): e0306196. https://doi.org/10.1371/journal.pone.0306196

About the Authors:

Md. Sagib

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Mathematics, Hajee Mohammad Danesh Science and Technology University, Dinajpur, Bangladesh

Md. Aslam Hossain

Roles: Conceptualization, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing

Affiliation: Department of Mathematics, Pabna University of Science and Technology, Pabna, Bangladesh

Bijan Krishna Saha

Roles: Data curation, Formal analysis, Methodology, Validation, Writing – original draft, Writing – review & editing

Affiliation: Department of Mathematics, University of Barishal, Barishal, Bangladesh

Kamruzzaman Khan

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliations: Department of Mathematics, Pabna University of Science and Technology, Pabna, Bangladesh, School of Science and Technology, University of New England, Armidale, NSW, Australia

ORICD: https://orcid.org/0000-0002-4531-288X

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Seadawy AR, Ali A, Raddadi MH. Exact and solitary wave solutions of conformable time fractional clannish rando walkers parabolic and Ablowitz-Kaup-Newell-Segur equation via modified mathematical methods. Results in Physics 2021;26:104374.

2. Kumar S, Ma WX, Kumar A. Lie symmetries, optimal system and group-invariant solutions of the (3+ 1)-dimensional generalized KP equation. Chinese J. Phys 2021;69:1–23.

3. Khater MMA, Elagen SK, Mousa AA, Shorbagy AAE, Alfalqi SH, Alzaidi JF, Lu D. Sub-10 -fs-Pulse propagation between analytical and numerical investigation. Results in Physics 2021;25:104133.

4. Unal M, Ekici M. The double expansion method and its applications for some nonlinear partial differential equaitons. J. Institute of Sci. and Tech 2021;11:599–608.

5. Alquran M. New symmetric bidirectional progressive surface-wave solutions to a generalized fourth-order nonlinear partial differential equation involving second-order time-derivative. J Ocean Eng Sci 2022.

6. Hosseini K, Salahshour S, Mirzazadeh M, Ahmadian A, Baleanu D, Khoshrang A. The (2+1)-dimensional Heisenberg ferromagnetic spin chain equation: its solitons and Jacobi elliptic function solutions. Eur Phys J Plus 2021;136:1–9.

7. Jaradat I, Alquran M. A variety of physical structures to the generalized equalwidth equation derived from wazwaz-benjamin-bona-mahony model. J Ocean Eng Sci 2022;7(3):244–7.

8. Khater MM. Abundant and accurate computational wave structures of the nonlinear fractional biological population model. International Journal of Modern Physics B. 2023 Jul 20;37(18):2350176.

9. Malik S, Kumar S, Das A. A (2+ 1)-dimensional combined KdV-mKdV equation: integrability, stability analysis and soliton solutions. Nonlinear Dynam 2022;1–13.

10. He XJ, Lü X. M-lump solution, soliton solution and rational solution to a (3+1)-dimensional nonlinear model. Math Comput Simulat 2022;197:327–40.

11. Yu Z-B, Zhu C, Zhao J-S, Zou L. Inverse scattering transform of the general three-component nonlinear schrödinger equation and its multisoliton solutions. Appl Math Lett 2022;128:107874.

12. Khater MM. Novel computational simulation of the propagation of pulses in optical fibers regarding the dispersion effect. International Journal of Modern Physics B. 2023 Apr 10;37(09):2350083.

13. Yin YH, Lü X, Ma WX. Backlund transformation, exact solutions and diverse interaction phenomena to a (3+1)-dimensional nonlinear evolution equation. Nonlin Dyn 2022;108:4181–94.

14. Islam MT, Akbar MA, Ahmad H, Ilhan OA, Gepreel KA. Diverse and novel soliton structures of coupled nonlinear schrodinger type equations through two competent techniques. Modern Physics Letters B 2022;36(11):2250004.

15. Ayad MS, Khaled AEA, Mohamed OME, Ali MA, Abdullah MHA. The -expansion method for solving a nonlinear PDE describing the nonlinear low pass electrical lines. J. Taiban Univ. Sci 2019;13:63–70.

16. Akbar MA, Ali NHM, Islam MT. Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics. AIMS Mathematics 2019;4:397–411.

17. Mohanty SK, Kravchenko OV, Dev AN. Exact traveling wave solutions of the Schamel Burgers’ equation by using generalized-improved and generalized (G′/G)-expansion methods. Results Phys 2020;33:105124.

18. Mohanty SK, Kumar S, Dev AN, Deka MK, Churikov DV, Kravchenko OV. An efficient technique of (G′/G)-expansion method for modified KdV and Burgers equations with variable coefficients. Results Phys 2022;37:105504.

19. Khater MM. Hybrid accurate simulations for constructing some novel analytical and numerical solutions of three-order GNLS equation. International Journal of Geometric Methods in Modern Physics. 2023 Aug 6;20(09):2350159.

20. Khater MM. In solid physics equations, accurate and novel soliton wave structures for heating a single crystal of sodium fluoride. International Journal of Modern Physics B. 2023 Mar 20;37(07):2350068.

21. Khater MM. Analyzing pulse behavior in optical fiber: Novel solitary wave solutions of the perturbed Chen–Lee–Liu equation. Modern Physics Letters B. 2023 Dec 10;37(34):2350177.

22. Khater MM. Computational simulations of propagation of a tsunami wave across the ocean. Chaos, Solitons & Fractals. 2023 Sep 1;174:113806.

23. Khater MM. Soliton propagation under diffusive and nonlinear effects in physical systems;(1+ 1)–dimensional MNW integrable equation. Physics Letters A. 2023 Aug 28;480:128945.

24. Khater MM. Horizontal stratification of fluids and the behavior of long waves. The European Physical Journal Plus. 2023 Aug 15;138(8):715.

25. Khater MM. Advancements in computational techniques for precise solitary wave solutions in the (1+ 1)-dimensional Mikhailov-Novikov-Wang equation. International Journal of Theoretical Physics. 2023 Jul 14;62(7):152.

26. Wang KJ. Dynamics of breather, multi-wave, interaction and other wave solutions to the new (3+ 1)-dimensional integrable fourth-order equation for shallow water waves. International Journal of Numerical Methods for Heat & Fluid Flow. 2023 Oct 31;33(11):3734–47.

27. Chen SJ, Lu X. Observation of resonant solitons and associated integrable properties for nonlinear waves. Chaos Solitons Fract 2022;163:112543.

28. Wang KJ, Shi F. Multi-soliton solutions and soliton molecules of the (2+ 1)-dimensional Boiti-Leon-Manna-Pempinelli equation for the incompressible fluid. Europhysics Letters. 2024.

29. Wang KJ. Soliton molecules, Y-type soliton and complex multiple soliton solutions to the extended (3+ 1)-dimensional Jimbo-Miwa equation. Physica Scripta. 2024 Jan 2;99(1):015254.

30. Benjamin TB, Bona JL, Mahony JJ. Model equations for long waves in nonlinear dispersive systems. Philos. Trans. Roy. Soc. London Ser. A 1972;272:47–48.

31. Abbasbandy S, Shirzadi A. The first integral method for modified Benjamin–Bona–Mahony equation. Communications in Nonlinear Science and Numerical Simulation 2010; 15(7):1759–1764.

32. Aslan I. Exact and explicit solutions to some nonlinear evolution equations by utilizing the (G′/G)-expansion method. Applied Mathematics and Computation 2009; 215(2):857–863.

33. Naher H, Abdullah FA, Bekir A. Some new traveling wave solutions of the modified Benjamin-Bona-Mahony equation via the improved (G′/G)-expansion method. New Trends in Mathematical Sciences 2015;3(1):78–89.

34. Khater MM, Nofal TA, Abu-Zinadah H, Lotayif MS, Lu D. Novel computational and accurate numerical solutions of the modified Benjamin–Bona–Mahony (BBM) equation arising in the optical illusions field. Alexandria Engineering Journal 2021;60(1):1797–1806.

35. Shapiro AG, Hamburger K. When Is an Illusion an Illusion? An Examination of Contrast Information in Grouping and Grid Phenomena. In A New Look at New Realism 2017; 81–102. Routledge.

36. Sun ED, Dekel R. ImageNet-trained deep neural networks exhibit illusion-like response to the Scintillating grid. Journal of Vision 2021;21(11):15–15. pmid:34677575

37. Khan K, Akbar MA. Application of exp-expansion method to find the exact solutions of modified Benjamin-Bona-Mahony equation. World Applied Sciences Journal 2013;24(10):1373–1377.

38. Noor MA, Noor KI, Waheed A, Al-Said EA. Some new solitonary solutions of the modified Benjamin–Bona–Mahony equation. Computers & Mathematics with Applications. 2011 Aug 1;62(4):2126–31.

39. Song M. Nonlinear wave solutions and their relations for the modified Benjamin–Bona–Mahony equation. Nonlinear Dyn 2015;80:431–446.

40. Tariq KU, Seadawy AR. On the soliton solutions to the modified Benjamin-Bona-Mahony and coupled Drinfel’d-SokolovWilson models and its applications. J. King Saud Univ. Sci. 2020;32:156–162.

41. Wang KJ. Diverse wave structures to the modified Benjamin–Bona–Mahony equation in the optical illusions field. Modern Physics Letters B. 2023 Apr 20;37(11):2350012.

42. Wang KJ. Variational principle and diverse wave structures of the modified Benjamin-Bona-Mahony equation arising in the optical illusions field. Axioms. 2022 Aug 31;11(9):445.

43. Ullah N, Asjad MI, Hussanan A, Akgül A, Alharbi WR, Algarni H, et al. Novel waves structures for two nonlinear partial differential equations arising in the nonlinear optics via Sardar-subequation method. Alexandria Engineering Journal 2023;71:105–113.

44. Khater MMA, Salama SA. Semi-analytical and numerical simulations of the modified Benjamin-Bona-Mahony model. J. Ocean. Eng. Sci. 2022;7:264–271.

45. Islam MT, Sarkar S, Alsaud H, Inc M. Distinct wave profiles relating to a coupled of Schrödinger equations depicting the modes in optics. Optical and Quantum Electronics. 2024 Mar;56(3):492.

46. Agrawal GP. Nonlinear fiber optics. In Nonlinear Science at the Dawn of the 21st Century 2000 Dec 21 (pp. 195–211). Berlin, Heidelberg: Springer Berlin Heidelberg.

47. Abdullahi Y, Isa AA, Dumitru B. Soliton solution and stability analysis for some conformable nonlinear partial differential equations in mathematical physics. Optical Quantum Electron 2018;50:190.

48. Islam SR, Ahmad H, Khan K, Wang H, Akbar MA, Awwad FA, et al. Stability analysis, phase plane analysis, and isolated soliton solution to the LGH equation in mathematical physics. Open Physics. 2023 Oct 30;21(1):20230104.

49. Barnes B, Fulford GR. Mathematical modelling with case studies: a differential equations approach using Maple and MATLAB. CRC Press 2011;25.

50. Islam SR, Khan K, Akbar MA. Optical soliton solutions, bifurcation and stability analysis of the Chen-Lee-Liu model. Results in Physics 2023; 106620.