Collocation for Singular Perturbation Problems I: First Order Systems with Constant Coefficients
- 1 June 1983
- journal article
- Published by Society for Industrial & Applied Mathematics (SIAM) in SIAM Journal on Numerical Analysis
- Vol. 20 (3), 537-557
- https://doi.org/10.1137/0720035
Abstract
The application of collocation methods for the numerical solution of singularly perturbed ordinary differential equations is investigated. Collocation at Gauss, Radau and Lobatto points is considered, for both initial and boundary value problems for first order systems with constant coefficients. Particular attention is paid to symmetric schemes for boundary value problems; these problems may have boundary layers at both interval ends. .br Our analysis shows that certain collocation schemes, in particular those based on Gauss or Lobatto points, do perform very well on such problems, provided that a fine mesh with steps proportional to the layers'' width is used in the layers only, and a coarse mesh, just fine enough to resolve the solution of the reduced problem, is used in between. Ways to construct appropriate layer meshes are proposed. Of all methods considered, the Lobatto schemes appear to be the most promising class of methods, as they essentially retain their usual superconvergence power for the smooth, reduced solution, whereas Gauss-Legendre schemes do not. .br We also investigate the conditioning of the linear systems of equations arizing in the discretization of the boundary value problem. For a row equilibrated version of the discretized system we obtain a pleasantly small bound on the maximum norm condition number, which indicates that these systems can be solved safely by Gaussian elimination with scaled partial pivoting.Keywords
This publication has 14 references indexed in Scilit:
- Analytical and Numerical Methods for Nonlinear Singular Singularly-Perturbed Initial Value ProblemsSIAM Journal on Applied Mathematics, 1980
- Solving boundary-value problems with a spline-collocation codeJournal of Computational Physics, 1980
- SOLVEBLOK: A Package for Solving Almost Block Diagonal Linear SystemsACM Transactions on Mathematical Software, 1980
- On the Boundary Value Problem for Systems of Ordinary Differential Equations with a Singularity of the Second KindSIAM Journal on Mathematical Analysis, 1980
- Stability Criteria for Implicit Runge–Kutta MethodsSIAM Journal on Numerical Analysis, 1979
- A collocation solver for mixed order systems of boundary value problemsMathematics of Computation, 1979
- An Adaptive Finite Difference Solver for Nonlinear Two-Point Boundary Problems with Mild Boundary LayersSIAM Journal on Numerical Analysis, 1977
- The numerical solution of boundary value problems for stiff differential equationsMathematics of Computation, 1977
- Collocation at Gaussian PointsSIAM Journal on Numerical Analysis, 1973
- Implicit Runge-Kutta processesMathematics of Computation, 1964