Abstract
We use the q-homotopy analysis Shehu transform method in this article to obtain analytical and numerical solutions to time fractional partial differential equations. We also give analytical solutions to two problems, as well as a comparison study in terms of absolute error with homotopy perturbation transform method, homotopy analysis transform method, and residual power series method to verify the suggested technique’s effectiveness and correctness. The numerical and graphical solutions achieved by the proposed method show that it is computationally accurate and may be used to obtain and investigate solutions to time fractional partial differential equations.
1. Introduction
Fractional calculus (FC), including integration and differentiation of arbitrary non-integer order, is the generalization of classical integration and differentiation (Oldham & Spanier, Citation1974). The beauty of FC is that fractional order derivatives and integrals are non-local. The purpose of using fractional models in differential equations in physical models (Ali, Osman, Baskonus, Elazabb, & İlhan, Citation2020; Cuahutenango-Barro et al., Citation2021; Khalid, Rehan, Nisar, & Osman, Citation2021) due to their non-local property. Recently, numerous models in biology (Dave, Khan, Purohit, & Suthar, Citation2021; Kumar, Kumar, Osman, & Samet, Citation2021; Mistry, Khan, & Suthar, Citation2020; Nisar et al., Citation2022), and dynamics (Bayones et al., Citation2021) can be modelled with the help of fractional order derivatives. In the past few years, many researchers have been established and applied numerous schemes to gain the solutions of fractional order differential equations such as (Akinyemi, Şenol, & Osman, Citation2022; Ramani et al., 2022, Ramani, Khan, & Suthar, Citation2019), HATM and RPSM (Wang, Wu, Ren, & Chen, Citation2019), reproducing kernel discretization method (Arqub, Osman, Abdel-Aty, Mohamed, & Momani, Citation2020), Chebyshev wavelet collocation method (Dhawan, Machado, Brzeziński, & Osman, Citation2021), Chebyshev collocation method (Ali, Abd El Salam, et al., Citation2020; Habenom & Suthar, Citation2020), Tikhonov regularization method (Djennadi, Shawagfeh, & Abu Arqub, Citation2021; Djennadi, Shawagfeh, Osman, Gómez-Aguilar, & Arqub, Citation2021), the Adams-Bashforth predictor-corrector method (Kumar, Chauhan, Osman, & Mohiuddine, Citation2021), Haar wavelet method (Amin, Mahariq, Shah, Awais, & Elsayed, Citation2021), reproducing kernel method (AL-Smadi, Arqub, & El-Ajou, Citation2014; Yildirim, Akgül, & Inc, Citation2021), Homotopy perturbation method (He, Citation1999; Pareek, Gupta, Agarwal, & Suthar, Citation2021).
In the present paper, we consider the initial value autonomous system of time fractional partial differential equations (TFPDEs) with proportional delay defined in (Sakar, Uludag, & Erdogan, Citation2016) (1) (1) where with represents the initial value, and denotes the differential operator. presents the Caputo derivative of
El-Tawil and Huseen (Citation2012, Citation2013) proposed a method namely q-homotopy analysis method (q-HAM) which is a more general method of HAM. The essential idea of this method is to introduce a homotopy parameter, say q, which varies from 0 to and a nonzero auxiliary parameter At q = 0, the system of equations usually has been reduced to a simplified form which normally admits a rather simple solution. As q gradually increases continuously towards the system goes through a sequence of deformations, and the solution at each stage is close to that at the previous stage of the deformation. Eventually at the system takes the original form of the equation and the final stage of the deformation gives the desired solution.
The Shehu transform (ST) is a generalization of the Laplace and the Sumudu integral transform (Watugala, Citation1998). Besides, the proposed integral transform is similar to natural transform (Khan & Khan, Citation2008). The ST become Laplace’s transform (Spiegel, Citation1965), when the variable μ = 1, and becomes the Yang’s integral transform (Yang, Citation2016) when the variable s = 1.
In this study, we used q-homotopy analysis Shehu transform method (q-HASTM) to gain the analytical solution of system (1). The proposed scheme, namely q-HASTM, is an elegant amalgamation of q-HAM and ST. Its superiority is its ability to adjust two strong computational methodologies for probing FDEs. By choosing proper we can control the convergence region of solution series in a large permissible domain. The advantage of q-HASTM in that it does not require linearization or discretization, shows little perturbations, has no restrictive assumptions, lessens mathematical computations significantly, offers non-local effect, promises a big convergence region and is free from obtaining difficult polynomials, integrations and physical parameters.
2. Definitions
Definition 2.1.
The fractional Caputo derivative of function is defined as (Diethelm & Ford, Citation2002) (2) (2)
In particular, if then
Lemma 2.1.
The Caputo derivative satisfies some following properties for
p = 0, for a constant;
Definition 2.2.
The ST (Maitama & Zhao, Citation2019b) of the function of exponential order is defined over the set of functions, by the following integral (3) (3)
It converges if the limit of the integral exists, and diverges if not.
Definition 2.3.
The ST of Caputo fractional derivative operator was introduced by Maitama and Zhao (Citation2019b) and expressed as (4) (4) here is the ST of
3. Basic idea of proposed technique
Consider a time fractional partial differential system (5) (5) where presents the Caputo derivative of and denotes the differential operator.
Taking the ST to both sides of Equation (5) and on simplifying, we get (6) (6)
Now, we define a non-linear operator as (7) (7)
In Equation (7), is an embedding parameter and is indicate the real function of and q. Liao (Citation1992, Citation1995) constructed zeroth-order deformation equation such as (8) (8) here S represents the ST, is nonzero auxiliary parameter, denoted an auxiliary function, and expresses the initial gauss of and is unknown function. Let q = 0 and q = 1 in Equation (8), we get (9) (9)
Thus, if q rises from 0 to the series solution varies from the initial guess to the solution Upon expanding with the help of Taylor’s series near to q, we have (10) (10) where (11) (11)
By proper choosing of and the series in Equation (10) converges at q = we will get (12) (12)
We define the vector as (13) (13)
First, differentiating Equation (8) p-times with respect to q, then evaluate at q = 0 and finally dividing by we have the so-called pth-order deformation equation (14) (14) where (15) (15) and (16) (16)
Taking the inverse ST to both sides of Equation (14) and with the aid of Equations (8) and (15), we get (17) (17)
Based on Equation (5), is define as (18) (18)
Finally, we compute by using Equation (17) for Hence the Mth order approximate solution of Equation (5) can be represented as (19) (19)
Moreover, for we get (20) (20)
The existence of the factor in the q-HASTM solution (20) allows for faster convergence than the standard HAM. Moreover, in the special case n = 1, the q-HASTM reduces to the standard homotopy analysis Shehu transform method (HASTM).
4. Convergence analysis
In this section, we investigate the convergence analysis of q-HASTM technique.
Theorem 4.1.
Let satisfy the Lipschitz condition with the Lipschitz constant δ. The solution derived with the aid of q-HASTM of the time fractional partial differential system (5) is unique, wherever , where
Proof.
The solution of the time fractional partial differential system (5) is presented as (21) (21) where (22) (22)
Now, let ν and be two different solutions of considered time fractional partial differential system, then we have (23) (23)
With the aid of the convolution theorem, we can obtain (24) (24)
Next, putting up the integral mean value theorem in use, it yields (25) (25)
It gives Because therefore, which implies that Hence the solution is unique.
Theorem 4.2
(Convergence theorem). Let us consider that X be a Banach space and there is a nonlinear mapping and assume that (26) (26)
Then in view of Banach’s fixed point theory, W has a fixed point. Furthermore, the sequence generated by the q-HASTM with an arbitrary selection of converges to the fixed point of W and (27) (27)
Proof.
Let us take a Banach space of all continuous functions on I with the norm expressed as
Now, we show that the sequence is a Cauchy sequence in the Banach space.
Now, making use of the convolution theorem for the ST, it gives (28) (28)
Next, by the application of the integral mean value theorem (Maitama & Zhao, Citation2019a), we obtain
Let m = n + 1, then we have (29) (29)
On using the triangular inequality, it yields
Because so then we have (30) (30)
But so as then Therefore, the sequence is Cauchy sequence in and so the sequence is convergent.
5. Numerical problem
In this section, we consider two numerical problems to prove the accuracy, and efficiency of our proposed method. All the numerical and graphical results for the following two problems are calculated by utilizing the software scilab-6.0.2.
Problem 1.
Consider the time-fractional generalized Burger’s equation (Sakar et al., Citation2016) as (31) (31) where presents the Caputo derivative of and subject to initial condition (32) (32)
By performing the ST on both sides of Equation (31) and with the help of Equation (32), we get (33) (33)
According to proposed scheme, we define the non-linear operator as (34) (34)
Form Equation (14) and choosing the pth order deformation equation is given as (35) (35) and by using of Equations (15) and (34) the valve of is given as (36) (36)
Operating the inverse ST to Equation (35) and by using Equation (36), we get (37) (37)
By putting p = 1, 2, 3 in Equation (37) and with Equation (32), we obtain the following results (38) (38) (39) (39)
In the same manner, we can get (40) (40) (41) (41)
Hence, the fourth order approximate solution of (31) is given as (42) (42)
In particular, when we take and then the solution converge to the exact solution of (31) very fastly (43) (43)
Problem 2.
Consider the TFPDEs as given in Sakar et al. (Citation2016) and Singh and Kumar (Citation2018) with proportional delay (44) (44) subject to initial condition (45) (45)
Applying ST to both sides of Equation (44) and on, we get (46) (46)
According to proposed technique, the nonlinear operator decomposed as following (47) (47)
Form Equation (14) and choosing the order deformation equation is given as (48) (48) and by using of Equations (15) and (47) the valve of is given as (49) (49)
Operating the inverse ST to Equation (48) and by using Equation (49), we get (50) (50)
By putting p = 1, 2, 3 in Equation (50) and with Equation (45), we obtain the following results (51) (51) (52) (52)
In the same manner, we can obtain (53) (53) (54) (54)
Similarly, we obtain next terms in the same manner. Hence, we get the fourth order approximate solution of (44) as follow (55) (55)
In particular, when we take and then the solution converge to the exact solution of (44) very fastly (56) (56)
6. Results and discussions
In the present study, we can analyzed from and that the numerical results gained by the proposed technique almost same to the schemes presented in the literature (Singh & Kumar, Citation2018; Wang et al., Citation2019). depicts the comparative analysis between approximate solution for distinct ζ and the exact solution. We also present the absolute error for Problem 1 at n = 1, represents the nature of obtained solution for (31) with the distinct ζ and compression between approximate solution and exact solution (43). presents the behaviour of approximate solution of (31) with n at and explores the -curves of the q-HASTM solution (31) with distinct values of n at and which helps to adjust the convergence region.
Moreover, represents the behaviour of approximate solution at distinct ζ in comparison with the exact solution and the absolute error for Problem 2. presents the behaviour of approximate solution with various ζ and exact solution of (44). explores nature of the q-HASTM solution with distinct ζ. explores the -curves of the q-HASTM solution (44) with distinct values of n at and which helps to adjust the convergence region.
7. Conclusions
In this article, we successfully implemented q-HASTM to find the analytical and numerical solutions of time-fractional partial differential equations (TFPDEs). We obtained the analytical and numerical solutions of two applications of TFPDEs to present the effectiveness and accuracy of proposed scheme. Moreover, the q-homotopy analysis transform method provided the convergent series solution with easily determinable components without using any perturbation, linearization or limiting assumption. If we assume and n = 1 in q-HASTM solution, then the q-HASTM solution presented an excellent agreement with the exact solution of TFPDEs. The numerical and graphical solutions obtained by q-HASTM are presented that the proposed technique is computationally very accurate and attractive technique to obtain and investigate the solutions of time fractional partial differential equations.
Authors’ contributions
The authors contributed equally and significantly in writing this paper. All authors read and approved the final manuscript.
Acknowledgements
The authors are grateful to the editor and reviewers for their thorough review and comments, which contributed to improving this article.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Akinyemi, L., Şenol, M., & Osman, M. S. (2022). Analytical and approximate solutions of nonlinear Schrödinger equation with higher dimension in the anomalous dispersion regime. Journal of Ocean Engineering and Science, 7(2), 143–154. doi:10.1016/j.joes.2021.07.006
- Ali, K. K., Abd El Salam, M. A., Mohamed, E., Samet, B., Kumar, S., & Osman, M. S. (2020). Numerical solution for generalized nonlinear fractional integro-differential equations with linear functional arguments using Chebyshev series. Advances in Difference Equations, 2020(1), 494. doi:10.1186/s13662-020-02951-z
- Ali, K. K., Osman, M. S., Baskonus, H. M., Elazabb, N. S., & İlhan, E. (2020). Analytical and numerical study of the HIV-1 infection of CD4+ T-cells conformable fractional mathematical model that causes acquired immunodeficiency syndrome with the effect of antiviral drug therapy. Mathematical Methods in the Applied Sciences, 1–17. doi:10.1002/mma.7022
- Al-Smadi, M., Arqub, O. A., & El-Ajou, A. (2014). A numerical iterative method for solving systems of first-order periodic boundary value problems. Journal of Applied Mathematics, 2014, 135465. doi:10.1155/2014/135465
- Amin, R., Mahariq, I., Shah, K., Awais, M., & Elsayed, F. (2021). Numerical solution of the second order linear and nonlinear integrodifferential equations using Haar wavelet method. Arab Journal of Basic and Applied Sciences, 28(1), 12–19. doi:10.1080/25765299.2020.1863561
- Arqub, O. A., Osman, M. S., Abdel-Aty, A.-H., Mohamed, A.-B. A., & Momani, S. (2020). A numerical algorithm for the solutions of ABC singular Lane–Emden type models arising in astrophysics using reproducing kernel discretization method. Mathematics, 8(6), 923. doi:10.3390/math8060923
- Bayones, F. S., Nisar, K. S., Khan, K. A., Raza, N., Hussien, N. S., Osman, M. S., & Abualnaja, K. M. (2021). Magneto-hydrodynamics (MHD) flow analysis with mixed convection moves through a stretching surface. AIP Advances, 11(4), 045001. doi:10.1063/5.0047213
- Cuahutenango-Barro, B., Taneco-Hernández, M. A., Lv, Y.-P., Gómez-Aguilar, J. F., Osman, M. S., Jahanshahi, H., & Aly, A. A. (2021). Analytical solutions of fractional wave equation with memory effect using the fractional derivative with exponential kernel. Results in Physics, 25, 104148. doi:10.1016/j.rinp.2021.104148
- Dave, S., Khan, A. M., Purohit, S. D., & Suthar, D. L. (2021). Application of green synthesized metal nano-particles in photo catalytic degradation of dyes and its mathematical modelling using Caputo Fabrizio fractional derivative without singular kernel. Journal of Mathematics, 2021, 9948422. doi:10.1155/2021/9948422
- Dhawan, S., Machado, J., Brzeziński, D. W., & Osman, M. S. (2021). A Chebyshev wavelet collocation method for some types of differential problems. Symmetry, 13(4), 536. doi:10.3390/sym13040536
- Diethelm, K., & Ford, N. J. (2002). Analysis of fractional differential equations. Journal of Mathematical Analysis and Applications, 265(2), 229–248. doi:10.1006/jmaa.2000.7194
- Djennadi, S., Shawagfeh, N., & Abu Arqub, O. (2021). A fractional Tikhonov regularization method for an inverse backward and source problems in the time-space fractional diffusion equations. Chaos, Solitons and Fractals, 150, 111127. doi:10.1016/j.chaos.2021.111127
- Djennadi, S., Shawagfeh, N., Osman, M. S., Gómez-Aguilar, J. F., & Arqub, O. A. (2021). The Tikhonov regularization method for the inverse source problem of time fractional heat equation in the view of ABC-fractional technique. Physica Scripta, 96(9), 094006. doi:10.1088/1402-4896/ac0867
- El-Tawil, M. A., & Huseen, S. N. (2012). The q-homotopy analysis method. International Journal of Applied Mathematics and Mechanics, 8, 51–75.
- El-Tawil, M. A., & Huseen, S. N. (2013). On convergence of the q-homotopy analysis method. International Journal of Contemporary Mathematical Sciences, 8, 481–497. doi:10.12988/ijcms.2013.13048
- Habenom, H., & Suthar, D. L. (2020). Numerical solution for the time-fractional Fokker–Planck equation via shifted Chebyshev polynomials of the fourth kind. Advances in Difference Equation, 315, 2020. doi:10.1186/s13662-020-02779-7
- He, J. (1999). Homotopy perturbation technique. Computer Methods in Applied Mechanics and Engineering, 178(3–4), 257–262 doi:10.1016/S0045-7825(99)00018-3
- Khalid, A., Rehan, A., Nisar, K. S., & Osman, M. S. (2021). Splines solutions of boundary value problems that arises in sculpturing electrical process of motors with two rotating mechanism circuit. Physica Scripta, 96(10), 104001. doi:10.1088/1402-4896/ac0bd0
- Khan, Z. H., & Khan, W. A. (2008). N-transform-properties and applications. NUST Journal of Engineering Sciences, 1(1), 127–133.
- Kumar, S., Chauhan, R. P., Osman, M. S., & Mohiuddine, S. A. (2021). A study on fractional HIV-AIDs transmission model with awareness effect. Mathematical Methods in the Applied Sciences, 1–15. doi:10.1002/mma.7838
- Kumar, S., Kumar, R., Osman, M. S., & Samet, B. (2021). A wavelet based numerical scheme for fractional order SEIR epidemic of measles by using Genocchi polynomials. Numerical Methods for Partial Differential Equations, 37(2), 1250–1268. doi:10.1002/num.22577
- Liao, S. J. (1992). The proposed homotopy analysis technique for the solution of nonlinear problems (PhD thesis). Shanghai Jiao Tong University, Shanghai.
- Liao, S. J. (1995). An approximate solution technique not depending on small parameters: A special example. International Journal of Non-Linear Mechanics, 30(3), 371–380. doi:10.1016/0020-7462(94)00054-E
- Maitama, S., & Zhao, W. (2019a). New Laplace-type integral transform for solving steady heat-transfer problem. Thermal Science, 25(1 Part A), 1–12. doi:10.2298/TSCI180110160M
- Maitama, S., & Zhao, W. (2019b). New integral transform: Shehu transform a generalization of Sumudu and Laplace transform for solving differential equations. International Journal of Analysis and Applications, 17(2), 167–190.
- Mistry, L., Khan, A. M., & Suthar, D. L. (2020). An epidemic slia mathematical model with Caputo Fabrizio fractional derivative. Test Engineering and Management, 83, 26374–26391.
- Nisar, K. S., Ciancio, A., Ali, K. K., Osman, M. S., Cattani, C., Baleanu, D., … Azeem, M. (2022). On beta-time fractional biological population model with abundant solitary wave structures. Alexandria Engineering Journal, 61(3), 1996–2008. doi:10.1016/j.aej.2021.06.106
- Oldham, K. B., & Spanier, J. (1974). The fractional calculus theory and applications of differentiation and integration to arbitrary order. New York: Academic Press.
- Pareek, N., Gupta, A., Agarwal, G., & Suthar, D. L. (2021). Natural transform along with HPM technique for solving fractional ADE. Advances in Mathematical Physics, 2021, 9915183. doi:10.1155/2021/9915183
- Ramani, P., Khan, A. M., Suthar, D. L. , Kumar, D. ( 2022). Approximate analytical solution for non-linear Fitzhugh–Nagumo equation of time fractional order through fractional reduced differential transform method. International Journal of Applied and Computational Mathematics, 8(1), 1–12. doi:https://doi.org/10.1007/s40819-022-01254-z
- Ramani, P., Khan, A. M., & Suthar, D. L. (2019). Revisiting analytical-approximate solution of time fractional Rosenau–Hyman equation via fractional reduced differential transform method. International Journal on Emerging Technologies, 10(2), 403–409.
- Sakar, M. G., Uludag, F., & Erdogan, F. (2016). Numerical solution of time-fractional nonlinear pdes with proportional delays by homotopy perturbation method. Applied Mathematical Modelling, 40(13–14), 6639–6649 doi:10.1016/j.apm.2016.02.005
- Singh, B. K., & Kumar, P. (2018). Homotopy perturbation transform method for solving fractional partial differential equations with proportional delay. SeMA Journal, 75(1), 111–125. doi:10.1007/s40324-017-0117-1
- Spiegel, M. R. (1965). Theory and problems of Laplace transform. Schaum’s outline series. New York: McGraw-Hill.
- Wang, L., Wu, Y., Ren, Y., & Chen, X. (2019). Two analytical methods for fractional partial differential equations with proportional delay. IAENG International Journal of Applied Mathematics, 49(1), 1–6.
- Watugala, G. K. (1998). Sumudu transform a new integral transform to solve differential equations and control engineering problems. Mathematical Engineering in Industries, 6(1), 319–329.
- Yang, X. J. (2016). A new integral transform method for solving steady heat-transfer problem. Thermal Science, 20(suppl. 3), 639–642. doi:10.2298/TSCI16S3639Y
- Yildirim, E. N., Akgül, A., & Inc, M. (2021). Reproducing kernel method for the solutions of non-linear partial differential equations. Arab Journal of Basic and Applied Sciences, 28(1), 80–86. doi:10.1080/25765299.2021.1891678