A parameter-uniform numerical method for a Sobolev problem with initial layer

The present study is concerned with the numerical solution, using finite difference method of a one-dimensional initial-boundary value problem for a linear Sobolev or pseudo-parabolic equation with initial jump. In order to obtain an efficient method, to provide good approximations with independence of the perturbation parameter, we have developed a numerical method which combines a finite difference spatial discretization on uniform mesh and the implicit rule on Shishkin mesh(S-mesh) for the time variable. The fully discrete scheme is shown to be convergent of order two in space and of order one expect for a logarithmic factor in time, uniformly in the singular perturbation parameter. Some numerical results confirming the expected behavior of the method are shown.      

On adaptive mesh for the initial boundary value singularly perturbed delay Sobolev problems

A uniform finite difference method on a B-mesh is applied to solve the initial-boundary value problem for singularly perturbed delay Sobolev equations. To solve the foresold problem, finite difference scheme on a special nonuniform mesh, whose solution converges point-wise independently of the singular perturbation parameter is constructed and analyzed. The present paper also aims at discussing the stability and convergence analysis of the method. An error analysis shows that the method is of second order convergent in the discrete maximum norm independent of the perturbation parameter. A numerical example and the simulation results show the effectiveness of our theoretical results.     

A Fitted Numerov Method for Singularly Perturbed Parabolic Partial Differential Equation with a Small Negative Shift Arising in Control Theory

In this paper, a fitted Numerov method is constructed for a class of singularly perturbed one-dimensional parabolic partial differential equations with a small negative shift in the temporal variable. Similar boundary value problems are associated with a furnace used to process a metal sheet in control theory. Here, the study focuses on the effect of shift on the boundary layer behavior of the solution via finite difference approach. When the shift parameter is smaller than the perturbation parameter, the shifted term is expanded in Taylor series and an exponentially fitted tridiagonal finite difference scheme is developed. The proposed finite difference scheme is unconditionally stable. When the shift parameter is larger than the perturbation parameter, a special type of mesh is used for the temporal variable so that the shift lies on the nodal points and an exponentially fitted scheme is developed. This scheme is also unconditionally stable. The applicability of the proposed methods is demonstrated by means of two examples.     

Numerical solution of singularly perturbed convection-diffusion problem using parameter uniform B-spline collocation method

This paper is concerned with a numerical scheme to solve a singularly perturbed convection-diffusion problem. The solution of this problem exhibits the boundary layer on the right-hand side of the domain due to the presence of singular perturbation parameter ε. The scheme involves B-spline collocation method and appropriate piecewise-uniform Shishkin mesh. Bounds are established for the derivative of the analytical solution. Moreover, the present method is boundary layer resolving as well as second-order uniformly convergent in the maximum norm. A comprehensive analysis has been given to prove the uniform convergence with respect to singular perturbation parameter. Several numerical examples are also given to demonstrate the efficiency of B-spline collocation method and to validate the theoretical aspects.     

High order parameter-robust numerical method for singularly perturbed reaction-diffusion problems

We present a high order parameter-robust finite difference method for singularly perturbed reaction-diffusion problems. The problem is discretized using a suitable combination of fourth order compact difference scheme and central difference scheme on generalized Shishkin mesh. The convergence analysis is given and the method is proved to be almost fourth order uniformly convergent in maximum norm with respect to singular perturbation parameter ε. Numerical experiments are conducted to demonstrate the theoretical results.     

A uniformly convergent B-spline collocation method on a nonuniform mesh for singularly perturbed one-dimensional time-dependent linear convection–diffusion problem

A numerical method is proposed for solving singularly perturbed one-dimensional parabolic convection–diffusion problems. The method comprises a standard implicit finite difference scheme to discretize in temporal direction on a uniform mesh by means of Rothe's method and B-spline collocation method in spatial direction on a piecewise uniform mesh of Shishkin type. The method is shown to be unconditionally stable and accurate of order O((Δx)²+Δt). An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. Several numerical experiments have been carried out in support of the theoretical results. Comparisons of the numerical solutions are performed with an upwind finite difference scheme on a piecewise uniform mesh and exponentially fitted method on a uniform mesh to demonstrate the efficiency of the method.     

A numerical method for singularly perturbed weakly coupled system of two second order ordinary differential equations with discontinuous source term

In this paper, a numerical method based on finite difference scheme and Shishkin mesh for singularly perturbed two second order weakly coupled system of ordinary differential equations with discontinuous source term is presented. An error estimate is derived to show that the method is uniformly convergent with respect to the singular perturbation parameter. Numerical results are presented to illustrate the theoretical results.     

Fourth-order finite-difference method for boundary value problems with two small parameters

We present a finite difference scheme for a class of linear singularly perturbed boundary value problems with two small parameters. The problem is discretized using a Bakhvalov-type mesh. It is proved under certain conditions that this scheme is fourth-order accurate and that its error does not increase when the perturbation parameter tends to zero. Numerical examples are presented which demonstrate computationally the fourth order of the method.     

Uniform Global Convergence of a Hybrid Scheme for Singularly Perturbed Reaction–Diffusion Systems

We consider a system of coupled singularly perturbed reaction–diffusion two-point boundary-value problems. A hybrid difference scheme on a piecewise-uniform Shishkin mesh is constructed for solving this system, which generates better approximations to the exact solution than the classical central difference scheme. Moreover, we prove that the method is third order uniformly convergent in the maximum norm when the singular perturbation parameter is small. Numerical experiments are conducted to validate the theoretical results.     

Optimal error estimate of upwind scheme on Shishkin-type meshes for singularly perturbed parabolic problems with discontinuous convection coefficients

This paper analyzes the implicit upwind finite difference scheme on Shishkin-type meshes (including the classical piecewise-uniform Shishkin mesh and the Bakhalov-Shishkin mesh) for a class of singularly perturbed parabolic convection-diffusion problems exhibiting strong interior layers. Suitable conditions on the mesh-generating functions are derived and are found to be sufficient for the convergence of the method, uniformly with respect to the perturbation parameter. Utilizing these conditions, it is shown that the method converges uniformly in the discrete supremum norm with an optimal error bound. Numerical results are presented to validate the theoretical results.     

A second order scheme for a time-dependent,singularly perturbed convection–diffusion equation

We consider a numerical scheme for a one-dimensional, time-dependent, singularly perturbed convection–diffusion problem. The problem is discretized in space by a standard finite element method on a Bakhvalov–Shishkin type mesh. The space error is measured in an L² norm. For the time integration, the implicit midpoint rule is used. The fully discrete scheme is shown to be convergent of order 2 in space and time, uniformly in the singular perturbation parameter.     

Finite-Difference Methods for a Class of Strongly Nonlinear Singular Perturbation Problems

The paper is concerned with strongly nonlinear singularly perturbed bound- ary value problems in one dimension.The problems are solved numerically by finite- difference schemes on special meshes which are dense in the boundary layers.The Bakhvalov mesh and a special piecewise equidistant mesh are analyzed.For the central scheme,error estimates are derived in a discrete L~1 norm.They are of second order and decrease together with the perturbation parameterε.The fourth-order Numerov scheme and the Shishkin mesh are also tested numerically.Numerical results showε-uniform pointwise convergence on the Bakhvalov and Shishkin meshes.     

Uniform convergence analysis of hybrid numerical scheme for singularly perturbed problems of mixed type

In this article, we consider a class of singularly perturbed mixed parabolic‐elliptic problems whose solutions possess both boundary and interior layers. To solve these problems, a hybrid numerical scheme is proposed and it is constituted on a special rectangular mesh which consists of a layer resolving piecewise‐uniform Shishkin mesh in the spatial direction and a uniform mesh in the temporal direction. The domain under consideration is partitioned into two subdomains. For the spatial discretization, the proposed scheme is comprised of the classical central difference scheme in the first subdomain and a hybrid finite difference scheme in the second subdomain, whereas the time derivative in the given problem is discretized by the backward‐Euler method. We prove that the method converges uniformly with respect to the perturbation parameter with almost second‐order spatial accuracy in the discrete supremum norm. Numerical results are finally presented to validate the theoretical results.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1931–1960, 2014     

On piecewise-uniform meshes for upwind- and central-difference operators for solving singularly perturbed problems

Numerical methods composed of upwind-difference operators onuniform meshes are shown analytically to be defective for solvingsingularly perturbed differential equations, in the sense that,as the mesh is refined, the error in the numerical approximationincreases until the mesh parameter has decreased to the sameorder of magnitude as the singular perturbation parameter. Itis also shown that the same is true for upwind-difference operatorson piecewise uniform meshes having a transition point that isdependent solely on the singular perturbation parameter. Itis then shown, with specific analytic examples, that both upwind-and central-difference operators on specially designed piecewise-uniformmeshes give numerical methods which do not suffer from thisdefect. Conditions are also given on the structure of a piecewiseuniform mesh that are necessary if the numerical method, composedof this mesh and an upwind-difference operator, is to be convergentuniformly with respect to the singular perturbation parameter.     

A singular perturbation approach to solve Burgers–Huxley equation via monotone finite difference scheme on layer-adaptive mesh

This paper deals with a numerical method for solving one-dimensional unsteady Burgers–Huxley equation with the viscosity coefficient ε. The parameter ε takes any values from the half open interval (0, 1]. At small values of the parameter ε, an outflow boundary layer is produced in the neighborhood of right part of the lateral surface of the domain and the problem can be considered as a non-linear singularly perturbed problem with a singular perturbation parameter ε. Using singular perturbation analysis, asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular components. We construct a numerical scheme that comprises of implicit-Euler method to discretize in temporal direction on uniform mesh and a monotone hybrid finite difference operator to discretize the spatial variable with piecewise uniform Shishkin mesh. To obtain better accuracy, we use central finite difference scheme in the boundary layer region. Shishkin meshes are refined in the boundary layer region, therefore stability constraint is satisfied by proposed scheme. Quasilinearization process is used to tackle the non-linearity and it is shown that quasilinearization process converges quadratically. The method has been shown to be first order uniformly accurate in the temporal variable, and in the spatial direction it is first order parameter uniform convergent in the outside region of boundary layer, and almost second order parameter uniform convergent in the boundary layer region. Accuracy and uniform convergence of the proposed method is demonstrated by numerical examples and comparison of numerical results made with the other existing methods.     

Fitted finite difference method for singularly perturbed delay differential equations

This paper deals with singularly perturbed initial value problem for linear second-order delay differential equation. An exponentially fitted difference scheme is constructed in an equidistant mesh, which gives first order uniform convergence in the discrete maximum norm. The difference scheme is shown to be uniformly convergent to the continuous solution with respect to the perturbation parameter. A numerical example is solved using the presented method and compared the computed result with exact solution of the problem.     

A compact finite difference scheme for 2D reaction–diffusion singularly perturbed problems

In this work we define a compact finite difference scheme of positive type to solve a class of 2D reaction–diffusion elliptic singularly perturbed problems. We prove that if the new scheme is constructed on a piecewise uniform mesh of Shishkin type, it provides better approximations than the classical central finite difference scheme. Moreover, the uniform parameter bound of the error shows that the scheme is third order convergent in the maximum norm when the singular perturbation parameter is sufficiently small. Some numerical experiments illustrate in practice the result of convergence proved theoretically.     

An ?-uniform hybrid scheme for singularly perturbed delay differential equations

In this paper, a singularly perturbed delay differential equation of first order has been considered. The problem is solved by using a hybrid scheme on a Shishkin mesh. The difference scheme is shown to converge to the continuous solution uniformly with respect to the perturbation parameter. Truncation errors are obtained. Finally, numerical experiments are carried out on a test problem, confirming the effectiveness of the proposed technique.     

A higher order difference method for singularly perturbed parabolic partial differential equations

This paper studies a higher order numerical method for the singularly perturbed parabolic convection-diffusion problems where the diffusion term is multiplied by a small perturbation parameter. In general, the solutions of these type of problems have a boundary layer. Here, we generate a spatial adaptive mesh based on the equidistribution of a positive monitor function. Implicit Euler method is used to discretize the time variable and an upwind scheme is considered in space direction. A higher order convergent solution with respect to space and time is obtained using the postprocessing based extrapolation approach. It is observed that the convergence is independent of perturbation parameter. This technique enhances the order of accuracy from first order uniform convergence to second order uniform convergence in space as well as in time. Comparative study with the existed meshes show the highly effective behavior of the present method.     

在非均匀网格上解椭圆—双曲型偏微分方程奇异摄动问题

本文考察了椭圆一双曲型偏微分方程奇异摄动问题(1.1),证明了迎风差分格式在一特殊的非均匀网格上是一阶一致收敛的.最后给出了一些数值结果.