On Λ -Fractional Viscoelastic Models

: Λ -Fractional Derivative ( Λ -FD) is a new groundbreaking Fractional Derivative (FD) introduced recently in mechanics. This derivative, along with Λ -Transform ( Λ -T), provides a reli-able alternative to fractional differential equations’ current solving. To put it straightforwardly, Λ -Fractional Derivative might be the only authentic non-local derivative that exists. In the present article, Λ -Fractional Derivative is used to describe the phenomenon of viscoelasticity, while the whole methodology is demonstrated meticulously. The fractional viscoelastic Zener model is studied, for relaxation as well as for creep. Interesting results are extracted and compared to other methodologies showing the value of the pre-mentioned method.


Introduction
Fractional Calculus (FC) is a modern and robust field of mathematics serving the need for non-local mathematical analysis. Eringen [1], in his book regarding nonlocal continuum field theories, explains the need for non-local mathematical analysis in the micro and nanoscale. More precisely, as Eringen wrote in his book: "There exists a large class of problems in classical physics and continuum mechanics (classical field theories) that fall outside their domain of applications. Fracture of solids, stress fields at the dislocation core and the tips of cracks, singularities present at the point of application of concentrated loads (forces, couples, heat, etc.), sharp corners and discontinuities in bodies, the failure in the prediction of short-wavelength behavior of elastic waves, and several decades of the viscosity rise in fluids flowing in microscopic channels are but a few major anomalies that defy classical treatment." Therefore, the need for developing unconventional methodologies for treating nonlocal phenomena is compulsory more than ever. Specifically, FC shows unique advantages when describing size effects, especially in materials with fractal characteristics or porous materials (Carpinteri et al. [2]). Moreover, there are a plethora of articles strongly connecting fractal geometry and Fractional Calculus (Butera et al. [3]), Tatom [4]). Thus, there are cases where Fractional Calculus seems irreplaceable compared to conventional non-local theories, especially when fractal structures are involved. We refer to physics (Hilfer [5], West et al. [6], Mechanics (Drapaca et al. [7], Di Paola et al. [8], Carpinteri et al. [9], Lazopoulos [10], Baleanu et al. [11], Sumelka [12], chemistry, biology (Magin [13]), control systems, economics, etc. Especially in mechanics, there are phenomena, such as viscoelasticity, where FC excels compared to other methodologies (Atanackovic et al. [14,15]). That is so because that mathematical theory is one of the few describing phenomena with fractal characteristics non-locally. Viscoelasticity is a privileged field in Fractional Calculus: It is the best-known phenomenon that Fractional Calculus applies, and a lot of fractional models were proposed over the years: Fractional Calculus models referred to Magin [13], Lazopoulos [16], Mainardi [17], Atanackovic [14] and Yang [18] are only a few of them. These models' problems are that they use the traditional approach; They replace in ODEs the conventional derivatives with fractional ones and solve the problem directly. This article will use a double-face derivative, the Λ-Derivative, and the dual Λ-Space to tackle this problem.
Moreover, Fractional Calculus is a suitable methodology for describing various phenomena, especially in the micro and nano regions where atomic interactions are strong; especially when a fractal structure is involved, FC is considered the only way forward. On the other hand, FDs not only represent the phenomena non-locally but also provide maximum flexibility when a mathematical model is designed. That holds because Fractional Derivatives are defined through the fractional-order γ (which can vary at will). Likewise, some papers enhance this flexibility even more, such as Olivar-Romero et al. [19], which uses fractional derivatives to describe a model that intermediates among wave, heat, and transport equations. Furthermore, Di Paola et al. [20] face situations in which the Boltzmann linear superposition principle does not apply in standard form since the fractional order is not constant in time.
Unfortunately, derivatives in Fractional Calculus are merely operators and not derivatives since none of them satisfies the criteria of the differential topology of a derivative, as described in Chillingworth [21]). Those criteria are: FC's great weakness was denounced single-mindedly in the past by various specialists in the field (Samko et al. [22], König et al. ([23,24]), Cresson et al. [25]). There were serious efforts to tackle this grave problem by designing fractional derivatives that follow Leibniz's rule (Jumarie [26,27], Yang [28]), without success. Therefore, as mentioned above, the problem is severe since it forbids the definition of a differential and, consequently, differential geometry.
It is repeatedly pointed out in every differential topology book (e.g., Gauld, D. [29] p. 95) that the Leibniz rule along with the chain rule is essential to define a differential, tangent space, and consequently, geometry.
Nevertheless, Fractional Calculus proceeded with describing phenomena in multiple scientific areas using derivatives that are operators and not derivatives, in the mathematical sense. Actually, in the book of Atangana [30], it is clearly stated in chapter one that " . . . We can therefore conclude that both the Riemann-Liouville and Caputo operators are not derivatives, and then they are not fractional derivatives, but fractional operators." This point of view is already shared by most of the fractional calculus community and emphasized on every occasion (e.g., conferences). Although that strategy offered many practical results, it lacked theoretical justification. The direct approach currently used to set up a Fractional Differential Equation (FDE) is to transform an Ordinary Differential Equation (ODE) by replacing the conventional derivatives with the appropriate fractional ones and solving it directly. Unfortunately, although posing an FDE is acceptable, solving it directly is utterly questionable since these FDs do not satisfy Equations (2) and (3); therefore, they cannot be considered proper derivatives at all.
On the other hand, Λ-Fractional Derivative is a newly defined derivative (Lazopoulos et al. [31]) that exhibits all the perquisites of a well-posed derivative in the dual Λ-Space (Λ-S). That mathematical concept introduces a transformation (Λ-Transformation (Λ-T)) that transfers the initial space (where the fractional problem is posed initially) to dual Λ-Space where the Λ-FD behaves as local, conventionally. Therefore, in Λ-Space, differentials corresponding to Λ-FDs exist, and fractional differential geometry may be generated. Solving the problem in the dual Λ-S, the results may be transferred back to the initial space. The beauty of this new mathematical concept is that it is local and non-local simultaneously, in two different Spaces: The initial Space (non-local) and Λ-Space (local). Thus Λ-Fractional Derivative has all the advantages of a classic local derivative (It is mainly consistent with all the requirements of differential Topology for a proper derivative), as well as the advantages of a non-local one (It is affected by the phenomenon non-locally). This central feature of Λ-FD provides great mathematical value to our study since this derivative might be the only proper derivative in the whole field of FC. Furthermore, Fractional Differential Equations (FDE) in the initial space are transformed in ODE in Λ-Space. This advantage is also revolutionary since no complex methodologies (i.e., application of Mittag-Leffler functions, etc.) are needed to solve these FDEs in the initial space.
In the present article, the basic concepts of Λ-Fractional Derivatives are introduced in the case of viscoelastic models (if the reader wants to know more on the application of Fractional Calculus on viscoelasticity then [18] by Yang et al. and [17] by Mainardi are considered the most valuable resources to start from). Afterward, the viscoelastic model's behavior is presented for the classical case (Zener's model). Then, using the Λ-Fractional derivative with the dual Λ-Space, Zener's model is scrutinized, and a specific example is shown. Finally, diagrams are presented, and conclusions are extracted.

Outline of Λ-Fractional Derivative, Λ-Transform, and Dual Λ-Space
Fractional Calculus is an interesting, demanding, and robust mathematical field serving the non-local mathematical analysis. Leibniz (who established the conventional local derivative) suggested whether a non-local derivative could be introduced. He proposed the derivative of half order in 1695. Many books have been written in that area ever since (Hilfer [5], Kilbas [32], Podlubny [33], Samko [22], Oldham [34]). Nevertheless, the methodology followed to solve Fractional Differential Equations is utterly questionable. The replacement of ordinary derivatives with fractional ones and solving them directly is a very fragile practice since the existing fractional derivatives fail to satisfy the Differential Topology requirements for a derivative. The question that arises then is how can we transform an ODE to an FDE and solve it more appropriately.
For an introduction to fractional analysis, the interested reader should be informed by references (Baleanu et al. [11], Sumelka [12], Magin [13], Atanackovic [14,15]). Nevertheless, a brief outline is presented in this chapter to make the present work more accessible.
Let Ω = [a,b] (−∞ < a < b < ∞) be a finite interval on the real axis R. The left and right Riemann-Liouville fractional integrals are defined by (Kilbas [32]): with γ (0 < γ ≤ 1) being the order of fractional integrals and Γ(n) = (n − 1)!, n integer. (Γ(γ) is called Euler's Gamma function). Fractional integrals are at our analysis's heart since the transformation of the equations from the initial to Λ-Space is based on them. Furthermore, since 0 < γ ≤ 1 applies, the Riemann-Liouville (RL) Fractional Derivatives are defined by (Kilbas [31]): where Equation (6) defines the left and Equation (7) the right fractional derivatives. Moreover, the fractional integrals with the corresponding Riemann-Liouville FDs are related by the equation: Riemann-Liouville Fractional Derivative is also essential to our methodology since Λ-Derivative is defined as the fraction of two such derivatives (see Lazopoulos [10,16] It is clear that RL a D γ x f (x) is the Riemann-Liouville Derivative of f(x), as described in FC (Equations (6) and (7)) and RL a I 1−γ x f (x) is the Riemann-Liouville fractional integral of real fractional dimension. In this article, 0 < γ ≤ 1 is considered. (See Samko et al. [7], Podlubny [33]).
Λ-Fractional Derivative is a new FD introduced by Lazopoulos et al. [31]. In that article, the Λ-Fractional Derivative and the dual Λ-Space and Λ-Transform, are defined. In this article, a variation of that method is proposed so that we can assure maximum mathematical validity.
When a problem is posed in Fractional Calculus, the first step is to find the appropriate Ordinary Differential Equation (ODE) for the phenomenon. The second step is to replace all integer derivatives with suitable fractional derivatives. Therefore, the fractional problem is posed. This is the procedure that almost every scientist in FC follows to solve problems in this field. Hence, up to this step, our methodology shows no difference. However, in contrast to mainstream methodology, we do not solve the Fractional Differential Equation (FDE) directly because fractional derivatives are not proper derivatives. So, we transform the FDE from the initial to Λ-Space in two steps.
In that space, the FDE posed in the initial space has taken the form of an ODE, undertaking all obstacles that FDs face: All perquisites from differential topology are satisfied, and there is a frame indifference, as well as all the advantages that an ODE possesses. We solve the equation in Λ-Space and afterward with the inverse Λ-Transformations Λ 1 −1 and Λ 2 −1 , we restore the solution F(X) of the equation in Λ-Space to f (x) in the initial Space ( Figure 1). It might be evident that Λ-Transform is following the structure of a Laplace transform. This is true, but the difference is that the Laplace transform only has two steps of evaluation, while the Λ-Transform has four ( Figure 1).
Riemann-Liouville Fractional Derivative is also essential to our methodology since Λ-Derivative is defined as the fraction of two such derivatives (see Lazopoulos [10,16]): It is clear that ( ) is the Riemann-Liouville Derivative of f(x), as described in FC (Equations (6) and (7)) and ( ) is the Riemann-Liouville fractional integral of real fractional dimension. In this article, 0 < γ ≤ 1 is considered. (See Samko et al. [7], Podlubny [33]). Λ-Fractional Derivative is a new FD introduced by Lazopoulos et al. [31]. In that article, the Λ-Fractional Derivative and the dual Λ-Space and Λ-Transform, are defined. In this article, a variation of that method is proposed so that we can assure maximum mathematical validity.
When a problem is posed in Fractional Calculus, the first step is to find the appropriate Ordinary Differential Equation (ODE) for the phenomenon. The second step is to replace all integer derivatives with suitable fractional derivatives. Therefore, the fractional problem is posed. This is the procedure that almost every scientist in FC follows to solve problems in this field. Hence, up to this step, our methodology shows no difference. However, in contrast to mainstream methodology, we do not solve the Fractional Differential Equation (FDE) directly because fractional derivatives are not proper derivatives. So, we transform the FDE from the initial to Λ-Space in two steps.
In that space, the FDE posed in the initial space has taken the form of an ODE, undertaking all obstacles that FDs face: All perquisites from differential topology are satisfied, and there is a frame indifference, as well as all the advantages that an ODE possesses. We solve the equation in Λ-Space and afterward with the inverse Λ-Transformations Λ1 −1 and Λ2 −1 , we restore the solution F(X) of the equation in Λ-Space to f(x) in the initial Space ( Figure 1). It might be evident that Λ-Transform is following the structure of a Laplace transform. This is true, but the difference is that the Laplace transform only has two steps of evaluation, while the Λ-Transform has four ( Figure 1). Along with the Λ-Transform, Λ-Fractional Derivative and Λ-Space come along. So, the non-local Λ-Fractional derivative in the initial space becomes local in the dual Λ-S where the differential may be defined, and differential geometry may be generated. Solving the problem in the dual Λ-Space, using conventional methods, the results may be transferred to the initial space. Along with the Λ-Transform, Λ-Fractional Derivative and Λ-Space come along. So, the non-local Λ-Fractional derivative in the initial space becomes local in the dual Λ-S where the differential may be defined, and differential geometry may be generated. Solving the problem in the dual Λ-Space, using conventional methods, the results may be transferred to the initial space.
Initially, the Λ-Transform dictates that a function f (x) R the initial space (x, f (x)) is transformed in Λ-Space in two steps: First, we have Λ 1 {f (x)}, which transforms f (x) from the initial to the intermediate space (x, F(x)) (F(x) R) as thus: ). Next, we complete the transformation with the final step Λ 2 {F(x)} by transforming variable x in the initial space to variable X R in Λ-Space: where and, It is evident from the procedure that the Rieman-Liouville integral plays a central role in Λ-Transform since it is the heart of this transform.
As far as the derivatives are concerned, we assume that the Λ-T transforms Λ-FD of any function f (x) of any order to itself in the intermediate space to the corresponding Derivative in Λ-Space. Therefore, we have: Later, we solve the Equation in Λ-Space (as an ODE) and afterward, with the inverse Λ-Transforms Λ 1 −1 and Λ 2 −1 , we restore the solution F(X) of the Equation in Λ-Space to f (x) in the initial space. In our case: And Further information may be found in Reference Lazopoulos [31] and References Samko [22], Podlubny [32], Oldham [33]. In the next paragraphs, viscoelasticity is studied using Λ-FD, where the pre-mentioned method will be applied and clarified explicitly and in every detail in the following parts of this article.

Models of Viscoelasticity-The Conventional Zener Model
There are many models of viscoelastic phenomena. The most representative is the viscoelastic model proposed by Zener. This model consists of the elastic springs with elastic moduli E 1 and E 2 and the viscous dashpot having a viscosity constant C. Then, the standard linear solid model called Zener viscoelastic model is shown in Figure 2.  (18) In that case, the creep compliance J(t), (constant stress σ), and relaxation G(t) (constant strain ε) functions take the form (Mainardi [17]): where Applying Equation (19) for the viscoelastic materials with E1 = 0.10 × 10 9 N/m 2 , E2 = 2.0 × 10 9 N/m 2 and C = 10 × 10 9 Ns/m 2 , and the initial conditions J(0) = 1/E2,G(0) = E2, the function of the compliance concerning time is shown in Figure 3   The equation that describes the viscoelastic deformation of the 3-parameter Zener model is: In that case, the creep compliance J(t), (constant stress σ), and relaxation G(t) (constant strain ε) functions take the form (Mainardi [17]): where Applying Equation (19) for the viscoelastic materials with E 1 = 0.10 × 10 9 N/m 2 , E 2 = 2.0 × 10 9 N/m 2 and C = 10 × 10 9 Ns/m 2 , and the initial conditions J(0) = 1/E 2 , G(0) = E 2 , the function of the compliance concerning time is shown in Figure 3 Axioms 2021, 10, x FOR PEER REVIEW 6 of 17 In that case, the creep compliance J(t), (constant stress σ), and relaxation G(t) (constant strain ε) functions take the form (Mainardi [17]): where Applying Equation (19) for the viscoelastic materials with E1 = 0.10 × 10 9 N/m 2 , E2 = 2.0 × 10 9 N/m 2 and C = 10 × 10 9 Ns/m 2 , and the initial conditions J(0) = 1/E2,G(0) = E2, the function of the compliance concerning time is shown in Figure 3  Moreover, the function of the relaxation modulus G(t) is shown in Figure 4 (Equation (20)):

The Λ-Fractional Zener Viscoelastic Model
Since the non-local action of the friction has been experimentally proven in viscoelastic models, Zener's viscoelastic model has been studied in the context of fractional analysis (Magin [13], Lazopoulos [16], Mainardi [17], Atanackovic [14]). We will first study the case of creep, and afterward, relaxation.

Creep
In this case, stress σ is constant (σ = σ0). Therefore, in the initial space (t, f(t)), the only acceptable equation is the Ordinary Differential Equation (ODE) for creep, which is extracted from the general viscoelasticity Equation, Equation (18). (C is the viscosity constant, and E1, E2 are the elastic constraints of the model): This is generally considered the according FDE of the creep ODE for 0 < γ ≤ 1. Up to this point, the procedure is the same as the common practice that most scientists follow. Nevertheless, this time Equation (23) is not solved directly; it is transformed into an ODE in Λ-Space.

The Λ-Fractional Zener Viscoelastic Model
Since the non-local action of the friction has been experimentally proven in viscoelastic models, Zener's viscoelastic model has been studied in the context of fractional analysis (Magin [13], Lazopoulos [16], Mainardi [17], Atanackovic [14]). We will first study the case of creep, and afterward, relaxation.

Creep
In this case, stress σ is constant (σ = σ 0 ). Therefore, in the initial space (t, f (t)), the only acceptable equation is the Ordinary Differential Equation (ODE) for creep, which is extracted from the general viscoelasticity Equation, Equation (18). (C is the viscosity constant, and E 1 , E 2 are the elastic constraints of the model): or where m = E 2 E 1 E 1 + E 2 , and b = CE 2 E 1 + E 2

Setup of Fractional Differential Equation in the Initial Space
To find the fractional version of Equation (22) This is generally considered the according FDE of the creep ODE for 0 < γ ≤ 1. Up to this point, the procedure is the same as the common practice that most scientists follow. Nevertheless, this time Equation (23) is not solved directly; it is transformed into an ODE in Λ-Space.

Λ 1 -Transform of Fractional Differential Equation from the Initial to the Intermediate Space
To transform it to the intermediate space, we apply the Λ 1 -Transform to Equation (23): And the equation that arises is: To execute the transform Λ 1 in Equation (25), we use the transform described in Equation (10): Especially when the Λ-fractional derivative of a function f(t) is transformed, then we assume that the result is the FD itself: And due to Equations (14) and (23), and considering a = 0: Consequently, the terms in Equation (25) are calculated as: Therefore, Equation (25) becomes: It is apparent that the last Equation (Equation (33)), although it reminds us Equation (22), is not similar due to the term t 1−γ Γ(2−γ) on the LHS. This last equation lies on the intermediate space (t, F(t)) (t > 0, F(t) R) and connects the FDE in the initial space with the ODE in Λ-Space.

Λ2-Transform of Fractional Differential Equation from the Intermediate to the Λ-Space
Through the Transform: and the inverse relation: t = (T · Γ(3 − γ)) 1 2−γ (35) We can then use Equation (31) and Λ 2 -Transform on Equation (33) to have: Which becomes: Or, where, Therefore, the initial right ODE for the specific problem (Equation (22)) in the initial space  Figure 2): And, Thus, the solution to the problem, in the initial space, is extracted. In Figure 5, the solving procedure is depicted in four steps. This process follows the example of a Laplace Transform. Only that two steps are required in the latter case, while in our case, the steps are 4. Or, where, Therefore, the initial right ODE for the specific problem (Equation (22) (23). It is directly connected with Equation (23), but it is slightly different in such a way that it can be transformed in Λ-Space.

) in the initial space (t, f(t)) is transformed to the according Fractional Differential Equation in the initial space Equation (23), and the intermediate equation in the intermediate space (t, F(t)) (Equation (33)), and finally to the according ODE (Equation (38)) in Λ-Space (T, F(T)). Equation (38) is a transformed version of Equation
Thus, the solution to the problem, in the initial space, is extracted. In Figure 5, the solving procedure is depicted in four steps. This process follows the example of a Laplace Transform. Only that two steps are required in the latter case, while in our case, the steps are 4.

Setup of Fractional Differential Equation in the Initial Space
To find the fractional version of Equation (42) This is generally considered the according FDE of the relaxation ODE for 0 < γ ≤ 1. As with creep, up to this point, the procedure is the same as the common practice that most scientists follow. Nevertheless, this time Equation (44) is not solved directly; it is transformed into an ODE in Λ-Space.

Λ 1 -Transform of Fractional Differential Equation from the Initial to the Intermediate Space
At this point, we must remind the two crucial transforms in the initial space (Equations (26) and (27)) Λ 1 -Transform of a function f(t): : We apply the Λ 1 -Transform to Equation (44): Therefore, the equation that arises is: As in creep, Λ 1 -Transform is given by the following premises: Therefore, the terms in Equation (48) become: Therefore, Equation (48) becomes: In the initial space (t, f (t)): In the intermediate space (t, F(t)): Again, Equation (55) or, Equation (57)  Finally, as in the case of creep, using the transform: And, The solution to the problem in the initial space is extracted. In Figure 6, the solving procedure is depicted in four steps. As in creep, it is also obvious that this process follows the example of a Laplace Transform. Only that two steps are required in the latter case, while in our case, the steps are four. It should be clarified at this point, just as with the case of Laplace Transform, that the transformed equations usually have no physical meaning. Λ-Transform, as well as Laplace Transform is a mathematical trick to accommodate the solution of the initial equation.

Application of Λ-Derivative to the Zener Model
Let us assume the Zener model depicted in Figure 2, as described in Section 3. As already mentioned, the values of the characteristic constants are: E1 = 0.10 × 10 9 N/m 2 , E2 = 2.0 × 10 9 N/m 2 and C = 10 × 10 9 Ns/m 2 , and the initial conditions J (0) = 1/E2, G(0) = E2. The resulting diagrams for creep compliances and relaxation functions for various γ are shown in the following sections.

Creep
The problem of creep was solved entirely in Section 4.1. The resulting ODE in Λ-Space was solved arithmetically (Equation (38)), as well as the inverse transform (Equations (39) and (40)). For the solution of the resulting ODE in Λ-space, we use a fifth-order Runge-Kutta method. Initial condition (T = 0) for the case of creep was taken as JΛ(0) = 1/E2, and for the case of relaxation, was taken as GΛ(0) = Ε2. for the numerical evaluation of this derivative, based on the so-called "L1-Method", as described in the book of Li et al. [35]. The occurring diagrams are: From the above figures, we may conclude that in creep, there is a small difference among the curves concerning γ close to 1, as expected, while that difference is more evident for γ much smaller than 1. (Additionaly, in Figure 7 we see the variation of the creep compliance versus time in Λ-Space.) The parameter γ expresses the extent that the microstructure of the material influences the viscoelastic response. Small γ designates a more porous material, which means a softer viscoelastic response. Therefore, it seems that the physical track of the phenomenon is utterly self-evident. More porous materials-softer response. Nevertheless, these results are not always compatible with other results of similar studies. For instance, the results on page 64 of Mainardi [17] are quite different. The point here is that in our case, the results seem more physically justified. Lastly, the units of measurement are unchanged in the initial space since the Λ-Transform (Figures 5 and  6) is structured in such a way that these units remain intact. Therefore, as it can be most clearly seen in Figures 8 and 9 we have Pa for stress and s for time.

Application of Λ-Derivative to the Zener Model
Let us assume the Zener model depicted in Figure 2, as described in Section 3. As already mentioned, the values of the characteristic constants are: E 1 = 0.10 × 10 9 N/m 2 , E 2 = 2.0 × 10 9 N/m 2 and C = 10 × 10 9 Ns/m 2 , and the initial conditions J (0) = 1/E 2 , G(0) = E 2 . The resulting diagrams for creep compliances and relaxation functions for various γ are shown in the following sections.

Creep
The problem of creep was solved entirely in Section 4.1. The resulting ODE in Λ-Space was solved arithmetically (Equation (38)), as well as the inverse transform (Equations (39) and (40)). For the solution of the resulting ODE in Λ-space, we use a fifthorder Runge-Kutta method. Initial condition (T = 0) for the case of creep was taken as J Λ (0) = 1/E 2 , and for the case of relaxation, was taken as G Λ (0) = E 2 .
For the transformation of the values of time and solution from Λ-space to real space via the Riemann-Liouville fractional derivative ( RL 0 D 1−γ t ), we constructed an algorithm for the numerical evaluation of this derivative, based on the so-called "L1-Method", as described in the book of Li et al. [35].
The occurring diagrams are: From the above figures, we may conclude that in creep, there is a small difference among the curves concerning γ close to 1, as expected, while that difference is more evident for γ much smaller than 1. (Additionaly, in Figure 7 we see the variation of the creep compliance versus time in Λ-Space.) The parameter γ expresses the extent that the microstructure of the material influences the viscoelastic response. Small γ designates a more porous material, which means a softer viscoelastic response. Therefore, it seems that the physical track of the phenomenon is utterly self-evident. More porous materials-softer response. Nevertheless, these results are not always compatible with other results of similar studies. For instance, the results on page 64 of Mainardi [17] are quite different. The point here is that in our case, the results seem more physically justified. Lastly, the units of measurement are unchanged in the initial space since the Λ-Transform (Figures 5 and 6) is structured in such a way that these units remain intact. Therefore, as it can be most clearly seen in Figures 8 and 9 we have Pa for stress and s for time.

Relaxation
The problem of relaxation was solved entirely in Section 4.2. The resulting ODE in Λ-Space was solved arithmetically (Equation (57)), as well as the inverse transform (Equa- Figure 9. Variation of the creep compliance versus time for γ = 1, 0.8, 0.7, 0.6 in the initial space.

Relaxation
The problem of relaxation was solved entirely in Section 4.2. The resulting ODE in Λ-Space was solved arithmetically (Equation (57)), as well as the inverse transform (Equations (59) and (60)). The occurring diagrams are: In relaxation, as in creep, the main attribute of the fractional phenomenon is still present. As γ decreases, the material's response becomes softer, signifying that friction and energy losses increase as γ decreases. Again, the corresponding results on page 64 of Mainardi [17] are quite different. The point here is that in our case, the results seem more physically justified. In Figure 10, we see the variation of the relaxation modulus versus time in Λ-Space. Moreover, from Figures 11 and 12 a problem emerges. Although we have an initial-value problem, nevertheless it does not comply with the initial conditions in Figures 11 and 12. This phenomenon is due to a problem in the structure of the Riemann-Liouville derivative, which has a singularity of the function at 0. So, for values of t which are very close to zero the value of RL fractional derivative exceeds the value of the function at 0 and, thus, we get this inconsistency. Lastly, as in creep, the units of measurement are again unchanged in the initial space since the Λ-Transform (Figures 5 and 6) is structured in such a way that these units remain intact. Therefore, as it can be most clearly seen in Figures 11 and 12, we have Pa for stress and s for time. Figure 9. Variation of the creep compliance versus time for γ = 1, 0.8, 0.7, 0.6 in the initial space.

Relaxation
The problem of relaxation was solved entirely in Section 4.2. The resulting ODE in Λ-Space was solved arithmetically (Equation (57)), as well as the inverse transform (Equations (59) and (60)). The occurring diagrams are: In relaxation, as in creep, the main attribute of the fractional phenomenon is still present. As γ decreases, the material's response becomes softer, signifying that friction and energy losses increase as γ decreases. Again, the corresponding results on page 64 of Mainardi [17] are quite different. The point here is that in our case, the results seem more physically justified. In Figure 10, we see the variation of the relaxation modulus versus time in Λ-Space. Moreover, from Figures 11 and 12 a problem emerges. Although we have an initial-value problem, nevertheless it does not comply with the initial conditions in Figures 11 and 12. This phenomenon is due to a problem in the structure of the Riemann-Liouville derivative, which has a singularity of the function at 0. So, for values of t which are very close to zero the value of RL fractional derivative exceeds the value of the function at 0 and, thus, we get this inconsistency. Lastly, as in creep, the units of measurement are again unchanged in the initial space since the Λ-Transform (Figures 5 and 6) is structured in such a way that these units remain intact. Therefore, as it can be most clearly seen in Figures 11 and 12, we have Pa for stress and s for time.

Conclusions
Λ-Fractional Derivative, a new innovative fractional derivative, was used to study viscoelasticity. The main advantage is that it makes the whole procedure mathematically proper, a feature that other fractional approaches lack. Thus, Λ-FD seems to be the only authentic non-local derivative. In contrast, the other so-called fractional "derivatives" are just operators since they cannot satisfy the differential topology requirements for a proper derivative. Hence, that derivative opens new horizons in fractional calculus, and generally non-local analysis since differentials and differential geometry could be defined. Finally, the Λ-Transform of a "Fractional Differential Equation" to an Ordinary Differential Equation in the dual Λ-Space simplifies the analysis since an ODE solution is much easier than a Fractional Differential Equation in the initial space. Nevertheless, it must be stressed out that the solution of any FDE in the initial space is mathematically questionable. In the present article, the fractional Zener viscoelastic model was studied. It seems that as the fractional-order γ decreases, the viscoelastic response becomes softer in creep and relaxation due to friction and energy losses. That means that as a porous microstructure takes over, the material exhibits a softer response. The results presented may express deep and

Conclusions
Λ-Fractional Derivative, a new innovative fractional derivative, was used to study viscoelasticity. The main advantage is that it makes the whole procedure mathematically proper, a feature that other fractional approaches lack. Thus, Λ-FD seems to be the only authentic non-local derivative. In contrast, the other so-called fractional "derivatives" are just operators since they cannot satisfy the differential topology requirements for a proper derivative. Hence, that derivative opens new horizons in fractional calculus, and generally non-local analysis since differentials and differential geometry could be defined. Finally, the Λ-Transform of a "Fractional Differential Equation" to an Ordinary Differential Equation in the dual Λ-Space simplifies the analysis since an ODE solution is much easier than a Fractional Differential Equation in the initial space. Nevertheless, it must be stressed out that the solution of any FDE in the initial space is mathematically questionable. In the present article, the fractional Zener viscoelastic model was studied. It seems that as the fractional-order γ decreases, the viscoelastic response becomes softer in creep and relaxation due to friction and energy losses. That means that as a porous microstructure takes over, the material exhibits a softer response. The results presented may express deep and profound natural intuition, and even though they do not always agree with other studies, they make the perspective of the Λ-FD method even more intriguing.