Abstract

Bagley-Torvik equation appears in viscoelasticity problems where fractional derivatives seem to play an important role concerning empirical data. There are several works treating this equation by using numerical methods and analytic formulations. However, the analytical solutions presented in the literature consider particular cases of boundary and initial conditions, with inhomogeneous term often expressed in polynomial form. Here, by using Laplace transform methodology, the general inhomogeneous case is solved without restrictions in boundary and initial conditions. The generalized Mittag-Leffler functions with three parameters are used and the solutions presented are expressed in terms of Wiman’s functions and their derivatives.

1. Introduction

Bagley-Torvik equations (BTE) firstly appeared in their seminal work [1] where they proposed to model viscoelastic behavior of geological strata, metals, and glasses by using fractional differential equations, showing that this approach is effective in describing structures containing elastic and viscoelastic components.

Initially, in [2], inhomogeneous BTE was studied with an analytical solution being proposed. Since then, there were several works to solve BTE, starting with numerical procedures for a reformulated BTE as a system of functional differential equations of order 1/2 [35].

Following a numerical way for solving BTE, a generalization of Taylor’s and Bessel’s collocation method [6, 7] and the use of evolutionary computation [8] provided acceptable solutions from engineering point of view.

Approximation techniques were successfully applied to BTE mainly by using enhanced homotopic perturbation methods [9, 10], fractional iteration techniques [11], and cubic polynomial spline functions [12].

Analytical exact solutions for BTE were obtained in [13] for the particular initial condition , considering the boundary condition given by . Besides, by using a modified generalized Laguerre spectral method for fractional differential equations, BTE was solved in [14] for some specific conditions.

The aim of this work is to provide an analytical solution for the most generic case of inhomogeneous BTE in terms of the derivatives of Wiman’s functions [2] considering that three-parameter generalized Mittag-Leffler function, introduced in [15], can be better explored.

The approach presented here consists of using the Laplace transform of Prabhakar’s function, in order to solve the most general BTE, carrying out a subsequent transformation of the achieved solution in terms of the Mittag-Leffler functions.

In the next section, the BTE solution by using the referred to method is presented, followed by a section with the solution expressed by using Wiman’s functions and their derivatives, motivated by the fact that Wiman’s functions can be implemented in computational software like , easing the solution of BTE in practical applications.

2. Laplace Transformation Solution of BTE: Mittag-Leffler Functions

The most general inhomogeneous BTE [1] is given bywhich, by using operational notation, can be expressed as

Here the domain of is considered to be the time; consequently it is assumed that . The Riemann-Liouville fractional integral is adopted and, therefore, Caputo’s fractional derivative can be simply denoted by [16, 17].

The method to be used to solve BTE is the traditional Laplace transform that, due to the fractional nature of the exponents, results from Prabhakar functions, that is, three-parameter Mittag-Leffler generalized functions [2, 15].

Considering Caputo’s fractional derivative [18], the Laplace transform of the fractional exponent term isimplying that the Laplace transform of the BTE is written as

Defining and and considering the binomial development,the Laplace transform of BTE’s solution is

In order to invert the Laplace transform to obtain the general BTE’s solution, is rewritten expanding the sums:

As shown in [2, 15], it is possible to define three-parameter Mittag-Leffler general functions aswith being the Pochhammer symbol. When , these functions are called Wiman’s functions.

Here, the main fact to be used is related to the inverse of the Laplace transform of the main terms that appear in (7), given by

Consequently, the first four terms of (7) can be inverted and written as

The inhomogeneous term is calculated by using convolution:and, consequently, the general solution of BTE is given by

In order to simplify the BTE general solution expressed by (12), a factor can be defined:with . Therefore the BTE solution is expressed as

Operator represents the integer order classical derivative, that is, , with the property .

Expression (14) is a general analytical solution of BTE in terms of three-parameter Mittag-Leffler functions, . However, the solution established in [2] for the special case with is given by using Wiman’s functions derivatives, that is, two-parameter Mittag-Leffler functions.

As Wiman’s functions and derivatives are part of several software packages, as , for instance, and three-parameter Mittag-Leffler are not, the results will be modified to be compatible with this fact, easing the practical approach of the problem.

3. BTE Solution via Wiman’s Functions and Derivatives

Here, the solution of BTE expressed by (14) is modified, in order to be written in terms of Wiman’s functions and their derivatives. The explicit forms of these functions are given by

For the sake of clearness, the work will be divided into two parts: homogeneous solution terms and inhomogeneous solution term.

3.1. Homogeneous Solution Terms

Taking the homogeneous part of (14), denoted by ,

If the and sums are uniformly convergent,

Observing that the last sum can be written as the th derivative of Wiman’s function, the homogeneous part of the solution is

3.2. Inhomogeneous Solution Term

Calling the inhomogeneous term of the BTE solution , its expression is

Considering the convergence of the integral and uniform convergence of the sums, their order can be changed and it can be written as

It can be noticed that the last sum of (20) can be written as the -order derivative of Wiman’s function and, consequently,

Then, definingit follows thataccording to the BTE analytical solution obtained in [2] for .

3.3. Complete BTE Solution

Considering that , the general analytical solution of BTE, in terms of Wiman’s function, iswith and .

4. Conclusion

A general analytical solution of BTE defined by Caputo’s fractional derivatives was obtained in terms of Wiman’s functions and their derivatives. The reasoning was conducted by using Laplace transform and Mittag-Leffler functions.

The final expression permits the calculations in practical cases as it is built considering Wiman’s functions that are part of the most usual numerical packages.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.