Stability and Error Analysis of an Efficient Numerical Method for Convection Dominated Parabolic PDEs with Jump Discontinuity in Source Function on Modified Layer-Adapted Mesh
{"title":"Stability and Error Analysis of an Efficient Numerical Method for Convection Dominated Parabolic PDEs with Jump Discontinuity in Source Function on Modified Layer-Adapted Mesh","authors":"Narendra Singh Yadav, Kaushik Mukherjee","doi":"10.1134/s0965542524030102","DOIUrl":null,"url":null,"abstract":"<h3 data-test=\"abstract-sub-heading\">Abstract</h3><p>This paper introduces a novel approach for analyzing efficient numerical solution of a class of singularly perturbed parabolic convection-diffusion PDEs having a finite jump discontinuity mostly in the source function at the interior of the spatial domain. These PDEs often appear in mathematical modeling of the semiconductor devices; and solutions of such problems usually exhibit a weak interior layer at one side of the point of discontinuity along with a boundary layer at one side of the spatial domain. The prime objectives of this paper are to overcome the theoretical challenge related to the monotonocity of the finite difference operator which utilizes one-sided second-order difference operators at the interface point; and also to establish higher-order numerical approximation in space, regardless of smaller and larger values of the parameter ε. Proving discrete maximum principle of the proposed difference operator is found to be challenging task on the standard Shishkin mesh adapted to both boundary and weak interior layers. We overcome this difficulty by constructing a special non-uniform mesh, called modified layer-adapted mesh; and hereby establish the stability as well as the parameter-uniform error estimate in the discrete supremum norm. Finally, the theoretical estimate is verified with numerical results for test examples with and without the exact solution. Moreover, numerical results are obtained for the semilinear PDEs and also compared with the implicit upwind scheme to exhibit the efficiency of the proposed algorithm.</p>","PeriodicalId":55230,"journal":{"name":"Computational Mathematics and Mathematical Physics","volume":"204 1","pages":""},"PeriodicalIF":0.7000,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational Mathematics and Mathematical Physics","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1134/s0965542524030102","RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
This paper introduces a novel approach for analyzing efficient numerical solution of a class of singularly perturbed parabolic convection-diffusion PDEs having a finite jump discontinuity mostly in the source function at the interior of the spatial domain. These PDEs often appear in mathematical modeling of the semiconductor devices; and solutions of such problems usually exhibit a weak interior layer at one side of the point of discontinuity along with a boundary layer at one side of the spatial domain. The prime objectives of this paper are to overcome the theoretical challenge related to the monotonocity of the finite difference operator which utilizes one-sided second-order difference operators at the interface point; and also to establish higher-order numerical approximation in space, regardless of smaller and larger values of the parameter ε. Proving discrete maximum principle of the proposed difference operator is found to be challenging task on the standard Shishkin mesh adapted to both boundary and weak interior layers. We overcome this difficulty by constructing a special non-uniform mesh, called modified layer-adapted mesh; and hereby establish the stability as well as the parameter-uniform error estimate in the discrete supremum norm. Finally, the theoretical estimate is verified with numerical results for test examples with and without the exact solution. Moreover, numerical results are obtained for the semilinear PDEs and also compared with the implicit upwind scheme to exhibit the efficiency of the proposed algorithm.
期刊介绍：
Computational Mathematics and Mathematical Physics is a monthly journal published in collaboration with the Russian Academy of Sciences. The journal includes reviews and original papers on computational mathematics, computational methods of mathematical physics, informatics, and other mathematical sciences. The journal welcomes reviews and original articles from all countries in the English or Russian language.
京公网安备 11010802042870号
京ICP备2023020795号-1
分享
求助内容：
应助结果提醒方式：
扫码关注我们
