Research Papers

Wavelets Galerkin Method for the Fractional Subdiffusion Equation

[+] Author and Article Information
M. H. Heydari

Faculty of Mathematics,
Yazd University,
P. O. Box 89195-741,
Yazd, Iran;
The Laboratory of Quantum
Information Processing,
Yazd University,
P. O. Box 89195-741,
Yazd, Iran
e-mail: heydari@stu.yazd.ac.ir

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received April 28, 2015; final manuscript received July 14, 2016; published online August 22, 2016. Assoc. Editor: Hiroshi Yabuno.

J. Comput. Nonlinear Dynam 11(6), 061014 (Aug 22, 2016) (7 pages) Paper No: CND-15-1114; doi: 10.1115/1.4034391 History: Received April 28, 2015; Revised July 14, 2016

The time fractional subdiffusion equation (FSDE) as a class of anomalous diffusive systems has obtained by replacing the time derivative in ordinary diffusion by a fractional derivative of order 0<α<1. Since analytically solving this problem is often impossible, proposing numerical methods for its solution has practical importance. In this paper, an efficient and accurate Galerkin method based on the Legendre wavelets (LWs) is proposed for solving this equation. The time fractional derivatives are described in the Riemann–Liouville sense. To do this, we first transform the original subdiffusion problem into an equivalent problem with fractional derivatives in the Caputo sense. The LWs and their fractional operational matrix (FOM) of integration together with the Galerkin method are used to transform the problem under consideration into the corresponding linear system of algebraic equations, which can be simply solved to achieve the solution of the problem. The proposed method is very convenient for solving such problems, since the initial and boundary conditions are taken into account, automatically. Furthermore, the efficiency of the proposed method is shown for some concrete examples. The results reveal that the proposed method is very accurate and efficient.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

The graphs of the approximate solution (left side) and absolute error (right side) for Example 1 when α=0.25 with m̂=48

Grahic Jump Location
Fig. 2

The graphs of the approximate solution (left side) and absolute error (right side) for Example 2 when α=0.50 with m̂=60

Grahic Jump Location
Fig. 3

The graphs of the approximate solution (left side) and absolute error (right side) for Example 3 when α=0.75 with m̂=48




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In