Next Article in Journal
A New Blow-Up Criterion to a Singular Non-Newton Polytropic Filtration Equation
Next Article in Special Issue
Elucidating the Effects of Ionizing Radiation on Immune Cell Populations: A Mathematical Modeling Approach with Special Emphasis on Fractional Derivatives
Previous Article in Journal
Hamilton–Jacobi Inequality Adaptive Robust Learning Tracking Controller of Wearable Robotic Knee System
Previous Article in Special Issue
A Novel Approach to Solving Fractional-Order Kolmogorov and Rosenau–Hyman Models through the q-Homotopy Analysis Transform Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing Both Fractional Porous Media and Heat Transfer Equations via Some Novel Techniques

1
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan
3
Department of Mechanical Engineering, Sejong University, Seoul 05006, Republic of Korea
4
Department of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
5
Department of Physics, Faculty of Science, Ain Shams University, Cairo 11566, Egypt
6
Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt
7
Research Center for Physics (RCP), Department of Physics, Faculty of Science and Arts, Al-Mikhwah, Al-Baha University, Al-Baha 1988, Saudi Arabia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work and are co-first authors.
Mathematics 2023, 11(6), 1350; https://doi.org/10.3390/math11061350
Submission received: 10 February 2023 / Revised: 6 March 2023 / Accepted: 6 March 2023 / Published: 10 March 2023

Abstract

:
It has been increasingly obvious in recent decades that fractional calculus (FC) plays a key role in many disciplines of applied sciences. Fractional partial differential equations (FPDEs) accurately model various natural physical phenomena and many engineering problems. For this reason, the analytical and numerical solutions to these issues are seriously considered, and different approaches and techniques have been presented to address them. In this work, the FC is applied to solve and analyze the time-fractional heat transfer equation as well as the nonlinear fractional porous media equation with cubic nonlinearity. The idea of solving these equations is based on the combination of the Yang transformation (YT), the homotopy perturbation method (HPM), and the Adomian decomposition method (ADM). These combinations give rise to two novel methodologies, known as the homotopy perturbation transform method (HPTM) and the Yang tranform decomposition method (YTDM). The obtained results show the significance of the accuracy of the suggested approaches. Solutions in various fractional orders are found and discussed. It is noted that solutions at various fractional orders lead to an integer-order solution. The application of the current methodologies to other nonlinear fractional issues in other branches of applied science is supported by their straightforward and efficient process. In addition, the proposed solution methods can help many plasma physics researchers in interpreting the theoretical and practical results.

1. Introduction

Fractional calculus (FC) has a significant advantage when it comes to simulating physical problems with whole-memory effects. In applied mathematics, fractional derivatives are a useful tool for relating non-integer-order derivatives as well as integrals. Numerous studies have been conducted in many fields of science, physics, and engineering concerning this branch due to its importance in modeling many realistic problems. In reality, the novel features of this distinguished branch appear in several disciplines, including fluid dynamics, chemistry, applied mathematics, physics of plasma, astrophysics, mathematical biology, control theory, image processing, controlled thermonuclear fusion, and many other fields. Fractional-order integrals and derivatives have become a significant new mathematical tool that helps in solving many problems that arise in engineering, physics, and nature. Podlubny [1] provided comprehensive details on arbitrary-order differential equations (DEs). A great book on the basics and details of FC and DEs was written by Miller and Ross [2]. There are also many research works on this very important topic in modeling several nonlinear phenomena in physics, especially the physics of plasmas. Singh [3] investigated a fractional-order blood alcohol model. Caputo presented the basic characteristics of the FC [4]. The topic of singular kernels connected to fractional derivatives was studied by Caputo, and Fabrizio [5]. New aspects of fractional derivatives were explained by Caputo, and Fabrizio [4].
Both nonlinear ordinary DEs (ODEs) and partial DEs (PDEs) are primarily applied to describe most physical phenomena, either simple or complex, that appear in reality [6,7,8] via mathematical modeling, such as many nonlinear phenomena that have been studied and analyzed in different models of plasma physics [9,10,11,12,13,14]. The mathematical modeling of numerous natural processes benefits greatly from the use of nonlinear fractional DEs (NLFDEs). The essential feature of fractional derivatives is their nonlocality, which emphasizes how the future state arises from both current and all previous states [15,16,17]. The investigation of nonlinear ODEs is an essential way of determining the behavior of the system, which has attracted the interest of many researchers (including scientists and engineers). Fractional PDEs are significant mathematical models that accurately represent a variety of complex processes in numerous scientific disciplines, such as fractional diffusion equations occurring in oil pollution [18], fractional Fornberg–Whitham equations [19], fractional Helmholtz equations [20], local fractional modified Zakharov–Kuznetsov (mZK) equations [21], fractional local damped wave and dissipative equations [22], fractional nonlinear Boussinesq Equations [23], time-fractional Navier–Stokes equations [24], fractional third-order Burgers and Korteweg–De Vries (KdV) equations [25], nonlinear Schrodinger (NLS) equations [26], Maxwell’s equations [27], fractional advection–dispersion equations [27], time-fractional Fisher’s equations [28], nonlinear biological predator–prey population systems [29], fractional-order Kaup–Kupershmidt equations [30], and fractional-order Gardner and Cahn–Hilliard equations [31].
Mathematicians have constructed many numerical and analytical strategies for solving and analyzing integer order DEs and FPDEs because of the aforementioned applications [32,33,34,35,36]. For example, Mehdi Dehghan et al., in [37], applied the homotopy analysis method (HAM) for solving the fractional KdV, Burgers, BBM–Burgers, coupled KdV, cubic Boussinesq, and Boussinesq-like equations. Mohamed et al. [38] applied Elzaki Transformation HPM (ETHPM) to determine a series form of solutions of different NLFDEs. Similarly, Mounirah Areshi et al. [39] applied the variational iteration transform method (VITM) to solve fractional-order Newell–Whitehead–Segel equations. An optimal HAM has been applied by Randhir Singh [40] to solve non-isothermal reaction–diffusion model equations with a spherical catalyst. With the help of the modified trial equation method, Hasan Bulut and Yusuf Pandir [41] solved the time-fractional Sharma–Tasso–Olever (STO) equation. Sheng Zhang et al. [42] used the modified extended Fan sub-equation technique to solve the (3+1)-dimensional Kadomstev–Petviashvili (KP) equation. Analytic results of the fractional-order diffusion equation have been investigated by applying the Yang transformation (YT) decomposition technique [43]. Furthermore, the numerical solutions of the multi-dimensional time-fractional telegraphic equation have been analyzed and discussed using the reduced differential transform method [44]. Moreover, Azzh Saad Alshehry et al. [45] developed natural transformation based on the decomposition method for analyzing a fractional Kuramoto–Sivashinsky equation.
In nonlinear flow studies concerning heat and mass transfer, diffusion, porous media flow, biology, and viscosity, as well as other related subjects, boundary-layer theory is a crucial component; for example, the porous media equation, sometimes known as a nonlinear heat equation, commonly exists, and more effort is needed to analyze them via simplified mathematical approaches [46,47,48,49]. Current research aims to present novel methods for solving the fractional porous media equation and the non-linear time-fractional heat transfer equation that features cubic non-linearity using the combination between YT, HPM, and ADM. Adomian introduced the ADM in 1980 as a powerful tool for resolving differential equations describing physical phenomena [50]. ADM is additionally altered using the YT. As a result, the Yang transformation decomposition method (YTDM), a novel technique, has been constructed.
The fundamental elements of both YT and ADM are included in the YTDM. This method is useful and effective for analyzing PDEs of initial and boundary value problems. He [51,52] created HPM by combining conventional homotopy and perturbation to address a variety of physical issues. It is important to note that the HPM is applied without any discretization, restricted assumptions, or transformation and is error-free concerning rounding. Due to the challenges raised by the nonlinear variables, the YT is completely unable to handle nonlinear equations. The well-known YT method is paired with the homotopy perturbation approach to develop a highly efficient method known as the homotopy perturbation transform method (HPTM) for dealing with numerous nonlinear issues. The advantages of the proposed techniques are their capability of combining two powerful approaches for obtaining exact and approximate analytical solutions to nonlinear equations. It is important to keep in mind that the suggested methods can perform better in general as they need less computational work than the alternative methods while keeping a high level of numerical precision; the size reduction amounts to an improvement in the performance of the approaches.
The following summarises this work’s structure: In Section 2, we present a few fundamental elements of calculus theory. In Section 3 and Section 4, the general solution is obtained using the HPTM and YTDM algorithms, and Section 5 provides the convergence analysis of these suggested techniques. In Section 6, we demonstrate the reliability and efficacy of both strategies. In Section 7, the conclusion is presented.

2. Preliminary Concepts

Definition 1.
The The Caputo definition of order ς is described as [39,53]
D ϑ ς ζ ( ß ) = 1 Γ ( ı ς ) 0 ß ( ß j ) ı ς 1 ζ ( ı ) ( j ) d j , ı 1 < ς ı , ı N .
Definition 2.
YT is defined by [43,54]
Y { ζ ( ß ) } = M ( u ) = 0 e ß u ζ ( ß ) d ß , ß > 0 ,
and its inverse reads
Y 1 { M ( u ) } = ζ ( ß ) .
Definition 3.
The YT inverse reads [43,54]
Y 1 [ Y u ] = ζ ( ß ) = 1 2 π ι ς ι ς + ι ζ 1 u e u ß u d u = Σ r e s i d u e s o f ζ 1 u e u ß u .
Definition 4.
The derivative of fractional order in YT reads [43,54]
Y { ζ ς ( ß ) } = M ( u ) u ς ı = 0 n 1 ζ ı ( 0 ) u ς ( ı + 1 ) , n 1 < ς n .

3. Fundamental Plan of HPTM

To discuss the main plan of the suggested scheme, we examine the following FPDE related to the Caputo derivative:
D ß ς ζ ( χ , ß ) = M 1 [ χ ] ζ ( χ , ß ) + N 1 [ χ ] ζ ( χ , ß ) , 0 < ς 1 ,
with initial solution
ζ ( χ , 0 ) = ξ ( χ ) .
Here, D ß ς = ς ß ς denotes the fractional derivative in a Caputo manner, and M 1 [ χ ] and N 1 [ χ ] represent the general differential operators.
Applying the YT, we obtain
Y [ D ß ς ζ ( χ , ß ) ] = Y [ M 1 [ χ ] ζ ( χ , ß ) + N 1 [ χ ] ζ ( χ , ß ) ] ,
1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y [ M 1 [ χ ] ζ ( χ , ß ) + N 1 [ χ ] ζ ( χ , ß ) ] .
The simplified form reads
M ( u ) = u ζ ( χ , 0 ) + u ς Y [ M 1 [ χ ] ζ ( χ , ß ) + N 1 [ χ ] ζ ( χ , ß ) ] .
By employing the inverse of YT, we have
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 [ u ς Y [ M 1 [ χ ] ζ ( χ , ß ) + N 1 [ χ ] ζ ( χ , ß ) ] ] .
In terms of HPM, the basic solution in a power series reads
ζ ( χ , ß ) = ı = 0 ϵ ı ζ ı ( χ , ß ) .
with homotopy parameter ϵ [ 0 , 1 ] .
The nonlinear term reads
N 1 [ χ ] ζ ( χ , ß ) = ı = 0 ϵ ı H ı ( ζ ) .
Accordingly, the following polynomials are obtained
H ı ( ζ 0 , ζ 1 , , ζ n ) = 1 Γ ( n + 1 ) D ϵ ı N 1 ı = 0 ϵ i ζ i ϵ = 0 ,
where D ϵ ı = ı ϵ ı .
Substituting Equations (10) and (11) into Equation (9), we obtain
ı = 0 ϵ ı ζ ı ( χ , ß ) = ζ ( χ , 0 ) + ϵ × Y 1 u ς Y { M 1 ı = 0 ϵ ı ζ ı ( χ , ß ) + ı = 0 ϵ ı H ı ( ζ ) } .
Equating similar coefficients of ϵ in the above equation, we have
ϵ 0 : ζ 0 ( χ , ß ) = ζ ( χ , 0 ) , ϵ 1 : ζ 1 ( χ , ß ) = Y 1 u ς Y ( M 1 [ χ ] ζ 0 ( χ , ß ) + H 0 ( ζ ) ) , ϵ 2 : ζ 2 ( χ , ß ) = Y 1 u ς Y ( M 1 [ χ ] ζ 1 ( χ , ß ) + H 1 ( ζ ) ) , . . . ϵ ı : ζ ı ( χ , ß ) = Y 1 u ς Y ( M 1 [ χ ] ζ ı 1 ( χ , ß ) + H ı 1 ( ζ ) ) , ı > 0 , ı N .
Finally, we approximate the analytical solution in the following series form:
ζ ( χ , ß ) = lim M ı = 1 M ζ ı ( χ , ß ) .

4. The Fundamental Idea of the YTDM

The principal scheme of the current approach is based on studying the FPDE associated with the Caputo derivative as follows:
D ß ς ζ ( χ , ß ) = M 1 ( χ , ß ) + N 1 ( χ , ß ) , 0 < ς 1 ,
with initial source
ζ ( χ , 0 ) = ξ ( χ ) .
Here, the fractional derivative in the Caputo manner is denoted by D ß ς = ς ß ς , and the general differential operators are represented by M 1 and N 1 .
Applying the YT, we obtain
Y [ D ß ς ζ ( χ , ß ) ] = Y [ M 1 ( χ , ß ) + N 1 ( χ , ß ) ] , 1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y [ M 1 ( χ , ß ) + N 1 ( χ , ß ) ] .
The simplified form reads
M ( u ) = u ζ ( χ , 0 ) + u ς Y [ M 1 ( χ , ß ) + N 1 ( χ , ß ) ] ,
By employing the inverse of the YT, we obtain
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 [ u ς Y [ M 1 ( χ , ß ) + N 1 ( χ , ß ) ] .
The series form solution of ζ ( χ , ß ) reads
ζ ( χ , ß ) = m = 0 ζ m ( χ , ß ) .
In Equation (19), the nonlinear term is put in as
N 1 ( χ , ß ) = m = 0 A m ( ζ ) .
with
A m ζ 0 , ζ 1 , ζ 2 , , ζ m = 1 m ! m m N 1 m = 0 m ζ m = 0 , m = 0 , 1 , 2 ,
The substitution of Equations (20) and (21) into Equation (19) yields
m = 0 ζ m ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y M 1 m = 0 ζ m ( χ , ß ) + m = 0 A m ( ζ ) .
Similarly,
ζ 0 ( χ , ß ) = ζ ( χ , 0 ) ,
ζ 1 ( χ , ß ) = Y 1 u ς Y { M 1 ( ζ 0 ) + A 0 } ,
Generally, for m 1 , we obtain
ζ m + 1 ( χ , ß ) = Y 1 u ς Y { M 1 ( ζ m ) + A m } .

5. Convergence Analysis

In this part, the suggested techniques of convergence analysis are discussed
Theorem 1.
Suppose that the exact solution to Equation (5) is G ( χ , ß ) , and let G ( χ , ß ) , G n ( χ , ß ) H and α ( 0 , 1 ) , where H represents the Hilbert space. The obtained solution q = 0 G q ( χ , ß ) converges G ( χ , ß ) if G q ( χ , ß ) G q 1 ( χ , ß ) q > A , i.e., for any ω > 0 A > 0 , such that | | G q + n ( χ , ß ) | | β , m , n N .
Proof. 
We take a sequence of q = 0 G q ( χ , ß ) .
C 0 ( χ , ß ) = G 0 ( χ , ß ) , C 1 ( χ , ß ) = G 0 ( χ , ß ) + G 1 ( χ , ß ) , C 2 ( χ , ß ) = G 0 ( χ , ß ) + G 1 ( χ , ß ) + G 2 ( χ , ß ) , C 3 ( χ , ß ) = G 0 ( χ , ß ) + G 1 ( χ , ß ) + G 2 ( χ , ß ) + G 3 ( χ , ß ) , C q ( χ , ß ) = G 0 ( χ , ß ) + G 1 ( χ , ß ) + G 2 ( χ , ß ) + + G q ( χ , ß ) ,
We must demonstrate that C q ( χ , ß ) forms a “Cauchy sequence” in order to achieve the desired outcome. Additionally, let us take
| | C q + 1 ( χ , ß ) C q ( χ , ß ) | | = | | G q + 1 ( χ , ß ) | | α | | G q ( χ , ß ) | | α 2 | | G q 1 ( χ , ß ) | | α 3 | | G q 2 ( χ , ß ) | | α q + 1 | | G 0 ( χ , ß ) | | .
For q , n N , we have
| | C q ( χ , ß ) C n ( χ , ß ) | | = | | G q + n ( χ , ß ) | | = | | C q ( χ , ß ) C q 1 ( χ , ß ) + ( C q 1 ( χ , ß ) C q 2 ( χ , ß ) ) + ( C q 2 ( χ , ß ) C q 3 ( χ , ß ) ) + + ( C n + 1 ( χ , ß ) C n ( χ , ß ) ) | | | | C q ( χ , ß ) C q 1 ( χ , ß ) | | + | | ( C q 1 ( χ , ß ) C q 2 ( χ , ß ) ) | | + | | ( C q 2 ( χ , ß ) C q 3 ( χ , ß ) ) | | + + | | ( C n + 1 ( χ , ß ) C n ( χ , ß ) ) | | α q | | G 0 ( χ , ß ) | | + α q 1 | | G 0 ( χ , ß ) | | + + α q + 1 | | G 0 ( χ , ß ) | | = | | G 0 ( χ , ß ) | | ( α q + α q 1 + α q + 1 ) = | | G 0 ( χ , ß ) | | 1 α q n 1 α q + 1 α n + 1 .
As 0 < α < 1 , and G 0 ( χ , ß ) are bound, we take β = 1 α / ( 1 α q n ) α n + 1 | | G 0 ( χ , ß ) | | , and we obtain
| | G q + n ( χ , ß ) | | β , q , n N .
Hence, { G q ( χ , ß ) } q = 0 makes a “Cauchy sequence” in H. It proves that the sequence { G q ( χ , ß ) } q = 0 is a convergent sequence with the limit lim q G q ( χ , ß ) = G ( χ , ß ) for G ( χ , ß ) H , which completes the proof. □
Theorem 2.
Let us assume that h = 0 k G h ( χ , ß ) is finite and that G ( χ , ß ) reflects the series solution that was found. Assuming α > 0 such that | | G h + 1 ( χ , ß ) | | | | G h ( χ , ß ) | | , the maximum absolute error is given by the following relation.
| | G ( χ , ß ) h = 0 k G h ( χ , ß ) | | < α k + 1 1 α | | G 0 ( χ , ß ) | | .
Proof. 
Suppose h = 0 k G h ( χ , ß ) is finite which implies that h = 0 k G h ( χ , ß ) < .
Let us consider
| | G ( χ , ß ) h = 0 k G h ( χ , ß ) | | = | | h = k + 1 G h ( χ , ß ) | | h = k + 1 | | G h ( χ , ß ) | | h = k + 1 α h | | G 0 ( χ , ß ) | | α k + 1 ( 1 + α + α 2 + ) | | G 0 ( χ , ß ) | | α k + 1 1 α | | G 0 ( χ , ß ) | | .
which completes the proof of theorem. □
Theorem 3.
The result of (16) is unique when 0 < ( φ 1 + φ 2 ) ( ß ς Γ ( ς + 1 ) ) < 1 .
Proof. 
Let H = ( C [ J ] , | | . | | ) with the norm | | ϕ ( ß ) | | = m a x ß J | ϕ ( ß ) | be Banach space, which is ∀-continuous function on J. Let I : H H be a non-linear mapping, where
ζ l + 1 C = ζ 0 C + Y 1 [ p ( ς , υ , ϖ ) Y [ M ( ζ l ( χ , ß ) ) + N ( ζ l ( χ , ß ) ) ] ] , l 0 .
Suppose that | M ( ζ ) M ( ζ * ) | < φ 1 | ζ ζ * | and | N ( ζ ) N ( ζ * ) | < φ 2 | ζ ζ * | , where ζ : = ζ ( χ , ß ) and ζ * : = ζ * ( χ , ß ) are are two different function values and φ 1 , φ 2 are Lipschitz constants.
| | I ζ I ζ * | | m a x t J | Y 1 [ p ( ς , υ , ϖ ) Y [ M ( ζ ) M ( ζ * ) ] + p ( ς , υ , ϖ ) Y [ N ( ζ ) N ( ζ * ) ] | ] m a x ß J [ φ 1 Y 1 [ p ( ς , υ , ϖ ) Y [ | ζ ζ * | ] ] + φ 2 Y 1 [ p ( ς , υ , ϖ ) Y [ | ζ ζ * | ] ] ] m a x t J ( φ 1 + φ 2 ) Y 1 [ p ( ς , υ , ϖ ) Y | ζ ζ * | ] ( φ 1 + φ 2 ) Y 1 [ p ( ς , υ , ϖ ) Y | | ζ ζ * | | ] = ( φ 1 + φ 2 ) ( ß ς Γ ( ς + 1 ) ) | | ζ ζ * | | ,
where I is a contraction as 0 < ( φ 1 + φ 2 ) ( ß ς Γ ( ς + 1 ) ) < 1 . From Banach’s fixed-point theorem, the result of (16) is unique. □
Theorem 4.
The result of (16) is convergent.
Proof. 
Let ζ m = r = 0 m ζ r ( χ , ß ) . To show that ζ m is a Cauchy sequence in H, let
| | ζ m ζ n | | = m a x ß J | r = n + 1 m ζ r | , n = 1 , 2 , 3 , m a x ß J Y 1 p ( ς , υ , ϖ ) Y r = n + 1 m ( M ( ζ r 1 ) + N ( ζ r 1 ) ) = m a x ß J Y 1 p ( ς , υ , ϖ ) Y r = n + 1 m 1 ( M ( ζ r ) + N ( ζ r ) ) m a x ß J | Y 1 [ p ( ς , υ , ϖ ) Y [ ( M ( ζ m 1 ) M ( ζ n 1 ) + N ( ζ m 1 ) N ( ζ n 1 ) ) ] ] | φ 1 m a x ß J | Y 1 [ p ( ς , υ , ϖ ) Y [ ( M ( ζ m 1 ) M ( ζ n 1 ) ) ] ] | + φ 2 m a x ß J | Y 1 [ p ( ς , υ , ϖ ) Y [ ( N ( ζ m 1 ) N ( ζ n 1 ) ) ] ] | = ( φ 1 + φ 2 ) ( ß ς Γ ( ς + 1 ) ) | | ζ m 1 ζ n 1 | | .
Let m = n + 1 ; then,
| | ζ n + 1 ζ n | | φ | | ζ n ζ n 1 | | φ 2 | | ζ n 1 ζ n 2 | | φ n | | ζ 1 ζ 0 | | ,
where φ = ( φ 1 + φ 2 ) ( ß ς Γ ( ς + 1 ) ) . Similarly, we have
| | ζ m ζ n | | | | ζ n + 1 ζ n | | + | | ζ n + 2 ζ n + 1 | | + + | | ζ m ζ m 1 | | , ( φ n + φ n + 1 + + φ m 1 ) | | ζ 1 ζ 0 | | φ n 1 φ m n 1 φ | | ζ 1 | | ,
As 0 < φ < 1 , we obtain 1 φ m n < 1 . Therefore,
| | ζ m ζ n | | φ n 1 φ m a x ß J | | ζ 1 | | .
Since | | ζ 1 | | < , | | ζ m ζ n | | 0 when n . As a result, ζ m is a Cauchy sequence in H, implying that the series ζ m is convergent. □

6. Numerical Problems

Example 1.
Here, we consider the following fractional porous media Equation [55]:
ς ζ ( χ , ß ) ß ς = χ ζ ( χ , ß ) ζ ( χ , ß ) χ ,
with the initial condition
ζ ( χ , 0 ) = χ
Applying the YT yields
Y ς ζ ß ς = Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ ,
The simplified form reads
1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ ,
M ( u ) = u ζ ( χ , 0 ) + u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ .
By employing the inverse of the YT, we obtain
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ , ζ ( χ , ß ) = χ + Y 1 u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ .
In terms of HPM, we have
ı = 0 ϵ ı ζ ı ( χ , ß ) = χ + Y 1 u ς Y χ ı = 0 ϵ ı H ı ( ζ ) .
Here, the nonlinear term reads
H 0 ( ζ ) = ζ 0 ζ 0 χ , H 1 ( ζ ) = ζ 0 ζ 1 χ + ζ 1 ζ 0 χ , H 2 ( ζ ) = ζ 0 ζ 2 χ + ζ 1 ζ 1 χ + ζ 2 ζ 0 χ , .
Now, by equating the coefficient of ϵ, we have
ϵ 0 : ζ 0 ( χ , ß ) = χ , ϵ 1 : ζ 1 ( χ , ß ) = Y 1 u ς Y χ ı = 0 ϵ ı H 0 ( ζ ) = ß ς Γ ( ς + 1 ) , ϵ 2 : ζ 2 ( χ , ß ) = Y 1 u ς Y χ ı = 0 ϵ ı H 1 ( ζ ) = 0 , .
Finally, we approximate the analytical solution in the form of the following series"
ζ ( χ , ß ) = ζ 0 ( χ , ß ) + ζ 1 ( χ , ß ) + ζ 2 ( χ , ß ) + , ζ ( χ , ß ) = χ + ß ς Γ ( ς + 1 ) + .
Now, in terms of YTDM and by applying the YT, we obtain
Y ς ζ ß ς = Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ ,
The simplified form reads
1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ ,
M ( u ) = u ζ ( χ , 0 ) + u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ .
By employing the inverse of the YT, we obtain
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ , ζ ( χ , ß ) = ( 15 8 sech 2 ( χ 2 ) ) + Y 1 u ς Y χ ζ ( χ , ß ) ζ ( χ , ß ) χ .
The series form solution of ζ ( χ , ß ) reads
ζ ( χ , ß ) = m = 0 ζ m ( χ , ß ) ,
with ζ ( χ , ß ) ζ ( χ , ß ) χ = m = 0 A m , showing the nonlinear term in terms of Adomian polynomial, and
m = 0 ζ m ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y χ m = 0 A m , m = 0 ζ m ( χ , ß ) = χ + Y 1 u ς Y χ m = 0 A m .
Here, the nonlinear terms read
A 0 = ζ 0 ζ 0 χ , A 1 = ζ 0 ζ 1 χ + ζ 1 ζ 0 χ , A 2 = ζ 0 ζ 2 χ + ζ 1 ζ 1 χ + ζ 2 ζ 0 χ , .
Similarly,
ζ 0 ( χ , ß ) = χ ,
for m = 0
ζ 1 ( χ , ß ) = ß ς Γ ( ς + 1 ) ,
for m = 1
ζ 2 ( χ , ß ) = 0 ,
Finally, we approximate the analytical solution in the form of the following series:
ζ ( χ , ß ) = m = 0 ζ m ( χ , ß ) = ζ 0 ( χ , ß ) + ζ 1 ( χ , ß ) + ζ 2 ( χ , ß ) + ,
ζ ( χ , ß ) = χ + ß ς Γ ( ς + 1 ) + .
The exact result at ς = 1 reads
ζ ( χ , ß ) = χ + ß .
Example 2.
Here, we consider the following fractional heat-transfer equation [55]
ς ζ ( χ , ß ) ß ς = 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) ,
with the initial condition
ζ ( χ , 0 ) = 1 + 2 χ χ 2 + χ + 1 .
Applying the YT yields
Y ς ζ ß ς = Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) ,
The simplified form reads
1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) ,
M ( u ) = u ζ ( χ , 0 ) + u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) .
By employing the inverse of the YT, we obtain
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) , ζ ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 + Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) .
In terms of HPM, we have
ı = 0 ϵ ı ζ ı ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 + Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 ı = 0 ϵ ı H ı ( ζ ) .
Here, the nonlinear terms read
H 0 ( ζ ) = 2 ζ 0 3 ( χ , ß ) , H 1 ( ζ ) = 6 ζ 0 2 ( χ , ß ) ζ 1 ( χ , ß ) , H 2 ( ζ ) = 6 ζ 0 ( χ , ß ) ζ 1 2 ( χ , ß ) + 6 ζ 0 2 ( χ , ß ) ζ 2 ( χ , ß ) , .
Now, by equating the coefficient of ϵ, we have
ϵ 0 : ζ 0 ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 , ϵ 1 : ζ 1 ( χ , ß ) = Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 ı = 0 ϵ ı H 0 ( ζ ) = 6 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 2 ß ς Γ ( ς + 1 ) , ϵ 2 : ζ 2 ( χ , ß ) = Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 ı = 0 ϵ ı H 1 ( ζ ) = 72 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 3 ß 2 ς Γ ( 2 ς + 1 ) , ϵ 3 : ζ 3 ( χ , ß ) = Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 ı = 0 ϵ ı H 2 ( ζ ) = ( 1296 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 4 + 432 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 216 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 · Γ ( 2 ς + 1 ) Γ 2 ( ς + 1 ) ) ß 3 ς Γ ( 3 ς + 1 ) , .
Finally, we approximate the analytical solution in the form of the following series
ζ ( χ , ß ) = ζ 0 ( χ , ß ) + ζ 1 ( χ , ß ) + ζ 2 ( χ , ß ) + ζ 3 ( χ , ß ) + ζ ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 + 6 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 2 ß ς Γ ( ς + 1 ) + 72 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 3 ß 2 ς Γ ( 2 ς + 1 ) + ( 1296 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 4 + 432 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 216 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 · Γ ( 2 ς + 1 ) Γ 2 ( ς + 1 ) ) ß 3 ς Γ ( 3 ς + 1 ) + .
Now, in terms of YTDM and applying the YT yields
Y ς ζ ß ς = Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) .
Upon simplification, we have
1 u ς { M ( u ) u ζ ( χ , 0 ) } = Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) ,
M ( u ) = u ζ ( χ , 0 ) + u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) .
By employing the inverse of YT, we obtain
ζ ( χ , ß ) = ζ ( χ , 0 ) + Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) , ζ ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 + Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 2 ζ 3 ( χ , ß ) .
The series form solution of ζ ( χ , ß ) reads
ζ ( χ , ß ) = m = 0 ζ m ( χ , ß ) ,
with 2 ζ 3 ( χ , ß ) = m = 0 A m , showing the nonlinear term in terms of Adomian polynomial
m = 0 ζ m ( χ , ß ) = ζ ( χ , 0 ) Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 m = 0 A m , m = 0 ζ m ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 Y 1 u ς Y 2 ζ ( χ , ß ) χ 2 m = 0 A m .
Here, the nonlinear terms read
A 0 = 2 ζ 0 3 ( χ , ß ) , A 1 = 6 ζ 0 2 ( χ , ß ) ζ 1 ( χ , ß ) , A 2 = 6 ζ 0 ( χ , ß ) ζ 1 2 ( χ , ß ) + 6 ζ 0 2 ( χ , ß ) ζ 2 ( χ , ß ) , .
Similarly,
ζ 0 ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 ,
for m = 0
ζ 1 ( χ , ß ) = 6 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 2 ß ς Γ ( ς + 1 ) ,
for m = 1
ζ 2 ( χ , ß ) = 72 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 3 ß 2 ς Γ ( 2 ς + 1 ) ,
for m = 2
ζ 3 ( χ , ß ) = 1296 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 4 + 432 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 216 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 · Γ ( 2 ς + 1 ) Γ 2 ( ς + 1 ) ß 3 ς Γ ( 3 ς + 1 ) .
Finally, we approximate the analytical solution in the form of the following series:
ζ ( χ , ß ) = m = 0 ζ m ( χ , ß ) = ζ 0 ( χ , ß ) + ζ 1 ( χ , ß ) + ζ 2 ( χ , ß ) + ζ 3 ( χ , ß ) + ,
ζ ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 + 6 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 2 ß ς Γ ( ς + 1 ) + 72 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 3 ß 2 ς Γ ( 2 ς + 1 ) + ( 1296 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 4 + 432 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 216 ( 1 + 2 χ ) 3 ( χ 2 + χ + 1 ) 5 · Γ ( 2 ς + 1 ) Γ 2 ( ς + 1 ) ) ß 3 ς Γ ( 3 ς + 1 ) + .
The obtained results agree with the result in Ref. [56]. For ς = 1 , the above can be rearranged as
ζ ( χ , ß ) = 1 + 2 χ χ 2 + χ + 1 6 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 2 ß + 36 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 3 ß 2 216 ( 1 + 2 χ ) ( χ 2 + χ + 1 ) 4 ß 3 + .
Numerical Simulation Studies
Here, we have presented the numerical simulations of the time-fractional heat transfer equation as well as the nonlinear fractional porous media equation with cubic nonlinearity of fractional order by applying the proposed methodologies. We used graphs and tables to show how the obtained solution behaves. Maple was used to complete all of the computational work for the problems stated. The behavior of the exact and suggested approaches’ solutions at ς = 1 is depicted in the graphs in Figure 1. Figure 2 displays the results of proposed methodologies at various fractional orders of ς = 1 , 0.75 , 0.50 , 0.25 , within the domain 1 χ , ß 0 for Example 1. Table 1 compares the accurate as well as proposed methods’ solutions with the aid of absolute error for Example 1. The approximate solutions ζ ( χ , ß ) for various fractional-order ς values are shown in Figure 3 and Figure 4, respectively. This illustrates how the graphical behavior varies on the order of the fractional derivative. Table 2 compares the approximations for the solutions to Example 2 for various values of χ and ß. The demonstrated plots help us to better understand the nature of the proposed equations when temporal-spatial variables vary in comparison with arbitrary order. The comparison of absolute errors demonstrates that our methods converge more quickly than other methods. Additionally, the graphical depiction demonstrates a good agreement between the exact solution and the suggested approaches solution.

7. Conclusions

In this work, both the nonlinear fractional porous media equation and time-fractional heat transfer equation have been solved by implementing HPTM and YTDM, an elegant combination of the hybrid Yang transformation (YT) with the homotopy perturbation method (HPM) and the Adomian decomposition method (ADM). Both He’s polynomials and Adomian polynomials have been applied to express the nonlinear terms in the targeted issues. The suggested hybrid methods feature simpler and clearer steps for determining the solutions to fractional problems. To make the offered approaches more understandable and to assess their applicability, various numerical problems were solved. The obtained results demonstrate a strong association between the suggested strategy and the actual solution. Plots of the fractional solutions depict the behavior of different dynamics of the specified physical phenomenon. The fractional solution rapidly approaches the integer-order ones. In conclusion, we recommend that these methods be applied to solving several other nonlinear fractional problems in various scientific fields because they are easy to implement and lead to the actual solution. The proposed methods can also address many evolution equations that govern different nonlinear phenomena in plasma physics [57,58,59,60].

Author Contributions

Conceptualization, W.A. and R.S.; methodology, R.S. and S.A.E.-T.; software, R.S. and S.A.E.-T.; validation, W.A., N.A.S., J.D.C. and S.M.E.I.; formal analysis, W.A., N.A.S., J.D.C. and S.M.E.I.; investigation, S.A.E.-T.; resources, R.S. and S.A.E.-T.; data curation, W.A., N.A.S., J.D.C. and S.M.E.I.; writing—original draft preparation, R.S.; writing—review and editing, W.A., S.M.E.I. and S.A.E.-T.; visualization, R.S., N.A.S. and J.D.C.; supervision, R.S. and S.A.E.-T.; project administration, S.A.E.-T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

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

Acknowledgments

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R157), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This study is supported via funding from Prince Sattam bin Abdulaziz University project number (PSAU/2023/R/1444). This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2017R1D1A1B05030422).

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this article.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999; Volume 198, p. 340. [Google Scholar]
  2. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Willey: New York, NY, USA, 1993. [Google Scholar]
  3. Singh, J. Analysis of fractional blood alcohol model with composite fractional derivative. Chaos Solitons Fractals 2020, 140, 110127. [Google Scholar] [CrossRef]
  4. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  5. Caputo, M.; Fabrizio, M. On the Singular Kernels for Fractional Derivatives. Some Applications to Partial Differential Equations. Progr. Fract. Differ. Appl. 2021, 7, 79–82. [Google Scholar]
  6. Lu, C.; Zhu, R.; Yu, F.; Jiang, X.; Liu, Z.; Dong, L.; Hua, Q.; Ou, Z. Gear rotational speed sensor based on FeCoSiB/Pb(Zr,Ti)O3 magnetoelectric composite. Meas. J. Int. Meas. Confed. 2021, 168, 108409. [Google Scholar] [CrossRef]
  7. Lei, W.; Hui, Z.; Xiang, L.; Zelin, Z.; Xu-Hui, X.; Evans, S. Optimal Remanufacturing Service Resource Allocation for Generalized Growth of Retired Mechanical Products: Maximizing Matching Efficiency. IEEE Access 2021, 9, 89655–89674. [Google Scholar] [CrossRef]
  8. Zhang, X.; Ma, F.; Dai, Z.; Wang, J.; Chen, L.; Ling, H.; Soltanian, M.R. Radionuclide transport in multi-scale fractured rocks: A review. J. Hazard. Mater. 2022, 424 Pt C, 127550. [Google Scholar] [CrossRef]
  9. Alharthi, M.R.; Alharbey, R.A.; El-Tantawy, S.A. Novel analytical approximations to the nonplanar Kawahara equation and its plasma applications. Eur. Phys. J. Plus 2022, 137, 1172. [Google Scholar] [CrossRef]
  10. Alhejaili, W.; Salas, A.H.; El-Tantawy, S.A. Approximate solution to a generalized Van der Pol equation arising in plasma oscillations. AIP Adv. 2022, 12, 105104. [Google Scholar] [CrossRef]
  11. El-Tantawy, S.A.; Salas, A.H.; Alyousef, H.A.; Alharthie, M.R. Novel approximations to a nonplanar nonlinear Schrödinger equation and modeling nonplanar rogue waves/breathers in a complex plasma. Chaos Solitons Fractals 2022, 163, 112612. [Google Scholar] [CrossRef]
  12. Aljahdaly, N.H.; El-Tantawy, S.A. Novel anlytical solution to the damped Kawahara equation and its application for modeling the dissipative nonlinear structures in a fluid medium. J. Ocean Eng. Sci. 2022, 7, 492–497. [Google Scholar] [CrossRef]
  13. Alyousef, H.A.; Salas, A.H.; Matoog, R.T.; El-Tantawy, S.A. On the analytical and numerical approximations to the forced damped Gardner Kawahara equation and modeling the nonlinear structures in a collisional plasma. Phys. Fluids 2022, 34, 103105. [Google Scholar] [CrossRef]
  14. Douanla, D.V.; Tiofack, C.G.L.; Alim; Aboubakar, M.; Mohamadou, A.; Albalawi, W.; El-Tantawy, S.A.; El-Sherif, L.S. Three-dimensional rogue waves and dust-acoustic dark soliton collisions in degenerate ultradense magnetoplasma in the presence of dust pressure anisotropy. Phys. Fluids 2022, 34, 087105. [Google Scholar] [CrossRef]
  15. Bai, B.; Zhou, R.; Yang, G.; Zou, W.; Yuan, W. The constitutive behavior and dissociation effect of hydrate-bearing sediment within a granular thermodynamic framework. Ocean Eng. 2023, 268, 113408. [Google Scholar] [CrossRef]
  16. Liu, K.; Yang, Z.; Wei, W.; Gao, B.; Xin, D.; Sun, C.; Wu, G. Novel detection approach for thermal defects: Study on its feasibility and application to vehicle cables. High Volt. 2022, 1–10. [Google Scholar] [CrossRef]
  17. Wang, X.; Lyu, X. Experimental study on vertical water entry of twin spheres side-by-side. Ocean Eng. 2021, 221, 108508. [Google Scholar] [CrossRef]
  18. Singh, J.; Alshehri, A.M.; Momani, S.; Hadid, S.; Kumar, D. Computational Analysis of Fractional Diffusion Equations Occurring in Oil Pollution. Mathematics 2022, 10, 3827. [Google Scholar] [CrossRef]
  19. Gupta, P.K.; Singh, M. Homotopy perturbation method for fractional Fornberg-Whitham equation. Comput. Math. Appl. 2011, 61, 250–254. [Google Scholar] [CrossRef]
  20. Abuasad, S.; Moaddy, K.; Hashim, I. Analytical treatment of two-dimensional fractional Helmholtz equations. J. King Saud Univ.-Sci. 2019, 31, 659–666. [Google Scholar] [CrossRef]
  21. Wang, K. A novel perspective to the local fractional Zakharov-Kuznetsov-modified equal width dynamical model on Cantor sets. Math. Methods Appl. Sci. 2022, 46, 622–630. [Google Scholar] [CrossRef]
  22. Dubey, V.P.; Singh, J.; Alshehri, A.M.; Dubey, S.; Kumar, D. A hybrid computational method for local fractional dissipative and damped wave equations in fractal media. Waves Random Complex Media 2022, 1–23. [Google Scholar] [CrossRef]
  23. Alyobi, S.; Shah, R.; Khan, A.; Shah, N.A.; Nonlaopon, K. Fractional Analysis of Nonlinear Boussinesq Equation under Atangana-Baleanu-Caputo Operator. Symmetry 2022, 14, 2417. [Google Scholar] [CrossRef]
  24. Jassim, H.K.; Zayir, M.Y. A Fractional Variational Iteration Approach for Solving Time-Fractional Navier-Stokes Equations. Math. Comput. Sci. 2022, 3, 41–47. [Google Scholar]
  25. Alderremy, A.A.; Aly, S.; Fayyaz, R.; Khan, A.; Shah, R.; Wyal, N. The analysis of fractional-order nonlinear systems of third order KdV and Burgers equations via a novel transform. Complexity 2022, 2022, 4935809. [Google Scholar] [CrossRef]
  26. Abdel-Salam, E.A.; Yousif, E.A.; El-Aasser, M.A. Analytical solution of the space-time fractional nonlinear Schrödinger equation. Rep. Math. Phys. 2016, 77, 19–34. [Google Scholar] [CrossRef]
  27. Huang, Q.; Huang, G.; Zhan, H. A finite element solution for the fractional advection-dispersion equation. Adv. Water Resour. 2008, 31, 1578–1589. [Google Scholar] [CrossRef]
  28. Zidan, A.M.; Khan, A.; Shah, R.; Alaoui, M.K.; Weera, W. Evaluation of time-fractional Fisher’s equations with the help of analytical methods. AIMS Math. 2022, 7, 18746–18766. [Google Scholar] [CrossRef]
  29. Singh, J.; Kilicman, A.; Kumar, D.; Swroop, R. Numerical study for fractional model of nonlinear predator-prey biological population dynamic system. Ther. Sci. 2018, 23, 366. [Google Scholar] [CrossRef]
  30. Shah, N.A.; Hamed, Y.S.; Abualnaja, K.M.; Chung, J.D.; Shah, R.; Khan, A. A comparative analysis of fractional-order kaup-kupershmidt equation within different operators. Symmetry 2022, 14, 986. [Google Scholar] [CrossRef]
  31. Kbiri, M.; Nonlaopon, K.; Zidan, A.M.; Khan, A.; Shah, R. Analytical investigation of fractional-order cahn-hilliard and gardner equations using two novel techniques. Mathematics 2022, 10, 1643. [Google Scholar] [CrossRef]
  32. Dahmani, Z.; Anber, A.; Gouari, Y.; Kaid, M.; Jebril, I. Extension of a Method for Solving Nonlinear Evolution Equations via Conformable Fractional Approach. In Proceedings of the 2021 International Conference on Information Technology, ICIT 2021, Guangzhou, China, 22–25 December 2021; pp. 38–42. [Google Scholar]
  33. Hammad, M.A.; Horani, M.A.; Shmasenh, A.; Khalil, R. Ruduction of order of fractional differential equations. J. Math. Comput. Sci. 2018, 8, 683–688. [Google Scholar]
  34. Hammad, M.A. Conformable Fractional Martingales and Some Convergence Theorems. Mathematics 2021, 10, 6. [Google Scholar] [CrossRef]
  35. Zhang, X.; Ma, F.; Yin, S.; Wallace, C.D.; Soltanian, M.R.; Dai, Z.; Lu, X. Application of upscaling methods for fluid flow and mass transport in multi-scale heterogeneous media: A critical review. Appl. Energy 2021, 303, 117603. [Google Scholar] [CrossRef]
  36. Bai, B.; Rao, D.; Chang, T.; Guo, Z. A nonlinear attachment-detachment model with adsorption hysteresis for suspension-colloidal transport in porous media. J. Hydrol. (Amst.) 2019, 578, 124080. [Google Scholar] [CrossRef]
  37. Dehghan, M.; Manafian, J.; Saadatm, A. Solving nonlinear fractional partial differential equations using the homotopy analysis method. Numer. Methods Partial. Differ. Equ. Int. J. 2010, 26, 448–479. [Google Scholar] [CrossRef]
  38. Mohamed, M.; Yousif, M.; Hamza, A.E. Solving nonlinear fractional partial differential equations using the Elzaki transform method and the homotopy perturbation method. Abstr. Appl. Anal. 2022, 2022, 4743234. [Google Scholar] [CrossRef]
  39. Areshi, M.; Khan, A.; Shah, R.; Nonlaopon, K. Analytical investigation of fractional-order Newell-Whitehead-Segel equations via a novel transform. AIMS Math. 2022, 7, 6936–6958. [Google Scholar] [CrossRef]
  40. Singh, R. Optimal homotopy analysis method for the non-isothermal reaction-diffusion model equations in a spherical catalyst. J. Math. Chem. 2018, 56, 2579–2590. [Google Scholar] [CrossRef]
  41. Bulut, H.; Pandir, Y. Modified trial equation method to the nonlinear fractional Sharma-Tasso-Olever equation. Int. J. Model. Optim. 2013, 3, 353. [Google Scholar] [CrossRef] [Green Version]
  42. Zhang, S.; Xia, T. A further improved extended Fan sub-equation method and its application to the (3+ 1)-dimensional Kadomstev-Petviashvili equation. Phys. Lett. A 2006, 356, 119–123. [Google Scholar] [CrossRef]
  43. Sunthrayuth, P.; Alyousef, H.A.; El-Tantawy, S.A.; Khan, A.; Wyal, N. Solving Fractional-Order Diffusion Equations in a Plasma and Fluids via a Novel Transform. J. Funct. Spaces 2022, 2022, 1899130. [Google Scholar] [CrossRef]
  44. Srivastava, V.K.; Awasthi, M.K.; Kumar, S. Analytical approximations of two and three dimensional time-fractional telegraphic equation by reduced differential transform method. Egypt. J. Basic Appl. Sci. 2014, 1, 60–66. [Google Scholar] [CrossRef] [Green Version]
  45. Saad Alshehry, A.; Imran, M.; Khan, A.; Shah, R.; Weera, W. Fractional View Analysis of Kuramoto-Sivashinsky Equations with Non-Singular Kernel Operators. Symmetry 2022, 14, 1463. [Google Scholar] [CrossRef]
  46. Arasteh, H.; Mashayekhi, R.; Toghraie, D.; Karimipour, A.; Bahiraei, M.; Rahbari, A. Optimal arrangements of a heat sink partially filled with multilayered porous media employing hybrid nanofluid. J. Therm. Anal. Calorim. 2019, 137, 1045–1058. [Google Scholar] [CrossRef]
  47. Barnoon, P.; Toghraie, D.; Dehkordi, R.B.; Afr, M. Two phase natural convection and thermal radiation of Non-Newtonian nanofluid in a porous cavity considering inclined cavity and size of inside cylinders. Int. Commun. Heat Mass Transf. 2019, 108, 104285. [Google Scholar] [CrossRef]
  48. Arasteh, H.; Mashayekhi, R.; Ghaneifar, M.; Toghraie, D.; Afrand, M. Heat transfer enhancement in a counter-flow sinusoidal parallel-plate heat exchanger partially filled with porous media using metal foam in the channels divergent sections. J. Therm. Anal. Calorim. 2020, 141, 1669–1685. [Google Scholar] [CrossRef]
  49. Martín Ortega, E.B.; Sastre, F.; Velázquez, A.; Baïri, A. Heat transfer enhancement around finned vertical antenna by means of porous media saturated with Water-Copper nanofluid. Case Stud. Ther. Eng. 2021, 28, 101555. [Google Scholar] [CrossRef]
  50. Adomian, G. Solving frontier problems of physics: The decomposition method. In Fundamental Theories of Physics; Cherruault, Y., Ed.; Kluwer Academic Publishers Group: Dordrecht, IL, USA, 1994; Volume 1. [Google Scholar]
  51. He, J.H. Homotopy perturbation method: A new nonlinear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar] [CrossRef]
  52. He, J.H. Comparison of homotopy perturbation method and homotopy analysis method. Appl. Math. Comput. 2004, 156, 527–539. [Google Scholar] [CrossRef]
  53. Sunthrayuth, P.; Ullah, R.; Khan, A.; Shah, R.; Kafle, J.; Mahariq, I.; Jarad, F. Numerical analysis of the fractional-order nonlinear system of Volterra integro-differential equations. J. Funct. Spaces 2021, 2021, 1537958. [Google Scholar] [CrossRef]
  54. Alaoui, M.K.; Fayyaz, R.; Khan, A.; Shah, R.; Abdo, M.S. Analytical investigation of NoyesField model for time-fractional BelousovZhabotinsky reaction. Complexity 2021, 2021, 3248376. [Google Scholar] [CrossRef]
  55. Ganji, D.D.; Sadighi, A. Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations. J. Comput. Appl. Math. 2007, 207, 24–34. [Google Scholar] [CrossRef] [Green Version]
  56. Yan, L.M. Modified Homotopy Perturbation Method Coupled with Laplace Transform for Fractional Heat Transfer and Porous Media Equations. Therm. Sci. 2013, 17, 1409–1414. [Google Scholar] [CrossRef] [Green Version]
  57. El-Tantawy, S.A.; El-Sherif, L.S.; Bakry, A.M.; Alhejaili, W.; Wazwaz, A. On the analytical approximations to the nonplanar damped Kawahara equation: Cnoidal and solitary waves and their energy. Phys. Fluids 2022, 34, 113103. [Google Scholar] [CrossRef]
  58. El-Tantawy, S.A.; Salas, A.H.; Alyouse, H.A.; Alharthi, M.R. Novel exact and approximate solutions to the family of the forced damped Kawahara equation and modeling strong nonlinear waves in a plasma. Chin. J. Phys. 2022, 77, 2454–2471. [Google Scholar] [CrossRef]
  59. El-Tantawy, S.A. Nonlinear dynamics of soliton collisions in electronegative plasmas: The phase shifts of the planar KdV-and mkdV-soliton collisions. Chaos Solitons Fractals 2016, 93, 162–168. [Google Scholar] [CrossRef]
  60. Ali, S.; Shohaib, M.; Masood, W.; Alyousef, H.A.; El-Tantawy, S.A. The attributes of the dust-acoustic solitary and periodic structures in the Saturn’s inner magnetosphere. Phys. Fluids 2023, 35, 023101. [Google Scholar] [CrossRef]
Figure 1. The visual representation of the proposed techniques and the specific solution at ς = 1 for ζ ( χ , ß ) are presented with a graphical layout.
Figure 1. The visual representation of the proposed techniques and the specific solution at ς = 1 for ζ ( χ , ß ) are presented with a graphical layout.
Mathematics 11 01350 g001
Figure 2. The proposed solution methods for ζ ( χ , ß ) have been visually represented in the form of a graphical layout for different values of ς .
Figure 2. The proposed solution methods for ζ ( χ , ß ) have been visually represented in the form of a graphical layout for different values of ς .
Mathematics 11 01350 g002
Figure 3. The visual representation of the proposed techniques at a value of ς = 0.7 and ς = 0.9 for the function ζ ( χ , ß ) is shown graphically.
Figure 3. The visual representation of the proposed techniques at a value of ς = 0.7 and ς = 0.9 for the function ζ ( χ , ß ) is shown graphically.
Mathematics 11 01350 g003
Figure 4. Visual representation of proposed techniques for finding ζ ( χ , ß ) at different values of ς .
Figure 4. Visual representation of proposed techniques for finding ζ ( χ , ß ) at different values of ς .
Mathematics 11 01350 g004
Table 1. A comparison between our approach and the actual outcome with a ς value of 1, including the calculation of the absolute error (AE).
Table 1. A comparison between our approach and the actual outcome with a ς value of 1, including the calculation of the absolute error (AE).
ß = 0.01 Exact ResultOur Techniques’ ResultAE of Our Techniques
χ ς = 1 ς = 1 ς = 1
11.010000000001.010000000000.0000000000  × 10 + 00
22.010000000002.010000000000.0000000000  × 10 + 00
33.010000000003.010000000000.0000000000  × 10 + 00
44.010000000004.010000000000.0000000000  × 10 + 00
55.010000000005.010000000000.0000000000  × 10 + 00
66.010000000006.010000000000.0000000000  × 10 + 00
77.010000000007.010000000000.0000000000  × 10 + 00
88.010000000008.010000000000.0000000000  × 10 + 00
99.010000000009.010000000000.0000000000  × 10 + 00
1010.0100000000010.010000000000.0000000000  × 10 + 00
Table 2. A solution is presented through the proposed method, which involves the use of different values of χ and ß for varying fractional orders ς in Example 2.
Table 2. A solution is presented through the proposed method, which involves the use of different values of χ and ß for varying fractional orders ς in Example 2.
χ ß ς = 0.50 ς = 0.75 ς = 1
10.250.9005770.7963500.750000
2 0.5773300.5831040.594023
3 0.4679020.4755400.483500
4 0.3905280.3954510.400145
5 0.3325370.3356390.338499
10.501.6126251.2559791.000000
2 0.5789170.5529330.539358
3 0.4381410.4392830.442876
4 0.3692070.3723290.376093
5 0.3186210.3211250.323822
10.752.8256670.3211251.750000
2 0.6553980.5959200.550291
3 0.4277560.4205330.416590
4 0.3550170.3551640.356413
5 0.3079950.3091670.310806
114.4485323.6923643.000000
2 0.79108880.7037480.626822
3 0.4320960.4168960.404642
4 0.3460440.3429910.341107
5 0.2996960.2992870.299452
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Albalawi, W.; Shah, R.; Shah, N.A.; Chung, J.D.; Ismaeel, S.M.E.; El-Tantawy, S.A. Analyzing Both Fractional Porous Media and Heat Transfer Equations via Some Novel Techniques. Mathematics 2023, 11, 1350. https://doi.org/10.3390/math11061350

AMA Style

Albalawi W, Shah R, Shah NA, Chung JD, Ismaeel SME, El-Tantawy SA. Analyzing Both Fractional Porous Media and Heat Transfer Equations via Some Novel Techniques. Mathematics. 2023; 11(6):1350. https://doi.org/10.3390/math11061350

Chicago/Turabian Style

Albalawi, Wedad, Rasool Shah, Nehad Ali Shah, Jae Dong Chung, Sherif M. E. Ismaeel, and Samir A. El-Tantawy. 2023. "Analyzing Both Fractional Porous Media and Heat Transfer Equations via Some Novel Techniques" Mathematics 11, no. 6: 1350. https://doi.org/10.3390/math11061350

APA Style

Albalawi, W., Shah, R., Shah, N. A., Chung, J. D., Ismaeel, S. M. E., & El-Tantawy, S. A. (2023). Analyzing Both Fractional Porous Media and Heat Transfer Equations via Some Novel Techniques. Mathematics, 11(6), 1350. https://doi.org/10.3390/math11061350

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