Chebyshev cardinal polynomials for delay distributed-order fractional fourth-order sub-diffusion equation
Introduction
The subject of fractional derivatives (which is a natural extension of the classical derivatives) is of interest to researchers because of its wide range of applications [1]. A very important point about fractional derivatives that is highly regarded is their non-local feature [2]. This property makes the future state of a dynamical system modeled by such derivatives dependent on its past and present states. In recent years, many fundamental problems in engineering and practical science have been modeled with good accuracy and finesse via fractional derivatives. For instance, they have been used in modeling the behavior of coronavirus [3], teletraffic problems [4], [5], dynamics of tuberculosis and HIV diseases [6], the vibration of a plate excited with supports movement [7], etc.
A serious point regarding fractional differential equations (i.e., equations that possess at least one term with a fraction derivative) is to find their analytic solution, which is often very arduous and in most cases, infeasible. In the face of such problems, the use of numerical methods is usually fruitful. Some well-known numerical approaches that have recently been applied for this category of problem include finite difference methods, meshless methods and spectral methods. For instance, see [8], [9], [10], [11], [12], [13], [14].
Various definitions of fractional derivatives have been provided in different texts. The most popular classical fractional derivatives are the Caputo and Riemann–Liouville derivatives [1]. As an extension of the traditional fractional derivatives, the distributed-order fractional derivatives are defined by integrating the classical cases over the order of the differentiation within a known domain [15], [16]. These derivatives play a middle role between the ordinary and fractional derivatives [17]. The differential equations produced by this kind of fractional derivatives can be propounded as generalizations of single and multiple fractional differential equations [16], [18] In recent years, fractional differential equations with distributed-order derivatives have been successfully applied to more accurately model many problems in the fields of signal processing [19], viscoelastic [20], electrochemistry [21], control [22], diffusion [23], etc.
Fractional diffusion equations of fourth-order model many important problems in engineering, such as beam vibration and groove formation on a metal surface [24]. There are several real models that are constructed via the diffusion-wave equations with distributed-order fractional derivatives. In [25], the authors used distributed-order fractional form of diffusion-wave equations to study waves in a viscoelastic rod. Some novel fractional differential equations based upon the Fokker–Planck equation with distributed-order fractional derivative are introduced in [26] for infiltration and absorption. In [27], the authors used the theory of distributed-order fractional derivatives regarding diffusion-wave equations in studying the radial groundwater flow and also investigated its application to pumping and slug tests.
Due to the importance of the fractional fourth-order sub-diffusion equations, we are interested in proposing a suitable technique to solve a delay version of such equations. So, we concentrate on the distributed-order fractional problem with , under the conditions and in which is a non-zero real number, , , and for are given smooth functions, and are given constants. Also, denotes the Caputo’s fractional differentiation of order relative to of . Moreover, the distribution function satisfies the following conditions [28]: Note that the linearized second-order numerical method is developed in [24] to solve the above problem. The implicit compact difference method is developed in [29] for solving the above fractional problem.
Since the distributed-order fractional differentiation is based upon the integration of a classical fractional differentiation, the polynomial basis functions can be useful in solving fractional problems with this type of derivative. Among basis functions of polynomial type, cardinal polynomials can be a good choice for solving these problems because their expansion coefficients are easily calculated. In recent years, cardinal polynomials have been extensively applied to solve diverse problems. Some of these uses can be found in [30], [31], [32], [33].
In this study, we utilize the Chebyshev cardinal polynomials (CCPs) [34] to generate a suitable and accurate numerical method for the above introduced delay distributed fractional problem. A derivative matrix is derived to compute the distributed-order fractional differentiation of these polynomials, and employed in making the provided method. The expressed methodology turns the problem solving into solving a simple algebraic problem by approximating the unknown function of the problem by these polynomials and exploiting the tau technique. Note that the interpolation property makes it easy to calculate nonlinear term.
The overview of this study is as follows: Required preparations are given in Section 2. A review of the CCPs with their properties is provided in Section 3. A matrix relationship is provided in Section 4 for distributed-order fractional differentiation of the CCPs. The generated methodology is explained in detail in Section 5. Several examples are investigated in Section 6. Eventually, the conclusion of this study is summarized in Section 7.
Section snippets
Definitions and preliminaries
Here, we provide some information that is required in this study.
Definition 2.1 Let is a differential function over and is a real number. The Caputo fractional differentiation of order of is given by Note that for , we have .[1]
Property 2.2 For and , we have [1]
Definition 2.3 An -point Gauss–Legendre quadrature formula can be defined over as follows [35]: where
Chebyshev cardinal polynomials
Here is a brief overview of CCPs along with some of their features.
Definition 3.1 For a given , the CCPs of order are generated on as [34] where and for .
Property 3.2 A continuous function on can be expressed by these polynomials as where and
Theorem 3.3 The first derivative of can
Distributed-order fractional derivative matrix
In the sequel, we introduce a matrix relationship for the distributed-order fractional derivative of the CCPs.
Lemma 4.1 Let ’s are the functions defined in (3.8). Then, for , we have
Proof Using the formula expressed in (3.1) and Property 2.2, the proof can be easily performed. □
Theorem 4.2 The Caputo fractional differentiation of order of in (3.8), can be computed as
Description of the method
Here, we establish a numerical method for the problem expressed in (1.1)–(1.3) by approximating the unknown in terms of the CCPs. To this end, we let where is an unknown matrix as From Corollary 3.4 and the above approximation, we get Using Theorem 4.3, the distributed-order fractional differentiation of (5.1), can be
Numerical examples
We here investigate some examples to show the reliability of our method. The following relations will be used to measure the values of error generated by the provided approach: where is the analytic solution and is the solution derived using the proposed scheme.
Example 1 Consider the equation where [24], [29]
Conclusion
In this paper, the Chebyshev cardinal polynomials were employed to make a suitable scheme for delay distributed-order time fractional fourth-order sub-diffusion equation. Because of this, some matrix relations regarding the distributed-order fractional derivative of these polynomials were derived and applied in the desired approach. The adopted technique was a matrix method to convert the main problem into a problem of algebraic equations. Some examples were considered for numerical evaluation
CRediT authorship contribution statement
M.H. Heydari: Conceptualization, Methodology, Software, Validation, Writing – original draft, Visualization, Supervision. M. Razzaghi: Conceptualization, Methodology, Software, Validation, Review and editing. J. Rouzegar: Conceptualization, Methodology, Software, Validation, Review and editing.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
References (35)
- et al.
SARS-COV-2 rate of spread in and across tissue, groundwater and soil: A meshless algorithm for the fractional diffusion equation
Eng Anal Bound Elem
(2022) Multi-fractional generalized cauchy process and its application to teletraffic
Physica A
(2020)- et al.
On analysis generalization of TB-HIV dynamics by a two-scale reduction process
Results Phys
(2021) - et al.
A fractional viscoelastic model for vibrational analysis of thin plate excited by supports movement
Mech Res Commun
(2020) - et al.
A new meshless method of solving 2D fractional diffusion-wave equations
Appl Math Lett
(2022) - et al.
Numerical solution of an inverse source problem for a time-fractional PDE via direct meshless local Petrov–Galerkin method
Eng Anal Bound Elem
(2022) Finite difference method for solving fractional differential equations at irregular meshes
Math Comput Simulation
(2022)- et al.
A novel finite difference technique with error estimate for time fractional partial integro-differential equation of Volterra type
J Comput Appl Math
(2022) - et al.
Lucas polynomials based spectral methods for solving the fractional order electrohydrodynamics flow model
Commun Nonlinear Sci Numer Simul
(2022) - et al.
Spectral method for the two-dimensional time distributed-order diffusion-wave equation on a semi-infinite domain
J Comput Appl Math
(2022)
On the formulation and numerical simulation of distributed order fractional optimal control
Commun Nonlinear Sci Numer Simul
A numerical method based on fractional-order generalized taylor wavelets for solving distributed-order fractional partial differential equations
Commun Nonlinear Sci Numer Simul
An algebraic fractional order differentiator for a class of signals satisfying a linear differential equation
Signal Process
Fractional differential equations in electrochemistry
Adv Eng Softw
On the formulation and numerical simulation of distributed-order fractional optimal control problems
Commun Nonlinear Sci Numer Simul
Numerical treatment of non-linear fourth-order distributed fractional sub-diffusion equation with time-delay
Commun Nonlinear Sci Numer Simul
Distributed-order fractional wave equation on a finite domain. stress relaxation in a rod
Internat J Engrg Sci
Cited by (5)
A hybrid method based on the Chebyshev cardinal functions/wavelets for time fractional coupled Klein–Gordon–Schrödinger equations
2023, Journal of Computational and Applied MathematicsA numerical method for distributed-order time fractional 2D Sobolev equation
2023, Results in PhysicsA direct transcription method for solving distributed-order fractional optimal control problems
2024, JVC/Journal of Vibration and Control