Exact solutions of stochastic fractional Korteweg de–Vries equation with conformable derivatives
Ghany Hossam A.3, Hyder Abd-Allah1, 2, †, Zakarya M1, 4, ‡
King Khalid University, College of Science, Department of Mathematics, P. O. Box 9004, 61413, Abha, Saudi Arabia
Department of Engineering Mathematics and Physics, Faculty of Engineering, Al-Azhar University, 11371, Cairo, Egypt
Department of Mathematics, Helwan University, Sawah Street (11282), Cairo, Egypt
Department of Mathematics, Faculty of Science, Al-Azhar University, 71524, Assiut, Egypt

 

† Corresponding author. E-mail: h.abdelghany@yahoo.com mohammed_zakaria1983@yahoo.com

Abstract

We deal with the Wick-type stochastic fractional Korteweg de–Vries (KdV) equation with conformable derivatives. With the aid of the Exp-function method, white noise theory, and Hermite transform, we produce a novel set of exact soliton and periodic wave solutions to the fractional KdV equation with conformable derivatives. With the help of inverse Hermite transform, we get stochastic soliton and periodic wave solutions of the Wick-type stochastic fractional KdV equation with conformable derivatives. Eventually, by an application example, we show how the stochastic solutions can be given as Brownian motion functional solutions.

1. Introduction

The nonlinear fractional differential equations (FDEs) are constructed by mathematical modeling of some complex physical phenomena. The study of such nonlinear physical models through wave solutions analysis corresponding to their FDEs has a dynamic role in applied sciences. The Korteweg–de Vries (KdV) equations have been applied to a broad variety of material science phenomena as a model for the development and communication of nonlinear waves. The KdV equation was introduced to describe shallow water waves of long wavelength with small amplitude. The soliton and periodic exact solutions of the KdV equation may be of significance in many physical contexts as collision-free hydromagnetic waves, stratified internal waves, particle acoustic waves, plasma physics, and so on.[30,38,42]

This work is devoted to investigate the stochastic fractional KdV equation with conformable derivatives

where and 0 < α ≤ 1. Here, R and S are non-zero integrable functions from to the Kondrative distribution space which was defined by Holden et al. in Ref. [23] as a Banach algebra with the Wick-product “⋄”. Equation 1 is the perturbation of the variable coefficients fractional KdV equation with conformable derivatives

where r and s are non-zero integrable functions on . Equation (2) is a general model which describes shallow water waves of small amplitude and long wavelength.[22] Moreover, if equation (2) is considered in a random environment, we have a random fractional KdV equation. In order to obtain the exact solutions of the random fractional KdV equation, we only consider it in a white noise environment, that is, we will discuss the Wick-type stochastic fractional KdV equation (1).

There are many studies done for the definition and properties of the conformable derivative. Conformable forms of the chain rule, Gronwalls inequality, exponential functions, Taylor power series expansions, integration by parts, and Laplace transform have been presented by Abdeljawad in Ref. [1]. Benkhettoua et al.[3] have presented the calculus of the conformable time-scale. The heat equation with conformable derivatives was investigated by Hammad and Khalil in Ref. [18]. Chung[8] used the conformable derivative and integral to study the fractional Newtonian mechanics. Moreover, the deterministic conformable partial differential equations (PDEs) became an important subject in mathematical physics. So, many scholars paid more attention to their approximate and analytical solutions. The existence and uniqueness theorems for the linear sequential differential equations with conformable derivatives were proved by Gokdogan et al. in Ref. [17]. Eslami and Rezazadeh[9] gave a set of analytical solutions to the Wu–Zhang system with conformable derivative via the first integral method. The stochastic traveling wave solutions for the fractional coupled KdV and two-dimensional (2D) KdV equations were obtained by the modified fractional sub-equation method in Refs. [13,15], respectively.

Many researchers have studied the subject of random traveling wave, which is a significant subject of stochastic partial differential equations (SPDEs). Wadati[37] first proposed and discussed the stochastic KdV equation and discussed the propagation of soliton of the KdV equation under the effect of Gaussian noise. Furthermore, Ghany and Hyder,[12,13,15,16] Chen and Xie,[57] Hyder and Zakarya,[26,27] and Hyder[24,25] investigated a wide class of Wick-type stochastic evolution equations by using different extension methods and white noise analysis. Recently, many research works have done to investigate the conformable PDEs and their exact solutions via various methods. In Ref. [46], a conformable sub-equation method was proposed to construct exact solutions of the space–time resonant nonlinear Schröinger equation. By using the generalized exponential rational function method, new periodic and hyperbolic soliton solutions were constructed to the conformable Ginzburg–Landau equation with the Kerr law nonlinearity.[10] Also, a family of exact solutions were obtained for the space–time conformable generalized Hirota–Satsuma-coupled KdV equation and coupled mKdV equation using the Atangana’s conformable derivative and conformable sub-equation method.[43] The analysis of the first integral method was given in Refs. [44,45] to construct exact solutions of the nonlinear PDEs described by beta-derivative. Moreover, new optical, dark, complex, and singular soliton solutions were obtained for some nonlinear PDEs with M-derivative.[2,11] The investigation of exact and approximate solutions of nonlinear evolution equations plays an important role in the study of nonlinear physical phenomena. Many helpful methods such as bilinear transformation method,[35] the modified Clarkson and Kruskal (CK) direct method,[39] the multi-scale expansion method,[33] the binary Bell polynomials method,[48] the Riemann–Hilbert method,[41] and the approximate symmetry method[28] have been presented in the recent literature. He and Wu[19] introduced a concise and straightforward method, called the Exp-function method, to obtain generalized soliton, periodic, and compacton-like solutions for many nonlinear PDEs.[21, 49] In this method, the exact solutions are obtained in the form of an exponential type rational function in which both the numerator and denominator are polynomials of exponential functions. The main merit of the Exp-function method over the others lies in the fact that a great variety of exact solutions to nonlinear PDEs can be derived easily by choosing the parameters that appear, and the degrees of the polynomials that are present, in the solution. The main feature of the exact solutions obtained by the Exp-function method is their reducibility.[32] Also, the proposed Exp-function method leads to both the generalized solitonary solutions and periodic solutions. Moreover, the solution procedure, by the help of Mathematica, is of utter simplicity, and can be easily extended to all kinds of nonlinear equations.

Our aim in this work is to obtain new stochastic soliton and periodic wave solutions for the variable coefficients fractional KdV equation and Wick-type stochastic fractional KdV equation with conformable derivatives. Using white noise theory and Hermite transform, the Wick-type stochastic fractional KdV equation with conformable derivatives can be transformed to a deterministic fractional KdV equation containing conformable derivatives. In view of the proposed Exp-function method, the exact solutions of the deterministic fractional KdV equation are constructed in the form of an exponential type rational function in which both the numerator and denominator are polynomials of exponential functions. The highest and lowest degrees of the polynomial sums in both the numerator and denominator are determined through a homogeneous balance between the highest nonlinear term and the linear term of the highest order derivative appearing in the deterministic fractional KdV equation. Using some symbolic computation and the software Mathematica, we can find soliton and periodic wave solutions for the variable coefficients fractional KdV equation with conformable derivatives. Under pronounced conditions, we can apply the inverse Hermite transform to obtain stochastic soliton and periodic wave solutions for the Wick-type stochastic fractional KdV equation with conformable derivatives. Finally, by an application example, we show how the stochastic solutions can be given as Brownian motion functional solutions. This paper is organized as follows. In Section 2, we recall the definitions and some properties of the conformable derivative and integral, some requisites from Gaussian white noise analysis, and the main steps for solving the conformable nonlinear PDEs. In Section 3, we use the Exp-function method, white noise theory, and Hermite transform to obtain new stochastic soliton and periodic wave solutions for the Wick-type stochastic fractional KdV equation with conformable derivatives. In Section 4, we give an example to show that the stochastic solutions can be given as Brownian motion functional solutions. Section 5 is devoted to conclusion.

2. Preliminaries
2.1. The conformable derivative and integral

In this subsection, we recall the definitions and some properties of the conformable derivative and integral.

Definition 1[4,29]Let f be a function from (0,∞) into . For α ∈ (0,1], we define the conformable derivative of f of order α as follows:

Definition 2[4,29] Let f be an α-conformable differentiable function for t ∈ (0,a), a > 0 and exists. Then, and the conformable integral of the function f beginning from a ≥ 0 is given by

where the integral in the right hand side is the classical improper Riemann integral and α ∈ (0,1].

The following theorems give some sustainable properties for the conformable derivative.

Theorem 1[4,29] Assume that α ∈ (0,1], f and g are α-conformable differentiable functions at t ∈ (0,∞) and f is differentiable (in the usual sense) with respect to t. Then,

,

,

,

,

,

where ′ denotes the usual derivative with respect to t.

Theorem 2[31] Assume that the function f is a differentiable and α-conformable differentiable function on (0,∞). Also, assume that g is a differentiable function defined on the range of f. Then,

2.2. Basic concepts from Gaussian white noise analysis

The Gaussian white noise analysis starts with the rigging , where is the Schwartz space of rapidly decreasing, infinite differentiable functions on and is the space of tempered distributions. From the Bochner–Minlos theorem,[23] we have a unique white noise measure μ on . Assume that are the Hermite functions, where hn(x) denotes the Hermite polynomials. It is well known that the collection forms an orthonormal basis for . Let α = (α1,…,αd) be d-dimensional multi-indices with , then the family of tensor products ξα: = ξ(α1,…,αd) = ξα1 ⊗ ··· ⊗ ξαd, constitutes an orthonormal basis for . Now, introduce an ordering in by

where

Using this ordering, we define . Let be the set of all sequences with and with compact support. For , we define

Let , the Kondrative space of stochastic test functions is defined by

and the Kondrative space of stochastic distributions is defined by

The family of seminorms ||f||1,k, produces a topology on and can be represented as the dual of under the action ⟨F,f⟩ = ∑α(bα,cα)α!, where , , and (bα,cα) is the usual scalar product on .

The Wick product of two distributions , with is defined by

Let with . The Hermite transform of F is defined by

where and , with and .

For , by the definition of Hermite transform, we have

for all z such that and exist. The multiplication on the right hand side of the above equality is the complex bilinear multiplication in which is defined by , where . Hence, the Hermite transform converts the Wick product into the usual product and the convergence in into pointwise and bounded convergence in a specific neighborhood of zero in . For more details about stochastic Kondrative spaces, Wick product, and Hermite transform, we refer the reader to Ref. [23].

In what follows, the stochastic distribution process (or is a measurable function u from into . Moreover, if the function u is continuous, differentiable, C1, Ck, etc., then the process u has the same properties, respectively. Now, for q < ∞, r > 0, consider the infinite-dimensional neighborhoods of zero in .[23] To investigate the stochastic conformable PDEs, we need the following results.

Lemma 1[14,23] Suppose X(t,ϖ) and Y(t,ϖ) are such that (i) for each (t,z) ∈ (a,b) × Nq(r) and (ii) is a bounded function for (t,z) ∈ (a,b) × Nq(r) and continuous with respect to t ∈ (a,b) for each zNq(r). Then X(t,ω) has an α-order conformable derivative and for each t ∈ (a,b),

Lemma 2[14,23] Let X(t) be an . Suppose there exist q < ∞, r > 0 such that

and is a continuous function with respect to t ∈ [a,b] for each zNq(r). Then the α-order conformable integral operator of X(t) exists and

Now, consider a conformable nonlinear PDE in the form

where u = u(x,t) is the unknown function and P is a polynomial in u and its conformable derivatives. To obtain travelling wave solutions for Eq. (9), we use the following wave transformation:

where a > 0,k,c are constants and θ is a nonzero function to be determined later. Hence, equation (10) converts Eq. (9) to a nonlinear ordinary differential equation

Subsequently, the transformed equation (11) can be solved by using the Exp-function method as in the following section.

3. Exact solutions for Eq. (1)

Applying Hermite transform to Eq. (1), we get the conformable deterministic equation

where . To obtain traveling wave solutions to Eq. (14), we introduce the transformations , , with

where k,ω are arbitrary constants and θ is a nonzero function to be determined. Hence, equation (14) can be converted to the following NODE:

According to the Exp-function method, equation (16) can be solved by expanding its general solution in terms of the exponential function as follows:

where c, d, p, and q are freely chosen positive integers and an, bm are functions to be determined. Alternatively, we can write Eq. (17) as follows:

By balancing the highest order linear and nonlinear terms in Eq. (18), we can determine the numbers c and p. By simple calculation, we get p = c. Similarly, by balancing the lowest order linear and nonlinear terms in Eq. (18), we can determine the numbers d and q. Hence, we can obtain d = q.

Now, for some special cases of the numbers p,c,d, and q, we can express the general solution of Eq. (16) as follows.

Assume that k is an imaginary number, the above soliton wave solutions can be transformed to periodic solutions[19,20] according to the transformations

Then, the solution (21) becomes

where and . In order to obtain real periodic solutions, the imaginary part in Eq. (30) must be vanishing. For this, we set ν = 0, then we have b0 = ± 2b1. So, we get the following periodic solutions to Eq. (14):

Following the same technique as presented for the solution (21), we can convert the soliton wave solutions (25) and (28) to periodic solutions. Obviously, we can find different soliton and periodic wave solutions for Eq. (14) by setting different values for the numbers p,c,d, and q.

The properties of the exponential and trigonometric functions yield that there exists a bounded open set , q < ∞, r > 0 such that the solution u(x,t,z) of Eq. (14) and all its derivatives which are involved in Eq. (14) are uniformly bounded for (x,t,z) ∈ D × Nq(r), continuous with respect to (x,t) ∈ D for all zNq(r), and analytic with respect to zNq(r), for all (x,t) ∈ D. From Theorem 3, there exists such that for all (x,t,z) ∈ D × Nq(r) and U(x,t) solves Eq. (1) in . Hence, by applying the inverse Hermite transform to Eqs. (21), (25), (28), (31), and (32), we obtain the exact solutions of Eq. (1) as follows:

Stochastic soliton wave solutions

where .

Stochastic periodic wave solutions

4. Example and numerical simulation

In this section, we provide specific application example and numerical simulation to demonstrate the effectiveness of our results and to justify the real contribution of these results. We observe that the solutions of Eq. (1) are strongly depend on the shape of the given functions R(t) and S(t). So, for dissimilar forms of R(t) and S(t), we can find dissimilar solutions of Eq. (1) which come from Eqs. (33)–(37). We illustrate this by giving the following example.

Assume that α = 1, R(t) = ρS(t), and S(t) = f(t) + σ Wt, where ρ and σ are arbitrary constants, f(t) is a bounded measurable function on , and Wt is the Gaussian white noise which is the time derivative (in the strong sense in ) of the Brownian motion Bt. The Hermite transform of Wt is given by .[23] Using the definition of , equations (33)–(37) yield the white noise functional solution of Eq. (1) as follows:

where Θ = a1k + ρ b1(γ2+ k3),

The numerical simulation of the solution 38 is given in Figs. 1 and 2 for a = 0, a1 = 0.01, b0 = b1 = 1, k = 0.01, ρ = 1.5, σ = −2.5, and f(t) = sin2t. Figure 1 presents the evolutional behaviors of Eq. 38 with the noise effect Bt = Random[0,1] × tan(1.7t), and figure 2 presents the behavior of Eq. 38 without the effect of stochastic term Bt = 0. From Figs. 1 and 2, it is concluded that the stochastic forcing term leads to the uncertainty of the wave amplitude.

Fig. 1. (a) The 3D plot of the solution (38) and (b) contour plot of the solution (38) under the noise effect.
Fig. 2. (a) The 3D plot of the solution (38) and (b) contour plot of the solution (38) without the noise effect.
5. Conclusion

Due to the fact that the stochastic models are more realistic than the deterministic models, we concentrate our study in this paper on the Wick-type stochastic fractional KdV equation with conformable derivatives. Besides that, we investigate and solve the deterministic fractional KdV equation with conformable derivatives. By using the Exp-function method, Hermite transform, and white noise theory, we produce a new set of exact soliton and periodic wave solutions for the variable coefficients and fractional KdV equation with conformable derivatives. With the aid of inverse Hermite transform, we obtain stochastic soliton and periodic wave solutions for the Wick-type stochastic fractional KdV equation with conformable derivatives. Furthermore, we show by an example how the stochastic solutions can be given as Brownian motion functional solutions. Moreover, if α = 1, then the stochastic solutions (33)–(37) give a new set of stochastic solutions for the Wick-type stochastic KdV equation with integer derivatives.[40] Moreover, the set of solutions (21), (25), (28), (31), and (32) gives a new set of exact solutions for the variable coefficients and deterministic KdV equation.[47] Note that, the schema proposed in this paper can be used for solving several nonlinear evolution equations in mathematical physics, both Wick-type stochastic and deterministic.

Moreover, our exact solutions can be compared with other exact solutions which are obtained by different methods. For example, if we set α = b1 = 1, a1 = σ0, b0 = 2, k = 2δ, a = 0, R(t) = r(t), S(t) = s(t), and , where σ0, σ1, δ are constants and r(t), s(t) are deterministic integrable functions on . Then, the exact solution 33 becomes

Equation (43) can be rewritten in the form

Equation (44) is just the solution (32) in Ref. [36], which was obtained by an auxiliary equation method.

Reference
[1] Abdeljawad T 2015 J. Comput. Appl. Math. 279 57
[2] Baskonus H M Gómez-Aguilar J F 2019 Mod. Phys. Lett. 33 1950251
[3] Benkhettoua N Hassania S Torres D F M 2016 J. King Saud Univ. Sci. 28 93
[4] Çenesiz Y Baleanu D Kurt A Tasbozan O 2017 Waves in Random and Complex Media 27 103
[5] Chen B Xie Y 2005 Chaos Soliton. Fract. 23 281
[6] Chen B Xie Y C 2006 J. Comput. Appl. Math. 197 345
[7] Chen B Xie Y C 2007 J. Comput. Appl. Math. 203 249
[8] Chung W S 2015 J. Comput. Appl. Math. 290 150
[9] Eslami M Rezazadeh H 2016 Calcolo 53 475
[10] Ghanbaria B Gómez-Aguilarb J F 2019 Revista Mexicana de FĂ­sica 65 73
[11] Ghanbaria B Gómez-Aguilarb J F 2019 Mod. Phys. Lett. 33 1950235
[12] Ghany H A Hyder A 2012 International Review of Physics 6 153
[13] Ghany H A Okb El Babb A S Zabel A M Hyder A 2013 Chin. Phys. 22 080501
[14] Ghany H A Hyder A 2014 Int. J. Math. Analysis 23 2199
[15] Ghany H A Hyder A 2014 Chin. Phys. 23 060503
[16] Ghany H A Hyder A 2014 Kuwait Journal of Science 41 75
[17] Gökdoğan A Ünal E Çelik E 2016 Miskolc Math. Notes 17 267
[18] Hammad M A Khalil R 2014 Int. J. Pure Appl. Math. 94 215
[19] He J H Wu X H 2006 Chaos Soliton. Fract. 30 700
[20] He J H Abdou M A 2007 Chaos Soliton. Fract. 34 1421
[21] He J H 2008 Int. J. Mod. Phys. 22 3487
[22] Hereman W 2009 Shallow water waves and solitary waves, Encyclopedia of Complexity and Systems Science Meyers R. A. Heibelberg Springer Verlag 1620 1536
[23] Holden H Øsendal B Ubøe J Zhang T 2010 Stochastic partial differential equations New York Springer Science+Business Media, LLC
[24] Hyder A 2017 Journal of Mathematical Sciences: Advances and Applications 45 1
[25] Hyder A 2018 Pioneer Journal of Advances in Applied Mathematics 24 39
[26] Hyder A Zakarya M 2016 Int. J. Pure Appl. Math. 109 539
[27] Hyder A Zakarya M 27 5
[28] Jiao X Y 2018 Chin. Phys. 27 100202
[29] Khalil R Al Horani M Yousef A Sababheh M A 2014 J. Comput. Appl. Math. 246 65
[30] Khusnutdinova K R Stepanyants Y A Tranter M R 2018 Phys. Fluids 30 022104
[31] Kumarab D Seadawy R A Joardare A K 2018 Chin. J. Phys. 56 75
[32] Kudryashov N A Loguinova N B Commun. Nonlinear Sci. Numer. Simul. 14 1881
[33] Liu X Z Yu J Lou Z M et al. 2019 Chin. Phys. 28 010201
[34] Pérez J E S Gómez-Aguilar J F Baleanu D. Tchier F. 2018 Entropy 20 384
[35] Qian C Rao J G Liu Y B et al. 2016 Chin. Phys. Lett. 33 110201
[36] Taogetusang Sirendaoerji 2006 Chin. Phys. 15 1143
[37] Wadati M 1983 J. Phys. Soc. Jpn. 52 2642
[38] Wazwaz A M 2018 Mathematical Methods in the Applied Sciences 41 80
[39] Hu X R Chen J C Chen Y 2015 Chin. Phys. Lett. 32 070201
[40] Xie Y C 2003 Phys. Lett. 310 161
[41] Xu S Q Geng X G 2018 Chin. Phys. 27 120202
[42] Yuliawati L Budhi W S Adytia D 2019 J. Phys.: Conf. Ser. 1127 012065
[43] Yépez-Martínez H Gómez-Aguilar J F 2019 Waves in Random and Complex Media 29 678
[44] Yépez-Martínez H Gómez-Aguilar J F Atangana A 2018 Math. Model. Nat. Phenom. 13 14
[45] Yépez-Martínez H Gómez-Aguilar J F Baleanu D 2018 Optik 155 357
[46] Yépez-Martínez H Gómez-Aguilar J F 2019 Waves in Random and Complex Media https://dx.doi.org/10.1080/17455030.2019.1603413
[47] Zhang S 2007 Phys. Lett. 365 448
[48] Zheng P Jia M 2018 Chin. Phys. 27 120201
[49] Zhu S D 2007 Int. J. Nonlinear Sci. Numer. Simul. 8 465