Abstract
The theory of Bessel functions is a rich subject due to its essential role in providing solutions for differential equations associated with many applications. As fractional calculus has become an efficient and successful tool for analyzing various mathematical and physical problems, the so-called fractional Bessel functions were introduced and studied from different viewpoints. This paper is primarily devoted to the study of developing two aspects. The starting point is to present a fractional Laplace transform via conformable fractional-order Bessel functions (CFBFs). We establish several important formulas of the fractional Laplace Integral operator acting on the CFBFs of the first kind. With this in hand, we discuss the solutions of a generalized class of fractional kinetic equations associated with the CFBFs in view of our proposed fractional Laplace transform. Next, we derive an orthogonality relation of the CFBFs, which enables us to study an expansion of any analytic functions by means of CFBFs and to propose truncated CFBFs. A new approximate formula of conformable fractional derivative based on CFBFs is provided. Furthermore, we describe a useful scheme involving the collocation method to solve some conformable fractional linear (nonlinear) multiorder differential equations. Accordingly, several practical test problems are treated to illustrate the validity and utility of the proposed techniques and examine their approximate and exact solutions. The obtained solutions of some fractional differential equations improve the analog ones provided by various authors using different techniques. The provided algorithm may be beneficial to enrich the Bessel function theory via fractional calculus.
1. Introduction
The theory of special functions is a critical branch of modern mathematical analysis. During the past three decades, several new classes of special functions have been proposed as solutions of fractional differential equations (FDEs). No other special functions have received such detailed treatment in readily available treatises as Bessel functions. The investigation of such functions is an important problem in fractional calculus, which has earned much attention as real-life problems can be analyzed well. Fractional calculus appears in many branches of science, such as medicine, material sciences, electromagnetics, and fluid mechanics (see [1–4]). Many applications have been performed through FDEs, and their solution techniques could be found, for example, in [5–11].
We trace the existing efforts regarding fractional order derivatives. Several definitions of the fractional order derivative have been introduced by many famous authors, such as Euler, Fourier, Letnikov, Laurent, Grünwald-Letnikov, Caputo, and Riemann-Liouville. Other definitions have also been provided by Kilbas et al. and Miller and Ross in [3, 10]. The most popular definitions considered frequently in the literature are derivatives by Riemann-Liouville, Caputo, and Grünwald-Letnikov. Interestingly, each definition of the arbitrary order derivative captures only a few properties of the classical integral derivative. However, a few drawbacks exist; for instance, does not fulfill the Riemann-Liouville definition. In Caputo’s definition, is assumed to be differentiable; otherwise, one cannot use such a definition. Moreover, Liouville’s theorem in the fractional setting does not hold. Therefore, it is clear that all definitions of fractional derivatives seem deficient regarding certain mathematical properties, such as the rules of product, quotient, and chain. For more details concerning other properties of these fractional order derivatives, see, for example, [4] and the references therein.
Due to the mentioned arguments, a new definition of the fractional order derivative is needed to achieve suitable mathematical properties. Khalil et al. [12] introduced a well-extended definition of the noninteger order derivative called the conformable fractional derivative (CFD). This definition is formulated as follows:
Definition 1 (see [12]). For the initial real value , the conformable fractional derivative of a real function is defined by the following:
The initial value can be zero, and if the limit exists, then is called -differentiable.
Along with the CFD’s Definition 1, if is differentiable, then (see [12]). Moreover, if Definition 1 holds for , we obtain . Additionally, we have ; hence, , which indicates that relies on some functions; it is unreasonable. Therefore, the CFD definition [12] does not need to hold for zero order. For the conformable fractional integral, we state the following definition as given in the following [12]:
Definition 2 (see [12]). Let Then, for any , the conformable fractional integral of order of is defined as follows:
Definition 1 depends entirely on the basic limit like the classical order derivative. Furthermore, Definition 1 fulfills various classical properties, such as the mean value theorem and the product, quotient, and chain rules. Moreover, this definition is provided with the Leibniz rule, in which other fractional derivatives can not achieve (see [13]). Growing attention has been paid to explore the conformable derivatives due to the enormous number of their meaningful applications in many fields of science. Abul-Ez et al. [14] introduced a comprehensive study on the conformable fractional Legendre polynomials. They presented the shifted conformable fractional Legendre polynomials and described an applicable scheme using the collocation method to solve some fractional differential equations (FDEs) in the sense of conformable derivative. Recently, the conformable fractional Gauss hypergeometric function and a class of conformable fractional differential equations through that function were treated in [15]. Further interesting ideas on the conformable derivative can be found in the work by [16–23].
Note that some authors have demonstrated that the conformable derivative is not the same as a fractional order derivative, but it is a first-order derivative multiplied by an additional factor (see for example [23]). Hence, Definition 1 seems to be a natural extension of the conventional order derivative to noninteger order loosing memory effect. In addition, a new approach for finding fractional operators was introduced by Antagan and Baleanu [24] with a nonsingular Mittag-Leffler kernel with a memory effect.
Returning to the purpose of the present work, we observe that Bessel functions are playing a significant role in investigating the solution of important differential equations (for example, see [25]). The theory of Bessel functions is usually used when solving problems related to information theory, nuclear physics, radiophysics, and hydrodynamics. Recently, as in [26–30], a resurgence of interest has occurred in the study of Bessel functions in the framework of fractional calculus theory. Along with the work in [26, 27], we employ conformable fractional order Bessel functions (CFBFs) to solve problems of a fractional nature. The study of a Bessel function of half-integer order led to discovering another interesting class of orthogonal polynomials called the Bessel polynomials. Many authors have used these polynomials. For example, Yüzbaşi et al. [31] solved linear integral, differential, and integro-differential equations, while Parand et al. [32] applied Bessel functions to solve nonlinear Lane-Emden equations. In [33], fractional optimal control problems were solved using the Bessel collocation method.
The present work proposes an approach to approximating the solution for some important linear and nonlinear FDEs in the conformable sense. The paper is designed with two objectives. The first is to establish some interesting properties of fractional Laplace-type integrals of functions via CFBFs. Then, we use the obtained results to establish the possible solutions of conformable fractional kinetic equations through CFBFs. The second objective is concerned with developing applications of the fundamental process of the proposed approach in terms of CFBFs. To achieve that, we derive an orthogonality relation, expand functions in terms of the truncated CFBFs, and effectively formulate a scheme involving the collocation method which employed to provide solutions of certain types of linear and nonlinear CFBEs.
The structure of this paper is organized as follows. The needed concepts and features of CFD are collected in Section 2. Next, Section 3 establishes useful properties of Laplace transforms in the sense of CFBFs, with some applications to solve a new type of conformable fractional kinetic equations. Section 4 is divided into four subsections. Section 4.1 provides essential results on orthogonality relations. A brief study on an expansion of any analytic function employing CFBFs is the subject of Section 4.2. In Section 4.3, we construct an algorithm for solving various kinds of problems using CFBFs through the collocation method. Section 4.4 presents the concepts that have been developed through previous subsections to solve some linear and nonlinear conformable fractional differential equations (CFDEs), including the nonlinear Riccati FDE. Concluding remarks are provided in Section 5.
2. Preliminaries and Basic Concepts
The Bessel equation is a special case of the Sturm-Liouville problem, and it can be written as [34]
In view of formula (1), the authors of [26] solved the following conformable fractional Bessel equation: around the regular singular point and introduced the CFBFs of the first kind as its solution such that
Moreover, they investigated in [26] some of its recurrence relations from which we may mention
In the following, we are about to recall some essential definitions and results which are needed in the sequel.
Definition 3. The Gauss hypergeometric function is defined by (see [35]) where stands for the familiar Pochhammar symbol which can be written in terms of Gamma function as
Definition 4 (see [36]). The function where and refer to its numerators and denominators, respectively, is called the Fox-Wright function, and it can be defined by the formula such that , where .
In particular, when in Definition 4, then the function immediately reduced to the generalized hypergeometric function (see [35]). Abdeljawad [16] defined the fractional Laplace transform in the conformable sense as follows:
Definition 5 (see [16]). For a real valued function , the conformable fractional Laplace transform of noninteger order is given by
The inverse fractional Laplace transform is the transformation of a fractional Laplace transform into a function of time. If , then is the inverse fractional Laplace transform of , and it can be written as
Remark 6. If , then (10) is the classical definition of the Laplace transform of integer order.
Furthermore, the author in [16] gave the following interesting results.
Lemma 7 [16]. For a real valued function satisfying , the following relations hold true: (1), where (2)(3)(4)
3. Fractional Laplace Transform of the CFBFs
In this section, we derive some new interesting fractional Laplace-type integrals of functions involving CFBFs. Then, as an application, we are going to employ the obtained results in order to find the possible solutions of the fractional kinetic equations in the conformable sense associated with CFBFs.
3.1. Fractional Laplace Integral Formulas
Theorem 8. Let , be the CFBFs; then,
Proof. Owing to the definition of CFBFs (5) and applying the conformable fractional Laplace transform operator of order as stated in Lemma 7, we have According to (3) of Lemma 7, one can see Using the identity , we obtain Therefore, the result is established.
Now, consider the Fox-Wright function defined in (9) to deduce the following important results.
Theorem 9. Let be the CFBFs. Then, the following relation is satisfied:
Proof. By combining (5) and (10), we get Due to (3) of Lemma 7, it follows that as required.
Theorem 10. Let be the CFBFs. Then,
Proof. As proceeded in the proof of Theorem 8 then, relying on Equation (5) and Lemma 7, we conclude that
Theorem 11. Let be the CFBFs. Then, the following identity holds:
Proof. The reduction of (5) and (10) yields which ends the proof.
The above obtained results provide the necessary tools which enable us to carry out the following interesting study.
3.2. Fractional Kinetic Equations Associated with the CFBFs and Their Solutions via Fractional Laplace Transform
One of the most important equations in mathematical physics and natural sciences is the kinetic equation, which describes the continuity of motion of substances. Therefore, many researchers investigated extensions and generalizations of this equation in the context of various fractional calculus operators. For such type of work, we refer for example to [37, 38]. We begin by briefly reviewing these previous efforts; then, we introduce our extended form of the fractional kinetic differential equation associated with the CFBFs. Assuming that denotes an arbitrary reaction which depends on time, refers to the destruction rate, and is the production rate on , Haubold and Mathai [37] characterized the FDE of the quantities and , by the formula where for . In the case where spatial fluctuation or inhomogeneities in the quantity is neglected, the authors in [37] handled the following equation: where the primary condition gives the number density of species at time and constant . Equation (24) is known as the standard kinetic equation. Alternatively, if the index is neglected, then, by integrating the standard kinetic Equation (24), one can get where is the standard integral operator. Equation Equation (25) has been extended to the fractional setting in the form (see [37]). where denotes the standard fractional integral operator in Riemann-Liouville sense (see [10, 39]).
Now, consider the conformable fractional kinetic equation in the form: where is the conformable fractional integral of order in the frame of Definition 2. Accordingly, we develop here a new type of generalization of the fractional kinetic differential equation in the conformable sense involving the fractional-order Bessel function in view of a fractional Laplace transform.
Remark 12. The solutions we are going to conclude for the conformable fractional kinetic equations will be determined through the generalized Mittag-Leffler function [36], which is defined as
Theorem 13. For and , the following conformable fractional equation has a solution in the form
Proof. Following Abdeljawad [16], then, in virtue of , we have where defined in (10) and is the conformable fractional integral operator (2). Acting by the conformable fractional Laplace transform on both sides of Equation (30) implies that Combining Equation (5) and Equation (31) leads to Considering Equation (10) with (34), we obtain Therefore, With the aid of the inverse Laplace transform (10), it follows that
Theorem 14. The solution of the following conformable fractional equation is given by
Proof. Operating the conformable fractional Laplace transform on both sides of Equation (38), we get In view of Equations (5) and (31), we obtain Hence, using (3) of Lemma 7 implies from which one can obtain Therefore, In view of the inverse fractional Laplace transform (10), it follows that Thus, the result is established.
4. Orthogonality of the CFBFs
Understanding the orthogonality relation of the CFBFs is mandatory to compute coefficients of series whose terms include the CFBFs. These series represent solutions of the FDEs as we will encounter in the application part of this section. Along with [40] and in view of the CFD definition (1), we introduce the following interesting results on orthogonality which will be useful in the current study.
4.1. An Orthogonal Relation of the CFBFs
Theorem 15. The orthogonality relation of the CFBFs is deduced over with respect to the weight function by the following: where is the familiar Kronker delta function and are distinct roots of
Proof. Since is a solution of the CFBE (4), it follows that which satisfies the more general equation It is convenient to reformulate (47) in the following way: Consequently, and satisfy the following CFDEs, respectively: Multiplying (49) by and (50) by and then subtracting the resulting equations produce In view of the conformable fractional integral formula (2) over , we obtain By performing integration by parts [16] on the right-hand side divided by the factor , one can conclude that Hence, according to the values of and , we consider the following two cases: (i)If and by hypothesis , then the right-hand side of (53) vanishes(ii)If , then the resulting integralcreates an interest to look at. In order to deduce its value, we take the limit of (53) as . As the right-hand side in (53) approaches the indeterminate form in the limit, we apply L’Hopital’s rule, which leads to Now, using the following recurrence relations of CFBFs [26] it follows that
As a special case of Theorem 15, the following result can be easily verified.
Corollary 16. The CFBFs are orthogonal over with respect to the weight function and
The results obtained in the current Section 4.1 treated the topic of orthogonal polynomials which pave the way to discuss a function representation via a series of the CFBFs.
4.2. Expansion of Functions via CFBFs
The classical theory of expressing analytic functions as expansions in terms of an arbitrary set of orthogonal polynomials can be described as the backbone of many topics in analysis. It was originated by several authors to whom we may mention Boas and Buck [41], Faber [42], and Whittaker and Gattegno [43] and later on in higher dimension by Abul-Ez et al. [29, 30, 44–46]. In the usual classical calculus, we found that not all functions have the Taylor power series representation around specific points, but this is not the case in the theory of conformable fractional calculus. This fact has been shown by Abdeljawad [16], where he also proposed the expansion of the fractional power series for an infinity -differentiable function through the fractional Taylor series. The expansion of a given real function in a series of Bessel functions is extremely useful in determining the solution of certain FDEs involving radial symmetry [40]. Related to the work of finding the expansion of a given function by means of Bessel polynomials in the higher-dimensional context, see the work given in [29, 30]. Using the orthogonality property (46), one can easily represent a given function over the interval by a series of Bessel functions such as where , and are determined by
As the topic of expansions of an arbitrary function, of either a real or a complex variable, into a series of polynomials has not been fully explored, we believe that several open problems remain untouched in particular in the framework of fractional calculus.
4.3. Applications
In this subsection, we intend to indicate the efficiency and applicability of the results developed in this study. Precisely, we construct a scheme which will be employed to solve some linear and nonlinear CFDEs. In this concern, we first define the truncated CFBFs of the first kind. Then, we introduce the noninteger derivative in the conformable context of an approximated function expanded in terms of the CFBFs. Several important examples including famous FDEs have been comprehensively treated, and their solutions have been compared to other existing methods in the literature to show the consistency and accuracy of our proposed method.
We begin by considering the truncated CFBFs of the first kind as follows: where is a positive integer such that and . For , we have
Figures 1 and 2 show the graphs of the truncated CFBFs when and by taking various values of and .
Let be a function defined over ; then, can be expanded in terms of CFBFs as follows:
Thus, the following truncated series for is supposed to be
Theorem 17. The noninteger derivative of order of the CFBFs in the conformable sense is given by where
Proof. The linearity of the conformable derivative (see [12, 16]) leads to as required.
Remark 18. If , where , then,
Theorem 19. Let be an approximated function given by means of the truncated formula of CFBFs (61). Then,
Proof. The induction of Theorem 17 and in view of the linearity property leads to the required result.
4.3.1. Proposed Scheme
(1) Linear Multiorder CFDEs. Suppose that the generalized linear multiorder CFDE is given in the form subject to the initial conditions where denotes to the CFD of order are known to be a continuous function, and are some constants.
Assume that the solution of the CFDE (69) can be given in the form
Substituting (71) into (69) and using Theorem 19, we have
Collocating Equation (72) at the points , we get
These equations give linear algebraic equations. Combining (71) and (70) produces -algebraic equations. Thus, we obtain a linear algebraic system of equations in the unknown . By solving this system, we obtain the solution of Equation (69) associated with the initial condition (70).
(2) Nonlinear Multiorder CFDEs. Consider the generalized nonlinear multiorder CFDE such that with the initial conditions where stands for the CFD of order is a nonlinear operator, and are given constants. A similar procedure as in establishing (69) can be used to approximate in terms of CFBFs. Thus, combining (71) with (74) and using Theorem 19 yield
Collocating Equation (76) at the points , thus, we have the following:
Equation (77) determines nonlinear equations. In virtue of (71) and the initial conditions (75), we get equations. Immediately, one can get a system of nonlinear equations in the unknown As usual, Newton’s iterative method can be used to solve this system. Thus, the solution of the nonlinear CFDE (74) with (75) can be deduced.
We can briefly clarify the achieved advantages of employing the proposed method as follows. The fractional-order Bessel functions approximate the fractional function with more accuracy. This feature has made the FBFs more effective than Bessel functions in solving the fractional problems. As the values of coefficients in Bessel polynomials are smaller than the coefficients of Chebyshev, Legendre, and Bernoulli polynomials, the computational error in the current method is less. Furthermore, our detailed treatment to these FDEs using the collocation method is aimed at encouraging the use of such approach which allows reaching the required solutions with ease and accuracy.
In the following subsection, we discuss several numerical examples to demonstrate the consequences of the above-mentioned features.
4.4. Illustrative Examples
In view of the above arguments, it is interesting to employ the provided techniques to solve some useful FDEs, as we shall see through the following examples.
Example 20. Consider the famous Bagley-Torvik equation subject to the initial conditions
This equation has been solved using Legendre polynomials in view of the conformable fractional sense in [14], as well as in [47–50] in the sense of Caputo derivative. The exact solution of (78) was given in ([47, 49]) as ; then for and in (71), one gets the approximated solution of (78) in the form
Owing to initial condition (79), then (80) gives
Using the collocation point and in view of (81), one can get
Hence, Therefore, we have , which coincides with the one given in ([47, 49]).
A similar procedure can be carried out as in Example 20, so that we may summarize the corresponding details for two different CFDEs as follows.
Remark 21. (i)The problem in Table 1 has been treated in [47–50] with the Caputo fractional derivative using various methods such as Legendre polynomials, homotopy perturbation method, and homotopy analysis method(ii)The problem in Table 2 has been manipulated in [51] by means of Legendre polynomials with conformable derivative and in [52] in view of Laguerre polynomials involving Caputo derivative
Example 22. Suppose that we have the following linear fractional differential equation in the form with the initial condition
The exact solution of (83) is (see [47]). This problem has been evaluated using various methods (see for example [47, 49, 53]) in the sense of the Caputo fractional derivative. Putting in (71) with , then using the presented technique, we computed the approximate solution of the problem (83) for the values . We indicate the approximate solution through Figure 3, while Figure 4 compares the obtained solution with the exact solution of (83). Moreover, the corresponding absolute errors of our approximate solution are displayed in Table 3. From Figure 4, it is clear that the approximate solutions converge to the exact solutions.
Example 23. It is well known that the theory of Bessel functions is connected with the Riccati equations. In fact, Bessel functions are defined as solutions of Bessel equations, which can be derived from the Riccati equations. This motivates us to consider the following the nonlinear Riccati fractional differential equation: subject to the initial condition
This problem has the exact solution (see [21]), and it has been treated by many authors using various methods such as the collocation method, operational matrix method, homotopy perturbation Pade technique, homotopy analysis method, modified homotopy perturbation (MHP) method, modified variational iteration method, B-spline operational matrix method, and polynomial least squares method (see, for example, [54, 55]. In all these methods, the authors used the Caputo fractional derivative.
Now, regarding the proposed method here, the approximated analytical solution of the initial value problem (85) has been computed for with various values of . Figure 5 shows the approximate solutions for and various values of . Table 4 illustrates both numerical and exact solutions of (85). The comparison between the obtained values of by the presented method and the modified Homotopy perturbation method (MHPM) given in [55] for and is shown in Table 5. Table 5 illustrates that our method is more accurate compared to other methods.
Example 24. Suppose that a system of fractional differential equations is given in the form with the primary conditions
We proceed the solution as follows. When , the exact solution of this system was given in [56, 57] as and Applying the same technique described above and by putting in (71) with , we have computed the approximate solution for the values . Figure 6 shows the graphs of approximate solution for various values of with Tables 6 and 7 display the corresponding absolute error to our approximate solutions for and various values of . Note that the absolute error tends to zero when the terms of CFBFs increase.
5. Conclusion
Exploring the multifaceted applications of Bessel functions in several fields of science with recognition of the growing impact conformable fractional calculus has in many applications, this paper exhibits further developments on the conformable fractional-order Bessel functions (CFBFs). The novelty of this research paper is determined by explaining the comparison between our obtained solutions (results) and the solutions in previously published articles by various authors. Such novelty may be described as follows. In the first part of the paper, some useful formulas concerning the properties of conformable fractional Laplace transforms are obtained. These formulas are successfully employed to obtain the solutions for a new type of fractional kinetic equations associated with the CFBFs in the conformable fractional sense which is generalized and developed in this study. These solutions are newly presented compared with those given by various authors (see [38, 58–61]).
In addition, an interesting orthogonal relation of the CFBFs is established. This gives rise to the discussion of functions expansions where we present a given function in a series of the CFBFs. Consequently, we investigate the analytical and approximate solutions of some linear and nonlinear CFDEs using a proposed scheme depending on the collocation method, involving CFBFs. Particular emphasis is paid to indicate that our approach, in some sense, is easily applicable and provides more accurate results with refined errors. The proposed scheme is used to approximate the solutions of linear, nonlinear CFDEs, and also systems of CFDEs. The numerical examples show the validity and efficiency of the method. Notably, it is demonstrated that the presented method works well and achieves better accuracy compared with exact solutions and with results obtained using other methods. Our scheme is beneficial in the way that using the fractional-order Bessel functions leads to more accurate results in approximating the fractional function. Because of this characteristic, FBFs are more effective than Bessel functions in treating fractional problems. Because Bessel polynomials have smaller coefficients than Chebyshev, Legendre, and Bernoulli polynomials, the computational error in the current approach is less. Importantly, the employment of the collocation method provides convenient and accurate solutions. Furthermore, all the given solutions are obtained differently, unlike those established by the authors in [47–55] and the references therein where other fractional derivatives such as Caputo derivative were used.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare there is no conflict of interest.
Authors’ Contributions
The authors have contributed equally to this work.
Acknowledgments
The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, Saudi Arabia, through research group program under grant R.G.P.1/86/42.