Mesh selection for a nearly singular boundary value problem
asymptotic error analysisboundary-value-problemscentral difference schemeconvergenceequidistribution principlemesh adaptationnumerical testsShishkin meshsingular perturbation
Linear boundary value problems for ordinary differential equations (34B05) Singular perturbations for ordinary differential equations (34E15) Numerical solution of boundary value problems involving ordinary differential equations (65L10) Finite difference and finite volume methods for ordinary differential equations (65L12) Stability and convergence of numerical methods for ordinary differential equations (65L20) Mesh generation, refinement, and adaptive methods for ordinary differential equations (65L50) Error bounds for numerical methods for ordinary differential equations (65L70)
The authors investigate the numerical solution of the model equation \(u_{xx}=\varepsilon^{-2}e^{-x/\varepsilon }\) (\(x\in (0,1)\)) with \(u(0)=1\) and \(u(1)=\alpha =e^{-1/\varepsilon }\) and several slightly more general problems when \(\varepsilon << 1\) using the standard central difference scheme on nonuniform grids. In particular, it is shown the error behavior in two limiting cases: (i) the total mesh point number \(N\) is fixed when the regularization parameter \(\varepsilon \to 0\); (ii) \(\varepsilon \) is fixed when \(N\to \infty \). Using a formal analysis, the authors show that a generalized version of a special piecewise uniform mesh and an adaptive grid based on the equidistribution principle share some common features. The optimal meshes give rates of convergence bounded by \(|\log (\varepsilon)|\) as \(\varepsilon\to 0\) and \(N\) is given, which are shown to be sharp by numerical tests.
This page was built for publication: Mesh selection for a nearly singular boundary value problem
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5960024)
