An efficient fourth order Hermite spline collocation method for time fractional diffusion equation describing anomalous diffusion in two space variables
- Resource Type
- Original Paper
- Authors
- Priyanka; Sahani, Saroj; Arora, Shelly
- Source
- Computational and Applied Mathematics. 43(4)
- Subject
- Diffusion equation
Caputo derivative
L1-scheme
Sobolev spaces
Orthogonal collocation
Bi-quintic Hermite splines
H~1m-norm
35C11
35G31
35K05
65D07
65M06
65M12
- Language
- English
- ISSN
- 2238-3603
1807-0302
Anomalous diffusion of particles in fluids is better described by the fractional diffusion models. A robust hybrid numerical algorithm for a two-dimensional time fractional diffusion equation with the source term is presented. The well-known L1 scheme is considered for semi-discretization of the diffusion equation. To interpolate the semi-discretized equation, orthogonal collocation with bi-quintic Hermite splines as the basis is chosen for the smooth solution. Quintic Hermite splines interpolate the solution as well as its first and second order derivatives. The technique reduces the proposed problem to an algebraic system of equations. Stability analysis of the implicit scheme is studied using H~1mO(h4)O(Δt)2-αΔtα-norm defined in Sobolev space. The optimal order of convergence is found to be of order H~1mO(h4)O(Δt)2-αΔtα in spatial direction and is of order H~1mO(h4)O(Δt)2-αΔtα in the temporal direction where h is the step size in space direction and H~1mO(h4)O(Δt)2-αΔtα is the step size in time direction and H~1mO(h4)O(Δt)2-αΔtα is the fractional order of the derivative. Numerical illustrations have been presented to discuss the applicability of the proposed hybrid numerical technique to the problems having fractional order derivative.