Three-Dimensional Fourth-Order Time-Fractional Parabolic Partial Differential Equations and Their Analytical Solution

Mussa, Yesuf Obsie; Gizaw, Ademe Kebede; Negassa, Ayana Deressa

doi:https://doi.org/10.1155/2021/5108202

Mathematical Problems in Engineering

On this page

Abstract Introduction Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Advances in Delay Differential Equations and Applications

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5108202 | https://doi.org/10.1155/2021/5108202

Three-Dimensional Fourth-Order Time-Fractional Parabolic Partial Differential Equations and Their Analytical Solution

Yesuf Obsie Mussa,¹Ademe Kebede Gizaw,¹and Ayana Deressa Negassa¹

Academic Editor: Chuangxia Huang

Received01 Aug 2020

Revised17 Dec 2020

Accepted09 Jan 2021

Published21 Jan 2021

Abstract

In this study, the fractional reduced differential transform method (FRDTM) is employed to solve three-dimensional fourth-order time-fractional parabolic partial differential equations with variable coefficients. The fractional derivative used in this study is in the Caputo sense. A few important lemmas which are essential to solve the problems using the proposed method are proved. The novelty of this method is that it uses appropriate initial conditions and finds the solution to the problems without any discretization, linearization, perturbation, or any restrictive assumptions. Two numerical examples are considered in order to validate the efficiency and reliability of the method. Furthermore, the FRDTM solution when α = 1 is compared with other analytical methods available in the existing literature. Computational results are shown in tables and graphs. The obtained results revealed that the method is capable and simple to solve fractional partial differential equations. The software used for the calculations in this study is Mathematica 7.

1. Introduction

In the last decade, fractional calculus has become an important mathematical tool in several branches of sciences and engineering due to its enormous numbers of applications ([1] and the references therein). What specifics fractional calculus is the existence of various approaches to define differentiation with the nonintegral order [2]. As a result, several definitions for fractional derivatives appear in the literature to present more accurate models for real life phenomena. Some of known fractional derivatives are Riemann–Liouville, modified Riemann–Liouville, Caputo, Hadmard, Erdélyi–Kober, Riesz, Grunwald–Letnikov, Marchaud, and others ([3] and the references therein). All known fractional derivatives satisfy one of the well-known properties of classical derivative, namely, the linear property. However, the other properties of classical derivative, such as the derivative of a constant is zero; the product rule, quotient rule, and the chain rule either do not hold or are too complicated for many fractional derivatives.

In recent years, fractional differential equations have received considerable attention owing to their applicability in different fields of sciences such as chemistry, biology, diffusion, control theory, rheology, viscoelasticity, and so on [4]. Consequently, the solution of FPDEs represents nowadays a vigorous research area for scientists, and finding approximate and exact solutions to FPDEs is an important task [5]. However, PDEs are commonly hard to tackle, and their fractional order types are more complicated [2, 6–10]. Therefore, several methods such as the homotopy perturbation method [11], sub-ODE method [12, 13], generalized tanh method [14], residual power series method [14, 15], and so on [16–21] are developed to obtain solutions of some nonlinear fractional differential equations. Most of these methods sometimes require complex and huge calculation in order to obtain approximate solutions. To overcome such difficulties and drawbacks, an alternative method, the so called the fractional reduced differential transform method (FRDTM), has been developed by Keskin and Oturanc [22]. FRDTM plays a vital role among all the listed methods because it takes small size computation, easy to implement as compared to other techniques [23]. Using this method, it is possible to find both exact and approximate solutions in a rapidly convergent power series form. FRDTM is a very reliable, efficient, and effective powerful computational technique for solving physical problems [24–26].

In the field of modern science and engineering, the fourth-order initial/boundary value problems are of great importance. For example, airplane wings, bridge slabs, floor systems, and window glasses are being modeled as plates subjected to different types of end supports which are successfully described in terms of fourth-order PDEs [27, 28]. Time-fractional NLPDEs come from classical NLPDEs by replacing its time derivative with a fractional derivative. Time-fractional partial differential equations have many applications in areas such as electrochemistry, diffusion processes, electromagnetic, material science, chaotic dynamics, and so forth [27, 29].

In this work, we consider the class of three-dimensional fourth-order time-fractional parabolic PDEs of the form

Subjected to the initial conditions,

And boundary conditionswhere , , , and are positive, and the functions for .

The remaining sections of this study are organized as follows: in Section 2, we give some fundamental definitions and lemmas associated with fractional calculus. In Section 3, some basic definitions and properties related to the three-dimensional fractional reduced differential transform method taking the fractional derivatives in the Caputo sense are presented, and some lemmas are proved. In Section 4, we present the formulation of the method. Section 5 is devoted to apply the method to solve three-dimensional time-fractional parabolic partial differential equation with variable coefficients. In Section 6, we present the numerical results and discussion. Tables and graphs are used to show the effectiveness, validity, and performance of the FRDTM for some values of . Finally, the conclusion is presented in Section 7.

2. Fractional Calculus

In this section, some basic definitions and lemmas associated with fractional calculus are presented. Some of these definitions are due to Riemann–Liouville and Caputo sense; for details, see [2, 24, 30–33].

Definition 1. The gamma function. The gamma function is simply a generalization of the factorial real arguments. The Gamma function can be defined as

Definition 2. The fractional derivative of in Caputo sense is defined asThe basic properties of the Caputo fractional derivative can be given by the following lemma.

Lemma 1. If and , then

Definition 3. For to be the smallest integer that exceeds , and m to be the smallest integer that exceeds, β the Caputo time-fractional derivative of order is defined asand the space fractional derivative operator with respect to x, y, and z of order is, respectively, defined as

3. Fractional Reduced Differential Transform Method (FRDTM)

In this section, the basic definition and properties related to the three-dimensional fractional reduced differential transform method taking the fractional derivatives in the Caputo sense are presented.

Definition 4. (see [33, 34]). If is analytic and continuously differentiable with respect to space variables x, y, and time variable t in the domain of interest, then the spectrum function is defined as

Definition 5. (see [33, 34]). The inverse FRDTM of is defined asThe inverse transform of the set of values of gives the n-terms approximate solution as follows:Therefore, the exact solution of the problem is given by

Lemma 2. If , then .

Proof. By Definition 3 and properties of FRDTM given in Table 1, we have

Lemma 3. If , then . , and are the constants.

Proof. By Definition 3 and properties of FRDTM given in Table 1, we have_._, and are the constants.

Lemma 4. If , then . , and are the constants.
The proof is similar to Lemma 3.
If , then Lemma 4 reduces to , where .

4. Solution of the Problem by FRDTM

The steps to find the general solution of the governing equation (1) subjected to initial conditions (2) by the proposed method are as follows:

Applying the properties of FRDTM given in Table 1 to equations (1) and (2), we obtain

Solving for from equation (15), we get

When k = 0,1,2,3,... by iterative calculations, we obtain the following successive terms:and so on.

Thus, the inverse transform of the set of values of giving the n-terms approximate solution and the FRDTM solution of equation (1) subjected to the initial condition (2) is given by

5. Illustrative Examples

Example 1. Consider the following partial differential equation in three space variablessubjected to the initial conditionsWhen , equation (20) reduces to the classical (or nonfractional) three-dimensional homogeneous parabolic partial differential equation [35–37].
Applying the properties of FRDTM given in Table 1 and Lemma 2 to equations (20) and (21), we getWhen k = 0,1,2,…, we haveand so on. Continuing in this way, the remaining steps of the FRDTM can be obtained. Then, by equation (19), we obtain the following approximate solution:The exact solution of the classical form of Example 1 is as in [35–37].

Example 2. Consider the following three-dimensional nonhomogeneous fourth-order parabolic partial differential equation:subjected to initial conditionsWhen , equation (25) reduces to the classical (or nonfractional) three-dimensional nonhomogeneous parabolic partial differential equation [37].
Applying properties of FRDTM given in Table 1 and Lemma 4 to equations (25) and (26), we getSolving for from equation (27), we haveNow, when k = 0, 1,3,…,and so on. Continuing in this way, the remaining steps of the FRDTM can be obtained. Then, by equation (19), we obtainThe exact solution of the classical form of Example 2 is as in [37].

6. Numerical Results and Discussion

In this section, the numerical values of the function u (x, y, z, t) of Examples 1 and 2 are computed for various values of the fractional order α, and the exact solutions for the variables x, y, z, and t are presented in tables and figures.

Tables 2 and 3 exhibit the behavior of FRDTM solutions for different values of fractional order α and the absolute errors (E) approximated in the 5^th and 6^th order at different values of the independent variables x, y, z, and time t of equations (20) and (25), respectively. As it can be seen from these tables, when the value of the fractional order α becomes close to 1, the obtained FRDTM solutions of Examples 1 and 2 are close to the exact solutions of its respective classical (nonfractional) equations obtained by different methods. In particular, when α = 1, the FRDTM solutions of the mentioned examples are in excellent agreement with the exact solutions of their corresponding classical (nonfractional) form, for details, see [35–37].

Figures 1(a)–1(d) show the evaluation results of the approximate analytical solution of equation (20). These figures also show the behavior of the approximate solution obtained by the proposed method at fixed values x = 0.6 and y = 0.7 and for different values of the fractional order α = 0.5, 0.7, 1 and the exact solution.

(a)

(b)

(c)

(d)

Figures 2(a) and 2(b) demonstrate the solutions u (x, y, z, t) verses t of equation (20) for a fixed value of x, y, and z. In Figure 2(a), the FRDTM solutions at α = 0.7, 0.8, 1 and the exact solution are compared for a fixed value x = y = z = 0.4. In Figure 2(b), the FRDTM solutions at α = 0.7, 0.8, 1 and the exact solution are compared for a fixed value x = y = z = 0.6.

(a)

(b)

Figures 3(a)–3(d) show the evaluation results of the approximate analytical solution of equation (25). These figures also depict the behavior of the approximate solutions obtained by the proposed method at fixed values and and for a different time-fractional order and the exact solution.

(a)

(b)

(c)

(d)

Figures 4(a) and 4(b) demonstrate the solutions verses of equation (25) for a fixed value of . In Figure 4(a), the FRDTM solutions at and exact solution are compared for a fixed value . In Figure 4(b), the FRDTM solutions at and exact solution are compared for a fixed value .

(a)

(b)

7. Conclusion

In this study, we successfully employed the FRDTM to find the approximate and exact solutions of three-dimensional fourth-order time-fractional parabolic partial differential equations with variable coefficients taking the fractional derivative in Caputo sense. It is found that the results we obtained in Examples 1 and 2 when are in excellent agreement with the exact solutions obtained by VIM [35], ADM [36], and RDTM [37]. Also, the illustrated examples demonstrate that the suggested technique can be observed as a fast, simple, and efficient tool for computing three-dimensional fourth-order time-fractional parabolic partial differential equations (PDEs). The main advantage of the method is that it can be applied to fractional PDE’s and obtain approximate and exact solutions of fractional PDEs without requiring any discretization, perturbation, or restrictive conditions. Moreover, the tables and graphs constructed in this work demonstrate the physical behavior of the FRDTM solutions of the problems for different values of . As a result, we conclude that the FRDTM is very powerful, straightforward, and effective to obtain approximate and exact solutions of a wide variety problems related to fractional PDEs applied in sciences such as mathematics, physics, and engineering.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

YOM and Ademe KG proposed the main idea of this study. YOM, Ademe KG, and Ayana DN equally contributed in making the first draft and revision and approved the final manuscript for submission.

Acknowledgments

The authors would like to thank Jimma University, College of Natural Sciences, and Department of Mathematics for providing the necessary resources during conducting this study.

References

Y. Chatibi, E. H. El Kinani, and A. Ouhadan, “Variational calculus involving nonlocal fractional derivative with Mittag-Leffler kernel,” Chaos, Solitons & Fractals, vol. 118, pp. 117–121, 2019.
View at: Publisher Site | Google Scholar
I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, CA, USA, 1999.
B. A. Tayyan and A. H. Sakka, “Lie symmetry analysis of some conformable fractional partial differential equations,” Arabian Journal of Mathematics, vol. 9, pp. 201–212, 2020.
View at: Publisher Site | Google Scholar
Y. Chatibi, E. H. El Kinani, and A. Ouhadan, “Lie symmetry analysis of conformable differential equations,” AIMS Mathematics, vol. 4, no. 4, pp. 1133–1144, 2019.
View at: Publisher Site | Google Scholar
M. Bishehniasar, S. Salahshour, A. Ahmadian, F. Ismail, and D. Baleanu, “An Accurate Approximate- Analytical Technique for Solving Time-Fractional Partial Differential equations,” Complexity, vol. 2017, Article ID 8718209, 12 pages, 2017.
View at: Publisher Site | Google Scholar
D. Baleanu, J. A. T. Machado, and A. C. J. Luo, Fractional Dynamics and Control, Springer Science & Business Media, New York, NY, USA, 2011.
S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009.
View at: Publisher Site | Google Scholar
A. Elsaid, “The variational iteration method for solving Riesz fractional partial differential equations,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 1940–1947, 2010.
View at: Publisher Site | Google Scholar
Y. Chatibi, E. H. El Kinani, and A. Ouhadan, “On the discrete symmetry analysis of some classical andfractional differential equations,” Mathematical Methods in the Applied Sciences, pp. 1–11, 2019.
View at: Publisher Site | Google Scholar
Y. Chatibi, E. H. El Kinani, and A. Ouhadan, “Lie symmetry analysis and conservation laws for the time fractional Black-Scholes equation,” International Journal of Geometric Methods in Modern Physics, vol. 17, no. 1, Article ID 2050010, 2020.
View at: Publisher Site | Google Scholar
J. H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, pp. 73–79, 2003.
View at: Publisher Site | Google Scholar
F. Tchier, M. Inc, and A. Yusuf, “Symmetry analysis, exact solutions and numerical approximations for the space-time Carleman equation in nonlinear dynamical systems,” European Physical Journal Plus, vol. 134, pp. 1–18, 2019.
View at: Publisher Site | Google Scholar
A. Yusuf, M. Inc, and M. Bayram, “Soliton solutions for Kudryashov-Sinelshchikov equation,” Sigma Journal of Engineering and Natural Sciences, vol. 37, pp. 439–444, 2019.
View at: Google Scholar
M. Inc, A. Yusuf, A. I. Aliyu et al., “Dark and singular optical Solitons for the conformable space time nonlinear Schrodinger equation with Kerr and power law nonlinearity,” Optik, vol. 162, pp. 65–75, 2018.
View at: Publisher Site | Google Scholar
M. Inc, A. Yusuf, A. I. Aliyu et al., “Lie symmetry analysis and explicit solutions for the time generalized Burgers-Huxley equation,” Optical and Quantum Electronics, vol. 50, pp. 1–16, 2018.
View at: Publisher Site | Google Scholar
A. Akgul, “Reproducing kernel method for fractional derivative with non-local and non-singular kernel,” in Fractional Derivatives with Mittag-Leffler Kernel, Series of Studies in Systems, Decision and Control, Springer, Cham, Switzerland, 2019.
View at: Google Scholar
A. Akgül, “A novel method for a fractional derivative with non-local and non-singular kernel,” Chaos, Solitons & Fractals, vol. 114, pp. 478–482, 2018.
View at: Publisher Site | Google Scholar
M. G. Sakar, O. Saldir, and A. Akgul, “A novel technique for fractional Bagley-Torvik equation,” Proceedings of the National Academy of Sciences, India, Section A: Physical Sciences, vol. 89, no. 2, pp. 1–7, 2018.
View at: Publisher Site | Google Scholar
M. Modanli and A. Akgül, “Numerical solution of fractional telegraph differential equations by theta-method,” The European Physical Journal Special Topics, vol. 226, no. 16-18, pp. 3693–3703, 2017.
View at: Publisher Site | Google Scholar
E. K. Akgul, “Solutions of the linear and nonlinear differential equations within the generalized fractional derivatives,” Chaos, vol. 29, no. 2, Article ID 023108, 2019.
View at: Publisher Site | Google Scholar
E. K. Akgül, B. Orcan, and A. Akgül, “Solving higher-order fractional differential equations by reproducing kernel Hilbert space method,” Journal of Advanced Physics, vol. 7, no. 1, pp. 98–102, 2018.
View at: Publisher Site | Google Scholar
Y. Keskin and G. Oturanc, “Reduced differential transform method: a new approach to fractional partial differential equations,” Nonlinear Science Letters A, vol. 1, pp. 61–72, 2010.
View at: Google Scholar
V. K. Srivastava, M. K. Awasthi, and M. Tamsir, “RDTM solution of Caputo time fractional-order hyperbolic Telegraph equation,” AIP Advances, vol. 3, Article ID 032142, 2013.
View at: Publisher Site | Google Scholar
V. K. Srivastava, M. K. Awasthi, and S. Kumar, “Analytical approximations of two and three dimensional time-fractional telegraphic equation by reduced differential transform method,” Egyptian Journal of Basic and Applied Sciences, vol. 1, no. 1, pp. 60–66, 2014.
View at: Publisher Site | Google Scholar
V. K. Srivastava, S. Kumar, M. K. Awasthi, and B. K. Singh, “Two-dimensional time fractional-order biological population model and its analytical solution,” Egyptian Journal of Basic and Applied Sciences, vol. 1, no. 1, pp. 71–76, 2014.
View at: Publisher Site | Google Scholar
V. K. Srivastava, N. Mishra, S. Kumar, B. K. Singh, and M. K. Awasthi, “Reduced differential transform method for solving (1 + n)-dimensional Burgers' equation,” Egyptian Journal of Basic and Applied Sciences, vol. 1, no. 2, pp. 115–119, 2014.
View at: Publisher Site | Google Scholar
D.M. Yadeta, A. K. Gizaw, and Y. O. Mussa, “Approximate analytical solution of one-dimensional Beam equations by using time-fractional reduced differential transform method,” Journal of Applied Mathematics, vol. 2020, Article ID 7627385, 13 pages, 2020.
View at: Publisher Site | Google Scholar
M. Amin, M. Abbas, M. K. Iqbal, and D. Baleanu, “Non-polynomial quintic spline for numerical solution of fourth-order time fractional partial differential equations,” Advances in Difference Equations, vol. 183, 2019.
View at: Publisher Site | Google Scholar
H. Tariq and G. Akram, “Quintic spline technique for time fractional fourth-order partial differential equation,” Numerical Methods for Partial Differential Equations, vol. 33, no. 2, pp. 445–466, 2017.
View at: Publisher Site | Google Scholar
B. K. Singh, “Fractional reduced differential transform method for numerical computation of a system of linear and nonlinear fractional partial differential equations,” International Journal of Open Problems in Computer Science and Mathematics, vol. 9, no. 3, pp. 20–38, 2016.
View at: Publisher Site | Google Scholar
M. S. Rawashdeh, “A reliable method for the space-time fractional Burgers and time-fractional CahnAllen equations via the FRDTM,” Advances in Difference Equations, vol. 99, p. 14, 2017.
View at: Publisher Site | Google Scholar
S. Abuasad, I. Hashim, and S. A. Abdul Karim, “Modified fractional reduced differential transform method for the solution of multiterm time-fractional diffusion equations,” Advances in Mathematical Physics, vol. 2019, Article ID 5703916, 14 pages, 2019.
View at: Publisher Site | Google Scholar
D. Lu, J. Wang, M. Arshad, Abdullah, and A. Ali, “Fractional reduced differential transform method for space-time fractional order Heat-like and wave-like partial differential equations,” Journal of Advanced Physics, vol. 6, no. 4, pp. 598–607, 2018.
View at: Publisher Site | Google Scholar
M. Arshad, D. Lu, and J. Wang, “(N + 1)-dimensional fractional reduced differential transform method for fractional order partial differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 48, pp. 509–519, 2017.
View at: Publisher Site | Google Scholar
J. Biazar and H. Ghazvini, “He’s variational iteration method for fourth-order parabolic equations,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 1047–1054, 2007.
View at: Publisher Site | Google Scholar
A.-M. Wazwaz, “Exact solutions for variable coefficients fourth-order parabolic partial differential equations in higher-dimensional spaces,” Applied Mathematics and Computation, vol. 130, no. 2-3, pp. 415–424, 2002.
View at: Publisher Site | Google Scholar
B. Ibis, “Application of reduced differential transformation method for solving fourth-order parabolic partial differential equations,” Journal of Mathematics and Computer Science, vol. 12, no. 02, pp. 124–131, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yesuf Obsie Mussa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1143

Downloads

816

Citations