# New Exact Soliton Solutions of the (3+1)-Dimensional Conformable Wazwaz–Benjamin–Bona–Mahony Equation via Two Novel Techniques

Authors: Mohammed K. A. Kaabar; Melike Kaplan; Zailan Siri
Journal: Journal of Function Spaces
 (2021)

---

## Abstract

In this work, the (3+1)-dimensional Wazwaz–Benjamin–Bona–Mahony equation is formulated in the sense of conformable derivative. Two novel methods of generalized Kudryashov and exp−φℵ are investigated to obtain various exact soliton solutions. All algebraic computations are done with the help of the Maple software. Graphical representations are provided in 3D and 2D profiles to show the behavior and dynamics of all obtained solutions at various parameters’ values and conformable orders using Wolfram Mathematica.


---

## Body

## 1. Introduction

Partial differential equations (PDEqs) have attracted a particular interest from researchers in the fields of natural sciences and engineering due to the applicability of these equations in modeling various scientific phenomena in interdisciplinary sciences such as mathematical physics, mechanics, signal and image processing, and chemistry. Most physical systems are not linear; therefore, nonlinear partial differential equations (NLPDEqs), particularly nonlinear evolution equations (NLEEqs) (see [1]), have inspired researchers to investigate the existence of exact solutions for such equations. Finding new exact solutions for NLPDEqs can significantly provide a good interpretation for the physical meaning and dynamics of these equations. Therefore, several research studies have recently been done on developing new methods for solving NLPDEqs exactly. Some of the most notable methods that have been applied to solve some interesting NLPDEqs are the methods of modified simple equation and extended simplest equation to solve (4+1)-dimensional nonlinear Fokas equation [2], the methods of exp-function and exp−Φℵ-expansion to solve Vakhnenko-Parkes equation [3], and the method of generalized Kudryashov to solve nonlinear Jaulent-Miodek hierarchy and (2+1)-dimensional Calogero-Bogoyavlenskii-Schiff equations [4]. To solve nonlinear integrable equations, a novel technique, known as Hirota bilinear method, was proposed by Ma and Ma et al. in [5, 6] in order to obtain new lump solutions for the investigated equations [7].One of the most interesting NLPDEqs is the Benjamin-Bona-Mahony equation (BBMEq), an extension of the Korteweg-de Vries equation (KdVEq), which is basically a model that represents the unidirectional propagation of long waves with small amplitude on the hydromagnetic and acoustic waves’ surface in shallow water channel [8, 9]. Consider the following (3+1)-dimensional modified form of BBMEq:(1)Ψt+Ψx+Ψ2Ψy−Ψyzt=0.Equation (1) was first proposed by Wazwaz [10] by formulating a new three-dimensional modified version of BBMEqs, known as the Wazwaz-Benjamin-Bona-Mahony equation (WBBMEq), via coupling or various generalized contexts or as combination of both of them. Higher dimensional problems with various applications can be described by the WBBMEq with its spatial and temporal variables [8, 9]. Therefore, there is a significant need to find exact solutions for WBBMEq to interpret their physical meaning and dynamics.Fractional differential equations (FDiEqs) are generalized forms of integer-order differential equations. FDiEqs have attracted the interests of researchers due to the ability of these equations in modeling various scientific phenomena better than the integer-order ones (refer to [11–13]). The behavior of some physical systems can be interpreted better than the integer-order ones due to the nonlocality of fractional derivatives, and many systems have memory effects. There are many definitions of fractional derivatives, but the most common ones are Riemann-Liouville and Caputo fractional derivatives where the properties of linearity are commonly shared among these derivatives, but other properties such as product rule, constant, chain rule, and quotient rule are not satisfied. FDiEqs are considered as a powerful tool in modeling scenarios (see [14–16]), but this tool comes with various challenges when dealing with FDiEqs due to the difficulty of obtaining exact or analytical solutions where the solutions can be very complicated or impossible to obtain them for certain cases. As a result, to overcome some of challenges associated with nonlocal fractional derivatives, a new generalized fractional derivative of local-type, named conformable derivative, was proposed by Khalil et al. which is basically a stretch for the usual limit-based definition where all usual derivative’s properties are satisfied [17]. Many research studies have discussed the mathematical analysis and applications of conformable derivative such as conformable Laplace’s equation [18] and generalized conformable mean value theorems [19].Seadawy et al. [20] and Bilal et al. [21] investigated various soliton solutions for conformable WBBMEq using the methods of simple ansatz and generalized exponential rational function, respectively.Inspired by all above studies, this work is mainly aimed at obtaining new exact solitary solutions for a version of WBBMEq formulated in the sense of conformable derivative (ComD) with the help of two novel techniques: generalized Kudryashov method and exp−φℵ method. The general fractional formulation of WBBMEq can be expressed as follows:(2)DtζΨ+DxζΨ+DyζΨ−Dxzt3ζΨ=0,where Dζ is the fractional operator of order: ζ∈0,1. The exact solutions of Eq. (2) have been investigated in some research works using the G′/G-expansion method [9], modified simple equation method [8], and Riccati-Bernoulli Sub-ODE method [8]. However, according to the best of our knowledge, none of previous research works has investigated Eq. (2) in the context of ComD via the methods of generalized Kudryashov and exp−φℵ.Therefore, all results in our work are new and worthy.This work is divided into the following sections: Essential notions about conformable derivative are presented, and the methodology of our two proposed approaches is discussed in Section2. In Section 3, the main results of our work are presented. In Section 4, the graphical comparisons of our obtained exact solutions are represented in both 2D and 3D plots for various values of parameters and ζ. The conclusion of our work is provided in Section 5.


## 2. Fundamental Preliminaries and Methodology

Some important notions about conformable derivative are introduced in this section. In addition, the methodology of two different approaches, namely, generalized Kudryashov method and exp−φℵ method, is also described, respectively. The conformable derivative can be defined as follows:Definition 1.
Given a functionΨ:0,∞⟶R∋∀t>0, the conformable derivative of order ζ∈0,1 of Ψ can be expressed as follows:(3)DtζΨt=limΩ⟶0Ψt+Ωt1−ζ−ΨtΩ.LetΨ be ζ-differentiable in some 0,ϱ, where ϱ>0, and the limit of DtζΨt exists as t⟶0+; then, from this definition, we obtain the following: (4)DtζΨ0=limt⟶0+DtζΨt.The theorem [17] below shows that Dtζ satisfies all usual limit-based derivative’s properties as follows:Theorem 2.
Forζ∈0,1, let functions Ψ and Φ be ζ-differentiable at a point t; then, we have the following:(a)
DtζΨtΦt=ΨtDtζΦt+ΦtDtζΨt(b)
DtζmΨt+wΦt=mDtζΨt+wDtζΦt,∀m,w∈R(c)
DtζΨt/Φt=ΦtDtζΨt−ΨtDtζΦt/Φ2t(d)
Dtζtk=ktk−ζ,∀k∈R(e)
IfΨt is assumed to be a differentiable function, thenDtζΨt=t1−ζdΨ/dt(f)
Dtζν=0,∀constantfunctionsΨt=νThe methodology of the generalized Kudryashov method (GKuM) and exp−φℵ method (ExpM) can be presented as follows:Consider the following form of NLEEq, with 4 independent variables:x,y,z, and t, formulated generally in the sense of fractional derivative:(5)TΨ,DtζΨ,DxζΨ,DyζΨ,DzζΨ,Dt2ζΨ,Dx2ζΨ,Dy2ζΨ,Dz2ζΨ,⋯=0;0<ζ≤1,where Ψ=Ψx,y,z,t is an unknown function and T is a polynomial of Ψ and its partial derivatives in which all of the nonlinear terms and highest-order derivatives are included in Eq. (5). First of all, to solve Eq. (5), we use the following traveling wave transformations for ComD:For ComD,(6)Ψx,y,z,t=Ψℵ;ℵ=pxζζ+qyζζ+γzζζ−δtζζ,where p,q,γ, and δ are all constants with the condition: p,q,γ,δ≠0, and δ is the wave speed.According to the above transformations, Eq. (5) is reduced to the following ODE. (7)LΨ,Ψ′,Ψ′′,Ψ′′′,⋯=0.The derivative with respect toℵ is represented by a prime. Equation (7) should be integrated term by term one or more times.


### 2.1. The Generalized Kudryashov Method

From GKuM, the obtained solution for the reduced equation is constructed by a polynomial inℏℵ as [4, 22]:(8)Ψℵ=∑k=0Jpkℏkℵ∑l=0Wqlℏlℵ,where pkk=0,1,⋯,j,bll=0,1,⋯,w are constants which are needed to be determined ∍pJ≠0,qW≠0, and L=Lℵ is the solution of the following equation:(9)dℏdℵ=ℏ2ℵ−ℏℵ.The solution of Eq. (9) can be expressed as follows:(10)ℏℵ=11+I1eℵ,where I1 is an integration constant.According to the homogeneous balance principle (HBPrp), the positive integers:J and W in Eq. (8) can be obtained with the help of Eq. (7). In addition, a polynomial, ℏ, can be determined by the substitution of Eq. (8) into Eq. (7) along with Eq. (9). Now, by equating all the coefficients of polynomial ℏ to 0 in order to construct a system of algebraic equations, this system is solved with the aid of the computer software such as Maple and Wolfram Mathematica in order to find the values of pkk=0,1,⋯,j,qll=0,1,⋯,w. At the end, all soliton-type solutions of the reduced Eq. (7) can be found by the substitution of these obtained values and Eq. (9) into Eq. (8).


### 2.2. The exp−φℵ Method

From ExpM [3], the obtained solution for the reduced equation is constructed by a polynomial in exp−Φℵ as follows:(11)Ψℵ=∑j=0wpjexp−ϕℵj,where pjpw≠0 are constants which are needed to be found, and φℵ satisfies the following auxiliary ODE:(12)φ′ℵ=exp−φℵ+ϑexpφℵ+χ.Note that Eq. (12) has distinct solutions which are expressed as follows:Case 1.
Whenχ2−4ϑ>0 and ϑ≠0, the hyperbolic function solutions are expressed as follows: (13)ϕ1ℵ=ln−χ2−4ϑtanhχ2−4ϑ/2ℵ+I−χ2ϑ.Case 2.
Whenχ2−4ϑ<0 and ϑ≠0, the trigonometric function solutions are expressed as follows:(14)ϕ2ℵ=ln4ϑ−χ2tan4ϑ−χ2/2ℵ+C−χ2ϑ.Case 3.
Whenχ2−4ϑ>0,ϑ=0 and χ≠0, the hyperbolic function solutions are expressed as follows:(15)ϕ3ℵ=−lnχcoshχℵ+I+sinhχℵ+I−1.Case 4.
Whenχ2−4ϑ=0,ϑ≠0 and χ≠0, the rational function solutions are expressed as follows:(16)ϕ4ℵ=ln−2χℵ+I+2χ2ℵ+I.Case 5.
Whenχ2−4ϑ=0,ϑ=0 and χ=0, we have the following:(17)φ5ℵ=lnℵ+I.From the above cases, the integration constant is represented byI. By the substitution of Eq. (11) into the reduced Eq. (7) and collecting all terms together that are in the same order of exp−φℵjj=0,1,2,⋯, the polynomial in terms of exp−φℵ is verified. Then, by equating all coefficients to 0, a set of algebraic equations is constructed for pjj=0,1,..m,χ,δ and ϑ. With the help of Maple, Walfram Mathematica, or MATLAB, the system can be solved to obtain diverse exact solutions for Eq. (6).


## 2.1. The Generalized Kudryashov Method

From GKuM, the obtained solution for the reduced equation is constructed by a polynomial inℏℵ as [4, 22]:(8)Ψℵ=∑k=0Jpkℏkℵ∑l=0Wqlℏlℵ,where pkk=0,1,⋯,j,bll=0,1,⋯,w are constants which are needed to be determined ∍pJ≠0,qW≠0, and L=Lℵ is the solution of the following equation:(9)dℏdℵ=ℏ2ℵ−ℏℵ.The solution of Eq. (9) can be expressed as follows:(10)ℏℵ=11+I1eℵ,where I1 is an integration constant.According to the homogeneous balance principle (HBPrp), the positive integers:J and W in Eq. (8) can be obtained with the help of Eq. (7). In addition, a polynomial, ℏ, can be determined by the substitution of Eq. (8) into Eq. (7) along with Eq. (9). Now, by equating all the coefficients of polynomial ℏ to 0 in order to construct a system of algebraic equations, this system is solved with the aid of the computer software such as Maple and Wolfram Mathematica in order to find the values of pkk=0,1,⋯,j,qll=0,1,⋯,w. At the end, all soliton-type solutions of the reduced Eq. (7) can be found by the substitution of these obtained values and Eq. (9) into Eq. (8).


## 2.2. The exp−φℵ Method

From ExpM [3], the obtained solution for the reduced equation is constructed by a polynomial in exp−Φℵ as follows:(11)Ψℵ=∑j=0wpjexp−ϕℵj,where pjpw≠0 are constants which are needed to be found, and φℵ satisfies the following auxiliary ODE:(12)φ′ℵ=exp−φℵ+ϑexpφℵ+χ.Note that Eq. (12) has distinct solutions which are expressed as follows:Case 1.
Whenχ2−4ϑ>0 and ϑ≠0, the hyperbolic function solutions are expressed as follows: (13)ϕ1ℵ=ln−χ2−4ϑtanhχ2−4ϑ/2ℵ+I−χ2ϑ.Case 2.
Whenχ2−4ϑ<0 and ϑ≠0, the trigonometric function solutions are expressed as follows:(14)ϕ2ℵ=ln4ϑ−χ2tan4ϑ−χ2/2ℵ+C−χ2ϑ.Case 3.
Whenχ2−4ϑ>0,ϑ=0 and χ≠0, the hyperbolic function solutions are expressed as follows:(15)ϕ3ℵ=−lnχcoshχℵ+I+sinhχℵ+I−1.Case 4.
Whenχ2−4ϑ=0,ϑ≠0 and χ≠0, the rational function solutions are expressed as follows:(16)ϕ4ℵ=ln−2χℵ+I+2χ2ℵ+I.Case 5.
Whenχ2−4ϑ=0,ϑ=0 and χ=0, we have the following:(17)φ5ℵ=lnℵ+I.From the above cases, the integration constant is represented byI. By the substitution of Eq. (11) into the reduced Eq. (7) and collecting all terms together that are in the same order of exp−φℵjj=0,1,2,⋯, the polynomial in terms of exp−φℵ is verified. Then, by equating all coefficients to 0, a set of algebraic equations is constructed for pjj=0,1,..m,χ,δ and ϑ. With the help of Maple, Walfram Mathematica, or MATLAB, the system can be solved to obtain diverse exact solutions for Eq. (6).


## 3. Solutions of the (3+1)-Dimensional Conformable WBBMEq

Exact soliton solutions of the proposed Eq. (2) are obtained in this section via GKuM and ExpM.(18)pγδΨ′′′+qΨ3′+−γ+pΨ′=0.By integrating Eq. (18) once with respect to ℵ, we obtain the following:(19)pγδΨ″+qΨ3−δ+pΨ=0.


### 3.1. Exact Solutions via the Generalized Kudryashov Method

According to the HBPrp, it is obvious to haveJ=W+1. Let us set W=1, and we obtain J=2. Thus, the solution can be written as follows: (20)Ψℵ=p0+p1ℏ+p2ℏ2q0+q1ℏ,where ℏ=ℏℵ is the solution of Eq. (9). As a result, by substituting Eq. (20) into Eq. (19) and using Eq. (9), we obtain system of algebraic equations by equating all coefficients of the functions ℏ0,ℏ1,ℏ2,ℏ3,ℏ4,ℏ5,ℏ6 to 0. Now, p0,p1,p2,q0, and q1 are all parameters.(21)ℏ6:qp23+2pγδp2q12=0,ℏ5:3qp1p22−3pγδp2q12+6pγδp2q0q1=0,ℏ4:−9pγδp2q0q1−δp2q12+pp2q12+6pγδp2q02+pγδp2q12+3qp0p22+3qp12p2=0,ℏ3:−δp1q12+2pγδp1q02+qp13−10pγδp2q02+pγδℏ1q0q1−pγδq12p0+pp1q12−2δp2q0q1+3pγδp2q0q1+6qp0p1p2+2pp2q0q1−2pγδq1p0q0=0,ℏ2:−2δp1q0q1−pγδp1q0q1+pp0q12+3qp02p2−δp0q12+3qp0p12−3pγδp1q02+2pp1q0q1+pp2q02−δp2q02+4pγδp2q02+pγδq12p0+3pγδq1p0q0=0,ℏ1:−pγδq1p0q0+2pp0q0q1+pγδp1q02−δp1q02+3qp02p1−2δp0q0q1+pp1q02=0,ℏ0:pp0q02+qp03=0.From the above set of algebraic equations, various cases are presented as follows:Case 1.
(22)p0=0,p1=±pq1−γ2q+pqγ,p2=±2pγq1pγ+2q−γ/2q+pqγ,q0=0,q1=q1,δ=2ppγ+2.Then, by the substitution of the obtained values into Eq. (20) with Eq. (10), the soliton-type solutions of the following WBBMEq in the sense of ComD are as follows: (23)Ψ1x,y,z,t=±1−I1coshℵ+sinhℵpγ−γ/qpγ+2qpγ+21+I1coshℵ+I1sinhℵ,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγ+2tς/ς for ComD. I1 is an arbitrary constant.Case 2.
(24)p0=0,p1=±2q0p2γ−q+pqγ,p2=∓2q0p2γ−q+pqγ,q0=q0,q1=−2q0,δ=p1−pγ.Then, by the substitution of the obtained values into Eq. (20) with Eq. (10), the soliton-type solutions of the following WBBMEq in the sense of ComD are as follows:(25)Ψ2x,y,z,t=±2pI12γ/q−1+pγcoshℵ+sinhℵ−1+I122coshℵsinhℵ+2cosh2ℵ−1,where ℵ=pxς/ς+qyς/ς+γzς/ς−p/1−pγtς/ς for ComD. I1 is a constant.


### 3.2. Exact Solutions via the exp−φℵ Method

The ExpM is a very helpful technique to construct the solution of Eq. (7) in the following form:(26)Ψξ=p0+p1exp−φℵ.Let us substitute Eq. (19) into Eq. (7) and collect the coefficient of each power of exp−Φℵj. Now, a set of algebraic equations of p,q,γ,δ,p0,p1,χ and ϑ is obtained by equating all coefficients to 0.(27)exp3ℵ:qp03δp0+pγδp1ϑχ+pp0=0,exp2ℵ:3qp02p1−δp1+2pγδp1ϑ+pγδp1χ2+pp1=0,expℵ:3pγδp1χ+3qp0p12=0,exp0ℵ:2pγδp1+qp13=0.With the help of software programs such as Maple or Mathematica, the solution is obtained as follows:(28)p0=±χpγ−pqγχ2−2q+4pqγϑ,p1=±2pγ−pqγχ2−2q+4pqγϑ,δ=2ppγχ2+2−4pγϑ.From all the above obtained values, the technique’s algorithm, and its auxiliary equations, different cases for the conformable WBBMEq are given as follows:Case 1.
Whenχ2−4ϑ>0 and ϑ≠0, the hyperbolic function solutions are expressed as follows:(29)Ψ1x,y,z,t=±pγ/q−pγχ2−2+4pγϑχ2+χtanhχ2−4ϑ/2ℵ+Iχ2−4ϑ−4ϑχ+tanhχ2−4ϑ/2ℵ+Iχ2−4ϑ,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 2.
Whenχ2−4ϑ<0 and ϑ≠0, the trigonometric function solutions are expressed as follows:(30)Ψ2x,y,z,t=±pγ/q−pγχ2−2+4pγϑ−χ2+χtan4ϑ−χ2/2ℵ+I4ϑ−χ2+4ϑ−χ+tan4ϑ−χ2/2ℵ+I4ϑ−χ2,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 3.
Whenχ2−4ϑ>0,ϑ=0 and χ≠0, the hyperbolic function solutions are expressed as follows: (31)Ψ3x,y,z,t=±χpγ/q−pγχ2−2+4pγϑcoshχℵ+I+sinhχℵ+I+1coshχℵ+I+sinhχℵ+I−1,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 4.
Whenχ2−4ϑ=0,ϑ≠0 and χ≠0, the rational function solutions are expressed as follows:(32)Ψ4x,y,z,t=±2χpγ/q−pγχ2−2+4pγϑχℵ+I+2,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 5.
Finally, whenχ2−4ϑ=0,ϑ=0 and χ=0, we have the following: (33)Ψ5x,y,z,t=±γ/q−pγχ2−2+4pγϑpχℵ+I+2ℵ+I,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.


## 3.1. Exact Solutions via the Generalized Kudryashov Method

According to the HBPrp, it is obvious to haveJ=W+1. Let us set W=1, and we obtain J=2. Thus, the solution can be written as follows: (20)Ψℵ=p0+p1ℏ+p2ℏ2q0+q1ℏ,where ℏ=ℏℵ is the solution of Eq. (9). As a result, by substituting Eq. (20) into Eq. (19) and using Eq. (9), we obtain system of algebraic equations by equating all coefficients of the functions ℏ0,ℏ1,ℏ2,ℏ3,ℏ4,ℏ5,ℏ6 to 0. Now, p0,p1,p2,q0, and q1 are all parameters.(21)ℏ6:qp23+2pγδp2q12=0,ℏ5:3qp1p22−3pγδp2q12+6pγδp2q0q1=0,ℏ4:−9pγδp2q0q1−δp2q12+pp2q12+6pγδp2q02+pγδp2q12+3qp0p22+3qp12p2=0,ℏ3:−δp1q12+2pγδp1q02+qp13−10pγδp2q02+pγδℏ1q0q1−pγδq12p0+pp1q12−2δp2q0q1+3pγδp2q0q1+6qp0p1p2+2pp2q0q1−2pγδq1p0q0=0,ℏ2:−2δp1q0q1−pγδp1q0q1+pp0q12+3qp02p2−δp0q12+3qp0p12−3pγδp1q02+2pp1q0q1+pp2q02−δp2q02+4pγδp2q02+pγδq12p0+3pγδq1p0q0=0,ℏ1:−pγδq1p0q0+2pp0q0q1+pγδp1q02−δp1q02+3qp02p1−2δp0q0q1+pp1q02=0,ℏ0:pp0q02+qp03=0.From the above set of algebraic equations, various cases are presented as follows:Case 1.
(22)p0=0,p1=±pq1−γ2q+pqγ,p2=±2pγq1pγ+2q−γ/2q+pqγ,q0=0,q1=q1,δ=2ppγ+2.Then, by the substitution of the obtained values into Eq. (20) with Eq. (10), the soliton-type solutions of the following WBBMEq in the sense of ComD are as follows: (23)Ψ1x,y,z,t=±1−I1coshℵ+sinhℵpγ−γ/qpγ+2qpγ+21+I1coshℵ+I1sinhℵ,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγ+2tς/ς for ComD. I1 is an arbitrary constant.Case 2.
(24)p0=0,p1=±2q0p2γ−q+pqγ,p2=∓2q0p2γ−q+pqγ,q0=q0,q1=−2q0,δ=p1−pγ.Then, by the substitution of the obtained values into Eq. (20) with Eq. (10), the soliton-type solutions of the following WBBMEq in the sense of ComD are as follows:(25)Ψ2x,y,z,t=±2pI12γ/q−1+pγcoshℵ+sinhℵ−1+I122coshℵsinhℵ+2cosh2ℵ−1,where ℵ=pxς/ς+qyς/ς+γzς/ς−p/1−pγtς/ς for ComD. I1 is a constant.


## 3.2. Exact Solutions via the exp−φℵ Method

The ExpM is a very helpful technique to construct the solution of Eq. (7) in the following form:(26)Ψξ=p0+p1exp−φℵ.Let us substitute Eq. (19) into Eq. (7) and collect the coefficient of each power of exp−Φℵj. Now, a set of algebraic equations of p,q,γ,δ,p0,p1,χ and ϑ is obtained by equating all coefficients to 0.(27)exp3ℵ:qp03δp0+pγδp1ϑχ+pp0=0,exp2ℵ:3qp02p1−δp1+2pγδp1ϑ+pγδp1χ2+pp1=0,expℵ:3pγδp1χ+3qp0p12=0,exp0ℵ:2pγδp1+qp13=0.With the help of software programs such as Maple or Mathematica, the solution is obtained as follows:(28)p0=±χpγ−pqγχ2−2q+4pqγϑ,p1=±2pγ−pqγχ2−2q+4pqγϑ,δ=2ppγχ2+2−4pγϑ.From all the above obtained values, the technique’s algorithm, and its auxiliary equations, different cases for the conformable WBBMEq are given as follows:Case 1.
Whenχ2−4ϑ>0 and ϑ≠0, the hyperbolic function solutions are expressed as follows:(29)Ψ1x,y,z,t=±pγ/q−pγχ2−2+4pγϑχ2+χtanhχ2−4ϑ/2ℵ+Iχ2−4ϑ−4ϑχ+tanhχ2−4ϑ/2ℵ+Iχ2−4ϑ,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 2.
Whenχ2−4ϑ<0 and ϑ≠0, the trigonometric function solutions are expressed as follows:(30)Ψ2x,y,z,t=±pγ/q−pγχ2−2+4pγϑ−χ2+χtan4ϑ−χ2/2ℵ+I4ϑ−χ2+4ϑ−χ+tan4ϑ−χ2/2ℵ+I4ϑ−χ2,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 3.
Whenχ2−4ϑ>0,ϑ=0 and χ≠0, the hyperbolic function solutions are expressed as follows: (31)Ψ3x,y,z,t=±χpγ/q−pγχ2−2+4pγϑcoshχℵ+I+sinhχℵ+I+1coshχℵ+I+sinhχℵ+I−1,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 4.
Whenχ2−4ϑ=0,ϑ≠0 and χ≠0, the rational function solutions are expressed as follows:(32)Ψ4x,y,z,t=±2χpγ/q−pγχ2−2+4pγϑχℵ+I+2,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.Case 5.
Finally, whenχ2−4ϑ=0,ϑ=0 and χ=0, we have the following: (33)Ψ5x,y,z,t=±γ/q−pγχ2−2+4pγϑpχℵ+I+2ℵ+I,where ℵ=pxς/ς+qyς/ς+γzς/ς−2p/pγχ2+2−4pγϑtς/ς for ComD. I is a constant.


## 4. The Graphical Comparisons of Solutions

To show the dynamics and behavior of our obtained solutions, various soliton solutions in Eqs. (23), (25), (29), (30), (31), (32), and (33) are graphically represented and compared in both 3D and 2D plots in Figures 1–7 for various parameters’ values and ζ. In this work, two techniques are employed in the sense of ComD. Therefore, all our results are new and worthy because studying WBBMEq is very important to investigate various nonlinear scientific phenomena. ComD is a local-type fractional derivative which is a generalized formulation of usual limit-based derivative. Since the measurements in physics are local, this can make ComD suitable for modeling many physical phenomena, and it is also efficient to work with ComD to obtain exact solutions for NLPDEqs although ComD does not have some of the essential properties to be categorized as a fractional derivative. From the authors’ opinion, formulating and solving NLPDEqs in the sense of ComD are always recommended since all local and nonlocal fractional derivatives have both advantages and disadvantages. Therefore, exploring new properties and definitions for all local and nonlocal fractional derivatives are helpful while working on modeling nonlinear scientific phenomena. In addition, studying the ComD’s definition is also interesting because any new mathematical definition deserves to be explored and investigated to show its validity and potentiality for more application.Figure 1
The plots of Eq. (23) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=−1; q=p=I1=1; t=15; y=z=0; ζ=0.50; ζ=0.80; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 2
The plots of Eq. (25) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=−1; q=p=I1=1; t=15; y=z=0; ζ=0.50; ζ=0.80; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 3
The plots of Eq. (29) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=1; ϑ=−1; q=p=χ=I=1; t=15; y=z=0; ζ=0.50; ζ=0.75; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 4
The plots of Eq. (30) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=1; ϑ=1; q=p=χ=I=1; t=15; y=z=0; ζ=0.50; ζ=0.75; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 5
The plots of Eq. (31) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=−1; ϑ=0; q=p=χ=I=1; t=15; y=z=0; ζ=0.50; ζ=0.75; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 6
The plots of Eq. (32) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=1; ϑ=0; q=p=χ=I=1; t=15; y=z=0; ζ=0.50; ζ=0.75; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)Figure 7
The plots of Eq. (33) are represented in 3D in (a), (c), and (e) and in 2D in (b), (d), and (f) for γ=−1; ϑ=0; q=p=I=1; χ=0; t=15; y=z=0; ζ=0.50; ζ=0.75; ζ=0.90 for ComD, respectively.
(a)(b)(c)(d)(e)(f)


## 5. Conclusion

Two novel techniques of generalized Kudryashov and exp−φℵ have been applied in this work to investigate exact soliton solutions of the (3+1)-dimensional conformable WBBMEq in the sense of ComD. The obtained solutions are new which imply that the studied techniques provide efficient results. All algebraic computations in this work have been done using the Maple software. Graphical representations have been provided for all obtained solutions at various parameters’ values and conformable orders with the help of Wolfram Mathematica. All in all, the studied methods can be potentially applied to solve various NLPDEqs that are apparent in many important nonlinear scientific phenomena in physics and engineering. Our results can be further extended in future research works into solving various classes of higher dimensional nonlinear partial differential equations which will provide a major contribution to soliton theory and mathematical physics.



---

*Source: 4659905-2021-07-22.xml*