A fast semi-discrete Kansa method to solve the two-dimensional spatiotemporal fractional diffusion equation

被引:26
作者
Sun, HongGuang [1 ]
Liu, Xiaoting [1 ]
Zhang, Yong [2 ,3 ]
Pang, Guofei [4 ]
Garrard, Rhiannon [2 ]
机构
[1] Hohai Univ, Coll Mech & Mat, State Key Lab Hydrol Water Resources & Hydraul En, Nanjing 210098, Jiangsu, Peoples R China
[2] Univ Alabama, Dept Geol Sci, Tuscaloosa, AL 35487 USA
[3] Hohai Univ, Coll Mech & Mat, Nanjing 210098, Jiangsu, Peoples R China
[4] Beijing Computat Sci Res Ctr, Beijing 100084, Peoples R China
基金
中国国家自然科学基金; 美国国家科学基金会;
关键词
Anomalous transport; Spatiotemporal FDE; Semi-discrete Kansa method; Vector fractional derivative; ADVECTION-DISPERSION EQUATIONS; FINITE-ELEMENT-METHOD; ANOMALOUS DIFFUSION; DIFFERENTIAL-EQUATIONS; HETEROGENEOUS MEDIA; SOLUTE TRANSPORT; BOUNDED DOMAINS; POROUS-MEDIA; TIME; APPROXIMATION;
D O I
10.1016/j.jcp.2017.05.012
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Fractional-order diffusion equations (FDEs) extend classical diffusion equations by quantifying anomalous diffusion frequently observed in heterogeneous media. Real-world diffusion can be multi-dimensional, requiring efficient numerical solvers that can handle long-term memory embedded in mass transport. To address this challenge, a semi-discrete Kansa method is developed to approximate the two-dimensional spatiotemporal FDE, where the Kansa approach first discretizes the FDE, then the Gauss-Jacobi quadrature rule solves the corresponding matrix, and finally the Mittag-Leffler function provides an analytical solution for the resultant time-fractional ordinary differential equation. Numerical experiments are then conducted to check how the accuracy and convergence rate of the numerical solution are affected by the distribution mode and number of spatial discretization nodes. Applications further show that the numerical method can efficiently solve two-dimensional spatiotemporal FDE models with either a continuous or discrete mixing measure. Hence this study provides an efficient and fast computational method for modeling super-diffusive, sub-diffusive, and mixed diffusive processes in large, two-dimensional domains with irregular shapes. (C) 2017 Elsevier Inc. All rights reserved.
引用
收藏
页码:74 / 90
页数:17
相关论文
共 67 条
[1]   Gaussian closure of one-dimensional unsaturated flow in randomly heterogeneous soils [J].
Amir, O ;
Neuman, SP .
TRANSPORT IN POROUS MEDIA, 2001, 44 (02) :355-383
[2]  
[Anonymous], 2014, RECENT ADV RADIAL BA
[3]  
[Anonymous], 2015, NUMERICAL METHODS FR
[4]  
Atluri SN, 2000, INT J NUMER METH ENG, V47, P537, DOI 10.1002/(SICI)1097-0207(20000110/30)47:1/3<537::AID-NME783>3.0.CO
[5]  
2-E
[6]   Adaptive low-rank approximation of collocation matrices [J].
Bebendorf, M ;
Rjasanow, S .
COMPUTING, 2003, 70 (01) :1-24
[7]   Application of a fractional advection-dispersion equation [J].
Benson, DA ;
Wheatcraft, SW ;
Meerschaert, MM .
WATER RESOURCES RESEARCH, 2000, 36 (06) :1403-1412
[8]   Anomalous transport in laboratory-scale, heterogeneous porous media [J].
Berkowitz, B ;
Scher, H ;
Silliman, SE .
WATER RESOURCES RESEARCH, 2000, 36 (01) :149-158
[9]   Physical pictures of transport in heterogeneous media: Advection-dispersion, random-walk, and fractional derivative formulations [J].
Berkowitz, B ;
Klafter, J ;
Metzler, R ;
Scher, H .
WATER RESOURCES RESEARCH, 2002, 38 (10)
[10]   Modeling non-Fickian transport in geological formations as a continuous time random walk [J].
Berkowitz, Brian ;
Cortis, Andrea ;
Dentz, Marco ;
Scher, Harvey .
REVIEWS OF GEOPHYSICS, 2006, 44 (02)