Next Article in Journal
Deep Contrastive Learning-Based Model for ECG Biometrics
Previous Article in Journal
Applicability of Deep Reinforcement Learning for Efficient Federated Learning in Massive IoT Communications
Previous Article in Special Issue
Differentiation of the Wright Functions with Respect to Parameters and Other Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Fading Memory Formalism with Mittag-Leffler-Type Kernels as A Generator of Non-Local Operators

Department of Chemical Engineering, University of Chemical Technology and Metallurgy, 8 Kliment Ohridsky, Blvd., 1756 Sofia, Bulgaria
Appl. Sci. 2023, 13(5), 3065; https://doi.org/10.3390/app13053065
Submission received: 27 January 2023 / Revised: 21 February 2023 / Accepted: 24 February 2023 / Published: 27 February 2023

Abstract

:
Transient heat conduction problems are systematically applied to the fading memory formalism with different Mittag-Leffler-type memory kernels. With such an approach, using various memories naturally results in definitions of various fractional operators. Six examples are given and interpreted from a common perspective, covering the most well-liked versions of the Mittag-Leffler function. The fading memory approach was used as a template and demonstrated that, if the constitutive equations are correctly built, it is also possible to directly determine where the hereditary terms are located in the models.

1. Introduction

1.1. Motivation of This Study

In recent years, there are “hot discussions” on definitions and applicability of fractional operators with different kernels [1,2,3,4,5,6,7]. In general, these discussions are from a mathematical point of view [4,8,9,10,11], because as a rule, the new operators are formulated constitutively, declaring their structures, with various kernels: correctly defined but depending only on the author’s imagination, without any physical backgrounds. The “mischance” in fractional modeling begins when operators with kernels not related to certain physical problems are implemented by replacements in already established integer-order models that make the results questionable.
Moreover, there are quests to see in the new operators what we already know without understanding that each operator has its physical basis and new properties. Moreover, there are attempts to “frame” all operators in already existing theories of fractional operator generalizations, not related to any real-world problems; This is an approach that, to some extent, is self-motivated and avoids the fact that fractional calculus is a branch of the applied mathematics; any steps to make it unified and generalized shrink the areas of applications and suppress its multifaceted nature, still developing in applied mathematical modeling.
Motivated by the idea to show some steps toward fractional modeling by applying constitutive equations with memories, obeying the causality principle, and having thermodynamic consistency, the following study was developed, using transient heat conduction with memory as an example.

1.2. Heat Conduction as a Transport Process with Memory

Transport phenomena with hereditaries that can be described by fractional operators (derivatives) are hot topics in modeling, and the approach to constructing physically adequate models is of primary importance. Models with memories can describe behaviors of elastic and viscoelastic solids, Non-Newtonian fluids, electromagnetic materials [12], as well as transport processes as flux models of diffusion and heat conduction [12,13,14,15].
It is well known that the earlier attempt to model transport processes with hereditaries comes from the works of Boltzmann [16] and Volterra [17]. We will skip the analysis of the Volterra constitutive equation but refer to [15] for a comprehensive analysis showing the differences between the Volterra approach and the fractional calculus modeling, which is the main task of this study. Precisely we address heat conduction as a physical problem where the constitutive equation of the flux can be modeled by fractional derivatives. Moreover, the main goal is to demonstrate that from the constitutive flux equation, following the fading memory principle (formalism) [18,19,20,21] and assuming different memory kernels, we may derive several operators widely used in fractional modeling. Hence, we will try to show that fractional operators appear logically from the correctly constructed constitutive equations rather than constituted as is generally done in mathematical papers- see the sequel for more details and reference quotations.
The first step in building the study’s exposition logically is to demonstrate the background in constitutive modeling of heat flux with a finite speed. Because of that, we will briefly present two examples: (i) How to create a constitutive equation for heat flux using the Caputo derivative [15] (Section 1.2.1) and (ii) constitutive equation formulation for heat conductors with fading memories [22] (Section 1.2.2).

1.2.1. Constitutive Fractional Heat Flux with Caputo Derivative [15]

Let us start with the energy balance (the First Law of Thermodynamics) in a rigid heat conductor
e x , t t = q x , t x + r
where e = ρ C p T is the internal energy, q x , t is the heat flux, and r is the heat generation (per unit volume). With the Fourier constitutive equation
q x , t = k T x , t T x , t x , k T x , t
where k T x , t is the thermal conductivity. This construction of the heat flux has no memory, propagates with infinite speed, and with the energy balance equation results in the parabolic Fourier model.
Attempts to model heat flux propagation with a finite speed were done by Cattaneo [23,24] and Vernotte [25] with slight modifications of the heat flux (via a Taylor series expansion of the heat flux, in time, and as a first-order approximation), namely
q x , t t = q x , t + k T x , t 1 τ q x , t x , τ > 0
with exponential memory (Jeffrey’s memory kernel) R t = exp t / τ , and a finite relaxation time τ , named the Cattaneo equation (see more details in Section 4.2).
As Fabrizio [15] demonstrated, replacing the time derivative in (2) by a fractional (Caputo) derivative, we get (in the original notations)
C D t α q x , t = q x , t + k T x , t T x , t x , α > 0 , t > 0
The thermodynamic restriction from the Second Law of Thermodynamics for a rigid heat conductor needs the following inequality to be obeyed [14]
η x , t t e x , t t T x , t + 1 T x , t 2 T x , t x
where the internal energy e x , t and the entropy η x , t are state functions.
In this framework, the following extension of (2) was proposed [15]
q x , t = C D t α q x , t x = k x 1 Γ 1 α t 1 t τ α t T x , τ x d τ = = k x α Γ 1 α t 1 t τ α + 1 T x , τ x T x , τ x d τ
The second version (in the second row) in (6) is an equivalent presentation of the Caputo derivative.
Finally, with the assumption that the internal energy e x , t depends only on T x , t (i.e., the density ρ and the heat capacity C p are temperature independent) and applying the energy balance equation the following model was obtained
ρ C p T x , t t = x k x α Γ 1 α t 1 t τ α + 1 T x , τ x T x , τ x d τ + r x , t
or in a more convenient form, for the presentation in this article, as
T x , t t = x a x 1 Γ 1 α t 1 t τ α d d τ T x , τ x d τ + r x , t , a x = k x ρ C p , t > 0
Remark 1.
It is worth noting and important for the entire discussion developed in this article that the time derivative in the energy balance equation (sometimes termed a continuity equation) is not replaced by a fractional derivative. The flux constitutive model is satisfying the causality principle [26,27,28], that between the cause (temperature gradient) and the result (output of the system), there is a time shift. This time shift in (7) is assured by the fractional Caputo derivative. This comment is of great importance for the model’s development in this article, where similar constructions of heat flux constitutive equations are developed with contributions of different fractional derivatives.
Remark 2.
The construction (6), in the light of the fractional diffusion model (we have a fractional diffusion model of thermal energy), does not consider a stationary state, that is, when the model reduces to the Fourier constitutive Equation (2) which is valid for long times where the relaxation of the heat flux vanishes. In addition, the analysis in [29] provides sufficient information on fractional modeling of transient heat conduction.

1.2.2. Constitutive Heat Flux with Linear Memory [22]

Following Giorgi and Gentili [22], the internal energy e x , t and heat flux q x , t , and following the fading memory principle [18,19,20,21] are described by the following linearized constitutive equations (with notations used in this work)
e x , t = e 0 + C p θ x , t + 0 R e s θ t s d s
q x , t = k 0 T x , t x 0 R q s θ t s d s
θ x , t = T x , t T 0 x , t , T 0 x , t = T x , 0 is temperature variation field. Further, R e t and R q t are relaxation functions (fading memories), and following the arguments of Day [21], the fading memory principle is equivalent to R e t , R q t L 1 R + and R q t .
In this context, the heat capacity C p t , and thermal conductivity k t are defined as [22]
C p t = C p 0 + 0 t R e s d s , k t = k 0 + 0 t R q s d s , C p 0 > 0
and C 0 = lim t C p t , and k 0 = lim t k t are equilibrium heat capacity, and the equilibrium thermal conductivity, that is when the memory integrals in (11) are vanishing.
At this point, we will discontinue the discussion of the thermodynamic principles of fading memory (it will be revisited in Section 3 devoted to model development with fractional operators), but it is critical to emphasize two principles in the constitutive model formulation:
1. The memory integral is an additive part (in the cases discussed here, a linear contribution) to the constitutive heat flux model. This gives us a great opportunity to apply fractional modeling through the use of adequately formulated relaxation functions (memory kernels), as we will do in the sequel to this article.
2. The energy balance equation remains unchanged since the First Law of Thermodynamics holds. This directly means that the time derivative of internal energy cannot be mechanistically replaced by fractional operators.
Remark 3.
We can see that we have two versions of constitutive heat flux formulations, but the energy balance equation remains unmodified since this is the First Law of Thermodynamics. Any changes in the balance equation, which is a very common step in fractional calculus publications, caused by replacing the time derivatives with fractional counterparts are violations of the First Law of Thermodynamics. The main task of mathematical modeling, as well as fractional modeling, is to solve correctly constructed models with physically interpretable results, which are often found in formalistic fractional models obtained by replacement.
Remark 4.
Before starting with the development of fractional operators with various Mittag-Leffler kernels, we have to say that there is no mathematics only for mathematics in this study. Fractional calculus is applied mathematics invoked by reality since the non-localities are natural phenomena in the transport processes. Because of that, the following exposition stresses the attention to the physical background of the modeled problems and their adequate interpretations.

1.3. The Main Task of This Article

The main task of this article is threefold:
1. To start the introduction of correctly formulated models of heat conductors with memory (actually, these are diffusion models, not diffusion models commonly modeled in heterogeneous media by application of fractional derivatives), applying the fading memory formalism.
2. To introduce different fractional operators in a logical way by using different relaxation functions (completely monotone Mittag-Leffler functions) and how these operators are related. From our point of view, this approach is of high importance for fractional modeling, where there are too many operators constituted without any physical background and then incorrectly applied in the modeling of real-world problems.
3. To analyze the developed constructions of diffusion (heat conduction) models and perform a comparative analysis regarding existing equations from the literature involving fractional operators based on various Mittag-Leffler-type kernels.

1.4. Paper Organization in the Sequel

The following part of this article is organized as follows: Section 2.1 provides the basic definitions and properties of the classical (one-parameter) Mittag-Leffler function (MLF) Section 2.1.1, two-parameter MLF Section 2.1.2, three-parameter MLF (Prabhakar function) Section 2.1.3, Prabhakar kernel Section 2.1.4. The main fractional operators based on Mittag-Leffler type kernels and their interrelations are briefly outlined in Section 2.2. The fading memory formalism with the concept for model build-up with different kernels are explained in Section 3 and Section 3.2. Constitutive heat flux models with different memories are developed in Section 4: Exponential memory Section 4.2, with Mittag-Leffler (one-parameter) memory kernel Section 4.3, with Prabhakar memory kernel Section 4.4. Additional experiments with some special well-kernels are developed in the Section 5: Rzanitsyn kernel Section 5.1, Miller–Ross kernel Section 5.2 and Rabotnov kernel Section 5.3, thus, demonstrating how new operators can be generated from a common viewpoint. An analysis of the developed models in parallel to some available equations from the literature is developed in Section 6.

2. Preliminaries on Mittag-Leffler Functions, Their Properties and Available Fractional Operators

To support the further development of the problems envisaged in this article, we need some preliminary information on the Mittag-Leffler functions and related operators to be briefly outlined.

2.1. Mittag-Leffler Functions and Related Kernels

2.1.1. One-Parameters Mittag-Leffler Function

The Mittag-Leffler function is defined as a power-series convergent in the whole complex plane [30,31]
E α z = k = 0 z k Γ α k + 1 , α > 0
and is an entire function of order 1 / α
It is a completely monotone (CM) function [30,32] that implies 1 m d m d z m E α z 0 , i.e., a negative function of class C for all m N , that it is a Bernstein function since its first derivative is CM. Our interest is oriented toward the function
E α t α = k = 0 1 k t α k Γ α k + 1 , 0 < α < 1 , t > 0
with a Laplace transform
L E α t α = s α 1 s α + 1 , L d d t E α t α = 1 s α + 1 , α > 0
For t 0 + the asymptotic expansion matches the stretched exponential [30]
E α t α = 1 t α Γ α + 1 + exp t α Γ α + 1 , t > 0
while for t the asymptote is a negative power-law t α
E α t α k = 0 1 k 1 t α k Γ 1 α k , t

2.1.2. Two-Parameters Mittag-Leffler Function

The two-parameter of the Mittag-Leffler type is defined as a series expansion
E α , β z = k = 0 z k Γ α k + β , z C , α > 0 , β > 0 , β C
For β = 1 we get the classical one-parameter Mittag-Leffler function E α , 1 z defined by (12)
From the definitions (12) and (17) it follows that [30,31,33]
E 1 , 1 z = k = 0 z k Γ k + 1 = k = 0 z k k ! = e z
E 1 , 2 = k = 0 z k Γ k + 2 = k = 0 z k ( k + 1 ) ! = 1 z k = 0 z k + 1 ( k + 1 ) = e z 1 z
The corresponding Laplace transforms are [30,31,33]
L t β 1 E α , β λ t α = s α β s α + λ = s β 1 + λ s α
L t α k + β 1 E α , β k λ t α = k ! s α β s α λ k + 1 , k = 0 , 1 , 2 ,
E α , β z is completely monotonic [34,35] for any α 0 , 1 , precisely if and only if α 0 , 1 and β α [36].
The integral of E α , β t is
0 t E α , β λ t α t β 1 d t = t β E α , β + 1 λ t α , β > 0
Applying the Riemann–Liouville fractional derivative to (17) we get [33]
0 D t μ t α k + β 1 E α , β k λ t α = t α k + β μ 1 E α , β μ k λ t α
When, k = 0 , λ = 1 and μ = m is an integer we get [33]
d m d t m t β 1 E α , β t α = t β m 1 E α , β m t α , m = 1 , 2 , 3 , . . . . .
For m = 1 we have
d d t t β 1 E α , β t α = t β 2 E α , β 1 t α
When β = 1 (one parameter Mittag-Leffler function), then (25) yields
d d t E α , 1 t α = E α , 0 t α t
In this context, the first derivative of E α t α , defined by (13), and the first derivative of E α , β t α are related as follows [30]
d d t E α , α t α = t 1 α E α , α t α = d d t E α t α
For λ 1 we have [33]
d d t E α , 1 λ t α = E α , 0 λ t α t
That is
0 t d d t E α , 1 λ t α = E α , 0 λ t α t

2.1.3. Three-Parameter Mittag-Leffler Function

Prabhakar [37] introduced the function (see also [38])
E α , β γ z = k = 0 γ k Γ α k + β z k k ! , α , β , γ C , Re α > 0
where γ k is the Pochhammer symbol
γ 0 = 1 , γ k = γ γ + 1 · · · γ + k 1 = Γ γ + k Γ γ , k = 1 , 2 , . . . .
For γ = 1 (30) reduces to (17). For γ = β = 1 we get the one-parameter Mittag-Leffler function (12).
The three-parameter Mittag-Leffler function E α , β γ z is completely monotone (CM) for 0 < α 1 and β α γ [39,40], precisely when 0 < α γ β 1 [40].
The differentiation of the generalized Mittag-Leffler function (30) [38] through the product z β 1 E α , β γ λ z α leads to
d m d z m z β 1 E α , β γ λ z α = z β m 1 E α , β m γ λ z α
For γ = 1 we have
d m d z m z β 1 E α , β λ z α = z β m 1 E α , β m λ z α
and for m = 1 the differentiation of E α , β m γ z is [41]
d d z z β 1 E α , β λ z α = z β 2 E α , β 1 λ z α
For β = 1 we get
d d t E α , 1 γ λ t α = E α , 0 γ λ t α t
and for γ = 1 we get (27).

2.1.4. Prabhakar Kernel

Prabhakar [37] studied a generalized Mittag-Leffler function (hereafter termed as Prabhakar kernel) (see also [40])
e α , β γ t = t β 1 E α , β γ t α , t 0
The Laplace transform of e α , β γ t , for some particular values of the parameters is [40,42,43]
L e α , β γ t = s α γ β s α + 1 γ , s > 0 , s α > 1
For λ 1 the Laplace transform and its inverse are [44]
L e α , β γ ± λ t α = s α γ β s λ γ , L 1 s α γ β s λ γ = t β 1 E α , β γ ± λ t α
The function e α , β γ t is locally integrable and completely monotone if the conditions 0 < α 1 and 0 < α γ β 1 are satisfied [45].
Asymptotically, for t we have [40]
e α , β γ t t β α γ 1 Γ β α γ , 0 < α γ < β 1 γ t β α γ 1 Γ β α γ , 0 < α γ = β 1
The integral of e α , β γ λ t for any t R + and λ C is [43]
0 t u β 1 E α , β γ λ u α d u = t β E α , β 1 γ λ t α
The m t h repeated integration of e α , β γ λ t = t β 1 E α , β γ λ t α [41]
I 0 m E α , β γ λ t α = t β + m 1 E α , β + m γ λ t α
In addition, for any m N the differentiation is is defined by (32) and in terms of time t is [43]
d m d t m t β 1 E α , β γ λ t α = t β m 1 E α , β m γ λ t α
For m = 1 the differentiation of t β 1 E α , β γ λ t α is (33) [41]
d d t t β 1 E α , β γ λ t α = t β 2 E α , β 1 γ λ t α

2.2. Fractional Operators with Mittag-Leffler-Type Kernels: Definitions and Interrelationships

Here, features of the main “artists” in this study employing different members of the Mittag-Leffler function family will be briefly encompassed. It is worth mentioning for a better understanding of the approach taken in the current work that all these fractional operators have been introduced constitutively, without proof of where they come from [38,40,41,42,43,45,46,47,48,49]. To be correct, the procedure is straightforward: we select a kernel, then insert it into constructions such as the Riemann–Liouville integral or the Caputo derivative [31,50] to obtain new operators. The sources contain no words that explain how these memories (kernels) are linked to certain physical models.
The three-parameter Mittag-Leffler function (30), as we mentioned, is CM’s entire function, and we focus the attention on the so-called Prabhakar kernel (36), which in particular reduces to [51]
e 1 , 1 1 λ t = exp λ t , e 1 , 1 0 λ t = t β 1 Γ β , e α , β γ λ = 0 , t = t β 1 Γ β

2.2.1. Prabhakar Integral

The Prabhakar integral with base point 0 is defined as [41,51] (in a Riemann–Liouville sense)
P I α , β , λ γ f t = 0 t e α , β γ λ , t τ f τ d τ , f t L 1 0 , 1
Moreover, P I α , β , λ γ f t , for f t L 1 0 , 1 can be presented as [51]
P I α , β , λ γ f t = k = 0 γ k k ! λ k I α k + β f t
That is P I α , β , λ γ f t is linear and bounded from L p 0 , 1 into L p 0 , 1 for any 1 p [51]. Recalling that for λ = 1 we get e 1 , t = t β 1 E α , β γ t α as CM function, if 0 < α γ β 1 as well as [41]
P I α , β , λ γ e α , μ ω λ t = e α , μ + β ω + γ λ t

2.2.2. Prabhakar Derivatives

The Prabhakar derivative in the Riemann–Liouville sense is defined as [41,49,51]
R L P D α , β , λ γ f t = d d t P I α , 1 β , λ γ f t
The Prabhakar derivative in the Caputo sense is defined as [41,49,51]
C P D α , β , λ γ f t = P I α , 1 β , λ γ d f t d t
In this context, for γ = 0 and λ = 0 the Prabhakar integral reduces to the Riemann–Liouville integral, [51].
P I α , β , 0 0 f t = 0 t e α , β γ 0 , t τ f τ d τ = 1 Γ β 0 t t τ β 1 f τ d τ = I β f t
The Laplace transform of C P D α , β , λ γ f t is [51]
L C P D α , β , λ γ f t = s β α γ s α λ γ L f t m = k = 0 1 s k 1 f k 0 +
With m defined as the integer part of β . For m = 0 (also for m = 1 ) we have [51]
L C P D α , β , λ γ t n = s β α γ s α λ γ Γ n + 1 s n + 1

2.2.3. Caputo Derivative

Considering the Caputo derivative [31,50]
C D t α = 1 Γ 1 α 0 t t τ α d f τ d τ d τ
and its Laplace transform
L C D t α t n = Γ n + 1 s n α + 1 , α > 0
we can see that for f t = t n , we get [51]
L C P D α , β , 0 γ t n = L C P D α , β , λ 0 f t = Γ n + 1 s n β + 1 = L C D β t n
That is, comparing the Laplace transforms (52) and (54), we see that the Caputo derivative (53) is a particular case of Prabhakar fractional derivative (49).

2.2.4. Atangana–Baleanu Derivative

Further, the Atangana–Baleanu derivative in the Caputo sense [52]
A B C D t α f t = M α 1 α 0 t E α α 1 α t τ α d f τ d τ d τ
With a normalization function M α such that M 0 = M 1 = 1 , and we may see that λ = α / 1 α , has a Laplace transform, when f t = t n , can be presented as [51]
L A B C D t α t n = M α 1 α 1 s α + α 1 α Γ n + 1 s n α + 1
Thus, we can see that A B C D t α f t is a particular case of C P D α , β , λ γ f t . Recall that E α , 1 1 z = e α , 1 1 z = k = 0 z k Γ α k + 1 (12).

2.2.5. Caputo–Fabrizio Derivative

Finally, let us consider the simplest case with the Caputo–Fabrizio derivative [53] (with M α = 1 [54]
C F D t α = 1 1 α 0 t exp α 1 α t τ d f τ d τ d τ
has a Laplace transform in case with f t = t n [51]
L C F D t α t n = Γ n + 1 s n 1 α + 1 α s
Then, considering L C P D α , β , λ γ t n for α = 1 , β = 0 , γ = 1 , and λ = α / 1 α we get [51].
L C P D 1 , 0 , α / 1 α 1 t n = Γ n + 1 s n 1 s + α 1 α = 1 α L C F D t α t n
That is, the Caputo–Fabrizio derivative is a particular case of the Prabhakar derivative of the Caputo-sense. Recall that exp z = E 1 , 1 1 z = e 1 , 1 1 z .
All these interrelationships between the derivatives of the Caputo sense naturally come from the interrelationships between the Mittag-Leffler function E α , β γ and the Prabhakar kernel e α , β γ when the parameters α , β and γ take different values as was discussed in Section 2.1.
In the sequel, we can see how these fractional derivatives logically appear through constitutive equations based on the fading memory formalism when the primary relaxation function (memory) is defined as the first derivative of a certain kernel based on the Mittag-Leffler type function.
After this brief exposition of interrelationships between fractional operators based on the Mittag-Leffler function, we highly recommend the analyses in [55,56,57,58,59] where more deep relations are developed. At this end, the naturally arising question is: Despite the mathematical correctness of the generalized operators, their applicability to real-world physical problems is still questionable, albeit some successful steps in this direction have been done [40,41,48,49,60,61,62].

2.3. The Fractional Order in Caputo–Fabrizio and Atangana–Baleanu Derivatives

Caputo–Fabrizio Fractional Operator: The Fractional Parameter

In the definition of the Caputo–Fabrizio operator, the stretched time t s is multiplied by a dimensional factor λ = α / 1 α which should have a dimension s 1 while actually, it should be dimensionless because physically α is a dimensionless parameter. By a nondimesionalization of the exponential function with the help of characteristic time of the relaxation process t 0 , namely
exp t s τ = exp t / t 0 s / t 0 τ / t 0 = exp t ¯ s ¯ τ ¯
where τ is the relaxation time and in the context used earlier τ = 1 / λ in the definition (58) (see also the explanations about the construction of (59) and has a dimension 1 / t i m e .
This nondimesionalization does not change the meaning of the exponential relaxation function but avoids any doubts about the definition of the fractional order α as [63].
1 α α = τ t 0 α = 1 1 + τ / t 0
The Atangana–Baleanu derivative of Caputo sense (ABC) with B α = 1 can be rescaled as [64]
A B C D a + α f t = 1 1 α 0 z d f s ¯ d s ¯ E α t ¯ s ¯ τ ¯ α d s ¯ , 1 α α = τ t 0 α = τ ¯ α , t ¯ = t t 0 , z = s / t 0
In detail, we have
A B C D a + α f t = 1 1 α 0 z d f s ¯ d s ¯ k = 0 1 Γ α k + 1 t ¯ s ¯ τ ¯ α k d s ¯
Therefore, the argument of the Mittag-Leffler kernel E α Z is Z = α / 1 α t s α and consequently we get [64]
1 α α k = τ t 0 k α = 1 1 + τ / t 0
The results are the same as established for the Caputo–Fabrizio operator [64,65,66].
Hence, the fractional order α introduced through λ = α / 1 α in the generalized formulation of the derivatives has a physical meaning and is strongly related to the characteristic time scales of the relaxation process modeled. In this context, the case with the Caputo–Fabrizio operator is more physically clear due to the large background of exponentially decaying processes.

3. Fading Memory Approach (Formalism) and Concept of Model Development

3.1. Fading Memory Concept

The Boltzmann linear superposition functional [16] for hereditary viscoelasticity with a time-dependent memory kernel (correlation function) R t , τ is
φ x , t = A 1 v x x , t + A 2 0 t R x , t τ v x τ d τ
The fading memory concept relating the flux to its gradient, for simple materials [67,68,69] is modeled by the following integrodifferential equation
j x , t = A 1 C x , t x A 2 t R t τ C x , τ x d τ
as a manifestation of the Boltzmann linear superposition functional 66). In (67) the transport coefficients A 1 and A 2 are diffusivities. The history value problem for (67) addresses the following integral [69]
d t = 0 R t τ C x , τ x d τ
allowing to give a function C x , t on < t 0 From (67) and (68) it follows that
x j x , t = A 1 2 C x , t x 2 A 2 0 t R t τ 2 C x , τ x 2 d τ + x d t
Since C x , t is a causal function (vanishing for t < 0 ) and considered only for 0 < t < we accept d t = 0 x d t = 0 , and therefore (69) can be rewritten as
x j x , t = A 1 C x , t x A 2 0 t R t τ C x , τ x d τ
Reducing (69) and setting the lower terminal of the memory integral to zero have deep physical meaning when applying hereditary integrals. The essence is well expressed by Hilfer in [70] that in fractional operators, the time is not the chronological time (instant time) but the intrinsic time of the process (the time of duration), starting at the point accepted as  t = 0 . Hence, in the context of the reduction of (69) to (70), there is no process before t = 0 .
The deep thermodynamic sense of the fading memory formulation is that the non-locality represented by the convolution term works for short times, while for long times, we get the first term in (70), i.e., the instant reaction of the system. Moreover, models constructed with the fading memory principle obey: the causality principle (through the convolution term) [26,27,28], thermodynamic consistency [71,72], and model observability (objectivity) [73,74,75].
Further, if A 1 = 0 , that is, the process modeled has no long-time asymptotically stable state, the result is a model close to the Continuum Time Random Walk (CTRW) without a stationary state [76,77], and applying the continuity equation to such flux, we get
C x , t t = x j x , t C x , t t = A 2 0 t R t τ C x , τ x d τ
To a greater extent, the model (71) mimics the master equation [77,78]; this is only a remark, and we will not elaborate on it further in this study but refer readers to [79,80,81] where this issue was developed in detail regarding Prabhakar-type memory.
Remark 5.
Fading memory concept in the case of heat conduction is well applicable to the so-called simple materials [21,82] where the flux is proportional to the temperature gradient (the term was coined after the work of Storm [82]), as in all examples studied here.

3.2. Model Build-Up Concept

Let us take a look at the flux definition of (67) where the relaxation (memory) function can be presented as M t = δ t + R t such that (67) expressed as a convolution product
j x , t = C x , t x M t = C x , t x A 1 δ t + A 2 R t
The two main requirements for R t are: it should be completely monotone, regardless of the fact that it is singular or non-singular, and that its first derivative d R t / d t = R t t , should also be completely monotone.
The two main requirements for R t are: it should be completely monotone, regardless of whether it is singular or non-singular at t 0 + , and that its first derivative d R t / d t = R t t should also be a completely monotone function.
Then, the main step in the construction of the flux constitutive model is the use R t t instead R t , that is
j x , t = C x , t x M t = C x , t x A 1 δ t + A 2 R t t
Let us consider the convolution term alone and use integration by parts to clarify things.
t R t t τ C x , τ x d τ = t C x , t x d R t τ = = C x , τ x d R t τ t t R t τ d d τ C x , τ x
The term = C x , t x d R t τ t is zero when we have: R t is a causal function and the initial conditions corresponding to a virgin material (medium) are C x , 0 = C 0 , 0 = C , t = C , t , C x x , 0 = C x x x , 0 = 0 (signaling problem). As a result, we get
t R t t τ C x , τ x d τ = t R t τ d d τ C x , τ x d τ
Denoting x C x , t = f x , t and setting the lower terminal at zero, we get a constriction of Caputo-type fractional derivative that is
C D t α f x , t = M α N α 0 t R t τ d f x , τ d τ d τ
where the pre-factor M α / N α that only takes into account fractional order, but it is determined by the type of memory kernel (see Section 2.2 and the developments in the sequel).
The flux can therefore be expressed as
j x , t = A 1 C x , t x A 2 C D t α C x , t x
In this case, the hereditary term is generically expressed as a fractional derivative of the Caputo type without any kernel-specific information. Then, applying the continuity equation C t = q x we get
C x , t t = A 1 C x , t x + A 2 C D t α C x , t x
For α = 1 , this model reduces to the local diffusion equation.
Remark 6.
Specific comment is required regarding the following exposition before moving on to the modeling method. Assuming that the kernel of the most straightforward scenario is built first, the model build-ups are demonstrated inductively. The model build-ups are illustrated deductively, which suggests that the kernel of the most straightforward scenario is constructed first. When the same approach is used to handle ever-more complex kernels, the resulting fractional operator structures gradually become more complex. The approach taken here is more instructive than the one frequently used in the literature, which starts with the most complex kernel and gradually creates simpler kernels by lowering the complexity parameters.

4. Heat Conduction Models with Mittag-Leffler-Type Memories: Examples

4.1. Heat Conduction with Infinite Flux Speed (Local in Time)

Here we start with simple heat conduction, local in time, which will allow us to see how the approach developed here upgrades it towards non-local versions expressed in terms of fractional operators. The energy conservation (continuity) equation is
ρ C p T t = q x
If the constitutive equation is (Fourier law)
q x , t = k T x , t x
we get the local in-time heat conduction model (Fourier model)
ρ C p T t = k 2 T x 2
or with constant density ρ = c o n s t . and heat capacity C p = c o n s t . as
T t = a 2 T x 2 , a = k ρ C p
where a = k ρ C p is the thermal diffusivity.

4.2. Heat Flux with Exponential Memory Kernel

Consider a virgin semi-infinite material (medium) with initial and boundary conditions (signaling problem)
T x , 0 = T , t = T , t = 0 , T x x , 0 = T x x x , 0 = 0 , T 0 , t = T s
These conditions will be valid for all examples developed in the sequel.
We start with the simplest case when the kernel is
e 1 , 1 1 λ t = exp λ t , d d t e 1 , 1 1 λ t = λ exp λ t
which will allow an easier understanding of the technology of the solution. This development follows the results of [66].
Hence, we consider an exponential memory (Jeffrey’s kernel) R t [83] with a finite relaxation time τ (and λ = 1 / τ )
R t s = exp t s τ , R t t = 1 τ exp t s τ
If we now apply the energy balance (79) (with ρ C p = c o n s t . ), the result is the Cattaneo equations [23] modeling only the relaxation part of the heat conduction (the RHS of (86))
T x , t t = a 2 τ 0 t exp t s τ T x , s x d s , a 2 = k 2 ρ C p
The Fourier law is the Cattaneo equation’s limit τ 0 resulting in (82)
Now, let us accept the construction of the relaxation kernel (73), namely
M t = k 1 δ t + k 2 R t t M t s = k 1 δ s + k 2 1 τ exp t s τ
In (87) δ s is the Dirac delta function, while a 1 and a 2 are the effective thermal conductivity and the elastic conductivity. Then, the heat flux can be expressed as
q x , t = k 1 T x , t x k 2 τ t e t s τ T x , s x d s
The energy balance (79) yields
T x , t t = a 1 2 T x , t x 2 + a 2 λ t e λ t s 2 T x , s x 2 d s , λ = 1 τ
Here a 1 = k 1 / ρ C p , a 2 = k 2 / ρ C p are thermal diffusivities: they have equal dimensions m 2 / s because the exponential second term in (87) is dimensionless.
Denoting F ( x , t ) = 2 T x , t / x 2 , for the sake of the simplicity of the expressions, and after integration by parts in the second term in (89) we get [66].
λ t e λ t s F x , s d s = e λ t s F x , s F x , t t + + λ t e λ t s F x , t F x , s d s
The first term in the RHS of (90), is zero, but the second one matches the definition of the Caputo–Fabrizio fractional derivative [53] (see also the general construction (75) of a Caputo-type derivative). Precisely, we have the formulation [53]
C F D t α = α 1 α 2 t f t f a s exp α 1 α t s d s , t > 0
This can be considered as a pro-Caputo (non-normalized) derivative denoted as P C D t λ . It can be expressed in two equivalent forms (following the notations used in [53] namely [66]
P C D t λ F x , t = λ t e λ t s F x , t F x , s d s = = λ t e λ t s d F x , s d t d s
where the rate constant λ 0 , controls the kernel. If we like to satisfy the conditions: α 0.1 1 / λ 0 , we get λ α = α / 1 α : the desired properties are obtained, namely
1 λ = 1 α α 0 , , α = 1 1 + 1 / λ = λ 1 + λ 0 , 1 , α 1 α 2 = λ 1 α
Further, following the definition of the Caputo–Fabrizio derivative [53] and considering the lower limit of integral at zero, we have [66]
C F D t α T x , t = N σ σ P C D t β T x , t = N σ σ 0 t e α 1 α t s d F x , s d t d s = = M α 1 α 0 t e α 1 α t s d F x , s d t
In the terms used here σ = 1 / λ = 1 α / α 0 , , while N α and M α are normalization functions [53] (see (58) and (76)).
Turning on (89) in terms of T x , t and taking into account the last expression of (94), we get [66]
T x , t t = a 1 2 T x , t x 2 + a 2 1 α C F D t α 2 T x , t x 2 , t > 0
Equation (95) models transient heat conduction with a damping term expressed through the Caputo–Fabrizio fractional derivative. For α = 1 we get the Fourier Equation (82).
If the Jeffrey kernel is only considered (in the present context (84)), for instance, then the equivalent form of (95) is (86) accounting only for the elastic part of heat diffusion and the result in terms of P C D t α is [66].
T x , t t = a 2 β 0 t exp β t s T x , s x d s T x , t t = a 2 α 1 α P C D t α T x , s x
As a result, we could see that the Caputo–Fabrizio derivative of a damping term was produced directly by the fading memory approach.
Remark 7.
It is worth noting and drawing attention to an incorrect interpretation of the Caputo–Fabrizio derivative based on its exponential kernel. As stated by Tarasov and Tarasova [84]: “Note that the memoryless property of the exponential distribution allows us to state that the differential operator (eq.6 in [84]) (i.e., the Caputo–Fabrizio derivative) cannot be used to describe processes with memory” (sic!).
As to the memoryless properties of the exponential function, this is a well-known fact from the students’ textbooks. However, the memory kernel should be a complete monotone, causal function, and thermodynamically consistent noting more. The exponential function obeys all of these requirements [14]. The memory properties of the non-local operator do not depend on whether the function used as a kernel has a memory or not, outside the context of the construction of the convolution integral.
The opinion about the Caputo–Fabrizio derivative in [84] could be interpreted either as a deep misunderstanding of the meaning of non-local operators or as a manipulative statement in the context of the confrontation [1,2,3,4,5,6,7] between some trends in fractional calculus (see the very beginning of the Introduction) rather than from a serious scientific perspective.
In addition, without continuing the discussion on the use of exponential memories, we refer to the book of Uchaikin [85] (in Section 2.1.3 Memory) where the same standpoint based on memoryless behavior of the exponential function is expressed.

4.3. Heat Flux with Mittag-Leffler (One-Parameter) Memory Kernel

Now, using the relation (26) we constitute a new relaxation (memory) function for the heat flux [86]
R M L = k 1 δ t + k 2 λ E α , 0 λ t α t = k 1 δ t + k 2 d d t E α , 1 λ t α
The second term in (97) is a monotonically decaying function, singular at t 0 , and thus, obeys all necessary properties to be used as a memory kernel. Then, the heat flux with memory can be expressed as [86]
q x , t = k 1 T x , t x k 2 λ t T x , s x d d s E α , 1 λ t s α d s
Now, the integration by parts of the integral (denoting F ( x , t ) = T x , t / x ) in (98) yields
t E α , 1 λ t s α F x , s d λ s = E α , 1 t s α F x , s F x , t t + + t E α , 1 λ t s α F x , s s d λ s
The first term in the RHS of (99) is zero, while the second one matches the definition of the Atangana–Baleanu derivative fractional derivative (56).
The coefficient λ , 0 can be mapped by (see the same relations in [66] (see (93) used earlier)
λ = α 1 α 0 , , α = 1 1 + 1 / λ = λ 1 + λ 0 , 1
Then, assuming the lower limit in the Stieltjes integral at zero (the causality principle), we have a new expression for the second term of the heat flux with a memory [86]
q e = a 2 0 t E α , 1 α 1 α t s α d d s T x , s x d s
We can recast the integral in (101) as the elastic component of the flux q e because it resembles how the ABC derivative is constructed. [86].
q e = a 2 1 α 1 1 α 0 t E α , 1 α 1 α t s α d d s T x , s x d s = = a 2 1 α A B C D t α T x , t x
Consequently
q e x = a 2 1 α 0 t E α , 1 α 1 α t s α d d s 2 T x , s x 2 d s = = a 2 1 α A B C D t α 2 T x , t x 2
Finally, the energy conservation equation yields [86]
T x , t t = a 1 2 T x , t x 2 + a 2 1 α A B C D t α 2 T x , t x 2
For α = 1 , it reduces to the Fourier equation.
Based on the fading memory formalisms for the heat flux, we found that the ABC derivative spontaneously manifests itself when a suitable memory kernel (97) is incorporated into the constitutive equation.

4.4. Heat Flux with Prabhakar Memory Kernel

After the two prior examples, which used relatively simple kernels, we now attempt to create a fractional derivative using the Prabhakar kernel and the fading memory formalism. Therefore, applying (36) and its first derivative (43) we get for the heat flux with relaxation function defined as
R P B = k 1 δ t + k 2 d d t e α , β γ λ t α = k 1 δ t + k 2 t β 2 E α , β 1 γ λ t α
Then, the heat flux can be formulated as
q x , t = k 1 T x , t x k 2 t T x , u x d d u e α , β γ t , λ u α d u
Assuming for the sake of simplicity λ = 1 , and applying integration by parts in the second term of (106), and what mutatis mutandis, we get
q x , t = k 1 T x , t x k 2 0 t e α , β γ t , u α d d u T x , u x d u
In general, the second term in (107) mimics the fractional derivative of the Caputo type.
Alternately, by using the Laplace transform on (106) and accounting for the convolution product in the second term, we obtain (see (38))
L q x , t = L k 1 T x , t x k 2 L d d t e α , β γ λ t α × L T x , t x
that is
q x , s = k 1 T x , s x k 2 s α γ β + 1 s + λ γ × T x , s x = = k 1 T x , s x k 2 s α γ β s + λ γ × s T x , s x
The inverse transformation L 1 q x , s yields
q x , t = k 1 T x , t x k 2 0 t e α , β γ t , λ u α d d u T x , u x d u
The second term in (110) has the construction of Caputo-type with the Prabhakar memory kernel (see (49)).
Assuming a normalization function N α = 1 α and setting λ = α / 1 α (110) can be reformulated as follows without loss of generality
q x , t = k 1 T x , t x k 2 1 α 1 1 α 0 t e α , β γ t , λ u α d d u T x , u x d u λ = α 1 α , t > 0
Now, applying the energy balance equation, we have
T x , t t = a 1 T x , t x + a 2 1 α 1 1 α 0 t e α , β γ t , λ u α d d u T x , u x d u , λ = α 1 α , t > 0
For α = 1 , the non-locality in (112) is absent, and it reduces to the Fourier model.
The construction
H C P D α , β γ f ( t ) = 1 1 α 0 t e α , β γ t , λ u α d d u f u d u , λ = α 1 α , 0 < α < 1 , t > 0
is a fractional operator (derivative) that uses a base point zero but does not exactly match the definition (49). The differences result from the fact that (49) is formulated constitutively, mathematically correct (but with no physics behind it), and similar to the traditional Caputo derivative. They also result from the definition of the Prabhakar fractional integral (45) (in a Riemann–Liouville sense) [38,41,43,45,49].
Hence, we may present (112) as
T x , t t = a 1 T x , t x + a 2 1 α H C P D α , β γ T ( x , t ) x
It reduces to the Fourier model for α = 1 .
The Laplace transform of H C P D α , β γ f ( t ) , with the assumption k = 0 1 s k 1 f k 0 + = 0 , for the sake of simplicity, yields
L H C P D α , β γ f t L e α , β γ λ t α × L f t s β α γ s + λ γ × f s
Comparing to (38) (see also (52)), we see that the Laplace transform is the same, which is a nice result.
Furthermore, as it was shown in Section 2.2, it is easy to see, and straightforwardly developed, that for some particular cases, the kernel of H C P D α , β γ f ( t ) reduces to memory functions of well-known fractional operators. For instance, we obtain the Atangana–Baleanu derivative for γ = β = 1 .
Additionally, it is evident that for γ = 1 we have e α , β , λ 1 λ t α , a kernel based on the two-parameter Mittag-Leffler function
e α , β , λ 1 = t β 1 E α , β λ t α
This enables the definition of a novel derivative (which has not yet been studied), namely
H C P D α , β 1 f t = 1 1 α 0 t t β 1 E α , β λ t α , τ d f τ d τ d τ
If λ = α / 1 α , it can be aligned to the ABC definition and reduced to it for β = 1 .

5. Experiments with Other Known Kernels

In this section, we clearly show how well-known memory kernels related to the Mittag-Leffler functions can be incorporated into constructions of fractional operators of Caputo-type.

5.1. Rzanitsyn Kernel

This is a composite function defined as [87,88,89] (singular at t 0 + )
M R α , t ; λ = t α 1 Γ α exp λ t , 0 < α < 1
which is a product of a particular case, the Sonine kernel t α 1 Γ α (known as the kernel of fractional integration or Gel’fand-Shilov distribution) [90,91] and the exponential decaying function.
It can be presented as a product of two particular versions of the Prabhakar kernel (see (44))
M R α , t ; λ = e α , β γ ( λ = 0 , t ) × e 1 , 1 1 ( λ t )
or in terms of the two-parameter Mittag-Leffler function as
M R = t α 1 Γ α E 1 , 1 1 λ t
The Laplace transform of M R is (see Equation (25) for the special case α = β = 1 )
L M R α , λ t = 1 Γ α 1 s + λ
The derivative d d t M R α , t ; λ is (recall the definition (25) with α = β = 1 ) is
d d t M R α , t ; λ = t 1 Γ α E 1 , 0 α , t ; λ
which is singular at t 0 + with a Laplace transform (see (121))
L d d t M R α , λ , t = 1 Γ α s s + λ
Assuming λ = α 1 α , we get
M R α , t = t α 1 Γ α E 1 , 1 1 α 1 α t , 0 < α < 1
d d t M R α , t = t 1 Γ α E 1 , 0 α , t ; α 1 α , 0 < α < 1
Now, we may write the heat flux constitutive equation as
q R x , t = k 0 T x , t x k 2 0 t t τ 1 Γ α E 1 , 0 α , t ; α 1 α T x , t x d τ
Then, what is mutatis mutandis, we get
q R x , t = k 0 T x , t x k 2 0 t M R α , t ; α 1 α d d τ T x , τ x d τ
and a heat diffusion equation
T x , t t = a 0 T x , t x + a 2 1 α 0 t M R α , t ; α 1 α d d τ T x , τ x d τ
With the Rzhanitsyn kernel, we, therefore, defined a fractional derivative of the Caputo type while assuming a normalization function N α = 1 α , namely
R C D t α f t = 1 1 α 0 t M R α , t ; α 1 α d d τ f τ d τ

5.2. Miller–Ross Kernel

The Miller–Ross function (MR) function is defined as
M M R α , λ t = t α k = 0 λ t k Γ α + k + 1 = k = 0 λ k t k + α Γ α + k + 1 t α
It is defined as the D t α e λ t = I t α e λ t as the α t h integral of the exponential function e λ t
It can be represented as the two-parameter Mittag-Leffler function by [31]
M M R α , λ t = t α E 1 , α + 1 λ t , 0 < α < 1
The Laplace transform of the MR function is
L M M R α , λ t α = 1 s α 1 1 λ s 1 , λ s < 1
The time derivative of the MR function is
d d t M M R α , λ t = d d t k = 0 λ k t k + α Γ α + k + 1 = t α 1 k = 0 k + α λ k t k + α Γ α + k + 1
and through the Laplace transform, with the rule L d d t f t = s L f t f 0 = s f s f 0
L d d t M M R α , λ t α = s 1 s α 1 1 + λ s 1 = s 2 α 1 λ s 1 , λ s < 1
because M M R α , λ t α t = 0 = 0 .
Then, if we define a relaxation function d d t M M R α , λ t = d d t t α E 1 , α + 1 λ t = R M R α , λ t such that R M R α , λ t and M M R α , λ t are monotonically decaying functions, we may write the heat flux constitutive equation as
q M R x , t = k 0 T x , t x k 2 0 t t τ α E 1 , α + 1 λ t τ α T x , τ x d τ
Then, following the integration by parts in the second term of (135), and applying all steps as in the preceding examples, as well as faster by the Laplace transform of q M R x , t
q M R x , s = k 0 T x , s x k 2 0 t s 1 s α 1 1 + λ s 1 T x , s x d s , 0 < α < 1
Then, a rearrangement of the second term in (136) yields
0 t s 1 s α 1 1 λ s 1 T x , s x d s = 0 t 1 s α 1 1 + λ s 1 s T x , s x
remembering that T x , 0 x = 0 , as the initial condition defined at the beginning. This allows expressing the hereditary term with the MR function as a memory kernel, namely
q M R x , t = k 0 T x , t x k 2 0 t M M R λ t τ α T x , t x d τ , 0 < α < 1 , t > 0
Now, without loss of generality, we may assume λ = α 1 α , as well as defining the normalization function N α = 1 α . Thus, we may formulate a Miller–Ross fractional operator of Caputo-type, namely
M R C D t α f t = 1 1 α 0 t M M R α 1 α t τ α d d τ f τ d τ .
Applying the energy conservation equation with the flux expressed as (138) we get
T x , t t = a 0 2 T x , t x 2 + 1 α a 2 M R C D t α T x , t x
For α = 1 , Equation (140) reduces to the local Fourier equation.
To further explain the recently developed result, it is important to draw your attention to the fact that a fractional integral with MR function is defined in [92] as (in a Riemann–Liouville style) (see also [93,94])
M R I 0 + α , λ f t = 0 t M R λ t τ α f τ d τ
and based on it, a fractional Liouville–Sonine derivative was proposed as [92]
L S M R D t α f t = d d t M R I 0 + α , λ f t = d d t 0 t M R λ t τ α f τ d τ

5.3. Rabotnov Kernel

The Rabotnov function is defined originally as [89,95]
M R a b α , λ t = t α 1 1 λ α k = 0 1 k λ t α k Γ α k + 1 , λ = 1 τ k , 1 < α < 0
and can be presented alternatively as [31]
M R a b α , λ = t α k = 0 λ k t k α + 1 Γ k + 1 α + 1 ] , 1 < α < 0
or in terms of two-parameter Mittag-Leffler functions as [31]
M R a b α , λ = t α E α + 1 , α + 1 λ t α + 1 , M R a b α , λ = t α 1 E α , α λ t α
For α = 1 the function M R a b 1 , λ reduces to the exponential function. However, when α 0 it goes to 1 2 δ t [89].
Its time derivative, applying (25) to the second version of (145), can be expressed as
d d t M R a b α , λ = t α 2 E α , α λ t α
The original Rabotnov function (143) is a growing in time function, conceived to model creep (deformation, extension) of viscoelastic and elastoplastic materials [89], and in such a case, the kernel in a hereditary integral should be
M R a b α , λ = t τ α k = 0 λ k t τ k 1 α Γ k + 1 1 α
A decaying (monotonic) version was used in a new fractional operator known as the Yang–Abdel-Aty–Cattani (YAC) derivative [92,96,97] (see also [98,99,100])
M Y A C λ , t α = k = 0 λ k t k + 1 α + 1 1 Γ k + 1 α + 1
Its Laplace transform is [98,99]
L M R a b λ t α = 1 s α + 1 1 1 λ s α + 1 , λ s α + 1 < 1
The time derivative of M R a b λ t α through the Laplace transform is
L d d t M R a b λ t α = s 1 s α + 1 1 1 λ s α + 1 = 1 s α 1 1 λ s α + 1 , λ s α + 1 < 1
Now, constructing the heat flux constitutive equation, we get
q Y A C x , t = k 1 T x , t x k 2 0 t d d t M R a b λ t α T x , τ x d τ
Then, applying the Laplace transform to both sides of (151) and taking into account (149), we have
q Y A C x , s = k 1 T x , s x k 2 0 t 1 s α + 1 1 1 λ s α + 1 s T x , s x d s
The inverse Laplace transform yields
q Y A C x , t = k 1 T x , t x k 2 0 t M R a b λ t τ α d d t T x , τ x d s
Consequently, the application of the energy balance equation results in
T x , t t = a 1 2 T x , t x 2 + a 2 0 t M R a b α 1 α t τ α d d t 2 T x , τ x 2 d τ
We may suggest again that λ = α 1 α and normalization function N α = 1 α . The result is a heat conduction equation
T x , t t = a 1 2 T x , t x 2 + a 2 1 α C R a b D t α 2 T x , t x 2
with damping terms expressed through fractional of Caputo-type with Rabotnov kernel, defined as
C R a b D t α f t = 1 1 α 0 t M R a b α 1 α t τ α d d t f τ d τ
This definition differs from the YAC derivative, where the approach uses a fractional integral of a Riemann–Liouville type [92,96,97]
R a b I t α f t = 0 t M R a b λ t τ α f τ d τ , f t L 0 ,
and a left-sided fractional derivative in the classical (termed Liouville–Caputo derivative) style is declared, without any derivation from a physical model, namely
R a b D t α f t = t M R a b λ t τ α f 1 τ d τ
and without any specification of the parameter λ .
Remark 8.
This section of the study can be summarized by pointing out that the new derivatives naturally appear when the constitutive equation for the heat flux is constructed using a proper Mittag-Leffler memory kernel and the fading memory formalism. The integer-order time derivative of the energy balance equation is unaffected by the new constitutive equation, and the fractional derivatives appear as damping terms.

6. Heat Conduction Model Analysis

Now, after the development of six models with various time-fractional derivatives as damping terms, we have to discuss:
What are the differences between the models using some of the fractional operators discussed in this article and those published in the literature?

6.1. A Trivial Example That Is Correct but Might Be Misleading

At the beginning of this discussion, we like to stress the attention on a trivial for the fractional calculus example. Let us consider again Equation (6), assuming that k x = k = c o n s t . (time and space independent thermal conductivity, for the sake of simplicity) precisely the first version of its RHS, namely
q x , t = k C D t α T x , t x
and represent it through its corresponding fractional integral of order 1 α as
q x , t = k I t 1 α T x , t x
This is a completely casual relationship, in a convolution form, due to the time shift and history (memory) of the gradient related to the time evolution of the flux. In addition, this formulation does not consider a long time term (that is, thermal conductivity k 1 = 0 , while k 2 0 ). Moreover, taking into account that C D t γ f = C I γ f and with γ = 1 α , we get an alternative form in terms of the left Caputo derivative, namely
q x , t = C D t α 1 k T x , t x
Taking into account that C D t γ f = C I γ f and with γ = 1 α we got (161) Further, since C D t γ is an operator left inverse of the fractional integral, we have C D t γ I γ = C D t γ C D t γ f = D t 0 f = f .
Thus, the application of the continuity equation yields
T x , t t = k I t 1 α 2 T x , τ x 2 x , t = k C D t α 1 2 x , t x 2
Applying to both sides of (162) the Caputo operator of order α 1 we get
α T x , t t α = a 2 T x , t x 2 , a = k ρ C p
This is the well-known time-fractional diffusion equation. It was correctly derived from a causal relationship (160). It is correct from both the mathematical and thermodynamic points of view.
However, when mathematical modeling skills are lacking, this simple equation may be misleading because it can be formalistically written based on the normal diffusion equation (when α = 1 ) by a simple fractional replacement of the integer-order time derivative. This is the well-known “replacement fractionalization,” which is generally an incorrect operation. To be precise, such a replacement means that the continuity equation is written as
α T x , t t α = k q x , t x T x , t = k I α q x , t x
which contradicts the First Law of Thermodynamics.
In the case presented here, the manipulations were correct mathematically, but the replacement is incorrect, and the starting point should be (160) as an element of the fading memory constitutive equation.
In a general sense, when a long time behavior is not considered (when the memory kernel goes to zero) as a model of the process (see the comments on CTRW models where stationary states are missing in the concept), we have to use the general constitutive equation
u x , t t = k j x , t x = x k 0 t R α , t τ d d τ u x , τ x d τ
The more general flux formulation with k 1 0 , when the long time behaviour should be taken into account, is
j = k 1 u x , t x k 2 0 t R α , t τ d d τ u x , τ x d τ
Then, after the application of the continuity equation, we have
u x , t t = k 1 2 u x , t x 2 + k 2 0 t R α , t τ d d τ 2 u x , τ x 2 d τ
In (165)–(167) u ( x , t ) could be concentration, temperature, or another variable (such as stress or deformation), and the non-locality is modeled by the memory kernel (relaxation function, also termed a correlation function), and it is strongly dependent on the physics of the modeled process. It would be a nice situation if the resulting fractional operator had semigroup properties, but this is not possible in some cases. The most important issue is the correct correspondence of the memory kernel to the physical relaxation process [64,65].
In the following section, we will examine some existing diffusion models from the literature to determine whether or not they confirm the results of the fading memory formalism.

6.2. Heat Conduction and Diffusion Models with Mittag-Leffler Kernels: An Analysis

Garra–Garrappa’s Non-Linear Heat Conduction Model

In the general context of this study, we address a model of a non-linear heat equation formulated as (Section 4 in [43])
D α , β , λ γ T x , t = K T 2 T x , t x 2 β T , K T = k 0 T ξ , ξ > 0 , t > 0 , x R
The derivative D α , β , λ γ was applied in two versions: a non-regularized (Riemann–Liouville type) derivative (169), and a regularized (Caputo-type) derivative (170), defined as (in the original notations established in [101]) (see also the definitions in Section 2.2.2),
0 D t α + λ γ f t d m d t m 0 t e α , α γ γ t u ; λ f u d u
C 0 D t α + λ γ f t 0 t e α , α γ γ t u ; λ f m u d u
with a Prabhakar integral (see also the definition (45))
0 I t α + λ γ f t e α , α γ γ t ; λ f t = 0 t e α , α γ γ t u ; λ f u d u
In general, this type of non-linearity in the diffusion term, when the equation is local in time, makes the heat conduction model a degenerate parabolic equation. It is well-known, since the first solutions in the 50s of the last century, that such degenerate parabolic problems have solutions moving with finite speeds, in contrast to the case with ξ = 0 where the flux speed is infinite (see [102] and the reference therein).
Now, for the sake of the coherence of the exposition, we return to the common symbol used here (let the definitions (169)–(171) remain, demonstrating that different versions of the symbols used are also possible; in general, there are various notations of these operators, so the reader should carefully understand the meaning in each particular case).
If we like to construct (derive) such a heat conduction model, we may define the flux-gradient relationship as
q x , t = 0 t e α , β γ λ , t τ k 0 T ξ T x , t x = = k 0 1 + ξ 0 t e α , β γ λ , t τ T ξ + 1 x , t x
If the energy balance equation is (1) is applied (neglecting the sink term for the sake of simplicity), then the heat conduction model would be
T x , t t = k 0 1 + ξ 0 t e α , β γ λ , t τ 2 T ξ + 1 x , τ x 2
or in a more compact form as
T x , t t = k 0 1 + ξ P I α , β γ 2 T ξ + 1 x , t x 2
In the local case when memory is δ t it will be
T x , t t = k 0 1 + ξ 2 T ξ + 1 x , t x 2
and the solutions are mainly approximate (see [102]) with convex temperature profiles (with almost infinite gradient close to the solution front for high values of ξ ), in contrast to the case for ξ = 0 with concave temperature distributions.
Now, the principle question is: how to transform (173) to (168)? Taking into account that the regularized Caputo-type derivative C D α , β , λ γ acts as left-inverse of the Prabhakar fractional integral [41,49], as well as that C D α , β , λ γ f t = R L C D α , β , λ γ f t at zero initial conditions (see eq. 5.11 in [41]), we may apply C D α , β , λ γ to both sides of (173) or (174). The result of this operation is
C D α , β , λ γ T x , t t = k 0 1 + ξ 2 T ξ + 1 x , t x 2
We can see that the left side of (176) is not the same as in (168). That is, the operation that was successful with the Caputo derivative in Section 6.1 does not work here. The origin is that (173) is a natural consequence of the constructive flux equation, while (168) is postulated and could be explained if the continuity equation is defined as
C D α , β , λ γ T x , t = q x , t x
that is an imitation of the First Law of Thermodynamics, where the causality principle is violated. Despite this, Equation (168) was solved elegantly even though a physical analysis on its basis is impossible.
Remark 9.
We can see that when the models are developed by the “copy-paste” technology, known as “replacement” in fractional modeling, the outcomes are non-physical. Especially for the Prabhakar derivative, it is important to stress the fact that all efforts were oriented towards “framing” in Kochubei’s “general fractional calculus”, [49], thus, attempting to make it “unipolar fractional calculus”. High merits, with good results, but physically inapplicable because Kochubei’s calculus is almost “sterile”, and not related to solutions of physical problems: Bright mathematics where the results are not provoked by nature and the consequence in fractional modeling is evident: nice mathematics, as it was done in [43], but with inexplicable results. We may consider all these studies, and others of the same style (products of replacement fractionalizations), as mathematical experiments, with elegant mathematics and still expecting relevant physical interpretations.
Remark 10.
Despite the strong verdict in the previous remarks, there are some good cases when the results are adequate and provide new information, such as modeling relaxations in anomalous dielectrics [60,101], diffusion with non-static stochastic resetting [103], generalized Langevin equation [104], anomalous diffusion [105], even though the models already exist, just the memory kernels are changed.
After this analysis, we will skip similar models with Miller–Ross [94] and Rabotnov derivatives [96,97,100], where the models are obtained by replacements.

7. Fading Memory Approach or Volterra Equations?

The results developed the fading memory formalism and the construction of the Stieltjes integral
0 t R t τ d f τ d τ d τ
lead to the formulation of normalized fractional operators of Caputo type, following the definition of Stieltjes integral that f τ is a function of bounded variation.
However, if the constitutive equations are of the Volterra type, such as (10), the hereditary term is
0 t R t τ f τ d τ
The construction (179) allows formulations of both fractional integrals and fractional derivatives of a Riemann–Liouville type (as in the studies devoted and analyzes here on the Prabhakar calculus)
I t α f τ 0 t R t τ f τ d τ
D t α f τ d m d t m 0 t R t τ f τ d τ , m = 1 , 2 , 3
The kernel R ( t ) in both formulations depends on the physics of the modeled relaxation process.
Taking into account that the “driving force,” i.e., the reason the heat flux to “flow” is the temperature gradient T / x (for simple materials [18,19,20,21,67,69,71]), then the construction (178) of the hereditary term is physically motivated. Hence, the physics of the process motivates the model, not the formal and voluntary use of fractional operators.
The application of Volterra equations in the process of generation of fractional operators with different kernels is beyond the scope of this study, but we refer to [89], where there are many examples of this approach applied in the early studies on fractional modeling in viscoelasticity.

8. Final Comments and Outcomes

Now, after this long study with many analyses and remarks, the main questions are
  • What the main achievements are?
  • How the systematic approach applying the fading memory formalism allows constructions of physically adequate fractional models?
Since developing this study, many explanations and specific remarks were made; the answers to these questions can be briefly outlined as
1. One of the principal achievements of this study is that it shows how different fractional operators appear in physical models, avoiding unmotivated replacement fractionalizations. The emphasis was on different versions of the Mittag-Leffler functions as memory kernels.
2. We demonstrated only six examples, but there is no limit to developing new ones if the constitutive equations are properly formulated.
3. We emphasize the importance of correct formulations of constitutive equations because they express physical laws relating to causes and effects in mathematical forms.
4. Therefore, the first step is to correctly formulate a constitutive equation with a hereditary term involving a memory kernel connected to the underlying physics of the modeled process. Then, the balance equation yields the final model (to which various initial and boundary conditions can be applied). However, the final models are not direct expressions of physical laws where hereditary terms can be inserted by replacements.
5. Some comments on existing heat conduction models confirm the standpoint of point 4. There are many articles with questionable results (see, for instance, the comments of Hilfer in [10] on physics of fractional models) where the final model is fractionalized by replacement, and the results that cannot be interpreted because the basic physics of the modeled process is violated.
We believe that this study will provide necessary and enough instructive information on how fractional modeling should be done. There is no need for high mathematics to formulate the constitutive equation correctly with a certain memory term and to get a correct model. High mathematics applied to incorrect models does not provide useful information and could be related to mathematical experiments, not to mathematical modeling of real-world phenomena.
The example of heat conduction used here is merely one to illustrate the general idea behind the physical basis of the formulation of fractional operators. This technique can be applied to a wide variety of other physical processes that exhibit memories and have properly formulated constitutive equations.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

There are no conflict of interests.

References

  1. Giusti, A. A comment on some new definitions of fractional derivative. Nonlinear Dyn. 2018, 93, 1757–1763. [Google Scholar] [CrossRef] [Green Version]
  2. Stynes, M. Fractional-order derivatives defined by continuous kernels are too restrictive. Appl. Math. Lett. 2018, 85, 22–26. [Google Scholar] [CrossRef]
  3. Ortigueira, M.D.; Machado, J.A.T. A critical analysis of the Caputo–Fabrizio operator. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 608–611. [Google Scholar] [CrossRef]
  4. Hilfer, R.; Luchko, Y. Desiderata for fractional derivatives and integrals. Mathematic 2019, 7, 149. [Google Scholar] [CrossRef] [Green Version]
  5. Diethelm, K.; Garrappa, R.; Giusti, A.; Stynes, M. Why fractional derivatives with nonsingular kernels should not be used. Fract. Calc. Appl. Anal. 2020, 23, 610–634. [Google Scholar] [CrossRef]
  6. Hanyga, A. A comment on a controversial issue: A generalized fractional derivative cannot have a regular kernel. Fract. Calc. Appl. Anal. 2020, 23, 211–223. [Google Scholar] [CrossRef] [Green Version]
  7. Diethelm, K.; Kiryakova, V.; Luchko, Y.; Machado, J.A.T.; Tarasov, V.A. Trends, directions for further research, and some open problems of fractional calculus. Nonlinear Dyn. 2022, 107, 3245–3270. [Google Scholar] [CrossRef]
  8. Baleanu, R.P.; Agarwal, R.P. Fractional calculus in the sky. Adv. Differ. Equ. 2021, 2021, 117. [Google Scholar] [CrossRef]
  9. Caputo, M.; Fabrizio, M. On the singular kernels for fractional derivatives. Some applications to partial differential equations. Prog. Fract. Differ. Appl. 2021, 7, 1–4. [Google Scholar]
  10. Hilfer, R. Mathematical and physical interpretations of fractional derivatives and integrals. In Handbook of Fractional Calculus with Applications; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019; Volume 1, pp. 47–85. [Google Scholar]
  11. Tateishi, A.A.; Ribeiro, H.V.; Sandev, T.; Petreska, I.; Lenzi, E.K. Quenched and annealed disorder mechanisms in comb models with fractional operators. Phys. Rev. E 2020, 101, 022135. [Google Scholar] [CrossRef] [Green Version]
  12. Amendola, G.; Fabrizio, M.; Golden, M. Thermodynamics of Materials with Memory; Springer: New York, NY, USA, 1964. [Google Scholar]
  13. Carillo, S. Some remarks on materials with memory: Heat conduction and viscoelasticity. J. Nonlinear Math. Phys. 2005, 12, 163–178. [Google Scholar] [CrossRef] [Green Version]
  14. Fabrizio, M.; Gentili, G.; Reynolds, D.W. On rigid heat conductors with memory. Ind. J. Eng. Sci. 1998, 36, 765–782. [Google Scholar] [CrossRef]
  15. Fabrizio, M. Fractional rheological models for thermomechanical systems. Dissipation and free energies. Frac. Calc. Appl. Anal. 2014, 17, 206–223. [Google Scholar] [CrossRef]
  16. Boltzmann, L. Zur Theorie der Elastischen Nachwirkung. Sitzungsber. Akad. Wiss. Wien. Mathem.-Naturwiss 1874, 70, 275–300. [Google Scholar]
  17. Volterra, V. Lecons sur la Theorie Mathematique de la Lutte Pour la Vie; Gauthier-Villars: Paris, France, 1931. [Google Scholar]
  18. Coleman, B.D.; Mizels, V.J. Thermodynamics and departure from Fourier’s law of heat conduction. Arch. Ration. Mech. Anal. 1963, 13, 245–261. [Google Scholar] [CrossRef]
  19. Gurtin, M.E. On the thermodynamics of materials with memory. Arch. Ration. Mech. Anal. 1968, 28, 40–50. [Google Scholar] [CrossRef]
  20. Nunziato, J.W. On heat conduction in materials with memory. Quart. J. Appl. Math. 1971, 29, 187–204. [Google Scholar] [CrossRef] [Green Version]
  21. Day, W.A. The Thermodynamics of Simple Materials with Fading Memory; Springer: Berlin/Heidelberg, Germany, 1972. [Google Scholar]
  22. Giorgi, C.; Gentili, G. Thermodynamic properties and stability for heat flux equation with linear memory. Quart. J. Appl. Math. 1993, 51, 343–362. [Google Scholar] [CrossRef] [Green Version]
  23. Cattaneo, C. On the conduction of heat. Atti Sem. Mat. Fis. Univ. Modena 1948, 3, 83–101. (In Italian) [Google Scholar]
  24. Cattaneo, C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation. Comp. Rend. Hebd. Séances Acad. Sci. Paris 1958, 247, 431–433. (In French) [Google Scholar]
  25. Vernotte, P. Paradoxes in the continuous theory of the heat equation. Comp. Rend. Hebd. Séances Acad. Sci. Paris 1958, 246, 3154–3155. (In French) [Google Scholar]
  26. Mittelstaedt, P.; Weingartner, P.A. Laws of Nature; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  27. Nussenzveig, H. Causality and Dispersion Relations; Academic Press: Cambridge, MA, USA, 1972. [Google Scholar]
  28. Lighthill, M.J. An Introduction to Fourier Analysis and Generalized Functions; Cambridge University Press: London, UK; New York, NY, USA, 1959. [Google Scholar]
  29. Fabrizio, M.; Giorgi, C.; Morro, A. Modeling of heat conduction via fractional derivatives. Heat Mass Transfer 2017, 53, 2785–2797. [Google Scholar] [CrossRef]
  30. Mainardi, F. Why the Mittag-Leffler function can be considered the queen of the fractional calculus? Entropy 2020, 22, 1359. [Google Scholar] [CrossRef] [PubMed]
  31. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA; Boston, MA, USA; New York, NY, USA; London, UK; Tokyo, Japan; Toronto, ON, Canada, 1999. [Google Scholar]
  32. Polard, H. The completely monotonic character of the Mittag-Leffler function. Bull. Am. Math. Soc. 1948, 52, 908–910. [Google Scholar]
  33. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  34. Miller, K.S.; Samko, S.G. A note on the complete monotonicity of the generalized Mittag-Leffler function. Real Anal. Exch. 1997, 23, 753–755. [Google Scholar] [CrossRef]
  35. Miller, K.S.; Samko, S.G. Completely monotonic functions. Integral Transform. Spec. Funct. 2001, 12, 389–402. [Google Scholar] [CrossRef]
  36. Schneider, W.R. Completely monotone generalized Mittag-Leffler functions. Expo. Math. Soc. 1996, 14, 3–16. [Google Scholar]
  37. Prabhakar, T.R. A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  38. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Generalized Mittag-Leffler function and generalized fractional calculus operators. Integral Trasform. Spec. Funct. 2004, 15, 31–49. [Google Scholar] [CrossRef]
  39. Górska, K.; Horzela, A.; Lattanzi, A.; Pogány, T.K. On complete monotonicity of three parameter Mittag-Leffler function. Appl. Anal. Discret. Math. 2021, 5, 118–128. [Google Scholar] [CrossRef]
  40. Mainardi, F.; Garrappa, R. On complete monotonicity of the Prabhakar function and non-Debye relaxation in dielectrics. J. Comp. Phys. 2015, 293, 70–80. [Google Scholar] [CrossRef] [Green Version]
  41. Giusti, A.; Colombaro, I.; Garra, R.; Garrappa, R.; Polito, F.; Popolizio, M.; Mainardi, F. A practical guide to Prabhakar fractional calculus. Frac. Calc. Appl. Anal. 2020, 23, 9–54. [Google Scholar] [CrossRef] [Green Version]
  42. Tomovski, Z.; Pogany, T.K.; Srivastava, H.M. Laplace type integral expressions for a certain three-parameter family of generalized Mittag-Leffler functions with applications involving complete monotonicity. J. Frankl. Inst. 2014, 351, 5437–5454. [Google Scholar] [CrossRef]
  43. Garra, R.; Garrappa, R. The Prabhakar of three parameter Mittag-Leffler function: Theory and application. Commun. Nonlinear Sci. Numer. Simulat. 2018, 56, 314–329. [Google Scholar] [CrossRef] [Green Version]
  44. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  45. Garra, R.; Gorenglo, R.; Polito, F.; Tomovski, Z. Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef] [Green Version]
  46. Mainardi, F.; Tomovski, Z. Some properties of Prabhakar-type fractional calculus operators. Fract. Differ. Calc. 2016, 6, 73–94. [Google Scholar]
  47. Srivastava, H.M.; Tomovski, Z. Fractional calculus with an integral operator containing a generalized Mittag-Leffler function in the kernel. Appl. Math. Comp. 2009, 211, 198–210. [Google Scholar] [CrossRef]
  48. Giusti, A.; Colombaro, I. Prabhakar-like fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 138–143. [Google Scholar] [CrossRef] [Green Version]
  49. Giusti, A. General fractional calculus and Prabhakar’ theory. Commun. Nonlinear Sci. Numer. Simul. 2020, 83, 1055114. [Google Scholar] [CrossRef] [Green Version]
  50. Caputo, M. Elasticita e Dissipazione, 1st ed.; Zanichelli Bologna: Bologna, Italy, 1969. [Google Scholar]
  51. Area, I.; Nieto, J.J. Fractional-order logistic differential equation with Mittag-Leffler-type kernel. Fractal Fract. 2021, 5, 273. [Google Scholar] [CrossRef]
  52. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  53. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  54. Losada, J.; Nieto, J.J. Properties of a New Fractional Derivative without Singular Kernel. Prog. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
  55. Abdeljawad, T.; Baleanu, D. Integration by parts and its applications to a new nonlocal fractional derivative with Mittag-Leffler nonsingular kernel. J. Nonlinear Sci. Appl. 2017, 10, 1098–1107. [Google Scholar] [CrossRef] [Green Version]
  56. Abdeljawad, T.; Baleanu, D. On fractional derivatives with generalized Mittag-Leffler kernels. Adv. Differ. Equ. 2018, 2018, 468. [Google Scholar] [CrossRef] [Green Version]
  57. Abdeljawad, T. Fractional operators with generalized Mittag-Leffler kernels and their iterated differintegrals. Chaos 2019, 29, 023102. [Google Scholar] [CrossRef]
  58. Fernandez, A.; Baleanu, D.; Srivastava, H.M. Series representation f0r fractional-calculus operators involving generalized Mittag-Leffler functions. Commun. Nonlinear Sci. Appl. Num. Simul. 2019, 67, 157–527. [Google Scholar]
  59. Fernandez, A.; Abdeljawad, T.; Baleanu, D. Relations between fractional models with three-parameter Mittag-Leffler kernels. Adv. Differ. Equ. 2020, 2020, 186. [Google Scholar] [CrossRef]
  60. Garrappa, R.; Maione, G. Fractional Prabhakar derivative and applications in anomalous dielectrics: A numerical approach. In Theory and Applications of Non-Integer Order Systems, Lecture Notes in Electrical Engineering 407; Babiarz, A., Czornik, A., Klamka, J., Niezabitowski, M., Eds.; Springer: Cham, Switzerland, 2017; pp. 429–439. [Google Scholar]
  61. Capelas de Oliveira, E.; Mainardi, F.; Vaz, J., Jr. Fractional models of anomalous relaxation based on the Kilbas and Saigo function. Meccanica 2014, 49, 2049–2060. [Google Scholar] [CrossRef]
  62. Skovaranek, T. The Mittag-Leffler fitting of the Phillips curve. Mathematics 2019, 7, 589. [Google Scholar] [CrossRef] [Green Version]
  63. Hristov, J. Derivatives with non-singular kernels: From the Caputo-Fabrizio definition and beyond: Appraising analysis with emphasis on diffusion models. In Frontiers in Fractional Calculus; Bhalekar, S., Ed.; Bentham Science Publishers: Sharjah, United Arab Emirates, 2018; pp. 269–342. [Google Scholar]
  64. Hristov, J. Linear viscoelastic responses and constitutive equations in terms of fractional operators with non-singular kernels: Pragmatic approach, Memory kernel correspondence requirement and analyses. Eur. Phys. J. Plus 2019, 134, 283. [Google Scholar] [CrossRef]
  65. Hristov, J. Response functions in linear viscoelastic constitutive equations and related fractional operators. Math. Model. Nat. Phenom. 2019, 14, 305. [Google Scholar] [CrossRef] [Green Version]
  66. Hristov, J. Transient heat diffusion with a non-singular fading memory: From the Cattaneo constitutive equation with Jeffrey’s kernel to the Caputo-Fabrizio time-fractional derivative. Therm. Sci. 2016, 20, 765–770. [Google Scholar] [CrossRef]
  67. Coleman, B.; Gurtin, M.E. Equipresence and constitutive equations for rigid heat conductors. Z. Angew. Math. Phys. 1967, 18, 188–208. [Google Scholar] [CrossRef]
  68. Gurtin, M.E.; Pipkin, A.C. A general theory of heat conduction with finite wave speeds. Arch. Ration. Mech. Anal. 1968, 31, 113–126. [Google Scholar] [CrossRef]
  69. Miller, R.K. An integrodifferential equation for rigid heat conductors with memory. J. Math. Anal. Appl. 1978, 20, 313–332. [Google Scholar] [CrossRef] [Green Version]
  70. Hilfer, R. Fractional Time evolution, In Applications of Fractional Calculus in Physics; Hilfer, R., Ed.; World Scientific: Singapore, 2000; pp. 89–116. [Google Scholar]
  71. Coleman, B. Thermodynamics of Materials with Memory; Springer: Wien, Austria, 1971. [Google Scholar]
  72. Fabrizio, M.; Morro, A. Thermodynamic restrictions on relaxation functions in linear viscoelasticity. Mech. Res. Commun. 1985, 12, 101–105. [Google Scholar] [CrossRef]
  73. Morro, A. A thermodynamic approach to rate equations in continuum physics. J. Phys. Sci. Appl. 2017, 7, 15–23. [Google Scholar]
  74. Morro, A. Thermodynamic consistency of objective rate equations. Mech. Res. Commun. 2017, 84, 72–76. [Google Scholar] [CrossRef]
  75. Morro, A.; Giorgi, C. Objective rate equations and memory properties in continuum physics. Math. Comp. Sim. 2020, 176, 243–253. [Google Scholar] [CrossRef]
  76. Bouchaud, J.-P.; Georges, A. Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications. Phys. Rep. 1990, 195, 127–293. [Google Scholar] [CrossRef]
  77. Zaslavsky, G.M. Chaos, fractional kinetics, and anomalous transport. Phys. Rep. 2002, 371, 461–580. [Google Scholar] [CrossRef]
  78. Ladman, U.; Montroll, E.W.; Shlesinger, M.F. Random walks and generalized master equations with integral degrees of freedom. Proc. Natl. Acad. Sci. USA 1977, 74, 430–433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Stanislavsky, A.; Weron, A. Transient diffusion with Prabhakar-type memory. J. Chem. Phys. 2018, 149, 044107. [Google Scholar] [CrossRef] [PubMed]
  80. Stanislavsky, A.; Weron, A. Control of the transient subdiffusion exponent at a short and long times. Phys. Rev. Res. 2019, 1, 023006. [Google Scholar] [CrossRef] [Green Version]
  81. Gajda, J.; Beghin, L. Prabhakar Levy processes. Stat. Probab. Lett. 2021, 178, 109162. [Google Scholar] [CrossRef]
  82. Storm, M.L. Heat conduction in simple metals. J. Appl. Phys. 1951, 22, 940–951. [Google Scholar] [CrossRef]
  83. Joseph, D.D.; Preciozi, L. Heat waves. Rev. Mod. Phys. 1989, 61, 41–73. [Google Scholar] [CrossRef]
  84. Tarasov, V.E.; Tarasova, V.V. Logistic equation with continuously distributed lag and applications in economics. Nonlinear Dyn. 2019, 97, 1313–1328. [Google Scholar] [CrossRef]
  85. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers, vol. 1. Background and Theory; Higher Education Press: Beijing, China; Springer: Berlin/Heidelberg, Germany, 2013; pp. 62–63. [Google Scholar]
  86. Hristov, J. Transient heat conduction with non-singular memory: Heat flux equation with a Mittag-Leffler memory naturally leads to ABC derivative. Therm. Sci. 2023; 27, in press. [Google Scholar] [CrossRef]
  87. Rzanitsyn, A.R. Some Questions in Mechanics of Systems Deformed with Time; Gostekhiizadat: Moscow, Russia, 1949. (In Russian) [Google Scholar]
  88. Selivanov, M.F. Effective properties of a linear viscoelastic composites. Int. J. Appl. Mech 2009, 45, 1084–1091. [Google Scholar] [CrossRef]
  89. Shitikova, M.V. Fractional operator viscoelastic models in dynamic problems of mechanics of solids: A review. Mech. Solids 2022, 57, 1–33. [Google Scholar] [CrossRef]
  90. Samko, S.G.; Cardoso, R.P. Integral equations of the first kind of Sonine type. Int. J. Math. Math. Sci. 2003, 57, 3609–3632. [Google Scholar] [CrossRef] [Green Version]
  91. Samko, S.G.; Cardoso, R.P. Sonine integral equations of the first kind. Frac. Calc. Appl. Anal. 2003, 6, 235–258. [Google Scholar] [CrossRef] [Green Version]
  92. Yang, X.-J. General Fractional Derivatives: Theory, Methods and Applications; CRC Press: New York, NY, USA, 2019. [Google Scholar]
  93. Feng, Y.-Y.; Yang, X.-Y.; Liou, J.-G.; Chen, Z.-Q. Rheological analysis of the general fractional-order viscoelastic model involving the Miller-Ross kernel. Acta Mech. 2021, 232, 3141–3148. [Google Scholar] [CrossRef]
  94. Feng, Y.-Y.; Liu, J.-G. Anomalous diffusion equation using a new general fractional derivative within the Miller-Ross kernel. Mod. Phys. Lett. B 2020, 34, 2050289. [Google Scholar] [CrossRef]
  95. Rabotnov, Y.N. Creep Problems in Structural Members; North-Holland: Amsterdam, The Netherlands, 1969. [Google Scholar]
  96. Yang, X.-J.; Abdel-Aty, M.; Cattani, C. A new general fractional-order derivative with Rabotnov fractional exponential kernel applied to the anomalous heat transfer. Therm. Sci. 2019, 23, 1677–1681. [Google Scholar] [CrossRef] [Green Version]
  97. Yang, X.-J.; Ragulskis, M.; Taha, T. A new general fractional-order derivative with Rabotnov fractional exponential kernel. Therm. Sci. 2019, 23, 3711–3718. [Google Scholar] [CrossRef]
  98. Kumar, S.; Ghosh, S.; Samet, B.; Goufo, E.F.D. An analysis of heat equations arises in diffusion process using new Yang-Abdel-Aty-Cattani fractional operator. Math. Model. Appl. Sci. 2020, 43, 6062–6080. [Google Scholar] [CrossRef]
  99. Kumar, S.; Ghosh, S.; Lotayif, M.S.M.; Samet, B. A model describing the velocity of a particle in Brownian motion by Robotnov function based fractional operator. Alex. Eng. J. 2020, 59, 1435–1449. [Google Scholar] [CrossRef]
  100. Sene, N. Fractional diffusion equation with new fractional operator. Alex. Eng. J. 2020, 59, 2921–2926. [Google Scholar] [CrossRef]
  101. Garrappa, R.; Mainardi, F.; Maione, G. Models of dielectric relaxation based on completely monotone functions. Frac. Calc. Appl. Anal. 2016, 19, 1105–1160. [Google Scholar] [CrossRef] [Green Version]
  102. Hristov, J. Integral solutions to transient nonlinear heat (mass) diffusion with a power-law diffusivity: A semi-infinite medium with fixed boundary conditions. Heat Mass Transf. 2016, 52, 635–655. [Google Scholar] [CrossRef]
  103. Dos Santos, A.F. Fractional Prabhakar derivative in diffusion equation with non-static stochastic resetting. Physics 2019, 1, 40–58. [Google Scholar] [CrossRef] [Green Version]
  104. Sandev, T. Generalized Langevin equation and the Prabhakar derivative. Mathematics 2017, 5, 66. [Google Scholar] [CrossRef] [Green Version]
  105. Sandev, T.; Deng, W.; Xu, P. Models for characterizing the transition among anomalous diffusions with different diffusion exponents. J. Phys. A Math. Theor. 2018, 51, 405002. [Google Scholar] [CrossRef] [Green Version]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hristov, J. The Fading Memory Formalism with Mittag-Leffler-Type Kernels as A Generator of Non-Local Operators. Appl. Sci. 2023, 13, 3065. https://doi.org/10.3390/app13053065

AMA Style

Hristov J. The Fading Memory Formalism with Mittag-Leffler-Type Kernels as A Generator of Non-Local Operators. Applied Sciences. 2023; 13(5):3065. https://doi.org/10.3390/app13053065

Chicago/Turabian Style

Hristov, Jordan. 2023. "The Fading Memory Formalism with Mittag-Leffler-Type Kernels as A Generator of Non-Local Operators" Applied Sciences 13, no. 5: 3065. https://doi.org/10.3390/app13053065

APA Style

Hristov, J. (2023). The Fading Memory Formalism with Mittag-Leffler-Type Kernels as A Generator of Non-Local Operators. Applied Sciences, 13(5), 3065. https://doi.org/10.3390/app13053065

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop