关于耦合型对流-扩散方程组的奇异摄动边值问题

本文考虑下述耦合型对流-扩散方程组的奇异摄动边值问题:本文提出两种方法:一种是初值化解法,用这种方法,原始问题转化成一系列没有扰动的一阶微分方程或方程组的初值问题,从而得到一个渐近展开式;第二种是边值化解法,用这种方法,原始问题转化成一组没有边界层现象的边值问题,从而可以得到精确解和使用经典的数值方法去得到具有关于摄动参数ε一致的高精度数值解.     

On Spurious Solutions of Singular Perturbation Problems

A study is made of several nonlinear boundary-value problems of singular perturbation type for which a straightforward application of boundary-layer theory leads to spurious solutions. It is shown that these problems can be treated successfully by a slight modification of the method of matched asymptotic expansions. The analysis leads to several novel features which are not present in routine singular perturbation problems.     

Helly’s Principle and Its Application to an Infinite-Horizon Optimal Control Problem

In this paper, a numerical method is presented to solve singularly-perturbed two-point boundary-value problems for second-order ordinary differential equations with a discontinuous source term. First, an asymptotic expansion approximation of the solution of the boundary-value problem is constructed using the basic ideas of the well-known WKB perturbation method. Then, some initial-value problems and terminal-value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial-value problems and terminal-value problems are singularly-perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples are provided to illustrate the method.     

Asymptotic and numerical analysis of singular perturbation problems: A survey

This paper is intended to be a brief survey of the asymptotic and numerical analysis of singular perturbation problems. The purpose is to find out what problems are treated and what numerical/asymptotic methods are employed, with an eye toward the goal of developing general methods to solve such problems. A summary of the results of some recent methods is presented, and this leads to conclusions and recommendations about what methods to use on singular perturbation problems. Finally, some areas of current research are indicated. A bibliography of about 130 items is provided.     

Numerical treatment of a mathematical model arising from a model of neuronal variability

In this paper, we describe a numerical approach based on finite difference method to solve a mathematical model arising from a model of neuronal variability. The mathematical modelling of the determination of the expected time for generation of action potentials in nerve cells by random synaptic inputs in dendrites includes a general boundary-value problem for singularly perturbed differential-difference equation with small shifts. In the numerical treatment for such type of boundary-value problems, first we use Taylor approximation to tackle the terms containing small shifts which converts it to a boundary-value problem for singularly perturbed differential equation. A rigorous analysis is carried out to obtain priori estimates on the solution of the problem and its derivatives up to third order. Then a parameter uniform difference scheme is constructed to solve the boundary-value problem so obtained. A parameter uniform error estimate for the numerical scheme so constructed is established. Though the convergence of the difference scheme is almost linear but its beauty is that it converges independently of the singular perturbation parameter, i.e., the numerical scheme converges for each value of the singular perturbation parameter (however small it may be but remains positive). Several test examples are solved to demonstrate the efficiency of the numerical scheme presented in the paper and to show the effect of the small shift on the solution behavior.     

Completeness theorem for singular differential pencils

A theorem completeness theorem of special vector functions induced by the products of the so-called Weyl solutions of a fourth-order differential equation and by their derivatives on the semiaxis is presented. We prove that such nonlinear combinations of Weyl solutions and their derivatives constitute a linear subspace of decreasing (at infinity) solutions of a linear singular differential system of Kamke type. We construct and study the Green function of the corresponding singular boundary-value problems on the semiaxis for operator pencils defining differential systems of Kamke type. The required completeness theorem is proved by using the analytic and asymptotic properties of the Green function, operator spectral theory methods, and analytic function theory.     

求解奇异摄动边值问题的精细积分法

提出了一种求解一端有边界层的奇异摄动边值问题的精细方法．首先将求解区域均匀离散,由状态参量在相邻节点间的精细积分关系式确定一组代数方程,并将其写成矩阵形式．代入边界条件后,该代数方程组的系数矩阵可化为块三对角形式,针对这一特性,给出了一种高效递推消元方法．由于在离散过程中,精细积分关系式不会引入离散误差,故所提出的方法具有极高的精度．数值算例充分证明了所提出方法的有效性．     

Asymptotic solution of a singularly perturbed set of functional-differential equations of Riccati type encountered in the optimal control theory

The boundary-value problem for the set of functional-differential equations with partial derivatives of Riccati type, associated with a singularly perturbed linear-quadratic optimal control problem with delay in state, is considered. The expression for a solution of the problem, which transforms it to the explicit singular perturbation form, is proposed. An asymptotic solution of this problem is constructed. Received August 7, 1997     

B-Spline Collocation Method for Nonlinear Singularly-Perturbed Two-Point Boundary-Value Problems

A B-spline collocation method is presented for nonlinear singularly-perturbed boundary-value problems with mixed boundary conditions. The quasilinearization technique is used to linearize the original nonlinear singular perturbation problem into a sequence of linear singular perturbation problems. The B-spline collocation method on piecewise uniform mesh is derived for the linear case and is used to solve each linear singular perturbation problem obtained through quasilinearization. The fitted mesh technique is employed to generate a piecewise uniform mesh, condensed in the neighborhood of the boundary layers. The convergence analysis is given and the method is shown to have second-order uniform convergence. The stability of the B-spline collocation system is discussed. Numerical experiments are conducted to demonstrate the efficiency of the method.     

A Method for Incorporating Transcendentally Small Terms into the Method of Matched Asymptotic Expansions

The essential ideas behind a method for incorporating exponentially small terms into the method of matched asymptotic expansions are demonstrated using an Ackerberg–O'Malley resonance problem and a spurious solutions problem of Carrier and Pearson. One begins with the application of the standard method of matched asymptotic expansions to obtain at least the leading terms in outer and inner (Poincaré-type) expansions; some, although not all, matching can be carried out at this stage. This is followed by the introduction of supplementary expansions whose gauge functions are transcendentally small compared to those in the standard expansions. Analysis of terms in these expansions allows the matching to be completed. Furthermore, the method allows for the inclusion of globally valid transcendentally small contributions to the asymptotic solution; it is well known that such terms may be numerically significant.     

椭圆-抛物偏微分方程奇异摄动问题的差分解法

本文讨论了含有小参数在高阶导数项的椭圆型方程奇异摄动问题的差分解法.当ε=0时椭圆型方程退化为抛物型方程.作者根据此问题解的边界层性质,构造了特殊的差分格式:研究了它的收敛性和解的渐近性态.最后给出一个数值例题.     

Slowly Varying Phase Planes and Boundary-Layer Theory

A method is presented that combines phase-plane techniques with the ideas of multiple scale and matched asymptotic expansions to explain the behavior of solutions to second-order, nonlinear, nonautonomous, singular boundary-value problems. It is shown that if one is willing to give up the detailed information provided by a procedure such as matched asymptotic expansions, then complete qualitative information can be obtained by the much simpler method given here. (“Complete” here means that the method provides a way of categorizing all possible solutions of such problems.) In addition, the similarities and differences between the present method and that of Melnikov, which has been useful in the study of dynamical systems, are noted.     

Asymptotic expansions using blow-up

The method of matched asymptotic expansions and geometric singular perturbation theory are the most common and successful approaches to singular perturbation problems. In this work we establish a connection between the two approaches in the context of the simple fold problem. Using the blow-up technique [5], [12] and the tools of geometric singular perturbation theory we derive asymptotic expansions of slow manifolds continued beyond the fold point. Our analysis explains the structure of the expansion and gives an algorithm for computing its coefficients.     

矩阵奇异值分解问题重分析的摄动法

本文提出了一般实矩阵奇异值分解问题重分析的摄动法.这是一种简捷、高效的快速重分析方法,对于提高各种需要反复进行矩阵奇异值分解的迭代分析问题的计算效率具有较重要的实用价值.文中导出了奇异值和左、右奇异向量的直到二阶摄动量的渐近估计算式.文末指出了将这种振动分析方法直接推广到一般复矩阵情况的途径.     

Numerical solution of singular control problems using multiple shooting techniques

Algorithms for calculating the junction points between optimal nonsingular and singular subarcs of singular control problems are developed. The algorithms consist in formulating appropriate initialvalue and boundary-value problems; the boundary-value problems are solved with the method of multiple shooting. Two examples are detailed to illustrate the proposed numerical methods.The author would like to thank Professor Dr. R. Bulirsch, who stimulated and encouraged this work, which is part of the author's dissertation.     

Asymptotic expansions using blow-up

The method of matched asymptotic expansions and geometric singular perturbation theory are the most common and successful approaches to singular perturbation problems. In this work we establish a connection between the two approaches in the context of the simple fold problem. Using the blow-up technique [5], [12] and the tools of geometric singular perturbation theory we derive asymptotic expansions of slow manifolds continued beyond the fold point. Our analysis explains the structure of the expansion and gives an algorithm for computing its coefficients.*Research supported by the Austrian Science Foundation under grant Y 42-MAT.Received: February 1, 2001; revised: November 22, 2002     

Uniform Global Convergence of a Hybrid Scheme for Singularly Perturbed Reaction–Diffusion Systems

We consider a system of coupled singularly perturbed reaction–diffusion two-point boundary-value problems. A hybrid difference scheme on a piecewise-uniform Shishkin mesh is constructed for solving this system, which generates better approximations to the exact solution than the classical central difference scheme. Moreover, we prove that the method is third order uniformly convergent in the maximum norm when the singular perturbation parameter is small. Numerical experiments are conducted to validate the theoretical results.     

差分方程奇异摄动问题的渐近解

本文对含小参数的差分方程奇异摄动问题构造了一种新的渐近方法.     

Modified deflation algorithm for the solution of singular problems. II. Nonlinear multipoint boundary value problems

Various kinds of iterative methods have been proposed for the solution of nonlinear multipoint boundary-value problems MPBVP's. However, it is necessary for these methods that the adjusting matrix, which corresponds to the Jacobian of nonlinear equations, is nonsingular at the solution. In this paper an algorithm for the singular solution of nonlinear MPBVP's, which is an extension of the modified deflation algorithm for the singular root of nonlinear algebraic equations developed by the author is presented. According to the present method, the singular solution can ultimately be reduced to the usual simple solution and both convergency and accuracy can greatly be improved. The effectiveness of the present method is shown by solving two illustrative examples.     

Innovative solution of a 2D elastic transmission problem¶

This article is concerned with a boundary-field equation approach to a class of boundary value problems exterior to a thin domain. A prototype of this kind of problems is the interaction problem with a thin elastic structure. We are interested in the asymptotic behavior of the solution when the thickness of the elastic structure approaches to zero. In particular, formal asymptotic expansions will be developed, and their rigorous justification will be considered. As will be seen, the construction of these formal expansions hinges on the solutions of a sequence of exterior Dirichlet problems, which can be treated by employing boundary element methods. On the other hand, the justification of the corresponding formal procedure requires an independence on the thickness of the thin domain for the constant in the Korn inequality. It is shown that in spite of the reduction of the dimensionality of the domain under consideration, this class of problems are, in general, not singular perturbation problems, because of appropriate interface conditions.