Abstract
In this paper, we consider a class of singularly perturbed convection-diffusion boundary-value problems with discontinuous convection coefficient which often occur as mathematical models for analyzing shock wave phenomena in gas dynamics. In general, interior layers appear in the solutions of this class of problems and this gives rise to difficulty while solving such problems using the classical numerical methods (standard central difference or standard upwind scheme) on uniform meshes when the perturbation parameter ε is small. To achieve better numerical approximation in solving this class of problems, we propose a new hybrid scheme utilizing a layer-resolving piecewise-uniform Shishkin mesh and the method is shown to be ε-uniformly stable. In addition to this, it is proved that the proposed numerical scheme is almost second-order uniformly convergent in the discrete supremum norm with respect to the parameter ε. Finally, extensive numerical experiments are conducted to support the theoretical results. Further, the numerical results obtained by the newly proposed scheme are also compared with the hybrid scheme developed in the paper [Z.Cen, Appl. Math. Comput., 169(1): 689-699, 2005]. It shows that the current hybrid scheme exhibits a significant improvement over the hybrid scheme developed by Cen, in terms of the parameter-uniform order of convergence.
Highlights
A class of singularly perturbed boundary-value problems (BVPs) with discontinuous convection coefficient is considered on the unit interval
The prime objective of this article is to develop a new hybrid scheme on a piecewise-uniform Shishkin mesh for solving the singularly perturbed BVPs of the form (1.1)–(1.2) so that the method is at least second-order uniformly convergent with respect to ε in the discrete supremum norm
(1.2), we describe the proposed hybrid numerical scheme which consists of a modified central difference scheme when ε > 2 a N −1 and a proper combination of the midpoint upwind scheme in the outer regions (0, ξ − σ1], [ξ + σ2, 1)
Summary
A class of singularly perturbed boundary-value problems (BVPs) with discontinuous convection coefficient is considered on the unit interval. The prime objective of this article is to develop a new hybrid scheme on a piecewise-uniform Shishkin mesh for solving the singularly perturbed BVPs of the form (1.1)–(1.2) so that the method is at least second-order uniformly convergent with respect to ε in the discrete supremum norm.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have