Robust higher-order numerical scheme for solving time-fractional singularly perturbed parabolic partial differential equations with large delay in time
This paper presents a robust higher-order numerical scheme for time-fractional singularly perturbed partial differential equations having large delay in time, where the time-fractional derivative term is taken in the Caputo sense with order α∈(0,1). The time-fractional singularly perturbed delay par...
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Elsevier
2025-07-01
|
| Series: | International Journal of Thermofluids |
| Subjects: | |
| Online Access: | http://www.sciencedirect.com/science/article/pii/S2666202725002642 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Summary: | This paper presents a robust higher-order numerical scheme for time-fractional singularly perturbed partial differential equations having large delay in time, where the time-fractional derivative term is taken in the Caputo sense with order α∈(0,1). The time-fractional singularly perturbed delay partial differential equation is the class of time-fractional delay partial differential equations in which the highest-order space derivative is multiplied by a very small positive parameter. The problem is discretized by the Alikhanov’s L2−1σ technique on a uniform mesh in the temporal direction and a hybrid finite difference scheme on a piecewise uniform mesh in the spatial direction. The hybrid scheme is composed of the mid-point upwind scheme in the outer region and central finite difference method in the boundary layer region. The uniform stability analysis and the bounds of the truncation error are performed. The convergence of the numerical scheme is proved in the maximum norm. It is shown that the proposed numerical scheme is parameter-uniformly convergent of order ON−γ+M−2ln2M, where N and M are the number of time and space mesh intervals, respectively and γ=3−α. To verify the applicability of the proposed numerical scheme, three examples are considered and the results agree with the theoretical concepts discussed. Moreover, comparisons are conducted with some of the methods available in literature and the present scheme outperforms from these methods. Finally, the numerical results verify that the proposed scheme is higher-order and ϵ-uniformly convergent. |
|---|---|
| ISSN: | 2666-2027 |