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.     

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.      

A Finite Difference Scheme on a Priori Adapted Meshes for a Singularly Perturbed Parabolic Convection-Diffusion Equation

A boundary value problem is considered for a singularly perturbed parabolic convection-diffusion equation;we construct a finite difference scheme on a priori (se-quentially) adapted meshes and study its convergence.The scheme on a priori adapted meshes is constructed using a majorant function for the singular component of the discrete solution,which allows us to find a priori a subdomain where the computed solution requires a further improvement.This subdomain is defined by the perturbation parameterε,the step-size of a uniform mesh in x,and also by the required accuracy of the discrete solution and the prescribed number of refinement iterations K for im- proving the solution.To solve the discrete problems aimed at the improvement of the solution,we use uniform meshes on the subdomains.The error of the numerical so- lution depends weakly on the parameterε.The scheme converges almostε-uniformly, precisely,under the condition N~(-1)=o(ε~v),where N denotes the number of nodes in the spatial mesh,and the value v=v(K) can be chosen arbitrarily small for suitable K.     

Finite element methods on piecewise equidistant meshes for interior turning point problems

We consider linear second order singularly perturbed two-point boundary value problems with interior turning points. Piecewise linear Galerkin finite element methods are constructed on various piecewise equidistant meshes designed for such problems. These methods are proved to be convergent, uniformly in the singular perturbation parameter, in a weighted energy norm and the usualL ² norm. Supporting numerical results are presented.     

A class of singularly perturbed semilinear differential equations with interior layers

In this paper singularly perturbed semilinear differential equations with a discontinuous source term are examined. A numerical method is constructed for these problems which involves an appropriate piecewise-uniform mesh. The method is shown to be uniformly convergent with respect to the singular perturbation parameter. Numerical results are presented that validate the theoretical results.

     

Generating Layer-Adapted Meshes Using Mesh Partial Differential Equations

We present a new algorithm for generating layer-adapted meshes for the finite element solution of singularly perturbed problems based on mesh partial differential equations (MPDEs). The ultimate goal is to design meshes that are similar to the well-known Bakhvalov meshes, but can be used in more general settings: specifically two-dimensional problems for which the optimal mesh is not tensor-product in nature. Our focus is on the efficient implementation of these algorithms, and numerical verification of their properties in a variety of settings. The MPDE is a nonlinear problem, and the efficiency with which it can be solved depends adversely on the magnitude of the perturbation parameter and the number of mesh intervals. We resolve this by proposing a scheme based on $h$-refinement. We present fully working FEniCS codes [Alnaes et al., Arch. Numer. Softw., 3 (100) (2015)] that implement these methods, facilitating their extension to other problems and settings.     

A posteriori error analysis for defect correction method for two parameter singular perturbation problems

Defect correction method is used for two parameter singular perturbation problem on Bakhvalov-Shishkin mesh. Use of defect correction method on Bakhvalov-Shishkin mesh gives a second order convergence. A posteriori error estimate is obtained. The numerical examples are given to establish the second order convergence in practice.     

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.     

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.     

四阶奇异摄动边值问题在自适应网格上的一致收敛分析

we study a difference scheme for the fourth-order singular pertur-bation differential equation on the Bakhvalov-Shishkin grid by Green‘‘s function.The method is shown to be uniformly convergent with respect to the perturbation parameter,of order N^-2 in the maxmum norm on Bakhvalov-Shishkin meshes.Numerical results support our 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 Curved Macro-Element for Cracks and Corners

This note provides an extension of the quarter-point type offinite elements which are widley used in the solution of ellipticdifferential equations in regions with re-entrant corners. Itis shown how curved two or three-dimensional elements can bedevised that have the desired singular behaviour and permitarbitrary mesh refinement in the neighbourhood of the re-entrantcorner thus facilitating the attainment of high resolutionsapproximations. It is also shown how macro-elements can be usedto estimate the most advantageous size to use for the singularelements.     

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.     

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 posteriori error estimation for convection dominated problems on anisotropic meshes

A singularly perturbed convection–diffusion problem in two and three space dimensions is discretized using the streamline upwind Petrov Galerkin (SUPG) variant of the finite element method. The dominant convection frequently gives rise to solutions with layers; hence anisotropic finite elements can be applied advantageously. The main focus is on a posteriori energy norm error estimation that is robust in the perturbation parameter and with respect to the mesh anisotropy. A residual error estimator and a local problem error estimator are proposed and investigated. The analysis reveals that the upper error bound depends on the alignment of the anisotropies of the mesh and of the solution. Hence reliable error estimation is possible for suitable anisotropic meshes. The lower error bound depends on the problem data via a local mesh Peclet number. Thus efficient error estimation is achieved for small mesh Peclet numbers. Altogether, error estimation approaches for isotropic meshes are successfully extended to anisotropic elements. Several numerical experiments support the analysis. Copyright © 2003 John Wiley & Sons, Ltd.     

分层网格上奇异摄动问题的一致NIPG分析

本文采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在分层网格上至多相差一个关于摄动参数对数因子的拟最优阶的一致收敛性,即在能量范数度量下其误差估计为O((log~2(1/e))/N),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

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.     

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 layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011