Abstract
Fractional differential equations can model various complex problems in physics and engineering, but there is no universal method to solve fractional models precisely. This paper offers a new hope for this purpose by coupling the homotopy perturbation method with Aboodh transform. The new hybrid technique leads to a simple approach to finding an approximate solution, which converges fast to the exact one with less computing effort. An example of the fractional casting-mold system is given to elucidate the hope for fractional calculus, and this paper serves as a model for other fractional differential equations.
1 Introduction
Fractional calculus has triggered much interest in both physics and mathematics [1, 2]. Traditional differential equations cannot accurately represent many physical problems, and the fractional partner can provide deeper insight into these complex physical phenomena with ease. In general, this newly developed field is for studying real-world applications in the fractal space, so most literature labeled it as the fractal–fractional calculus [3–5] or the local fractional calculus on the Cantor set [6]. A continuum medium, e.g., water or air, becomes a fractal space (porous medium) when we observe it on a molecule’s scale. Any phenomena arising in molecules’ perturbation have to be modeled by the fractal–fractional model [7]. As an example, we consider a nanoparticle’s motion in the air, which is stochastic and difficult to be modeled by the traditional differential equation; however, if the air is considered as a fractal space on a molecule’s scale, its motion is determinate and can be modeled by the fractal–fractional model. So, we need two scales for a porous medium; one is large enough so that the continuum assumption works, and the other is small enough so that the porosity can be measured, as pointed out by Ji-Huan He that “seeing with a single scale is always unbelieving” [8]. Another example is the motion of the Moon, which is naturally periodic; however, if we measure its motion at an extremely far distance, its motion becomes stochastic, and the Heisenberg-like uncertainty principle works for the Moon [9]. He and Qian showed that the fractal diffusion process in water depends on the fractal dimensions [10], and other scientists also discussed the fractional advection–reaction–diffusion process [11] and the fractal diffusion–reaction process [12]. A cocoon’s air/moisture permeability and its thermal property can best be revealed by the fractal–fractional model [13, 14], and the fractal micro-electromechanical systems show even more amazing properties [15–18].
Fractional calculus is a good and reliable tool for scientists and engineers but a mixed blessing for practical applications because an intractable problem arises; that is, fractional models are extremely difficult to be solved. Researchers have been racing to test various analytical methods which were originally proposed to solve traditional differential equations. Though there are many famous analytical methods in the literature, for example, the homotopy perturbation method [19–23] and its various modifications [24–26], the decomposition method [27], the variational iteration method [28–30], the exp-function method [31], and the differential transform method [32], so far, there is not a universal approach to solving exactly fractional differential equations, and this paper offers a new hope for this purpose by coupling the homotopy perturbation method [33, 34] and the Aboodh transform [35].
The homotopy perturbation method (HPM) was first proposed by Chinese mathematician Prof. Ji-Huan He in the later 1990s [33]; it is mathematically simple and physically insightful. The method is equally suitable for linear or non-linear, homogeneous or inhomogeneous, and initial and/or boundary value problems. The solution is expressed in an infinite series and typically converges to the exact solution. The HPM is now considered a matured tool for almost all kinds of problems, and many researchers have used this method for an accurate insight into the solution properties of a complex problem [36–38].
The Aboodh transform (AT) was proposed by Aboodh [35] and derived from the classical Fourier integral. This transform is now considered a simple technique for solving linear differential equations but is unable to solve non-linear ones. By coupling AT with the HPM, one has the capability to solve linear and non-linear problems, and a lot of literature works have been witnessed to utilize this coupling for solving various types of problems. Using AT–HPM, Manimegalai et al. [39] solved strongly non-linear oscillators with great success. Jani and Singh [40] found it had obvious advantages over the decomposition method, Yasmin [41] revealed the dynamic behavior of the fractional convection–reaction–diffusion process, and Jani and Singh [42] extended it to the soliton theory.
Though much work was achieved, in this study, we will show that AT–HPM is a universal tool for fractional calculus. As an example, we consider the time-fractional casting-mold system which is used in manufacturing various medical equipment, ranging from injections to the COVID-19 tool-kit [43]. The significant findings reveal that AT–HPM is an accurate and effective approach that reduces the computational work with fast convergence ratio.
2 Aboodh transform-based homotopy perturbation technique
This section is divided into two sections. In the first section, the methodology will be proposed, and the convergence of the suggested technique will be discussed in the second section.
2.1 Methodology
In this section, we give a brief introduction to the Aboodh transform [35] and homotopy perturbation method [33, 34].
If is a continuous piecewise function of time , then the Aboodh transform of is that can be expressed as follows [35]:where and are positive and can be finite or infinite. is considered a function of the exponential order, which assures the convergence of the integrand. is the kernel of the transform, and is the transform variable. Table 1 includes the Aboodh transformation of some elementary functions helpful for this manuscript. This table can also be used for inverse Aboodh transform.
TABLE 1
Aboodh transform of some elementary functions.
The Aboodh transform of the partial derivative of time can be obtained using the following formula:where is the independent variable. Now, suppose the general system of PDEs is expressed aswhere is the linear operator, are the non-linear operators, are the dependent variables, and are the inhomogeneous functions. We assume the initial conditions aswhere and are known functions of the independent variable . The methodology composed of initially applying the Aboodh transform to both sides of the system of equations written in Eq. 3 and then employing the given initial conditions expressed in Eq. 4, thus yielding
By employing the differential characteristic of Aboodh transform, we can express Eq. 3 asand after using the initial conditions, we haveorwhere and denote the terms arising from the initial condition. According to the standard homotopy perturbation method [33, 34], the solution and can be expanded into an infinite series aswhere is the embedding parameter. Also, the non-linear terms and can be written aswhere and are He’s polynomials [44] and can be generated by the recursive formula
By substituting Eqs 7, 8 in Eq. 6, we get
Comparing the coefficients of like powers of , we have
We can obtain the best approximation for the solution as
2.2 Convergence analysis
To show that the series solution of the system in Eq. 14 converges to the solution of Eq. 3, we are to prove the sufficient condition of the convergence, and the following theorem will help us.
Theorem: We assume that X and Y are Banach spaces and is a non-linear contractive mapping such that
Then, according to Banach’s fixed point theorem,
Mhas a unique fixed point
, that is,
Supposing that the sequence in Eq.
14can be written as
and considering that
, we have
(i)
(ii)
Proof: (i) By the principle of mathematical induction, for , we have
Assuming as an induction hypothesis, we get
By employing the definition of , we have
(ii) As ,
Hence, the given statement is proved.
3 Numerical examples
In this section, three examples are presented to illustrate the idea explained in Section 2. First, we will study the method for a homogeneous linear system of PDEs. Second, the analytical solution will be obtained for an inhomogeneous linear system of PDEs. Finally, the inhomogeneous non-linear system of PDEs will be examined.
3.1 The system of homogeneous linear PDEs
We consider the following linear system:with initial conditions
By employing the Aboodh transform method, we have
Using the initial conditions given in Eq. 16, we reachor
The Aboodh transform-based homotopy perturbation method considers a series solution given by
By using the aforestated equation, the system of equations in Eq. 19 gets the form
By comparing like powers of from the aforestated equation, we obtain
Hence, the series solution by using Eq. 14 can be expressed asor in a closed form aswhich is the exact solution of Eq. 15.
3.2 The system of inhomogeneous linear PDEs
Suppose the following inhomogeneous linear system of PDEs:with initial conditions
Applying the Aboodh transform on each side of the equations in Eq. 24 and then putting on the given initial conditions, we obtainor
By using the Aboodh transform-based homotopy perturbation method, the series solution is expressed by
The system of equations in Eq. 27 gets the following form after employing the aforestated equation:
By comparing the coefficient of like powers of , we have
Therefore, the solution in the form of an infinite series by using Eq. 14 can be expressed asor in its convergent form aswhich is the exact solution of Eq. 24.
3.3 The system of inhomogeneous non-linear PDEs
Suppose the following inhomogeneous non-linear system of PDEs:with initial conditions
Employing the Aboodh transform on each side of the equations in Eq. 36 and then applying the given initial conditions give
Taking the inverse Aboodh transform on each side, we obtain
According to the Aboodh transform-based homotopy perturbation method, the solution functions and are series solutions, and inserting these series into both sides of each equation of the system yieldswhere the non-linear terms and are denoted by He’s polynomials and , respectively. A few He’s polynomials are
By comparing the coefficient of like powers of , we have
Therefore, the solution in the form of an infinite series by using Eq. 14 can be expressed asor in its convergent form aswhich is the exact solution of Eq. 36.
4 Time-fractional casting-mold system
Now, we turn back to a time-fractional casting-mold system which models the temperature distribution in the casting and molding processes. For this, two heat conduction equations are used with initial and Dirichlet boundary conditions [45]. The mathematical model is depicted as follows:where are parameters, are functions of time and space that represent the temperature on casting and molding plates, respectively, and is the fractal dimension. For more details on the modeling aspect of the aforementioned model, readers can see [45].
It is necessary to point out that Eq. 48 was originally studied in [45], where the series solution was presented and no closed-form solution was formulated. Our aim here is to overcome the main shortcomings in [45] and to offer a totally new hope for numerical approximation. To this end, applying the Aboodh transform in the aforementioned system, we have
Now, by inverse Aboodh transformation, we obtainwhich can further be written as
According to the standard HPM [33, 34], the solution and can be expanded into a finite series as
By substituting Eq. 52 in Eq. 51, the solution can be written as
Equating coefficients of powers of , we yield the following:
The approximate solution can be obtained as
4.1 Example
We consider the system expressed in Eq. 48 for the case . By utilizing Eqs 54–56, we have
By employing Eq. 57, the solution can be written as
The expressions are similar to those obtained by the fractional complex transform [46–49]. In the closed form, we obtainwhere is the Mittag-Leffler function [50]. One can check that Eq. 59 is an exact solution of Eq. 48 for the said parameters.
4.2 Results and discussion
This section is devoted to test the applicability and validity of the suggested technique for the time-fractional casting-mold system over the series-based solution of the same model.
Figures 1, 2 present the errors of the series solutions obtained by the HPM [45] for the fractal dimension β = 1. It is observed that for all the parameters and for both casting and molding processes, the errors grow exponentially for the case of a series solution [45] and can be reduced by adding more terms in the solution. On the other hand, the suggested solution has the exact solution, and there is no chance of error even for a larger range of t. Therefore, based on these findings, we can say that the proposed technique is more effective than the previous method [45].
FIGURE 1
FIGURE 2
5 Conclusion
The Aboodh transform-based homotopy perturbation method is successfully employed to solve traditional differential equations and fractional differential equations successfully. This approach has been shown to have the potential to solve both linear and non-linear problems. For a linear system, the exact solution is predicted, while for a non-linear system, with the help of He’s polynomials, a series solution is obtained, which converges fast to the exact one. So, the method pushes the progress of non-linear science and will make a “big change” to increase the number of practical applications, and this paper serves as a model for other applications.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding author.
Author contributions
Conceptualization: HT and NA; methodology: NA and YY; validation: NA and YY; writing—original draft preparation: HT and YY; writing—review and editing: HT, NA, and YY; supervision: HT and YY; and funding acquisition: YY. All authors read and agreed to the published version of the manuscript.
Funding
The study was supported by the Natural Science Foundation of Henan Province (No. 222300420507); National Natural Science Foundation of China (No. 12171193), Key Scientific Research Project of High Education Institutions of Henan Province (No. 23A110019), Science and Technology Research Projects of Henan Province (No. 182102110292), Basic and Frontier Technology Research Project of Henan Province (Nos. 12300410398 and 132300410084), and Zhumadian Key Laboratory of Statistical Computing and Data Modeling [No. (2022)12].
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
| continuous piecewise function | |
| Aboodh transform operator | |
| independent variable | |
| Linear operator | |
| dependent variable | |
| purturbation parameter | |
| inverse Aboodh transform operator | |
| mapping from X to Y | |
| parameter | |
| temperature at casting plate | |
| Mittag-Leffler function | |
| fractal dimension | |
| time | |
| transformed variable | |
| functions of independent variables | |
| Nonlinear operators | |
| functions of variable | |
| He’s polynomials | |
| Banach spaces | |
| fixed point | |
| elements of Banach space | |
| temperature at molding plate | |
| Gamma function | |
| parameters of casting and molding |
References
1.
HongBJ. Exact solutions for the conformable fractional coupled nonlinear Schrodinger equations with variable coefficients. J Low Frequency Noise, Vibration Active Control (2023) 146134842211354. 10.1177/14613484221135478
2.
HeJHHeCHSaeedT. A fractal modification of Chen-Lee-Liu equation and its fractal variational principle. Int J Mode Phys B (2021) 35(21):2150214. 10.1142/s0217979221502143
3.
WangKLWeiCF. Fractal soliton solutions for the fractal-fractional shallow water wave equation arising in ocean engineering. Alexandria Eng J (2023) 65(2023):859–65. 10.1016/j.aej.2022.10.024
4.
FengGQNiuJY. An analytical solution of the fractal toda oscillator. Results Phys (2023) 44:106208. 10.1016/j.rinp.2023.106208
5.
HeJHQieNHeCH. Solitary waves travelling along an unsmooth boundary. Results Phys (2021) 24:104104. 10.1016/j.rinp.2021.104104
6.
YangXJSrivastavaHMHeJHBaleanuD. Cantor-type cylindrical-coordinate method for differential equations with local fractional derivatives. Phys Lett A (2013) 377(28-30):1696–700. 10.1016/j.physleta.2013.04.012
7.
HeJH. Fractal calculus and its geometrical explanation. Results Phys (2018) 10:272–6. 10.1016/j.rinp.2018.06.011
8.
HeJH. Seeing with a single scale is always unbelieving: From magic to two-scale fractal. Therm Sci (2021) 25(2):1217–9. 10.2298/tsci2102217h
9.
HeJH. Frontier of modern textile engineering and short remarks on some topics in physics. Int J Nonlinear Sci Numer Simulation (2010) 11(7):555–63. 10.1515/ijnsns.2010.11.7.555
10.
HeJHQianMY. A fractal approach to the diffusion process of red ink in a saline water. Therm Sci (2022) 26(3B):2447–51. 10.2298/tsci2203447h
11.
DaiDDBanTTWangYLZhangW. The piecewise reproducing kernel method for the time variable fractional order advection-reaction-diffusion equations. Therm Sci (2021) 25(2B):1261–8. 10.2298/tsci200302021d
12.
LinLQiaoY. Fractal diffusion-reaction model for a porous electrode. Therm Sci (2021) 25(2):1305–11. 10.2298/tsci191212026l
13.
LiuFJZhangTHeCHTianD. Thermal oscillation arising in a heat shock of a porous hierarchy and its application. Facta Universitatis Ser Mech Eng (2022) 20(3):633–45. 10.22190/fume210317054l
14.
XueRJLiuFJ. A Fractional model and its application to heat prevention coating with cocoon-like hierarchy. Therm Sci (2022) 26(3):2493–8. 10.2298/tsci2203493x
15.
TianDAinQTAnjumNHeCHChengB. Fractal N/MEMS: From pull-in instability to pull-in stability. Fractals (2021) 29:2150030. 10.1142/s0218348x21500304
16.
TianDHeCH. A fractal micro-electromechanical system and its pull-in stability. J Low Frequency Noise Vibration Active Control (2021) 40(3):1380–6. 10.1177/1461348420984041
17.
HeJHYangQHeCHLiHBBuheE. Pull-in stability of a fractal MEMS system and its pull-in plateau. Fractals (2023) 30. 10.1142/S0218348X22501857
18.
HeCH. A variational principle for a fractal nano/microelectromechanical (N/MEMS) system. Int J Numer Methods Heat Fluid Flow (2023) 33(1):351–9. 10.1108/hff-03-2022-0191
19.
AnjumNHeJH. Higher-order homotopy perturbation method for conservative nonlinear oscillators generally and microelectromechanical systems’ oscillators particularly. Int J Mod Phys B (2020) 34:2050313. 10.1142/S0217979220503130
20.
HeCHEl-DibYO. A heuristic review on the homotopy perturbation method for non-conservative oscillators. J Low Frequency Noise, Vibration Active Control (2022) 41(2):572–603. 10.1177/14613484211059264
21.
AnjumNHeJH. Homotopy perturbation method for N/MEMS oscillators. Math Methods Appl Sci (2020). 10.1002/mma.6583
22.
HeJHJiaoMLHeCH. Homotopy perturbation method for fractal Duffing oscillator with arbitrary conditions. Fractals (2023) 30. 10.1142/S0218348X22501651
23.
HeJHMoatimidGMZekryMH. Forced nonlinear oscillator in a fractal space. Facta Universitatis Ser Mech Eng (2022) 20(1):001–20. 10.22190/fume220118004h
24.
LiXXHeCH (2019). Homotopy perturbation method coupled with the enhanced perturbation method, J Low Frequency Noise Vibration Active Control38, 1399–403. 10.1177/1461348418800554
25.
AnjumNHeJHAinQTTianD. Li-He’s modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system. Facta Universitatis Ser Mech Eng (2021) 19(4):601–12. 10.22190/fume210112025a
26.
HeJHEl-DibYO. The enhanced homotopy perturbation method for axial vibration of strings. Facta Universitatis Ser Mech Eng (2021) 19(4):735–50. 10.22190/fume210125033h
27.
WazwazAM. The decomposition method applied to systems of partial differential equations and to the reaction–diffusion Brusselator model. Appl Math Comput (2000) 110:251–64. 10.1016/s0096-3003(99)00131-9
28.
WangSQHeJH. Variational iteration method for solving integro-differential equations. Phys Lett A (2007) 367(3):188–91. 10.1016/j.physleta.2007.02.049
29.
WangSQ. A variational approach to nonlinear two-point boundary value problems. Comput Math Appl (2009) 58(11):2452–5. 10.1016/j.camwa.2009.03.050
30.
ShenYYHuangXXKwakKYangBWangS. Subcarrier-pairing-based resource optimization for OFDM wireless powered relay transmissions with time switching scheme. IEEE Trans Signal Process (2016) 65(5):1130–45. 10.1109/tsp.2016.2628351
31.
ChenQLSunZQ. The exact solution of the non-linear Schrodinger equation by the exp-function method. Therm Sci (2021) 25(3B):2057–62. 10.2298/tsci200301088c
32.
GüzelNKurulayM. Solution of shiff systems by using differential transform method. Journal of Science and Technology of Dumlupinar University (2008) 16:49–60.
33.
HeJH. Homotopy perturbation technique. Comp Methods Appl Mech Eng (1999) 178:257–62. 10.1016/s0045-7825(99)00018-3
34.
HeJHHeCHAlsolamiAA. A good initial guess for approximating nonlinear oscillators by the homotopy perturbation method. Facta Universitatis, Ser Mech Eng (2023). 10.22190/FUME230108006H
35.
AboodhKS. Application of new transform “Aboodh transform” to partial differential equations. Glob J Pure Appl Math (2014) 10(2):249–54.
36.
PekerHACuhaFA. Application of Kashuri Fundo transform and homotopy perturbation methods to fractional heat transfer and porous media equations. Therm Sci (2022) 26(4):2877–84. 10.2298/tsci2204877p
37.
AnjumNHeJH. Two modifications of the homotopy perturbation method for nonlinear oscillators. J Appl Comput Mech (2020) 2020:2482. 10.22055/JACM.2020.34850.2482
38.
NadeemMLiFQ. He-Laplace method for nonlinear vibration systems and nonlinear wave equations. J Low Frequency Noise, Vibration Active Control (2019) 38(3-4):1060–74. 10.1177/1461348418818973
39.
ManimegalaiKZephania C FSBeraPKBeraPDasSKSilT. Study of strongly nonlinear oscillators using the Aboodh transform and the homotopy perturbation method. Eur Phys J Plus (2019) 134:462–71. 10.1140/epjp/i2019-12824-6
40.
JaniHPSinghTR. Aboodh transform homotopy perturbation method for solving fractional-order Newell-Whitehead-Segel equation. Math Methods Appl Sci (2022). 10.1002/mma.8886
41.
YasminH. Application of Aboodh homotopy perturbation transform method for fractional-order convection–reaction–diffusion equation within caputo and atangana–baleanu operators. Symmetry (2023) 15(2):453. 10.3390/sym15020453
42.
JaniHPSinghTR. A robust analytical method for regularized long wave equations. Iranian J Sci Technol Trans A: Sci (2022) 46(6):1667–79. 10.1007/s40995-022-01380-9
43.
XuRHYangLBQinZ. Design, manufacture, and testing of customized sterilizable respirator. J Mech Behav Biomed Mater (2022) 131:105248. 10.1016/j.jmbbm.2022.105248
44.
GhorbaniA. Beyond adomian polynomials: He polynomials. Chaos Solitons Fractals (2009) 39:1486–92. 10.1016/j.chaos.2007.06.034
45.
LuoXKNadeemMAsjadMIAbdoMS. A computational approach for the calculation of temperature distribution in casting-mould heterogeneous system with fractional order. Comput Math Methods Med (2022) 2022:1–10. 10.1155/2022/3648277
46.
LiZBHeJH. Fractional complex transform for fractional differential equations. Math Comput Appl (2010) 15(5):970–3. 10.3390/mca15050970
47.
HeJHElaganSKLiZB. Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus. Phys Lett A (2012) 376(4):257–9. 10.1016/j.physleta.2011.11.030
48.
AinQTHeJHAnjumNAliM. The fractional complex transform: A novel approach to the time-fractional schrödinger equation. Fractals (2021) 28(7):2050141. 10.1142/s0218348x20501418
49.
HeJHEl-DibYO. A tutorial introduction to the two-scale fractal calculus and its application to the fractal Zhiber-Shabat Oscillator. Fractals (2021) 29:2150268. 10.1142/s0218348x21502686
50.
HauboldHJMathaiAMSaxenaRK. Mittag-leffler functions and their applications. J Appl Math (2011) 2011:1–51. 10.1155/2011/298628
Summary
Keywords
homotopy perturbation method, Aboodh transform, He’s polynomials, fractional differential equation, two-scale fractal theory
Citation
Tao H, Anjum N and Yang Y-J (2023) The Aboodh transformation-based homotopy perturbation method: new hope for fractional calculus. Front. Phys. 11:1168795. doi: 10.3389/fphy.2023.1168795
Received
18 February 2023
Accepted
03 April 2023
Published
27 April 2023
Volume
11 - 2023
Edited by
Hamid M. Sedighi, Shahid Chamran University of Ahvaz, Iran
Reviewed by
Szabolcs Fischer, Széchenyi István University, Hungary
Santanu Saha Ray, National Institute of Technology Rourkela, India
Magaji Adamu Yunbunga, Abubakar Tafawa Balewa University, Nigeria
Updates
Copyright
© 2023 Tao, Anjum and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Naveed Anjum, xsnaveed@yahoo.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.