Analysis of a higher-order scheme for multi-term time-fractional integro-partial differential equations with multi-term weakly singular kernels

Santra, S (2024) Analysis of a higher-order scheme for multi-term time-fractional integro-partial differential equations with multi-term weakly singular kernels. In: Numerical Algorithms .

PDF Num_Alg_24.pdf - Published Version Restricted to Registered users only Download (1MB) | Request a copy

Abstract

This work is focused on developing a hybrid numerical method that combines a higher-order finite difference method and multi-dimensional Hermite wavelets to address two-dimensional multi-term time-fractional integro-partial differential equations with multi-term weakly singular kernels having bounded and unbounded time derivatives at the initial time t=0. Specifically, the multi-term fractional operators are discretized using a higher-order approximation designed by employing different interpolation schemes based on linear, quadratic, and cubic interpolation leading to O(N-(4-α1)) accuracy on a suitably chosen nonuniform mesh and O(N-α1) accuracy on a uniformly distributed mesh. The weakly singular integral operators are approximated by a modified numerical quadrature, which is a combination of the composite trapezoidal approximation and the midpoint rule. The effects of the exponents of the weakly singular kernels over fractional orders are analyzed in terms of accuracy over uniform and nonuniform meshes for the solution having both bounded and unbounded time derivatives. The stability of the proposed semi-discrete scheme is derived based on L�-norm for uniformly distributed temporal mesh. Further, we employ the uniformly distributed collocation points in spatial directions to estimate the tensor-based wavelet coefficients. Moreover, the convergence analysis of the fully discrete scheme is carried out based on L2-norm leading to O(N-α1) accuracy on a uniform mesh. It also highlights the higher-order accuracy over nonuniform mesh. Additionally, we discuss the convergence analysis of the proposed scheme in the context of the multi-term time-fractional diffusion equations involving time singularity demonstrating a O(N-(4-α1)) accuracy on a nonuniform mesh with suitably chosen grading parameter. Note that the scheme reduces to O(N-α1) accuracy on a uniform mesh. Several tests are performed on numerous examples in L�- and L2-norm to show the efficiency of the proposed method. Further, the solutions� nature and accuracy in terms of absolute point-wise error are illustrated through several isosurface plots for different regularities of the exact solution. These experiments confirm the theoretical accuracy and guarantee the convergence of approximations to the functions having time singularity, and the higher-order accuracy for a suitably chosen nonuniform mesh. © The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature 2024.

Item Type:	Journal Article
Publication:	Numerical Algorithms
Publisher:	Springer
Additional Information:	The copyright for this article belongs to the Publisher.
Department/Centre:	Division of Interdisciplinary Sciences > Computational and Data Sciences
Date Deposited:	20 Sep 2024 11:36
Last Modified:	20 Sep 2024 11:43
URI:	http://eprints.iisc.ac.in/id/eprint/86279

Actions (login required)

View Item