Parameter-uniform numerical methods for some linear and nonlinear singularly perturbed convection diffusion boundary turning point problems

Both linear and nonlinear singularly perturbed two point boundary value problems are examined in this paper. In both cases, the problems have a boundary turning point and are of convection-diffusion type. Parameter-uniform numerical methods composed of monotone finite difference operators and piecewise-uniform Shishkin meshes, are constructed and analyzed for both the linear and the nonlinear class of problems. Numerical results are presented to illustrate the theoretical parameter-uniform error bounds established.     

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.     

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.     

GUARANTEED BOUNDS FOR THE SOLUTION OF THE WAVE EQUATION

Abstract

Numerical solution of the wave equation in the form of close lower and upper bounds provides a secure a posteriori error estimate that can be used for efficient accuracy control. The method considered in this paper uses some monotone properties of the differential operator in the wave equation to construct bounds for the solution in the form of trigonometric polynomials of x. Aspects of the numerical implementation, the accuracy of the computed bounds and some numerical examples are discussed.     

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.     

Scaled discrete derivatives of singularly perturbed elliptic problems

Numerical approximations to the solution of a singularly perturbed elliptic convection–diffusion problem in two space dimensions are generated using a monotone finite difference operator on a tensor product of piecewise‐uniform Shishkin meshes. The bilinear interpolants of these numerical approximations are parameter‐uniformly convergent to the solution of the continuous problem, in the pointwise maximum norm. In this article, discrete approximations to the first derivatives of the solution are shown to be globally first‐order (up to logarithmic factors) uniformly convergent, when the errors are scaled within the analytical layers of the continuous problem. Numerical results are presented to illustrate the theoretical error bounds established in an appropriated weighted C¹–norm. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 225–252, 2015     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

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.

     

A posteriori error control in low-order finite element discretisations of incompressible stationary flow problems

Computable a posteriori error bounds and related adaptive mesh-refining algorithms are provided for the numerical treatment of monotone stationary flow problems with a quite general class of conforming and nonconforming finite element methods. A refined residual-based error estimate generalises the works of Verfürth; Dari, Duran and Padra; Bao and Barrett. As a consequence, reliable and efficient averaging estimates can be established on unstructured grids. The symmetric formulation of the incompressible flow problem models certain nonNewtonian flow problems and the Stokes problem with mixed boundary conditions. A Helmholtz decomposition avoids any regularity or saturation assumption in the mathematical error analysis. Numerical experiments for the partly nonconforming method analysed by Kouhia and Stenberg indicate efficiency of related adaptive mesh-refining algorithms.

     

多边形网格上扩散方程新的单调格式

 本文在星形多边形网格上, 构造了扩散方程新的单调有限体积格式.该格式与现有的基于非线性两点流的单调格式的主要区别是, 在网格边的法向流离散模板中包含当前边上的点, 在推导离散法向流的表达式时采用了定义于当前边上的辅助未知量, 这样既可适应网格几何大变形, 同时又兼顾了当前网格边上物理量的变化. 在光滑解情形证明了离散法向流的相容性.对于具有强各向异性、非均匀张量扩散系数的扩散方程, 证明了新格式是单调的, 即格式可以保持解析解的正性. 数值结果表明在扭曲网格上, 所构造的格式是局部守恒和保正的, 对光滑解有高于一阶的精度, 并且, 针对非平衡辐射限流扩散问题, 数值结果验证了新格式在计算效率和守恒精度上优于九点格式.     

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.

     

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 modified accelerated monotone iterative method for finite difference reaction-diffusion-convection equations

This paper is concerned with monotone algorithms for the finite difference solutions of a class of nonlinear reaction-diffusion-convection equations with nonlinear boundary conditions. A modified accelerated monotone iterative method is presented to solve the finite difference systems for both the time-dependent problem and its corresponding steady-state problem. This method leads to a simple and yet efficient linear iterative algorithm. It yields two sequences of iterations that converge monotonically from above and below, respectively, to a unique solution of the system. The monotone property of the iterations gives concurrently improving upper and lower bounds for the solution. It is shown that the rate of convergence for the sum of the two sequences is quadratic. Under an additional requirement, quadratic convergence is attained for one of these two sequences. In contrast with the existing accelerated monotone iterative methods, our new method avoids computing local maxima in the construction of these sequences. An application using a model problem gives numerical results that illustrate the effectiveness of the proposed method.     

Accelerated monotone iterative methods for finite difference equations of reaction-diffusion

This paper is concerned with numerical methods for a finite difference system of reaction-diffusion-convection equation under nonlinear boundary condition. Various monotone iterative methods are presented, and each of these methods leads to an existence-comparison theorem as well as a computational algorithm for numerical solutions. The monotone property of the iterations gives improved upper and lower bounds of the solution in each iteration, and the rate of convergence of the iterations is either quadratic or nearly quadratic depending on the property of the nonlinear function. Application is given to a model problem from chemical engineering, and some numerical results, including a test problem with known analytical solution, are presented to illustrate the various rates of convergence of the iterations. Received November 2, 1995 / Revised version received February 10, 1997     

A block monotone domain decomposition algorithm for a semilinear convection–diffusion problem

This paper deals with discrete monotone iterative algorithms for solving a nonlinear singularly perturbed convection–diffusion problem. A block monotone domain decomposition algorithm based on a Schwarz alternating method and on block iterative scheme is constructed. This monotone algorithm solves only linear discrete systems at each iterative step of the iterative process and converges monotonically to the exact solution of the nonlinear problem. The rate of convergence of the block monotone domain decomposition algorithm is estimated. Numerical experiments are presented.     

Approximation of Derivative for a Singularly Perturbed Second-Order ODE of Robin Type with Discontinuous Convection Coefficient and Source Term

In this paper, a singularly perturbed Robin type boundary value problem for second-order ordinary differential equation with discontinuous convection coefficient and source term is considered. A robust-layer-resolving numerical method is proposed. An e-uniform global error estimate for the numerical solution and also to the numerical derivative are established. Numerical results are presented, which are in agreement with the theoretical predictions.AMS subject classifications: 65L10, CR G1.7     

Error and stability of monotone method for numerical solutions of fourth-order semilinear elliptic boundary value problems

This paper is concerned with the error and stability analysis of the monotone method for numerical solutions of fourth-order semilinear elliptic boundary value problems. A comparison result among the various monotone sequences is given. The global error is analyzed, and some sufficient conditions are formulated to guarantee a geometric rate of convergence. The stability of the monotone method is proved. Some numerical results are presented.     

Hamilton-Jacobi方程的单调有限元格式

A monotone finite element scheme is obtained by applying the finite element method to the viscosity equation of the Hamilton-Jacobi equation on unstructured meshes. Under some constraints, we show that this scheme is monotone and its numerical solution converges to the viscosity solution of the Hamilton-Jacobi equa-tion. Numerical examples test the stability and the convergence of this scheme.     

Monotonous enclosures for the Thomas-Fermi equation in the isolated neutral atom case

Please address any correspondence to: C. Grossmann, Department of Mathematics, Kuwait University, PO Box 5969, Safat 13060, Kuwait A numerical method for the generation of enclosures for thesolution of the Thomas–Fermi equation on a semi-infiniteinterval is proposed. The method is based on the monotone discretizationprinciple and on available global bounds for the solution. Theconvergence of the new method on refined and extended gridsis investigated where available bounds are used to increasethe local step size with increasing arguments.     

Uniform quadratic convergence of monotone iterates for nonlinear singularly perturbed parabolic problems

This paper deals with a monotone iterative method for solving nonlinear singularly perturbed parabolic problems. Monotone sequences, based on the method of upper and lower solutions, are constructed for a nonlinear difference scheme which approximates the nonlinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. The monotone sequences possess quadratic convergence rate. An analysis of uniform convergence of the monotone iterative method to the solutions of the nonlinear difference scheme and to the continuous problem is given. Numerical experiments are presented.