Next Article in Journal / Special Issue
Numerical Solution of Fractional Differential Equations: A Survey and a Software Tutorial
Previous Article in Journal
Acknowledgement to Reviewers of Mathematics in 2017
Previous Article in Special Issue
Best Approximation of the Fractional Semi-Derivative Operator by Exponential Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Storage and Dissipation of Energy in Prabhakar Viscoelasticity

1
Department of Information and Communication Technologies, Universitat Pompeu Fabra, C/Roc Boronat 138, 08018 Barcelona, Spain
2
Department of Physics & Astronomy, University of Bologna, Via Irnerio 46, 40126 Bologna, Italy
3
I.N.F.N., Sezione di Bologna, via B. Pichat 6/2, I-40127 Bologna, Italy
4
Arnold Sommerfeld Center, Ludwig-Maximilians-Universität, Theresienstraße 37, 80333 München, Germany
*
Author to whom correspondence should be addressed.
Mathematics 2018, 6(2), 15; https://doi.org/10.3390/math6020015
Submission received: 12 December 2017 / Revised: 11 January 2018 / Accepted: 14 January 2018 / Published: 23 January 2018
(This article belongs to the Special Issue Fractional Calculus: Theory and Applications)

Abstract

:
In this paper, after a brief review of the physical notion of quality factor in viscoelasticity, we present a complete discussion of the attenuation processes emerging in the Maxwell–Prabhakar model, recently developed by Giusti and Colombaro. Then, taking profit of some illuminating plots, we discuss some potential connections between the presented model and the modern mathematical modelling of seismic processes.

1. Introduction

The linear theory of viscoelasticity, despite its apparent simplicity, keeps having a striking impact in geophysics, theoretical mechanics and biophysics; see, for example, [1,2,3,4,5,6,7]. Besides, fractional calculus [8,9,10] has proven itself to be one of the fundamental languages for describing processes involving memory effects, like the ones that are typically featured by viscoelastic systems. Concerning the latter, it is also worth remarking on the pivotal role of the notion of complete monotonicity, which was first (implicitly) hinted at by Gross [11] in 1953, and then brought to light by Molinari [12], in 1973. These seminal studies were then followed by many other authors; see, for example, [10,13,14].
A simple generalization of the well-known fractional Maxwell model of linear viscoelasticity was first introduced by Giusti and Colombaro in [15]. In this paper the authors provide an extension of the classical model by replacing the Caputo fractional derivative with the Prabhakar one in the constitutive equation. Concretely, if we denote with σ , ε , the stress and the strain for a given system, respectively, and we further assume that these functions are both causal such that σ , ε A C 1 0 , + , then the constitutive equation of the Maxwell–Prabhakar model [15] reads
σ ( t ) + a C D α , β , Ω γ σ ( t ) = b C D α , β , Ω γ ε ( t ) ,
where a and b are two suitable real constants and provided that α , β , γ , Ω R , α > 0 and 0 < β < 1 .
Here, C D α , β , Ω γ represents the regularized Prabhakar derivative [16], which is defined by
a C D α , β , Ω γ f ( t ) = a E α , m β , Ω γ f ( m ) ( t ) ,
where
a E μ , ν , λ ρ f ( t ) = a t ( t τ ) ν 1 E μ , ν ρ λ ( t τ ) μ f ( τ ) d τ
denotes the Prabhakar fractional integral [17] and E μ , ν ρ ( t ) represents the Prabhakar function [18,19,20,21]. Furthermore, it is important to stress that the function t ν 1 E μ , ν ρ ( t μ ) , for t > 0 , is locally integrable and completely monotone provided that 0 < μ 1 and 0 < μ ρ ν 1 ; see, for example, [22,23].
It is important to remark that the Prabhakar fractional calculus has been attracting much attention in the mathematical community [15,17,18,23,24,25,26], particularly because of its connection with the theoretical description of the Havriliak–Negami model [18,24,27,28]. Moreover, this growing interest in Prabhakar’s calculus is also reflected by the increasing literature on the recently proposed Maxwell–Prabhakar model, which was also kindly referred to as the Giusti–Colombaro model in [29].
In this paper we wish to analyze the important phenomena of storage and dissipation of energy in linear viscoelastic media, with particular regard for the class of models emerging from the constitutive equation in Equation (1). In viscoelasticity, as well as in electrical engineering, the process of dissipation of energy is usually accounted for in terms of a dimensionless parameter, called the quality factor, that is roughly defined as the ratio of the peak of energy stored in the system under a cycle of forced harmonic oscillation to the total rate of change of the energy, per cycle, by damping processes. Therefore, the aim of this paper is to compute and discuss the quality factor for the model defined in Equation (1).

2. Storage and Dissipation of Energy in Linear Viscoelasticity

In this section we wish to review the general theory, concerning the theoretical foundations, that leads to the definition of quality factor for a viscoelastic system. In order to do so, we will mimic the arguments presented in [10,30], unifying these formulations according to the notations employed in this paper.
Let us consider a quiescent viscoelastic body for t < 0 . Then, under the hypothesis of sufficiently well-behaved causal histories, its constitutive equation in the creep representation reads
ε ( t ) = 0 t J ( t τ ) d σ ( τ ) = σ ( 0 + ) J ( t ) + 0 t J ( t τ ) σ ˙ ( τ ) d τ ,
where d σ ( τ ) represents the Riemann–Stieltjes measure and J ( t ) is the so-called creep compliance of the system, that in the Laplace domain is given by
ε ˜ ( s ) = s J ˜ ( s ) σ ˜ ( s ) .
In order to consider the harmonic behavior of a linear viscoelastic material, we should assume that a sufficient amount of time has elapsed since the original perturbation so that the effect of initial conditions could be considered negligible. So, let us consider some harmonic excitation of the material, which can be described in terms of the complex exponential representation, that is
σ ( t ; ω ) = χ exp i ω t , ω > 0 , < t < , χ C .
Clearly, a similar argument can be presented in terms of the relaxation representation, however we will only focus on the creep one for sake of brevity.
If we now plug (5) into (3) we get
ε ( t ; ω ) = i ω J ^ ( ω ) ( χ exp i ω t ) = i ω J ^ ( ω ) σ ( t ; ω ) ,
where J ^ ( ω ) stands for the Fourier transform of J ( t ) that, the latter being a causal function, ultimately reads
J ^ ( ω ) = 0 exp ( i ω t ) J ( t ) d t .
Moreover, if we denote J ( ω ) = i ω J ^ ( ω ) , as in [10], the constitutive equation, in the creep representation, for a viscoelastic body subject to a harmonic stress excitation reduces to
ε ( t ; ω ) = J ( ω ) σ ( t ; ω ) .
The time rate of change of energy in the system is then given by
W ˙ = σ R ( t ; ω ) ε ˙ R ( t ; ω ) ,
where the subscript R indicates that we are considering the real part of the corresponding function.
Now, one can easily solve (7) for σ ( t ; ω ) , then taking the real part of the resulting equation gives
σ R ( t ; ω ) = ε R ( t ; ω ) J R ( ω ) ε I ( t ; ω ) J I ( ω ) | J ( ω ) | 2 ,
where we denoted σ σ R + i σ I , ε ε R + i ε I and J ( ω ) J R ( ω ) i J I ( ω ) for future convenience. Besides, from (5) and (7) it is also easy to see that
ε I ( t ; ω ) = ε ˙ R ( t ; ω ) ω .
Hence, if we plug (9) into (8) and recall the result in (10), then after some simple manipulations W ˙ reduces to
W ˙ = t 1 2 J R ( ω ) | J ( ω ) | 2 ε R 2 + 1 ω J I ( ω ) | J ( ω ) | 2 ε ˙ R 2 ,
where we omitted the explicit dependence on t and ω in the strain for sake of clarity.
Thus, it is easy to see that the total rate of change of energy over one cycle is accounted for by the integral over the cycle of the second term on the right-hand side of (11), namely
Δ E Cycle = t t + T W ˙ ( τ ) d τ = π J I ( ω ) | J ( ω ) | 2 | χ | 2 ,
with T = 2 π / ω the period of the cycle.
Due to the second law of thermodynamics, which requires that the total amount of energy dissipated increases with time, one can further infer that J I ( ω ) 0 .
However, despite defining a boundary term, the first piece of the right-hand side of (11) carries a very important physical meaning. Indeed, it tells us that the peak energy stored during a cycle is given by
P max = 1 2 J R ( ω ) | J ( ω ) | 2 | χ | 2 .
One can now define the specific attenuation factor, or quality factor (Q-factor), as a normalized non-dimensional quantity defined by
Q 1 1 2 π Δ E / Cycle P max .
Then, taking profit of the previous discussion it is easy to see that
Q 1 = 1 2 π Δ E / Cycle P max = J I ( ω ) J R ( ω ) = i ω J ^ ( ω ) i ω J ^ ( ω ) ,
recalling that J I ( ω ) = J ( ω ) .
If we combine the fact that the Fourier transform is equivalent to evaluating the bilateral Laplace transform with imaginary argument s = i ω , together with the assumption for which J ( t ) is a causal function, one can conclude that J ^ ( ω ) = J ˜ ( s ) | s = i ω . This argument ultimately leads us to a very useful expression for the Q-factor, namely
Q 1 ( ω ) = s J ˜ ( s ) | s = i ω s J ˜ ( s ) | s = i ω ,
where we shall consider some positive real frequencies ω .

3. Quality Factor in Prabhakar-Like Viscoelasticity

Let us now compute the Q-factor for the Maxwell–Prabhakar model. Recalling that the Laplace transform of the Prabhakar integral kernel is given by
L t β 1 E α , β γ ( λ t α ) = s β 1 λ s α γ ,
where t R , α , β , γ , λ C and Re ( β ) > 0 , then it is easy to see that the creep compliance, in the Laplace domain, for a system described in terms of Equation (3) is therefore given by
s J ˜ ( s ) = a b + 1 b s β 1 Ω s α γ .
Then, if we apply the replacement s = i ω , the latter turns into
J ( ω ) = i ω J ^ ( ω ) = a b + 1 b ( i ω ) β 1 Ω ( i ω ) α γ .
Let us define an auxiliary variable z ( ω ) = 1 Ω ( i ω ) α . Then, considering ω R + , one can easily rewrite s = i ω = ω exp ( i π / 2 ) , where | s | = ω , that allows us to recast this new variable in the exponential representation, that is
z ( ω ) = | z ( ω ) | exp i θ z ( ω ) ,
with
| z ( ω ) | 2 = 1 + Ω 2 ω 2 α 2 Ω ω α cos α π 2 ,
θ z ( ω ) = arctan Ω sin α π / 2 ω α Ω cos α π / 2 .
Then, plugging z ( ω ) into Equation (18), one can easily infer that
J ( ω ) = a b + 1 b ω β | z ( ω ) | γ exp i β π 2 + γ θ z ( ω ) ,
from which we can conclude that
J ( ω ) = a b + 1 b ω β | z ( ω ) | γ cos β π 2 + γ θ z ( ω ) ,
and
J ( ω ) = 1 b ω β | z ( ω ) | γ sin β π 2 + γ θ z ( ω ) .
Hence, the quality factor for a Maxwell–Prabhakar viscoelastic body is given by
Q 1 ( ω ) = sin β π 2 + γ θ z ( ω ) a ω β | z ( ω ) | γ + cos β π 2 + γ θ z ( ω ) .

4. Quality Factor for Some Specific Realizations of the Maxwell–Prabhakar Model

In this section we discuss the quality factor for different choices of the parameters of the discussed model. Specifically, we will focus our discussion on four cases corresponding to two well-known classical viscoelastic models and the viscoelastic analogue of the Havriliak–Negami model for dielectric relaxation.

4.1. Fractional Maxwell Model

As argued in [15], it is easy to see that Equation (1) naturally reduces to the fractional Maxwell model, namely
σ ( t ) + a C D ν σ ( t ) = b C D ν ε ( t ) ,
where C D ν represents the Caputo fractional derivative, provided that the parameters are chosen according to one of these two configurations:
(i)
γ = 0 , a > 0 , b > 0 , β = ν , Ω R ;
(ii)
γ R , a > 0 , b > 0 , β = ν , Ω = 0 .
Here, it is trivial to infer that if γ = 0 then neither | z ( ω ) | nor θ z ( ω ) enter in the expression for the Q-factor, whereas if Ω = 0 it is easy to see that | z ( ω ) | = 1 and θ z ( ω ) = 0 . Hence, we get
Q 1 ( ω ) = sin ν π / 2 a ω ν + cos ν π / 2 .
Furthermore, it is also worth remarking that if we set ν = 1 (ordinary limit), we explicitly recover the Q-factor for the (ordinary) Maxwell model, that is, Q 1 ( ω ) = ( a ω ) 1 (see, for example, [10]).

4.2. Fractional Voigt Model

Again, following the analysis presented in [15], one has that Equation (1) reduces to the fractional Voigt model, that is
σ ( t ) = M ε ( t ) + B C D ν ε ( t ) ,
by setting γ = 1 , a = 0 , b = B < 0 , α = β = ν , Ω = M / B < 0 .
Now, if we plug this choice of parameters into Equation (26), we get
Q 1 ( ω ) = sin ν π 2 + θ z ( ω ) cos ν π 2 + θ z ( ω ) = tan ν π 2 + θ z ( ω ) ,
with
θ z ( ω ) = arctan M B sin ν π / 2 ω ν + ( M / B ) cos ν π / 2 .
The latter, in the limit for ν = 1 , reduces to
Q 1 ( ω ) = tan π 2 + θ z ( ω ) = 1 tan θ z ( ω ) = B ω M .
which corresponds to the quality factor for the (ordinary) Voigt model (see, for example, [10]).

4.3. Havriliak–Negami Model

The Havriliak–Negami relaxation is an empirical model which was first introduced in order to describe the dielectric relaxation of certain types of polymers [24,25,31].
Now, it is very well known that a viscoelastic system can usually be mapped onto a class of electrical ladder networks and vice versa; see, for example, [32,33].
Following this line of thought, it is easy to see that the constitutive equation for the Havriliak–Negami viscoelastic model is given by [15]
σ ( t ) + a C D α , α γ , λ γ σ ( t ) = b C D α , α γ , λ γ ε ( t ) ,
with β = α γ , Ω = λ , with λ > 0 , and 0 < α , γ < 1 .
Then, following the procedure presented in Section 3, one finds that
Q 1 ( ω ) = sin γ α π 2 + θ z ( ω ) a ω α γ ω 2 α + λ 2 + 2 λ ω α cos α π 2 γ / 2 + cos γ α π 2 + θ z ( ω ) ,
where
θ z ( ω ) = arctan λ sin α π / 2 ω α + λ cos α π / 2 .

5. Discussion and Conclusions

Attenuation effects represent one of the main fields of study in modern seismology, and consequently the specific attenuation factor (or Q-factor, for simplicity) embodies one of the key ingredients in geophysical sciences. Indeed, were it not for the damping capabilities of the soil, the energy of past earthquakes would still be resonating within the earth’s interior.
In this paper, after a thorough review of the physical definition and meaning of the quality factor Q, we have provided an analysis of the storage and dissipation of energy in viscoelastic material of the Maxwell–Prabhakar class, namely, the one featured by a constitutive relation given by Equation (1). Specifically, in Section 3 we have computed the Q-factor for the general Maxwell–Prabhakar model, for which some interesting configurations of the parameters are shown in Figure 1. In Section 4 we have further provided some explicit realizations of the Maxwell–Prabhakar theory, namely the fractional Maxwell model, the fractional Voigt model and the viscoelastic equivalent of the Havriliak–Negami model for dielectric relaxation, shown in Figure 2 and Figure 3.
Let us pay particular attention to the cases displayed in Figure 1. Indeed, as argued in [34], there exists much experimental evidence supporting the theses for which the Q-factor of homogeneous materials is substantially independent of the frequency. In this respect, it is worth noting that the Maxwell–Prabhakar class shows a very slow varying (almost constant) behaviour of the quality factor for low frequencies, for certain choices of the parameters of the model (see Figure 1). This is quite consistent with the results for the dumping of long-period teleseismic body waves and surface waves. Furthermore, for high frequencies, the model shows a power law behaviour, namely, Q ( ω ) ω β as ω , which appears to be consistent with the experimental results concerning the attenuation of the coda of high-frequency teleseismic waves in Earth’s upper mantle [35]. Furthermore, it is rather easy to prove that some very well-known constant-Q models are nothing but some specific realizations of the model defined in Equation (1). Indeed, for example, the renowned Kjartansson model [36,37] can be obtained from (18) by setting γ = 0 , β 2 η with η ( 0 , 1 / 2 ) and a = 0 , which is nothing but the Scott–Blair model [36]. In view of the last few comments, we believe that the Maxwell–Prabhakar model of viscoelasicity can potentially provide some stimulating new insights into the mathematical modelling of seismic processes and therefore is worthy of further studies.

Acknowledgments

The work of I.C. and A.G. has been carried out in the framework of the activities of the National Group of Mathematical Physics (GNFM, INdAM). Moreover, the work of A.G. has been partially supported by GNFM/INdAM Young Researchers Project 2017 “Analysis of Complex Biological Systems”. Besides, the work of S.V. has been partially supported by the Interdepartmental Center “Luigi Galvani” for integrated studies of Bioinformatics, Biophysics and Biocomplexity of the University of Bologna.

Author Contributions

The authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Colombaro, I.; Giusti, A.; Mainardi, F. On transient waves in linear viscoelasticity. Wave Motion 2017, 74, 191–212. [Google Scholar] [CrossRef]
  2. Colombaro, I.; Giusti, A.; Mainardi, F. On the propagation of transient waves in a viscoelastic Bessel medium. Z. Angew. Math. Phys. 2017, 68, 62–74. [Google Scholar] [CrossRef]
  3. Colombaro, I.; Giusti, A.; Mainardi, F. A class of linear viscoelastic models based on Bessel functions. Meccanica 2017, 52, 825–832. [Google Scholar] [CrossRef]
  4. Garra, R.; Mainardi, F.; Spada, G. A generalization of the Lomnitz logarithmic creep law via Hadamard fractional calculus. Chaos Solitons Fractals 2017, 102, 333–338. [Google Scholar] [CrossRef]
  5. Giusti, A. On infinite order differential operators in fractional viscoelasticity. Fract. Calc. Appl. Anal. 2017, 20, 854–867. [Google Scholar] [CrossRef]
  6. Giusti, A.; Mainardi, F. A dynamic viscoelastic analogy for fluid-filled elastic tubes. Mecanica 2016, 51, 2321–2330. [Google Scholar] [CrossRef]
  7. Mainardi, F. Fractional Calculus: Some basic problems in continuum and statistical mechanics. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: New York, NY, USA; Wien, Austria, 1997. [Google Scholar]
  8. Giusti, A. A comment on some new definitions of fractional derivative. arXiv, 2017; arXiv:1710.06852. [Google Scholar]
  9. 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: New York, NY, USA; Wien, Austria, 1997. [Google Scholar]
  10. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity; Imperial College Press: London, UK, 2010. [Google Scholar]
  11. Gross, B. Mathematical Structure of the Theories of Viscoelasticity; Hermann & Cie: Paris, France, 1953. [Google Scholar]
  12. Molinari, A. Viscoélasticité linéaire et functions complétement monotones. Journal de Mécanique 1973, 12, 541–553. [Google Scholar]
  13. Mainardi, F.; Turchetti, G. Positivity constraints and approximation methods in linear viscoelasticity. Lettere al Nuovo Cimento 1979, 26, 38–40. [Google Scholar] [CrossRef]
  14. Hanyga, A. Wave propagation in linear viscoelastic media with completely monotonic relaxation moduli. Wave Motion 2013, 50, 909–928. [Google Scholar] [CrossRef]
  15. Giusti, A.; Colombaro, I. Prabhakar-like fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 138–143. [Google Scholar] [CrossRef]
  16. D’Ovidio, M.; Polito, F. Fractional Diffusion-Telegraph Equations and their Associated Stochastic Solutions. Teoriya Veroyatnostei i ee Primeneniya 2017, 62, 692–718. [Google Scholar] [CrossRef]
  17. Garra, R.; Gorenflo, R.; Polito, F.; Tomovski, Z. Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef]
  18. Garra, R.; Garrappa, R. The Prabhakar or three parameter Mittag-Leffler function: Theory and application. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 314–329. [Google Scholar] [CrossRef]
  19. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  20. Paneva-Konovska, J. From Bessel to Multi-Index Mittag Leffler Functions: Enumerable Families, Series in Them and Convergence; World Scientific Publishing: London, UK, 2016. [Google Scholar]
  21. 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]
  22. Capelas de Oliveira, E.; Mainardi, F.; Vaz, J., Jr. Models based on Mittag-Leffler functions for anomalous relaxation in dielectrics. Eur. Phys. J. Spec. Top. 2011, 193, 161–171. [Google Scholar] [CrossRef]
  23. Mainardi, F.; Garrappa, R. On complete monotonicity of the Prabhakar function and non-Debye relaxation in dielectrics. J. Comput. Phys. 2015, 293, 70–80. [Google Scholar] [CrossRef]
  24. Garrappa, R. Grünwald-Letnikov operators for fractional relaxation in Havriliak-Negami models. Commun. Nonlinear Sci. Numer. Simul. 2016, 38, 178–191. [Google Scholar] [CrossRef]
  25. Garrappa, R.; Mainardi, F.; Maione, G. Models of dielectric relaxation based on completely monotone functions. Fract. Calc. Appl. Anal. 2016, 19, 1105–1160. [Google Scholar] [CrossRef]
  26. Sandev, T. Generalized Langevin equation and the Prabhakar derivative. Mathematics 2017, 5, 66. [Google Scholar] [CrossRef]
  27. Hanyga, A.; Seredyńska, M. On a Mathematical Framework for the Constitutive Equations of Anisotropic Dielectric Relaxation. J. Stat. Phys. 2008, 131, 269–303. [Google Scholar] [CrossRef]
  28. Seredyńska, M.; Hanyga, A. Relaxation, dispersion, attenuation, and finite propagation speed in viscoelastic media. J. Math. Phys. 2010, 51, 092901. [Google Scholar] [CrossRef]
  29. Ding, X.; Zhang, G.; Zhao, B.; Wang, Y. Unexpected viscoelastic deformation of tight sandstone: Insights and predictions from the fractional Maxwell model. Sci. Rep. 2017, 7, 11336. [Google Scholar] [CrossRef] [PubMed]
  30. Borcherdt, R. Viscoelastic Waves in Layered Media; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  31. Havriliak, S.; Negami, S. A complex plane representation of dielectric and mechanical relaxation processes in some polymers. Polymer 1967, 8, 161–210. [Google Scholar] [CrossRef]
  32. Giusti, A.; Mainardi, F. On infinite series concerning zeros of Bessel functions of the first kind. Eur. Phys. J. Plus. 2016, 131, 206–212. [Google Scholar] [CrossRef]
  33. Gross, B.; Fuoss, R. Ladder structures for representation of viscoelastic systems. J. Polym. Sci. 1956, 19, 39–50. [Google Scholar] [CrossRef]
  34. Knopoff, L. Q. Rev. Geophys. 1964, 2, 625–660. [Google Scholar] [CrossRef]
  35. Shito, A.; Karato, A.; Park, J. Frequency dependence of Q in Earth’s upper mantle inferred from continuous spectra of body waves. Geophys. Res. Lett. 2004, 31, L12603. [Google Scholar] [CrossRef]
  36. Carcione, J.; Cavallini, F.; Mainardi, F.; Hanyga, A. Time-domain modeling of constant-Q Seismic waves using fractional derivatives. Pure Appl. Geophys. 2002, 159, 1719–1736. [Google Scholar] [CrossRef]
  37. Kjartansson, E. Constant Q-wave propagation and attenuation. J. Geophys. Res. Solid Earth 1979, 84, 4737–4748. [Google Scholar] [CrossRef]
Figure 1. Q-factor of the Maxwell–Prabhakar model for a = 0.01 , Ω = 10 . Model dependence on α = 0.1 , 0.2 , , 1 , β = 0.5 , γ = 0.3 (a); model dependence on β = 0.1 , 0.2 , , 1 , α = 0.3 , γ = 0.3 (b); model dependence on γ = 0.1 , 0.2 , , 1 , β = 0.5 , α = 0.3 (c). Parameter value increasing from blue to red. In all panels the dashed line corresponds to α = 0.3 , β = 0.5 , γ = 0.3 .
Figure 1. Q-factor of the Maxwell–Prabhakar model for a = 0.01 , Ω = 10 . Model dependence on α = 0.1 , 0.2 , , 1 , β = 0.5 , γ = 0.3 (a); model dependence on β = 0.1 , 0.2 , , 1 , α = 0.3 , γ = 0.3 (b); model dependence on γ = 0.1 , 0.2 , , 1 , β = 0.5 , α = 0.3 (c). Parameter value increasing from blue to red. In all panels the dashed line corresponds to α = 0.3 , β = 0.5 , γ = 0.3 .
Mathematics 06 00015 g001
Figure 2. Q-factor of fractional Maxwell model with a = 0.1 (a), fractional Voigt model with M = 10 3 , B = 1 (b). Parameter value increasing from blue to red.
Figure 2. Q-factor of fractional Maxwell model with a = 0.1 (a), fractional Voigt model with M = 10 3 , B = 1 (b). Parameter value increasing from blue to red.
Mathematics 06 00015 g002
Figure 3. Q-factor of Havriliak–Negami model for dielectric relaxation with a = 1 , λ = 1 and α = 0.1 , 0.2 , , 1 , γ = 0.5 (a), or a = 1 , λ = 1 and γ = 0.1 , 0.2 , , 1 , α = 0.5 (b). Parameter value increasing from blue to red. Decreasing a, the widths of the curves become wider, while increasing λ , the peak results are shifted to larger ω for lower α .
Figure 3. Q-factor of Havriliak–Negami model for dielectric relaxation with a = 1 , λ = 1 and α = 0.1 , 0.2 , , 1 , γ = 0.5 (a), or a = 1 , λ = 1 and γ = 0.1 , 0.2 , , 1 , α = 0.5 (b). Parameter value increasing from blue to red. Decreasing a, the widths of the curves become wider, while increasing λ , the peak results are shifted to larger ω for lower α .
Mathematics 06 00015 g003

Share and Cite

MDPI and ACS Style

Colombaro, I.; Giusti, A.; Vitali, S. Storage and Dissipation of Energy in Prabhakar Viscoelasticity. Mathematics 2018, 6, 15. https://doi.org/10.3390/math6020015

AMA Style

Colombaro I, Giusti A, Vitali S. Storage and Dissipation of Energy in Prabhakar Viscoelasticity. Mathematics. 2018; 6(2):15. https://doi.org/10.3390/math6020015

Chicago/Turabian Style

Colombaro, Ivano, Andrea Giusti, and Silvia Vitali. 2018. "Storage and Dissipation of Energy in Prabhakar Viscoelasticity" Mathematics 6, no. 2: 15. https://doi.org/10.3390/math6020015

APA Style

Colombaro, I., Giusti, A., & Vitali, S. (2018). Storage and Dissipation of Energy in Prabhakar Viscoelasticity. Mathematics, 6(2), 15. https://doi.org/10.3390/math6020015

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