Research Article  Open Access
Constantin Bota, Bogdan Căruntu, Mădălina Sofia Paşca, Dumitru Ţucu, Marioara Lăpădat, "Least Squares Differential Quadrature Method for the Generalized Bagley–Torvik Fractional Differential Equation", Mathematical Problems in Engineering, vol. 2020, Article ID 4806387, 7 pages, 2020. https://doi.org/10.1155/2020/4806387
Least Squares Differential Quadrature Method for the Generalized Bagley–Torvik Fractional Differential Equation
Abstract
In this paper, the least squares differential quadrature method for computing approximate analytical solutions for the generalized Bagley–Torvik fractional differential equation is presented. This new method is introduced as a straightforward and accurate method, fact proved by the examples included, containing a comparison with previous results obtained by using other methods.
1. Introduction
Even though the fractional calculus is a relatively old concept, dating back from the time of the mathematicians L’Hospital and Leibniz in the 17^{th} century, it has seen a significant development in the last decades only. The fractional differential equations have numerous applications in biology, physics, chemistry, and engineering, which is why many scientists are concerned about finding more effective ways to solve them. Unfortunately, finding analytical solutions for fractional differential equations is often difficult or even impossible, which is why in recent years various numerical and approximate methods have been developed. One of the fractional differential equations with wide applicability in engineering is the Bagley–Torvik equation. This equation was introduced by Bagley and Torvik in [1] in 1984 to “model a motion of a rigid plate immersed in a Newtonian fluid” and in [2] for modeling the damping properties of the polymers. In engineering (construction, biotechnology, and chemistry), condensing polymers are widely used, and in their study, the Bagley–Torvik equation is the one employed [3, 4].
Bagley–Torviktype fractional differential equations have been studied both numerically and analytically in numerous articles. Among the methods used to solve this equation, we mention the following: numerical methods [5], Adomian decomposition method [6], discrete spline methods [7], Haar wavelet method [8, 9], homotopy perturbation method [10], sinc collocation method [11, 12], cubic spline method [13], quadratic spline solution [14], Bspline collocation method [15, 16], Chebyshev collocation method [17], hybrid functions approximation [18], harmonic wavelets [19], predictorcorrector method of Adams type [20], spectral methods [21], fractional natural decomposition method [22], and finite element method [23].
The motivation of this paper is to introduce a new method for obtaining analytical approximate solutions for fractional differential equations. This method is obtained as a combination of the differential quadrature method (DQM) ([24]) and the least squares method (LSM), and we named it the least squares differential quadrature method (LSDQM). We applied this new method to find solutions of the following generalized Bagley–Torvik fractional differential equation [25]:where , , together with the conditions:where , , , and are constants (for , we obtain boundary conditions, while for , we obtain initial conditions). , and are given such that problems (1) and (2) satisfy the existence and uniqueness conditions for a continuous solution [25].
denotes Caputo’s fractional derivative:
In the next section, we will present the least squares differential quadrature method (LSDQM), and in the third section, some numerical examples are presented followed by the conclusions.
2. The Least Squares Differential Quadrature Method
First, we will consider a numerical mesh of the interval by means of a partition consisting of points: . In order to find an approximate solution of the generalized Bagley–Torvik problems (1) and (2), we will compute the values of a certain functional (the functional (12) introduced in the following) at the point .
To equation (1), we attach the following operator:
We denote by an approximate solution of equation (1). By replacing in the exact solution with this approximate solution, we obtain the following reminder:
Definition 1. We call an approximate solution of problems (1) and (2) related to the partition an approximate polynomial solution, which satisfies the following relations:
Definition 2. We consider the sequence of polynomials:We call the sequence of polynomials convergent to the solution of problems (1) and (2) ifTaking into account the above definitions, we will compute approximate polynomial solutions for a Bagley–Torvik problem of the types (1) and (2) by taking the steps described in the following algorithm:(i)Step 1. we will compute an approximate solution of the type where the sequence of polynomials must also satisfy the boundary conditions: and the constants are calculated in the following steps.(ii)Step 2. From the boundary conditions, we obtain and as functions of . Thus, the expression of from now on will be a function of and only.(iii)Step 3. We attach to problems (1) and (2) the functional :(iv)Step 4. By minimizing the functional (12), we obtain the values of the coefficients which give the minimum of (12).(v)Step 5. After computing , we return to the boundary conditions and obtain the final value of and .(vi)Step 6. Finally, we replace the values , calculate in the expression of , and denote by the analytical approximate polynomial solution by LSDQM of problems (1) and (2).The following convergence theorem is satisfied:
Theorem 1. The sequence of polynomials satisfies the following relation:
Proof. Let be an exact solution of problems (1) and (2), which means from hypothesis that there exists a sequence of polynomials converging to (according to the Weierstrass theorem on Polynomial Approximation, see [24]), i.e., . Taking into account Definition 2, this means that satisfiesDenoting by the values of the coefficients which give the minimum of the functional (12), it follows that . If is the approximate solution by LSDQM (with and computed by using the initial conditions), we obtainHence, .
We conclude that .
2.1. Error Estimation
We can rewrite equation (1) as
We denote with the following operator:
Thus, equation (1) becomes
Using (for simplicity) the notation: , in (4), we obtainwith as an approximate solution for problems (1) and (2), which means as satisfiestogether with the conditions:
We define the error function in the following way: , where is the exact solution for problems (1) and (2) and we obtain (using (18) and (19)) the differential equation for the error function:with the conditions:
The problem for the error function becomes
Solving (24) together with (23) in the same manner as described above for problems (1) and (2), we obtain the approximation . We will be able to determine the absolute maximum error:
In this manner, we can estimate the error without knowing the exact solution of the initial problems (1) and (2).
Remark 1. We remark that the above theorem proves the convergence of LSDQM since as the degree of the polynomial increases the remainder corresponding to the approximation tends to zero. This fact is illustrated in the second example.
3. Numerical Examples
In this section, we will present some examples of problems including generalized Bagley–Torvik fractional differential equations together with boundary or initial conditions and problems solved using the least squares differential quadrature method (LSDQM).
3.1. Example 1
We consider the generalized Bagley–Torvik equation [26, 27]:together with the initial conditions: where , , and . The exact solution of this problem is .
Approximate solutions for this problem with absolute errors larger than were proposed by ElMesiry et al. in [26] and with absolute errors larger than by Li and Zhao in [27].
We applied LSDQM to find an approximate analytical solution on the interval for (26).
In Step 1, we choose an approximate solution of the following type:
In Step 2, by using the initial conditions, we find and , and the approximate solution becomes
In Step 3, the corresponding reminder (5) is
We consider the simplest possible partition : , and we attach to the equation the functional (too long to be included here).
In Step 4, in order to find the minimum of this functional, we can compute the stationary points by equating to zero its first derivative (the computations are performed using the software Wolfram Mathematica). We obtain , and it is easy to show (by means of the second derivative) that this stationary point is indeed the minimum.
In the case of this problem, Step 5 is not needed since we already have the final values of and .
Finally, in Step 6, we replace , , and in the initial expression of , thus obtaining the exact solution of the equation: .
We remark that, in the case of this example, we could choose any partition since the particular form of the remainder (5) assures the fact that the stationary point is always no matter how large is.
3.2. Example 2
We consider the fractional Bagley–Torvik boundary value problem [21]:where , , and is the generalized hypergeometric function.
The exact solution of this problem is .
We remark the fact that if the exact solutions of the problem is a polynomial function (as in the case of the first example), then, naturally, in applying the first step of our method, we have chosen an approximate polynomial solution (10) of the same degree as the exact solution. Since, in this example, the exact solution is not a polynomial, we computed approximate polynomial solutions of several degrees, and for each degree, we calculated the maximal value of the absolute error on the whole interval for . The results presented in Table 1 together with a comparison with previous results from [21] clearly illustrate the convergence of the method.

The approximate polynomial solutions computed by LSDQM using an equidistant partition are
3.3. Example 3
Choosing in (1) , , , , , and , we obtain the problem [5, 28]:
This equation is actually a Riccatitype equation, and for , the exact solution of the problem is .
Following the LSDQM steps presented in the previous section, we computed approximate polynomial solutions of the 9th degree for problems (41) and (42) choosing in turn , , , and . The solutions are computed on the interval using again an equidistant partition . The results are presented in Figure 1.
For example, the approximate polynomial solutions corresponding to is
The maximal absolute error corresponding to this approximation is .
We remark that the results presented in Figure 1 are in good agreement with previously obtained results [5].
3.4. Example 4
While LSDQM, as most methods of this type, works best on small intervals, it can also be successfully employed in the case of large intervals. In order to illustrate this feature, we consider the Bagley–Torvik problem [6]:where .
The exact solution of this problem is not known.
An approximate solution obtained by using the Adomian decomposition method was proposed in [6] by V. DaftardarGejji and H. Jafari.
For , using an equidistant partition on the interval , we obtain the approximate analytical solution, presented in Figure 2:
We remark that our result is in very good agreement with the one in [6].
4. Conclusions
In this paper, the application of the least squares differential quadrature method (LSDQM) to the Bagley–Torvik fractional differential equation is presented.
Due to the fact that the method is relative straightforward, the approximations may be obtained in a quick and simple manner. The numerical examples included clearly illustrate the accuracy of the method by means of comparisons with solutions previously computed by other methods. If the exact solution of the problem is itself a polynomial, then usually LSDQM is able to find the exact solution, as in the case of the first example. If not, in most of the cases, the solutions provided by LSDQM are not only more precise, but being of the simplest form (i.e., polynomial solutions), they are much easier to use in subsequent computation than many of the solutions provided by other methods.
Data Availability
All the necessary data are included in the paper.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
References
 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 Site  Google Scholar
 R. L. Bagley and P. J. Torvik, “Fractional calculusa different approach to the analysis of viscoelastically damped structures,” AIAA Journal, vol. 21, no. 5, pp. 741–748, 1983. View at: Publisher Site  Google Scholar
 T. Aleroev, S. Erokhin, and E. Kekharsaeva, “Modeling of deformationstrength characteristics of polymer concrete using fractional calculus,” IOP Conf. Series: Materials Science and Engineering, vol. 365, Article ID 032004, 2018. View at: Publisher Site  Google Scholar
 B. H. E. Ibrahim, Q. Dong, and Z. Fan, “Existence for boundary value problems of twoterm Caputo fractional differential equations,” The Journal of Nonlinear Sciences and Applications, vol. 10, no. 2, pp. 511–520, 2017. View at: Publisher Site  Google Scholar
 H. Jafari, H. Tajadodi, H. Nazari, and C. M. Khalique, “Numerical solution of nonlinear Riccati differential equations with fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 11, pp. 179–182, 2010. View at: Publisher Site  Google Scholar
 V. DaftardarGejji and H. Jafari, “Adomian decomposition: a tool for solving a system of fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 301, no. 2, pp. 508–518, 2005. View at: Publisher Site  Google Scholar
 W. K. Zahra and M. Van Daele, “Discrete spline methods for solving two point fractional bagleytorvik equation,” Applied Mathematics and Computation, vol. 296, pp. 42–56, 2017. View at: Publisher Site  Google Scholar
 M. U. Rehman and R. A. Khan, “A numerical method for solving boundary value problems for fractional differential equations,” Applied Mathematical Modeling, vol. 36, no. 3, pp. 894–907, 2012. View at: Google Scholar
 S. Micula and C. Cattani, “On a numerical method based on wavelets for fredholmhammerstein integral equations of the second kind,” Mathematical Methods in the Applied Sciences, vol. 41, no. 18, pp. 9103–9115, 2018. View at: Publisher Site  Google Scholar
 Y.G. Wang, H.F. Song, and D. Li, “Solving twopoint boundary value problems using combined homotopy perturbation method and green’s function method,” Applied Mathematics and Computation, vol. 212, no. 2, pp. 366–376, 2009. View at: Publisher Site  Google Scholar
 S. Alkan, K. Yildirim, and A. Secer, “An efficient algorithm for solving fractional differential equations with boundary conditions,” Open Physics, vol. 14, no. 1, pp. 6–14, 2016. View at: Publisher Site  Google Scholar
 S. Yüzbasi, “Numerical solution of the BagleyTorvik equation by the bessel collocation method,” Mathematical Methods in the Applied Science, vol. 36, pp. 300–312, 2013. View at: Google Scholar
 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
 W. K. Zahra and S. M. Elkholy, “Quadratic spline solution for boundary value problem of fractional order,” Numerical Algorithms, vol. 59, no. 3, pp. 373–391, 2012. View at: Publisher Site  Google Scholar
 K. Thula and P. Roul, “A highorder Bspline collocation method for solving nonlinear singular boundary value problems arising in engineering and applied science,” Mediterranean Journal of Mathematics, vol. 15, 2018. View at: Google Scholar
 H. Allouche, A. Tazdayte, and K. Tigma, “Highly accurate method for solving singular boundaryvalue problems via padé approximation and twostep quartic Bspline collocation,” Mediterranean Journal of Mathematics, vol. 19, 2019. View at: Publisher Site  Google Scholar
 H. Jafari, A. Babaei, and S. Banihashemi, “A novel approach for solving an inverse reactiondiffusionconvection problem,” Journal of Optimization Theory and Applications, vol. 183, no. 2, pp. 688–704, 2019. View at: Publisher Site  Google Scholar
 S. Mashayekhi and M. Razzaghi, “Numerical solution of the fractional bagleytorvik equation by using hybrid functions approximation,” Mathematical Methods in the Applied Sciences, vol. 39, no. 3, pp. 353–365, 2016. View at: Publisher Site  Google Scholar
 C. Cattani, “A review on harmonic wavelets and their fractional extension,” Journal of Advanced Engineering and Computation, vol. 2, no. 4, pp. 224–238, 2018. View at: Publisher Site  Google Scholar
 K. Diethelm and J. Ford, “Numerical solution of the bagleytorvik equation,” BIT Numerical Mathematics, vol. 42, no. 3, pp. 490–507, 2002. View at: Google Scholar
 E. H. Doha, A. H. Bhrawy, and S. S. EzzEldien, “Efficient chebyshev spectral methods for solving multiterm fractional orders differential equations,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5662–5672, 2011. View at: Publisher Site  Google Scholar
 M. S. Rawashdeh and H. AlJammal, “Numerical solutions for systems of nonlinear fractional ordinary differential equations using the FNDM,” Mediterranean Journal of Mathematics, vol. 13, no. 6, pp. 4661–4677, 2016. View at: Publisher Site  Google Scholar
 M. Badr, A. Yazdani, and H. Jafari, “Stability of a finite volume element method for the timefractional advectiondiffusion equation,” Numerical Methods for Partial Differential Equations, vol. 34, no. 5, pp. 1459–1471, 2018. View at: Publisher Site  Google Scholar
 C. Shu, Differential Quadrature and its Application in Engineering, Springer, New York, NY, USA,, 2000.
 S. Stanek, “Twopoint boundary value problems for the generalized bagleytorvik fractional differential equation,” Central European Journal of Mathematics, vol. 11, no. 3, pp. 574–593, 2013. View at: Google Scholar
 A. E. M. ElMesiry, A. M. A. ElSayed, and H. A. A. ElSaka, “Numerical methods for multiterm fractional (arbitrary) orders differential equations,” Applied Mathematics and Computation, vol. 160, no. 3, pp. 683–699, 2005. View at: Publisher Site  Google Scholar
 Y. Li and W. Zhao, “Haar wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations,” Applied Mathematics and Computation, vol. 216, no. 8, pp. 2276–2285, 2010. View at: Publisher Site  Google Scholar
 B. Căruntu, C. Bota, M. S. Paşca, and M. Lăpădat, “Approximate solutions for Riccati differential equation of fractional order using the Least Squares differential quadrature method,” Scientific Bulletin of the ”Politehnica” University of Timisoara, vol. 64, no. 78, pp. 18–28, 2019. View at: Google Scholar
Copyright
Copyright © 2020 Constantin Bota 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.