Numerical treatment of singularly perturbed two point boundary value problems

We propose a method 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 underlying interval. This is a practical method and can be easily implemented on a computer. The original problem is divided into inner and outer region differential equation systems. The reduced problem is solved to obtain the terminal boundary condition. Then, a new inner region problem is created and solved as a two point boundary value problem (TPBVP). In turn, the outer region problem is also solved as a TPBVP. Both these TPBVPs are efficiently treated by employing a slightly modified classical finite difference scheme coupled with discrete invariant imbedding algorithm to obtain the numerical solutions. The stability of some recurrence relations involved in the algorithm is investigated. The proposed method is iterative on the terminal point. Some numerical examples are included, and the computational results are compared with exact solutions. It is observed that the accuracy predicted can always be achieved with very little computational effort.     

Numerical treatment of singularly perturbed two point boundary value problems exhibiting boundary layers

A boundary value method for solving a class of nonlinear singularly perturbed two point boundary value problems with a boundary layer at one end is proposed. Using singular perturbation analysis the method consists of solving two problems; namely, a reduced problem and a boundary layer correction problem. We use Pade’ approximation to obtain the solution of the latter problem and to satisfy the condition at infinity. Numerical examples will be given to illustrate the method.     

Numerical solution of singularly perturbed fifth order two point boundary value problem

We will consider Adomain decomposition method and the homotopy method to solve a fifth order singularly perturbed BVP arising in viscoelastic flows. The success and pitfalls of the methods will be investigated. Numerical testing will be provided to show the efficiency of the methods proposed. Comparison with the work of others will also be done.     

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.     

Global convergence method for singularly perturbed boundary value problems

We give uniformly convergent splines difference scheme for singularly perturbed boundary value problems(1)$-\epsilon u^{″}+p(x)u^{\prime }+q(x)u=f(x),u(a)={\alpha }_0,u(b)={\alpha }_1,$by using splines fitted with delta sequence due to the very stiff nature of the problem under consideration. We prove the $\mathrm{O}(\min (h^2,{\epsilon }^2))$ order of uniform convergence with respect to small parameter $\epsilon$ at nodes on uniform mesh and $\mathrm{O}(\min (h,\epsilon ))$ order of uniform global convergence with respect to the approximate solution given by $S(x)={\sum }_{i=1}^N{S}_{{\Delta }_i}(x)H(x_i-x)$ where H is the Heaviside function, which is the approximation for the closed form of the exact solution.     

奇摄动非线性边值问题

The singularly perturbed nonlinear boundary value problems are considered. Using the stretched variable and the method of boundary layer correction,the formal asymptotic expansion of solution is obtained. And then the uniform validity of solution is proved by using the differential inequalities.     

Exponential B-spline collocation method for self-adjoint singularly perturbed boundary value problems

We present an exponential B-spline collocation method for self-adjoint singularly perturbed boundary value problem. The convergence analysis is given and the method is shown to have second order uniform convergence. Numerical experiments are conducted to demonstrate the efficiency of the method.     

Numerical solution of second-order nonlinear singularly perturbed boundary-value problems by initial-value methods

Nonlinear singularly perturbed boundary-value problems are considered, with one or two boundary layers but no turning points. The theory of differential inequalities is used to obtain a numerical procedure for quasilinear and semilinear problems. The required solution is approximated by combining the solutions of suitable auxiliary initial-value problems easily deduced from the given problem. From the numerical results, the method seems accurate and solutions to problems with extremely thin layers can be obtained at reasonable cost.This work was supported by CNR, Rome, Italy (Progetto Finalizzato Sistemi Informatici e Calcolo Parallelo, Sottoprogetto 1).     

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.     

The Lie-group shooting method for solving nonlinear singularly perturbed boundary value problems

A new computational method for solving the second-order nonlinear singularly perturbed boundary value problems (SPBVPs) is provided in this paper. In order to overcome a highly singular behavior very near to the boundary as being not easy to treat by numerical method, we adopt a coordinate transformation from an x-domain to a t-domain via a rescaling technique, which can reduce the singularity within the boundary layer. Then, we construct a Lie-group shooting method (LGSM) to search a missing initial condition through the finding of a suitable value of a parameter r ∈ [0, 1]. Moreover, we can derive a closed-form formula to express the initial condition in terms of r, which can be determined properly by an accurate matching to the right-boundary condition. Numerical examples are examined, showing that the present approach is highly efficient and accurate.     

Piecewise shooting reproducing kernel method for linear singularly perturbed boundary value problems

In this letter, a new numerical method is proposed for solving second order linear singularly perturbed boundary value problems with left layers. Firstly a piecewise reproducing kernel method is proposed for second order linear singularly perturbed initial value problems. By combining the method and the shooting method, an effective numerical method is then proposed for solving second order linear singularly perturbed boundary value problems. Two numerical examples are used to show the effectiveness of the present method.     

The Galerkin method for singularly perturbed boundary value problems on adaptive grids

Voronezh. Translated from Sibirskii Matematicheskii Zhurnal, Vol. 31, No. 5, pp. 138–148, September–October, 1990.     

Reproducing kernel method for singularly perturbed turning point problems having twin boundary layers

A numerical method is proposed for solving singularly perturbed turning point problems exhibiting twin boundary layers based on the reproducing kernel method (RKM). The original problem is reduced to two boundary layers problems and a regular domain problem. The regular domain problem is solved by using the RKM. Two boundary layers problems are treated by combining the method of stretching variable and the RKM. The boundary conditions at transition points are obtained by using the continuity of the approximate solution and its first derivatives at these points. Two numerical examples are provided to illustrate the effectiveness of the present method. The results compared with other methods show that the present method can provide very accurate approximate solutions.     

Boundary integral methods for singularly perturbed boundary value problems

In this paper we consider boundary integral methods appliedto boundary value problems for the positive definite Helmholtz-typeproblem –U + ²U = 0 in a bounded or unbounded domain,with the parameter real and possibly large. Applications arisein the implementation of space–time boundary integralmethods for the heat equation, where is proportional to 1/(t),and t is the time step. The corresponding layer potentials arisingfrom this problem depend nonlinearly on the parameter and havekernels which become highly peaked as , causing standard discretizationschemes to fail. We propose a new collocation method with arobust convergence rate as . Numerical experiments on a modelproblem verify the theoretical results.     

一类具有无穷边界值的二次奇摄动边值问题

研究一类具有无穷边界值的二次奇摄动Robin边值问题解的存在性与解的渐进行为,重点关注边界值的奇异程度对解的边界层行为的影响;同时将所得的结果与Chang及Howes的结果(带正常边界值)进行比较.研究表明:(1)当边界值大小为O(1/)时,得到的边界层大小为O( ln ),这比Chang及Howes带正常边界值的情形提高了O(ln )量级;(2)增大边界值的奇性至O(1/ r),这里r >1,边界层大小的量级不变,依然为O( ln );(3)若要使得边界层大小为O(1),则边界值的大小需为O(e?1/).最后给出一个算例验证得到的结果.     

A novel method for a class of parameterized singularly perturbed boundary value problems

In this paper a novel approach is presented for solving parameterized singularly perturbed two-point boundary value problems with a boundary layer. By the boundary layer correction technique, the original problem is converted into two non-singularly perturbed problems which can be solved using traditional numerical methods, such as Runge–Kutta methods. Several non-linear problems are solved to demonstrate the applicability of the method. Numerical experiments indicate the high accuracy and the efficiency of the new method.     

Convergence of a Bakhvalov grid adaptation method for singularly perturbed boundary value problems

In this paper, a Galerkin finite element method for non-self-adjoint boundary value problems on Bakhvalov grids is considered. Using the Galerkin projectionmethod, the convergence of a sequence of computational grids with an unknown boundary of the boundary layer is proved. Some numerical examples are presented.     

Parameter-uniform numerical methods for singularly perturbed mixed boundary value problems using grid equidistribution

In this paper, we present the analysis of an upwind scheme for obtaining the solution of a convection-diffusion two-point boundary value problem with Robin boundary conditions. The solution is obtained on a suitable nonuniform mesh which is formed by equidistributing the arc-length monitor function. It is shown that the discrete solution obtained by the upwind scheme converges uniformly with respect to the perturbation parameter. Numerical results are presented that demonstrate the sharpness of the theoretical estimates.