Approximation of a system of singularly perturbed reaction-diffusion parabolic equations in a rectangle

The Dirichlet problem for a system of singularly perturbed reaction-diffusion parabolic equations in a rectangle is considered. The higher order derivatives of the equations are multiplied by a perturbation parameter ?², where ? takes arbitrary values in the interval (0, 1]. When ? vanishes, the system of parabolic equations degenerates into a system of ordinary differential equations with respect to t. When ? tends to zero, a parabolic boundary layer with a characteristic width ? appears in a neighborhood of the boundary. Using the condensing grid technique and the classical finite difference approximations of the boundary value problem, a special difference scheme is constructed that converges ?-uniformly at a rate of O(N ^?2ln² N + N ₀ ^?1 , where \(N = \mathop {\min }\limits_s N_s \), N ^s + 1 and N ₀ + 1 are the numbers of mesh points on the axes x _s and t, respectively.     

Computer difference scheme for a singularly perturbed convection-diffusion equation

The Dirichlet problem for a singularly perturbed ordinary differential convection-diffusion equation with a perturbation parameter ? (that takes arbitrary values from the half-open interval (0, 1]) is considered. For this problem, an approach to the construction of a numerical method based on a standard difference scheme on uniform meshes is developed in the case when the data of the grid problem include perturbations and additional perturbations are introduced in the course of the computations on a computer. In the absence of perturbations, the standard difference scheme converges at an \(\mathcal{O}\) (δ _st ) rate, where δ _st = (? + N ^?1)^?1 N ^?1 and N + 1 is the number of grid nodes; the scheme is not ?-uniformly well conditioned or stable to perturbations of the data. Even if the convergence of the standard scheme is theoretically proved, the actual accuracy of the computed solution in the presence of perturbations degrades with decreasing ? down to its complete loss for small ? (namely, for ? = \(\mathcal{O}\) (δ^?2max _i,j |δa _i ^j | + δ^?1 max _{i, j} |δb _i ^j |), where δ = δ _st and δa _i ^j , δb _i ^j are the perturbations in the coefficients multiplying the second and first derivatives). For the boundary value problem, we construct a computer difference scheme, i.e., a computing system that consists of a standard scheme on a uniform mesh in the presence of controlled perturbations in the grid problem data and a hypothetical computer with controlled computer perturbations. The conditions on admissible perturbations in the grid problem data and on admissible computer perturbations are obtained under which the computer difference scheme converges in the maximum norm for ? ∈ (0, 1] at the same rate as the standard scheme in the absence of perturbations.     

Grid approximation of singularly perturbed parabolic convection-diffusion equations subject to a piecewise smooth initial condition

A boundary value problem for a singularly perturbed parabolic convection-diffusion equation on an interval is considered. The higher order derivative in the equation is multiplied by a parameter ? that can take arbitrary values in the half-open interval (0, 1]. The first derivative of the initial function has a discontinuity of the first kind at the point x ₀. For small values of ?, a boundary layer with the typical width of ? appears in a neighborhood of the part of the boundary through which the convective flow leaves the domain; in a neighborhood of the characteristic of the reduced equation outgoing from the point (x ₀, 0), a transient (moving in time) layer with the typical width of ?^1/2 appears. Using the method of special grids that condense in a neighborhood of the boundary layer and the method of additive separation of the singularity of the transient layer, special difference schemes are designed that make it possible to approximate the solution of the boundary value problem ?-uniformly on the entire set $\bar G$ , approximate the diffusion flow (i.e., the product ?(?/?x)u(x, t)) on the set $\bar G^ * = \bar G\backslash \{ (x_0 ,0)\} $ , and approximate the derivative (?/?x)u(x, t) on the same set outside the m-neighborhood of the boundary layer. The approximation of the derivatives ?²(?²/?x ²)u(x, t) and (?/?t)u(x, t) on the set $\bar G^ * $ is also examined.     

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.     

Analysis of a DG method for singularly perturbed convection-diffusion problems

In this article, we studied a discontinuous Galerkin finite element method for convection-diffusion-reaction problems with singular perturbation. Our approach is highly flexible by allowing the use of discontinuous approximating function on polytopal mesh without imposing extra conditions on the convection coefficient. A priori error estimate is devised in a suitable energy norm on general polytopal mesh. Numerical examples are provided.     

Numerical method for a singularly perturbed convection-diffusion problem with delay

This paper deals with the singularly perturbed boundary value problem for a linear second-order delay differential equation. For the numerical solution of this problem, we use an exponentially fitted difference scheme on a uniform mesh which is accomplished by the method of integral identities with the use of exponential basis functions and interpolating quadrature rules with weight and remainder term in integral form. It is shown that one gets first order convergence in the discrete maximum norm, independently of the perturbation parameter. Numerical results are presented which illustrate the theoretical results.     

A posteriori error estimate for finite volume approximations to singularly perturbed nonlinear convection-diffusion equations

Summary. This paper is devoted to the study of a posteriori and a priori error estimates for the scalar nonlinear convection diffusion equation . The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the -norm in the situation, where the diffusion parameter is smaller or comparable to the mesh size. Numerical experiments underline the theoretical results. Received February 25, 1999 / Revised version received July 6, 1999 / Published online August 2, 2000     

Approximation of systems of singularly perturbed elliptic reaction-diffusion equations with two parameters

In a rectangle, the Dirichlet problem for a system of two singularly perturbed elliptic reaction-diffusion equations is considered. The higher order derivatives of the ith equation are multiplied by the perturbation parameter ? _i ² (i = 1, 2). The parameters ?_i take arbitrary values in the half-open interval (0, 1]. When the vector parameter ? = (?₁, ?₂) vanishes, the system of elliptic equations degenerates into a system of algebraic equations. When the components ?₁ and (or) ?₂ tend to zero, a double boundary layer with the characteristic width ?₁ and ?₂ appears in the vicinity of the boundary. Using the grid refinement technique and the classical finite difference approximations of the boundary value problem, special difference schemes that converge ?-uniformly at the rate of O(N ^?2ln² N) are constructed, where N = min N _s, N _s + 1 is the number of mesh points on the axis x _s.     

Effective solution of a singularly perturbed system of nonlinear differential equations

The authors consider a boundary-value problem for a singularly-perturbed system of three nonlinear ordinary differential equations arising in modeling the interaction of physical fields in a semiconducting diode. They suggest a higher-accuracy method for solution of this problem based on the combination of the Newton operator method and the continuation method with respect to a parameter. In this case, the Frechét derivative, the Green function of the corresponding differential equation, and the initial approximation are found in explicit analytical form. The numerical implementation confirms the higher effectivity and the superexponential rate of convergence of the proposed method. __________ Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 36, Suzdal Conference-2004, Part 2, 2005.     

A steplike contrast structure in a singularly perturbed system of elliptic equations

A singularly perturbed boundary value problem for a system of elliptic equations in a two-dimensional region is considered. The asymptotics and existence of a solution with an internal transition layer are studied. The asymptotics is justified by the method of differential inequalities.     

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.     

Conditioning of finite difference schemes for a singularly perturbed convection-diffusion parabolic equation

In the case of the boundary value problem for a singularly perturbed convection-diffusion parabolic equation, conditioning of an ε-uniformly convergent finite difference scheme on a piecewise uniform grid is examined. Conditioning of a finite difference scheme on a uniform grid is also examined provided that this scheme is convergent. For the condition number of the scheme on a piecewise uniform grid, an ε-uniform bound O(δ ₁ ^?2 lnδ ₁ ^?1 + δ ₀ ^?1 ) is obtained, where δ₁ and δ₀ are the error components due to the approximation of the derivatives with respect to x and t, respectively. Thus, this scheme is ε-uniformly well-conditioned. For the condition number of the scheme on a uniform grid, we have the estimate O(ε^?1δ ₁ ^?2 + δ ₀ ^?1 ); this scheme is not ε-uniformly well-conditioned. In the case of the difference scheme on a uniform grid, there is an additional error due to perturbations of the grid solution; this error grows unboundedly as ε → 0, which reduces the accuracy of the grid solution (the number of correct significant digits in the grid solution is reduced). The condition numbers of the matrices of the schemes under examination are the same; both have an order of O(ε^?1δ ₁ ^?2 + δ ₀ ^?1 ). Neither the matrix of the ε-uniformly convergent scheme nor the matrix of the scheme on a uniform grid is ε-uniformly well-conditioned.     

Application of Sinc-Galerkin method to singularly perturbed parabolic convection-diffusion problems

We develop a numerical algorithm for solving singularly perturbed one-dimensional parabolic convection-diffusion problems. The method comprises a standard finite difference to discretize in temporal direction and Sinc-Galerkin method in spatial direction. The convergence analysis and stability of proposed method are discussed in details, it is justifying that the approximate solution converges to the exact solution at an exponential rate. we know that the conventional methods for these problems suffer due to decreasing of perturbation parameter, but the Sinc method handel such difficulty as singularity. This scheme applied on some test examples, the numerical results illustrate the efficiency of the method and confirm the theoretical behavior of the rates of convergence.     

Dynamics of a singularly perturbed system of two differential equations with delay

Theoretical and Mathematical Physics - We study a two-dimensional singularly perturbed system with delay, which is a simplification of models used in laser physics. We analyze several cases of a...     

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.     

Numerical approximations to a singularly perturbed convection-diffusion problem with a discontinuous initial condition

A singularly perturbed parabolic problem of convection-diffusion type with a discontinuous initial condition is examined. An analytic function is identified which matches the discontinuity in the initial condition and also satisfies the homogenous parabolic differential equation associated with the problem. The difference between this analytical function and the solution of the parabolic problem is approximated numerically, using an upwind finite difference operator combined with an appropriate layer-adapted mesh. The numerical method is shown to be parameter-uniform. Numerical results are presented to illustrate the theoretical error bounds established in the paper.

     

Inverse singularly perturbed problems of the convection-diffusion type in quadrangular curvilinear domains

We have constructed an algorithm for the asymptotic approximation of the solutions of inverse singularly perturbed boundary-value problems of the convection-diffusion type with unknown diffusion coefficient, depending on the coordinates of a quadrangular curvilinear domain of filtration. The case of sufficient smoothness and consistency of the overdetermination, initial, and boundary conditions is considered. Unlike the construction of an algorithm for the solution of similar problems in doubly connected domains, here, in the corresponding relations, there appear corrections taking into account the influence of “lateral sources of pollution.” With the help of this algorithm, we have carried out a computer experiment, the results of which confirm the well-known fact of “strong sensitivity” of the model to assignment of the overdetermination condition. In particular, we have revealed the specific character of influence of this condition on the required diffusion coefficient depending on the filtration velocity.     

Uniform mesh approximation to nonsmooth solutions of a singularly perturbed convection-diffusion equation in a rectangle

Here nonsmooth solutions of a differential equation are treated as solutions for which the compatibility conditions are not required to hold at the corner points of the domain and hence corner singularities can occur. In the present paper, we drop the compatibility conditions at three of the four vertices of a rectangle. At the remaining vertex, from which a characteristic (inclined) of the reduced equation issues, we impose compatibility conditions providing the C ^3,λ -smoothness of the desired solution in a neighborhood of that vertex as well as additional conditions leading to the smoothness of solutions of the reduced equation occurring in the regular component of the solution of the considered problem. Under our assumptions and for a sufficient smoothness of the coefficients of the equation and its right-hand side, we show that the classical five-point upwind approximation on a Shishkin piecewise uniform mesh preserves the accuracy specific for the smooth case; i.e., the mesh solution uniformly (with respect to a small parameter) converges in the L _∞ ^h -norm to the exact solution at the rate O(N ⁻¹ ln² N), where N is the number of mesh nodes in each of the coordinate directions.     

Asymptotic behaviour of three-dimensional singularly perturbed convection-diffusion problems with discontinuous data

We consider three singularly perturbed convection-diffusion problems defined in three-dimensional domains: (i) a parabolic problem −?(uxx+uyy)+ut+v₁ux+v₂uy=0 in an octant, (ii) an elliptic problem −?(uxx+uyy+uzz)+v₁ux+v₂uy+v₃uz=0 in an octant and (iii) the same elliptic problem in a half-space. We consider for all of these problems discontinuous boundary conditions at certain regions of the boundaries of the domains. For each problem, an asymptotic approximation of the solution is obtained from an integral representation when the singular parameter ?→⁰⁺. The solution is approximated by a product of two error functions, and this approximation characterizes the effect of the discontinuities on the small ?− behaviour of the solution and its derivatives in the boundary layers or the internal layers.