In this work, we theoretically and numerically discuss the time fractional subdiffusion-normal transport equation, which depicts a crossover from sub-diffusion (as $t\rightarrow 0$) to normal diffusion (as $t\rightarrow \infty$). Firstly, the well-posedness and regularities of the model are studied by using the bivariate Mittag-Leffler function. Theoretical results show that after introducing the first-order derivative operator, the regularity of the solution can be improved in substance. Then, a numerical scheme with high-precision is developed no matter the initial value is smooth or non-smooth. More specifically, we use the contour integral method (CIM) with parameterized hyperbolic contour to approximate the temporal local and non-local operators, and employ the standard Galerkin finite element method for spacial discretization. Rigorous error estimates show that the proposed numerical scheme has spectral accuracy in time and optimal convergence order in space. Besides, we further improve the algorithm and reduce the computational cost by using the barycentric Lagrange interpolation. Finally, the obtained theoretical results as well as the acceleration algorithm are verified by several 1-D and 2-D numerical experiments, which also show that the numerical scheme developed in this paper is effective and robust.
翻译:在这项工作中,我们从理论上和数字上讨论了时间分数子子子扩散-正常运输等式,它描述了从次扩散(作为$t\rightrow 0.0美元)到正常扩散(作为$t\rightrow \ infty$)的交叉跨度。首先,模型的精度和规律性是通过使用双变量Mittag-Leffler函数来研究的。理论结果显示,在引入第一级衍生器操作器后,解决方案的规律性在实质上是可以改进的。然后,一个具有高精确度的数值方法,无论初始值是平滑还是非滑动。更具体地说,我们使用带有参数化超偏差的等离子集成法集成法(CIM)到正常扩散(作为参数化的超偏差等同度)到与当地和非当地操作器的接近,并使用标准加勒金定元素元素法用于慢离子化。严格错误估计表明,拟议的数字方法在时间和最佳融合顺序上是可以改进的。此外,我们还改进了算法,并降低了计算成本,通过使用巴中中心拉格-D中心的计算法和数级之间测算法,最后通过数法显示了数率的加速率。