Higher order numerical methods for singular perturbation problems

Philosophiae Doctor - PhDIn recent years, there has been a great interest towards the higher order numerical methods for singularly perturbed problems. As compared to their lower order counterparts, they provide better accuracy with fewer mesh points. Construction and/or implementation of direct higher order methods is usually very complicated. Thus a natural choice is to use some convergence acceleration techniques, e.g., Richardson extrapolation, defect correction, etc. In this thesis, we will consider various classes of problems described by singularly perturbed ordinary and partial differential equations. For these problems, we design some novel numerical methods and attempt to increase their accuracy as well as the order of convergence. We also do the same for existing numerical methods in some instances. We find that, even though the Richardson extrapolation technique always improves the accuracy, it does not perform equally well when applied to different methods for certain classes of problems. Moreover, while in some cases it improves the order of convergence, in other cases it does not. These issues are discussed in this thesis for linear and nonlinear singularly perturbed ODEs as well as PDEs. Extrapolation techniques are analyzed thoroughly in all the cases, whereas the limitations of the defect correction approach for certain problems is indicated at the end of the thesis.South Afric

A fitted numerical method for parabolic turning point singularly perturbed problems with an interior layer

The objective of this paper is to construct and analyzea fitted operator finite difference method (FOFDM) forthe family of time-dependent singularly perturbed parabolicconvection–diffusion problems. The solution to the problemswe consider exhibits an interior layer due to the presence ofa turning point. We first establish sharp bounds on the solu-tion and its derivatives. Then, we discretize the time variableusing the classical Euler method. This results in a system ofsingularly perturbed interior layer two-point boundary valueproblems. We propose a FOFDM to solve the system above

Using the Kellogg-Tsan Solution Decomposition in NumericalMethods for Singularly Perturbed Convection-Diffusion Problems

The linear one-dimensional singularly perturbed convection-diffusion problem is solved numerically by a second-order method that is uniform in the perturbation parameter . The method uses the Kellogg-Tsan decomposition of the continuous solution. This increases the accuracy of the numerical results and simplifies the proof of their -uniformit

On the design and implementation of a hybrid numerical method for singularly perturbed two-point boundary value problems

>Magister Scientiae - MScWith the development of technology seen in the last few decades, numerous solvers have been developed to provide adequate solutions to the problems that model different aspects of science and engineering. Quite often, these solvers are tailor-made for specific classes of problems. Therefore, more of such must be developed to accompany the growing need for mathematical models that help in the understanding of the contemporary world. This thesis treats two point boundary value singularly perturbed problems. The solution to this type of problem undergoes steep changes in narrow regions (called boundary or internal layer regions) thus rendering the classical numerical procedures inappropriate. To this end, robust numerical methods such as finite difference methods, in particular fitted mesh and fitted operator methods have extensively been used. While the former consists of transforming the continuous problem into a discrete one on a non-uniform mesh, the latter involves a special discretisation of the problem on a uniform mesh and are known to be more accurate. Both classes of methods are suitably designed to accommodate the rapid change(s) in the solution. Quite often, finite difference methods on piece-wise uniform meshes (of Shishkin-type) are adopted. However, methods based on such non-uniform meshes, though layer-resolving, are not easily extendable to higher dimensions. This work aims at investigating the possibility of capitalising on the advantages of both fitted mesh and fitted operator methods. Theoretical results are confirmed by extensive numerical simulations

A uniformly accurate finite elements method for singular perturbation problems

We consider piecewise polynomial finite elements method for a singular perturbation problem. The finite elements method of Griffiths for a problem with non-constant coefficients was adapted by introducing piecewise polynomial approximation. We generate the tridiagonal difference schemes which are second order accurate in uniform norm

An exponentially fitted finite difference scheme for a class of singularly perturbed delay differential equations with large delays

AbstractThis paper deals with singularly perturbed boundary value problem for a linear second order delay differential equation. It is known that the classical numerical methods are not satisfactory when applied to solve singularly perturbed problems in delay differential equations. In this paper we present an exponentially fitted finite difference scheme to overcome the drawbacks of the corresponding classical counter parts. The stability of the scheme is investigated. The proposed scheme is analyzed for convergence. Several linear singularly perturbed delay differential equations have been solved and the numerical results are presented to support the theory