Parameter-uniform finite difference schemes for singularly perturbed parabolic diffusion-convection-reaction problems

In this paper, parameter-uniform numerical methods for a class of singularly perturbed parabolic partial differential equations with two small parameters on a rectangular domain are studied. Parameter-explicit theoretical bounds on the derivatives of the solutions are derived. The solution is decomposed into a sum of regular and singular components. A numerical algorithm based on an upwind finite difference operator and an appropriate piecewise uniform mesh is constructed. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are given to illustrate the parameter-uniform convergence of the numerical approximations.

     

Fitted mesh method for singularly perturbed reaction-convection-diffusion problems with boundary and interior layers

A robust numerical method for a singularly perturbed secondorder ordinary differential equation having two parameters with a discontinuous source term is presented in this article. Theoretical bounds are derived for the derivatives of the solution and its smooth and singular components. An appropriate piecewise uniform mesh is constructed, and classical upwind finite difference schemes are used on this mesh to obtain the discrete system of equations. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are provided to illustrate the convergence of the numerical approximations.     

Approximation of derivative in a system of singularly perturbed convection-diffusion equations

In this paper, a numerical method for a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative is presented. Parameter-uniform error bounds for the numerical solution and also to numerical derivative are established. Numerical results are provided to illustrate the theoretical results.     

A technique to prove parameter-uniform convergence for a singularly perturbed convection–diffusion equation

A priori parameter explicit bounds on the solution of singularly perturbed elliptic problems of convection–diffusion type are established. Regular exponential boundary layers can appear in the solution. These bounds on the solutions and its derivatives are obtained using a suitable decomposition of the solution into regular and layer components. By introducing extensions of the coefficients to a larger domain, artificial compatibility conditions are not imposed in the derivation of these decompositions.     

Numerical analysis of a strongly coupled system of two singularly perturbed convection–diffusion problems

A system of two coupled singularly perturbed convection–diffusion ordinary differential equations is examined. The diffusion term in each equation is multiplied by a small parameter, and the equations are coupled through their convective terms. The problem does not satisfy a conventional maximum principle. Its solution is decomposed into regular and layer components. Bounds on the derivatives of these components are established that show explicitly their dependence on the small parameter. A numerical method consisting of simple upwinding and an appropriate piecewise-uniform Shishkin mesh is shown to generate numerical approximations that are essentially first order convergent, uniformly in the small parameter, to the true solution in the discrete maximum norm.      

Parameter-uniform numerical method for singularly perturbed convection-diffusion problem on a circular domain

A linear singularly perturbed elliptic problem, of convection-diffusion type, posed on a circular domain is examined. Regularity constraints are imposed on the data in the vicinity of the two characteristic points. The solution is decomposed into a regular and a singular component. A priori parameter-explicit pointwise bounds on the partial derivatives of these components are established. By transforming to polar co-ordinates, a monotone finite difference method is constructed on a piecewise-uniform layer-adapted mesh of Shishkin type. Numerical analysis is presented for this monotone numerical method. The numerical method is shown to be parameter-uniform. Numerical results are presented to illustrate the theoretical error bounds established.     

A parameter‐uniform scheme for the parabolic singularly perturbed problem with a delay in time

In this paper, a parameter‐uniform numerical scheme for the solution of singularly perturbed parabolic convection–diffusion problems with a delay in time defined on a rectangular domain is suggested. The presence of the small diffusion parameter ? leads to a parabolic right boundary layer. A collocation method consisting of cubic B ‐spline basis functions on an appropriate piecewise‐uniform mesh is used to discretize the system of ordinary differential equations obtained by using Rothe's method on an equidistant mesh in the temporal direction. The parameter‐uniform convergence of the method is shown by establishing the theoretical error bounds. The numerical results of the test problems validate the theoretical error bounds.     

An elliptic singularly perturbed problem with two parameters I: Solution decomposition

In this paper we consider a singularly perturbed elliptic problem with two small parameters posed on the unit square. Its solution may have exponential, parabolic and corner layers. We give a decomposition of the solution into regular and layer components and derive pointwise bounds on the components and their derivatives. The estimates are obtained by the analysis of appropriate problems on unbounded domains.     

A-posteriori error analysis and adaptive finite element methods for singularly perturbed convection-diffusion equations

For singularly perturbed one-dimensional convection-diffusion equations, finite element approximations are constructed based on a so-called approximate symmetrization of the given unsymmetric problem. Local a-posteriori error estimates are established with respect to an appropriate energy norm where the bounds are proved to be realistic. The local bounds, called error indicators, provide a basis for a self-adaptive mesh refinement. For a model problem numerical results are presented showing that the adaptive method detects and resolves the boundary layer.     

Convergence Analysis of Weighted Difference Approximations on Piecewise Uniform Grids to a Class of Singularly Perturbed Functional Differential Equations

We consider the numerical approximation of a singularly perturbed time delayed convection diffusion problem on a rectangular domain. Assuming that the coefficients of the differential equation be smooth, we construct and analyze a higher order accurate finite difference method that converges uniformly with respect to the singular perturbation parameter. The method presented is a combination of the central difference spatial discretization on a Shishkin mesh and a weighted difference time discretization on a uniform mesh. A?priori explicit bounds on the solution of the problem are established. These bounds on the solution and its derivatives are obtained using a suitable decomposition of the solution into regular and layer components. It is shown that the proposed method is $L_{2}^{h}$ -stable. The analysis done permits its extension to the case of adaptive meshes which may be used to improve the solution. Numerical examples are presented to demonstrate the effectiveness of the method. The convergence obtained in practical satisfies the theoretical predictions.     

A Parameter-Uniform Finite Difference Method for a Coupled System of Convection-Diffusion Two-Point Boundary Value Problems

A system of m (≥2) linear convection-diffusion two-point boundary value problems is examined,where the diffusion term in each equation is multiplied by a small parameterεand the equations are coupled through their convective and reactive terms via matrices B and A respectively.This system is in general singularly perturbed. Unlike the case of a single equation,it does not satisfy a conventional maximum princi- ple.Certain hypotheses are placed on the coupling matrices B and A that ensure exis- tence and uniqueness of a solution to the system and also permit boundary layers in the components of this solution at only one endpoint of the domain;these hypotheses can be regarded as a strong form of diagonal dominance of B.This solution is decomposed into a sum of regular and layer components.Bounds are established on these compo- nents and their derivatives to show explicitly their dependence on the small parameterε.Finally,numerical methods consisting of upwinding on piecewise-uniform Shishkin meshes are proved to yield numerical solutions that are essentially first-order conver- gent,uniformly inε,to the true solution in the discrete maximum norm.Numerical results on Shishkin meshes are presented to support these theoretical bounds.     

The necessity of Shishkin decompositions

In the present paper, we study model singularly perturbed convection-diffusion problems with exponential boundary layers. It has been believed for some time that only a complete splitting of the exact solution into regular and layer parts provides the information necessary for the study of the uniform convergence properties of numerical methods for these problems on layer-adapted grids (such as Shishkin meshes). In the present paper, we give new proofs of uniform interpolation error estimates for linear and bilinear interpolation; these proofs are based on the older a priori bounds derived by Kellogg and Tsan [1].     

Numerical solution of the high thermal loss problem presented by a fractional differential equation

Three different numerical methods are used to solve singularly perturbed Able Volterra integral equation as presented by a fractional differential equation. Convergence and stability analysis together with the results of these methods are compared and contrasted when applied to the high thermal loss problem as an example of singularly perturbed Able Volterra integral equation.     

一类四阶非线性奇摄动方程的边界层问题

本文研究了一类四阶非线性奇摄动方程的边界层问题,利用在左右边界层的两次匹配,得出了原问题解的一致有效的渐近表达式.这个结果是奇摄动理论在研究高阶微分方程中的一个应用.     

A new fitted operator finite difference method to solve systems of evolutionary reaction-diffusion equations

Abstract

In recent years, fitted operator finite difference methods (FOFDMs) have been developed for numerous types of singularly perturbed ordinary differential equations. The construction of most of these methods differed though the final outcome remained similar. The most crucial aspect was how the difference operator was designed to approximate the differential operator in question. Very often the approaches for constructing these operators had limited scope in the sense that it was difficult to extend them to solve even simple one-dimensional singularly perturbed partial differential equations. However, in some of our most recent work, we have successfully designed a class of FOFDMs and extended them to solve singularly perturbed time-dependent partial differential equations. In this paper, we design and analyze a robust FOFDM to solve a system of coupled singularly perturbed parabolic reaction-diffusion equations. We use the backward Euler method for the semi-discretization in time. An FOFDM is then developed to solve the resulting set of boundary value problems. The proposed method is analyzed for convergence. Our method is uniformly convergent with order one and two, respectively, in time and space, with respect to the perturbation parameters. Some numerical experiments supporting the theoretical investigations are also presented.     

Some Fredholm equations with boundary layer resonance

Certain singularly perturbed differential equations which exhibit boundary layer resonance are difficult to solve by the application of standard asymptotic methods. After reformulation as a singularly perturbed integral equation and treatment by a recently developed asymptotic methodology, the desired solution is obtained in a straightforward manner.     

A singular approach to a class of impulsive differential equation

In this paper, a singular approach to study the solutions of an impulsive differential equation from a qualitative and quantitative point of view is proposed. In the approach, a suitable singular perturbation term is introduced and a singularly perturbed system with infinite initial values is defined, in which, the reduced problem of the singularly perturbed system is exactly the impulsive differential equation under consideration. Then the boundary layer function method is applied to construct the uniformly valid asymptotic solutions to the singularly perturbed system. Based on the continuous asymptotic solution, the discontinuous solutions of the impulsive differential equation are described and approximated. An example, namely, a classical Lotka-Volterra prey-predator model with one pulse is carried out to illustrate the main results.     

A Shishkin mesh for a singularly perturbed Riccati equation

In this paper a singularly perturbed Riccati initial value problem is examined. Parameter explicit bounds on the solution and its derivatives are given. A numerical method composed of an implicit difference operator and a piecewise-uniform Shishkin mesh is constructed. A theoretical parameter independently bound on the errors in the numerical approximations is established. Numerical results are presented which are in agreement with the theoretical error bound.     

A nonasymptotic method for singular perturbation problems

In this paper, an approximate method for the numerical integration of singularly perturbed two-point boundary-value problems with a boundary layer on the left end of the underlying interval is presented. The method is distinguished by the following fact: the original second-order differential equation is replaced by an approximate first-order differential equation with a small deviating argument and is solved efficiently by employing the Simpson rule, coupled with the discrete invariant imbedding algorithm. The proposed method is iterative on the deviating argument. Several numerical examples have been solved to demonstrate the applicability of the method.     

A numerical method for a system of equations with a small parameter and a point source on an infinite interval

A method for solving a boundary-value problem on an infinite interval is considered for a linear system of second-order ordinary differential equations with a small parameter at the highest derivatives and a point source. The question is addressed of reduction of this problem to a finite interval. A mesh, condensing in the boundary layer, is used for numerical solution of a system of singularly perturbed equations on a finite interval.