Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2015, Article ID 591715, 4 pages
Research Article

Analytical Solution of General Bagley-Torvik Equation

Escola Politécnica da Universidade de São Paulo, Avenida Prof. Luciano Gualberto, Travessa 3, No. 158, 05508-900 São Paulo, SP, Brazil

Received 21 October 2015; Accepted 18 November 2015

Academic Editor: Ivan D. Rukhlenko

Copyright © 2015 William Labecca et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


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.


  1. P. J. Torvik and R. L. Bagley, “On the appearance of the fractional derivative in the behavior of real materials,” Journal of Applied Mechanics, vol. 51, no. 2, pp. 294–298, 1984. View at Publisher · View at Google Scholar · View at Scopus
  2. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999. View at MathSciNet
  3. A. H. Bhrawy, T. M. Taha, and J. A. T. Machado, “A review of operational matrices and spectral techniques for fractional calculus,” Nonlinear Dynamics, vol. 81, no. 3, pp. 1023–1052, 2015. View at Publisher · View at Google Scholar · View at MathSciNet
  4. K. Diethelm and N. J. Ford, “Numerical solution of the Bagley-Torvik equation,” BIT. Numerical Mathematics, vol. 42, no. 3, pp. 490–507, 2002. View at Google Scholar · View at MathSciNet
  5. Z. H. Wang and X. Wang, “General solution of the Bagley-Torvik equation with fractional-order derivative,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 5, pp. 1279–1285, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  6. Y. Çenesiz, Y. Keskin, and A. Kurnaz, “The solution of the Bagley-Torvik equation with the generalized TAYlor collocation method,” IEEE Engineering and Applied Mathematics, vol. 347, no. 2, pp. 452–466, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  7. Ş. Yüzbaşı, “Numerical solution of the Bagley-Torvik equation by the Bessel collocation method,” Mathematical Methods in the Applied Sciences, vol. 36, no. 3, pp. 300–312, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  8. M. A. Raja, J. A. Khan, and I. M. Qureshi, “Solution of fractional order system of Bagley-Torvik equation using evolutionary computational intelligence,” Mathematical Problems in Engineering, vol. 2011, Article ID 675075, 18 pages, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. F. Abidi and K. Omrani, “The homotopy analysis method for solving the Fornberg-Whitham equation and comparison with Adomian's decomposition method,” Computers & Mathematics with Applications, vol. 59, no. 8, pp. 2743–2750, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. M. Zolfaghari, R. Ghaderi, A. SheikholEslami et al., “Application of the enhanced homotopy perturbation method to solve the fractional-order Bagley–Torvik differential equation,” Physica Scripta, vol. 2009, article T136, Article ID 014032, 2009. View at Publisher · View at Google Scholar
  11. T. Mekkaoui and Z. Hammouch, “Approximate analytical solutions to the Bagley-Torvik equation by the fractional iteration method,” Annals of the University of Craiova-Mathematics and Computer Science Series, vol. 39, no. 2, pp. 251–256, 2012. View at Google Scholar · View at MathSciNet
  12. W. K. Zahra and S. M. Elkholy, “Cubic spline solution of fractional Bagley-Torvik equation,” Electronic Journal of Mathematical Analysis and Applications, vol. 1, no. 2, pp. 230–241, 2013. View at Google Scholar
  13. D. Baleanu, A. H. Bhrawy, and T. M. Taha, “A modified generalized Laguerre spectral method for fractional differential equations on the half line,” Abstract and Applied Analysis, vol. 2013, Article ID 413529, 12 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  14. S. Staněk, “Two-point boundary value problems for the generalized Bagley-Torvik fractional differential equation,” Central European Journal of Mathematics, vol. 11, no. 3, pp. 574–593, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  15. T. R. Prabhakar, “A singular integral equation with a generalized Mittag-Leffler function in the kernel,” Yokohama Mathematical Journal, vol. 19, pp. 7–15, 1971. View at Google Scholar
  16. K. M. Kolwankar and A. D. Gangal, “Fractional differentiability of nowhere differentiable functions and dimensions,” Chaos, vol. 6, no. 4, pp. 505–513, 1996. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  17. T. Kisela, Fractional differential equations and their applications [M.S. thesis], University of l’Aquila & Brno University of Technology, Brno, Czech Republic, 2008.
  18. A. M. O. Anwar, F. Jarad, D. Baleanu, and F. Ayaz, “Fractional Caputo heat equation within the double Laplace transform,” Romanian Journal of Physics, vol. 58, no. 1-2, pp. 15–22, 2013. View at Google Scholar · View at MathSciNet