An approach to singular perturbation problems insoluble by asymptotic methods

Singular perturbation problems not amenable to solution by asymptotic methods require special treatment, such as the method of Carrier and Pearson. Rather than devising special methods for these problems, this paper suggests that there may be a uniform way to solve singular perturbation problems, which may or may not succumb to asymptotic methods. A potential mechanism for doing this is the author's boundary-value technique, a nonasymptotic method, which previously has only been applied to singular perturbation problems that lend themselves to asymptotic techniques. Two problems, claimed by Carrier and Pearson to be insoluble by asymptotic methods, are solved by the boundary-value method.     

Second-order numerical method for domain optimization problems

This paper is concerned with a second-order numerical method for shape optimization problems. The first variation and the second variation of the objective functional are derived. These variations are discretized by introducing a set of boundary-value problems in order to derive the second-order numerical method. The boundary-value problems are solved by the conventional finite-element method.The authors would like to express their thanks to Mr. T. Masanao, who was an undergraduate student, for his cooperation and comments. They also thank Professor Y. Sakawa of Osaka University for his encouragement.A part of this paper was presented at the IFIP Conference on Control of Boundaries and Stabilization, Clermont-Ferrand, France, 1988.     

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.     

A steepest-descent method for optimization of mechanical systems

Steepest-descent optimal control techniques have been used extensively for dynamic systems in one independent variable and with a full set of initial conditions. This paper presents an extension of the steepest-descent technique to mechanical design problems that are described by boundary-value problems with one or more independent variables. The method is illustrated by solving finite-dimensional problems, problems with distribution of design over one space dimension, and problems with distribution of design over two space dimensions.     

Asymptotic Initial-Value Method for Singularly-Perturbed Boundary-Value Problems for Second-Order Ordinary Differential Equations

A computational method is presented to solve a class of nonturning-point singularly-perturbed two-point boundary-value problems for second-order ordinary differential equations with a small parameter multiplying the highest derivative, subject to Dirichlet-type boundary conditions. In this method, first we construct a zeroth order asymptotic expansion for the solution of the given boundary-value problem. Then, this problem is integrated to get an equivalent initial-value problem for first-order ordinary differential equations. This initial-value problem is solved by either a classical method or a fitted operator method after approximating some of the terms in the differential equations by using the zeroth order asymptotic expansion. This method is effective and easy to implement. An error estimate is derived for the numerical solution. Examples are given to illustrate the method.     

Fast solution of elliptic control problems

Elliptic control problems with a quadratic cost functional require the solution of a system of two elliptic boundary-value problems. We propose a fast iterative process for the numerical solution of this problem. The method can be applied to very special problems (for example, Poisson equation for a rectangle) as well as to general equations (arbitrary dimensions, general region). Also, nonlinear problems can be treated. The work required is proportional to the work taken by the numerical solution of a single elliptic equation.     

A collocation/quadrature-based Sturm–Liouville problem solver

We present a computational method for solving a class of boundary-value problems in Sturm–Liouville form. The algorithms are based on global polynomial collocation methods and produce discrete representations of the eigenfunctions. Error control is performed by evaluating the eigenvalue problem residuals generated when the eigenfunctions are interpolated to a finer discretization grid; eigenfunctions that produce residuals exceeding an infinity-norm bound are discarded. Because the computational approach involves the generation of quadrature weights and arrays for discrete differentiation operations, our computational methods provide a convenient framework for solving boundary-value problems by eigenfunction expansion and other projection methods.     

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.     

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 model trust-region modification of Newton's method for nonlinear two-point boundary-value problems

The method of quasilinearization for nonlinear two-point boundary-value problems is Newton's method for a nonlinear differential operator equation. A model trust-region approach to globalizing the quasilinearization algorithm is presented. A double-dogleg implementation yields a globally convergent algorithm that is robust in solving difficult 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 J. E. Dennis, Jr., and Professor R. A. Tapia.     

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.     

A spectral collocation method based on integrated Chebyshev polynomials for two-dimensional biharmonic boundary-value problems

This paper reports a new spectral collocation method for numerically solving two-dimensional biharmonic boundary-value problems. The construction of the Chebyshev approximations is based on integration rather than conventional differentiation. This use of integration allows: (i) the imposition of the governing equation at the whole set of grid points including the boundary points and (ii) the straightforward implementation of multiple boundary conditions. The performance of the proposed method is investigated by considering several biharmonic problems of first and second kinds; more accurate results and higher convergence rates are achieved than with conventional differential methods.     

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).     

Parametric Cubic Spline Approach to the Solution of a System of Second-Order Boundary-Value Problems

We use parametric cubic spline functions to develop a numerical method for computing approximations to the solution of a system of second-order boundary-value problems associated with obstacle, unilateral, and contact problems. We show that the present method gives approximations which are better than those produced by other collocation, finite-difference, and spline methods. A numerical example is given to illustrate the applicability and efficiency of the new method.     

Parametric Quintic-Spline Approach to the Solution of a System of Fourth-Order Boundary-Value Problems

In this paper, we use parametric quintic splines to derive some consistency relations which are then used to develop a numerical method for computing the solution of a system of fourth-order boundary-value problems associated with obstacle, unilateral, and contact problems. It is known that a class of variational inequalities related to contact problems in elastostatics can be characterized by a sequence of variational inequations, which are solved using some numerical method. Numerical evidence is presented to show the applicability and superiority of the new method over other collocation, finite difference, and spline methods.     

Explicit solutions of boundary-value problems for (2 + 1)-dimensional integrable systems

Two nonlinear integrable models with two space variables and one time variable, the Kadomtsev-Petviashvili equation and the two-dimensional Toda chain, are studied as well-posed boundary-value problems that can be solved by the inverse scattering method. It is shown that there exists a multitude of integrable boundary-value problems and, for these problems, various curves can be chosen as boundary contours; besides, the problems in question become problems with moving boundaries. A method for deriving explicit solutions of integrable boundary-value problems is described and its efficiency is illustrated by several examples. This allows us to interpret the integrability phenomenon of the boundary condition in the traditional sense, namely as a condition for the availability of wide classes of solutions that can be written in terms of well-known functions.     

Conjugate residual methods for almost symmetric linear systems

This paper concerns the use of conjugate residual methods for the solution of nonsymmetric linear systems arising in applications to differential equations. We focus on an application derived from a seismic inverse problem. The linear system is a small perturbation to a symmetric positive-definite system, the nonsymmetries arising from discretization errors in the solution of certain boundary-value problems. We state and prove a new error bound for a class of generalized conjugate residual methods; we show that, in some cases, the perturbed symmetric problem can be solved with an error bound similar to the one for the conjugate residual method applied to the symmetric problem. We also discuss several applications for special distributions of eigenvalues.This work was supported in part by the National Science Foundation, Grants DMS-84-03148 and DCR-81-16779, and by the Office of Naval Research, Contract N00014-85-K-0725.     

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.     

Shooting Methods for Numerical Solution of Stochastic Boundary-Value Problems

Abstract

In the present investigation, numerical methods are developed for approximate solution of stochastic boundary-value problems. In particular, shooting methods are examined for numerically solving systems of Stratonovich boundary-value problems. It is proved that these methods accurately approximate the solutions of stochastic boundary-value problems. An error analysis of these methods is performed. Computational simulations are given.