Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Particle–Antiparticle Asymmetry in Relativistic Deformed Kinematics
Next Article in Special Issue
The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform
Previous Article in Journal
Particle-Based Dynamic Water Drops with High Surface Tension in Real Time
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation of Fractional-Order Swift–Hohenberg Equations via a Novel Transform

1
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Mathematics and Statistics, College of Science, Taif University, Taif 21944, Saudi Arabia
3
Department of Mathematics, College of Science, King Khalid University, Abha 61413, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Al-Azhar University, Assiut 71524, Egypt
5
Department of Mathematics, Abdul Wali University Mardan, Mardan 23200, Pakistan
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(7), 1263; https://doi.org/10.3390/sym13071263
Submission received: 21 June 2021 / Revised: 6 July 2021 / Accepted: 8 July 2021 / Published: 14 July 2021

Abstract

:
In this paper, the Elzaki transform decomposition method is implemented to solve the time-fractional Swift–Hohenberg equations. The presented model is related to the temperature and thermal convection of fluid dynamics, which can also be used to explain the formation process in liquid surfaces bounded along a horizontally well-conducting boundary. In the Caputo manner, the fractional derivative is described. The suggested method is easy to implement and needs a small number of calculations. The validity of the presented method is confirmed from the numerical examples. Illustrative figures are used to derive and verify the supporting analytical schemes for fractional-order of the proposed problems. It has been confirmed that the proposed method can be easily extended for the solution of other linear and non-linear fractional-order partial differential equations.

1. Introduction

The concept of Fractional Calculus (FC) is old, which arises from the nth derivative notation used by Leibniz in their publication. So, L’Hopital asks Leibniz what the result would be if the order is non-integer [1]. Riemann and Liouville defined the concept of fractional order differentiation in the 19th century. Later on, researchers began to research FC and found that fractional-order models are more suitable than integer-order models for some real-world problems [2,3,4]. FC is an efficient and powerful tool for describing memory and hereditary properties in various materials and processes.
Fractional Differential Equations (FDEs) have significantly gained much attention from researchers due to providing fractional modeling of different phenomena in nature. Due to this reason, the implementation of FDEs to model different physical systems and processes has been increased, for example, colored noise [5], economics [6], oscillation of earthquake [7], and bioengineering [8]. The other applications are control theory [9], rheology [10], visco-elastic materials [11], signal processing [12], damping method [13], polymers [14], and so on. The scheme consisting of integer partial differential equations and fractional-order partial differential equations with the fractional Caputo derivative has a well-designed symmetry structure. This problem is utilized to analyze dispersive wave phenomena in different areas of applied science, like quantum mechanics and plasma physics. Nonlinear phenomena play a crucial role in applied mathematics and physics; we know that most engineering problems are non-linear, and solving them analytically is difficult. In physics and mathematics, obtaining exact or approximate solutions to nonlinear FPDEs is still a significant problem that requires new methods to discover exact or approximate solutions.
Because of the above fact, researchers have developed numerous numerical and analytical techniques for the solution of FPDEs [15,16]. In [17], A.A. Alderremy et al. used Modified Reduced Differential Transform Method (MRDTM) to solve the fractional nonlinear Newell–Whitehead–Segel equation. M.S. Rawashdeh and H. Al-Jammal [18] implemented the fractional natural decomposition method (FNDM) for finding approximate analytical solutions to systems of nonlinear PDEs. Certain analytical solutions of the fractional-order diffusion equations were found by K. Shah et al. [19], who used Natural Transform Method (NTM). To obtain the approximate and exact solutions of space and time-fractional Burgers equations with initial conditions [20], M. Inc implemented a variational iteration method.
In [21], using the approximate analytical method, travelling wave solutions for Korteweg-de Vries equations having fractional-order were discussed. Similarly, in [22], F.A. Alawad et al. solved space-time fractional telegraph equations using a new technique of the Laplace variational iteration method. H. Jafari et al. [23] found the approximate solution of the nonlinear gas dynamic equation by implementing homotopy analysis method. To obtain a series form solution of time-fractional coupled Burgers equations. P. Veereshaa and D.G. Prakash used a reliable technique q-homotopy analysis transform method (q-HATM) [24]. In [25], L. Yan used the iterative Laplace transform method, which combines two methods, the iterative method, and the Laplace transform method, to obtain the numerical solutions of fractional Fokker–Planck equations.
However, we used a new technique formed by the combination of Elzaki transform [26] and the Adomian decomposition method [27,28] known as the Elzaki Transform Decomposition Method (ETDM). The Elzaki transformation is renowned for handling linear ordinary differential equations, linear partial differential equations, and integral equations, as seen in [29,30,31]. In contrast, the Adomian decomposition method [27,28] is a well-known method for handling linear and nonlinear, homogeneous and nonhomogeneous differential and partial differential equations, integro-differential, and FDEs series form solution.
In this paper, we aim to solve Swift–Hohenberg (S-H) equation with the help of ETDM. The S-H equation was first introduced and derived from the equations for thermal convection by J. Swift and P. Hohenberg [32]. The general form of the S-H equation is
δ ρ τ δ = b ρ 1 + 2 τ 2 ρ N ( ρ ) , b R , 0 < δ 1 ,
where ρ is a scalar function, b is the real constant, and N ( ρ ) is a nonlinear term. The S-H equation has many applications in engineering and science, such as physics, biology, laser study fluid, and hydro-dynamics [33,34,35]. The S-H equation plays an important role in pattern formation theory in fluid layers confined between horizontal well-conducting boundaries [36]. This equation has many applications in the modeling pattern formation and its different issues, including the selection of pattern, effects of noise on bifurcations, the dynamics of defects, and spatiotemporal chaos [37,38,39,40].

2. Preliminaries

In this subsection, we recall some simple and most significant concepts concerning fractional calculus.
Definition 1
([41,42,43]). Abel-Riemann (A-R) described D δ operator of the δ order as
D δ ρ ( ψ ) = d n d ψ n ρ ( ψ ) , δ = n ; 1 Γ ( n δ ) d d ψ n 0 ψ ρ ( T ) ( ψ τ ) δ n + 1 d τ , n 1 < δ < n ,
where δ R + , n Z + and
D δ ρ ( ψ ) = 1 Γ ( δ ) 0 ψ ( ψ τ ) δ 1 ρ ( τ ) d τ , 0 < δ 1 .
Definition 2
([42,43]). The fractional order A-R integral operator J δ is given as
J δ ρ ( ψ ) = 1 Γ ( δ ) 0 ψ ( ψ τ ) δ 1 ρ ( τ ) d τ , τ > 0 , δ > 0 .
By Podlubny [42] we may have
J δ τ n = Γ ( n + 1 ) Γ ( n + δ + 1 ) τ n + δ ,
D δ τ n = Γ ( n + 1 ) Γ ( n δ + 1 ) τ n δ .
Definition 3
([41,42,44,45]). In the Caputo manner, the operator D δ with the order δ is given as
D δ ρ ( ψ ) = 1 Γ ( n δ ) 0 ψ ρ n ( τ ) ( ψ τ ) δ n + 1 d τ , n 1 < δ < n ; d n d τ n ρ ( ψ ) , δ = n
with the following properties:
( a ) D τ δ J τ δ f ( τ ) = f ( τ ) ; ( b ) J τ δ D τ δ f ( τ ) = f ( τ ) k = 0 n f k ( 0 + ) τ k k !
for τ > 0 and n 1 < δ n , n N .
Definition 4
([46]). The Mittag–Leffler function ψ is defined as
E δ ( ψ ) = k = 0 ψ k Γ ( δ k + 1 ) , δ > 0 .
For a f ( t ) function, the ET or modified Sumudu transform definition is given as
E [ f ( τ ) ] = F ( r ) = r 0 h ( τ ) e τ r d τ , r > 0 .
The transformation of Elzaki is a very useful and powerful tool for solving the integral equation that can not be solved by the Sumudu transformation method.
The following ET transformations of partial derivatives, which can be obtained by using integration by parts, may be used in (8):
  • E [ f ( ψ , τ ) τ ] = 1 r F ( ψ , r ) r f ( ψ , 0 ) .
  • E [ 2 f ( ψ , τ ) τ 2 ] = 1 r 2 F ( ψ , r ) f ( ψ , 0 ) r f ( ψ , 0 ) τ .
  • E [ f ( ψ , τ ) ψ ] = d d ψ F ( ψ , r ) .
  • E [ 2 f ( ψ , τ ) ψ 2 ] = d 2 d ψ 2 F ( ψ , r ) .
Theorem 1
([47]). Let U ( s ) be the Laplace transform of f ( τ ) then ET F ( r ) of f ( τ ) is specified as
F ( r ) = r U ( 1 r ) .
Theorem 2
([47]). If F ( r ) is the ET of the f ( τ ) function, then
E [ D δ f ( τ ) ] = F ( r ) r δ k = 0 n 1 r k δ + 2 f ( k ) ( 0 ) , n 1 < δ n .

3. Idea of ETDM

The ETDM solution for fractional partial differential equations is described in this section.
D τ δ ρ ( ψ , τ ) + G ¯ 1 ρ ( ψ , τ ) + N 1 ρ ( ψ , τ ) P 1 ( ψ , τ ) = 0 , n 1 < δ n
with initial conditions
ρ ( ψ , 0 ) = g 1 ( ψ ) ,
where D τ δ = δ τ δ is the fractional derivative in Caputo sense having order δ , G ¯ 1 and N 1 are linear and non-linear functions, respectively, and source operator is P 1 .
By applying Elzaki transform on both sides of (13), we obtain
E [ D τ δ ρ ( ψ , τ ) ] + E [ G ¯ 1 ρ ( ψ , τ ) + N 1 ρ ( ψ , τ ) P 1 ( ψ , τ ) ] = 0 .
By Elzaki transform property of differentiation, we get
E [ ρ ( ψ , τ ) ] = s δ k = 0 m 1 s 2 + k δ k ρ ( ψ , τ ) k τ | τ = 0 + s δ E [ P 1 ( ψ , τ ) ] s δ E { G ¯ 1 ρ ( ψ , τ ) + N 1 ρ ( ψ , τ ) } ] .
ETDM determines the solution of the infinite sequence of ρ ( ψ , τ )
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) .
The decomposition of nonlinear terms by Adomian polynomials N 1 is defined as
N 1 ρ ( ψ , τ ) = m = 0 A m .
The Adomian polynomials can represent all forms of nonlinearity as
A m = 1 m ! m m N 1 k = 0 k ψ k , k = 0 k τ k = 0 .
Putting (17) and (19) into (16), gives
E m = 0 ρ m ( ψ , τ ) = s δ k = 0 m 1 s 2 + k δ k ρ ( ψ , τ ) k τ | τ = 0 + s δ E { P 1 ( ψ , τ ) } s δ E G ¯ 1 ρ ( m = 0 ψ m , m = 0 τ m ) + m = 0 A m .
By applying inverse Elzaki to (20), we obtain
m = 0 ρ m ( ψ , τ ) = E s δ k = 0 m 1 s 2 + k δ k ρ ( ψ , τ ) k τ | τ = 0 + s δ E { P 1 ( ψ , τ ) } s δ E G ¯ 1 ρ m = 0 ψ m , m = 0 τ m + m = 0 A m .
The following terms are described as
ρ 0 ( ψ , τ ) = E s δ k = 0 m 1 s 2 + k δ k ρ ( ψ , τ ) k τ | τ = 0 + s δ E + { P 1 ( ψ , τ ) } , ρ 1 ( ψ , τ ) = E s δ E + { G ¯ 1 ρ ( ψ 0 , τ 0 ) + A 0 } ,
for m 1 , is determined as
ρ m + 1 ( ψ , τ ) = E s δ E + { G ¯ 1 ρ ( ψ m , τ m ) + A m } .

4. Existence and Uniqueness Results for ETDM

In what follows, we will demonstrate that the sufficient conditions assure the existence of a unique solution. Our desired existence of solutions in the case of SDM follows by [40].
Theorem 3.
(Uniqueness theorem): Equation (23) has a unique solution whenever 0 < ϵ < 1 , where ϵ = ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) τ ( δ 1 ) δ ! .
Proof. 
Assume that M = C [ I ] , . represents all continuous mappings on the Banach space, defined on I = [ 0 , T ] having the norm . . For this we introduce a mapping W : M M , we have
ρ n + 1 ( ψ , τ ) = ρ ( ψ , τ ) + E 1 [ ω ξ δ E [ L ρ n ( ψ , τ ) + R [ ρ n ( ψ , τ ) ] + N [ ρ n ( ψ , τ ) ] ] ] , n 0 ,
where L ρ ( ψ , τ ) 3 ρ ( ψ , τ ) ψ 2 and R ρ ( ψ , τ ) ρ ( ψ , τ ) ψ . Now assume that L ρ ( ψ , τ ) and M ρ ( ψ , τ ) are also Lipschitzian with | R ρ R ρ ˇ | < L ˇ 1 | ρ ρ ˇ | and | L ρ L ρ ˇ | < L ˇ 2 | ρ ρ ˇ | , where L ˇ 1 and L ˇ 2 are Lipschitz constant, respectively, and ρ , ρ ˇ are various values of the mapping.
W ρ W ρ ˇ = max τ I E 1 [ s δ E [ L ρ ( ψ , τ ) + R ρ ( ψ , τ ) + N ρ ( ψ , τ ) ] ] E 1 [ s δ E [ L ρ ˇ ( ψ , τ ) + R ρ ˇ ( ψ , τ ) + N ρ ˇ ( ψ , τ ) ] ] max τ I E 1 [ s δ E [ L ρ ( ψ , τ ) L ρ ˇ ( ψ , τ ) ] ] + E 1 [ s δ E [ R ρ ( ψ , τ ) R ρ ˇ ( ψ , τ ) ] ] + E 1 [ s δ E [ N Φ ( ψ , τ ) N ρ ˇ ( ψ , τ ) ] ] max τ I L ˇ 1 E 1 [ s δ E | ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) | ] + L ˇ 2 E 1 [ s δ E | ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) | ] + L ˇ 3 E 1 [ s δ E | ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) | ]
max τ I ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) E 1 [ s δ E | ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) | ] ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) E 1 [ s δ E ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) ] = ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) τ ( δ 1 ) δ ! ρ ( ψ , τ ) ρ ˇ ( ψ , τ ) .
Under the assumption 0 < ϵ < 1 , the mapping is contraction. Thus, by Banach contraction fixed point theorem, there exists a unique solution to (13). Therefore, this completes the proof. □
Theorem 4
(Convergence Analysis). The general form solution of (13) will be convergent.
Proof. 
Suppose S ^ n be the n t h partial sum, that is S ^ n = m = 0 n ρ m ( ψ , τ ) . Firstly, we show that { S ^ n } is a Cauchy sequence in Banach space in M . Taking into consideration a new representation of Adomian polynomials we obtain
R ¯ ( S ^ n ) = H ˇ n + p = 0 n 1 H ˇ p , N ¯ ( S ^ n ) = H ˇ n + c = 0 n 1 H ˇ c .
Now
S ^ n S ^ q = max τ I | S ^ n S ^ q | = max τ I | m = q + 1 n ρ ˇ ( ψ , τ ) | , ( m = 1 , 2 , 3 , ) max τ I E 1 s δ E m = q + 1 n L ρ n 1 ( ψ , τ ) + E 1 s δ E m = q + 1 n R ρ n 1 ( ψ , τ ) + E 1 s δ E m = q + 1 n H ˇ n 1 ( ψ , τ ) = max τ I E 1 s δ E m = q n 1 L ρ n ( ψ , τ ) + E 1 s δ E m = q n 1 R ρ n ( ψ , τ ) + E 1 s δ E m = q n 1 H ˇ n ( ψ , τ ) max τ I E 1 s δ E m = q n 1 L ( S ^ n 1 ) L ( S ^ q 1 ) + E 1 s δ E m = q n 1 R ( S ^ n 1 ) R ( S ^ q 1 ) + E 1 s δ E m = q n 1 N ( S ^ n 1 ) N ( S ^ q 1 )
max τ I E 1 s δ E L ( S ^ n 1 ) L ( S ^ q 1 ) + E 1 s δ E R ( S ^ n 1 ) R ( S ^ q 1 ) + E 1 s δ E N ( S ^ n 1 ) N ( S ^ q 1 ) L ˇ 1 max τ I E 1 s δ E ( S ^ n 1 ) ( S ^ q 1 ) + L ˇ 2 max τ I E 1 s δ E ( S ^ n 1 ) ( S ^ q 1 ) + L ˇ 3 max τ I E 1 s δ E ( S ^ n 1 ) ( S ^ q 1 ) = ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) τ ( δ 1 ) δ ! S ^ n 1 S ^ q 1 .
Consider n = q + 1 ; then
S ^ q + 1 S ^ q ϵ S ^ q S ^ q 1 ϵ 2 S ^ q 1 S ^ q 2 ϵ q S ^ 1 S ^ 0 ,
where ( L ˇ 1 + L ˇ 2 + L ˇ 3 ) τ ( δ 1 ) δ ! . Analogously, from the triangular inequality we have
S ^ n S ^ q S ^ q + 1 S ^ q + S ^ q + 2 S ^ q + 1 + + S ^ n S ^ n 1 ϵ q + ϵ q + 1 + + ϵ n 1 S ^ 1 S ^ 0 ϵ q 1 ϵ n q ϵ ρ 1 ,
since 0 < ϵ < 1 , we have ( 1 ϵ n q ) < 1 , then
S ^ n S ^ q ϵ q 1 ϵ max τ I ρ 1 .
However, | ρ 1 | < (since ρ ( ψ , τ ) is bounded). Thus, as q , then S ^ n S ^ q 0 . Hence, { S ^ 1 } is a Cauchy sequence in K . As a result, the series n = 0 ρ n is convergent and this completes the proof. □
Theorem 5
([40] (Error estimate)). The maximum absolute truncation error of the series solution (13) to (23) is computed as
max τ I | ρ ( ψ , τ ) n = 1 q ρ n ( ψ , τ ) | ϵ q 1 ϵ max τ I ρ 1 .

5. Numerical Examples:

Example 1.
Consider the following linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 + 3 ρ ( ψ , τ ) ψ 3 = 0 , 0 < δ 1 ,
with initial condition
ρ ( ψ , 0 ) = exp ( ψ ) .
Taking Elzaki transformation of (27), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 3 ρ ( ψ , τ ) ψ 3 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 3 ρ ( ψ , τ ) ψ 3 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 3 ρ ( ψ , τ ) ψ 3 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 3 ρ ( ψ , τ ) ψ 3 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution:
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) , m = 0 ρ m ( ψ , τ ) = exp ψ + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 3 ρ ( ψ , τ ) ψ 3 .
Thus, by comparing both sides of (30), we have
ρ 0 ( ψ , τ ) = exp ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( b 4 ) exp ψ 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( b 4 ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ ;
for m = 2 ,
ρ 3 ( ψ , τ ) = ( b 4 ) 3 exp ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ ;
for m = 3 ,
ρ 4 ( ψ , τ ) = ( b 4 ) 4 exp ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 3 ) are easy to get. Thus, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + ρ 3 ( ψ , τ ) + ρ 4 ( ψ , τ ) + , ρ ( ψ , τ ) = exp ( ψ ) + ( b 4 ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 4 ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ + ( b 4 ) 3 exp ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ + ( b 4 ) 4 exp ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ + .
Exact solution for (27) at δ = 1 is
ρ ( ψ , τ ) = exp ( ψ ) E δ ( ( b 4 ) τ δ ) .
The Figure 1 shows the ETDM graph for Example 1 at various fractional order.
Example 2.
Consider the following linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 + 4 ρ ( ψ , τ ) ψ 4 = 0 , 0 < δ 1
with initial condition
ρ ( ψ , 0 ) = sin ( ψ ) .
Taking Elzaki transformation of (35), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) , m = 0 ρ m ( ψ , τ ) = sin ( ψ ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Here we will discuss the following two cases.
Case 1: ( b = 0 ). According to (35), we have
ρ 0 ( ψ , τ ) = sin ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = 0 ;
for m = 1 ,
ρ 2 ( ψ , τ ) = 0 ;
for m = n ,
ρ n ( ψ , τ ) = 0 .
Hence the solution to (35) in this case is
ρ ( ψ , τ ) = sin ( ψ ) .
Case 2: ( b 0 ) . According to (37), we have
ρ 0 ( ψ , τ ) = sin ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( ( b 1 ) + 1 ) sin ( ψ ) 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( ( b 1 ) 2 + ( b 1 ) + 1 ) sin ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ ;
for m = 2 ,
ρ 3 ( ψ , τ ) = ( ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 ) sin ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ ;
for m = 3 ,
ρ 4 ( ψ , τ ) = ( ( b 1 ) 4 + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 ) sin ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 3 ) are easy to get. Thus, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + ρ 3 ( ψ , τ ) + ρ 4 ( ψ , τ ) + , ρ ( ψ , τ ) = sin ( ψ ) + ( ( b 1 ) + 1 ) sin ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 1 ) 2 + ( b 1 ) + 1 sin ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 sin ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ + ( b 1 ) 4 + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 sin ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ + .
Exact solution for (33) at δ = 1 is
ρ ( ψ , τ ) = 1 b sin ( ψ ) lim N n = 0 N 1 ( b 1 ) n + 1 1 Γ ( n δ + 1 ) τ n δ .
The Figure 2 shows ETDM solution graph at different fractional order for Example 2.
Example 3.
Consider the following linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 + 4 ρ ( ψ , τ ) ψ 4 = 0 , 0 < δ 1
with initial condition
ρ ( ψ , 0 ) = cos ( ψ ) .
Taking Elzaki transformation of (41), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) , m = 0 ρ m ( ψ , τ ) = cos ( ψ ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 .
Here we will discuss the following two cases.
Case 1: ( b = 0 ). According to Equation (44), we have
ρ 0 ( ψ , τ ) = cos ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = 0 ;
for m = 1 ,
ρ 2 ( ψ , τ ) = 0 ;
for m = n ,
ρ n ( ψ , τ ) = 0 .
Hence the solution to (39) in this case is
ρ ( ψ , τ ) = cos ( ψ ) .
Case 1: ( b 0 ) . According to (43), we have
ρ 0 ( ψ , τ ) = cos ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( ( b 1 ) + 1 ) cos ( ψ ) 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( ( b 1 ) 2 + ( b 1 ) + 1 ) cos ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ ;
for m = 2 ,
ρ 3 ( ψ , τ ) = ( ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 ) cos ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ ;
for m = 3 ,
ρ 4 ( ψ , τ ) = ( ( b 1 ) 4 + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 ) cos ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 3 ) are easy to get. Thus, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + ρ 3 ( ψ , τ ) + ρ 4 ( ψ , τ ) + , ρ ( ψ , τ ) = cos ( ψ ) + ( b 1 ) + 1 cos ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 1 ) 2 + ( b 1 ) + 1 cos ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 cos ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ + ( b 1 ) 4 + ( b 1 ) 3 + ( b 1 ) 2 + ( b 1 ) + 1 cos ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ + .
Exact solution for (41) at δ = 1 is
ρ ( ψ , τ ) = 1 b cos ( ψ ) lim N n = 0 N 1 ( b 1 ) n + 1 1 Γ ( n δ + 1 ) τ n δ .
The Figure 3 shows ETDM solution graph at different fractional order for Example 3.
Example 4.
Consider the following linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 σ 3 ρ ( ψ , τ ) ψ 3 + 4 ρ ( ψ , τ ) ψ 4 = 0 , 0 < δ 1
with initial condition
ρ ( ψ , 0 ) = exp ( ψ ) .
Taking Elzaki transformation of (47), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) , m = 0 ρ m ( ψ , τ ) = exp ( ψ ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 .
Thus, by comparing both sides of (49), we have
ρ 0 ( ψ , τ ) = exp ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( b 4 + σ ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( b 4 + σ ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ ;
for m = 2 ,
ρ 3 ( ψ , τ ) = ( b 4 + σ ) 3 exp ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ ;
for m = 3 ,
ρ 4 ( ψ , τ ) = ( b 4 + σ ) 4 exp ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 3 ) are easy to get. Thus, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + ρ 3 ( ψ , τ ) + ρ 4 ( ψ , τ ) + , ρ ( ψ , τ ) = exp ( ψ ) + ( b 4 + σ ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 4 + σ ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ + ( b 4 + σ ) 3 exp ( ψ ) 1 Γ ( 3 δ + 1 ) τ 3 δ + ( b 4 + σ ) 4 exp ( ψ ) 1 Γ ( 4 δ + 1 ) τ 4 δ + .
Exact solution for (47) at δ = 1 is
ρ ( ψ , τ ) = exp ( ψ ) E δ ( ( b 4 + σ ) τ δ ) .
The Figure 4 shows ETDM solution graph at different fractional order for Example 4.
Example 5.
Consider the following non-linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 + 4 ρ ( ψ , τ ) ψ 4 ρ 2 ( ψ , τ ) + ρ ( ψ , τ ) ψ 2 = 0 , 0 < δ 1 ,
with initial condition
ρ ( ψ , 0 ) = exp ( ψ ) .
Taking Elzaki transformation of (53), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution:
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) ,
where the Adomian polynomials ρ 2 = m = 0 A m and ( ρ ψ ) 2 = m = 0 B m and the nonlinear terms have been characterised. (53) can be rewritten in the form using certain terms
m = 0 ρ m ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + m = 0 A m m = 0 B m , m = 0 ρ m ( ψ , τ ) = exp ( ψ ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 4 ρ ( ψ , τ ) ψ 4 + m = 0 A m m = 0 B m .
The decomposition of nonlinear terms by Adomian polynomials is defined as, according to (20),
A 0 = ρ 0 2 , A 1 = 2 ρ 0 ρ 1 , B 0 = ( ρ 0 ψ ) 2 , B 1 = 2 ρ 0 ψ ρ 1 ψ .
Thus, by comparing both sides of (56), we have
ρ 0 ( ψ , τ ) = exp ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( b 4 ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( b 4 ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 1 ) are easy to obtain. So, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + ρ ( ψ , τ ) = exp ( ψ ) + ( b 4 ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 4 ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ +
Exact solution for (51) at δ = 1 is
ρ ( ψ , τ ) = exp ( ψ ) E δ ( ( b 4 ) τ δ ) .
The Figure 5 shows the ETDM graph for Example 5 at various fractional order.
Example 6.
Consider the following non-linear time-fractional S-H equation
δ ρ ( ψ , τ ) τ δ + ( 1 b ) ρ ( ψ , τ ) + 2 2 ρ ( ψ , τ ) ψ 2 σ 3 ρ ( ψ , τ ) ψ 3 + 4 ρ ( ψ , τ ) ψ 4 ρ 2 ( ψ , τ ) + ρ ( ψ , τ ) ψ 2 = 0 , 0 < δ 1
with initial condition
ρ ( ψ , 0 ) = exp ( ψ ) .
Taking Elzaki transformation of (59), we obtain
E δ ρ τ δ = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 , 1 s δ E ρ ( ψ , τ ) s 2 δ ρ ( ψ , 0 ) = E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
The above algorithm’s simplified form is
E ρ ( ψ , τ ) = s 2 ρ ( ψ , 0 ) + s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
Using inverse Elzaki transformation, we get
ρ ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 + ρ 2 ( ψ , τ ) ρ ( ψ , τ ) ψ 2 .
Assume that the unknown ρ ( ψ , τ ) function, in infinite series form, has the following solution
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) ,
where the Adomian polynomials ρ 2 = m = 0 A m and ( ρ ψ ) 2 = m = 0 B m and the nonlinear terms have been characterised. (61) can be rewritten in the form using certain terms
m = 0 ρ m ( ψ , τ ) = ρ ( ψ , 0 ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 + m = 0 A m m = 0 B m , m = 0 ρ m ( ψ , τ ) = exp ( ψ ) + E s δ E ( 1 b ) ρ ( ψ , τ ) 2 2 ρ ( ψ , τ ) ψ 2 + σ 3 ρ ( ψ , τ ) ψ 3 4 ρ ( ψ , τ ) ψ 4 + m = 0 A m m = 0 B m .
The decomposition of nonlinear terms by Adomian polynomials is defined as, according to (20),
A 0 = ρ 0 2 , A 1 = 2 ρ 0 ρ 1 , B 0 = ( ρ 0 ψ ) 2 , B 1 = 2 ρ 0 ψ ρ 1 ψ .
Thus, by comparing both sides of (63), we have
ρ 0 ( ψ , τ ) = exp ( ψ ) ;
for m = 0 ,
ρ 1 ( ψ , τ ) = ( b 4 + σ ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ ;
for m = 1 ,
ρ 2 ( ψ , τ ) = ( b 4 + σ ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ .
Similarly, the remaining ETDM solution elements ρ m ( m 1 ) are easy to get. Thus, we define the sequence of alternatives as
ρ ( ψ , τ ) = m = 0 ρ m ( ψ , τ ) = ρ 0 ( ψ , τ ) + ρ 1 ( ψ , τ ) + ρ 2 ( ψ , τ ) + , ρ ( ψ , τ ) = exp ( ψ ) + ( b 4 + σ ) exp ( ψ ) 1 Γ ( δ + 1 ) τ δ + ( b 4 + σ ) 2 exp ( ψ ) 1 Γ ( 2 δ + 1 ) τ 2 δ + .
Exact solution for (58) at δ = 1 is
ρ ( ψ , τ ) = exp ( ψ ) E δ ( ( b 4 + σ ) τ δ ) .
The Figure 6 shows the ETDM graph for Example 6 at various fractional order.

6. Results and Discussion

In this paper, ETDM is implemented to solve time-fractional Swift–Hohenberg equations. The results, we get by using suggested technique are explain with the help of its graphical representation. Figure 1 show the 3D and 2D graph at different values of δ . The ETDM solution graph are plotted at b = 5 in the domain 4 ψ 4 . In Figure 2, ETDM solutions graphs at ( a ) δ = 0.25 , ( b ) δ = 0.50 and ( a ) δ = 0.75 are plotted in which we fix b = 1 in the given domain 6 ψ 6 . In Figure 3, the ETDM solutions at fractional orders are drawn. The graph (a) represent the solution of Example 3 at ( a ) δ = 0.25 , ( b ) δ = 0.50 , while graph (c) is the plotted at δ = 0.75 . The given figures are plotted at b = 1 with ψ ranges from 0 to 1. In Figure 4, the ETDM solutions are plotted at various fractional order for b = 5 , σ = 1 with 0 ψ 5 . The graph (a) represent the solution of Example 4 at ( a ) δ = 0.25 , ( b ) δ = 0.50 , while graph (c) is the plotted at δ = 0.75 . The solution in Figure 5 are calculated at different fractional-orders. It is observed that the solutions at various fractional-orders are converges to he solution of integer-order solution as fractional-orders approaches to an integer-order. The graphs are plotted at b = 5 having 0 ψ 1 . In Figure 6, the same graphical representation have been made at b = 5 , σ = 1 and 1 ψ 1 .

7. Conclusions

An efficient analytical technique is used to solve time-fractional Swift–Hohenberg equations. We take the linear and nonlinear Swift–Hohenberg equations with different initial conditions to illustrate the effectiveness of such a method. The results we get are displayed by solution graph for each problem. The present method has simple, accurate, and straightforward implementation to solve fractional-order Swift–Hohenberg equations. In conclusion, the suggested approach is considered a sophisticated tool for the solution of other fractional-order differential equations.

Author Contributions

Conceptualization, K.N., A.M.Z., A.M.A. and R.S.; investigation, K.N., A.M.Z., A.M.A. and R.S.; methodology, K.N., A.M.Z., Y.S.H. and R.S.; validation, K.N., A.M.Z. and R.S.; Formal Analysis, K.N., A.M.Z., A.K. and R.S.; Resources, A.M.Z., A.M.A. and R.S.; Data Curation, R.S.; Writing—Original Draft Preparation, K.N. and R.S.; Writing—Review and Editing, K.N. and R.S.; Project Administration, K.N.; Funding Acquisition, K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The numerical data used to support the findings of this study are included within the article.

Acknowledgments

One of the co-authors (A. M. Zidan) extends their appreciation to the Deanship of Scientific Research at King Khalid University, Abha 61413, Saudi Arabia, for funding this work through research groups program under grant number R.G.P.1/30/42. This Research was supported by Taif University Researchers Supporting Project Number (TURSP-2020/96), Taif University, Taif, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baleanu, D.; Guvenc, Z.B.; Machado, J.T. New Trends in Nanotechnology and Fractional Calculus Applications; Springer: New York, NY, USA, 2010. [Google Scholar]
  2. Baleanu, D.; Machado, J.A.; Luo, A.C. Fractional Dynamics and Control; Springer Science & Business Media: New York, NY, USA, 2011. [Google Scholar]
  3. Liu, Q.; Xu, Y.; Kurths, J. Active vibration suppression of a novel airfoil model with fractional order viscoelastic constitutive relationship. J. Sound Vib. 2018, 432, 50–64. [Google Scholar] [CrossRef]
  4. Xu, Y.; Li, Y.; Liu, D. A method to stochastic dynamical systems with strong nonlinearity and fractional damping. Nonlinear Dyn. 2016, 83, 2311–2321. [Google Scholar] [CrossRef]
  5. Xu, Y.; Li, Y.; Liu, D.; Jia, W.; Huang, H. Responses of Duffing oscillator with fractional damping and random phase. Nonlinear Dyn. 2013, 74, 745–753. [Google Scholar] [CrossRef]
  6. Caputo, M. Linear models of dissipation whose Q is almost frequency independent II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  7. Ford, N.J.; Simpson, A.C. The numerical solution of fractional differential equations: Speed versus accuracy. Numer. Algorithms 2001, 26, 333–346. [Google Scholar] [CrossRef]
  8. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  9. Shah, N.A.; Dassios, I.; Chung, J.D. A decomposition method for a fractional-order multi-dimensional telegraph equation via the Elzaki transform. Symmetry 2021, 13, 8. [Google Scholar] [CrossRef]
  10. Ryzhkov, S.V.; Kuzenov, V.V. New realization method for calculating convective heat transfer near the hypersonic aircraft surface. Z. Angew. Math. Phys. 2019, 70, 1–9. [Google Scholar] [CrossRef]
  11. Saadeh, R.; Qazza, A.; Burqan, A. A new integral transform: ARA transform and its properties and applications. Symmetry 2020, 12, 925. [Google Scholar] [CrossRef]
  12. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Fract. Differ. Appl. 2015, 2, 731–785. [Google Scholar]
  13. Losada, J.; Nieto, J.J. Properties of the new fractional derivative without singular kernel. Fract. Differ. Appl. 2015, 2, 87–92. [Google Scholar]
  14. Sunthrayuth, P.; Zidan, A.M.; Yao, S.-W.; Shah, R.; Inc, M. The Comparative Study for Solving Fractional-Order Fornberg–Whitham Equation via ρ-Laplace Transform. Symmetry 2021, 5, 784. [Google Scholar] [CrossRef]
  15. Baleanu, D.; Mustafa, O.G. On the global existence of solutions to a class of fractional differential equations. Comput. Math. Appl. 2010, 59, 35–41. [Google Scholar] [CrossRef] [Green Version]
  16. Yousef, F.; Alquran, M.; Jaradat, I.; Momani, S.; Baleanu, D. Ternary-fractional differential transform schema: Theory and application. Adv. Differ. Equ. 2019, 2019, 197. [Google Scholar] [CrossRef]
  17. Bokhari, A.; Baleanu, D.; Belgacem, R. Application of Shehu transform to Atangana-Baleanu derivatives. Int. J. Math. Comput. Sci. 2019, 20, 101–107. [Google Scholar] [CrossRef] [Green Version]
  18. He, J.H.; Ji, F.Y. Two-scale mathematics and fractional calculus for thermodynamics. Therm. Sci. 2019, 21, 2131–2133. [Google Scholar] [CrossRef]
  19. Wang, K.L.; Yao, S.W.; Yang, H.W. A fractal derivative model for snow’s thermal insulation property. Therm. Sci. 2019, 23, 2351–2354. [Google Scholar] [CrossRef]
  20. Kakutani, T.; Ono, H. Weak non-linear hydromagnetic waves in a cold collision-free plasma. J. Phys. Soc. Japan 1969, 26, 1305–1318. [Google Scholar] [CrossRef]
  21. Yang, X.J.; Srivastava, H.M.; Machado, J.A. A new fractional derivative without singular kernel: Application to the modelling of the steady heat flow. Therm. Sci. 2016, 20, 753–756. [Google Scholar] [CrossRef]
  22. Yang, X.J. Fractional derivatives of constant and variable orders applied to anomalous relaxation models in heat-transfer problems. Therm. Sci. 2017, 21, 1161–1171. [Google Scholar] [CrossRef] [Green Version]
  23. Singh, J.; Kumar, D.; Kumar, S. A new fractional model of nonlinear shock wave equation arising in flow of gases. Nonlinear Eng. 2014, 3, 43–50. [Google Scholar] [CrossRef]
  24. Naeem, M.; Zidan, A.M.; Nonlaopon, K.; Syam, M.I.; Al-Zhour, Z.; Shah, R. A New Analysis of Fractional-Order Equal-Width Equations via Novel Techniques. Symmetry 2021, 5, 886. [Google Scholar] [CrossRef]
  25. Paolo, D.B.; Fattorusso, L.; Versaci, M. Electrostatic field in terms of geometric curvature in membrane MEMS devices. Commun. Appl. Ind. Math. 2017, 8, 165–184. [Google Scholar]
  26. Yong, L.; Wang, H.; Chen, X.; Yang, X.; You, Z.; Dong, S.; Gao, J. Shear property, high-temperature rheological performance and low-temperature flexibility of asphalt mastics modified with bio-oil. Constr. Build. Mater. 2018, 174, 30–37. [Google Scholar]
  27. Meerschaert, M.M.; Tadjeran, C. Finite difference approximations for two-sided space-fractional partial differential equations. Appl. Numer. Math. 2006, 56, 80–90. [Google Scholar] [CrossRef]
  28. Ray, S.S.; Bera, R.K. Analytical solution of the Bagley Torvik equation by Adomian decomposition method. Appl. Math. Comput. 2005, 168, 398–410. [Google Scholar] [CrossRef]
  29. Jiang, Y.; Ma, J. High-order finite element methods for time-fractional partial differential equations. J. Comput. Appl. Math. 2011, 235, 3285–3290. [Google Scholar] [CrossRef] [Green Version]
  30. Odibat, Z.; Momani, S.; Erturk, V.S. Generalized differential transform method: Application to differential equations of fractional order. Appl. Math. Comput. 2008, 197, 467–477. [Google Scholar] [CrossRef]
  31. Arikoglu, A.; Ozkol, I. Solution of fractional differential equations by using differential transform method. Chaos Solitons Fractals 2007, 34, 1473–1481. [Google Scholar] [CrossRef]
  32. Zhang, X.; Zhao, J.; Liu, J.; Tang, B. Homotopy perturbation method for two dimensional time-fractional wave equation. Appl. Math. Model. 2014, 38, 5545–5552. [Google Scholar] [CrossRef]
  33. Prakash, A. Analytical method for space-fractional telegraph equation by homotopy perturbation transform method. Nonlinear Eng. 2016, 5, 123–128. [Google Scholar] [CrossRef]
  34. Dhaigude, C.; Nikam, V. Solution of fractional partial differential equations using iterative method. Fract. Calc. Appl. Anal. 2012, 15, 684–699. [Google Scholar] [CrossRef]
  35. Safari, M.; Ganji, D.D.; Moslemi, M. Application of He’s variational iteration method and Adomian’s decomposition method to the fractional KdV-Burgers-Kuramoto equation. Comput. Math. Appl. 2009, 58, 2091–2097. [Google Scholar] [CrossRef] [Green Version]
  36. Liao, S.J. The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problems. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 1992. [Google Scholar]
  37. Liao, S. Homotopy analysis method: A new analytical technique for nonlinear problems. Commun. Nonlinear Sci. Numer. Simulat. 1997, 2, 95–100. [Google Scholar] [CrossRef]
  38. Liao, S. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. 2004, 147, 499–513. [Google Scholar] [CrossRef]
  39. Abbasbandy, S.; Hashemi, M.S.; Hashim, I. On convergence of homotopy analysis method and its application to fractional integro-differential equations. Quaest. Math. 2013, 36, 93–105. [Google Scholar] [CrossRef]
  40. Kumar, D.; Singh, J.; Baleanu, D. A fractional model of convective radial fins with temperature-dependent thermal conductivity. Rom. Rep. Phys. 2017, 69, 103. [Google Scholar]
  41. Kumar, D.; Agarwal, R.P.; Singh, J. A modified numerical scheme and convergence analysis for fractional model of Lienards equation. J. Comput. Appl. Math. 2018, 339, 405–413. [Google Scholar] [CrossRef]
  42. Hang, X.; Cang, J. Analysis of a time fractional wave-like equation with the homotopy analysis method. Phys. Lett. A 2008, 372, 1250–1255. [Google Scholar]
  43. Dehghan, M.; Manafian, J.; Saadatmandi, A. The solution of the linear fractional partial differential equations using the homotopy analysis method. Z. Nat. A 2010, 65, 935–949. [Google Scholar] [CrossRef] [Green Version]
  44. Goufo, E.F.D.; Pene, M.K.; Mwambakana, J.N. Duplication in a model of rock fracture with fractional derivative without singular kernel. Open Math. 2015, 13, 839–846. [Google Scholar] [CrossRef] [Green Version]
  45. Jafari, H.; Das, S.; Tajadodi, H. Solving a multi-order fractional differential equation using homotopy analysis method. J. King Saud Univ. Sci. 2011, 23, 151–155. [Google Scholar] [CrossRef] [Green Version]
  46. Diethelm, K.; Ford, N.J. Multi-order fractional differential equations and their numerical solution. Appl. Math. Comput. 2004, 154, 621–640. [Google Scholar] [CrossRef]
  47. Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. 2006, 316, 753–763. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a). the solution of ETDM at different fractional-order δ. (b). the graph show that the close relation with each other.
Figure 1. (a). the solution of ETDM at different fractional-order δ. (b). the graph show that the close relation with each other.
Symmetry 13 01263 g001
Figure 2. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Figure 2. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Symmetry 13 01263 g002
Figure 3. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Figure 3. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Symmetry 13 01263 g003
Figure 4. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Figure 4. (a). δ = 0.25. (b). δ = 0.50. (c). δ = 0.75.
Symmetry 13 01263 g004
Figure 5. (a). the solution of ETDM at different fractional-order δ. (b). the graph show that the close relation with each other.
Figure 5. (a). the solution of ETDM at different fractional-order δ. (b). the graph show that the close relation with each other.
Symmetry 13 01263 g005
Figure 6. (a). the solution of ETDM at different fractional-order δ (b). the graph show that the close relation with each other.
Figure 6. (a). the solution of ETDM at different fractional-order δ (b). the graph show that the close relation with each other.
Symmetry 13 01263 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nonlaopon, K.; Alsharif, A.M.; Zidan, A.M.; Khan, A.; Hamed, Y.S.; Shah, R. Numerical Investigation of Fractional-Order Swift–Hohenberg Equations via a Novel Transform. Symmetry 2021, 13, 1263. https://doi.org/10.3390/sym13071263

AMA Style

Nonlaopon K, Alsharif AM, Zidan AM, Khan A, Hamed YS, Shah R. Numerical Investigation of Fractional-Order Swift–Hohenberg Equations via a Novel Transform. Symmetry. 2021; 13(7):1263. https://doi.org/10.3390/sym13071263

Chicago/Turabian Style

Nonlaopon, Kamsing, Abdullah M. Alsharif, Ahmed M. Zidan, Adnan Khan, Yasser S. Hamed, and Rasool Shah. 2021. "Numerical Investigation of Fractional-Order Swift–Hohenberg Equations via a Novel Transform" Symmetry 13, no. 7: 1263. https://doi.org/10.3390/sym13071263

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop