Compact Difference Schemes and Layer Resolving Grids for Numerical Modeling of Problems with Boundary and Interior Layers

• Published in: Numerical analysis and applications Vol. 12; no. 1; pp. 37 - 50
• Main Authors: Liseikin, V. D., Paasonen, V. I.
• Format: Journal Article
• Language: English
• Published: Moscow Pleiades Publishing 2019; Springer Nature B.V
• Subjects: Boundary layers; Coordinate transformations; Formulations; Interpolation; Mathematical models; Mathematics; Mathematics and Statistics; Numerical Analysis; Singularity (mathematics); compact scheme; boundary layer; interior layer; equation with a small parameter; adaptive grid; high-order scheme; layer resolving grid
• ISSN: 1995-4239; 1995-4247
• DOI: 10.1134/S199542391901004X

A combination of two approaches to numerically solving second-order ODEs with a small parameter and singularities, such as interior and boundary layers, is considered, namely, compact high-order approximation schemes and explicit generation of layer resolving grids. The generation of layer resolving grids, which is based on estimates of solution derivatives and formulations of coordinate transformations eliminating the solution singularities, is a generalization of a method for a first-order scheme developed earlier. This paper presents formulas of the coordinate transformations and numerical experiments for first-, second-, and third-order schemes on uniform and layer resolving grids for equations with boundary, interior, exponential, and power layers of various scales. Numerical experiments confirm the uniform convergence of the numerical solutions performed with the compact high-order schemes on the layer resolving grids. By using transfinite interpolation or numerical solutions to the Beltrami and diffusion equations in a control metric based on coordinate transformations eliminating the solution singularities, this technology can be generalized to the solution of multi-dimensional equations with boundary and interior layers.