A robust adaptive method for a quasi-linear one-dimensional convection-diffusion problem

A quasi-linear conservative convection-diffusion two-point boundary value problem is considered. To solve it numerically, an upwind finite difference scheme is applied. The mesh used has a fixed number (N + 1) of nodes and is initially uniform, but its nodes are moved adaptively using a simple algorithm of de Boor based on equidistribution of the arc-length of the current computed piecewise linear solution. It is proved for the first time that a mesh exists that equidistributes the arc-length along the polygonal solution curve and that the corresponding computed solution is first-order accurate, uniformly in epsilon, where epsilon is the diffusion coefficient. In the case when the boundary value problem is linear, if N is sufficiently large independently of epsilon, it is shown that after O(ln(1/epsilon)/lnN) iterations of the algorithm, the piecewise linear interpolant of the computed solution achieves first-order accuracy in the L-infinity[0,1] norm uniformly in epsilon. Numerical experiments are presented that support our theoretical results.