Next Article in Journal
A Fractal Approach to Nonlinear Topographical Features of Healthy and Keratoconus Corneas Pre- and Post-Operation of Intracorneal Implants
Previous Article in Journal
Possibilities of Use for Fractal Techniques as Parameters of Graphic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Biswas–Milovic Equation in Random Case Study

by
Abdulwahab Almutairi
Department of Mathematics, College of Arts and Sciences Onaizah, Qassim University, Buraidah 51452, Saudi Arabia
Fractal Fract. 2022, 6(11), 687; https://doi.org/10.3390/fractalfract6110687
Submission received: 18 October 2022 / Revised: 4 November 2022 / Accepted: 8 November 2022 / Published: 19 November 2022
(This article belongs to the Section Mathematical Physics)

Abstract

:
We apply two mathematical techniques, specifically, the unified solver approach and the exp ( φ ( ξ ) ) -expansion method, for constructing many new solitary waves, such as bright, dark, and singular soliton solutions via the fractional Biswas–Milovic (FBM) model in the sense of conformable fractional derivative. These solutions are so important for the explanation of some practical physical problems. Additionally, we study the stochastic modeling for the fractional Biswas–Milovic, where the parameter and the fraction parameters are random variables. We consider these parameters via beta distribution, so the mathematical methods that were used in this paper may be called random methods, and the exact solutions derived using these methods may be called stochastic process solutions. We also determined some statistical properties of the stochastic solutions such as the first and second moments. The proposed techniques are robust and sturdy for solving wide classes of nonlinear fractional order equations. Finally, some selected solutions are illustrated for some special values of parameters.

1. Introduction

Nonlinear fractional differential equations (NFDEs) play vital roles in many interesting applications in chemical engineering, fluid mechanics, biology, electromagnetic theory, physics, and others [1,2,3,4]. These equations are powerful instruments for depicting real-world problems more accurately than the classical integer-order equations. Thus, the investigation of solitary wave solutions for NFDEs becomes very useful in scientific research. Recently, many researches have proposed and developed various numerical and analytical methods for solving NFDEs. Shah et al. applied the time-fractional Caputo and Caputo–Fabrizio fractional derivatives to the Chua type nonlinear chaotic system [5]. Alshehry et al. presented the Laplace residual-power-series method (LRPSM), a powerful new technique for solving fractional partial differential equations [6].
For the above reasons, recently, several efficient mathematical approaches have been proposed to obtain solutions of NFDEs, such as the ( G G ) expansion method [7], the fractional sub-equation method [8], the first integral method [9], the tanh-sech method [10], the unified solver method [11], the exponential function method [12] and others [13,14,15,16]. Recently, the dynamical behavior of the exact traveling wave solutions and their phase portrait analysis were extensively studied using the dynamical system theory. Zhu et al. [17] studied the exact traveling wave solutions and bifurcations of the fractional Klein–Gordon equation and the fractional generalized Hirota–Satsuma coupled KdV system. Specifically, they applied this technique for the first time to NFDEs. Based on these powerful results, some other authors followed the same technique, for example, see [18,19,20]. Li et al. [18] studied the dynamical behavior of a time–space fractional Phi-4 equation using the bifurcation method of a planar dynamical system via conformable fractional derivative. Liu et al. [19] investigated the dynamical behavior and bifurcation of solutions of the traveling wave system for a generalized (3 + 1)-dimensional time-fractional gCH-KP equation. Li and Han [20] investigated the bifurcation and new exact solutions for the (2 + 1)-dimensional conformable time-fractional Zoomeron equation.
The Biswas–Milovic equation studied in this paper is given by [21]
i ( q m ) t + a ( q m ) x x + b F ( q 2 ) q m = 0 , i = 1 ,
where x and t are two independent variables, and q is a complex valued function. The coefficients a and b are constants with a b > 0 , and m 1 is a parameter. Moreover, F is a real-valued algebraic function, which is necessary to have smoothness of the complex function F ( q 2 ) : C C . Equation (1) is not integrable, in general. The non-integrability is not necessarily related to the nonlinear term in it. This equation emerges in the study of long-distance optical communications and all-optical ultra fast switching devices. Furthermore, this equation has been indicated to manage the evolution of a wave packet in a weakly nonlinear and dispersive medium and has eventuated diverse fields, such as plasma, nonlinear optics, and water waves.
To highlight numerous complex phenomena in various fields of nature, such as biology, economy, engineering, chemical engineering, signal processing, solid state physics, and electromagnetic theory, it is crucial to take into account random problem effects. Practical considerations require that the stochastic type perturbations be considered. More focus has been placed in recent years on the impact of noise on the spread of these solton solutions. Therefore, we must take into consideration the stochastic effect.
The novelties of this paper are mainly exhibited in four aspects. First, we present a general form of a new fractional Biswas–Milovic equation (FBM). Second, we use a new method, the so-called unified solver method [22], in order to solve the FBM. Many other nonlinear models developing in applied science and new physics can be solved using the solver as a box solver. In comparison to previous methods, this solver has some advantages, such as avoiding difficult and time-consuming calculations and presenting important solutions explicitly. Moreover, we use the exp-function method [23,24] in order to solve the proposed equation. Indeed, this method presents powerful solutions in vital applications in natural science [25,26]. Third, we obtain new types of exact analytical solutions. Comparing our results with other results, one can see that our results are new and more extensive. Fourth, we study the stochastic modeling for the fractional Biswas–Milovic through the parameter where the fraction parameters are random variables. Thus, we must discuss our fractional problem under beta distribution [27].
The remainder of this paper is arranged as follows. In Section 2, some preliminaries, notions of local fractional calculus, and the beta random distribution method are introduced. Section 3 presents the description of the two methods applied in this work, namely, the unified solver method and the exp-function method. In Section 4, some exact solutions of the fractional BM equation are presented, using the two methods. Stochastic study with beta distribution is studied in Section 5. Conclusions are in Section 6.

2. Preliminaries

Here, we give some fundamental notions of fractional calculus theory, which turn out to be very useful to complete this article in a unified way. There are various types of fractional derivatives [2], such as G r u ¨ n w a l d –Letnikov, Caputo local fractional derivative [28], Abel–Riemann fractional derivative [29] and He’s fractional derivative [15,16]. Khalil et al. [30] presented a new vital definition of fractional derivatives, the so-called conformable fractional derivative. This definition preserves various advantages that cannot be satisfied by the known fractional derivatives, such as Rolle’s theorem, mean value theorem, the product of two functions, and the chain rule [31]. The conformable fractional derivative has attracted important attention due to its simplicity. Thus, many studies has been done on it by many scientists. First, we introduce some properties of the conformable fractional derivative. Second, we give the beta random distribution.
Definition 1
([30]). Let a function ϕ : ( 0 , ) R , then the conformable fraction derivative of ϕ of order α is
D t α ( ϕ ) ( t ) = lim ι 0 ϕ ( t + ι t 1 α ) ϕ ( t ) ι , t > 0 , 0 < α 1 .
The conformable fractional derivative satisfies:
(i)
D t α ( a χ + b ϕ ) = a D t α ( χ ) + b D t α ( ϕ ) , a , b R ,
(ii)
D t α ( t n ) = n t n α , n R ,
(iii)
D t α ( χ ϕ ) = χ D t α ( ϕ ) + ϕ D t α ( χ ) ,
(vi)
D t α ( χ ϕ ) = ϕ D t α ( χ ) χ D t α ( ϕ ) ϕ 2 ,
(v)
If ϕ is differentiable, thus D t α ( ϕ ) ( t ) = t 1 α d ϕ d t .
Theorem 1
([30]). Let χ , ϕ : ( 0 , ) R be differentiable and also α-differentiable, then:
D t α ( χ ϕ ) ( t ) = t 1 α ϕ ( t ) ψ ( ϕ ( t ) ) .

Beta Distribution

If we have a real random variable X defined on the probability space ( Ω , F , P ) that has the probability density function f ( x ) , we can defined the statistical moments as follows:
Remark 1.
1.
The first moment:
E X = x f x d x
2.
The second moment:
E X 2 = x 2 f x d x .
Therefore, if X has a beta distribution, then the probability density function defined as
f ( x ) = ( x a ) α 1 ( b x ) β 1 B ( α , β ) ( b a ) α + β 1 a x b ; α , β > 0
is called a beta random variable. Here, α and β are the shape parameters, a and b are the lower and upper bounds, respectively, of the distribution. In addition, B ( α , β ) is the beta function, and if we have a = 0 and b = 1 , it is called the standard beta distribution.
Additionally, E X = α α + β , where, E denotes the expectation value operator. In this work, we will deal with the standard beta distribution.
Any NFDEs in two independent variables x and t can be expressed as follows:
Λ ( ϕ , D t α ϕ , D x α ϕ , D t α D x α ϕ , D x α D x α ϕ , . . . ) = 0 ,
1 α > 0 .
Utilizing the traveling wave transformation
ϕ ( x , t ) = Φ ( ξ ) , ξ = x α α + w t α α
converts Equation (5) to the following ODE:
Λ ( Φ , Φ , Φ , Φ , . . . ) = 0 .
If we have any NFDEs in two independent variables x and t, then we develop it only in the beta random distribution because the range of the beta random variable is 1 α > 0 . Therefore, the random case of (5) is given when α has a beta random distribution. We use the beta random distribution traveling wave transformation (6), where k, λ are non zero constants and α has beta random distribution. Therefore, we can complete the beta random distribution methods as in the deterministic α where 0 < α 1 . In this paper, the beta random distribution method is considered through two methods. When we use the beta random variable as in the random problem itself or when we use beta random distribution traveling wave transformation, the proposed methods, namely the unified solver and the exp ( φ ( ξ ) ) -expansion methods, are called random methods. Additionally, the solutions produced by these random methods are stochastic process solutions.

3. Description of the Methods

We briefly give the unified solver [22] and the exp [ φ ( ξ ) ] -expansion techniques [23,24].
Various NFDEs (5) reduce to the following ODE:
Λ 1 Φ + Λ 2 Φ 3 + Λ 3 Φ = 0 ,
where Λ 1 , Λ 2 , and Λ 3 are constants that depend on the proposed equation’s constants and the wave transformations’ speed.

3.1. Unified Solver Method

In light of the unified solver [22], the solutions of Equation (8) are:
(i)
Rational solutions: (when Λ 3 = 0 )
Φ 1 , 2 ( x , t ) = Λ 2 2 Λ 1 ξ + ε 1 .
(ii)
Trigonometric solutions: (when Λ 3 Λ 1 < 0 )
Φ 3 , 4 ( x , t ) = ± Λ 3 Λ 2 t a n Λ 3 2 Λ 1 η + ε
and
Φ 5 , 6 ( x , t ) = ± Λ 3 Λ 2 c o t Λ 3 2 Λ 1 η + ε .
(iii)
Hyperbolic solutions: (when Λ 3 Λ 1 > 0 )
Φ 7 , 8 ( x , t ) = ± Λ 3 Λ 2 t a n h Λ 3 2 Λ 1 η + ε
and
Φ 9 , 10 ( x , t ) = ± Λ 3 Λ 2 c o t h Λ 3 2 Λ 1 η + ε .
Here, ε is an arbitrary constant.

3.2. The Exp [ φ ( ξ ) ] -Expansion Method

In view of the exp [ φ ( ξ ) ] -expansion approach [23,24], the solution of Equation (8) is rewritten in the following polynomial form of e x p [ φ ( ξ ) ] :
Φ ( ξ ) = B 0 + B 1 exp [ φ ( ξ ) ]
B 0 and B 1 0 are constants and φ ( ξ ) satisfies
φ ( ξ ) = exp [ φ ( ξ ) ] + ν exp [ φ ( ξ ) ] + λ .
The solutions of Equation (15) are:
1.
At λ 2 4 ν > 0 & ν 0 ,
φ ( ξ ) = l n λ 2 4 ν t a n h λ 2 4 ν 2 ξ + C λ 2 ν ,
2.
At λ 2 4 ν < 0 & ν 0 ,
φ ( ξ ) = l n 4 ν λ 2 t a n 4 ν λ 2 2 ξ + C λ 2 ν ,
3.
At λ 2 4 ν > 0 & λ 0 & ν = 0 ,
φ ( ξ ) = l n λ exp λ ξ + C 1 ,
4.
At λ 2 4 ν = 0 & λ 0 & ν 0 ,
φ ( ξ ) = l n 2 λ ξ + C + 2 λ 2 ξ + C ,
5.
At λ 2 4 ν = 0 & λ = 0 & ν = 0 ,
φ ( ξ ) = l n ξ + C .
Here, C is an arbitrary constant.
Setting Equation (14) and Equation (15) into Equation (8) and adding all phrases with the same power exp m φ ( ξ ) , m = 0 , 1 , 2 , 3 . Then, putting them with zero yields algebraic equations. Solving these equations provides the values of B 0 , B 1 . Then, we obtain the solutions (14) that produce the exact solutions of Equation (7).

4. Application

In this section, we are concerned with the fractional type of BM equation [32]
i D t α q m + γ D x 2 α q m + δ q 2 q m = 0 ,
where 1 α > 0 , m 1 . There are many approaches that have been applied to provide solutions for Equation (21), such as the four trigonometric analytical methods [33], Adomian decomposition method [34], the extended fractional sinh–Gordon equation expansion approach [32], etc. Indeed, most standard papers considered the same methods to study the fractional differential equations in a deterministic sense. In contrast to these papers, we considered the stochastic modeling for the fractional Biswas–Milovic through the parameter, and the fraction parameters are random variables. In probability theory, this means that the only random distribution for these parameters is the beta distribution, so the mathematical methods used in this paper may be called random methods, and the exact solutions found using these methods may be called stochastic process solutions. In addition, we determined some of statistical properties of the stochastic solutions as the first and the second moments. To perform this procedure, we apply the unified solver method and the exp-function method. We consider the traveling wave transformation:
q ( x , t ) = u ( ξ ) e i ψ , ξ = x α α ν t α α , ψ = k x α α + w t α α + θ ,
where ν , k, θ , and w denote the speed of the wave, frequency, phase constant, and wave number, respectively. Substituting (22) into (21) gives [32]
a u + b u 3 ( w + a k 2 ) u = 0
and
ν = 2 a k .
Now we apply the unified solver method and exp-function methods for Equation (24).

4.1. On Solving Equation (21) Using the Unified Solver Method

In view of the unified solver approach introduced in [22], the solutions for Equation (24) are:

4.1.1. Rational Solutions

The rational solution of Equation (24) is
u 1 ( x , t ) = b 2 a x α α ν t α α + μ 1 .
Using Equations (22) and (25), the solutions of Equation (21) take the form:
q 1 ( x , t ) = e i k x α α + w t α α + θ b 2 a x α α ν t α α + μ 1 .

4.1.2. Trigonometric Solutions

The trigonometric solutions of Equation (24) are
u 2 , 3 ( x , t ) = ± ( w + a k 2 ) b t a n ( w + a k 2 ) 2 a x α α ν t α α + μ
and
u 4 , 5 ( x , t ) = ± ( w + a k 2 ) b c o t ( w + a k 2 ) 2 a x α α ν t α α + μ .
Using Equations (22), (28) and (43), the solutions of Equation (21) take the forms:
q 2 , 3 ( x , t ) = ± ( w + a k 2 ) b e i k x α α + w t α α + θ t a n ( w + a k 2 ) 2 a x α α ν t α α + μ
and
q 4 , 5 ( x , t ) = ± ( w + a k 2 ) b e i k x α α + w t α α + θ c o t ( w + a k 2 ) 2 a x α α ν t α α + μ .

4.1.3. Hyperbolic Solutions

The hyperbolic solutions of Equation (24) are
u 6 , 7 ( x , t ) = ± w + a k 2 b t a n h ( w + a k 2 ) 2 a x α α ν t α α + μ
and
u 8 , 9 ( x , t ) = ± w + a k 2 b c o t h ( w + a k 2 ) 2 a x α α ν t α α + μ .
Using Equations (22), (31) and (32), the solutions of Equation (21) take the forms:
q 6 , 7 ( x , t ) = ± w + a k 2 b e i k x α α + w t α α + θ t a n h ( w + a k 2 ) 2 a x α α ν t α α + μ
and
q 8 , 9 ( x , t ) = ± w + a k 2 b e i k x α α + w t α α + θ c o t h ( w + a k 2 ) 2 a x α α ν t α α + μ .

4.2. On Solving Equation (21) Using the Exp ( φ ( ξ ) ) -Expansion Method

In view of the exp ( φ ( ξ ) ) -expansion technique [23,24], the solutions of Equation (24) have the following solution:
u = A 0 + A 1 exp ( φ ) ,
where A 0 and A 1 are constants to be determined, such that A 1 0 . It is easy to see that
u = A 1 2 exp ( 3 φ ) + 3 λ exp ( 2 φ ) + ν λ + 2 ν exp ( φ ) + λ 2 exp ( φ ) ,
u 3 = A 1 3 exp ( 3 φ ) + 3 A 0 A 1 2 exp ( 2 φ ) + 3 A 0 2 A 1 exp ( φ ) + A 0 3
Substituting u, u , u 3 into Equation (24) and then equating the coefficients of exp ( φ ) to zero gives a system of algebraic equations. Solving this system yields
u ( ξ ) = ± c 2 b λ + 2 exp ( φ ( ξ ) ) .
Hence, the solutions of Equation (24) are:
Case 1. At λ 2 4 ν > 0 , ν 0 ,
u ˜ 1 , 2 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν tanh λ 2 4 ν 2 x α α ν t α α + K + λ ,
u ˜ 3 , 4 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν coth λ 2 4 ν 2 x α α ν t α α + K + λ .
Using Equations (22), (39) and (40), the solutions of Equation (21) take the forms:
q ˜ 1 , 2 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν e i k x α α + w t α α + θ tanh λ 2 4 ν 2 x α α ν t α α + K + λ ,
q ˜ 3 , 4 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν e i k x α α + w t α α + θ coth λ 2 4 ν 2 x α α ν t α α + K + λ .
Case 2. At λ 2 4 ν < 0 , ν 0 ,
u ˜ 5 , 6 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 tan 4 ν λ 2 2 x α α ν t α α + K λ ,
u ˜ 7 , 8 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 cot 4 ν λ 2 2 x α α ν t α α + K λ .
Using Equations (22), (43) and (44), the solutions of Equation (21) take the forms:
q ˜ 5 , 6 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 e i k x α α + w t α α + θ tan 4 ν λ 2 2 x α α ν t α α + K λ ,
q ˜ 7 , 8 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 e i k x α α + w t α α + θ cot 4 ν λ 2 2 x α α ν t α α + K λ .
Case 3. At λ 2 4 ν > 0 , ν = 0 , λ 0
u ˜ 9 , 10 ( ξ ) = ± c 2 b 1 + 2 λ exp λ x α α ν t α α + K 1 .
Using Equations (22) and (47), the solutions of Equation (21) take the forms:
q ˜ 9 , 10 ( ξ ) = ± c 2 b e i k x α α + w t α α + θ 1 + 2 λ exp λ x α α ν t α α + K 1 .

5. Discussion in Some Stochastic Cases under Beta Random Distribution

In this section, we consider Equation (21) in a stochastic sense, where α is a beta random variable. Namely, we study the statistical properties for the stochastic processes solutions as follows:

5.1. Stochastic Solutions of (21) via the Unified Solver Method

Here, we can write the stochastic solutions of problem (21) by using the unified solver method to show the statistical properties of the stochastic processes solutions. As in the deterministic case, we find the random relations as follows:

5.1.1. Rational Stochastic Solutions

s q 1 ( x , t ) = e i k x α α + w t α α + θ b 2 a x α α ν t α α + μ 1 .
The mean value E [ s q 1 ] and the second moment E [ s q 1 2 ] are depicted in Figure 1a,b, respectively.

5.1.2. Trigonometric Stochastic Solutions

s q 2 , 3 ( x , t ) = ± ( w + a k 2 ) b e i k x α α + w t α α + θ t a n ( w + a k 2 ) 2 a x α α ν t α α + μ
The mean value E [ s q 3 ] and the second moment E [ s q 3 2 ] are depicted in Figure 2a,b, respectively.
s q 4 , 5 ( x , t ) = ± ( w + a k 2 ) b e i k x α α + w t α α + θ c o t ( w + a k 2 ) 2 a x α α ν t α α + μ .
The mean value E [ s q 5 ] and the second moment E [ s q 5 2 ] are depicted in Figure 3a,b, respectively.

5.1.3. Hyperbolic Stochastic Solutions

s q 6 , 7 ( x , t ) = ± w + a k 2 b e i k x α α + w t α α + θ t a n h ( w + a k 2 ) 2 a x α α ν t α α + μ
The mean value E [ s q 7 ] and the second moment E [ s q 7 2 ] are depicted in Figure 4a,b, respectively.
And
s q 8 , 9 ( x , t ) = ± w + a k 2 b e i k x α α + w t α α + θ c o t h ( w + a k 2 ) 2 a x α α ν t α α + μ .
The mean value E [ s q 9 ] and the second moment E [ s q 9 2 ] are depicted in Figure 5a,b, respectively.

5.2. Stochastic Solutions of (21) via the Exp ( φ ( ξ ) ) -Expansion Method

Here, we can write the stochastic solutions of problem (21) using the exp ( φ ( ξ ) ) -expansion method and try to find the statistical properties of the stochastic processes solutions. As in the deterministic case, we can find the random relations as follows:
Case 1. At λ 2 4 ν > 0 , ν 0 ,
s q ˜ 1 , 2 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν e i k x α α + w t α α + θ tanh λ 2 4 ν 2 x α α ν t α α + K + λ ,
where α is a beta random variable.
The mean value E s q ˜ 1 and the second moment E s q ˜ 1 2 are depicted in Figure 6a,b, respectively.
Case 2. At λ 2 4 ν < 0 , ν 0 ,
s q ˜ 3 , 4 ( ξ ) = ± c 2 b 1 4 ν λ 2 4 ν e i k x α α + w t α α + θ coth λ 2 4 ν 2 x α α ν t α α + K + λ ,
where α is a beta random distribution.
The mean value E s q ˜ 3 and the second moment E s q ˜ 3 2 are depicted in Figure 7a,b, respectively.
Case 3. At λ 2 4 ν > 0 , ν = 0 , λ 0
s q ˜ 5 , 6 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 e i k x α α + w t α α + θ tan 4 ν λ 2 2 x α α ν t α α + K λ ,
where α is a beta random variable.
The mean value E s q ˜ 5 and the second moment E s q ˜ 5 2 are depicted in Figure 8a,b, respectively.
And
s q ˜ 7 , 8 ( ξ ) = ± c 2 b 1 + 4 ν 4 ν λ 2 e i k x α α + w t α α + θ tan 4 ν λ 2 2 x α α ν t α α + K λ .
where α is a beta random variable.
The mean value E s q ˜ 7 and the second moment E s q ˜ 7 2 are depicted in Figure 9a,b, respectively.
s q ˜ 9 , 10 ( ξ ) = ± c 2 b e i k x α α + w t α α + θ 1 + 2 λ exp λ x α α ν t α α + K 1 ,
where α is a beta random variable. The mean value E s q ˜ 9 and the second moment E s q ˜ 9 2 are depicted in Figure 10a,b, respectively.

5.3. The Influence of Randomness

Here, we show the effect of the beta random variable on the behavior of solutions. We present a number of graphs for various values of beta distribution parameters. As illustrated in Figure 1, Figure 5 and Figure 7, some are more dispersive than the others. Additionally, Figure 8 and Figure 9 show a near stable solution, whereas in Figure 4, Figure 6 and Figure 10, the surface becomes more planar after minor adjustments.

6. Conclusions

We have successfully applied the unified solver and the exp ( φ ( ξ ) ) -expansion approaches to extract some solitary waves through deterministic and beta distribution cases for the (stochastic) fractional Biswas–Milovic (FBM) equation in the sense of conformable fractional derivatives. Specifically, some new random solutions for the FBM equation have successfully been obtained. These solutions may have important significance for the explanation of some practical physical phenomena. The graphs of some solutions are illustrated for suitable coefficients. Additionally, the first and the second moments are obtained with graphical representation for the beta distribution case. This method can be applied for other stochastic NFDEs with complex valued solutions appearing in applied sciences.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The researcher would like to thank the Deanship of Scientific Research, Qassim University for funding the publication of this project.

Conflicts of Interest

The author declares that he has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  3. Abdelrahman, M.A.E.; AlKhidhr, H. Closed-form solutions to the conformable space-time fractional simplified MCH equation and time fractional Phi-4 equation. Results Phys. 2020, 18, 103294. [Google Scholar] [CrossRef]
  4. Abdelrahman, M.A.E.; Sohaly, M.A.; Alharbi, Y.F. Fundamental stochastic solutions for the conformable fractional NLSE with spatiotemporal dispersion via exponential distribution. Phys. Scr. 2021, 96, 125223. [Google Scholar] [CrossRef]
  5. Shah, N.A.; Ahmed, I.; Asogwa, K.K.; Zafar, A.A.; Weera, W.; Akgül, A. Numerical study of a nonlinear fractional chaotic Chua’s circuit. Aims Math. 2022, 8, 1636–1655. [Google Scholar] [CrossRef]
  6. Alshehry, A.S.; Shah, R.; Shah, N.A.; Dassios, I. A reliable technique for solving fractional partial differential equation. Axioms 2022, 11, 574. [Google Scholar] [CrossRef]
  7. Shang, N.; Zheng, B. Exact solutions for three fractional partial differential equations by the ( G G ) Method. Int. J. Appl. Math. 2013, 43, 1–6. [Google Scholar]
  8. Tong, B.; He, Y.; Wei, L.; Zhang, X. A generalized fractional sub-equation method for fractional differential equations with variable coefficients. Phys. Lett. 2012, 376, 2588–2590. [Google Scholar] [CrossRef]
  9. Lu, B. The first integral method for some time fractional differential equations. J. Math. Anal. Appl. 2012, 395, 684–693. [Google Scholar] [CrossRef] [Green Version]
  10. Ray, S.S.; Sahoo, S. A novel analytical method with fractional complex transform for new exact solutions of time-fractional fifth-order Sawada-Kotera equation. Rep. Math. Phys. 2015, 75, 63–72. [Google Scholar] [CrossRef]
  11. Abdelrahman, M.A.E.; Hassan, S.Z.; Alomair, R.A.; Alsaleh, D.M. The new wave structures to the fractional ion sound and Langmuir waves equation in plasma physics. Fractal Fract. 2022, 6, 227. [Google Scholar] [CrossRef]
  12. Bekir, A.; Guner, O.; Cevikel, A.C. The fractional complex transform and exp-function methods for fractional differential equations. Abstr. Appl. Anal. 2013, 2013, 426–462. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, W.; Chen, K. The functional variable method for finding exact solutions of some nonlinear time-fractional differential equations. Pramana J. Phys. 2013, 81, 377–384. [Google Scholar] [CrossRef]
  14. Abdelrahman, M.A.E.; Hassan, S.Z.; Alomair, R.A.; Alsaleh, D.M. Fundamental solutions for the conformable time fractional Phi-4 and space-time fractional simplified MCH equations. AIMS Math. 2021, 6, 6555–6568. [Google Scholar] [CrossRef]
  15. Wang, K.J.; Li, G.; Liu, J.H.; Wang, G.D. Solitary waves of the fractal regularized long-wave equation traveling along an unsmooth boundary. Fractals 2022, 30, 2250008. [Google Scholar] [CrossRef]
  16. Wang, K.J.; Wang, G.D.; Zhu, H.W. A new perspective on the study of the fractal coupled Boussinesq-Burger equation in shallow water. Fractals 2021, 29, 2150122. [Google Scholar] [CrossRef]
  17. Zhu, W.; Xia, Y.; Zhang, B.; Bai, Y. Exact traveling wave solutions and bifurcations of the time-fractional differential equations with applications. Int. J. Bifurc. Chaos 2019, 29, 1950041. [Google Scholar] [CrossRef]
  18. Li, Z.; Han, T.; Huang, C. Bifurcation and new exact traveling wave solutions for time-space fractional Phi-4 equation. Aip Adv. 2020, 10, 115113. [Google Scholar] [CrossRef]
  19. Liu, Z.; Zhang, K.; Li, M. Exact traveling wave solutions and bifurcation of a generalized (3+1)-dimensional time-fractional Camassa-Holm-Kadomtsev-Petviashvili equation. J. Funct. Spaces 2020, 2020, 4532824. [Google Scholar] [CrossRef]
  20. Li, Z.; Han, T. Bifurcation and exact solutions for the (2+1)-dimensional conformable time-fractional Zoomeron equation. Adv. Contin. Discret. Model. 2020, 2020, 656. [Google Scholar] [CrossRef]
  21. Biswas, A.; Milovic, D. Bright and dark solitons of the generalized nonlinear Schrödinger’s equation. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1473–1484. [Google Scholar] [CrossRef]
  22. Abdelrahman, M.A.E.; AlKhidhr, H. A robust and accurate solver for some nonlinear partial differential equations and tow applications. Phys. Scr. 2020, 95, 065212. [Google Scholar] [CrossRef]
  23. He, J.H.; Wu, X.H. Exp-function method for nonlinear wave equations. Chaos Solitons Fractals 2006, 30, 700–708. [Google Scholar] [CrossRef]
  24. Aminikhad, H.; Moosaei, H.; Hajipour, M. Exact solutions for nonlinear partial differential equations via Exp-function method. Numer. Methods Partial. Differ. Equ. 2009, 26, 1427–1433. [Google Scholar]
  25. Shakeel, M.; Attaullah; Alaoui, M.K.; Zidan, A.M.; Shah, N.A.; Weera, W. Closed-form solutions in a magneto-electro-elastic circular rod via generalized exp-function method. Mathematics 2022, 10, 3400. [Google Scholar] [CrossRef]
  26. Rani, A.; Shakeel, M.; Alaoui, M.K.; Zidan, A.M.; Shah, N.A.; Junsawang, P. Application of the exp(-φ(ξ))-expansion method to find the soliton solutions in biomembranes and nerves. Mathematics 2022, 10, 3372. [Google Scholar] [CrossRef]
  27. Evans, M.; Hastings, N.; Peacock, B. Statistical Distributions; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  28. Wang, K.J.; Wang, G.D.; Shi, F. Abundant exact traveling wave solutions to the local fractional (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equation. Fractals 2022. [Google Scholar] [CrossRef]
  29. Kimeu, J.M. Fractional Calculus: Definitions and Applications; Western Kentucky University: Bowling Green, KY, USA, 2009. [Google Scholar]
  30. Khalil, R.; Horani, M.A.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. [Google Scholar] [CrossRef]
  31. Abdeljawad, T. On conformable fractional calculus. J. Comput. Appl. Math. 2015, 279, 57–66. [Google Scholar] [CrossRef]
  32. Foroutan, M.; Kumar, D.; Manafian, J.; Hoque, A. New explicit soliton and other solutions for the conformable fractional Biswas–Milovic equation with Kerr and parabolic nonlinearity through an integration scheme. Optik 2018, 170, 190–202. [Google Scholar] [CrossRef]
  33. Ahmadian, S.; Darvishi, M.T. Fractional version of (1+1)-dimensional Biswas-Milovic equation and its solutions. Optik 2016, 127, 10135–10147. [Google Scholar] [CrossRef]
  34. Zaidan, L.I.; Darvishi, M.T. Semi-analytical solutions for different kinds of fractional Biswas-Milovic equation. Optik 2017, 136, 403–410. [Google Scholar] [CrossRef]
Figure 1. The first moment for the random solution s q 1 in (49); α is a beta random variable and 2 t , x 2 .
Figure 1. The first moment for the random solution s q 1 in (49); α is a beta random variable and 2 t , x 2 .
Fractalfract 06 00687 g001
Figure 2. The first and second moments for the random solution s q 3 in (50); α is a beta distribution random variable and 2 t , x 2 .
Figure 2. The first and second moments for the random solution s q 3 in (50); α is a beta distribution random variable and 2 t , x 2 .
Fractalfract 06 00687 g002
Figure 3. The first and second moments for the random solution s q 5 in (51); α is a beta random variable and 2 t , x 2 .
Figure 3. The first and second moments for the random solution s q 5 in (51); α is a beta random variable and 2 t , x 2 .
Fractalfract 06 00687 g003
Figure 4. The first and second moments for the random solution s q 7 in (52); α is a beta random variable and 2 t , x 2 .
Figure 4. The first and second moments for the random solution s q 7 in (52); α is a beta random variable and 2 t , x 2 .
Fractalfract 06 00687 g004
Figure 5. The first and second moments for the random solution s q 9 in (53); α is a beta random variable and 2 t , x 2 .
Figure 5. The first and second moments for the random solution s q 9 in (53); α is a beta random variable and 2 t , x 2 .
Fractalfract 06 00687 g005
Figure 6. The first and second moments for the random solution s q ˜ 1 in (54); α is a beta random variable and 2 t , x 2 , the first moment is on the upper part and the second moment is on the lower part.
Figure 6. The first and second moments for the random solution s q ˜ 1 in (54); α is a beta random variable and 2 t , x 2 , the first moment is on the upper part and the second moment is on the lower part.
Fractalfract 06 00687 g006
Figure 7. The first and second moments for the random solution s q ˜ 3 in (55); α is a beta distribution random variable and 2 t , x 2 .
Figure 7. The first and second moments for the random solution s q ˜ 3 in (55); α is a beta distribution random variable and 2 t , x 2 .
Fractalfract 06 00687 g007
Figure 8. The first and second moments for the random solution s q ˜ 5 in (56); α is a beta random variable and 2 t , x 2 .
Figure 8. The first and second moments for the random solution s q ˜ 5 in (56); α is a beta random variable and 2 t , x 2 .
Fractalfract 06 00687 g008
Figure 9. The first and second moments for the random solution s q ˜ 7 in (57); α is a beta distribution random variable and 2 t , x 2 .
Figure 9. The first and second moments for the random solution s q ˜ 7 in (57); α is a beta distribution random variable and 2 t , x 2 .
Fractalfract 06 00687 g009
Figure 10. The first and second moments for the random solution s q ˜ 9 in (58); α is a beta distribution random variable and 2 t , x 2 .
Figure 10. The first and second moments for the random solution s q ˜ 9 in (58); α is a beta distribution random variable and 2 t , x 2 .
Fractalfract 06 00687 g010aFractalfract 06 00687 g010b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almutairi, A. Fractional Biswas–Milovic Equation in Random Case Study. Fractal Fract. 2022, 6, 687. https://doi.org/10.3390/fractalfract6110687

AMA Style

Almutairi A. Fractional Biswas–Milovic Equation in Random Case Study. Fractal and Fractional. 2022; 6(11):687. https://doi.org/10.3390/fractalfract6110687

Chicago/Turabian Style

Almutairi, Abdulwahab. 2022. "Fractional Biswas–Milovic Equation in Random Case Study" Fractal and Fractional 6, no. 11: 687. https://doi.org/10.3390/fractalfract6110687

APA Style

Almutairi, A. (2022). Fractional Biswas–Milovic Equation in Random Case Study. Fractal and Fractional, 6(11), 687. https://doi.org/10.3390/fractalfract6110687

Article Metrics

Back to TopTop