A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction-diffusion problems

A coupled system of two singularly perturbed linear reaction–diffusiontwo-point boundary value problems is examined. The leading termof each equation is multiplied by a small positive parameter,but these parameters may have different magnitudes. The solutionsto the system have boundary layers that overlap and interact.The structure of these layers is analysed, and this leads tothe construction of a piecewise-uniform mesh that is a variantof the usual Shishkin mesh. On this mesh central differencingis proved to be almost first-order accurate, uniformly in bothsmall parameters. Supporting numerical results are presentedfor a test problem.     

An efficient and uniformly convergent scheme for one-dimensional parabolic singularly perturbed semilinear systems of reaction-diffusion type

Numerical Algorithms - In this work we are interested in the numerical approximation of the solutions to 1D semilinear parabolic singularly perturbed systems of reaction-diffusion type, in the...     

A uniformly convergent continuous-discontinuous Galerkin method for singularly perturbed problems of convection-diffusion type

In this paper, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving singularly perturbed convection-diffusion problems. On Shishkin mesh with linear elements, a rate O(N^-1lnN) in an associated norm is established, where N is the number of elements. Numerical experiments complement the theoretical results. Moreover, a rate O(N^-2ln²N) in a discrete L^∞ norm, and O(N^-2) in L² norm, are observed numerically on the Shishkin mesh.     

A second order uniformly convergent numerical scheme for parameterized singularly perturbed delay differential problems

In this work, we propose a hybrid difference scheme for solving parameterized singularly perturbed delay differential problems. A unified error analysis framework for the proposed hybrid scheme is given that allows to conclude uniform convergence of \(\mathcal {O}(N^{-2}\ln ^{2} N)\) on Shishkin meshes and \(\mathcal {O}(N^{-2})\) on Bakhvalov meshes, where N is the number of mesh intervals in the domain. Numerical results are included to confirm the theoretical estimates.     

An almost fourth order uniformly convergent difference scheme for a semilinear singularly perturbed reaction-diffusion problem

Summary. This paper is concerned with a high order convergent discretization for the semilinear reaction-diffusion problem: , for , subject to , where . We assume that on , which guarantees uniqueness of a solution to the problem. Asymptotic properties of this solution are discussed. We consider a polynomial-based three-point difference scheme on a simple piecewise equidistant mesh of Shishkin type. Existence and local uniqueness of a solution to the scheme are analysed. We prove that the scheme is almost fourth order accurate in the discrete maximum norm, uniformly in the perturbation parameter . We present numerical results in support of this result. Received February 25, 1994     

A high order uniformly convergent alternating direction scheme for time dependent reaction–diffusion singularly perturbed problems

In this work we design and analyze an efficient numerical method to solve two dimensional initial-boundary value reaction–diffusion problems, for which the diffusion parameter can be very small with respect to the reaction term. The method is defined by combining the Peaceman and Rachford alternating direction method to discretize in time, together with a HODIE finite difference scheme constructed on a tailored mesh. We prove that the resulting scheme is ε-uniformly convergent of second order in time and of third order in spatial variables. Some numerical examples illustrate the efficiency of the method and the orders of uniform convergence proved theoretically. We also show that it is easy to avoid the well-known order reduction phenomenon, which is usually produced in the time integration process when the boundary conditions are time dependent. This research has been partially supported by the project MEC/FEDER MTM2004-01905 and the Diputación General de Aragón.     

A higher order uniformly convergent method for singularly perturbed parabolic turning point problems

In this article, we study numerical approximation for a class of singularly perturbed parabolic (SPP) convection-diffusion turning point problems. The considered SPP problem exhibits a parabolic boundary layer in the neighborhood of one of the sides of the domain. Some a priori bounds are given on the exact solution and its derivatives, which are necessary for the error analysis. A numerical scheme comprising of implicit finite difference method for time discretization on a uniform mesh and a hybrid scheme for spatial discretization on a generalized Shishkin mesh is proposed. Then Richardson extrapolation method is applied to increase the order of convergence in time direction. The resulting scheme has second-order convergence up to a logarithmic factor in space and second-order convergence in time. Numerical experiments are conducted to demonstrate the theoretical results and the comparative study is done with the existing schemes in literature to show better accuracy of the proposed schemes.     

An almost fourth order uniformly convergent domain decomposition method for a coupled system of singularly perturbed reaction-diffusion equations

We consider a system of M(≥2) singularly perturbed equations of reaction-diffusion type coupled through the reaction term. A high order Schwarz domain decomposition method is developed to solve the system numerically. The method splits the original domain into three overlapping subdomains. On two boundary layer subdomains we use a compact fourth order difference scheme on a uniform mesh while on the interior subdomain we use a hybrid scheme on a uniform mesh. We prove that the method is almost fourth order ε-uniformly convergent. Furthermore, we prove that when ε is small, one iteration is sufficient to get almost fourth order ε-uniform convergence. Numerical experiments are performed to support the theoretical results.     

An almost second order uniformly convergent scheme for a singularly perturbed initial value problem

In this paper a nonlinear singularly perturbed initial problem is considered. The behavior of the exact solution and its derivatives is analyzed, and this leads to the construction of a Shishkin-type mesh. On this mesh a hybrid difference scheme is proposed, which is a combination of the second order difference schemes on the fine mesh and the midpoint upwind scheme on the coarse mesh. It is proved that the scheme is almost second-order convergent, in the discrete maximum norm, independently of singular perturbation parameter. Numerical experiment supports these theoretical results.     

Uniformly convergent additive schemes for 2d singularly perturbed parabolic systems of reaction-diffusion type

Numerical Algorithms - In this work, we consider parabolic 2D singularly perturbed systems of reaction-diffusion type on a rectangle, in the simplest case that the diffusion parameter is the same...     

A second-order hybrid finite difference scheme for a system of singularly perturbed initial value problems

A system of coupled singularly perturbed initial value problems with two small parameters is considered. The leading term of each equation is multiplied by a small positive parameter, but these parameters may have different magnitudes. The solution of the system has boundary layers that overlap and interact. The structure of these layers is analyzed, and this leads to the construction of a piecewise-uniform mesh that is a variant of the usual Shishkin mesh. On this mesh a hybrid finite difference scheme is proved to be almost second-order accurate, uniformly in both small parameters. Numerical results supporting the theory are presented.     

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     

Parallel initial-value algorithms for singularly perturbed boundary-value problems

Initial-value methods for linear and semilinear singularly perturbed boundary-value problems are examined with a view to designing and implementing algorithms on parallel architectures. Practical experiments on a CRAY Y-MP 8/432 multiprocessor have been performed, showing the reliability and performance of several proposed parallel schemes.This work was supported by CNR, Rome, Italy (Progetto Finalizzato Sistemi Informatici e Calcolo Parallelo, Sottoprogetto 1).The authors wish to thank Dr. A. Papini, who carried out most of the computations reported in this work.     

A high order HODIE finite difference scheme for 1D parabolic singularly perturbed reaction-diffusion problems

This paper deals with the numerical approximation of the solution of 1D parabolic singularly perturbed problems of reaction-diffusion type. The numerical method combines the standard implicit Euler method on a uniform mesh to discretize in time and a HODIE compact fourth order finite difference scheme to discretize in space, which is defined on a priori special meshes condensing the grid points in the boundary layer regions. The method is uniformly convergent having first order in time and almost fourth order in space. The analysis of the uniform convergence is made in two steps, splitting the contribution to the error from the time and the space discretization. Although this idea has been previously used to prove the uniform convergence for parabolic singularly perturbed problems, here the proof is based on a new study of the asymptotic behavior of the exact solution of the semidiscrete problems obtained after the time discretization by using the Euler method. Some numerical results are given corroborating in practice the theoretical results.     

Initial-value methods for second-order singularly perturbed boundary-value problems

In this paper, we present a numerical method for solving linear and nonlinear second-order singularly perturbed boundary-value-problems. For linear problems, the method comes from the well-known WKB method. The required approximate solution is obtained by solving the reduced problem and one or two suitable initial-value problems, directly deduced from the given problem. For nonlinear problems, the quasilinearization method is applied. Numerical results are given showing the accuracy and feasibility of the proposed method.This work was supported in part by the Consiglio Nazionale delle Ricerche (Contract No. 86.02108.01 and Progetto Finalizzatto Sistemi Informatia e Calcolo Paralello, Sottoprogetto 1), and in part by the Ministero della Pubblica Istruzione, Rome, Italy.     

New initial-value method for singularly perturbed boundary-value problems

An initial-value method is given for second-order singularly perturbed boundary-value problems with a boundary layer at one endpoint. The idea is to replace the original two-point boundary value problem by two suitable initial-value problems. The method is very easy to use and to implement. Nontrivial text problems are used to show the feasibility of the given method, its versatility, and its performance in solving linear and nonlinear singularly perturbed problems.This work was supported in part by the Consiglio Nazionale delle Ricerche, Contract No. 86.02108.01, and in part by the Ministero della Pubblica Istruzione.     

Layer-adapted meshes for a linear system of coupled singularly perturbed reaction-diffusion problems

Niall Madden ^{We consider a system of 2 one-dimensional singularly perturbedreaction–diffusion equations coupled at the zero-orderterm. The second derivative of each equation is multiplied bya distinct small parameter. We show how to decompose the solutionto the problem into regular and layer parts. Properties of thediscretized operator are established using discrete Green'sfunctions. We prove that a central difference scheme on certainlayer-adapted meshes converges independently of the perturbationparameters. Supporting numerical examples confirm our theoreticalresults.}     

Finite-element approximation for singularly perturbed nonlinear elliptic boundary-value problems

Numerical solution of a two-dimensional nonlinear singularly perturbed elliptic partial differential equation ∈ Δu = f(x, u), 0 < x, y < 1, with Dirichlet boundary condition is discussed here. The modified Newton method of third-order convergence is employed to linearize the nonlinear problem in place of the standard Newton method. The finite-element method is used to find the solution of the nonlinear differential equation. Numerical results are provided to demonstrate the usefulness of the method.     

A coupled system of singularly perturbed parabolic reaction-diffusion equations

In this paper systems with an arbitrary number of singularly perturbed parabolic reaction-diffusion equations are examined. A numerical method is constructed for these systems which involves an appropriate layer-adapted piecewise-uniform mesh. The numerical approximations generated from this method are shown to be uniformly convergent with respect to the singular perturbation parameters. Numerical experiments supporting the theoretical results are given.     

A uniformly convergent method for a singularly perturbed semilinear reaction--diffusion problem with multiple solutions

This paper considers a simple central difference scheme for a singularly perturbed semilinear reaction--diffusion problem, which may have multiple solutions. Asymptotic properties of solutions to this problem are discussed and analyzed. To compute accurate approximations to these solutions, we consider a piecewise equidistant mesh of Shishkin type, which contains points. On such a mesh, we prove existence of a solution to the discretization and show that it is accurate of order , in the discrete maximum norm, where the constant factor in this error estimate is independent of the perturbation parameter and . Numerical results are presented that verify this rate of convergence.