Next Article in Journal
Implementation of Analytical Techniques for the Solution of Nonlinear Fractional Order Sawada–Kotera–Ito Equation
Next Article in Special Issue
Fractional Maclaurin-Type Inequalities for Multiplicatively Convex Functions
Previous Article in Journal
Some Estimates of k-Fractional Integrals for Various Kinds of Exponentially Convex Functions
Previous Article in Special Issue
Fractional Gradient Methods via ψ-Hilfer Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Subordination Principle for Generalized Fractional Zener Models

Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. G. Bonchev, Block 8, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(4), 298; https://doi.org/10.3390/fractalfract7040298
Submission received: 28 February 2023 / Revised: 22 March 2023 / Accepted: 27 March 2023 / Published: 29 March 2023
(This article belongs to the Special Issue New Trends on Generalized Fractional Calculus)

Abstract

:
The fractional Zener constitutive law is frequently used as a model of solid-like viscoelastic behavior. In this work, a class of linear viscoelastic models of Zener type, which generalize the fractional Zener model, is studied by the use of Bernstein functions technique. We prove that the corresponding relaxation moduli are completely monotone functions under appropriate thermodynamic restrictions on the parameters. Based on this property, we study the propagation function and establish the subordination principle for the corresponding Zener-type wave equation, which provides an integral representation of the solution in terms of the propagation function and the solution of a related classical wave equation. The analytical findings are supported by numerical examples.

1. Introduction

The application of fractional derivatives in time in rheological constitutive equations has a long history [1]. One of the advantages of fractional time-derivatives in linear viscoelasticity is that with their help the damping behavior of viscoelastic media is described by using less parameters, compared with integer-order models. For a review of the main aspects of wave propagation in linear homogeneous viscoelastic media and the simplest and most used fractional constitutive models, we refer to [2].
In linear viscoelasticity, the rheological properties of a viscoelastic medium are described through a linear constitutive relation between stress σ and strain ε . Following [2], we restrict our considerations to the uniaxial case, in which σ = σ ( x , t ) and ε = ε ( x , t ) , and consider systems quiescent for all times prior to some starting time t = 0 .
The fractional Zener constitutive law is extensively used as a model of solid-like viscoelastic behavior [2,3]. The fractional Zener stress–strain relation reads
( 1 + a D t α ) σ ( x , t ) = ( 1 + b D t α ) ε ( x , t ) ,
where D t α denotes the Riemann–Liouville fractional derivative [4] of order α ( 0 , 1 ) and a > 0 , b > 0 . Constitutive Equation (1) is introduced for the first time in 1971 by Caputo and Mainardi [5,6]. Another early extensive study of model (1) can be found in [7]. A necessary and sufficient condition for physical acceptability of model (1) is the restriction a b . Under this condition, the corresponding relaxation modulus is a completely monotone function, see, e.g., [2].
Recall that a functions f ( t ) is said to be completely monotone if it is real-valued and infinitely differentiable on R + and
( 1 ) n f ( n ) ( t ) 0 , t > 0 , n = 0 , 1 , 2 ,
Completely monotonic relaxation moduli are commonly encountered in polymer physics and other branches of rheology and mechanics [2,8,9]. For a discussion on the importance of the complete monotonicity property of the relaxation moduli in linear viscoelasticity, we refer to [9].
Wave equations in viscoelastic media described by the fractional Zener model are studied in many works, see, e.g., [3,10] and the references cited there. Various types of distributed-order generalizations of model (1) are proposed and studied along with the corresponding wave equations: see, e.g., [11,12,13,14,15,16]. The complete monotonicity of the relaxation modulus for Zener-type constitutive equations with multiple time-derivatives and distributed-order time-derivatives with power-type weight function is established in [17] under appropriate constraints on the parameters. The proofs are based on a representation of the relaxation modulus as a Laplace transform of a non-negative spectral function. For other applications of this method, we refer to [2,18].
Recently, equations with general fractional differential operators have gained increasing interest, see the pioneering work [19], the review [20], as well as [21,22,23,24,25], to mention only few of many publications from the last years.
In the article [26], published very recently, a generalization of model (1) in the context of general fractional calculus is proposed and studied. In this work, the fractional derivative D t α in the constitutive Equation (1) is replaced by a general convolutional derivative. For two specific cases of the kernel of this derivative, it is proven that the condition a < b is sufficient for the dissipation inequality to hold and thus, for the thermodynamic compatibility of the corresponding viscoelastic models.
Motivated by the work [26], in the present paper, we study the following generalized fractional Zener model
1 + a D t ( k ) σ ( x , t ) = 1 + b D t ( k ) ε ( x , t ) .
Here, a > 0 and b > 0 are given parameters, and D t ( k ) is a generalized convolutional derivative in the Riemann–Liouville sense, defined by
( D t ( k ) f ) ( t ) = d d t 0 t k ( t τ ) f ( τ ) d τ , t > 0 ,
where k is a locally integrable memory kernel, k L l o c 1 ( R + ) . For the kernel k ( t ) , we assume in addition that its Laplace transform k ^ ( s ) exists for all s > 0 and
k ^ ( s ) SF , lim s + s k ^ ( s ) = + ,
where SF denotes the class of Stieltjes functions. The precise definition of this class is given in the next section, see (6). The basic particular examples of operators D t ( k ) , defined in this way, are the first-order derivative (in this case, k ^ ( s ) = 1 ), the Riemann–Liouville fractional derivative of order α ( 0 , 1 ) (corresponding to k ^ ( s ) = s α 1 ), as well as linear combinations with positive coefficients of such derivatives.
Let us note that assumptions (5) are slightly weaker than those required in the definition of the general fractional derivative, introduced in [19] and studied in detail in [20,22], where additional assumptions on the limiting behavior of k ^ ( s ) are imposed. We also point out that the assumption k ^ ( s ) SF , which is essential also in the definition of the general fractional derivative in [19], allows the use of the convenient Bernstein functions technique developed in [27].
In this work, based on the Bernstein functions technique, we establish the complete monotonicity of the relaxation modulus for model (3), provided the constraint a b is satisfied. This result is stronger than the thermodynamic compatibility results proven in [26], since the dissipation inequality is satisfied for systems with completely monotone relaxation moduli. The complete monotonicity property implies also the subordination principle for the generalized fractional Zener wave equation. The propagation function is obtained as the solution of an equation governing wave propagation in a viscoelastic half-space with the generalized fractional Zener model (3) and its properties are studied. The relation between the propagation function and the subordination kernel is established. Some examples of application of subordination principle are given. Numerical results for specific cases of generalized derivatives are presented and discussed.
The rest of the paper is organized as follows. Section 2 contains preliminaries and properties of Bernstein functions and related classes of functions, which are necessary for the study. In Section 3, we establish the complete monotonicity property of the relaxation modulus and discuss the relation of this property with thermodynamic compatibility conditions. Propagation of waves in generalized fractional Zener medium is discussed in Section 4. A subordination principle is formulated in Section 5 and applied to two basic problems. Numerical examples for some specific kernels are discussed in Section 6. The last Section 7 contains concluding remarks.

2. Preliminaries

This section contains preliminaries on the following functions classes: completely monotone functions ( CMF ), Stieltjes functions ( SF ), Bernstein functions ( BF ) and complete Bernstein function ( CBF ). Their definitions and properties, necessary for our proofs, are given in the terminology of [27].
We use the following notations for the Laplace transform of a function
L { f ( t ) } ( s ) = f ^ ( s ) = 0 e s t f ( t ) d t .
The Laplace transform of a function u ( x , t ) with respect to the time variable t, where x is considered as a parameter, is denoted by u ^ ( x , s ) or L { u } ( x , s ) .
We proceed with the definitions of the four special classes of functions. The definition of the class CMF of completely monotone functions is already given in the Section 1, see (2). The characterization of the class CMF is given by the Bernstein’s theorem: a function is completely monotone if and only if it can be represented as the Laplace transform of a non-negative (generalized) function.
The class BF of Bernstein functions consists of all non-negative functions f ( t ) defined and infinitely differentiable on R + , such that f ( t ) CMF .
The class of Stieltjes functions ( SF ) consists of all functions φ defined on R + which admit the representation (see [19])
φ ( s ) = φ 0 + 0 e s τ f ( τ ) d τ , s > 0 ,
where φ 0 0 , f CMF L l o c 1 ( R + ) , and the Laplace transform of f exists for any s > 0 . The constant φ 0 can be determined from the identity (see, e.g., [28], Theorem 2.6) φ 0 = lim s + φ ( s ) .
The class CBF of complete Bernstein functions consists of all function φ defined on R + , such that φ ( s ) / s SF , s > 0 .
The simplest examples of functions from the classes CBF and SF are given as:
if   α [ 0 , 1 ]   then   s α CBF   and   s α SF .
The following properties are satisfied for s > 0 :
(P1)
The set CMF is a convex cone: λ 1 f 1 + λ 2 f 2 CMF for all f 1 , f 2 CMF and λ 1 , λ 2 0 . Moreover, f 1 · f 2 CMF for all f 1 , f 2 CMF . The set CMF is closed under pointwise limits;
(P2)
The sets SF , BF and CBF are convex cones. They are closed under pointwise limits;
(P3)
Any Stieltjes function is completely monotone. Any complete Bernstein function is a Bernstein function;
(P4)
Let φ BF . Then, the function φ ( s ) / s is completely monotone;
(P5)
CMF BF CMF , where ∘ denotes composition of functions from the corresponding sets;
(P6)
Every function φ SF , such that lim s + φ ( s ) = 0 , is the Laplace transform of a function f ( t ) , which is locally integrable on R + and completely monotone;
(P7)
φ SF if and only if s φ ( s ) CBF ;
(P8)
Let φ 0 . Then, φ ( s ) CBF if and only if ( φ ( s ) ) 1 SF ;
(P9)
Let φ , ψ CBF . Assume μ 1 , μ 2 ( 0 , 1 ) are such that μ 1 + μ 2 1 . Then,
φ μ 1 ( s ) · ψ μ 2 ( s ) CBF ;
(P10)
Any function φ from the sets CBF or SF admits an analytic extension to the complex plane cut along the negative real axis C \ ( , 0 ] , which satisfies
| arg φ ( z ) | | arg z | , z C \ ( , 0 ] .
Moreover, φ ( z ) * = φ ( z * ) , where * denotes the complex conjugate, and
φ ( z ) · z 0 for φ SF ; φ ( z ) · z 0 for φ CBF .
For details on the proofs, see [27,28], Theorem 2.6.
The Mittag–Leffler function is defined as follows:
E α , β ( z ) = k = 0 z k Γ ( α k + β ) ,
where α > 0 , β R and z C . If λ > 0 and 0 < α β 1 , then the function of Mittag–Leffler type t β 1 E α , β ( λ t α ) is completely monotone for t > 0 . The Laplace transform of this function is given by the identity
L t β 1 E α , β ( λ t α ) ( s ) = s α β s α + λ , t > 0 .

3. Completely Monotone Relaxation Modulus

Consider the generalized Zener constitutive Equation (3). We assume that ε ( x , t ) = 0 and σ ( x , t ) = 0 for all x R and t < 0 . Moreover, we suppose that ε C 1 ( R + ) .
Let us note that the assumption k ^ ( s ) SF is equivalent to the following representation of the kernel k ( t )
k ( t ) = k 0 δ ( t ) + k 1 ( t ) ,
where k 0 0 , δ ( · ) denotes the Dirac delta function and k 1 ( t ) L l o c 1 ( R + ) is a completely monotone function.
The relaxation modulus G ( t ) in a linear viscoelastic model is defined as the stress response to a unit step of strain. Specifically, for a quiescent system at t = 0 , the following equation is satisfied [2]:
σ ( x , t ) = 0 t G ( t τ ) ε τ ( x , τ ) d τ .
Applying formally the Laplace transform to (3) and taking into account the Laplace transform identity L { D t ( k ) u } ( s ) = s k ^ ( s ) , we deduce
1 + a s k ^ ( s ) σ ^ ( x , s ) = 1 + b s k ^ ( s ) ε ^ ( x , s ) .
This, together with the identity σ ^ ( x , s ) = G ^ ( s ) s ε ^ ( x , s ) (which follows from (11)), yields
G ^ ( s ) = 1 + b s k ^ ( s ) s ( 1 + a s k ^ ( s ) ) .
We next discuss the properties of the relaxation modulus G ( t ) , based on the Laplace transform pair (12).
In a thermodynamically compatible model, the function G ( t ) is non-negative and non-increasing. We prove next that these properties can be satisfied only if a b . Indeed, since the rate of relaxation is non-increasing, G ( t ) 0 , it follows L { G } ( s ) 0 for s > 0 . Therefore, taking into account identity (12), we obtain for s > 0
L d G d t ( s ) = s G ^ ( s ) lim t 0 + G ( t ) = 1 + b s k ^ ( s ) 1 + a s k ^ ( s ) b a = a b a ( 1 + a s k ^ ( s ) ) 0 .
Here, we have used that the initial value theorem for the Laplace transform implies
lim t 0 + G ( t ) = lim s + s G ^ ( s ) = lim s + 1 + b s k ^ ( s ) 1 + a s k ^ ( s ) = b a
due to the assumption s k ^ ( s ) + as s + . The inequality in (13) implies a b , since the denominator is a positive function. In this way, we have proved that condition a b is necessary for viscoelastic model (3) to be physically acceptable.
We prove next that the restriction a b together with the assumptions (5) on the kernel k ( t ) implies the complete monotonicity of the relaxation modulus for viscoelastic model (3).
Theorem 1.
Assume that 0 < a b and that conditions (5) on the kernel k ( t ) are satisfied. Then, there exists a locally integrable and completely monotone function G ( t ) , t 0 , which Laplace transform is given by the identity (12).
Proof. 
For the proof, we use the Bernstein functions technique. For s > 0 , it holds
s G ^ ( s ) = 1 + b s k ^ ( s ) 1 + a s k ^ ( s ) = 1 + ( b a ) 1 ( s k ^ ( s ) ) 1 + a CBF .
Here, we use the assumption k ^ ( s ) SF , which is equivalent to ( s k ^ ( s ) ) 1 SF (see properties (P7) and (P8)). Since a > 0 , this implies ( s k ^ ( s ) ) 1 + a SF . By (P8), this is equivalent to ( s k ^ ( s ) ) 1 + a 1 CBF , which, by noting that b a 0 , yields s G ^ ( s ) CBF . Applying again (P7), it follows G ^ ( s ) SF . Moreover, (12) and assumption s k ^ ( s ) + as s + imply lim s + G ^ ( s ) = 0 . Therefore, according to (P6), there exists a locally integrable function G ( t ) CMF , which Laplace transform is given by (12). □
By the uniqueness of Laplace transform, it follows that the function G ( t ) from the above theorem coincides with the relaxation modulus of constitutive model (3). In this way, we proved that the generalized Zener constitutive Equation (3) defines a viscoelastic system with a completely monotonic relaxation modulus. It is known (see, e.g., [8]) that systems with completely monotonic relaxation moduli obey the dissipation inequality, that is, for any x R and T > 0 , the following holds:
0 T σ ( x , t ) ε t ( x , t ) d t 0 .
Corollary 1.
Under the conditions of Theorem 1, the dissipation inequality (16) is satisfied for model (3).
Proof. 
For completeness, we give a direct proof of (16). Taking into account (11), the dissipation inequality (16) is equivalent to
0 T 0 t G ( t τ ) ε τ ( x , τ ) d τ ε t ( x , t ) d t 0 .
This is satisfied when G ( t ) is a function of positive type, that is,
0 T 0 t G ( t τ ) u ( x , τ ) d τ u ( x , t ) d t 0 for every u L 2 ( 0 , T ) .
We will use the fact that every completely positive function is a function of positive type, see, e.g., [29], Theorem 3.1. A function G ( t ) is called completely positive if it satisfies the following conditions: G L 1 ( 0 , T ) and there exist γ 0 and a non-negative and non-increasing function K L 1 ( 0 , T ) , such that
γ G ( t ) + 0 t K ( t τ ) G ( τ ) d τ = 1 , t [ 0 , T ] .
It is known that any completely monotone function is completely positive, see, e.g., [29]. For completeness, we will check directly that the relaxation modulus G ( t ) for the generalized Zener model (3) is a completely positive function. To this end, let us consider the function
F ( s ) = b a b ( 1 + b s k ^ ( s ) ) , s > 0 .
Since s k ^ ( s ) is a complete Bernstein function and b a 0 , property (P8) implies F ( s ) SF . Moreover, lim s F ( s ) = 0 . Therefore, according to (P6), there exists a locally integrable function K ( t ) CMF , such that K ^ ( s ) = F ( s ) . Then, identity (19) is satisfied with K ( t ) defined in this way and γ = a / b . This can be easily checked by applying Laplace transform since
γ G ^ ( s ) + K ^ ( s ) G ^ ( s ) = 1 / s .
Therefore, the dissipation inequality (16) is satisfied for model (3). □
Let us note that conditions for thermodynamic compatibility are considered in different form in [7]. In our notation, these conditions read
{ i ω G ^ ( i ω ) } 0 , { i ω G ^ ( i ω ) } 0 , ω > 0 .
For completely monotone functions G ( t ) , conditions (20) are always fulfilled. This follows from the fact that in this case, G ^ ( s ) SF , which is equivalent to s G ^ ( s ) CBF (see (P7)), and then properties (P10) imply (20).
Let us point out that since s k ^ ( s ) CBF , the limit lim s 0 + s k ^ ( s ) exists and is real (see [27], Theorem 6.2), which implies
lim t + G ( t ) = lim s 0 + s G ^ ( s ) = lim s 0 + 1 + b s k ^ ( s ) 1 + a s k ^ ( s ) > 0 .
This limit confirms that the class of models (3) govern solid-like behavior, see, e.g., [2].

4. Generalized Fractional Zener Wave Equation

The one-dimensional equation governing wave propagation in a viscoelastic medium with the generalized Zener constitutive law is derived by combining constitutive Equation (3), the one-dimensional equation of motion
x σ ( x , t ) + f ( x , t ) = 2 t 2 u ( x , t )
and the kinematic equation
ε ( x , t ) = x u ( x , t ) ,
where u ( x , t ) denotes the particle displacement, f ( x , t ) is the body force, x R , and t > 0 . Combining the three equations, we derive formally the following generalized wave equation of Zener type
1 + a D t ( k ) 2 t 2 u ( x , t ) = 1 + b D t ( k ) 2 x 2 u ( x , t ) + 1 + a D t ( k ) f ( x , t ) .
It is worth noting that in the limiting case a = b , Equation (24) is the classical wave equation, which is well studied. Therefore, in what follows, we consider only coefficients satisfying
0 < a < b .
Assuming the initial conditions
lim t 0 + u ( x , t ) = v ( x ) and lim t 0 + t u ( x , t ) = 0 ,
we rewrite the wave Equation (24) as a Volterra integral equation as follows. Applying Laplace transform to (24) and taking into account the initial conditions (25), we obtain
u ^ ( x , s ) = v ( x ) 1 s + 1 + b s k ^ ( s ) s 2 ( 1 + a s k ^ ( s ) ) L 2 u x 2 ( x , s ) + 1 s 2 f ^ ( x , s ) .
Denote by ξ ( t ) the function with Laplace transform
ξ ^ ( s ) = 1 + b s k ^ ( s ) s 2 ( 1 + a s k ^ ( s ) ) .
According to (12), the function ξ ( t ) is well defined and
ξ ( t ) = 0 t G ( τ ) d τ , t 0 .
Taking the inverse Laplace transform in (26), we derive the following Volterra integral equation
u ( x , t ) = v ( x ) + 0 t ξ ( t τ ) 2 x 2 u ( x , τ ) d τ + h ( x , t ) ,
where h ( x , t ) = 0 t ( t τ ) f ( x , τ ) d τ . Let us note that G ( t ) CMF implies ξ ( t ) BF .
In what follows, we study Equation (28) as a model for wave propagation in a viscoelastic medium with the generalized Zener constitutive law (3). We consider different possibilities for the domain of the spatial operator.

4.1. Propagation Function

Let us consider first the basic problem of wave propagation of impact waves in a viscoelastic half-space with constitutive Equation (3):
w ( x , t ) = 0 t ξ ( t τ ) 2 x 2 w ( x , τ ) d τ , x , t > 0 ,
w ( 0 , t ) = H ( t ) , w 0 a s x , t > 0 ,
where the kernel ξ is defined in (27) and H ( t ) is the Heaviside unit step function. The solution w ( x , t ) is often referred to as propagation function, since it represents the propagation in time of a disturbance at x = 0 .
It is worth noting that in the limiting case a = b , corresponding to the classical wave equation, the propagation function is w ( x , t ) = H ( t x ) .
Let us also note that the initial conditions
lim t 0 + w ( x , t ) = lim t 0 + t w ( x , t ) = 0 ,
follow from Equation (29), taking into account representation (27) for the kernel ξ ( t ) .
By applying Laplace transform with respect to the temporal variable in (29) and (30), we obtain the following ODE in x (with s considered as a parameter):
g ( s ) w ^ ( x , s ) = 2 w ^ x 2 ( x , s ) , w ^ ( 0 , s ) = 1 / s , w ^ ( x , s ) 0 a s x ,
where
g ( s ) = 1 ξ ^ ( s ) = s 2 1 + a s k ^ ( s ) 1 + b s k ^ ( s ) .
The solution of problem (32) is given by
w ^ ( x , s ) = 1 s exp x g ( s ) .
The following property is important in the study of the behavior of function w ( x , t ) .
Proposition 1.
If 0 < a b and k ( t ) satisfies conditions (5), then for the characteristic function g ( s ) , defined in (33), it holds
g ( s ) CBF .
Proof. 
Applying (33) and (12), we obtain the relation
g ( s ) = s G ^ ( s ) .
Since s CBF and 1 / G ^ ( s ) CBF (here, we use G ^ ( s ) SF , and property (P8)), the function g ( s ) is a product of two functions from the set CBF . Then, (35) follows by applying property (P9) with μ 1 = μ 2 = 1 / 2 . □
According to (35) g ( s ) CBF BF . Then, property (P5) implies
exp x g ( s ) CMF
with respect to s > 0 (and x > 0 considered as a parameter). Therefore, the function w ^ is a product of two functions, completely monotone for s > 0 , which yields w ^ ( x , s ) CMF . This implies by Bernstein’s theorem w ( x , t ) 0 . Moreover, the propagation function w ( x , t ) is non-decreasing in t and non-increasing in x. This can be established in the same way as Theorem 2.2 in [30] by proving that the Laplace transforms of the functions w / t and w / x are completely monotone functions. Indeed, since 1 / s CMF , and g ( s ) CBF BF , by property (P4) also g ( s ) / s CMF . This, together with (37), yields
L w x ( x , s ) = x w ^ ( x , s ) = g ( s ) s exp x g ( s ) CMF
as products of two completely monotone functions. Further, since lim t 0 + w ( x , t ) = 0 , it follows from (34) and (37)
L w t ( x , s ) = s w ^ ( x , s ) lim t 0 + w ( x , t ) = exp x g ( s ) CMF .

4.2. Propagation Speed

Applying the general formula for the propagation speed c of a disturbance in a viscoelastic medium [2,31], we obtain for the considered model (3) with conditions (5)
c = lim s s g ( s ) = b a .
This can be seen by inverting the Laplace transform in representation (34). Therefore, the propagation speed is finite for model (3) and w ( x , t ) 0 for x > t b / a .
From general theory (see, e.g., [31], Section 5.4) in the case of finite propagation speed c, a jump discontinuity of the propagation function w ( x , t ) at the planar surface x = c t exists only if
κ = lim s g ( s ) s / c < .
Moreover, this discontinuity is damped exponentially with exponent κ , i.e., w ( x , t ) exp ( κ x ) as t x / c , and x < c t . In our specific case, it holds that
κ = lim s s 1 + a s k ^ ( s ) 1 + b s k ^ ( s ) 1 / 2 a b 1 / 2 = lim s a b 1 / 2 s 1 + ( a s k ^ ( s ) ) 1 1 + ( b s k ^ ( s ) ) 1 1 / 2 1 .
We find the leading terms in the last expression as s by taking into account the expansion ( 1 + z ) γ 1 + γ z as z 0 , γ R , for z = ( a s k ^ ( s ) ) 1 and z = ( b s k ^ ( s ) ) 1 . In this way, we deduce
1 + ( a s k ^ ( s ) ) 1 1 + ( b s k ^ ( s ) ) 1 1 / 2 1 + 1 2 a s k ^ ( s ) 1 2 b s k ^ ( s ) , s .
Therefore,
κ = lim s a b 1 / 2 1 a 1 b 1 2 k ^ ( s ) .
Therefore, a jump at the wave front exists only when lim s k ^ ( s ) > 0 , i.e., only when k 0 0 in the kernel representation (10). This is, for instance, the case with the composite derivative D t ( k ) = d d t + λ D t α , where D t α denotes the Riemann–Liouville fractional derivative of order α ( 0 , 1 ) and λ 0 . For λ = 0 , the classical Zener model is recovered. On the other hand, if lim s k ^ ( s ) = 0 ( k 0 = 0 in the kernel representation (10) and thus, k ( t ) = k 1 ( t ) CMF ), the wave front is smooth. In this case, the coexistence of finite wave speed and absence of jump discontinuity at the wave front is observed.
In the next theorem, we summarize the most important findings in this section concerning the propagation function.
Theorem 2.
Assume that 0 < a b and conditions (5) on the kernel k ( t ) are satisfied. Then, the solution w ( x , t ) of problem (29) and (30) is a non-negative function, non-increasing in x and non-decreasing in t . Moreover, w ( x , t ) 0 for x > t b / a . If lim s k ^ ( s ) > 0 , a jump discontinuity is observed at x = t b / a , otherwise, the function is smooth.
Plots of the propagation functions w ( x , t ) for specific kernels at different values of the parameters are given in Figure 1 and Figure 2. The observed behavior in the figures is in agreement with the established properties of the propagation function.

5. Subordination Principle

Subordination results for different types of generalized wave equations have been established recently, see, e.g., [30,32,33,34]. Solution representations in the context of the principle of subordination are useful, since they provide an integral decomposition of the solution of such equations in terms of the solutions of two simpler problems.
In this section, we formulate the subordination principle for equations governing wave propagation in media with the generalized Zener constitutive law (3). For the sake of generality, we formulate the results in an abstract setting, using the terminology of [31].

5.1. General Formulation

Let A be a closed operator on a Banach space X. Assume A is a generator of a strongly continuous cosine family of operators Cos ( t ) . This means that the second-order abstract Cauchy problem
d 2 u d t 2 = A u ( t ) , t > 0 ; u ( 0 ) = v X , lim t 0 d u d t = 0 .
is well posed, with solution operator Cos ( t ) . The Laplace transform of the cosine family Cos ( t ) , generated by the operator A, is given by the identity
0 e s t Cos ( t ) d t = s ( s 2 A ) 1 ,
where ( z A ) 1 denotes the resolvent operator of the operator A.
The problem for the abstract Volterra integral equation
u ( t ) = 0 t ξ ( t τ ) A u ( τ ) d τ + f ( t ) , t > 0 ,
where ξ L l o c 1 ( R + ) is a scalar kernel and f-an X-valued function, is said to be well posed if and only if there exists a solution operator S ( t ) , i.e., a strongly continuous family of bounded linear operators in X which commutes with A and such that the function u ( t ) = S ( t ) v satisfies the equation
u ( t ) = v + 0 t ξ ( t τ ) A u ( τ ) d τ , t > 0 , v X .
In this case, the solution of problem (43) is found by the variation of parameters formula
u ( t ) = d d t 0 t S ( t τ ) f ( τ ) d τ .
The solution operator S ( t ) obeys the Laplace transform identity
0 e s t S ( t ) d t = 1 s ( 1 ξ ^ ( s ) A ) 1 .
For detailed theory of abstract Volterra integral equations, we refer to [31].
Let ξ ( t ) be the kernel defined by (27). The corresponding function g ( s ) is given in (33). Based on the property (35), the following subordination theorem can be established in an analogous way as Theorem 3.1 in [30].
Theorem 3.
Let 0 < a < b and the kernel k ( t ) satisfy conditions (5). Let operator A generate a bounded cosine family Cos ( t ) . Assume the kernel ξ ( t ) is defined by (27). Then, problem (43) is well posed and the corresponding solution operator S ( t ) is related to Cos ( t ) as follows
S ( t ) = 0 t b / a ϕ ( t , τ ) Cos ( τ ) d τ , t > 0 .
The subordination function ϕ ( t , τ ) satisfies the Laplace transform identity
ϕ ^ ( s , τ ) = g ( s ) s exp τ g ( s ) ,
where the function g ( s ) is defined in (33). Moreover, ϕ ( t , τ ) 0 , ϕ ( t , τ ) = 0 for τ > t b / a , and 0 ϕ ( t , τ ) d τ = 1 .
Proof. 
For completeness, we give a sketch of the proof. The well-posedness of the Volterra equation follows by applying Theorem 4.3(ii) in [31]. According to this theorem, it suffices to prove that ξ ( t ) is a creep function (i.e., it is non-negative, non-decreasing and concave) and that the function ξ 1 ( t ) = ξ ˙ ( t ) lim t ξ ( t ) / t is log-convex. Indeed, representation (27), together with complete monotonicity of G ( t ) , implies that ξ ( t ) is a Bernstein function, therefore, it is a creep function. Moreover, (27) yields ξ 1 ( t ) = G ( t ) G ( + ) . This is a completely monotone function, which implies that ξ 1 ( t ) is log-convex (see, e.g., [28], Theorem 2.8).
Let us check the next integral decomposition (47). Assume function ϕ is defined as in (48). Application of Laplace transform implies by the use of (42)
0 e s t 0 ϕ ( t , τ ) Cos ( τ ) d τ d t = 0 ϕ ^ ( s , τ ) Cos ( τ ) d τ = g ( s ) s 0 exp τ g ( s ) Cos ( τ ) d τ = g ( s ) s ( g ( s ) A ) 1 .
Comparing (46)–(49), it follows by the uniqueness of the Laplace transform that
S ( t ) = 0 ϕ ( t , τ ) Cos ( τ ) d τ , t > 0 .
Let us note that (34) and (48) imply the relation
ϕ ( t , τ ) = x w ( x , t ) x = τ ,
where w ( x , t ) is the propagation function. Therefore, ϕ ( t , τ ) = 0 for τ > t b / a and (50) implies (47). We already proved the non-negativity of w / x in the previous section, which implies the non-negativity of ϕ ( t , τ ) . Moreover, (51) implies
0 ϕ ( t , τ ) d τ = 0 w x ( x , t ) d x = w ( t , 0 ) w ( t , ) = 1 ,
where the boundary conditions (30) are taken into account. □

5.2. Applications

To demonstrate the applicability of Theorem 3, two examples are given next.
The first example is the following Cauchy problem for the generalized Zener-type wave equation
u ( x , t ) = v ( x ) + 0 t ξ ( t τ ) 2 x 2 u ( x , τ ) d τ , x R , t > 0 ,
lim | x | u ( x , t ) = 0 , t > 0 ,
where the kernel ξ is defined in (27).
Let X = L p ( R ) , 1 p < . Define the operator A by means of ( A u ) ( x ) = u ( x ) , with domain D ( A ) = { u X : u , u X , u ( ± ) = 0 } . Then, A generates a bounded cosine family given by the d’Alembert formula
( Cos ( t ) v ) ( x ) = 1 2 v ( x + t ) + v ( x t ) .
Inserting (54) in the subordination Formula (47), we obtain for the solution to problem (52) and (53)
u ( x , t ) = ( S ( t ) v ) ( x ) = 0 c t ϕ ( t , τ ) ( Cos ( τ ) v ) ( x ) d τ = 1 2 c t c t ϕ ( t , | σ | ) v ( x σ ) d σ ,
where ϕ is the subordination function from Theorem 3 and c = b / a .
As a second example, let us consider the inhomogeneous Zener wave equation in integral form (28) on a finite interval, x [ 0 , 1 ] , subject to boundary and initial conditions
u ( 0 , t ) = u ( 1 , t ) = 0 , t > 0 , u 0 ( x ) = 0 , x [ 0 , 1 ] .
Let us set X = L 2 ( [ 0 , 1 ] ) and define the operator A by ( A u ) ( x ) = u ( x ) , x [ 0 , 1 ] with domain D ( A ) = { u X : u , u X , u ( 0 ) = u ( 1 ) = 0 } . The cosine family Cos ( t ) is defined in this case by the solution of the following problem for the classical wave equation
2 u t 2 = 2 u x 2 , 0 < x < 1 , t > 0 , u ( 0 , t ) = u ( 1 , t ) = 0 , t > 0 ; u ( x , 0 + ) = v ( x ) , lim t 0 + u t = 0 , 0 < x < 1 .
Therefore, if v ( x ) admits the eigenfunction expansion v ( x ) = n = 1 v n sin ( n π x ) then
( Cos ( t ) v ) ( x ) = n = 1 v n sin ( n π x ) cos ( n π t ) .
Let f ( x , t ) admits the expansion f ( x , t ) = n = 1 f n ( t ) sin ( n π x ) . From the subordination identity (47) and the variation of parameters formula (45), we obtain for the solution of the Zener-type wave Equation (28) with initial-boundary conditions (56)
u ( x , t ) = n = 1 sin ( n π x ) 0 t S n ( t τ ) f n ( τ ) d τ ,
where
S n ( t ) = 0 t 0 c τ ϕ ( τ , σ ) cos ( n π σ ) d σ d τ .
Here, ϕ is the subordination function from Theorem 3 and c = b / a .

6. Numerical Results

In this section, we present numerical examples. The function w ( x , t ) and the subordination kernel ϕ ( t , τ ) are computed for three types of kernels k ( t ) . The numerical computations are based on integral representations of these functions, obtained as follows.
Integral representation for the propagation function w ( x , t ) can be derived by inverting its Laplace transform (34). Following the method of our previous work [30], Theorem 2.5, which uses solely property (35) and g ( s ) + as s + (and not the specific form of g ( s ) ), we have the following integral representation for the propagation function:
w ( x , t ) = 1 2 + 1 π 0 exp ( x K + ( r ) ) sin ( r t x K ( r ) ) d r r , x , t > 0 ,
where
K ± ( r ) = 1 2 A 2 ( r ) + B 2 ( r ) 1 / 2 ± A ( r ) 1 / 2 , A ( r ) = g ( i r ) , B ( r ) = g ( i r ) .
In the case of characteristic function g ( s ) , given in (33), we obtain for the functions A ( r ) and B ( r )
A ( r ) = r 2 ( 1 + a p ( r ) ) ( 1 + b p ( r ) ) + a b q ( r ) 2 ( 1 + b p ( r ) ) 2 + ( b q ( r ) ) 2 , B ( r ) = r 2 ( b a ) q ( r ) ( 1 + b p ( r ) ) 2 + ( b q ( r ) ) 2 ,
where p ( r ) = { s k ^ ( s ) } and q ( r ) = { s k ^ ( s ) } , for s = i r .
Further, the relation (51) and the integral representation (59) imply the following integral representation for ϕ ( t , τ ) :
ϕ ( t , τ ) = 1 π 0 exp τ K + ( r ) K + ( r ) sin r t τ K ( r ) + K ( r ) cos r t τ K ( r ) d r r , t , τ > 0 .
Here, K ± ( r ) are the functions defined in (60) and (61).
Next, we list some examples of kernels k ( t ) in the generalized derivative D t ( k ) , which satisfy conditions (5). These conditions guarantee that the results obtained in this work for constitutive model (3) with 0 < a < b are satisfied. The specific kernels are used in numerical experiments, the results of which are discussed in this section.
Since the numerical experiments are based on the integral representations (59) for the propagation function w ( x , t ) and (62) for the subordination function ϕ ( t , τ ) , for each kernel, we provide also the corresponding Laplace transforms and the expressions for the functions p ( r ) and q ( r ) , which appear in (61).
For the sake of brevity, we use the notation
ω α ( t ) = t α 1 Γ ( α ) , α > 0 , t > 0 .
If α = 0 , we set ω 0 ( t ) = δ ( t ) , t > 0 , where δ ( t ) is the Dirac delta function.
Example 1.
Two-term composite derivative: D t ( k ) = D t α + λ D t β , where D t α and D t β are Riemann–Liouville fractional derivatives with orders satisfying 0 < β < α 1 and λ 0 . In this case,
k ( t ) = ω 1 α ( t ) + λ ω 1 β ( t ) , k ^ ( s ) = s α 1 + λ s β 1 , 0 < β < α 1 , λ 0 ,
and
p ( r ) = r α cos ( α π / 2 ) + λ r β cos ( β π / 2 ) , q ( r ) = r α sin ( α π / 2 ) + λ r β sin ( β π / 2 ) .
Example 2.
Generalized derivative D t ( k ) with a singular Mittag–Leffler-type kernel:
k ( t ) = t β 1 E α , β ( t α ) , k ^ ( s ) = s α β s α + 1 , 0 < α β < 1 .
In this case,
p ( r ) = r 2 α β + 1 sin ( β π / 2 ) + r α β + 1 sin ( ( β α ) π / 2 ) r 2 α + 2 r α cos ( α π / 2 ) + 1 , q ( r ) = r 2 α β + 1 cos ( β π / 2 ) + r α β + 1 cos ( ( β α ) π / 2 ) r 2 α + 2 r α cos ( α π / 2 ) + 1 .
Example 3.
Distributed-order fractional derivative with uniform distribution over the interval [ 0 , 1 ] . In this case,
k ( t ) = 0 1 ω 1 α ( t ) d α , k ^ ( s ) = 0 1 s α 1 d α = s 1 s ln s ,
and
p ( r ) = π r / 2 ln r ln 2 r + π 2 / 4 , q ( r ) = r ln r + π / 2 ln 2 r + π 2 / 4 .
Another example of kernel, satisfying conditions (5), is the following
k ( t ) = ω 1 α ( t ) e μ t , k ^ ( s ) = ( s + μ ) α 1 ,
where 0 < α < 1 and μ > 0 . In fact, in [26], the generalized fractional Zener model is studied for two types of kernels: (65) and (67). Generalizations of the Mittag–Leffler kernel (65) can also be considered, such as the completely monotone multinomial Mittag–Leffler-type functions, see, e.g., [35,36].
In the present work, we treat numerically only the cases in Examples 1–3. All mentioned kernels admit decomposition (10) and their Laplace transforms satisfy conditions (5).
In Figure 1 and Figure 2, plots of propagation functions are given, while Figure 3 presents plots of subordination functions for different kernels k ( t ) . The numerical results provide a numerical verification of the analytical findings in this paper.
In all figures, except Figure 2b, the parameters a and b are taken a = 0.2 , b = 1.0 , which corresponds to a constant propagation speed, c = b / a = 5 . Therefore, in these figures, w ( x , t ) vanishes for x > t 5 and ϕ ( t , τ ) vanishes for τ > t 5 .
In Figure 1a, the propagation function w ( x , t ) is plotted for the kernel (64) with λ = 0 , which corresponds to the fractional Zener model (1) when α ( 0 , 1 ) . The function w ( x , t ) is plotted as a function of x for t = 0.4 and for different values of the parameter α . The limiting case α = 1 corresponds to the classical integer-order Zener model, where a jump discontinuity is observed at the wave front. The step function given with thinner line corresponds to α 0 , recovering the classical elastic wave equation. For α ( 0 , 1 ) , we observe, as expected, the coexistence of a finite wave speed and smooth wave front.
In Figure 1b, the propagation function w ( x , t ) is plotted for the composite kernel (64) with α = 1 , λ = 1 and β = 0.5 . The function w ( x , t ) is plotted as a function of x for different times t. Since α = 1 , it holds κ > 0 , see (40). Therefore, a jump discontinuity at the wave front is expected and observed. The thinner line is the graph of the function exp ( κ x ) and it marks the height of the jump.
Figure 2a corresponds to the Mittag–Leffler-type kernel (65) with α = 0.1 and different values of β , t = 0.4 . Figure 2b corresponds to the distributed-order kernel (66), a takes different values in the interval ( 0 , 1 ) , b = 1 , t = 0.3 . In Figure 2, the propagation function w ( x , t ) is plotted as a function of x. Since lim s k ^ ( s ) = 0 for the considered kernels, the propagation function is smooth at the wave front.
In Figure 3, the subordination function ϕ is plotted as a function of τ for two types of kernels. Figure 3a corresponds to kernel (64) with α = 0.75 and β = 0.25 , λ = 1 . Different values of t are considered. Figure 3b corresponds to the Mittag–Leffler-type kernel (65) with α = 0.3 and different β , α β < 1 ; t = 0.4 . The plotted functions in Figure 3 are probability densities as functions of τ .

7. Concluding Remarks

In this work, we study the Zener-type model of a viscoelastic body in the context of general fractional calculus, proposed recently in [26]. The model is studied for a general class of kernels by applying the Bernstein functions technique. This technique allows us to prove that the thermodynamic constraint a < b , used in [26], implies not only the dissipation inequality, but also guarantees the complete monotonicity property of the relaxation modulus G ( t ) , which is a stronger result. Based on this property, the propagation function, which is characteristic for the wave propagation in a viscoelastic medium, is studied and the subordination principle is established for wave equations with generalized Zener-type constitutive law. The subordination identity provides an integral decomposition of the solution into two parts. The first part is a probability density, related to the propagation function and depending only on the constitutive model. The second part is the solution of a related classical wave equation and it depends only on the problem geometry. The analytical findings are supported by numerical work.

Author Contributions

Conceptualization, E.B. and I.B.; methodology and formal analysis, E.B.; software and visualization, I.B.; writing—original draft preparation and revision, E.B. and I.B. All authors have read and agreed to the published version of the manuscript.

Funding

The first author (E.B.) is supported by Grant No BG05M2OP001-1.001-0003, financed by the Science and Education for Smart Growth Operational Program and co-financed by the European Union through the European structural and Investment funds. The second author (I.B.) is supported by the Bulgarian National Science Fund under Grant FNI KP-06-H22/2.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mainardi, F. An historical perspective on fractional calculus in linear viscoelasticity. Fract. Calc. Appl. Anal. 2012, 15, 712–717. [Google Scholar] [CrossRef] [Green Version]
  2. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity, 2nd ed.; World Scientific (Europe): London, UK, 2022. [Google Scholar] [CrossRef]
  3. Atanacković, T.M.; Pilipović, S.; Stanković, B.; Zorica, D. Fractional Calculus with Applications in Mechanics: Vibrations and Diffusion Processes; John Wiley & Sons: London, UK, 2014. [Google Scholar]
  4. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
  5. Caputo, M.; Mainardi, F. A new dissipation model based on memory mechanism. Pure Appl. Geophys. 1971, 91, 134–147. [Google Scholar] [CrossRef]
  6. Caputo, M.; Mainardi, F. Linear models of dissipation in anelastic solids. Riv. Nuovo Cimento (Ser. II) 1971, 1, 161–198. [Google Scholar] [CrossRef]
  7. Bagley, R.L.; Torvik, P.J. On the fractional calculus model of viscoelastic behavior. J. Rheol. 1986, 30, 137–148. [Google Scholar] [CrossRef]
  8. Hanyga, A. Physically acceptable viscoelastic models. In Trends in Applications of Mathematics to Mechanics; Hutter, K., Wang, Y., Eds.; Shaker Verlag GmbH: Aachen, Germany, 2005; pp. 125–136. [Google Scholar]
  9. Hanyga, A. Wave propagation in linear viscoelastic media with completely monotonic relaxation moduli. Wave Motion 2013, 50, 909–928. [Google Scholar] [CrossRef] [Green Version]
  10. Konjik, S.; Oparnica, L.; Zorica, D. Waves in fractional Zener type viscoelastic media. J. Math. Anal. Appl. 2010, 365, 259–268. [Google Scholar] [CrossRef] [Green Version]
  11. Atanacković, T.M. On a distributed derivative model of a viscoelastic body. C. R. Méc. 2003, 331, 687–692. [Google Scholar] [CrossRef]
  12. Rossikhin, Y.; Shitikova, M. Analysis of rheological equations involving more than one fractional parameter by the use of the simplest mechanical systems based on these equations. Mech. Time-Depend. Mat. 2001, 5, 131–175. [Google Scholar] [CrossRef]
  13. Rossikhin, Y.; Shitikova, M. Analysis of the viscoelastic rod dynamics via models involving fractional derivatives or operators of two different orders. Shock Vib. Digest 2004, 36, 3–26. [Google Scholar] [CrossRef]
  14. Atanacković, T.; Konjik, S.; Oparnica, L.; Zorica, D. Thermodynamical restrictions and wave propagation for a class of fractional order viscoelastic rods. Abstr. Appl. Anal. 2011, 2011, 975694. [Google Scholar] [CrossRef] [Green Version]
  15. Konjik, S.; Oparnica, L.; Zorica, D. Distributed-order fractional constitutive stress–strain relation in wave propagation modeling. Z. Angew. Math. Phys. 2019, 70, 51. [Google Scholar] [CrossRef] [Green Version]
  16. Atanacković, T.M.; Janev, M.; Pilipović, S. Restrictions on parameters in distributed order fractional linear constitutive equations. Appl. Math. Model. 2022, 110, 99–111. [Google Scholar] [CrossRef]
  17. Bazhlekova, E.; Bazhlekov, I. Complete monotonicity of the relaxation moduli of distributed-order fractional Zener model. AIP Conf. Proc. 2018, 2048, 050008. [Google Scholar] [CrossRef]
  18. 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: Vienna, Austria, 1997; pp. 223–276. [Google Scholar]
  19. Kochubei, A. General fractional calculus, evolution equations, and renewal processes. Integr. Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  20. Luchko, Y.; Yamamoto, M. The general fractional derivative and related fractional differential equations. Mathematics 2020, 8, 2115. [Google Scholar] [CrossRef]
  21. Sandev, T.; Tomovski, Z.; Dubbeldam, J.; Chechkin, A. Generalized diffusion-wave equation with memory kernel. J. Phys. A Math. Theor. 2019, 52, 015201. [Google Scholar] [CrossRef] [Green Version]
  22. Luchko, Y. Operational calculus for the general fractional derivative and its applications. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  23. Górska, K.; Horzela, A.; Lenzi, E.; Pagnini, G.; Sandev, T. Generalized Cattaneo (telegrapher’s) equations in modeling anomalous diffusion phenomena. Phys. Rev. E 2020, 102, 023719. [Google Scholar] [CrossRef]
  24. Dos Santos, M. Non-Gaussian distributions to random walk in the context of memory kernels. Fractal Fract. 2018, 2, 20. [Google Scholar] [CrossRef] [Green Version]
  25. Bazhlekov, I.; Bazhlekova, E. A predictor-corrector numerical approach to equations with general fractional derivative. Int. J. Appl. Math. 2022, 35, 693–709. [Google Scholar] [CrossRef]
  26. Atanacković, T.M.; Pilipović, S. Zener model with General Fractional Calculus: Thermodynamical restrictions. Fractal Fract. 2022, 6, 617. [Google Scholar] [CrossRef]
  27. Schilling, R.L.; Song, R.; Vondraček, Z. Bernstein Functions: Theory and Applications; De Gruyter: Berlin, Germany, 2010. [Google Scholar]
  28. Gripenberg, G.; Londen, S.O.; Staffans, O. Volterra Integral and Functional Equations; Cambridge University Press: Cambridge, UK, 1990. [Google Scholar]
  29. Clément, P.; Mitidieri, E. Qualitative properties of solutions of Volterra equations in Banach spaces. ISR J. Math. 1988, 64, 1–24. [Google Scholar] [CrossRef]
  30. Bazhlekova, E.; Bazhlekov, I. Subordination approach to multi-term time-fractional diffusion-wave equations. J. Comput. Appl. Math. 2018, 339, 179–192. [Google Scholar] [CrossRef] [Green Version]
  31. Prüss, J. Evolutionary Integral Equations and Applications; Birkhäuser: Basel, Switzerland, 1993. [Google Scholar]
  32. Bazhlekova, E. Subordination in a class of generalized time-fractional diffusion-wave equations. Fract. Calc. Appl. Anal. 2018, 21, 869–900. [Google Scholar] [CrossRef]
  33. Bazhlekova, E.; Bazhlekov, I. Transition from diffusion to wave propagation in fractional Jeffreys-type heat conduction equation. Fractal Fract. 2020, 4, 32. [Google Scholar] [CrossRef]
  34. Górska, K. Integral decomposition for the solutions of the generalized Cattaneo equation. Phys. Rev. E 2021, 104, 024113. [Google Scholar] [CrossRef]
  35. Bazhlekova, E. Completely monotone multinomial Mittag-Leffler type functions and diffusion equations with multiple time-derivatives. Fract. Calc. Appl. Anal. 2021, 24, 88–111. [Google Scholar] [CrossRef]
  36. Bazhlekova, E.; Pshenichnov, S. Wave propagation in viscoelastic half-space with memory functions of Mittag-Leffler type. Int. J. Appl. Math. 2021, 34, 423–440. [Google Scholar] [CrossRef]
Figure 1. Propagation function w ( x , t ) as a function of x for a = 0.2 , b = 1.0 and kernel (64); (a) λ = 0 , t = 0.4 and different values of α ; marked with a thinner line—the limiting elastic case α 0 ; (b) λ = 1 , α = 1 and β = 0.5 , marked with a thinner line—the function exp ( κ x ) .
Figure 1. Propagation function w ( x , t ) as a function of x for a = 0.2 , b = 1.0 and kernel (64); (a) λ = 0 , t = 0.4 and different values of α ; marked with a thinner line—the limiting elastic case α 0 ; (b) λ = 1 , α = 1 and β = 0.5 , marked with a thinner line—the function exp ( κ x ) .
Fractalfract 07 00298 g001
Figure 2. Propagation function w ( x , t ) as a function of x for (a) a = 0.2 , b = 1.0 and Mittag–Leffler-type kernel (65) with α = 0.1 and different values of β , t = 0.4 ; (b) distributed-order kernel (66), different values of a ( 0 , 1 ) , b = 1 , t = 0.3 .
Figure 2. Propagation function w ( x , t ) as a function of x for (a) a = 0.2 , b = 1.0 and Mittag–Leffler-type kernel (65) with α = 0.1 and different values of β , t = 0.4 ; (b) distributed-order kernel (66), different values of a ( 0 , 1 ) , b = 1 , t = 0.3 .
Fractalfract 07 00298 g002
Figure 3. Subordination kernel ϕ ( t , τ ) as a function of τ for a = 0.2 , b = 1 and (a) kernel (64) with α = 0.75 , β = 0.25 , λ = 1 , and different values of t; (b) the Mittag–Leffler-type kernel (65) with α = 0.3 and different β , t = 0.4 .
Figure 3. Subordination kernel ϕ ( t , τ ) as a function of τ for a = 0.2 , b = 1 and (a) kernel (64) with α = 0.75 , β = 0.25 , λ = 1 , and different values of t; (b) the Mittag–Leffler-type kernel (65) with α = 0.3 and different β , t = 0.4 .
Fractalfract 07 00298 g003
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

Bazhlekova, E.; Bazhlekov, I. Subordination Principle for Generalized Fractional Zener Models. Fractal Fract. 2023, 7, 298. https://doi.org/10.3390/fractalfract7040298

AMA Style

Bazhlekova E, Bazhlekov I. Subordination Principle for Generalized Fractional Zener Models. Fractal and Fractional. 2023; 7(4):298. https://doi.org/10.3390/fractalfract7040298

Chicago/Turabian Style

Bazhlekova, Emilia, and Ivan Bazhlekov. 2023. "Subordination Principle for Generalized Fractional Zener Models" Fractal and Fractional 7, no. 4: 298. https://doi.org/10.3390/fractalfract7040298

APA Style

Bazhlekova, E., & Bazhlekov, I. (2023). Subordination Principle for Generalized Fractional Zener Models. Fractal and Fractional, 7(4), 298. https://doi.org/10.3390/fractalfract7040298

Article Metrics

Back to TopTop