Approximation and application of the Riesz-Caputo fractional derivative of variable order with fixed memory

In this paper, the Riesz-Caputo fractional derivative of variable order with fixed memory is considered. The studied non-integer differential operator is approximated by means of modified basic rules of numerical integration. The three proposed methods are based on polynomial interpolation: piecewise constant, piecewise linear, and piecewise quadratic interpolation. The errors generated by the described methods and the experimental rate of convergence are reported. Finally, an application of the Riesz-Caputo fractional derivative of space-dependent order in continuum mechanics is depicted.


Introduction
In the last decades, fractional calculus has become a powerful mathematical tool to model a lot of physical processes. Because of the nature of fractional derivatives (nonlocal operators), they are especially used in describing nonlocal and multiscale phenomena. In this regard, we refer the reader works on the following topics: viscoelastic materials [6,23,28], heat conduction [24,40,54,56], diffusion processes in complex systems [12,30,32], nonlocal constitutive laws [13,19,46], microscopic interaction forces [18,25], nonlocal stress-strain constitutive relations [3], nonlocal beams [2,8,10,41,47]. A fairly broad overview of various aspects of fractional calculus applications can be found in the review paper by Failla and Zingales [20].
Notwithstanding this branch of mathematics is not new (the creation of fractional calculus dates back to the birth of the classical theory of differential calculus [27]) there are still remain a lot of problems that require attention. In particular, the branch of fractional calculus, that is devoted to differential equations containing simultaneously both the left and right fractional derivatives, is a research area where exists a large number of an unsolved/unexplored issues [15,21,29]. No exact solutions, for many of the aforementioned equations, have been found till now. For this reason, many researchers focus their attention on approximate solutions of equations containing various types of combinations of fractional derivatives. For this reason, many researchers focus their attention on approximate solutions of equations containing various types of combinations of fractional derivatives [1,7,11,14,22,53].
On the other hand, certain types of phenomena in physics are better depicted when the order of the fractional operator is variable (it depends on the spatial or time variable) [26,43]. Promising approaches have been presented in papers where authors studied: diffusion processes in an inhomogeneous and heterogeneous medium [16,33,44,[49][50][51]55], control theory [17,34,35], mechanics [4,31,36,38,39,42]. More details of the available applications of a variable order fractional calculus in the area of scientific and engineering modelling are presented in the following review papers: [37,52].
In this paper we focus on the Riesz-Caputo fractional operator of variable order with fixed memory and its application in a strong form of the space-fractional continuum (1D case), which is just a straightforward extension of the results presented in [45]. Such approach allows, when adopted to the mechanical applications, for clear physical interpretation and moreover gives smooth passage to classical approach as a limit case [5,48]. We propose three methods, based on polynomial interpolation, to approximation of the studied fractional operator. All presented numerical sachems were implemented in Python using library for real and complex floating-point arithmetic with arbitrary precision (http:// mpmath.org/). where n À 1\aðxÞ\n. In following, based on our previous work [9] where the numerical schemes for differential operators of constant order a were developed, we elaborate the algorithms for the Riesz-Caputo derivative of variable order. For discretisation, we generate the grid of ðN þ 1Þ -equidistant nodes: Moreover, we introduce the following notation:

Method I: Approximation based on piecewise constant interpolation
where where and

Method III: Approximation based on piecewise quadratic interpolation
where with Fig. 1 Approximations of the Riesz-Caputo derivative of functions x 2 , exp x, sin x and cos x for space-dependent order a 3 Numerical analysis: errors and experimental rates of convergence In this subsection, we report on a few numerical examples for checking the accuracy of the presented methods. Algebraic problems given by Eqs. (4), (7), and (10) were built using original procedures developed by the authors in Python. Figure 1 shows graphs of approximations of the Riesz-Caputo derivative of functions x 2 , exp x, sin x, and cos x for space-dependent orders a x ð Þ ¼ 1 2 þ x 2 and a x ð Þ ¼ 1 À x 2 . The calculations presented on the plots have been performed for N ¼ 2000 and ' ¼ 0:2.
In Tables 1, 2, and 3 we present the absolute errors of numerical calculations carried out using the proposed methods (Method I, Method II, and Method III) for the following functions: (Table 1), Table 3).
The reported errors have been determined by comparing numerical results with exact values (using series representation of the Riesz-Caputo derivative [9]). We take a particular value of the Riesz-Caputo operator and compare it with a value received from the series at the same specific point. In this case, for each node x i the order of derivative has the form a x i ð Þ and is an Analysing results presented in tables, one can observe that the convergence is the fastest for the piecewise quadratic interpolation (scheme 10), a little slower in case of the piecewise linear interpolation (scheme 7) and the slowest for the piecewise constant interpolation (scheme 4). Additionally, the errors generated by Method I are the biggest one, while for Method III are the smallest one.

Application of Riesz-Caputo derivative of variable order for the space-fractional continuum mechanics -numerical study
In this section we present the application of the Riesz-Caputo fractional derivative of variable order in continuum mechanics. We consider one-dimensional tension problem of the space-fractional continua in the following form [46]: with Dirichlet's boundary conditions where b E ¼ 0:1 and l ¼ 1. In Eq. (18), x denotes spatial variable, ' is the length scale, f denotes the  (left side) and a x ð Þ ¼ 1 À x displacements, b is the body force, and E stands for the Young modulus. We start from introducing the spatial discretization for the analysed one-dimensional continuum body that is presented in Fig. 2 Next, after applying the backward difference formula, for the first order derivative, in Eq. (15) we obtain: Finally, the following system of linear equations equations can be formulated: where The system (17) (containing N þ 1 equations) has been solved for two types of function a, namely a x ð Þ ¼ 1 2 þ x 2 and a x ð Þ ¼ 1 À x 2 and the following parameters: N ¼ 200 and ' 2 f0:01; 0:1; 0:2; 0:4g. The obtained results are presented in Fig. 3.
The analysis of Fig. 3 allows us to draw the conclusion that the deformation of the space-fractional continuum body solely depends on the length scale ' and the order a. Moreover, the introduction of variable order of fractional continua makes the model more flexible, thus in perspective validation of this model with experimental data should be easier.

Conclusions
The approximation of the Riesz-Caputo fractional derivative of variable order and with fixed memory can be elaborated using different approaches. Herein, three original methods are considered based on the polynomial interpolation: piecewise constant (Method I), piecewise linear (Method II), and piecewise quadratic (Method III) interpolation. As presented, the convergence is the fastest for Method III, a little slower in case of Method II and the slowest for Method I. Moreover, the errors generated by Method I are the biggest one, while for Method III are the smallest one. Finally, the presented schemes can be directly applied for modelling of scale effect is terms of spacefractional continuum mechanics, making its ability to mimic experimental results higher.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.