Jacobi convolution polynomial for Petrov-Galerkin scheme and general fractional calculus of arbitrary order over finite interval

Pavan Pranjivan Mehta,Gianluigi Rozza
2024-11-12
Abstract:Recently, general fractional calculus was introduced by Kochubei (2011) and Luchko (2021) as a further generalisation of fractional calculus, where the derivative and integral operator admits arbitrary kernel. Such a formalism will have many applications in physics and engineering, since the kernel is no longer restricted. We first extend the work of Al-Refai and Luchko (2023) on finite interval to arbitrary orders. Followed by, developing an efficient Petrov-Galerkin scheme by introducing Jacobi convolution polynomials as basis functions. A notable property of this basis function, the general fractional derivative of Jacobi convolution polynomial is a shifted Jacobi polynomial. Thus, with a suitable test function it results in diagonal stiffness matrix, hence, the efficiency in implementation. Furthermore, our method is constructed for any arbitrary kernel including that of fractional operator, since, its a special case of general fractional operator.
Numerical Analysis,Mathematical Physics,Analysis of PDEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop generalized fractional - order calculus on a finite interval and propose an efficient Petrov - Galerkin numerical method based on Jacobi convolution polynomial for solving generalized fractional - order differential equations. Specifically, the paper mainly focuses on the following aspects: 1. **Development of generalized fractional - order calculus**: The paper first extends the work of Al - Refai and Luchko on a finite interval to an arbitrary order, introducing a new class of kernel functions, which enables the generalized fractional - order calculus to be applied to a wider range of physical and engineering problems. These kernel functions are no longer limited to the power - law form, but can be in any form, thus increasing the flexibility and applicability of the model. 2. **Introduction of Jacobi convolution polynomial**: In order to improve the efficiency of the numerical method, the paper introduces Jacobi convolution polynomial as a basis function. An important property of this type of basis function is that its generalized fractional - order derivative is a shifted Jacobi polynomial. This makes the stiffness matrix a diagonal matrix after choosing appropriate test functions, thus greatly improving the computational efficiency. 3. **Construction of Petrov - Galerkin scheme**: Based on the above basis functions, the paper constructs an efficient Petrov - Galerkin numerical scheme for solving generalized fractional - order differential equations. This scheme is applicable not only to generalized fractional - order differential equations, but also to traditional fractional - order differential equations, since the latter is a special case of the former. Through these methods, the paper aims to provide a theoretical framework and numerical tools to solve complex fractional - order differential equation problems encountered in physics and engineering, especially those problems that require accurate solutions within a finite interval.