Convergence of linear multistep methods for two-parameter singular perturbation problems

1. IntroductionIt is well known that the IVPs of ODEs in singular perturbation form may be consideredas a special class of stab problems. But it is regrettable that they can't be satisfactorily covered by B-theory (of. [2-8]) because of their very special structure. In the recede ten yeaxslmany importallt and illteresting results on the convergence of linear multistep methods,Runge-Kutta methods, Rosenbrock methods and general linear methods applied to oneparameter singular perturbation pr…     

Convergence of parallel multistep hybrid methods for singular perturbation problems

Parallel multistep hybrid methods (PHMs) can be implemented in parallel with two processors, accordingly have almost the same computational speed per integration step as BDF methods of the same order with the same stepsize. But PHMs have better stability properties than BDF methods of the same order for stiff differential equations. In the present paper, we give some results on error analysis of A(α)-stable PHMs for the initial value problems of ordinary differential equations in singular perturbation form. Our convergence results are similar to those of linear multistep methods (such as BDF methods), i.e. the convergence orders are equal to their classical convergence orders, and no order reduction occurs. Some numerical examples also confirm our results.     

Numerical solution of singular perturbation problems by a terminal boundary-value technique

In this paper, we present a new approach for numerically solving linear singularly perturbed two-point boundary-value problems in ordinary differential equations with a boundary layer on the left end of the interval. The original problem is divided into outer and inner region problems. A terminal boundary condition in implicit form is introduced. Then, the outer region problem is solved as a two-point boundary-value problem (TPBVP), and an explicit terminal boundary condition is obtained. In turn, the inner region problem is modified and solved as a TPBVP using the explicit terminal boundary condition. The proposed method is iterative on the terminal point of the inner region. Three numerical examples have been solved to demonstrate the applicability of the method.     

High order finite volume methods for singular perturbation problems

In this paper we establish a high order finite volume method for the fourth order singular perturbation problems.In conjunction with the optimal meshes,the numerical solutions resulting from the method have optimal convergence order.Numerical experiments are presented to verify our theoretical estimates.     

Initial-value technique for a class of nonlinear singular perturbation problems

An initial-value technique, which is simple to use and easy to implement, is presented for a class of nonlinear, singularly perturbed two-point boundary-value problems with a boundary layer on the left end of the underlying interval. It is distinguished by the following fact: The original second-order problem is replaced by an asymptotically equivalent first-order problem and is solved as an initial-value problem. Numerical experience with several examples is described.     

An exponentially fitted special second-order finite difference method for solving singular perturbation problems

An exponentially fitted special second-order finite difference method is presented for solving singularly perturbed two-point boundary value problems with the boundary layer at one end (left or right) point. A fitting factor is introduced in a tri-diagonal finite difference scheme and is obtained from the theory of singular perturbations. Thomas Algorithm is used to solve the system and its stability is investigated. To demonstrate the applicability of the method, we have solved several linear and non-linear problems. From the results, it is observed that the present method approximates the exact solution very well.     

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.     

Convergence of one-leg methods for singular perturbation problems with delays

This paper is concerned with the error behaviour of one-leg methods applied to some classes of one-parameter multiple stiff singularly perturbed problems with delays. We obtain convergence results of A-stable one-leg methods with linear interpolation procedure. Numerical experiments further confirm our theoretical analysis.     

Numerical integration of a class of singular perturbation problems

In this paper, we discuss an approximate method for the numerical integration of a class of linear, singularly perturbed two-point boundary-value problems in ordinary differential equations with a boundary layer on the left end of the underlying interval. This method requires a minimum of problem preparation and can be implemented easily on a computer. We replace the original singular perturbation problem by an approximate first-order differential equation with a small deviating argument. Then, we use the trapezoidal formula to obtain the three-term recurrence relationship. Discrete invariant imbedding algorithm is used to solve a tridiagonal algebraic system. The stability of this algorithm is investigated. The proposed method is iterative on the deviating argument. Several numerical experiments have been included to demonstrate the efficiency of the method.The authors wish to express their sincere thanks to Dr. S. M. Roberts for his comments and valuable suggestions.     

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.     

A nonasymptotic method for general linear singular perturbation problems

The nonasymptotic method developed in Ref. 1 has been extended for solving general linear singularly perturbed two-point boundary-value problems. Firstly, we discuss problems with a right-hand boundary layer. Secondly, we discuss problems with an interior layer. Finally, we discuss problems with two boundary layers. Numerical experience with the method for some model problems is also reported to confirm the theoretical analysis.     

Variable step size initial value algorithm for singular perturbation problems using locally exact integration

We consider quasilinear singular perturbation problems of the form εy^″+p(x)y^′+q(x,y)=h(x),x[0,1];y(0)=,y(1)=β with a boundary layer at one end point. The original problem is reduced to an asymptotically equivalent linear first order initial-value problem (IVP). Then, a variable step size initial value algorithm is applied to solve this (IVP). The algorithm is based on the locally exact integration of quadratic linearized problem coefficients on a non-uniform mesh. Two term-recurrence relation with controlled step size is obtained. Several problems are solved to demonstrate the applicability and efficiency of the algorithm. It is observed that the present method approximates the exact solution very well.     

The layer-resolving transformation and mesh generation for quasilinear singular perturbation problems

The relationship is analyzed between layer-resolving transformations and mesh-generating functions for numerical solution of singularly perturbed boundary-value problems. The analysis is carried out for one-dimensional quasilinear problems without turning points, which are discretized by first-order finite-difference schemes. It is proved that if a general layer-resolving function is used to generate the discretization mesh, then the numerical solution converges uniformly in the perturbation parameter.     

A finite differences MATLAB code for the numerical solution of second order singular perturbation problems

We show the main features of the MATLAB code HOFiD_UP for solving second order singular perturbation problems. The code is based on high order finite differences, in particular on the generalized upwind method. Within its simplicity, it uses order variation and continuation for solving any difficult nonlinear scalar problem. Several numerical tests on linear and nonlinear problems are considered. The best performances are reported on problems with perturbation parameters near the machine precision, where most of the codes for two-point BVPs fail.     

A local monotonicity formula for an inhomogenous singular perturbation problem and applications

In this paper we prove a local monotonicity formula for solutions to an inhomogeneous singularly perturbed diffusion problem of interest in combustion. This type of monotonicity formula has proved to be very useful for the study of the regularity of limits u of solutions of the singular perturbation problem and of ∂{u > 0}, in the global homogeneous case. As a consequence of this formula we prove that u has an asymptotic development at every point in ∂{u > 0} where there is a nonhorizontal tangent ball. These kind of developments have been essential for the proof of the regularity of ∂{u > 0} for Bernoulli and Stefan free boundary problems. We also present applications of our results to the study of the regularity of ∂{u > 0} in the stationary case including, in particular, its regularity in the case of energy minimizers. We present as well a regularity result for traveling waves of a combustion model that relies on our monotonicity formula and its consequences.The fact that our results hold for the inhomogeneous problem allows a very wide applicability. Indeed, they may be applied to problems with nonlocal diffusion and/or transport. The research of the authors was partially supported by Fundación Antorchas Project 13900-5, Universidad de Buenos Aires grant X052, ANPCyT PICT No 03-13719, CONICET PIP 5478. The authors are members of CONICET.     

A variational approach to singular quasilinear elliptic problems with sign changing nonlinearities

We obtain positive solutions of singular p-Laplacian problems with sign changing nonlinearities using variational methods.     

Geometric mesh FDM for self-adjoint singular perturbation boundary value problems

A numerical method based on finite difference method with variable mesh is given for second order singularly perturbed self-adjoint two point boundary value problems. The original problem is reduced to its normal form and the reduced problem is solved by FDM taking variable mesh(geometric mesh). The maximum absolute errors max_i|y(x_i)-y_i|, for different values of parameter , number of points N, and the mesh ratio r, for three examples have been given in tables to support the efficiency of the method.     

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.     

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.     

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

利用匹配渐近展开法,讨论了一类四阶非线性方程的具有两个边界层的奇摄动边值问题．引进伸长变量,根据边界条件与匹配原则,在一定的可解性条件下,给出了外部解和左右边界层附近的内层解,得到了该问题的二阶渐近解,并举例说明了这类非线性问题渐近解的存在性．