Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)



Numerical analysis of boundary-layer problems in ordinary differential equations

Author: W. D. Murphy
Journal: Math. Comp. 21 (1967), 583-596
MSC: Primary 65.61
MathSciNet review: 0225496
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: We categorize some of the finite-difference methods that can be used to treat the initial-value problem for the boundary-layer differential equation

$\displaystyle (1){\text{ }}\mu y' = f(y,x);y(0) = {y^0}.$

These methods take the form

$\displaystyle (2){\text{ }}\sum\limits_{i = 0}^k {{\alpha _i}{Y_{n + i}} = {h^{... ...amma }}} \sum\limits_{i = 0}^k {{\beta _i}f} ({Y_{n + i}},{x_{n + i}}) + {R_n},$

where $ {\alpha _\nu }$ and $ {\beta _\nu }(\nu = 0,1, \cdots ,k)$ denote real constants which do not depend upon $ n,{R_n}$ is the round-off error, $ \mu = {h^r},0 < \gamma < 1$, and $ h$ is the mesh size. We define a new kind of stability called $ \mu $stability and prove that under certain conditions $ \mu $stability implies convergence of the difference method. We investigate $ \mu $stability and the optimal methods which it allows, i.e., methods of maximum accuracy.

The idea of relating $ \mu $ to $ h$ allows us to study the nature of the difference equation for very small $ \mu $. We can, however, look at this in another way. Given a differential equation in the form of Eq. (1) we ask how can we choose $ h$ so that the associated difference equation will give an accurate approximation. If $ \mu $ is sufficiently small, choose $ h$ by the formula $ h = {\mu ^{1/\gamma }}$ where $ 0 < \gamma < 1$.

References [Enhancements On Off] (What's this?)

  • [1] Germund Dahlquist, Convergence and stability in the numerical integration of ordinary differential equations, Math. Scand. 4 (1956), 33–53. MR 0080998
  • [2] Peter Henrici, Discrete variable methods in ordinary differential equations, John Wiley & Sons, Inc., New York-London, 1962. MR 0135729
  • [3] T. E. Hull and W. A. J. Luxemburg, Numerical methods and existence theorems for ordinary differential equations, Numer. Math. 2 (1960), 30–41. MR 0114017
  • [4] T. E. Hull and A. C. R. Newbery, Error bounds for a family of three-point integration procedures, J. Soc. Indust. Appl. Math. 7 (1959), 402–412. MR 0136079
  • [5] T. E. Hull and A. C. R. Newbery, Integration procedures which minimize propagated errors, J. Soc. Indust. Appl. Math. 9 (1961), 31–47. MR 0120770
  • [6] W. D. Murphy, "Numerical analysis of boundary layer problems," AEC Research and Development Report NYO-1480-63, New York University.
  • [7] A. B. Vasiĺeva, "Asymptotic behavior of solutions to certain problems involving nonlinear differential equations containing a small parameter multiplying the highest derivatives," Russian Math. Surveys, v. 18, 1963, no. 3, pp. 13-84.

Similar Articles

Retrieve articles in Mathematics of Computation with MSC: 65.61

Retrieve articles in all journals with MSC: 65.61

Additional Information

Article copyright: © Copyright 1967 American Mathematical Society