Acceso abierto

Dispersive modified Benjamin-Bona-Mahony and Kudryashov-Sinelshchikov equations: non-topological, topological, and rogue wave solitons


Cite

Introduction

Nonlinear differential equations arise in modeling several physical phenomena. These equations help handle the dynamics of the specific model [1,2,3]. Particularly interesting circumstances are modeled using nonlinear partial differential equations (NPDEs). Nonlinear wave phenomena are studied in a wide range of scientific and engineering fields, such as optical fibers, computational fluid dynamics, plasma physics, solid-state physics chemical dynamics, biological, and chemical-physical science, geochemistry, and shallow waves [4,5,6]. Nonlinear wave processes involving dissipation, dispersion, responses, convection, and propagation are crucial in understanding nonlinear wave equations [7,8,9]. Many approaches have been applied in recent years to examine nonlinear problems. Some of these are the Jacobi elliptic function method [10,11,12], the sub-equation approach [13], the exponential rational function method [14], new extended direct algebraic method [15] as well as symmetry approaches [16,17,18,19,20,21], and others [22,23,24,25,26]. Solitons, which are nonlinear confined solitary waves that preserve their shape while traveling at a constant speed, are a part of numerous nonlinear physical systems. The theory of soliton propagation has been utilized for many important oceanic phenomena, including that of the dispersion of nonlinear shallow as well as deep sea waves. In addition to gravity waves in a related domain and nonlinear wave motion rescaling, the wave motion includes wave propagation in dispersive media such as gas bubbles, liquid flow in lakes, rivers, the ocean, and fluid flow in elastic tubes. This type of evolution equation may have a substantial impact on studies of fluid flow and ocean wave motion. It is a common practice to use nonlinear evolution equations to simulate ocean and atmospheric dynamics. Ocean engineering is concerned with the size and frequency of the waves in the ocean as well as their temperature and density (such as tsunamis and wind waves). The nonlinear evolution equations [27,28,29,30] can be used to describe nonlinear wave propagation in seas. One of the main manifestations of nonlinear wave phenomena is the spread of nonlinear waves in a dynamical system, which has attracted increasing attention. In this study, the rogue wave, topological, and non-topological soliton solutions were derived utilizing the solitary wave ansatz method. We categorically acknowledge the reliability of this approach, which has recently been successfully applied. The improved simple equation method was utilized by Khan et al. [31] to identify accurate traveling wave solutions for the coupled Klein-Gordon equations and the NDMBBM equation. To create novel, accurate solutions to the NDMBBM problem and the mKdV-Burgers equation, Filiz and Arzu [32] used the consistent Riccati expansion approach. Kumar et al. [33] used the sinh-Gordon function method to validate the stability analysis and discover numerical and fresh closed-form trigonometric function solutions of the NDMBBM equation. They carried out numerical simulations of the discovered solutions using the proper parameters. The exact and numerical approximations of the nonlinear NDMBBM equation were also tackled using the finite forward difference method. To evaluate the stability of the numerical technique, they used the Fourier-Von Neumann analysis. In a table format, they presented the outcomes for both the numerical and exact answers.

The rest of the paper is organized as follows. Section 2 presents the ansatz approach to NDMBBM equation. Section 3 describe the ansatz approach to KS equation. The models, which includes the development of topological, non-topological, and rogue wave soliton solutions are also introduced. Section 4 introduces the physical interpretation of the results. Finally, Section 5 reports the novelties of this paper through the discussion and conclusions.

Ansatz approach to NDMBBM equation

The nonlinear NDMBBM equation is defined as follows Yt+YxαY2Yx+Yxxx=0, {Y_t} + {Y_x} - \alpha {Y^2}{Y_x} + {Y_{xxx}} = 0, where α is a nonzero real constant. The nonlinear NDMBBM equation was originally formulated to model surface wave propagation in nonlinear dispersive media, but it also describes hydromagnetic waves in plasma, acoustic waves in unique form crystals, and acoustic gravitational waves in compressible fluids [34].

Solitons and their topological characteristics provide an intriguing subject of study. Dark solitons, labeled as topological solitons, and bright solitons, termed non-topological solitons, showcase unique traits. This topological classification enhances our grasp of these phenomena, finding applications in diverse physical systems. The interplay between topology and solitons opens exciting avenues for investigating and controlling complex structures across various scientific domains.

Non-topological (bright) soliton solutions

The following non-topological wave ansatz [35] has been used in this section to find non-topological solutions for the equation (1), Y(x,t)=Ecoshp(κ), Y(x,t) = \frac{E}{{\mathop {\cosh }\nolimits^p (\kappa )}}, in which κ = C(xvt). In the aforementioned equations, the variables E and C are referred to as unconstrained parameters. E stands for amplitude of soliton, C is the inverse width, and v denotes the speed. Later, the homogeneous balance principle is used to determine how to evaluate the exponent p. From (2), it is possible to follow Yn(x,t)=Encoshpn(κ), {Y^n}(x,t) = \frac{{{E^n}}}{{\mathop {\cosh }\nolimits^{pn} (\kappa )}}, so that (Y(x,t))t=ECpvtanh(κ)coshp(κ), {(Y(x,t))_t} = ECpv\frac{{\tanh (\kappa )}}{{\mathop {\cosh }\nolimits^p (\kappa )}}, (Y(x,t))x=ECptanh(κ)coshp(κ), {(Y(x,t))_x} = - ECp\frac{{\tanh (\kappa )}}{{\mathop {\cosh }\nolimits^p (\kappa )}}, and (Y(x,t))xxx=EC3ptanh(κ)p2coshp(κ)+(p+1)(p+2)coshp+2(κ) {(Y(x,t))_{xxx}} = E{C^3}p\tanh (\kappa )\left( { - \frac{{{p^2}}}{{\mathop {\cosh }\nolimits^p (\kappa )}} + \frac{{(p + 1)(p + 2)}}{{\mathop {\cosh }\nolimits^{p + 2} (\kappa )}}} \right) \cdot Eq (1) undergoes ECpvtanh(κ)coshp(κ)ECptanh(κ)coshp(κ)+α(Ecoshp(κ))2ECptanh(κ)coshp(κ)+EC3ptanh(κ)p2coshp(κ)+(p+1)(p+2)coshp+2(κ)=0. \begin{array}{*{20}{r}}{ECpv\frac{{\tanh (\kappa )}}{{\mathop {\cosh }\nolimits^p (\kappa )}} - ECp\frac{{\tanh (\kappa )}}{{\mathop {\cosh }\nolimits^p (\kappa )}} + \alpha {{(\frac{E}{{\mathop {\cosh }\nolimits^p (\kappa )}})}^2}ECp\frac{{\tanh (\kappa )}}{{\mathop {\cosh }\nolimits^p (\kappa )}} + }\\{E{C^3}p\tanh (\kappa )\left( { - \frac{{{p^2}}}{{\mathop {\cosh }\nolimits^p (\kappa )}} + \frac{{(p + 1)(p + 2)}}{{\mathop {\cosh }\nolimits^{p + 2} (\kappa )}}} \right) = 0.}\end{array} Through the division of ECp tanh(κ) and subsequent manipulation, we reach the following result v1C2p2coshp(κ)+(αE2cosh3p(κ))+C2(p+1)(p+2)coshp+2(κ)=0. \frac{{v - 1 - {C^2}{p^2}}}{{\mathop {\cosh }\nolimits^p (\kappa )}} + (\frac{{\alpha {E^2}}}{{\mathop {\cosh }\nolimits^{3p} (\kappa )}}) + \frac{{{C^2}(p + 1)(p + 2)}}{{\mathop {\cosh }\nolimits^{p + 2} (\kappa )}} = 0. By applying the homogeneous balancing principle, we may assume that 3p = p + 2 and C2(p + 1)(p + 2) ≠ 0.

This results in p = 1. The following system is obtained from (7) as v1C2p2=0,αE2+C2(p+1)(p+2)=0. \begin{array}{*{20}{c}}{v - 1 - {C^2}{p^2} = 0,}\\{\alpha {E^2} + {C^2}(p + 1)(p + 2) = 0.}\end{array} We solve the aforementioned system of equations for p = 1 and obtain v1C2=0,αE2+6C2=0. v - 1 - {C^2} = 0,\;\;\;\alpha {E^2} + 6{C^2} = 0. For α < 0, we can compute the velocity of the soliton, v=1C2. v = 1 - {C^2}. Subsequently, we ascertain the amplitude (E) of the soliton. E=C6α E = C\sqrt { - \frac{6}{\alpha }} \cdot Hence, the NDMBBM equation (1) possesses a bright soliton solution Y(x,t)=C6αcosh(C(x(1C2)t)) Y(x,t) = \frac{{C\sqrt { - \frac{6}{\alpha }} }}{{\cosh (C(x - (1 - {C^2})t))}} \cdot

Rogue wave solutions

In this section, the rogue wave solutions for (1) have been generated using the wave ansatz [35]. Y(x,t)=Esinhp(κ), Y(x,t) = \frac{E}{{\mathop {\sinh }\nolimits^p (\kappa )}}, in which κ = C(xvt). In the given equations, the variables E and C are denoted as unconstrained parameters. Here, E represents the wave’s amplitude, C signifies their inverse width, and v indicates their speed. Subsequently, the homogeneous balance principle is employed to establish the evaluation of the exponent p. The progression can be traced from (13) (Y(x,t))t=ECpvcoth(κ)sinhp(κ), {(Y(x,t))_t} = ECpv\frac{{\coth (\kappa )}}{{\mathop {\sinh }\nolimits^p (\kappa )}}, (Y(x,t))x=ECpcoth(κ)sinhp(κ), {(Y(x,t))_x} = - ECp\frac{{\coth (\kappa )}}{{\mathop {\sinh }\nolimits^p (\kappa )}}, and (Y(x,t))xxx=EC3pcoth(κ)p2sinhp(κ)+(p+1)(p+2)sinhp+2(κ). {(Y(x,t))_{xxx}} = E{C^3}p\coth (\kappa )\left( { - \frac{{{p^2}}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + \frac{{(p + 1)(p + 2)}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}}} \right). In consideration of this, equation (1) can be expressed as ECpvcoth(κ)sinhp(κ)ECpcoth(κ)sinhp(κ)+α(Esinhp(κ))2ECpcoth(κ)sinhp(κ)+EC3pcoth(κ)p2sinhp(κ)+(p+1)(p+2)sinhp+2(κ)=0. \begin{array}{*{20}{r}}{ECpv\frac{{\coth (\kappa )}}{{\mathop {\sinh }\nolimits^p (\kappa )}} - ECp\frac{{\coth (\kappa )}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + \alpha {{(\frac{E}{{\mathop {\sinh }\nolimits^p (\kappa )}})}^2}ECp\frac{{\coth (\kappa )}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + }\\{E{C^3}p\coth (\kappa )\left( { - \frac{{{p^2}}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + \frac{{(p + 1)(p + 2)}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}}} \right) = 0.}\end{array} Upon dividing by ECp tanh(κ), we obtain vsinhp(κ)1sinhp(κ)+(Esinhp(κ))2αsinhp(κ)+C2p2sinhp(κ)+(p+1)(p+2)sinhp+2(κ)=0. \begin{array}{*{20}{r}}{\frac{v}{{\mathop {\sinh }\nolimits^p (\kappa )}} - \frac{1}{{\mathop {\sinh }\nolimits^p (\kappa )}} + {{(\frac{E}{{\mathop {\sinh }\nolimits^p (\kappa )}})}^2}\frac{\alpha }{{\mathop {\sinh }\nolimits^p (\kappa )}}}\\{ + {C^2}\left( { - \frac{{{p^2}}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + \frac{{(p + 1)(p + 2)}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}}} \right) = 0.}\end{array} v1C2p2sinhp(κ)+αE2sinh3p(κ)+C2(p+1)(p+2)sinhp+2(κ)=0. \frac{{v - 1 - {C^2}{p^2}}}{{\mathop {\sinh }\nolimits^p (\kappa )}} + \frac{{\alpha {E^2}}}{{\mathop {\sinh }\nolimits^{3p} (\kappa )}} + \frac{{{C^2}(p + 1)(p + 2)}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}} = 0. Following the homogeneous balancing principle, we assume that αE2 ≠ 0 and C2 ≠ (p + 1)(p + 2) ≠ 0 possess equal powers. Consequently, we deduce that p = 1. The subsequent system is derived from (19) v1C2p2=0,αE2+C2(p+1)(p+2)=0. \begin{array}{*{20}{c}}{v - 1 - {C^2}{p^2} = 0,}\\{\alpha {E^2} + {C^2}(p + 1)(p + 2) = 0.}\end{array} Upon solving the aforementioned set of equations for p = 1, we reach v1C2=0,αE2+6C2=0. v - 1 - {C^2} = 0,\;\;\;\alpha {E^2} + 6{C^2} = 0. For α < 0, we determine the speed of the wave v=1C2. v = 1 - {C^2}. Subsequently, we ascertain the amplitude E of the soliton. E=C6α E = C\sqrt { - \frac{6}{\alpha }} \cdot Therefore, the NDMBBM equation (1) features a rogue wave solution represented by Y(x,t)=C6αsinh(C(x(1C2)t)) Y(x,t) = \frac{{C\sqrt { - \frac{6}{\alpha }} }}{{\sinh (C(x - (1 - {C^2})t))}} \cdot

Topological (dark) soliton solution

In this section, topological solutions to equation (1) have been derived by employing the topological wave ansatz [35] Y(x,t)=Etanhy(κ), Y(x,t) = E\, \mathop {\tanh }\nolimits^y (\kappa ), in which κ = D(xvt). In the given equations, the variables E and D are denoted as unconstrained parameters. Here, E represents the amplitude of solitons, D signifies their inverse width, and v indicates their velocity. Subsequently, the homogeneous balance principle is employed to establish the evaluation of the exponent y. The progression can be traced from (25) (Y(x,t))t=EDyv(tanhy1(κ)tanhy+1(κ)), {(Y(x,t))_t} = - EDyv(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa )), (Y(x,t))x=EDy(tanhy1(κ)tanhy+1(κ)), {(Y(x,t))_x} = EDy(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa )), and (Y(x,t))xxx=ED3y((y1)(y2)tanhy3(κ)(3y23y+2)tanhy1(κ)+(3y2+3y+2)tanhy+1(κ)(y+1)(y+2)tanhy3(κ)). \begin{array}{*{20}{r}}{{{(Y(x,t))}_{xxx}} = E{D^3}y((y - 1)(y - 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa ) - (3{y^2} - 3y + 2)\mathop {\tanh }\nolimits^{y - 1} (\kappa )}\\{ + \;(3{y^2} + 3y + 2)\mathop {\tanh }\nolimits^{y + 1} (\kappa ) - (y + 1)(y + 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa )).}\end{array} Eq (1) consequently transforms into EDyv(tanhy1(κ)tanhy+1(κ))+EDy(tanhy1(κ)tanhy+1(κ))αEDy(Etanhy(κ))2(tanhy1(κ)tanhy+1(κ))+ED3y(y1)(y2)tanhy3(κ)(3y23y+2)tanhy1(κ)+(3y2+3y+2)tanhy+1(κ)(y+1)(y+2)tanhy3(κ)=0. \begin{array}{*{20}{r}}{ - EDyv(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa )) + EDy(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ - \alpha EDy{{(E\, \mathop {\tanh }\nolimits^y (\kappa ))}^2}(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ + E{D^3}y\left( {(y - 1)(y - 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa ) - (3{y^2} - 3y + 2)\mathop {\tanh }\nolimits^{y - 1} (\kappa )} \right.}\\{\left. { + (3{y^2} + 3y + 2)\mathop {\tanh }\nolimits^{y + 1} (\kappa ) - (y + 1)(y + 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa )} \right) = 0.}\end{array} Upon division by EDy, we get v(tanhy1(κ)tanhy+1(κ))+(tanhy1(κ)tanhy+1(κ))α(Etanhy(κ))2(tanhy1(κ)tanhy+1(κ))+D2(y1)(y2)tanhy3(κ)(3y23y+2)tanhy1(κ)+(3y2+3y+2)tanhy+1(κ)(y+1)(y+2)tanhy3(κ)=0. \begin{array}{*{20}{r}}{ - v(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa )) + (\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ - \alpha {{(E\, \mathop {\tanh }\nolimits^y (\kappa ))}^2}(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ + {D^2}\left( {(y - 1)(y - 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa ) - (3{y^2} - 3y + 2)\mathop {\tanh }\nolimits^{y - 1} (\kappa )} \right.}\\{\left. { + (3{y^2} + 3y + 2)\mathop {\tanh }\nolimits^{y + 1} (\kappa ) - (y + 1)(y + 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa )} \right) = 0.}\end{array} We further simplify it in the form v(tanhy1tanhy+1(κ))+v(tanhy1(κ)tanhy+1(κ))α(Etanhy(κ))2(tanhy1(ξ)tanhy+1(κ))+D2(y1)(y2)tanhy3(κ)(3y23y+2)tanhy1(κ)+(3y2+3y+2)tanhy+1(κ)(y+1)(y+2)tanhy3(κ)=0. \begin{array}{*{20}{r}}{ - v(\mathop {\tanh }\nolimits^{y - 1} - \mathop {\tanh }\nolimits^{y + 1} (\kappa )) + v(\mathop {\tanh }\nolimits^{y - 1} (\kappa ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ - \alpha {{(E\, \mathop {\tanh }\nolimits^y (\kappa ))}^2}(\mathop {\tanh }\nolimits^{y - 1} (\xi ) - \mathop {\tanh }\nolimits^{y + 1} (\kappa ))}\\{ + {D^2}\left( {(y - 1)(y - 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa ) - (3{y^2} - 3y + 2)\mathop {\tanh }\nolimits^{y - 1} (\kappa )} \right.}\\{\left. { + (3{y^2} + 3y + 2)\mathop {\tanh }\nolimits^{y + 1} (\kappa ) - (y + 1)(y + 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa )} \right) = 0.}\end{array} By separating different powers of tanh, we obtain tanhy1(κ)(1vD2(3y23y+2))+tanhy+1(κ)(v1+D2(3y2+3y+2))αE2tanh3y1(κ)+αE2tanh3y+1(κ)+D2(y1)(y2)tanhy3(κ)D2(y+1)(y+2)tanhy+3(κ)=0. \begin{array}{*{20}{r}}{\mathop {\tanh }\nolimits^{y - 1} (\kappa )(1 - v - {D^2}(3{y^2} - 3y + 2)) + \mathop {\tanh }\nolimits^{y + 1} (\kappa )(v - 1 + {D^2}(3{y^2} + 3y + 2))}\\{ - \alpha {E^2}\mathop {\tanh }\nolimits^{3y - 1} (\kappa ) + \alpha {E^2}\mathop {\tanh }\nolimits^{3y + 1} (\kappa ) + {D^2}(y - 1)(y - 2)\mathop {\tanh }\nolimits^{y - 3} (\kappa )}\\{ - {D^2}(y + 1)(y + 2)\mathop {\tanh }\nolimits^{y + 3} (\kappa ) = 0.}\end{array} Given that αE2 ≠ 0, D2(y + 1)(y + 2) ≠ 0, and D2(y − 1)(y − 2) ≠ 0, we assume their powers to be equal. By setting 3y − 1 = y − 3, we find y = −1. However, since y cannot be negative, we explore alternatives. Considering 3y + 1 = y + 3, we arrive at y = 1. Equation (32) for y = 1 implies (1v2D2)+tanh2(κ)(v1+8D2αE2)+tanh4(κ)(αE26D2)=0. (1 - v - 2{D^2}) + \mathop {\tanh }\nolimits^2 (\kappa )(v - 1 + 8{D^2} - \alpha {E^2}) + \mathop {\tanh }\nolimits^4 (\kappa )(\alpha {E^2} - 6{D^2}) = 0. The system derived from (33) is as follows 1v2D2=0,v1+8D2E2α=0,αE26D2=0. \begin{array}{*{20}{r}}{1 - v - 2{D^2} = 0,}\\{v - 1 + 8{D^2} - {E^2}\alpha = 0,}\\{\alpha {E^2} - 6{D^2} = 0.}\end{array} The speed of the soliton, for arbitrary values of D and α > 0, is v=12D2. v = 1 - 2{D^2}. Subsequently, we determine the amplitude E of the soliton as E=±D6α E = \pm D\sqrt {\frac{6}{\alpha }} \cdot Therefore, the NDMBBM equation (1) features the dark soliton solution given by Y1(x,t)=D6αtanh(D(x(12D2)t)), {Y_1}(x,t) = D\sqrt {\frac{6}{\alpha }} \tanh (D(x - (1 - 2{D^2})t)), and Y2(x,t)=D6αtanh(D(x(12D2)t)). {Y_2}(x,t) = - D\sqrt {\frac{6}{\alpha }} \tanh (D(x - (1 - 2{D^2})t)).

Ansatz approach to KS equation

The KS equation is given by Yt+αYYx+βYxxx+γ(YYx)xx+σYxYxx=0. {Y_t} + \alpha Y{Y_x} + \beta {Y_{xxx}} + \gamma {(Y{Y_x})_{xx}} + \sigma {Y_x}{Y_{xx}} = 0. The equation under consideration is defined by the constants α, β, γ, and σ, as proposed by Kudryashov and Sinelshchikov to elucidate pressure waves in a mixture of liquid and gas bubbles, accounting for fluid viscosity and heat transfer [36]. In our exploration, we delve into topological, non-topological, and rogue wave solutions for this equation. While Akram et al. [37] previously identified kink, periodic, dark, and bright soliton solutions, our current investigation focuses on dark, bright, and rogue wave solutions. The addition of rogue wave solutions contributes to the richness of the KS equation’s theoretical framework.

Non-topological (bright) soliton solutions

In this section, we have employed the following non-topological wave ansatz to derive solutions for equation (39). The chosen ansatz takes the form Y(x,t)=Ecoshy(κ), Y(x,t) = \frac{E}{{\mathop {\cosh }\nolimits^y (\kappa )}}, in which κ = C(xvt). In the provided equations, the variables E and C are denoted as unconstrained parameters. Upon substituting equation (40) into equation (39) and carrying out some manipulation, we reach vC2β(y2+6y+4)C2β(y+1)(y+2)cosh2(κ)+C2Ey(y+2)+C2γE(y2+6y+4)coshy(κ)C2Ey(y+1)+C2γE(y+1)(y+2)coshy+2(κ)=0. \begin{array}{*{20}{r}}{v - {C^2}\beta ({y^2} + 6y + 4) - \frac{{{C^2}\beta (y + 1)(y + 2)}}{{\mathop {\cosh }\nolimits^2 (\kappa )}} + \frac{{{C^2}Ey(y + 2) + {C^2}\gamma E({y^2} + 6y + 4)}}{{\mathop {\cosh }\nolimits^y (\kappa )}}}\\{ - \frac{{{C^2}Ey(y + 1) + {C^2}\gamma E(y + 1)(y + 2)}}{{\mathop {\cosh }\nolimits^{y + 2} (\kappa )}} = 0.}\end{array} Applying the homogeneous balancing principle, we assume that y = 2 and C2(y + 1)(y + 2) ≠ 0. Upon setting y = 2, we obtain equation (41) as v20C2β12C2β+8C2E+20C2γEcosh2(κ)6C2E+12C2γEcosh4(κ)=0. \begin{array}{*{20}{l}}{v - 20{C^2}\beta - \frac{{12{C^2}\beta + 8{C^2}E + 20{C^2}\gamma E}}{{\mathop {\cosh }\nolimits^2 (\kappa )}} - \frac{{6{C^2}E + 12{C^2}\gamma E}}{{\mathop {\cosh }\nolimits^4 (\kappa )}} = 0.}\end{array} Consequently, we obtain the following system of equations v20C2β=0,12C2β+8C2E+20C2γE=0,6C2E+12C2γE=0. \begin{array}{*{20}{r}}{v - 20{C^2}\beta = 0,}\\{12{C^2}\beta + 8{C^2}E + 20{C^2}\gamma E = 0,}\\{6{C^2}E + 12{C^2}\gamma E = 0.}\end{array} The speed of the soliton, for arbitrary values of C and β, can be calculated as v=20C2β. v = - 20{C^2}\beta . Subsequently, we determine the amplitude of the soliton, denoted as E E=6β. E = - 6\beta . Therefore, the bright soliton solution for the KS equation (39) is given by Y(x,t)=6βcosh2(C(x+20C2βt)) Y(x,t) = - \frac{{6\beta }}{{\mathop {\cosh }\nolimits^2 (C(x + 20{C^2}\beta t))}} \cdot

Rogue wave solutions

In this section, the rogue wave solutions for equation (39) have been generated using the following wave ansatz Y(x,t)=Esinhp(κ), Y(x,t) = \frac{E}{{\mathop {\sinh }\nolimits^p (\kappa )}}, in which ξ = C(xvt). In the provided equations, the variables E and C are denoted as unconstrained parameters. Substituting equation (47) into equation (39) and performing some manipulation, we obtain v+C2βp2+C2β(p+1)(p+2)cosh2(κ)+EαC2Ep2(3γ+σ)+EC2γp2sinhp(κ)+EγC2(p+1)(p+2)C2Ep(p+1)(3γ+σ)sinhp+2(κ)=0. \begin{array}{*{20}{r}}{v + {C^2}\beta {p^2} + \frac{{{C^2}\beta (p + 1)(p + 2)}}{{\mathop {\cosh }\nolimits^2 (\kappa )}} + \frac{{E\alpha - {C^2}E{p^2}(3\gamma + \sigma ) + E{C^2}\gamma {p^2}}}{{\mathop {\sinh }\nolimits^p (\kappa )}}}\\{ + \frac{{E\gamma {C^2}(p + 1)(p + 2) - {C^2}Ep(p + 1)(3\gamma + \sigma )}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}} = 0.}\end{array} Applying the homogeneous balancing principle, we assume that p = 2. Consequently, equation (48) transforms into v+4C2β+12C2β+Eα4C2Eσ8EC2γsinh2(κ)6EC2(γ+σ)sinhp+2(κ)=0. \begin{array}{*{20}{l}}{v + 4{C^2}\beta + \frac{{12{C^2}\beta + E\alpha - 4{C^2}E\sigma - 8E{C^2}\gamma }}{{\mathop {\sinh }\nolimits^2 (\kappa )}} - \frac{{6E{C^2}(\gamma + \sigma )}}{{\mathop {\sinh }\nolimits^{p + 2} (\kappa )}} = 0.}\end{array} By comparing powers of sinh κ, we obtain the following system from equation (49) v+4C2β=0,12C2β+Eα4C2Eσ8EC2γ=0,6EC2(γ+σ)=0. v + 4{C^2}\beta = 0, 12{C^2}\beta + E\alpha - 4{C^2}E\sigma - 8E{C^2}\gamma = 0, 6E{C^2}(\gamma + \sigma ) = 0. The speed of the wave, for arbitrary values of C and α, is obtained as v=4C2β. v = 4{C^2}\beta . Subsequently, we determine the amplitude of the rogue wave E=12C2βα+4C2σ E = - \frac{{12{C^2}\beta }}{{\alpha + 4{C^2}\sigma }} \cdot Therefore, the rogue wave solution for the KS equation (39) is given by Y(x,t)=12C2βα+4C2σcsch2(C(x+4C2βt)). Y(x,t) = - \frac{{12{C^2}\beta }}{{\alpha + 4{C^2}\sigma }}{{\mathop{\rm csch}\nolimits} ^2}(C(x + 4{C^2}\beta t)).

Topological (dark) soliton solution

By employing the topological wave ansatz explained subsequently, this section has successfully established dark soliton solutions for equation (39). We follow the ansatz Y(x,t)=Etanhp(κ), Y(x,t) = E\, \mathop {\tanh }\nolimits^p (\kappa ), in which κ = C(xβt). In the given equations, the parameters E and C are labeled as unconstrained variables. Upon substituting equation (54) into equation (39) and conducting some manipulations, we deduce q = 2 along with the corresponding parameter values v=4C2β. v = 4{C^2}\beta . Subsequently, we determine the amplitude of the soliton as E=16C2σ12γ+164C2γ+32C2σ E = \frac{{16{C^2}\sigma }}{{12\gamma + 164{C^2}\gamma + 32{C^2}\sigma }} \cdot Therefore, the KS equation (39) features a dark soliton solution given by Y(x,t)=16C2σ12γ+164C2γ+32C2σtanh2(C(x+4C2βt)). Y(x,t) = \frac{{16{C^2}\sigma }}{{12\gamma + 164{C^2}\gamma + 32{C^2}\sigma }}\mathop {\tanh }\nolimits^2 (C(x + 4{C^2}\beta t)).

Physical interpretation

The comprehension of solutions to Nonlinear Partial Differential Equations (NPDEs) is essential for grasping the underlying physical phenomena. Solitons, maintaining their shape and constant velocity during propagation, play a crucial role in mathematics and physics. Distinguishing between topological (dark) and non-topological (bright) solitons is relevant to the modulation instability of the carrier wave train. Topological solitons arise when the carrier wave is unstable to long-wave modulations, while non-topological solitons emerge when the carrier wave is modulationally stable. Rogue waves, also known as freak or killer waves, have evolved from marine folklore to recognized phenomena. These waves, surpassing twice the size of surrounding waves, are highly unpredictable and often appear unexpectedly from different directions than the prevailing wind and waves. The exploration of rogue waves contributes to a deeper understanding of extreme events in fluid dynamics. The study features Figures 1 and 4 presenting non-topological soliton solutions, Figures 2 and 5 showcasing rogue wave solutions, and Figures 3 and 6 illustrating topological soliton solutions. These visual depictions provide valuable insights into the physical structures and behaviors of these solitons and rogue waves, enhancing our comprehension of their real-world manifestations.

Fig. 1

Non-topological (bright) soliton nature of surface wave propagation by (12) with C = 1, α = −1 and t = 0.1, 0.2, 0.3.

Fig. 2

Rogue wave solutions of the equation (24) with C = 0.2, α = −1 and t = 0, 1, 2.

Fig. 3

Topological (dark) soliton nature of surface wave propagation by (37) with D = 4, α = 1 and t = 2, 4, 6.

Fig. 4

Non-topological (bright) soliton solutions of the equation (46) with C = 0.05, β = 1 and t = 0.2, 0.4, 0.6.

Fig. 5

Rogue wave solutions of the equation (53) with C = 1, α = 1, σ = 1, β = 1 and t = 0, 1, 2.

Fig. 6

Topological (dark) soliton solutions of the equation (57) with C = 1, σ = 1, γ = 1, β = −0.5 and t = 0.2, 0.4, 0.6.

Discussion and conclusions

The study of transmissions within dynamic systems is currently a subject of emerging fascination. Utilizing the ansatz approach for the NDMBBM equation and KS equation, we have developed novel elucidations of rogue wave solitons, topological, and non-topological solitons, showing the efficacy of this recently successful method. Existing literature, exemplified by Khan et al. [31], has discussed characteristics of trigonometric and hyperbolic type solutions with dark and bright solitons. Moreover, hyperbolic solutions have been reported by [32, 33]. However, our findings stand out as we discovered dark, bright, and rogue wave solutions, setting them apart from earlier research. Akram et al. [37] identified kink, periodic, dark, and bright soliton solutions, but our study uniquely includes dark, bright, and rogue wave solutions, making the discussion of rogue wave solutions intriguing in the context of the NDMBBM and KS equations. We plan to apply the same ansatz to explore nonlinear evolution equations in mathematical physics for future investigations.

Declarations
Conflict of interests

The authors declare that there is no conflict of interest regarding the publication of this paper.

Author’s contributions

A.H.-Methodology, Writing-Review, Software and Editing and Supervision; H.A.-Resources, Writing-Original Draft and Methodology; F.D.Z-Validation, Supervision, Conceptualization and Formal analysis; N.A.-Data curation, Investigation, plotting and Visualization. The paper has been submitted with the knowledge and consent of all authors.

Funding

Not applicable.

Acknowledgement

The second author extends their gratitude to the Abdus Salam School of Mathematical Sciences for their valuable support during the research.

Data availability statement

All data that support the findings of this study are included within the article.

Using of AI tools

The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.

eISSN:
2956-7068
Idioma:
Inglés
Calendario de la edición:
2 veces al año
Temas de la revista:
Computer Sciences, other, Engineering, Introductions and Overviews, Mathematics, General Mathematics, Physics