Initial-Value Technique for Singularly Perturbed Boundary-Value Problems for Second-Order Ordinary Differential Equations Arising in Chemical Reactor Theory

An initial-value technique is presented for solving singularly perturbed two-point boundary-value problems for linear and semilinear second-order ordinary differential equations arising in chemical reactor theory. In this technique, the required approximate solution is obtained by combining solutions of two terminal-value problems and one initial-value problem which are obtained from the original boundary-value problem through asymptotic expansion procedures. Error estimates for approximate solutions are obtained. Numerical examples are presented to illustrate the present technique.     

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.     

Numerical experiments using Sukhanov's initial-value method for nonlinear two-point boundary-value problems,II

In a recent paper, Sukhanov derived a new method for transforming a nonlinear two-point boundary-value problem into an initial-value problem. Sukhanov's equations involve only the solution of ordinary differential equations and not partial differential equations. An earlier paper by the authors presented their interpretation of Sukhanov's method. An alternative method is presented in this paper. Numerical results are given.     

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

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

微分方程初值问题泰勒级数法的实现

提出了用微分变换来求解常微分方程初值问题的一个方法,该方法能通过迭代获得问题解析解的高阶Taylor级数的展开式,从而实现了高阶泰勒级数方法.     

Introduction to the theory of functional differential equations and their applications. Group approach

In this work, we give an introduction to the theory of nonlinear functional differential equations of pointwise type on a finite interval, semi-axis, or axis. This approach is based on the formalism using group peculiarities of such differential equations. For the main boundary-value problem and the Euler-Lagrange boundary-value problem, we consider the existence and uniqueness of the solution, the continuous dependence of the solution on boundary-value and initial-value conditions, and the “roughness” of functional differential equations in the considered boundary-value problems. For functional differential equations of pointwise type we also investigate the pointwise completeness of the space of solutions for given boundary-value conditions, give an estimate of the rank for the space of solutions, describe types of degeneration for the space of solutions, and establish conditions for the “smoothness” of the solution. We propose the method of regular extension of the class of ordinary differential equations in the class of functional differential equations of pointwise type. __________ Translated from Sovremennaya Matematika. Fundamental’nye Napravleniya (Contemporary Mathematics. Fundamental Directions), Vol. 8, Functional Differential Equations, 2004.     

Volteera型积分微分方程组边值问题的奇摄动

该文研究二阶积分微分方程组边值问题奇摄动，在适当的条件下，利用渐近分析方法和对角化技巧，还得解的存在性和给出解的渐近展开式与相应的余项估计．然后，应用这些结果到三阶常微分方程组边值问题的奇摄动，最后也得到解的一致有效的渐近展开式．     

Booster Method for Singularly-Perturbed One-Dimensional Reaction-Diffusion Neumann Problems

A numerical method for singularly-perturbed self-adjoint boundary-value problems for second-order ordinary differential equations subject to Neumann boundary conditions is proposed. In this method (booster method), an asymptotic approximation is incorporated into a finite-difference scheme to improve the numerical solution. Uniform error estimates are derived for this method when implemented in known difference schemes. Numerical examples are presented to illustrate the present method.     

Booster Method for Singularly-Perturbed One-Dimensional Convection-Diffusion Neumann Problems

An improved numerical method for singularly-perturbed two-point boundary-value problems for second-order ordinary differential equations subject to Neumann-type boundary conditions is proposed. In this method, an asymptotic approximation is incorporated into a finite-difference scheme to improve the numerical solution. Uniform error estimates are derived when implemented in known difference schemes. Numerical results are presented in support of the proposed method.     

On solvability of non-linear semi-periodic boundary-value problem for system of hyperbolic equations

We study a non-linear semi-periodic boundary-value problem for a system of hyperbolic equations with mixed derivative. At that, the semi-periodic boundary-value problem for a system of hyperbolic equations is reduced to an equivalent problem, consisting of a family of periodic boundary-value problems for ordinary differential equations and functional relation. When solving a family of periodic boundary-value problems of ordinary differential equations we use the method of parameterization. This approach allowed to establish sufficient conditions for the existence of an isolated solution of non-linear semi-periodic boundary-value problem for a system of hyperbolic equations.     

一类广义非线性强阻尼扰动发展方程的行波解

研究了一类非线性强阻尼广义扰动发展方程问题.它们在数学、力学、物理学等领域中广泛出现.首先,引入一个行波变换,把相应的偏微分方程问题转化为行波方程问题并求出原典型问题的精确解.再用小参数方法和引入伸长变量构造了问题的渐近解.最后, 用泛函分析的不动点理论证明了原非线性强阻尼广义扰动发展方程初值问题渐近行波解的存在性,并证明渐近解具有较高的精度和一致有效性.该文求得的渐近解是一个解析展开式, 所以它还可继续进行解析运算, 而单纯用数值模拟的方法是不行的.     

Numerical Applications of the Scaling Concept

We present a survey of recent developments in the applications of the scaling concept to numerical analysis. In addition, we report on some relevant topics not covered in existing surveys. Therefore, the present work updates and complements the existing surveys on the subject concerned.Applications of the scaling concept are useful in the numerical treatment of both ordinary and partial differential problems. Applications to boundary-value problems governed by ordinary differential equations are mainly related to their transformation into initial-value problems. Within this context, special emphasis is placed on systems of governing equations, eigenvalue, and free boundary-value problems. An error analysis for a truncated boundary formulation of the Blasius problem is also reported. As far as initial-value problems governed by ordinary differential equations are concerned, we discuss the development of adaptive mesh methods. Applications to partial differential problems considered herein are related to the construction of finite-difference schemes for conservations laws, the solution structure of the Riemann problem, rescaling schemes and adaptive schemes for blow-up problems.In writing this paper, our aim was to promote further and more important numerical applications of the scaling concept. Meanwhile, the pertinent bibliography is highlighted and is available on internet as the BIB file sc-gita.bib from the anonymous ftp area at the URL ftp://dipmat.unime.it/pub/papers/fazio/surveys.     

Symmetry analysis of initial-value problems

Symmetry analysis is a powerful tool that enables the user to construct exact solutions of a given differential equation in a fairly systematic way. For this reason, the Lie point symmetry groups of most well-known differential equations have been catalogued. It is widely believed that the set of symmetries of an initial-value problem (or boundary-value problem) is a subset of the set of symmetries of the differential equation. The current paper demonstrates that this is untrue; indeed, an initial-value problem may have no symmetries in common with the underlying differential equation. The paper also introduces a constructive method for obtaining symmetries of a particular class of initial-value problems.     

The solution of high-order nonlinear ordinary differential equations by Chebyshev Series

By the use of the Chebyshev series, a direct computational method for solving the higher order nonlinear differential equations has been developed in this paper. This method transforms the nonlinear differential equation into the matrix equation, which corresponds to a system of nonlinear algebraic equations with unknown Chebyshev coefficients, via Chebyshev collocation points. The solution of this system yields the Chebyshev coefficients of the solution function. An algorithm for this nonlinear system is also proposed in this paper. The method is valid for both initial-value and boundary-value problems. Several examples are presented to illustrate the accuracy and effectiveness 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.     

Invariant imbedding and the resolvent of Fredholm integral equations with semi-degenerate kernels

In many branches of astrophysics, physics, biology, and nuclear engineering, the underlying functional equation is a Fredholm integral equation. In this paper, it is shown that Fredholm integral equations with semi-degenerate kernels can be reduced to initial-value problems for systems of ordinary differential equations using an interesting formula for the Fredholm resolvent. Semi-degenerate kernels are encountered in many applications in the foregoing fields. This procedure facilitates the computational solution of the two-point boundary-value problem by both analog and digital computers.     

On Warner's algorithm for the solution of boundary-value problems for ordinary differential equations

The paper discusses the solution of boundary-value problems for ordinary differential equations by Warner's algorithm. This shooting algorithm requires that only the original system of differential equations is solved once in each iteration, while the initial conditions for a new iteration are evaluated from a matrix equation. Numerical analysis performed shows that the algorithm converges even for very bad starting values of the unknown initial conditions and that the number of iterations is small and weakly dependent on the starting point. Based on this algorithm, a general subroutine can be realized for the solution of a large class of boundary-value problems.     

An inexact Newton method for nonlinear two-point boundary-value problems

The method of quasilinearization for nonlinear two-point boundary-value problems is an application of Newton's method to a nonlinear differential operator equation. Since the linear boundary-value problem to be solved at each iteration must be discretized, it is natural to consider quasilinearization in the framework of an inexact Newton method. More importantly, each linear problem is only a local model of the nonlinear problem, and so it is inefficient to try to solve the linear problems to full accuracy. Conditions on size of the relative residual of the linear differential equation can then be specified to guarantee rapid local convergence to the solution of the nonlinear continuous problem. If initial-value techniques are used to solve the linear boundary-value problems, then an integration step selection scheme is proposed so that the residual criteria are satisfied by the approximate solutions. Numerical results are presented that demonstrate substantial computational savings by this type of economizing on the intermediate problems.This work was supported in part by DOE Contract DE-AS05-82-ER13016 and NSF Grant RII-89-17691 and was part of the author's doctoral thesis at Rice University. It is a pleasure to thank the author's thesis advisors, Professor R. A. Tapia and Professor J. E. Dennis, Jr.     

Asymptotic analysis of some classes of ordinary differential equations with large high-frequency terms

A complete asymptotic expansion is constructed for solutions of the Cauchy problem for nth order linear ordinary differential equations with rapidly oscillating coefficients, some of which may be proportional to ω ^n/2, where ω is oscillation frequency. A similar problem is solved for a class of systems of n linear first-order ordinary differential equations with coefficients of the same type. Attention is also given to some classes of first-order nonlinear equations with rapidly oscillating terms proportional to powers ω ^d . For such equations with d ∈ (1/2, 1], conditions are found that allow for the construction (and strict justification) of the leading asymptotic term and, in some cases, a complete asymptotic expansion of the solution of the Cauchy problem.     

Extrapolation methods in Lie groups

Summary. The numerical solution of differential equations on Lie groups by extrapolation methods is investigated. The main principles of extrapolation for ordinary differential equations are extended on the general case of differential equations in noncommutative Lie groups. An asymptotic expansion of the global error is given. A symmetric method is given and quadratic asymptotic expansion of the global error is proved. The theoretical results are verified by numerical experiments. Received September 27, 1999 / Revised version received February 14, 2000 / Published online April 5, 2001