- 1Department of Mathematics, Faculty of Science, Al-Balqa’ Applied University, Al-Salt, Jordan
- 2Department of Basic Engineering Sciences, College of Engineering, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
In this paper, we introduce a series solution to a class of hyperbolic system of time-fractional partial differential equations with variable coefficients. The fractional derivative has been considered by the concept of Caputo. Two expansions of matrix functions are proposed and used to create series solutions for the target problem. The first one is a fractional Laurent series, and the second is a fractional power series. A new approach, via the residual power series method and the Laplace transform, is also used to find the coefficients of the series solution. In order to test our proposed method, we discuss four interesting and important applications. Numerical results are given to authenticate the efficiency and accuracy of our method and to test the validity of our obtained results. Moreover, solution surface graphs are plotted to illustrate the effect of fractional derivative arrangement on the behavior of the solution.
1 Introduction
Many natural phenomena have been modeled through partial differential equations (PDEs), especially in physics, engineering, chemistry, and biology, as well as in humanities [1, 2]. A wide range of PDEs can be classified under the name of hyperbolic PDEs that have the following general form [2–6]:
subject to the following initial condition:
The equations of compressible fluid flow and the Euler equations are examples of PDEs that can be reduced to hyperbolic PDEs when the effects of viscosity and heat conduction are neglected [6]. In addition, many mathematical models are appearing as hyperbolic systems of PDEs that have the following general form:
subject to
where
In recent decades, many mathematical models have been reformulated using the concept of fractional calculus because they are found to reflect the phenomenon that has been modeled in a more precise and realistic way by replacing the ordinary derivative with a fractional derivative (FD) of the model. The concept of fractional calculus dates back to the 17th century [7, 8] and has recently gained considerable interest because of its extensive use in widespread fields, for instance, engineering, biological, chemical, and applied physics such as in nonlinear oscillation, waves, and diffusion as we mentioned [7–13]. In fact, from that date until now, there are many definitions of the FD. The most popular definition is the Caputo FD that is denoted and defined as [7, 8]
where
For more details about the properties of the two previous definitions, readers can refer to the references [7–12]. The most useful properties that we need in this research are
As mentioned, the definition of Caputo is one of the most important definitions of the FD, since it was and still is an appropriate and effective tool in the modeling of many natural phenomena in all sciences and fields. For example, but not limited to, the definition of Caputo has recently been used to construct a mathematical model to illustrate the impacts of deforestation on wildlife species [13], in a fractional investigation of bank data [14], to model the spread of hookworm infection [15], and newly to model and analyze the dynamics of novel coronavirus (COVID-19) [16].
It is difficult to find exact solutions (ESs) for the fractional differential and integral equations; for this reason, analytical and numerical methods are usually used to find approximate solutions (ASs) for those equations. In recent decades, many methods have been used to find analytical and numerical solutions for fractional differential and integral equations such as the variational iteration method [17], the Adomian decomposition method [18], the homotopy perturbation method [17], the homotopy analysis method, the fractional transform method [19], Green’s function method [20], and other methods [21, 22].
In the last five years, the residual power series method (RPSM) has achieved an advanced rank among the methods used to find ASs for many fractional differential and integral equations. It has been used in determining ESs and ASs for many equations such as homogeneous and non-homogeneous time- and space-fractional telegraph equation [23], time-fractional Boussinesq-type and space-fractional Klein–Gordon–type equations [24], fractional multi-pantograph system [25], space- and time-fractional linear and nonlinear KdV–Burgers equation [26], multi-energy groups of neutron diffusion equations [27], and other equations. The RPSM is characterized by its ease and speed in finding solutions for equations in the form of a power series. In fact, the RPSM is a mechanism for finding the coefficients of the fractional power series (FPS) without having to find a recurrence relation that we normally obtain from equating the corresponding coefficients in the series. The RPSM is a good alternate proceeding for gaining analytic solutions for fractional PDEs.
Despite all the features we mentioned about the RPSM, we will present in this paper a major modification to the method. We use the concept of limit at infinity instead of the concept of FD in determining the coefficients of the power series solution (SS). As is well known, finding an FD manually is not easy and sometimes takes tens of minutes when it is calculated by software packages, while calculating the limit is much easier than calculating the FD manually and faster by compute. Indeed, the RPSM determines the coefficients of the power SS of the differential or integral equations, whereas the proposed technique determines the coefficients of the expansion that represents the Laplace transform (LT) of the solution. Therefore, we do not need FDs during the transaction-finding process. To be able to implement the new method, we need to provide two expansions of the matrix functions, one to represent the solution of the equation and the other to represent the LT of the solution. Moreover, the convergence of the introduced expansions is studied. In fact, the proposed method called the Laplace-RPSM (L-RPSM) was first introduced by the authors in [28] and used for introducing exact and approximate SSs to the linear and nonlinear neutral FDEs. El-Ajou [29] then adapted the new method in creating solitary solutions for the nonlinear time-fractional partial differential equations (T-FPDEs).
Several articles are interested in providing ASs to T-FPDEs of hyperbolic type, such as the Caputo time-fractional–order hyperbolic telegraph equation [30], hyperbolic T-FPDEs [31–35], the time-fractional diffusion equation [36], fractional advection–dispersion flow equations [37], and other hyperbolic equations. However, a limited number of research studies provided analytical and numerical solutions for hyperbolic systems of T-FPDEs. Kochubei [38] presented a numerical–analytical solution for homogeneous hyperbolic fractional systems, and Hendy et al. [39] introduced a solution for two-dimensional fractional systems that was provided by a separate contrast scheme. Therefore, more research is needed in providing analytical and numerical solutions for such systems due to their importance in many applications as mentioned above.
The present work aims to apply the L-RPSM to construct ASs of a hyperbolic system of T-FPDEs with variable coefficients in the sense of Caputo’s FD, which are given in the form of the following model:
subject to
where
This paper is organized as follows: In Section 2, the analysis of matrix FPS is prepared. In Section 3, the construction of FPS solution to a hyperbolic system of T-FPDEs with variable coefficients in the sense of Caputo’s FD is presented using the L-RPSM. In Section 4, application models and graphical and numerical simulations are performed in order to illustrate the capability and the simplicity of the proposed method. Finally, the conclusion is presented in Section 5.
2 Preliminaries of Matrix FPS
Here, we present some definitions and theories regarding matrix analysis and the matrix FPS, which are playing a central role in constructing the L-RPSM solution to a hyperbolic system of T-FPDEs with variable coefficients.
Definition 2.1. The R-LFIO of order
Definition 2.2. Caputo’s time-FD operator of order
Lemma 2.1. If
1.
2.
Definition 2.3 Let
Definition 2.4 For
is called a bivariate matrix FPS around
Theorem 2.1. Assume that
Proof. Of the operator definition in Eqs 6, 13 we have
Theorem 2.2. Assume that
Proof. From Theorem 2.1, it suffices to demonstrate that
According to Lemma 2.1, it is easy to show that the formula is correct for
Thus, the proof of Theorem 2.2 has been completed. Let us call the series (Eq. 17) the bivariate fractional matrix Taylor’s formula (BFMTF) of the matrix function
Theorem 2.3. (The Remainder Estimation Theorem) Assume that
Proof. The definition of the remainder of the BFMTF of
According to Theorem 2.2, the remainder can be expressed as
So, for
Thus, the proof is completed. Note that when
which can be applied throughout this work. Finally, it is worth to mention that if
which is the bivariate classical matrix Taylor’s formula of a matrix function.
Lemma 2.2.Let
1.
2.
3.
Corollary 2.1. Let
where
Theorem 2.4. Let
Proof. As it is assumed in the text of the theorem, suppose that
As in Eq. 19, the remainder of the FME in Eq. 23 is
Multiplying Eq. 26 by
Thus, it follows that
Finally, for
Thus, we reach the end of the proof.
3 Applying the L-RPSM to the Hyperbolic System of T-FPDEs
In this section, we construct an AS to the hyperbolic system of T-FPDEs with variable coefficients given in Eqs 11, 12 by using the L-RPSM. To achieve it, firstly, we apply the LT on both sides of Eq. 11, and use Lemma 2.2, and Eq. 12; then, we have
where
where
Of course, treating with a finite series is acceptable more than an infinite series. For this reason, the L-RPSM deals with a finite series while calculating coefficients of the SS. So, we express the kth truncated series (kth TS) of
To apply the L-RPSM for determining the coefficients
and the kth residual matrix function (
The main idea of the L-RPSM can be shown in the following clear facts related to the RMF and
1.
2.
3.
where
4. Thus,
5.
where
6. Using the following fact determines the unknown coefficients
Now, to find
Multiply Eq. 38 by
Compute the limit to Eq. 39 as
Similarly, to determine the second unknown coefficient in Eq. 32,
Multiply Eq. 41 by
In general, to determine the nth unknown coefficient in Eq. 32,
This procedure can be repeated for the required number of FME coefficients representing the solution of Eq. 30. Therefore, the kth approximation of the solution of Eq. 30 can be represented as the following finite series:
If we act the inverse LT on both sides of Eq. 44, then we obtain the kth approximation of the solution of the initial value problem (IVP) (Eqs 11, 12), which takes the following expression:
4 Applications and Numerical Simulations
To test our proposed method, we present in this section four interesting and important applications. The first three applications are prepared so that the ES is already known, while the last application is prepared without knowing the solution in advance to test the predictability of the solution or obtain a suitable AS.
Application 4.1. Consider the following homogeneous hyperbolic system of T-FPDEs with variable coefficients:
subject to
where
and the ES is
To obtain an FME solution for this application using the L-RPSM, transform Eq. 46 to the Laplace space as follows:
Let the solution of Eq. 48 have a form of the FME as in Eq. 31. According to the condition in Eq. 47, the first coefficient of the FME in Eq. 31,
and the kth
To find the first unknown coefficient
Multiply Eq. 51 by
Take the limit at infinity for Eq. 52 and use Eq. 37 to get
Similarly, we can obtain the second unknown coefficient in Eq. 49,
Multiply Eq. 54 by
If we repeat the previous procedure for
So, the ES for Eq. 48 will be as follows:
If we apply the inverse LT on Eq. 57, then the SS of the IVP (Eqs 46, 47) will take the following form:
which coincides with the ES.
Application 4.2. Consider the following non-homogeneous hyperbolic system of T-FPDEs with variable coefficients:
subject to
where
According to the construction in Section 3, the LT of Eq. 59 can be represented by
The kth TS of the FME of the solution of Eq. 61 has the following form:
and the kth
So, to set the first unknown coefficient of Eq. 62, substitute
and use the result in Eq. 37 to obtain
Again, substitute
According to the fact in Eq. 37, we have
Similarly, anyone can check that
Thus, the FME solution of the IVP (Eqs 59, 60) would be as follows:
which is identical to the ES
Application 4.3. Consider the following non-homogeneous hyperbolic system of T-FPDEs with variable coefficients:
subject to
where
and the ES is
where
Mathematica 7 software has been used through a low-RAM PC for obtaining all numerical calculations and symbolism. Since the Mittag-Leffler function is an infinite expansion, it was difficult to perform the calculations using the Mittag-Leffler function as it is. For this, the fifth truncated series of the expansion in Eq. 73 was used throughout the calculations.
Like the previous applications, transform Eq. 70 to the Laplace space using the initial condition in Eq. 71 to read as follows:
Let the solution of the algebraic equation (74) has an FME as in Eq. 31. Then, the kth TS of the FME of
and the kth
Now, to determine
where
Using the same previous approach, we find the following vector coefficients of Eq. 75:
So, the fifth AS of Eq. 74 can be written as follows:
Transforming the AS in Eq. 80 to the t-space by the inverse LT, we get the fifth approximation of the solution of the IVP (Eqs 70, 71) as follows:
Obviously, there is a pattern between the terms of Eq. 81 that gives us the ES as in Eq. 72.
The mathematical behavior of the solution of the IVP (Eqs 70, 71) is illustrated next by plotting the three-dimensional space figures of the fifth approximation of the two components of the vector solution in Eq. 81 for different values of α and a fixed value of
FIGURE 1. Surface graphs of the fifth AS of
Figures 1C,D show that the fifth AS of the IVP (Eq. 70, 71) is excellent compared to the ES, as well as in the previous cases, which have not been documented in order not to increase the numbers of graphs. It is known that, by increasing the number of terms in the series, the accuracy of the solution increases and, thus, the error of solution reduces; therefore, we can reduce the error of the solution by calculating more coefficients of the FME solution as in Eq. 31.
In the next application, the ES is unknown. Therefore, we are trying to find the ES or an appropriate approximation of the solution.
Application 4.4. Consider the following non-homogeneous hyperbolic system of T-FPDEs with variable coefficients:
subject to
where
Similar to the previous applications, the LT of Eq. 82 is given by
the kth TS of the expansion of the solution of Eq. 84 is given as
and the kth
According to the fact in Eq. 37, we can create, successively, the following first eight coefficients of the expansion in Eq. 31 for this application:
Thus, the seventh AS of Eq. 84 has the following expression:
so the seventh AS of the IVP (Eqs 82, 83) can be expressed as follows:
To test the AS in Eq. 89, we need to find the norm of residual error vector (
and the Frobenius norm is chosen for error analysis and defined by
Tables 1, 2 show the values of
5 Conclusion
We have found that the ES for the hyperbolic system of T-FPDEs with variable coefficients is available if the solution is a linear combination of power functions or if it is a composite of an elementary function and a power function. In case the ES is not available, a good approximation of the solution can be obtained. The L-RPSM is an effective, accurate, easy, and speed technique in obtaining the values of coefficients for the SS. Through this work, we have presented a solution that may be missing for this kind of problem and we have opened the way for researchers to provide other ways to solve this class of equations. Moreover, the newly proposed technique can be used to construct many types of the ordinary or partial DEs of fractional order such as Lane–Emden, Boussinesq, KdV–Burgers,
Data Availability Statement
All datasets generated for this study are included in the article/supplementary files.
Author Contributions
The idea of this work, implementation, and output of this form was carried out by both authors.
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.
Acknowledgments
The first author expresses his appreciation and thanks to Al-Balqa Applied University for granting him a sabbatical leave that made this research possible.
References
1. Burger M, Caffarelli L, Markowich PA. Partial Differential Equation Models in the Socio-Economic Sciences. Phil Trans R Soc A (2014) 372:20130406. doi:10.1098/rsta.2013.0406
2. Mattheij R, Rienstra S, Boonkkamp J. Partial Differential Equations: Modeling, Analysis, Computation, Technische Universiteit Eindhoven Eindhoven. Netherlands: Siam (2005). doi:10.1137/1.9780898718270
3. Evans LC. Partial Differential Equations. American Mathematical Society (2002). doi:10.1090/gsm/019
4. Leray J. Hyperbolic Differential Equations. Princeton, NJ: Institute for Advanced Study (1953). doi:10.11948/2017095
5. Sanchez Y, Vickers J. Generalised Hyperbolicity in Spacetimes with Lipschitz Regularity. J Math Phys (2017) 58:022502. doi:10.1063/1.4975216
6. Toro EF. Notions on Hyperbolic Partial Differential Equations. In: Riemann Solvers and Numerical Methods for Fluid Dynamics. Heidelberg, Berlin: Springer (2009). doi:10.1007/b79761
7. Caputo M. Linear Models of Dissipation Whose Q Is Almost Frequency Independent--II. Geophys J Int (1967) 13:529–39. doi:10.1111/j.1365-246x.1967.tb02303.x
8. Mainardi F. Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press. London, UK: Imperial College (2010). doi:10.1142/p614
9. Kilbas AA, Srivastava HM, Trujillo JJ. Theory and Applications of Fractional Differential Equations. Netherlands: Elsevier, Amsterdam (2006).
10. Magin RL, Ingo C, Colon-Perez L, Triplett W, Mareci TH. Characterization of Anomalous Diffusion in Porous Biological Tissues Using Fractional Order Derivatives and Entropy. Microporous Mesoporous Mater (2013) 178:39–43. doi:10.1016/j.micromeso.2013.02.054
11. Tarasov VE. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles. Fields and Media. Berlin, Germany: Springer (2011).
13. Qureshi S, Yusuf A. Mathematical Modeling for the Impacts of Deforestation on Wildlife Species Using Caputo Differential Operator. Chaos, Solitons & Fractals (2019) 126:32–40. doi:10.1016/j.chaos.2019.05.037
14. Li Z, Liu Z, Khan MA. Fractional Investigation of Bank Data with Fractal-Fractional Caputo Derivative. Chaos, Solitons & Fractals (2020) 131:109528. doi:10.1016/j.chaos.2019.109528
15. Mustapha UT, Qureshi S, Yusuf A, Hincal E. Fractional Modeling for the Spread of Hookworm Infection under Caputo Operator. Chaos, Solitons & Fractals (2020) 137:109878. doi:10.1016/j.chaos.2020.109878
16. Ali A, Alshammari FS, Islam S, Khan MA, Ullah S. Modeling and Analysis of the Dynamics of Novel Coronavirus (COVID-19) with Caputo Fractional Derivative. Results Phys (2021) 20:103669. doi:10.1016/j.rinp.2020.103669
17. Momani S, Odibat Z. Comparison between the Homotopy Perturbation Method and the Variational Iteration Method for Linear Fractional Partial Differential Equations. Comput Maths Appl (2007) 54:910–9. doi:10.1016/j.camwa.2006.12.037
18. Daftardar-Gejji V, Bhalekar S. Solving Multi-Term Linear and Non-linear Diffusion-Wave Equations of Fractional Order by Adomian Decomposition Method. Appl Maths Comput (2008) 202:113–20. doi:10.1016/j.amc.2008.01.027
19. Das S, Gupta PK. Homotopy Analysis Method for Solving Fractional Hyperbolic Partial Differential Equations. Int J Comput Maths (2011) 88:578–88. doi:10.1080/00207161003631901
20. Momani S, Odibat ZM. Fractional Green Function for Linear Time-Fractional Inhomogeneous Partial Differential Equations in Fluid Mechanics. J Appl Math Comput (2007) 24:167–78. doi:10.1007/bf02832308
21. El-Ajou A. Taylor’s Expansion for Fractional Matrix Functions: Theory and Applications. J Math Comput Sci (2020) 21(2):1–17. doi:10.22436/jmcs.021.01.01
22. Oqielat Ma. N, El-Ajou A, Al-Zhour Z, Alkhasawneh R, Alrabaiah H. Series Solutions for Nonlinear Time-Fractional Schrödinger Equations: Comparisons between Conformable and Caputo Derivatives. Alexandria Eng J (2020) 59(4):2101–14. doi:10.1016/j.aej.2020.01.023
23. El-Ajou A, Oqielat Ma. N, Al-Zhour Z, Momani S. A Class of Linear Non-homogenous Higher Order Matrix Fractional Differential Equations: Analytical Solutions and New Technique. Fract Calc Appl Anal (2020) 23(2):356–77. doi:10.1515/fca-2020-0017
24. El-Ajou A, Al-Smadi M, Oqielat Ma. N, Momani S, Hadid S. Smooth Expansion to Solve High-Order Linear Conformable Fractional PDEs via Residual Power Series Method: Applications to Physical and Engineering Equations. Ain Shams Eng J (2020) 11(4):1243–54. doi:10.1016/j.asej.2020.03.016
25. El-Ajou A, Oqielat Ma. N, Al-Zhour Z, Momani S. Analytical Numerical Solutions of the Fractional Multi-Pantograph System: Two Attractive Methods and Comparisons. Results Phys (2019) 14:102500. doi:10.1016/j.rinp.2019.102500
26. El-Ajou A, Al-Zhour Z, Oqielat Ma., Momani S, Hayat T. Series Solutions of Nonlinear Conformable Fractional KdV-Burgers Equation with Some Applications. Eur Phys J Plus (2019) 134(8):402. doi:10.1140/epjp/i2019-12731-x
27. Shqair M, El-Ajou A, Nairat M. Analytical Solution for Multi-Energy Groups of Neutron Diffusion Equations by a Residual Power Series Method. Mathematics (2019) 7(7):633. doi:10.3390/math7070633
28. Eriqat T, El-Ajou A, Oqielat Ma. N, Al-Zhour Z, Momani S. A New Attractive Analytic Approach for Solutions of Linear and Nonlinear Neutral Fractional Pantograph Equations. Chaos, Solitons & Fractals (2020) 138:109957. doi:10.1016/j.chaos.2020.109957
29. El-Ajou A. Adapting the Laplace Transform to Create Solitary Solutions for the Nonlinear Time-Fractional Dispersive PDEs via a New Approach. Eur Phys J Plus (2021) 136:229. doi:10.1140/epjp/s13360-020-01061-9
30. Srivastava VK, Awasthi MK, Tamsir M. RDTM Solution of Caputo Time Fractional-Order Hyperbolic Telegraph Equation. AIP Adv (2013) 3(3):032142. doi:10.1063/1.4799548
31. Abbas S, Benchohra M. Fractional Order Partial Hyperbolic Differential Equations Involving Caputo Derivative. Stud Univ Babes-bolyai Math (2012) 57:469–79.
32. Akilandeeswari A, Balachandran K, Annapoorani N. Solvability of Hyperbolic Fractional Partial Differential Equations. J App Anal Comp (2017) 7:1570–85. doi:10.11948/2017095
33. Ashyralyev A, Dal F, Pinar Z. On the Numerical Solution of Fractional Hyperbolic Partial Differential Equations. Math Probl Eng (2009) 2009:1–11. doi:10.1155/2009/730465
34. Khan H, Shah R, Baleanu D, Kumam P, Arif M. Analytical Solution of Fractional-Order Hyperbolic Telegraph Equation, Using Natural Transform Decomposition Method. Electronics (2019) 8(9):1015. doi:10.3390/electronics8091015
35. Modanl M. Two Numerical Methods for Fractional Partial Differential Equation with Nonlocal Boundary Value Problem. Adv Differ Equ (2018) 333:333. doi:10.1186/s13662-018-1789-2
36. Lin Y, Xu C. Finite Difference/spectral Approximations for the Time-Fractional Diffusion Equation. J Comput Phys (2007) 225:1533–52. doi:10.1016/j.jcp.2007.02.001
37. Meerschaert MM, Tadjeran C. Finite Difference Approximations for Fractional Advection-Dispersion Flow Equations. J Comput Appl Maths (2004) 172:65–77. doi:10.1016/j.cam.2004.01.033
38. Kochubei AN. Fractional-hyperbolic Systems. Fract Calc Appl Anal (2013) 16:860–73. doi:10.2478/s13540-013-0053-4
39. Hendy AS, Macías-Díaz JE, Serna-Reyes AJ. On the Solution of Hyperbolic Two-Dimensional Fractional Systems via Discrete Variational Schemes of High Order of Accuracy. J Comput Appl Maths (2019) 354:612–22. doi:10.1016/j.cam.2018.10.059
Keywords: hyperbolic systems, power series, analytical–numerical methods, fractional derivatives, Laplace transform
Citation: El-Ajou A and Al-Zhour Z (2021) A Vector Series Solution for a Class of Hyperbolic System of Caputo Time-Fractional Partial Differential Equations With Variable Coefficients. Front. Phys. 9:525250. doi: 10.3389/fphy.2021.525250
Received: 08 January 2020; Accepted: 30 April 2021;
Published: 28 May 2021.
Edited by:
Feliz Minhós, University of Evora, PortugalReviewed by:
Yang Liu, Inner Mongolia University, ChinaNdolane Sene, Cheikh Anta Diop University, Senegal
Sania Qureshi, Mehran University of Engineering and Technology, Pakistan
Copyright © 2021 El-Ajou and Al-Zhour. 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: Ahmad El-Ajou, YWpvdTQ0QGJhdS5lZHUuam8=; Zeyad Al-Zhour, emFsemhvdXJAaWF1LmVkdS5zYQ==