A uniform numerical method for quasilinear singular perturbation problems without turning points (Q1115126)

Consider the singularly perturbed boundary value problem \((1)\quad Tu:=- \epsilon u''-b(x,u)u'+c(x,u)=0,\) \(x\in I=[0,1]\), \(Bu:=(u(0),u(1))=(U_ 0,U_ 1),\) where the functions b, c and the numbers \(U_ 0\), \(U_ 1\) are given, \(\epsilon\) is a small parameter, \(0<\epsilon \leq \epsilon_ 0\). Suppose that \(b,c\in C^ 2\) (I\(\times R)\), \(b(x,u)\geq \beta >0\), \(b_ x(x,u)+c_ u(x,u)\geq \gamma\), \(\gamma\leq 0\), \(x\in I\), \(\beta^ 2+4\epsilon \gamma >0\). Then (T,B) is inverse monotone and there exists a unique solution \(u\in C^ 4(I)\) of the problem (1). It has a boundary layer near the origin. To approximate u uniformly in \(\epsilon\), the problem (1) is first transformed by introducing a special new independent variable. Then the Engquist-Osher type finite-difference scheme is applied. Stability uniform in \(\epsilon\) is proved in the discrete \(L^ 1\)-norm. First order consistency uniform in \(\epsilon\) and first order convergence uniform in \(\epsilon\) in the discrete \(L^ 1\)-norm are proved. Numerical examples illustrate the pointwise covergence uniform in \(\epsilon\).