Fitted non-polynomial spline method for singularly perturbed differential difference equations with integral boundary condition

The aim of this paper is to present fitted non-polynomial spline method for singularly perturbed differential-difference equations with integral boundary condition. The stability and uniform convergence of the proposed method are proved. To validate the applicability of the scheme, two model problems are considered for numerical experimentation and solved for different values of the perturbation parameter, ε and mesh size, h. The numerical results are tabulated in terms of maximum absolute errors and rate of convergence and it is observed that the present method is more accurate and uniformly convergent for h ≥ ε where the classical numerical methods fails to give good result and it also improves the results of the methods existing in the literature

A seventh order numerical method for singular perturbed differential-difference equations with negative shift

In this paper, a seventh order numerical method is presented for solving singularly perturbed differential-difference equations with negative shift. In recent papers the term negative shift has been used for delay. Such problems are associated with expected first exit time problem of the membrane, potential in models for neuron and in variational problems in control theory. In the numerical treatment for such type of boundary value problems, we first use Taylor approximation to tackle terms containing small shifts which converts into a singularly perturbed boundary value problem. This two point boundary value problem is transformed into general first order ordinary differential equation system. A discrete approximation of a seventh order compact difference scheme is employed for the first order system and solved by using the boundary conditions. Several numerical examples are solved and compared with exact solution. We also present least square errors, maximum errors and observed that the present method approximates the exact solution very well

Hybrid Algorithm for Singularly Perturbed Delay Parabolic Partial Differential Equations

This study aims at constructing a numerical scheme for solving singularly perturbed parabolic delay differential equations. Taylor’s series expansion is applied to approximate the shift term. The obtained result is approximated by using the implicit Euler method in the temporal discretization on a uniform step size with the hybrid numerical scheme consisting of the midpoint upwind method in the outer layer region and the cubic spline method in the inner layer region on a piecewise uniform Shishkin mesh in the spatial discretization. The constructed scheme is an ε−uniformly convergent accuracy of order one. Some test examples are considered to testify the theoretical investigations

A nonstandard fitted operator finite difference method for two-parameter singularly perturbed time-delay parabolic problems

In this article, a class of singularly perturbed time-delay two-parameter second-order parabolic problems are considered. The presence of the two small parameters attached to the derivatives causes the solution of the given problem to exhibit boundary layer(s). We have developed a uniformly convergent nonstandard fitted operator finite difference method (NSFOFDM) to solve the considered problems. The Crank-Nicolson scheme with a uniform mesh is used for the discretization of the time derivative, while for the spatial discretization, we have applied a fitted operator finite difference method following the nonstandard methodology of Mickens. Moreover, the solution bounds of the governing equation are shown by asymptotic analysis. The convergence of the proposed numerical scheme is investigated using truncation error and the barrier function approach. The study shows that our proposed scheme is uniformly convergent independent of the perturbation parameters, quadratically in time, and linearly in space. Numerical experiments are carried out, and the results are presented in tables and graphically

A numerical approach for a two-parameter singularly perturbed weakly-coupled system of 2-D elliptic convection–reaction–diffusion PDEs

In this work, we consider the numerical approximation of a two dimensional elliptic singularly perturbed weakly-coupled system of convection–reaction–diffusion type, which has two different parameters affecting the diffusion and the convection terms, respectively. The solution of such problems has, in general, exponential boundary layers as well as corner layers. To solve the continuous problem, we construct a numerical method which uses a finite difference scheme defined on an appropriate layer-adapted Bakhvalov–Shishkin mesh. Then, the numerical scheme is a first order uniformly convergent method with respect both convection and diffusion parameters. Numerical results obtained with the algorithm for some test problems are presented, which show the best performance of the proposed method, and they also corroborate in practice the theoretical analysis

A hybrid approximation scheme for 1-D singularly perturbed parabolic convection-diffusion problems

Our study is concerned with a hybrid spectral collocation approach to solving singularly perturbed 1-D parabolic convection-diffusion problems. In this approach, discretization in time is carried out with the help of Taylor series expansions before the spectral based on novel special polynomials is applied to the spatial operator in the time step. A detailed error analysis of the presented technique is conducted with regard to the space variable. The advantages of this attempt are presented through comparison of our results in the model problems obtained by this technique and other existing schemes

A hybrid approximation scheme for 1-D singularly perturbed parabolic convection-diffusion problems

Our study is concerned with a hybrid spectral collocation approach to solving singularly perturbed 1-D parabolic convection-diffusion problems. In this approach, discretization in time is carried out with the help of Taylor series expansions before the spectral based on novel special polynomials is applied to the spatial operator in the time step. A detailed error analysis of the presented technique is conducted with regard to the space variable. The advantages of this attempt are presented through comparison of our results in the model problems obtained by this technique and other existing schemes

On the numerical integration of singularly perturbed Volterra integro-differential equations

Magister Scientiae - MScEfficient numerical approaches for parameter dependent problems have been an inter- esting subject to numerical analysts and engineers over the past decades. This is due to the prominent role that these problems play in modeling many real life situations in applied sciences. Often, the choice and the e ciency of the approaches depend on the nature of the problem to solve. In this work, we consider the general linear first-order singularly perturbed Volterra integro-differential equations (SPVIDEs). These singularly perturbed problems (SPPs) are governed by integro-differential equations in which the derivative term is multiplied by a small parameter, known as "perturbation parameter". It is known that when this perturbation parameter approaches zero, the solution undergoes fast transitions across narrow regions of the domain (termed boundary or interior layer) thus affecting the convergence of the standard numerical methods. Therefore one often seeks for numerical approaches which preserve stability for all the values of the perturbation parameter, that is "numerical methods. This work seeks to investigate some "numerical methods that have been used to solve SPVIDEs. It also proposes alternative ones. The various numerical methods are composed of a fitted finite difference scheme used along with suitably chosen interpolating quadrature rules. For each method investigated or designed, we analyse its stability and convergence. Finally, numerical computations are carried out on some test examples to con rm the robustness and competitiveness of the proposed methods

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

Parameter-uniformly convergent numerical scheme for singularly perturbed delay parabolic differential equation via extended B-spline collocation

This paper presents a parameter-uniform numerical method to solve the time dependent singularly perturbed delay parabolic convection-diffusion problems. The solution to these problems displays a parabolic boundary layer if the perturbation parameter approaches zero. The retarded argument of the delay term made to coincide with a mesh point and the resulting singularly perturbed delay parabolic convection-diffusion problem is approximated using the implicit Euler method in temporal direction and extended cubic B-spline collocation in spatial orientation by introducing artificial viscosity both on uniform mesh. The proposed method is shown to be parameter uniform convergent, unconditionally stable, and linear order of accuracy. Furthermore, the obtained numerical results agreed with the theoretical results