Asymptotic initial-value method for second-order singular perturbation problems of reaction-diffusion type with discontinuous source term (Q946312)

The authors consider the boundary value problem (BVP) \[ -\varepsilon u''(x)+a(x)u(x)=-f(x),\quad x \in (0,d) \cup (d,1),\quad u(0)=p,\quad u(1)=q, \] where \(\varepsilon\) is a small parameter, \(a(x)\geq\alpha>0\) a smooth function, \(f(x)\) is a sufficient smooth function on \([0,\;d)\) and \((d,\;1]\) which has a discontinuity at \(d\). The solution of this singular BVP problem has, as a rule, two boundary and one interior layers. It is shown, that the solution can be approximated by the solution of the degenerated equation (with \(\varepsilon=0\)), i.e. \(u_0(x)=f(x)/a(x)\), and solutions of four auxiliary problems for equations of first order: two initial value problems on \((0,\;d)\) and two terminal value problems on \((d,\;1)\). For the numerical solution of the auxiliary problems the fitted mesh method is used. Two numerical examples are presented.