
The objective of this work is to provide the method of getting the closed-form solitary wave solution of the fractional (3+1)-generalized nonlinear wave equation that characterizes the behavior of liquids with gas bubbles. The same phenomena are evident in science, engineering, and even in the field of physics. This is done by employing the Riccati-Bernoulli sub-ode in a systematic manner as applied to the Bäcklund transformation in the study of this model. New soliton solutions, in the forms of soliton, are derived in the hyperbolic and trigonometric functions. The used software is the computational software Maple, which makes it possible to perform all the necessary calculations and the check of given solutions. The result of such calculations is graphical illustrations of the steady-state characteristics of the system and its dynamics concerning waves and the inter-relationships between the parameters. Moreover, the contour plots and the three-dimensional figures describe the essential features, helping readers understand the physical nature of the model introduced in this work.
Citation: Musawa Yahya Almusawa, Hassan Almusawa. Dark and bright soliton phenomena of the generalized time-space fractional equation with gas bubbles[J]. AIMS Mathematics, 2024, 9(11): 30043-30058. doi: 10.3934/math.20241451
[1] | Naeem Ullah, Muhammad Imran Asjad, Jan Awrejcewicz, Taseer Muhammad, Dumitru Baleanu . On soliton solutions of fractional-order nonlinear model appears in physical sciences. AIMS Mathematics, 2022, 7(5): 7421-7440. doi: 10.3934/math.2022415 |
[2] | Aly R. Seadawy, Bayan Alsaedi . Contraction of variational principle and optical soliton solutions for two models of nonlinear Schrödinger equation with polynomial law nonlinearity. AIMS Mathematics, 2024, 9(3): 6336-6367. doi: 10.3934/math.2024309 |
[3] | Khalid Khan, Amir Ali, Muhammad Irfan, Zareen A. Khan . Solitary wave solutions in time-fractional Korteweg-de Vries equations with power law kernel. AIMS Mathematics, 2023, 8(1): 792-814. doi: 10.3934/math.2023039 |
[4] | Imran Siddique, Khush Bukht Mehdi, Sayed M Eldin, Asim Zafar . Diverse optical solitons solutions of the fractional complex Ginzburg-Landau equation via two altered methods. AIMS Mathematics, 2023, 8(5): 11480-11497. doi: 10.3934/math.2023581 |
[5] | Ghazala Akram, Maasoomah Sadaf, Mirfa Dawood, Muhammad Abbas, Dumitru Baleanu . Solitary wave solutions to Gardner equation using improved tan$ \left(\frac{\Omega(\Upsilon)}{2}\right) $-expansion method. AIMS Mathematics, 2023, 8(2): 4390-4406. doi: 10.3934/math.2023219 |
[6] | Ghazala Akram, Saima Arshed, Maasoomah Sadaf, Hajra Mariyam, Muhammad Nauman Aslam, Riaz Ahmad, Ilyas Khan, Jawaher Alzahrani . Abundant solitary wave solutions of Gardner's equation using three effective integration techniques. AIMS Mathematics, 2023, 8(4): 8171-8184. doi: 10.3934/math.2023413 |
[7] | Naher Mohammed A. Alsafri, Hamad Zogan . Probing the diversity of kink solitons in nonlinear generalised Zakharov-Kuznetsov-Benjamin-Bona-Mahony dynamical model. AIMS Mathematics, 2024, 9(12): 34886-34905. doi: 10.3934/math.20241661 |
[8] | Wafaa B. Rabie, Hamdy M. Ahmed, Taher A. Nofal, Soliman Alkhatib . Wave solutions for the (3+1)-dimensional fractional Boussinesq-KP-type equation using the modified extended direct algebraic method. AIMS Mathematics, 2024, 9(11): 31882-31897. doi: 10.3934/math.20241532 |
[9] | Mohammad Alqudah, Safyan Mukhtar, Haifa A. Alyousef, Sherif M. E. Ismaeel, S. A. El-Tantawy, Fazal Ghani . Probing the diversity of soliton phenomena within conformable Estevez-Mansfield-Clarkson equation in shallow water. AIMS Mathematics, 2024, 9(8): 21212-21238. doi: 10.3934/math.20241030 |
[10] | Islam Samir, Hamdy M. Ahmed, Wafaa Rabie, W. Abbas, Ola Mostafa . Construction optical solitons of generalized nonlinear Schrödinger equation with quintuple power-law nonlinearity using Exp-function, projective Riccati, and new generalized methods. AIMS Mathematics, 2025, 10(2): 3392-3407. doi: 10.3934/math.2025157 |
The objective of this work is to provide the method of getting the closed-form solitary wave solution of the fractional (3+1)-generalized nonlinear wave equation that characterizes the behavior of liquids with gas bubbles. The same phenomena are evident in science, engineering, and even in the field of physics. This is done by employing the Riccati-Bernoulli sub-ode in a systematic manner as applied to the Bäcklund transformation in the study of this model. New soliton solutions, in the forms of soliton, are derived in the hyperbolic and trigonometric functions. The used software is the computational software Maple, which makes it possible to perform all the necessary calculations and the check of given solutions. The result of such calculations is graphical illustrations of the steady-state characteristics of the system and its dynamics concerning waves and the inter-relationships between the parameters. Moreover, the contour plots and the three-dimensional figures describe the essential features, helping readers understand the physical nature of the model introduced in this work.
The fractional differential equations have been the focus of many contributions, taking into account their application in various areas such as finance, engineering, biology and physics, control theory, system identification, and signal processing [1,2,3,4,5]. Their applications include; social sciences, dietary supplements, climate, and economics, among others [6,7,8,9,10,11]. Hence, the exact solutions to fractional differential equations are important. To solve these equations, numerous analytical and numerical approaches have been reported, in the literature such as the modified F-expansion systematic [12], the extended tanh-coth method [13], the mapping method [14], the (G'/G)-expansion technique [15], etc. Rayleigh [16] did the first work on bubble dynamics. Since many subjects and industries deal with bubbly liquids, including medical science and engineering, bubbly liquids have received significant attention. This collection of recent research articles covers a broad range of topics in mathematical physics, nonlinear dynamics, and applied mathematics. They delve into Gaussian traveling wave solutions for a Schrödinger equation with logarithmic nonlinearity [17], while their subsequent work explores exact solutions for a nonlinear fourth-order time-fractional partial differential equation [18]. Focus on soliton solutions to the (2+1) dimensional Chaffee-Infante equation [19]. Investigate synchronization patterns in memristive neuron maps [20], alongside studies of ADHD brain networks and chimera states in FitzHugh-Nagumo oscillators [21,22]. The present digital integrators are based on trigonometric quadrature rules, contributing to advancements in industrial electronics [23]. It has been found out that the propagation of linear acoustic waves in isothermal bubbly liquids, assuming that the radius of the bubbles is uniform, is governed by a fourth-order linear partial differential equation. A generalized (3+1)-dimensional nonlinear wave equation is one such model used to characterize a liquid containing gas bubbles:
∂∂x(∂f∂t+(p1)∂f∂t+(p2)∂3f∂x3+(p3)f∂f∂x)+(p4)∂2f∂z2+(p5)∂2f∂y2=0. | (1.1) |
In the context of the equation, f(x,y,z,t) stands for the amplitude of a wave, where p1 is the bubble liquid viscosity and p2 is the bubble liquid dispersion, p3 is the bubble liquid nonlinearity, p4 is the z transverse perturbation, and p5 is the transverse perturbation. Equation (1.1) can be used for the description of different nonlinear processes occurring in liquids with gas bubbles. Therefore, this equation has received a lot of attention in research. For example, within the frame of linear superposition, soliton solutions of N-soliton waves have been obtained [24]. Furthermore, lump-stripe solitons and rogue wave-stripe solutions have been built [25], and lump, stripe periodic, and multi solitons also have been found [26]. Soliton solutions have been obtained by using the technique that is now known as the Hirota bilinear form and has resulted in a new form of a generalized exponential rational function of Hirota [27,28]. Additionally, the use of the modified Kudryashov method and Nuccis reduction has been made to investigate the behavior of solitary waves [29]. However, in the past studies, the consideration of fractional derivatives on Eq (1.1) has not been thoroughly investigated, which hinders the ability to simulate the detailed kinetics of wave in systems having memory as well as non-local interactions. This has limited the modeling of wave interactions in near-neighbourhood systems with memory and non-locality. To this end, in this study, we will extend the above analysis to the generalized fractional nonlinear wave equation (GFNWE). By means of this, we obtain a more reliable and extensive model for nonlinear wave interactions that is significant for the description of various processes in fluid dynamics, plasma physics, and other complex systems. In this study, we extend this work by examining the generalized fractional nonlinear wave equation (GFNWE) as follows:
Dαx(Dαtf+(p1)Dαxf+(p2)D3αxf+(p3)fDαxf)+(p4)D2αzf+(p5)D2αyf=0. | (1.2) |
Further, the operator integrating α-derivatives of powers agrees exactly to the idea of conformable fractional derivatives [30].
DαϕW(ϕ)=limi→0W(ϕ+i(ϕ)1−α)−W(ϕ)i,0<α≤1. | (1.3) |
{Dαϕϕk=kϕk−α.Dαϕ(k1η(ϕ)±k2t(ϕ))=k1Dαϕ(η(ϕ))±m2Dαϕ(t(ϕ)).Dαϕ[f∘g]=ϕ1−αg(ϕ)Dαϕf(g(ϕ)). | (1.4) |
The Riccati-Bernoulli sub-ode method is a new effective analytical tool for obtaining the solitary wave solutions of PDEs and FPDEs. This method works on the basis of altering the original FPDE into a first-order ODE, that is nonlinear. Subsequent to this transformation, there is the use of a series-form solution, which results in the derivation of non-linear algebraic equations. Solving the algebraic equations as mentioned above helps in the determination of the solitary wave solutions of the concerned FPDE only [31,32,33,34,35,36]. This approach is specifically efficient so far as it concerns other analytical techniques utilized for deriving soliton solutions. The Riccati-Bernoulli sub-ode method allows the identification of other more wide-spread and multifaceted subfamilies of soliton solutions, thus, proving its vitality in describing various and complicated aspects of solitary waves in various systems [37,38,39].
The second part of the paper describes the overall approach used; a detailed explanation is given in Section 3 when the solution for the generalized fractional nonlinear wave equation is discussed. The analysis of the results is provided in Section 4 along with graphical representations of the findings. Lastly, Section 5 provides the conclusion of the study where the major observations and recommendations will be made.
Since it is crucial to provide the reader with a clear understanding of the methods used in this research, a slightly deeper explanation of the process will be given below. We begin by considering a general class of nonlinear partial differential equations (PDEs) in the following form:
P1(R1,Dαt(R1),Dαq1(R1),Dαq2(R1),R1Dαq1(R1),…)=0,0<α≤1, | (2.1) |
where R1=R(t,q1,q2,q3,…,qk) is a function of (t,q1,q2,q3,…,qk) and its partial derivatives. This transformation changes Eq (2.1) into a nonlinear ordinary differential equation (ODE) of the following form:
Q1(F,F′(ϕ),F″(ϕ),FF′(ϕ),…)=0. | (2.2) |
Let us suppose that Eq (2.2) has the following solution:
G(ϕ)=n∑j=−nkjg(ϕ)j, | (2.3) |
where kj are constants and g(ϕ) is obtained from the Bäcklund transformation, g(ϕ)=−ΩE2+E1Z(ϕ)E1+E2Z(ϕ). Where, (Ω), (E1), and (E2) are constants such that E2≠0 and Z(ϕ) are solutions of the following ODE.
dZdϕ=Ω+Z(ϕ)2. | (2.4) |
The Ricatti equation (2.4) possesses the following general solutions [40].
Z(ϕ)={−√−Ωtanh(√−Ωϕ),as Ω<0,−√−Ωcoth(√−Ωϕ),as Ω<0,Z(ϕ)=−1ψ,as Ω=0,Z(ϕ)={√Ωtan(√Ωϕ),as Ω>0,−√Ωcot(√Ωϕ),as Ω>0. | (2.5) |
To eliminate the homogeneous term, which is the ratio of the largest nonlinear term and the highest order derivative on the right-hand side of Eq (2.2), a positive integer n is derived as given in Eq (2.3). Therefore, while presenting the balance number of a processor, it is pertinent to mention that it can be arrived at using the following method [41].
D[dmFdψm]=n+m,D[FJdmFdψm]w=nJ+w(m+n). | (2.6) |
Subsequently, the function obtained from Eq (2.3) is plugged into Eq (2.2) or in the expression that we obtain after integrating Eq (2.2). Subsequently, all terms with g(ϕ) are grouped, and the coefficients of the polynomial are set to zero. This process brings out a system of algebraic equations with the use of (ki) and other variables applicable in the sequence. These algebraic equations are further solved using the Maple a computational software. Last of all, the solutions of Eq (1.1) are derived as single wave solutions and are displayed.
In this section, various solutions are discussed for the fractional (3+1)-dimensional generalized wave model of liquids containing gas bubbles, which has been defined in Eq (1.1). According to the Riccati-Bernoulli sub-ode technique and Bäcklund transformation process, solitary wave solutions are attained. The wave transformation employed in this context is as follows:
f(x,y,z,t)=F(ψ),whereψ=μ1xαα+μ2yαα+μ3zαα−ωtαα. | (3.1) |
Where, μ1,μ2,μ3 and ω are unknown constants. This transformation changes Eq (1.2) into a nonlinear ordinary differential equation (ODE) of the following form:
μ41p3(d2Fdψ2)+μ21p1(F(ψ)2)+F(ψ)(μ21p3−μ1ω+μ22p4+μ23p5)=0. | (3.2) |
This makes it possible for us to reduce and solve the system balancing equations and other multiple wave structures that are complex in a new approach that is adopted. Used in a systematic way, such integration allows extracting particular traits of the primary constituents of waves, next to which conventional concepts are employed. By this distinct integration, the characteristics of interactions and behaviors of waves can be observed. With a view to furthering the solution, Eqs (2.3) and (2.4) are integrated into the governing Eq (3.2). Thus, we gather coefficients of the function Z(ϕ) in detail and obtain an analogous system of equations. This system provides a basis for the analysis of the objects and future quantitative or qualitative analysis and represents the key characteristics of the wave structures.
Z1(ϕ)=μ12p1k−22E28+6μ14p3k−2E28Ω2=0,Z2(ϕ)=−2μ12p1k−2E28k−1Ω−2μ14p3k−1E28Ω3=0,Z3(ϕ)=k−2E28μ22p4Ω2+μ12p1k−12E28Ω2+k−2E28μ32p5Ω2+k−2E28μ12p3Ω2−k−2E28μ1ωΩ2+8μ14p3k−2E28Ω3+2μ12p1k−2E28k0Ω2=0,Z4(ϕ)=−2μ12p1k−1E28k0Ω3+k−1E28μ1ωΩ3−k−1E28μ22p4Ω3−2μ12p1k−2E28k1Ω3−2μ14p3k−1E28Ω4−k−1E28μ32p5Ω3−k−1E28μ12p3Ω3=0,Z5(ϕ)=μ12p1k02Ω4E28+k0μ32p5Ω4E28+k0μ12p3Ω4E28+k0μ22p4Ω4E28+2μ14p3k2E28Ω6+2μ14p3k−2E28Ω4−k0μ1ωΩ4E28+2μ12p1k−2E28k2Ω4+2μ12p1k−1E28k1Ω4=0,Z6(ϕ)=−2μ12p1k−1E28k2Ω5−2μ12p1k0k1Ω5E28+k1Ω5E28μ1ω−2μ14p3k1E28Ω6−k1Ω5E28μ32p5−k1Ω5E28μ22p4−k1Ω5E28μ12p3=0,Z7(ϕ)=μ12p1k12Ω6E28+k2Ω6E28μ12p3+k2Ω6E28μ22p4+k2Ω6E28μ32p5+2μ12p1k0k2Ω6E28+8μ14p3k2E28Ω7−k2Ω6E28μ1ω=0,Z8(ϕ)=−2μ12p1k1Ω7E28k2−2μ14p3k1E28Ω7=0,Z9(ϕ)=6μ14p3k2E28Ω8+μ12p1k22Ω8E28=0. | (3.3) |
This gives us the algebraic equations by setting Z(ϕ)=0. The solutions of this system of algebraic equations obtained from Maple are:
Set 1.
k0=k0,k1=0,k−1=0,k−2=−3/8p1k02μ12p3,k2=−6μ12p3p1,Ω=1/4p1k0μ12p3,μ1=μ1,μ2=μ2,μ3=√−−μ1ω+p4μ22+μ12p3+4μ12p1k0p5. | (3.4) |
Set 2.
k0=k0,k1=0,k−1=0,k−2=0,k2=−6μ12p3p1,Ω=−1/2p1k0μ12p3,μ1=μ1,μ2=μ2,μ3=√−−2μ12p1k0−μ1ω+p4μ22+μ12p3p5. | (3.5) |
Solution Set 1
Set 1. For (Ω<0), μ3=√−−μ1ω+p4μ22+μ12p3+4μ12p1k0p5 and ψ=μ1xαα+μ2yαα+μ3zαα−ωtαα, the following set of solutions for Eq (2.6) are obtained.
f1=−3/8p1k02(E1−E2√−Ωtanh(√−Ωψ))2μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2+k0−6μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2p1(E1−E2√−Ωtanh(√−Ωψ))2, | (3.6) |
or
f2=−3/8p1k02(E1−E2√−Ωcoth(√−Ωψ))2μ12p3(−ΩE2−E1√−Ωcoth(√−Ωψ))2+k0−6μ12p3(−ΩE2−E1√−Ωcoth(√−Ωψ))2p1(E1−E2√−Ωcoth(√−Ωψ))2. | (3.7) |
Solution Set 2
Set 1. For (Ω>0), the following set of solutions for Eq (2.6) are obtained.
f3=−3/8p1k02(E1+E2√Ωtan(√Ωψ))2μ12p3(−ΩE2+E1√Ωtan(√Ωψ))2+k0−6μ12p3(−ΩE2+E1√Ωtan(√Ωψ))2p1(E1+E2√Ωtan(√Ωψ))2, | (3.8) |
or
f4=−3/8p1k02(E1−E2√Ωcot(√Ωψ))2μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2+k0−6μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2p1(E1−E2√Ωcot(√Ωψ))2. | (3.9) |
Solution Set 3
Set 2. For (Ω<0) and μ3=√−−2μ12p1k0−μ1ω+p4μ22+μ12p3p5., the following set of solutions for Eq (2.6) are obtained.
f5=k0−6μ12p3(−ΩE2−E1√−Ωtanh(√−Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))2p1−1(E1−E2√−Ωtanh(√−Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))−2, | (3.10) |
or
f6=k0−6μ12p3(−ΩE2−E1√−Ωcoth(√−Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))2p1−1(E1−E2√−Ωcoth(√−Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))−2. | (3.11) |
Solution Set 4
Set 2. For (Ω>0), the following set of solutions for Eq (2.6) are obtained.
f7=k0−6μ12p3(−ΩE2+E1√Ωtan(√Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))2p1−1(E1+E2√Ωtan(√Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))−2, | (3.12) |
or
f8=k0−6μ12p3(−ΩE2−E1√Ωcot(√Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))2p1−1(E1−E2√Ωcot(√Ω(μ1xαα+μ2yαα+μ3zαα−ωtαα)))−2. | (3.13) |
The use of the Riccati-Bernoulli sub-ODE method with the Bäcklund transformation gave us exact solutions, which provided wider applicability in comparison with other approaches, such as the sine-cosine or the homotopy perturbation. It also allows for a powerful solution in handling of nonlinearity and fractional derivatives; this reduces the number of steps to be followed since FPDE is transformed into an algebraic statement. This results more stable and accurate waveforms compared to recent approaches applicable to fluid dynamics and plasma physics, as seen from the comparison. Using the figures of 3D surface plots and contour plots, it is easy to understand the physical changes of solitary wave solutions for the generalized fractional nonlinear wave equation with p1=ω=0.001, μ1=0.05, μ2=0.05, μ3=0.0315, p2=2, and p3=3. Both f1 and f8 solitary wave solutions depict dark solitary waves, which consist of regions of low amplitude enclosed by regions of high intensity, thereby illustrating how dark solitary waves create depressions in the media through which they travel. On the other hand, the bright solitary wave solutions f4 and f6 display the maximum peak in the amplitudes, and they represent enhanced regional intensity of the wave. Figure 1, 3D and contour representation of the solution, f1(x,y,z,t). Figure 2, 3D and contour representation of the solution, f4(x,y,z,t). Figure 3, 3D and contour representation of the solution, f6(x,y,z,t). Figure 4, 3D and contour representation of the solution, f8(x,y,z,t). Table 1, Comparison of the current approach with the alternative approach, specifically modified Kudryashov method [29]. These representations added emphasis on the spatial relationships as well as the changes in amplitudes which, showing how the wave behaves with its surroundings and how it advances. Taking into account these plots, one obtains more information about the actual behavior of wave solutions and some characteristics such as their stability and some features of their interactions with other solutions underlying physical processes modeled by the equations introduced above. The solutions are trigonometric and hyperbolic, and thus come with a rich structure, indicative of the wide range of physical processes that can be modeled. These types of solutions are of special importance for the description of liquid phases with gaseous bubbles, which are often met in different scientific and engineering practices. The knowledge of such waves in the field of fluid dynamics is significant to describe the multiphase flow as the key to its behavior consists of the interaction between the liquid and gaseous phases. For instance, accurate descriptions of bubbly flows, in the oil and gas industry are critical in enhancing the extraction processes and the stability of pipelines.
Case Ⅰ: Ω<0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√−Ωtanh(√−Ωψ))2μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2+k0−6μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2p1(E1−E2√−Ωtanh(√−Ωψ))2. |
Case Ⅰ: g2(g3l21+g4l22+g5l23−l1v)<0 (modified Kudryashov) method |
f(x,y,z,t)=3(g3l21+g4l22+g5l23−l1v)g1l21×csch2[√−g3l21−g4l22−g5l23+l1v4g2l41(l1x+l2y+l3z−vt)]. |
Case Ⅱ: Ω>0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√Ωcot(√Ωψ))2μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2+k0−6μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2p1(E1−E2√Ωcot(√Ωψ))2. |
Case Ⅱ: g2(g3l21+g4l22+g5l23−l1v)>0 (modified Kudryashov) method |
f(x,y,z,t)=g3l21+g4l22+g5l23−l1vg1l21×(−3csch2[√g3l21+g4l22+g5l23−l1v4g2l41(l1x+l2y+l3z−vt)]−2). |
The nonlinear wave equation derived in this work as a fractional derivative can be a useful model for liquids containing gaseous bubbles that are observed for their complicated wave interaction with a given environment. This makes the fractional derivatives important for the description of the memory and nonlocality in the fluid and necessary for the description of bubble, containing liquids. The solitary wave solutions obtained from the model reflect localized waveform structure, which can travel without loss of form or amplitude, hence resembling pressure waves in such media. These solutions offer information on wave behaviour, dispersion, and stability, all of which are critical to the behaviour of gas bubbles in diverse liquid systems, be it industrial or observed in natural occurrences. Moreover, these wave solutions are important when discussing the problem of wave carrying in the bubbly medium, a field important in the sphere of physics. Such models are employed to estimate the processes affecting sound waves in the media containing bubble gas; it is crucial in various fields, starting with underwater acoustics and ending with ultrasonography. In underwater acoustics, these models are used in sonar systems formulation so that the sound can be accurately used to find objects in bubbly water and in medical ultrasound to improve the techniques used in diagnosis.
In the next steps, potential future studies may be aimed at the further development of the described analytical approach with the inclusion of more elaborate models of a higher dimensionality and different forms of nonlinearity. In the same way, the enhancement of these theoretical models by the experimental confirmation would be the superior achievement, especially concerning the real-life applications in fluid dynamics of industrial processes, advanced ultrasounds, and diagnostics in medicine. The investigation of the above-mentioned relations between fractional order parameters and physical characteristics might help to expand the knowledge of wave processes in the media and to create more detailed models in the future. These advancements could make a breath of fresh air in the areas of engineering and technology that require the control and change of wave behavior.
In this study, considering the concept of the Bäcklund transformation integrated with the Riccati-Bernoulli sub-ode method, the solitary wave solutions for the generalized fractional nonlinear wave equation have been derived and analyzed. The mentioned approach allowed us to obtain a wide set of analytical solutions in the form of hyperboles, trigonometric functions, and rational expressions. These solutions were accurately and completely depicted in terms of space by both the 3D surface plot and the contour plot regarding the wave behavior. Another innovation has been where wave systems are nonlinear, and the use of 3D plus contour visualization has been able to provide a clearer picture of the various behaviors of the wave systems, including applications in fluid dynamics, chemical kinetics, and biological morphogenesis. The fact that the outlined methodology allows for the depiction of the diffusion reaction processes with sufficient accuracy to capture the intricate details of the interaction between the processes underlies one of the primary perceived values derived in this work the solutions obtained are practical in nature. This paper shows that both the Bäcklund transformation and the sub-ode method based on the Riccati-Bernoulli equation can be used as a powerful instruments when studying the generalized fractional nonlinear wave equation. What has been gained by reducing the equation to a form that admits an ordinary differential equation solution is an asset that is both mathematically sound to the extent that the analytical method is valid or reliable and qualitatively useful in charting the nature of solitary waves.
Musawa Yahya Almusawa: Conceptualization, Formal analysis, Investigation, Project administration, Validation, Visualization, Funding; Hassan Almusawa: Data curation, Resources, Validation, Software, Project administration, Writing-review & editing. All authors have read and agreed to the published version of the manuscript.
The authors gratefully acknowledge the funding of the Deanship of Graduate Studies and Scientific Research, Jazan University, Saudi Arabia, through project number: RG24-L02.
The authors declare that they have no conflicts of interest.
[1] | K. S. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, New York: Wiley, 1993. |
[2] | V. Daftardar-Gejji, Fractional calculus and fractional differential equations, Singapore: Springer, 2019. https://doi.org/10.1007/978-981-13-9227-6 |
[3] |
S. Mukhtar, M. Sohaib, I. Ahmad, A numerical approach to solve volume-based batch crystallization model with fines dissolution unit, Processes, 7 (2019), 453. https://doi.org/10.3390/pr7070453 doi: 10.3390/pr7070453
![]() |
[4] |
M. Bouloudene, M. A. Alqudah, F. Jarad, Y. Adjabi, T. Abdeljawad, Nonlinear singular P-Laplacian boundary value problems in the frame of conformable derivative, Discret. Contin. Dyn. Syst. S, 14 (2021), 3497–3528. https://doi.org/10.3934/dcdss.2020442 doi: 10.3934/dcdss.2020442
![]() |
[5] |
T. Abdeljawad, On conformable fractional calculus, J. Comput. Appl. Math., 279 (2015), 57–66. https://doi.org/10.1016/j.cam.2014.10.016 doi: 10.1016/j.cam.2014.10.016
![]() |
[6] |
H. Yasmin, A. S. Alshehry, A. H. Ganie, A. M. Mahnashi, R. Shah, Perturbed Gerdjikov-Ivanov equation: Soliton solutions via Bäcklund transformation, Optik, 298 (2024), 171576. https://doi.org/10.1016/j.ijleo.2023.171576 doi: 10.1016/j.ijleo.2023.171576
![]() |
[7] |
A. Zafar, M. Raheel, A. M. Mahnashi, A. Bekir, M. R. Ali, A. S. Hendy, Exploring the new soliton solutions to the nonlinear M-fractional evolution equations in shallow water by three analytical techniques, Results Phys., 54 (2023), 107092. https://doi.org/10.1016/j.rinp.2023.107092 doi: 10.1016/j.rinp.2023.107092
![]() |
[8] |
W. Hamali, J. Manafian, M. Lakestani, A. M. Mahnashi, A. Bekir, Optical solitons of M-fractional nonlinear Schrödinger's complex hyperbolic model by generalized Kudryashov method, Opt. Quantum Electron., 56 (2024), 7. https://doi.org/10.1007/s11082-023-05602-1 doi: 10.1007/s11082-023-05602-1
![]() |
[9] |
I. Ahmad, H. Seno, An epidemic dynamics model with limited isolation capacity, Theory Biosci., 142 (2023), 259–273. https://doi.org/10.1007/s12064-023-00399-9 doi: 10.1007/s12064-023-00399-9
![]() |
[10] |
N. Iqbal, M. B. Riaz, M. Alesemi, T. S. Hassan, A. M. Mahnashi, A. Shafee, Reliable analysis for obtaining exact soliton solutions of (2+1)-dimensional Chaffee-Infante equation, AIMS Mathematics, 9 (2024), 16666–16686. https://doi.org/10.3934/math.2024808 doi: 10.3934/math.2024808
![]() |
[11] |
A. S. Alshehry, S. Mukhtar, A. M. Mahnashi, Analytical methods in fractional biological population modeling: Unveiling solitary wave solutions, AIMS Mathematics, 9 (2024), 15966–15987. https://doi.org/10.3934/math.2024773 doi: 10.3934/math.2024773
![]() |
[12] |
W. W. Mohammed, F. M. Al-Askar, C. Cesarano, M. El-Morshedy, Solitary wave solutions of the fractional-stochastic quantum Zakharov-Kuznetsov equation arises in quantum magneto plasma, Mathematics, 11 (2023), 488. https://doi.org/10.3390/math11020488 doi: 10.3390/math11020488
![]() |
[13] |
W. W. Mohammed, C. Cesarano, F. M. Al-Askar, Solutions to the (4+1)-dimensional time-fractional Fokas equation with M-truncated derivative, Mathematics, 11 (2023), 194. https://doi.org/10.3390/math11010194 doi: 10.3390/math11010194
![]() |
[14] |
W. W. Mohammed, F. M. Al-Askar, C. Cesarano, The analytical solutions of the stochastic mKdV equation via the mapping method, Mathematics, 10 (2022), 4212. https://doi.org/10.3390/math10224212 doi: 10.3390/math10224212
![]() |
[15] |
F. M. Al-Askar, C. Cesarano, W. W. Mohammed, The analytical solutions of stochastic-fractional Drinfeld-Sokolov-Wilson equations via (G'/G)-expansion method, Symmetry, 14 (2022), 2105. https://doi.org/10.3390/sym14102105 doi: 10.3390/sym14102105
![]() |
[16] | Y. Xie, I. Ahmad, T. I. Ikpe, E. F. Sofia, H., Seno, What influence could the acceptance of visitors cause on the epidemic dynamics of a reinfectious disease?: A mathematical model, Acta Biotheor., 72 (2024), 3. |
[17] |
Y. Kai, Z. X. Yin, On the Gaussian traveling wave solution to a special kind of Schrödinger equation with logarithmic nonlinearity, Mod. Phys. Lett. B, 36 (2022), 2150543. https://doi.org/10.1142/S0217984921505436 doi: 10.1142/S0217984921505436
![]() |
[18] |
Y. Kai, S. Q. Chen, K. Zhang, Z. X. Yin, Exact solutions and dynamic properties of a nonlinear fourth-order time-fractional partial differential equation, Wave Random Complex, 2022. https://doi.org/10.1080/17455030.2022.2044541 doi: 10.1080/17455030.2022.2044541
![]() |
[19] |
C. Y. Zhu, M. Al-Dossari, S. Rezapour, B. Gunay, On the exact soliton solutions and different wave structures to the (2+1) dimensional Chaffee-Infante equation, Results Phys., 57 (2024), 107431. https://doi.org/10.1016/j.rinp.2024.107431 doi: 10.1016/j.rinp.2024.107431
![]() |
[20] |
Z. Wang, F. Parastesh, H. Natiq, J. Li, X. Xi, M. Mehrabbeik, Synchronization patterns in a network of diffusively delay-coupled memristive Chialvo neuron map, Phys. Lett. A, 514-515 (2024), 129607. https://doi.org/10.1016/j.physleta.2024.129607 doi: 10.1016/j.physleta.2024.129607
![]() |
[21] |
X. Xi, J. Li, Z. Wang, H. Tian, R. Yang, The effect of high-order interactions on the functional brain networks of boys with ADHD, Eur. Phys. J. Spec. Top., 233 (2024), 817–829. https://doi.org/10.1140/epjs/s11734-024-01161-y doi: 10.1140/epjs/s11734-024-01161-y
![]() |
[22] |
Z. Wang, M. Chen, X. Xi, H. Tian, R. Yang, Multi-chimera states in a higher order network of FitzHugh-Nagumo oscillators, Eur. Phys. J. Spec. Top., 233 (2024), 779–786. https://doi.org/10.1140/epjs/s11734-024-01143-0 doi: 10.1140/epjs/s11734-024-01143-0
![]() |
[23] |
T. A. A. Ali, Z. Xiao, H. Jiang, B. Li, A class of digital integrators based on trigonometric quadrature rules, IEEE Trans. Ind. Electron., 71 (2024), 6128–6138. https://doi.org/10.1109/TIE.2023.3290247 doi: 10.1109/TIE.2023.3290247
![]() |
[24] |
G. Shen, J. Manafian, D. T. N. Huy, K. S. Nisar, M. Abotaleb, N. D. Trung, Abundant soliton wave solutions and the linear superposition principle for generalized (3+1)-D nonlinear wave equation in liquid with gas bubbles by bilinear analysis, Results Phys., 32 (2022), 105066. https://doi.org/10.1016/j.rinp.2021.105066 doi: 10.1016/j.rinp.2021.105066
![]() |
[25] |
M. Wang, B. Tian, Y. Sun, Z. Zhang, Lump, mixed lump-stripe and rogue wave-stripe solutions of a (3+1)-dimensional nonlinear wave equation for a liquid with gas bubbles, Comput. Math. Appl., 79 (2020), 576–587. https://doi.org/10.1016/j.camwa.2019.07.006 doi: 10.1016/j.camwa.2019.07.006
![]() |
[26] |
Y. R. Guo, A. H. Chen, Hybrid exact solutions of the (3+1)-dimensional variable-coefficient nonlinear wave equation in liquid with gas bubbles, Results Phys., 23 (2021), 103926. https://doi.org/10.1016/j.rinp.2021.103926 doi: 10.1016/j.rinp.2021.103926
![]() |
[27] |
X. J. Zhou, O. A. Ilhan, F. Y. Zhou, S. Sutarto, J. Manafian, M. Abotaleb, Lump and interaction solutions to the (3+1)-dimensional variable-coefficient nonlinear wave equation with multidimensional binary Bell polynomials, J. Funct. Space, 2021 (2021), 4550582. https://doi.org/10.1155/2021/4550582 doi: 10.1155/2021/4550582
![]() |
[28] |
Y. H. Zhao, T. Mathanaranjan, H. Rezazadeh, L. Akinyemi, M. Inc, New solitary wave solutions and stability analysis for the generalized (3+1)-dimensional nonlinear wave equation in liquid with gas bubbles, Results Phys., 43 (2022), 106083. https://doi.org/10.1016/j.rinp.2022.106083 doi: 10.1016/j.rinp.2022.106083
![]() |
[29] | A. Akbulut, A. H. Arnous, M. S. Hashemi, M. Mirzazadeh, Solitary waves for the generalized nonlinear wave equation in (3+1) dimensions with gas bubbles using Nucci's reduction, enhanced and modified Kudryashov algorithms, J. Ocean. Eng. Sci., 2022, In press. https://doi.org/10.1016/j.joes.2022.07.002 |
[30] | E. M. E. Zayed, K. A. E. Alurrfi, The modified extended tanh-function method and its applications to the generalized KdV-mKdV equation with any-order nonlinear terms, Int. J. Environ. Eng. Sci. Tech. Res., 1 (2013), 165–170. |
[31] |
S. Alshammari, M. M. Al-Sawalha, R. Shah, Approximate analytical methods for a fractional-order nonlinear system of Jaulent-Miodek equation with energy-dependent Schrodinger potential, Fractal Fract., 7 (2023), 140. https://doi.org/10.3390/fractalfract7020140 doi: 10.3390/fractalfract7020140
![]() |
[32] |
A. A. Alderremy, R. Shah, N. Iqbal, S. Aly, K. Nonlaopon, Fractional series solution construction for nonlinear fractional reaction-diffusion Brusselator model utilizing Laplace residual power series, Symmetry, 14 (2022), 1944. https://doi.org/10.3390/sym14091944 doi: 10.3390/sym14091944
![]() |
[33] |
M. M. Al-Sawalha, R. Shah, A. Khan, O. Y. Ababneh, T. Botmart, Fractional view analysis of Kersten-Krasil-shchik coupled KdV-mKdV systems with non-singular kernel derivatives, AIMS Mathematics, 7 (2022), 18334–18359. https://doi.org/10.3934/math.20221010 doi: 10.3934/math.20221010
![]() |
[34] |
E. M. Elsayed, R. Shah, K. Nonlaopon, The analysis of the fractional-order Navier-Stokes equations by a novel approach, J. Funct. Space, 2022 (2022), 8979447. https://doi.org/10.1155/2022/8979447 doi: 10.1155/2022/8979447
![]() |
[35] |
M. Alqhtani, K. M. Saad, R. Shah, W. Weera, W. M. Hamanah, Analysis of the fractional-order local Poisson equation in fractal porous media, Symmetry, 14 (2022), 1323. https://doi.org/10.3390/sym14071323 doi: 10.3390/sym14071323
![]() |
[36] |
M. Naeem, H. Rezazadeh, A. A. Khammash, R. Shah, S. Zaland, Analysis of the fuzzy fractional-order solitary wave solutions for the KdV equation in the sense of Caputo-Fabrizio derivative, J. Math., 2022 (2022), 3688916. https://doi.org/10.1155/2022/3688916 doi: 10.1155/2022/3688916
![]() |
[37] |
M. A. E. Abdelrahman, M. A. Sohaly, Solitary waves for the modified Korteweg-de Vries equation in deterministic case and random case, J. Phys. Math., 8 (2017), 1000214. https://doi.org/10.4172/2090-0902.1000214 doi: 10.4172/2090-0902.1000214
![]() |
[38] |
M. A. E. Abdelrahman, M. A. Sohaly, Solitary waves for the nonlinear Schrödinger problem with the probability distribution function in the stochastic input case, Eur. Phys. J. Plus, 132 (2017), 339. https://doi.org/10.1140/epjp/i2017-11607-5 doi: 10.1140/epjp/i2017-11607-5
![]() |
[39] |
X. F. Yang, Z. C. Deng, Y. Wei, A Riccati-Bernoulli sub-ODE method for nonlinear partial differential equations and its application, Adv. Differ. Equ., 2015 (2015), 117. https://doi.org/10.1186/s13662-015-0452-4 doi: 10.1186/s13662-015-0452-4
![]() |
[40] |
Y. Jawarneh, H. Yasmin, A. H. Ganie, M. M. Al-Sawalha, A. Ali, Unification of adomian decomposition method and zz transformation for exploring the dynamics of fractional kersten-krasil'shchik coupled kdv-mkdv systems, AIMS Mathematics, 9 (2024), 371–390. https://doi.org/10.3934/math.2024021 doi: 10.3934/math.2024021
![]() |
[41] |
S. Bibi, S. T. Mohyud-Din, R. Ullah, N. Ahmed, U. Khan, Exact solutions for STO and (3+1)-dimensional KdV-ZK equations using G′/G-expansion method, Results Phys., 7 (2017), 4434–4439. https://doi.org/10.1016/j.rinp.2017.11.009 doi: 10.1016/j.rinp.2017.11.009
![]() |
Case Ⅰ: Ω<0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√−Ωtanh(√−Ωψ))2μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2+k0−6μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2p1(E1−E2√−Ωtanh(√−Ωψ))2. |
Case Ⅰ: g2(g3l21+g4l22+g5l23−l1v)<0 (modified Kudryashov) method |
f(x,y,z,t)=3(g3l21+g4l22+g5l23−l1v)g1l21×csch2[√−g3l21−g4l22−g5l23+l1v4g2l41(l1x+l2y+l3z−vt)]. |
Case Ⅱ: Ω>0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√Ωcot(√Ωψ))2μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2+k0−6μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2p1(E1−E2√Ωcot(√Ωψ))2. |
Case Ⅱ: g2(g3l21+g4l22+g5l23−l1v)>0 (modified Kudryashov) method |
f(x,y,z,t)=g3l21+g4l22+g5l23−l1vg1l21×(−3csch2[√g3l21+g4l22+g5l23−l1v4g2l41(l1x+l2y+l3z−vt)]−2). |
Case Ⅰ: Ω<0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√−Ωtanh(√−Ωψ))2μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2+k0−6μ12p3(−ΩE2−E1√−Ωtanh(√−Ωψ))2p1(E1−E2√−Ωtanh(√−Ωψ))2. |
Case Ⅰ: g2(g3l21+g4l22+g5l23−l1v)<0 (modified Kudryashov) method |
f(x,y,z,t)=3(g3l21+g4l22+g5l23−l1v)g1l21×csch2[√−g3l21−g4l22−g5l23+l1v4g2l41(l1x+l2y+l3z−vt)]. |
Case Ⅱ: Ω>0 Present method |
f(x,y,z,t)=−3/8p1k02(E1−E2√Ωcot(√Ωψ))2μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2+k0−6μ12p3(−ΩE2−E1√Ωcot(√Ωψ))2p1(E1−E2√Ωcot(√Ωψ))2. |
Case Ⅱ: g2(g3l21+g4l22+g5l23−l1v)>0 (modified Kudryashov) method |
f(x,y,z,t)=g3l21+g4l22+g5l23−l1vg1l21×(−3csch2[√g3l21+g4l22+g5l23−l1v4g2l41(l1x+l2y+l3z−vt)]−2). |