Next Article in Journal
Recovering Density and Speed of Sound Coefficients in the 2D Hyperbolic System of Acoustic Equations of the First Order by a Finite Number of Observations
Next Article in Special Issue
Some Properties of the Kilbas-Saigo Function
Previous Article in Journal
Acknowledgment to Reviewers of Mathematics in 2020
Previous Article in Special Issue
Approximation of CDF of Non-Central Chi-Square Distribution by Mean-Value Theorems for Integrals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Some Applications of the Wright Function in Continuum Physics: A Survey

Department of Mathematics and Computer Science, Faculty of Science and Technology, Jan Dlugosz University in Czestochowa, al. Armii Krajowej 13/15, 42-200 Czestochowa, Poland
Mathematics 2021, 9(2), 198; https://doi.org/10.3390/math9020198
Submission received: 23 December 2020 / Revised: 15 January 2021 / Accepted: 17 January 2021 / Published: 19 January 2021
(This article belongs to the Special Issue Special Functions with Applications to Mathematical Physics)

Abstract

:
The Wright function is a generalization of the exponential function and the Bessel functions. Integral relations between the Mittag–Leffler functions and the Wright function are presented. The applications of the Wright function and the Mainardi function to description of diffusion, heat conduction, thermal and diffusive stresses, and nonlocal elasticity in the framework of fractional calculus are discussed.

1. Introduction

The fractional calculus (the theory of integrals and derivatives of non-integer order) has attracted considerable interest of researchers and has many applications in physics, chemistry, rheology, geology, hydrology, medicine, engineering, finance, etc. (see, for example, West–Bologna–Grigolini [1], Magin [2], Povstenko [3], Tarasov [4], Povstenko [5], Uchaikin [6], Atanacković–Pilipović–Stanković–Zorica [7], Herrmann [8], Povstenko [9], Datsko–Gafiychuk–Podlubny [10], West [11], Skiadas [12], Tarasov [13], Kumar–Singh [14], Su [15] and references therein). The Mittag–Leffler functions and the Wright function appear in solutions of various types of equations with fractional operators. The Mittag–Leffler function in one parameter E α z was introduced in [16,17]. The generalized Mittag–Leffler function in two parameters E α , β z was considered in [18,19]. A comprehensive treatment of properties of the Mittag–Leffler functions can be found in Erdélyi–Magnus–Oberhettinger–Tricomi [20], Gorenflo–Mainardi [21], Podlubny [22], Kilbas–Srivastava–Trujillo [23], Gorenflo–Kilbas–Mainardi–Rogosin [24]. Numerical algorithms for calculation of the Mittag–Leffler functions were proposed in [25] and implemented in [26]. The Wright function was presented in [27,28] and later on discussed by Erdélyi–Magnus–Oberhettinger–Tricomi [20], Gorenflo–Mainardi [21], Podlubny [22], Kilbas–Srivastava–Trujillo [23], Gorenflo–Kilbas–Mainardi–Rogosin [24], Luchko [29], among others. Numerical algorithms for calculating the Wright function were suggested in [30].
In 1996, Mainardi [31,32] solved the diffusion-wave equation with the Caputo fractional derivative of the order α
α T t α = a 2 T x 2 , 0 < α 2 ,
on a real line (the Cauchy problem) and a half-line (the signaling problem). The solutions were obtained in terms of the Mainardi function M z ; α 2 [33], where
z = | x | a t α / 2
is the similarity variable, a can be treated as the generalized thermal diffusivity coefficient.
Equation (1) in the limiting case α 0 corresponds to the Helmholtz equation (localized diffusion); the subdiffusion regime is characterized by the values 0 < α < 1 . For 1 < α < 2 , the diffusion-wave Equation (1) interpolates between the diffusion equation ( α = 1 ) and the wave equation ( α = 2 ).
Applications of fractional calculus to viscoelasticity have been studied by many authors. The historical notes and the extensive bibliography on this subject can be found in the book of Mainardi [34]. According to the Scott–Blair stress-strain law, the dependence between the stress σ ( x , t ) and the strain ϵ ( x , t ) can be written as [34,35]
σ ( x , t ) = ρ a ν ϵ ( x , t ) t ν , 0 ν 1 .
The constitutive Equation (3) characterizes a viscoelastic material intermediate between a perfectly elastic solid (the Hooke law for the value ν = 0 ) and a perfectly viscous fluid (the Newton law when ν = 1 ) with the corresponding interpretations of the coefficient a in terms of the elasticity constant or the kinematic viscosity. The relation (3) leads to the evolution Equation (1) with α = 2 ν .
The book [36] presents a picture of the state-of-the-art for solutions of the diffusion-wave equation with one, two, and three space variables in Cartesian, cylindrical, and spherical coordinates under different kinds of boundary conditions.
In the present survey article, we briefly discuss the properties of the Mittag–Leffler functions and Wright function and present the integral relations between the Mittag–Leffler functions and the Wright function. The applications of the Wright function and the Mainardi function to the description of diffusion, heat conduction, thermal and diffusive stresses, and nonlocal elasticity in the framework of fractional calculus are reviewed.

2. Mathematical Preliminaries

2.1. Integrals and Derivatives of Fractional Order

The Riemann–Liouville integral of fractional order α is defined as [21,22,23]:
I α f ( t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 f ( τ ) d τ , α > 0 ,
where Γ ( α ) is the gamma function.
The Riemann–Liouville derivative of fractional order α has the form
D R L α f ( t ) = d n d t n 1 Γ ( n α ) 0 t ( t τ ) n α 1 f ( τ ) d τ , n 1 < α < n ,
whereas the Caputo fractional derivative is written as
D C α f ( t ) d α f ( t ) d t α = 1 Γ ( n α ) 0 t ( t τ ) n α 1 d n f ( τ ) d τ n d τ , n 1 < α < n .
The fractional operators have the following Laplace transform rules:
L I α f ( t ) = 1 s α f * ( s ) ,
L D R L α f ( t ) = s α f * ( s ) k = 0 n 1 D k I n α f ( 0 + ) s n 1 k , n 1 < α < n ,
L d α f d t α = s α f * ( s ) k = 0 n 1 f ( k ) ( 0 + ) s α 1 k , n 1 < α < n .
Here, the asterisk denotes the transform, and s is the Laplace transform variable.

2.2. Mittag–Leffler Functions

The Mittag–Leffler function in one parameter α
E α ( z ) = k = 0 z k Γ ( α k + 1 ) , α > 0 , z C ,
can be considered as the extension of the exponential function e z = E 1 ( z ) , whereas the generalized Mittag–Leffler function in two parameters α and β is defined by the series representation
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , α > 0 , β > 0 , z C .
In the general case, the parameters α and β can be treated as complex numbers with some limitations on their real parts [24], but we restrict ourselves to positive values of α and β .
The following recurrence relations [20,24]
E α , β ( z ) = 1 Γ ( β ) + z E α , α + β ( z ) .
E α , β ( z ) = β E α , β + 1 ( z ) + α z d E α , β + 1 ( z ) d z
are valid for the Mittag–Leffler functions.
For investigation of the convergence of integrals containing the Mittag–Leffler functions, their asymtotic representations for large negative values of argument are useful. For x , we have
E α ( x ) 1 Γ ( 1 α ) x ,
E α , 2 ( x ) 1 Γ ( 2 α ) x ,
E α , α ( x ) 1 Γ ( α ) x 2 ,
E α , β ( x ) 1 Γ ( β α ) x .
The essential role of the Mittag–Leffler functions in fractional calculus is connected with the formula for the inverse Laplace transform (see Gorenflo–Mainardi [21], Podlubny [22], Kilbas–Srivastava–Trujillo [23], Gorenflo–Kilbas–Mainardi–Rogosin [24]):
L 1 s α β s α + b = t β 1 E α , β ( b t α ) .

2.3. Wright Function and Mainardi Function

The Wright function is a generalization of the exponential function and the Bessel functions and is defined as [27,28] (see also refs. [20,21,22,23,24,31,32,37,38,39])
W ( α , β ; z ) = k = 0 z k k ! Γ ( α k + β ) , α > 1 , β C , z C .
The Wright function satisfies the recurrence equations [20]
α z W ( α , α + β ; z ) = W ( α , β 1 ; z ) + ( 1 β ) W ( α , β ; z ) ,
d W ( α , β ; z ) d z = W ( α , α + β ; z ) .
The Mainardi function M ( α ; z ) [22,31,32,33] is a particular case of the Wright function
M ( α ; z ) = W ( α , 1 α ; z ) = k = 0 ( 1 ) k z k k ! Γ [ α k + ( 1 α ) ] , 0 < α < 1 , z C .
The Wright function and the Mainardi function appear in formulae for the inverse Laplace transform (see Mainardi [31,32], Stanković [40], Gajić–Stanković [41]):
L 1 exp λ s α = α λ t α + 1 M α ; λ t α , 0 < α < 1 , λ > 0 ,
L 1 s α 1 exp λ s α = 1 t α M α ; λ t α , 0 < α < 1 , λ > 0 ,
L 1 s β exp λ s α = t β 1 W α , β ; λ t α , 0 < α < 1 , λ > 0 .

2.4. The Integral Transform Relations between the Mittag–Leffler Function and Wright Function

The Laplace transform of the Wright function is expressed in terms of the Mittag–Leffler function [20,22,23]
L W ( α , β ; t ) = 1 s E α , β 1 s , α > 0 , β > 0 ,
and [37]
L W ( α , β ; t ) = E α , β α s , 1 < α < 0 , β > 0 ,
whereas, for the Mainardi function, the corresponding relation takes the form
L M α ; t = E α s , 0 < α < 1 .
The Mittag–Leffler functions and the Wright function are related by the Fourier cosine transform (Povstenko [36,42]):
0 E α ( ξ 2 ) cos ( x ξ ) d ξ = π 2 M α 2 ; x , 0 < α < 2 , x > 0 ,
0 E α , 2 ( ξ 2 ) cos ( x ξ ) d ξ = π 2 W α 2 , 2 α 2 ; x , 0 < α < 2 , x > 0 ,
0 E α , α ( ξ 2 ) cos ( x ξ ) d ξ = π 2 W α 2 , α 2 ; x , 0 < α < 2 , x > 0 ,
0 E α , β ( ξ 2 ) cos ( x ξ ) d ξ = π 2 W α 2 , β α 2 ; x , 0 < α < 2 , β > 0 , x > 0 ,
as well as by the Fourier sine transform
0 ξ E α ( ξ 2 ) sin ( x ξ ) d ξ = π 2 W α 2 , 1 α ; x , 0 < α < 2 , x > 0 ,
0 ξ E α , 2 ( ξ 2 ) sin ( x ξ ) d ξ = π 2 W α 2 , 2 α ; x , 0 < α < 2 , x > 0 ,
0 ξ E α , α ( ξ 2 ) sin ( x ξ ) d ξ = π 2 W α 2 , 0 ; x = α π 4 x M α 2 ; x , 0 < α < 2 , x > 0 ,
0 ξ E α , β ( ξ 2 ) sin ( x ξ ) d ξ = π 2 W α 2 , β α ; x , 0 < α < 2 , β > 0 , x > 0 .
Due to (16), we can also obtain for E α , α ξ 2
0 ξ 2 E α , α ( ξ 2 ) cos ( x ξ ) d ξ = π 2 W α 2 , α 2 ; x , 0 < α < 2 , x > 0 ,
0 ξ 3 E α , α ( ξ 2 ) sin ( x ξ ) d ξ = π 2 W α 2 , α ; x , 0 < α < 2 , x > 0 .
The equations presented above allow us to obtain additional integral relations between the Mittag–Leffler functions and the Wright function, which can be helpful when solving problems in polar or cylindrical coordinates using the Hankel transform of order zero. Taking into account the integral representation of the Bessel function J 0 ( x ) (Watson [43], Abramowitz-Stegun [44])
J 0 ( x ) = 1 π 0 π cos x sin θ d θ ,
J 0 ( x ) = 2 π 0 sin x cosh t d t , x > 0 ,
J 0 ( x ) = 2 π 1 sin x t t 2 1 d t , x > 0 ,
we get
0 E α ξ 2 J 0 r ξ d ξ = 1 2 0 π M α 2 ; r sin θ d θ , 0 < α 2 , r > 0 ,
0 E α , 2 ξ 2 J 0 r ξ d ξ = 1 2 0 π W α 2 , 2 α 2 ; r sin θ d θ , 0 < α 2 , r > 0 ,
0 E α , α ξ 2 J 0 r ξ d ξ = 1 2 0 π W α 2 , α 2 ; r sin θ d θ , 0 < α 2 , r > 0 ,
0 E α , β ξ 2 J 0 r ξ d ξ = 1 2 0 π W α 2 , β α 2 ; r sin θ d θ , 0 < α 2 , r > 0 .
Similarly,
0 E α ξ 2 J 0 r ξ ξ d ξ = 0 W α 2 , 1 α ; r cosh t d t , 0 < α 2 , r > 0 ,
0 E α , 2 ξ 2 J 0 r ξ ξ d ξ = 0 W α 2 , 2 α ; r cosh t d t , 0 < α 2 , r > 0 ,
0 E α , α ξ 2 J 0 r ξ ξ d ξ = 0 W α 2 , 0 ; r cosh t d t , 0 < α 2 , r > 0 ,
0 E α , β ξ 2 J 0 r ξ ξ d ξ = 0 W α 2 , β α ; r cosh t d t , 0 < α 2 , r > 0 ,
and
0 E α ξ 2 J 0 r ξ ξ d ξ = 1 W α 2 , 1 α ; r t 1 t 2 1 d t , 0 < α 2 , r > 0 ,
0 E α , 2 ξ 2 J 0 r ξ ξ d ξ = 1 W α 2 , 2 α ; r t 1 t 2 1 d t , 0 < α 2 , r > 0 ,
0 E α , α ξ 2 J 0 r ξ ξ d ξ = 1 W α 2 , 0 ; r t 1 t 2 1 d t = α r 2 0 M α 2 ; r 1 + u 2 d u , 0 < α 2 , r > 0 ,
0 E α , β ξ 2 J 0 r ξ ξ d ξ = 1 W α 2 , β α ; r t 1 t 2 1 d t , 0 < α 2 , r > 0 .
In addition,
0 E α , α ξ 2 J 0 r ξ ξ 2 d ξ = 1 2 0 π W α 2 , α 2 ; r sin θ d θ , 0 < α 2 , r > 0 ,
0 E α , α ξ 2 J 0 r ξ ξ 3 d ξ = 0 W α 2 , α ; r cosh t d t = 1 W α 2 , α ; r t 1 t 2 1 d t , 0 < α 2 , r > 0 .

3. Applications of the Wright Function

3.1. Fractional Heat Conduction in Nonhomogeneous Media under Perfect Thermal Contact

Time-fractional heat conduction in two joint half-lines was considered by Povstenko [36,45,46]. In the general case, the heat conduction equation with the Caputo derivative of the order 0 < α 2 in one half-line
α T 1 t α = a 1 2 T 1 x 2 , x > 0 ,
and the corresponding equation with the Caputo derivative of the order 0 < β 2 in another half-line
β T 2 t β = a 2 2 T 2 x 2 , x < 0 ,
were treated under the boundary conditions of perfect thermal contact which state that two bodies must have the same temperature at the contact point and the heat fluxes through the contact point must be the same:
T 1 ( x , t ) | x = 0 + = T 2 ( x , t ) | x = 0 ,
k 1 D R L 1 α T 1 ( x , t ) x | x = 0 + = k 2 D R L 1 β T 2 ( x , t ) x | x = 0 , 0 < α 2 , 0 < β 2 .
In the condition (59), k 1 and k 2 are the generalized thermal conductivities of two bodies; the Riemann–Liouville fractional derivative of the negative order D R L α ( f ( t ) ) is understood as the Riemann–Liouville fractional integral I α ( f ( t ) .
Here, we present the fundamental solution to the first Cauchy problem with the initial condition
t = 0 : T 1 = p 0 δ ( x ϱ ) , x > 0 , ϱ > 0 ,
for the case α = β (for details see Povstenko [46]):
T 1 ( x , t ) = p 0 2 a 1 t α / 2 M α 2 ; | x ϱ | a 1 t α / 2 + ε 1 ε + 1 M α 2 ; x + ϱ a 1 t α / 2 , x 0 ,
T 2 ( x , t ) = ε p 0 ( ε + 1 ) a 1 t α / 2 M α 2 ; | x | a 2 t α / 2 + ρ a 1 t α / 2 , x 0 ,
where
ε = k 1 a 2 k 2 a 1 .
For the corresponding problem with uniform initial temperature T 0 in one of half-lines [45], in the particular case α = β , we have:
T 1 = T 0 T 0 ( 1 + ε ) W α 2 , 1 ; x a 1 t α / 2 , x > 0 ,
T 2 = ε T 0 ( 1 + ε ) W α 2 , 1 ; | x | a 2 t α / 2 , x < 0 .
The time-fractional heat conduction equations with the Caputo derivatives in a semi-infinite medium composed of a region 0 < x < L and a region L < x < under the boundary conditions of perfect thermal contact at x = L and the insulated boundary condition at x = 0 with uniform initial temperature in a layer were investigated in [47]. The approximate solution of the considered problem for small values of time is obtained based on Tauberian theorems for the Laplace transform. For α = β , this solution reads
T 1 T 0 T 0 1 + ε W α 2 , 1 ; L x a 1 t α / 2 , 0 x L ,
T 2 ε T 0 1 + ε W α 2 , 1 ; x L a 2 t α / 2 , L x < .
Fractional heat conduction in an infinite medium with a spherical inclusion when a sphere 0 r < R is at the initial uniform temperature T 0 and a matrix R < r < is at a zero initial temperature was considered by Povstenko [36,48]. In the case of perfect thermal contact at the boundary r = R ,
r = R : T 1 ( r , t ) = T 2 ( r , t ) ,
k 1 D R L 1 α T 1 ( r , t ) r = k 2 D R L 1 β T 2 ( r , t ) r , 0 < α 2 , 0 < β 2 ,
the approximate solution for small values of time has the following form (we present only the solution for α = β ):
T 1 ( r , t ) T 0 R T 0 k 2 ( k 2 k 1 ) r W α 2 , 1 ; R r a 1 t α / 2 W α 2 , 1 ; R + r a 1 t α / 2 + c R T 0 r 0 t ( t τ ) α / 2 1 τ α / 2 [ M α 2 ; R r a 1 τ α / 2 M α 2 ; R + r a 1 τ α / 2 ] E α / 2 , α / 2 b ( t τ ) α / 2 d τ ,
T 2 ( r , t ) R T 0 k 1 ( k 2 k 1 ) r W α 2 , 1 ; r R a 2 t α / 2 + c R T 0 r 0 t ( t τ ) α / 2 1 τ α / 2 ×   M α 2 ; r R a 2 τ α / 2 E α / 2 , α / 2 b ( t τ ) α / 2 d τ ,
where
b = ( k 2 k 1 ) a 1 a 2 R k 1 a 1 + k 2 a 2 , c = k 1 k 2 a 1 + a 2 ( k 2 k 1 ) k 1 a 1 + k 2 a 2 .
It should be mentioned that, for the classical heat conduction, the method of analysis of the solution for small values of time was described by Luikov [49] and Özişik [50]. In the case of fractional diffusion equation, the decay rate at large values of time was analyzed by Sakamoto–Yamamoto [51].

3.2. Fractional Heat Conduction in Nonhomogeneous Media under Nonperfect Thermal Contact

Near the interface between two solids, a transition region arises whose state differs from the state of contacting media owing to different conditions of material–particle interaction. The transition region has its own physical, mechanical, and chemical properties, and processes occurring in it differ from those in the bulk. Small thickness of the intermediate region between two solids allows us to reduce a three-dimensional problem to a two-dimensional one for median surface endowed with equivalent physical properties. There are several approaches to reducing three-dimensional equations to the corresponding two-dimensional equations for the median surface. For example, introducing the mixed coordinate system ( ξ , η , z ) , where ξ and η are the curvilinear coordinates in the median surface and z is the normal coordinate, the linear or polynomial dependence of the considered functions on the normal coordinate can be assumed. This assumption is often used in the theory of elastic shells.
For the classical heat conduction equation, which is based on the conventional Fourier law, the reduction of the three-dimensional problem to the simplified two-dimensional one was pioneered by Marguerre [52,53] and later on developed by many authors. In this case, the assumption of linear or polynomial dependence of temperature on the normal coordinate or more general operator method were used. An extensive literature on this subject can be found, for example, in [9]. For time-fractional heat conduction, the reduction of the three-dimensional equation to the two-dimensional one was carried out by Povstenko [9,54,55].
A solution to the problem (56), (57) with uniform initial temperature in one of half-lines under conditions of nonperfect thermal contact was obtained in [56]. In the particular case α = β , the solution reads
T 1 = T 0 T 0 ( 1 + ε ) W α 2 , 1 ; x a 1 t α / 2 + T 0 ( 1 ε ) 2 ( 1 + ε ) 0 t ( t τ ) α / 2 1 τ α / 2 × M α 2 ; x a 1 τ α / 2 E α / 2 , α / 2 b Σ t τ α / 2 d τ , x > 0 ,
T 2 = ε T 0 ( 1 + ε ) W α 2 , 1 ; | x | a 2 t α / 2 + T 0 ( 1 ε ) 2 ( 1 + ε ) 0 t ( t τ ) α / 2 1 τ α / 2 × M α 2 ; | x | a 2 τ α / 2 E α / 2 , α / 2 b Σ t τ α / 2 d τ , x < 0 ,
where ε is defined by (63),
b Σ = k 1 a 2 + k 2 a 1 C Σ a 1 a 2 ,
C Σ is the reduced heat capacity of the median surface of the transition region. When C Σ 0 , the solutions (73), (74) coincide with the solutions (64), (65).

3.3. Fractional Heat Conduction under Time-Harmonic Impact

Ångström [57] was the first to investigate the standard parabolic heat conduction equation under time-harmonic impact. An extensive review of literature in this field in the case of classical diffusion equation can be found in the book by Mandelis [58].
Fractional heat conduction with a source varying harmonically in time was studied by Povstenko [59]. Equation (1) with a source term
α T t α = a 2 T x 2 + Q 0 δ x e i ω t , 0 < α 2 ,
was solved in the domain < x < under zero initial conditions. Temperature is expressed as
T ( x , t ) = Q 0 2 a 0 t τ α / 2 1 W α 2 , α 2 ; | x | a τ α / 2 e i ω ( t τ d τ .
The corresponding problem in the central symmetric case
α T t α = a 2 T r 2 + 2 r T r + Q 0 δ r 4 π r 2 e i ω t , 0 < r < , 0 < α 2 ,
has the solution
T ( x , t ) = α Q 0 8 π a 3 / 2 0 t 1 τ 1 + α / 2 M α 2 ; r a τ α / 2 e i ω ( t τ d τ .

3.4. Fractional Nonlocal Elasticity

Nonlocal continuum physics assumes integral constitutive equations. In the nonlocal theory of the continuum mechanics, stresses at the reference point x of an elastic solid at time t depend not only on the strains at this point at this time, but also on strains at all the points x of a body and all the times prior to and at time t:
t x , t , ϵ L , ϵ T = 0 t V γ | x x | , t t , ϵ L , ϵ T σ x , t d v ( x ) d t ,
σ x , t = 2 μ e x , t + λ tr e x , t I ,
where t and σ are the nonlocal and classical stress tensors, x and x are the reference and running points, e the linear strain tensor, λ and μ are Lamé constants, I stands for the unit tensor. The volume integral in (80) is over the region occupied by the solid. The time-non-locality describes memory effects, distributed lag (distributed time delay), and frequency dispersion; the space-non-locality deals with the long-range interaction. The weight function (the non-locality kernel) γ | x x | , t t , ϵ L , ϵ T depends on two basic non-locality parameters (see Eringen [60]): the characteristic length ratio
ϵ L = Internal characteristic length External characteristic length
and the characteristic time ratio
ϵ T = Internal characteristic time External characteristic time .
When ϵ T 0 , the memory effects are eliminated; for ϵ L 0 the space-non-locality disappears.
In the pioneering works by Podstrigach [61,62], a new nontraditional thermodynamic pair (the chemical potential tensor φ and the concentration tensor c ) was introduced (see also [63,64]). The tensor character of the chemical potential means that, for solids, the work of bringing the substance into a point in a body depends on the direction. In this case, the diffusion equation, split into the mean and deviatoric parts, has the form
ρ ( tr c ) t = 3 a Δ ( tr φ ) ,
ρ ( dev c ) t = 2 a 1 Δ ( dev φ ) ,
where ρ is the mass density, and a and a 1 are the corresponding diffusion coefficients.
Starting from interrelated equations describing elasticity and diffusion, Podstrigach [65] eliminated the chemical potential tensor from the constitutive equation for the stress tensor and obtained the stress–strain relation containing spatial and time derivatives. In the infinite medium, this relation can be integrated using the Fourier and Laplace integral transforms, and the final result, written for the mean and deviatoric parts, has the nonlocal integral form:
tr σ = 3 K c tr e + 3 K φ K c p 0 t γ ( p ) x x , y y , z z , t t × tr e x , y , z , t d x d y d z d t ,
dev σ = 2 μ c dev e + 2 μ φ μ c q 0 t γ ( q ) x x , y y , z z , t t × dev e x , y , z , t d x d y d z d t .
Here, K c , K φ , μ c , μ φ , p, and q are material constants (for details, see [42,65]). The kernel γ ( p ) ( x , y , z , t ) has the following form:
γ ( p ) ( x , y , z , t ) = p 2 t 5 / 2 3 p r 2 2 t exp p r 2 4 t ,
where r = x 2 + y 2 + z 2 ; the kernel γ ( q ) x , y , z , t is obtained from the kernel γ ( p ) x , y , z , t substituting p by q.
The results of Podstrigach [65] were generalized by Povstenko [42] for the case of fractional diffusion equations
ρ α ( tr c ) t α = 3 a Δ ( tr φ ) ,
ρ α ( dev c ) t α = 2 a 1 Δ ( dev φ ) .
The kernel γ ( p ) ( x , y , z , t ) in the fractional generalization of the constitutive Equation (84) for the mean part of the stress tensor is expressed in terms of the Wright function:
γ ( p ) ( x , y , z , t ) = π p 2 2 t α + 1 r W α 2 , α ; p r t α / 2 .
The kernel γ ( q ) ( x , y , z , t ) in the fractional generalization of the constitutive Equation (85) for the deviatoric part of the stress tensor is obtained by substituting p with q.
In the case of only space-non-locality, the constitutive equation for the stress tensor reads
t x , ϵ L = V γ | x x | , ϵ L σ x d v ( x ) .
The space-nonlocal elasticity reduces to the classical theory of elasticity in the long wave-length limit and to the atomic lattice theory in the short wave-length limit. Several versions of nonlocal elasticity based on various assumptions were proposed by different authors (see, for example, Podstrigach [65], Eringen [66,67], Kunin [68,69] and references therein).
In the case of space-nonlocal constitutive Equation (90), the nonlocal kernel γ | x x | , ϵ L is a delta sequence and in the classical elasticity limit ϵ L 0 becomes the Dirac delta function. For example, slightly changing the notation, the nonlocal kernel γ | x x | , τ can be considered as the Green function of the Cauchy problem for the diffusion operator (see Eringen [67,70]):
γ x , τ τ a Δ γ x , τ = 0 ,
τ = 0 : γ x , τ = δ ( x ) ,
which results in the kernel
γ x , τ = 1 ( 2 π a τ ) n exp | x | 2 4 a τ
for n = 1 , 2 , 3 space variables. In this case, the nonlocal stress tensor is a solution of the corresponding Cauchy problem:
t x , τ τ a Δ t x , τ = 0 ,
τ = 0 : t x , τ = σ x .
It should be emphasized that in, the formal sense, τ in the initial-value problems (91), (92) and (94), (95) looks like time, but in fact τ is a non-locality parameter related to the space-non-locality characteristic ratio ϵ L .
In the paper [71], the nonlocal kernel γ | x x | , τ was considered as the Green function of the Cauchy problem for the fractional diffusion operator
α γ x , τ τ α a Δ γ x , τ = 0 , 0 < α 1 ,
τ = 0 : γ x , τ = δ ( x ) .
In the framework of this approach, instead of the Cauchy problem (94)–(95), we obtain
α t x , τ τ α a Δ t x , τ = 0 , 0 < α 1 ,
τ = 0 : t x , τ = σ x .
In the case of one spatial coordinate, the nonlocal kernel takes the form
γ x , τ = 1 2 a τ α / 2 M α 2 ; | x | a τ α / 2 , 0 α 1 ,
and in the central symmetric case
γ r , τ = 1 4 π a τ α r W α 2 , 1 α ; r a τ α / 2 , 0 α 1 .

4. Conclusions

In this survey, we have reviewed the main applications of the Wright function and the Mainardi function in continuum physics based essentially on the author’s works. We have presented the integral relations between the Mittag–Leffler functions and the Wright function, which can be useful when solving fractional differential equations. We have restricted ourselves to the standard Mittag–Leffler functions and Wright function. The interested reader is referred to publications on further generalizations of the Mittag–Leffler functions [24,72,73,74,75] and of the Wright function [24,29,76,77,78].

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The author would like to thank the anonymous reviewers for their helpful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. West, B.J.; Bologna, M.; Grigolini, P. Physics of Fractal Operators; Springer: New York, NY, USA, 2003. [Google Scholar]
  2. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Redding, CT, USA, 2006. [Google Scholar]
  3. Povstenko, Y. Fractional heat conduction equation and associated thermal stresses. J. Therm. Stress. 2005, 28, 83–102. [Google Scholar] [CrossRef]
  4. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin, Germany, 2010. [Google Scholar]
  5. Povstenko, Y. Fractional Cattaneo-type equations and generalized thermoelasticity. J. Therm. Stress. 2011, 34, 97–114. [Google Scholar] [CrossRef]
  6. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin, Germany, 2013. [Google Scholar]
  7. Atanacković, T.M.; Pilipović, S.; Stanković, B.; Zorica, D. Fractional Calculus with Applications in Mechanics: Vibrations and Diffusion Processes; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  8. Herrmann, R. Fractional Calculus: An Introduction for Physicists, 2nd ed.; World Scientific: Singapore, 2014. [Google Scholar]
  9. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  10. Datsko, B.; Gafiychuk, V.; Podlubny, I. Solitary travelling auto-waves in fractional reaction–diffusion systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 23, 378–387. [Google Scholar] [CrossRef] [Green Version]
  11. West, B.J. Fractional Calculus View of Complexity: Tomorrow’s Science; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  12. Skiadas, C.H. (Ed.) Fractional Dynamics, Anomalous Transport and Plasma Science; Springer: Cham, Switzerland, 2018. [Google Scholar]
  13. Tarasov, V.E. (Ed.) Handbook of Fractional Calculus with Applications. Volume 4: Application in Physics. Part A; Walter de Gruyter: Berlin, Germany, 2019. [Google Scholar]
  14. Kumar, D.; Singh, J. (Eds.) Fractional Calculus in Medical and Health Science; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  15. Su, N. Fractional Calculus for Hydrology, Soil Science and Geomechanics: An Introduction to Applications; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  16. Mittag–Leffler, G.M. Sur la nouvelle fonction Eα(x). C. R. Acad. Sci. Paris Ser. II 1903, 137, 554–558. [Google Scholar]
  17. Mittag–Leffler, G.M. Sopra la funzione Eα(x). Rend. Accad. Lincei Ser. V 1904, 13, 3–5. [Google Scholar]
  18. Humbert, P. Quelques résultats relatifs a‘ la fonction de Mittag–Leffler. C. R. Acad. Sci. Paris 1953, 236, 1467–1468. [Google Scholar]
  19. Humbert, P.; Agarwal, R.P. Sur la fonction de Mittag–Leffler et quelques-unes de ses généralisations. Bull. Sci. Math. 1953, 77, 180–185. [Google Scholar]
  20. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1955; Volume 3. [Google Scholar]
  21. Gorenflo, R.; Mainardi, F. Fractional calculus: Integral and differential equations of fractional order. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: Wien, Austria, 1997; pp. 223–276. [Google Scholar]
  22. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  23. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  24. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag–Leffler Functions, Related Topics and Applications, 2nd ed.; Springer: New York, NY, USA, 2020. [Google Scholar]
  25. Gorenflo, R.; Loutchko, J.; Luchko, Y. Computation of the Mittag–Leffler function and its derivatives. Fract. Calc. Appl. Anal. 2002, 5, 491–518. [Google Scholar]
  26. Podlubny, I. Mittag–Leffler Function; Calculates the Mittag–Leffler Function with Desired Accuracy, MATLAB Central File Exchange, File ID 8738. Available online: www.mathworks.com/matlabcentral/fileexchange/8738 (accessed on 16 November 2020).
  27. Wright, E.M. On the coefficients of power series having exponential singularities. J. Lond. Math. Soc. 1933, 8, 71–79. [Google Scholar] [CrossRef]
  28. Wright, E.M. The asymptotic expansion of the generalized Bessel function. Proc. Lond. Math. Soc. Ser. II 1935, 38, 257–270. [Google Scholar] [CrossRef]
  29. Luchko, Y. The Wright function and its applications. In Handbook of Fractional Calculus with Applications. Volume 1. Basic Theory; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany, 2019; pp. 241–268. [Google Scholar]
  30. Luchko, Y. Algorithms for evaluation of the Wright function for the real arguments’ values. Fract. Calc. Appl. Anal. 2008, 11, 57–75. [Google Scholar]
  31. Mainardi, F. The fundamental solutions for the fractional diffusion-wave equation. Appl. Math. Lett. 1996, 9, 23–28. [Google Scholar] [CrossRef] [Green Version]
  32. Mainardi, F. Fractional relaxation-oscillation and fractional diffusion-wave phenomena. Chaos Solitons Fractals 1996, 7, 1461–1477. [Google Scholar] [CrossRef]
  33. Mainardi, F.; Tomirotti, M. On a special function arising in the time fractional diffusion-wave equation. In Transform Methods & Special Functions, Sofia’ 94; Rusev, P., Dimovski, I., Kiryakova, V., Eds.; Science Culture Technology Publishing: Singapore, 1995; pp. 171–183. [Google Scholar]
  34. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  35. Mainardi, F.; Tomirotti, M. Seismic pulse propagation with constant Q and stable probability distributions. Ann. Geofis. 1997, 40, 1311–1328. [Google Scholar]
  36. Povstenko, Y. Linear Fractional Diffusion-Wave Equation for Scientists and Engineers; Birkhäuser: New York, NY, USA, 2015. [Google Scholar]
  37. Gorenflo, R.; Luchko, Y.; Mainardi, F. Analytical properties and applications of the Wright function. Fract. Calc. Appl. Anal. 1999, 2, 383–414. [Google Scholar]
  38. Gorenflo, R.; Luchko, Y.; Mainardi, F. Wright functions as scale-invariant solutions of the diffusion-wave equation. J. Comput. Appl. Math. 2000, 118, 175–191. [Google Scholar] [CrossRef] [Green Version]
  39. Mainardi, F.; Pagnini, G. The Wright functions as solutions of the time-fractional diffusion equation. Appl. Math. Comput. 2003, 141, 51–62. [Google Scholar] [CrossRef] [Green Version]
  40. Stanković, B. On the function of E. M. Wright. Publ. Inst. Math. 1970, 10, 113–124. [Google Scholar]
  41. Gajić, L.; Stanković, B. Some properties of Wright’s function. Publ. Inst. Math. 1976, 20, 91–98. [Google Scholar]
  42. Povstenko, Y. Generalized theory of diffusive stresses associated with the time-fractional diffusion equation and nonlocal constitutive equations for the stress tensor. Comput. Math. Appl. 2019, 78, 1819–1825. [Google Scholar] [CrossRef]
  43. Watson, G.N. A Treatise on the Theory of Bessel Functions, 2nd ed.; Cambridge University Press: Cambridge, UK, 1944. [Google Scholar]
  44. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphics and Mathematical Tables; Dover: New York, NY, USA, 1972. [Google Scholar]
  45. Povstenko, Y. Fractional heat conduction in infinite one-dimensional composite medium. J. Therm. Stress. 2013, 36, 351–363. [Google Scholar] [CrossRef]
  46. Povstenko, Y. Fundamental solutions to time-fractional heat conduction equations in two joint half-lines. Cent. Eur. J. Phys. 2013, 11, 1284–1294. [Google Scholar] [CrossRef] [Green Version]
  47. Povstenko, Y. Fractional heat conduction in a semi-infinite composite body. Comm. Appl. Industr. Math. 2014, 6, 1–13. [Google Scholar] [CrossRef]
  48. Povstenko, Y. Fractional heat conduction in an infinite medium with a spherical inclusion. Entropy 2013, 15, 4122–4133. [Google Scholar] [CrossRef] [Green Version]
  49. Luikov, A.V. Analytical Heat Diffusion Theory; Academic Press: New York, NY, USA, 1968. [Google Scholar]
  50. Özişik, M.N. Heat Conduction; John Wiley: New York, NY, USA, 1980. [Google Scholar]
  51. Sakamoto, K.; Yamamoto, M. Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems. J. Math. Anal. Appl. 2011, 382, 426–447. [Google Scholar] [CrossRef] [Green Version]
  52. Marguerre, K. Thermo-elastische Platten–Gleichungen. Z. Angew. Math. Mech. 1935, 15, 369–372. [Google Scholar] [CrossRef]
  53. Marguerre, K. Temperaturverlauf und Temperaturspannumgen in platten- und schalenformigen Körpern. Ing. Arch. 1937, 8, 216–228. [Google Scholar] [CrossRef]
  54. Povstenko, Y. Fractional thermoelasticity of thin shells. In Shell Structures; Pietraszkiewicz, W., Górski, J., Eds.; CRC Press: Boca Raton, FL, USA, 2014; Volume 3, pp. 141–144. [Google Scholar]
  55. Povstenko, Y. Generalized boundary conditions for the time-fractional advection diffusion equation. Entropy 2015, 17, 4028–4039. [Google Scholar] [CrossRef] [Green Version]
  56. Povstenko, Y.; Kyrylych, T. Fractional heat conduction in solids connected by thin intermediate layer: Nonperfect thermal contact. Contin. Mech. Thermodyn. 2019, 31, 1719–1731. [Google Scholar] [CrossRef] [Green Version]
  57. Ångström, A.J. Neue Methode, das Wärmeleitungsvermögen der Körper zu bestimmen. Ann. Phys. Chem. 1861, 114, 513–530. [Google Scholar] [CrossRef] [Green Version]
  58. Mandelis, A. Diffusion-Wave Fields: Mathematical Methods and Green Functions; Springer: New York, NY, USA, 2001. [Google Scholar]
  59. Povstenko, Y. Fractional heat conduction in a space with a source varying harmonically in time and associated thermal stresses. J. Therm. Stress. 2016, 39, 1442–1450. [Google Scholar] [CrossRef]
  60. Eringen, A.C. Vistas of nonlocal continuum physics. Int. J. Engng. Sci. 1992, 30, 1551–1565. [Google Scholar] [CrossRef]
  61. Pidstryhach, Y.S. Differential equations of the diffusion theory of deformation of a solid. Dopovidi Ukr. Acad. Sci. 1963, 3, 336–340. (In Ukrainian) [Google Scholar]
  62. Podstrigach, Y.S. Diffusion theory of the anelasticity of metals. J. Appl. Mech. Tech. Phys. 1965, 6, 56–60. [Google Scholar] [CrossRef]
  63. Podstrigach, Y.S.; Povstenko, Y. Introduction to Mechanics of Surface Phenomena in Deformable Solids; Naukova Dumka: Kiev, Ukraine, 1985. (In Russian) [Google Scholar]
  64. Povstenko, Y. From the chemical potential tensor and concentration tensor to nonlocal continuum theories. J. Math. Sci. 2020, 249, 389–403. [Google Scholar] [CrossRef]
  65. Podstrigach, Y.S. On a nonlocal theory of solid body deformation. Internat. Appl. Mech. 1967, 3, 44–46. [Google Scholar] [CrossRef]
  66. Eringen, A.C. Linear theory of nonlocal elasticity and dispersion of plane waves. Int. J. Engng. Sci. 1972, 10, 425–435. [Google Scholar] [CrossRef]
  67. Eringen, A.C. Nonlocal Continuum Field Theories; Springer: New York, NY, USA, 2002. [Google Scholar]
  68. Kunin, I.A. Elastic Media with Microstructure I: One-Dimensional Models; Springer: Berlin, Germany, 1982. [Google Scholar]
  69. Kunin, I.A. Elastic Media with Microstructure II: Three-Dimensional Models; Springer: Berlin, Germany, 1983. [Google Scholar]
  70. Eringen, A.C. On differential eqations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  71. Povstenko, Y. Fractional nonlocal elasticity and solutions for straight screw and edge dislocations. Fiz. Mesomekhanika 2020, 23, 35–44. [Google Scholar]
  72. Al-Bassam, M.A.; Luchko, Y.F. On generalized fractional calculus and its application to the solution of integro-diferential equations. J. Fract. Calc. 1995, 7, 69–88. [Google Scholar]
  73. Kiryakova, V. The multi-index Mittag–Leffler functions as an important class of special functions of fractional calculus. Comp. Math. Appl. 2010, 59, 1885–1895. [Google Scholar] [CrossRef] [Green Version]
  74. Luchko, Y. Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation. J. Math. Anal. Appl. 2011, 374, 538–548. [Google Scholar] [CrossRef] [Green Version]
  75. Li, Z.; Liu, Y.; Yamamoto, M. Initial-boundary value problems for multi-term time-fractional diffusion equations with positive constant coefficients. Appl. Math. Comput. 2015, 257, 381–397. [Google Scholar] [CrossRef] [Green Version]
  76. Kilbas, A.A.; Saigo, M.; Trujillo, J.J. On the generalized Wright function. Fract. Calc. Appl. Anal. 2002, 5, 437–460. [Google Scholar]
  77. Kilbas, A.A. Fractional calculus of the generalized Wright function. Fract. Calc. Appl. Anal. 2005, 8, 113–126. [Google Scholar]
  78. Khan, N.U.; Usman, T.; Aman, M. Some properties concerning the analysis of generalized Wright function. J. Comput. Appl. Math. 2020, 376, 112840. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Povstenko, Y. Some Applications of the Wright Function in Continuum Physics: A Survey. Mathematics 2021, 9, 198. https://doi.org/10.3390/math9020198

AMA Style

Povstenko Y. Some Applications of the Wright Function in Continuum Physics: A Survey. Mathematics. 2021; 9(2):198. https://doi.org/10.3390/math9020198

Chicago/Turabian Style

Povstenko, Yuriy. 2021. "Some Applications of the Wright Function in Continuum Physics: A Survey" Mathematics 9, no. 2: 198. https://doi.org/10.3390/math9020198

APA Style

Povstenko, Y. (2021). Some Applications of the Wright Function in Continuum Physics: A Survey. Mathematics, 9(2), 198. https://doi.org/10.3390/math9020198

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