Grid equidistribution for reaction-diffusion problems in one dimension (Q2583229)

This paper deals with an adaptive numerical method for the linear reaction-diffusion two-point boundary value problem: \(-\varepsilon^2u''(x)+b(x)u(x)=f(x)\), \(x\in(0,1)\), \(u(0) = u(1) =0\), where \(b,f\in [0,1]\) and \(0<\beta <b(x)\leq \overline b\) on \([0,1]\). Some facts about the solution \(u(x)\) and stability bounds for the computed solution \(u_N\) (on an arbitrary mesh) are listed. These bounds are used to provide strong heuristic evidence for the choice of monitor function to be equidistributed in an adaptive algorithm. Some numerical results are presented.