Interpolation on the Bakhvalov mesh in the presence of an exponential boundary layer

I.A. Blatova∗ , N.A. Zadorinb∗∗

aPovolzhskiy State University of Telecommunications and Informatics, Samara, 443010 Russia

bSobolev Institute of Mathematics, Siberian Branch, Russian Academy of Sciences, Novosibirsk, 630090 Russia

E-mail: ∗blatow@mail.ru, ∗∗nik-zadorin@yandex.ru

Received September 26, 2019

Full text pdf

DOI: 10.26907/2541-7746.2019.4.497-508

For citation: Blatov I.A., Zadorin N.A. Interpolation on the Bakhvalov mesh in the presence of an exponential boundary layer. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2019, vol. 161, no. 4, pp. 497–508. doi: 10.26907/2541-7746.2019.4.497-508. (In Russian)

 

Abstract

Interpolation of the function of one variable with large gradients in the region of the exponential boundary layer was studied. The interpolated function corresponds to solution of a boundary value problem for an ordinary differential equation of the second order with the small parameter ε before the highest derivative. Applying classical polynomial interpolation formulas on a uniform mesh to this function can lead to unacceptable errors. In the paper, the error of the piecewise linear interpolation formula on the Bakhvalov mesh condensing in the region of the boundary layer was estimated. The Bakhvalov mesh is used in a number of works when constructing difference schemes for singularly perturbed problems; therefore, estimating the error of interpolation formulas on this mesh is of interest. An error estimate of the order of O(1/N 2) was obtained uniformly with respect to the parameter ε, where N is the number of mesh nodes. The problem of computing the derivative of the function with large gradients given in the nodes of the Bakhvalov mesh was investigated. The classical difference formula with two nodes was considered obtained by differentiating the linear interpolant studied above. An estimate of the relative error of the order of O(1/N ), uniform in the parameter ε, was obtained. The results of the numerical experiments consistent with the obtained error estimates were presented. Numerical comparison of the errors obtained during the interpolation and numerical differentiation on the Bakhvalov mesh with errors on the Shishkin mesh and on the uniform mesh was carried out.

Keywords: function of one variable, boundary layer, Bakhvalov mesh, piecewise linear interpolation, numerical differentiation, ε -uniform error estimation

Acknowledgments. The study was supported by the Russian Foundation for Basic Research (project no. 19-31-60009).

References

  1. Linß T. Layer-Adapted Meshes for Reaction-Convection-Diffusion Problems. Berlin, Springer, 2010. 233 p.
  2. Zadorin A.I. Interpolation method for the boundary-layer problem. Sib. Zh. Vychisl. Mat., 2007, vol. 10, no. 3, pp. 267–275. (In Russian)
  3. Bakhvalov N.S. The optimization of methods of solving boundary value problems with a boundary layer. USSR Comput. Math. Math. Phys., 1969, vol. 9, no. 4, pp. 139–166. doi: 10.1016/0041-5553(69)90038-X.
  4. Shishkin G.I. Setochnye approksimatsii singulyarno vozmushchennykh ellipticheskikh i parabolicheskikh uravnenii [Grid Approximation of Singularly Perturbed Elliptic and Parabolic Equations]. Yekaterinburg, Ural. Otd. Ross. Akad. Nauk, 1992. 233 p. (In Russian)
  5. Zadorin A.I. Lagrange interpolation and Newton-Cotes formulas for functions with boundary layer components on piecewise-uniform grids. Numer. Anal. Appl., 2015, vol. 8, no. 3, pp. 235–247. doi: 10.1134/S1995423915030040.
  6. Zadorin A.I. Analysis of numerical differentiation formulas in a boundary layer on a Shishkin grid. Numer. Anal. Appl., 2018, vol. 11, no. 3, pp. 193–203. doi: 10.1134/S1995423918030011.
  7. Linß T. The necessity of Shishkin decompositions. Appl. Math. Lett., 2001, vol. 14, no. 7, pp. 891–896. doi: 10.1016/S0893-9659(01)00061-1.
  8. Zadorin A.I., Zadorin N.A. Spline interpolation on a uniform grid for functions with a boundary-layer component. Comput. Math. Math. Phys., 2010, vol. 50, no. 2, pp. 211–223. doi: 10.1134/S0965542510020028.
  9. Zadorin A.I., Zadorin N.A. Interpolation formula for functions with a boundary layer component and its application to derivatives calculation. Sib. Electron. Math. Rep., 2012, vol. 9, pp. 445–455.
  10. Blatov I.A., Zadorin N.A. Analysis of the interpolation formula exact on the boundary layer component of the interpolated function. Nauka Mir, 2015, vol. 1, no. 2, pp. 13–17. (In Russian)

 

The content is available under the license Creative Commons Attribution 4.0 License.