Гибридные Схемы На Специальных Сетках Для Решения Задач Со Степенными Слоями

Ранее авторами предложены две гибридные схемы для численного решения задач с экспоненциальными слоями, имеющие на явно задаваемых сетках, генерируемых специальным координатным преобразованием, второй порядок точности и сходящиеся равномерно относительно малого параметра к точным решениям краевых задач для обыкновенных дифференциальных уравнений второго порядка с малым параметром диффузии. В данной работе показано, что эти схемы применимы также для решения задач с иными типами слоев (степенных первого и второго рода и гибридных, представляющих собой сочетания простых) и проявляют при их моделировании ровно те же качества. Проведен анализ погрешности схем внутри и вне слоев при умеренных и стремящихся к нулю значениях малого параметра. Теоретические результаты подтверждены численными экспериментами на тестовых задачах со слоями различных типов, приведены апостериорные оценки ошибок решения и реальных порядков точности, подтверждающие равномерный характер сходимости со вторым порядком. The paper addresses the development of numerical algorithms based on layer-resolving grids and special hybrid difference schemes for solving one-dimensional stationary problems with a small diffusion coefficient that have boundary and interior layers. It describes coordinate transformations and corresponding layer-resolving grids, and analyses the convergence of numerical solutions by using the hybrid schemes on the layer-resolving grids. Theoretical results are confirmed by numerical experiments. The layers are classified by the qualitative behaviour of the corresponding singular functions, in particular by their derivatives. At present, basic types of layers (exponential, power-of-type-1 and type-2, and logarithmic) and their combinations (hybrid types of layers) are known for solving problems with small parameters. Of course, solutions of real problems, in particular to the Navier – Stokes equations, may have new types of layers that have not yet been discovered. The relevance and novelty of the paper are concluded on two powerful factors for increasing effici- encies of algorithms, namely, adaptive grids with special rules for grid clustering in the zones of layers which provide uniform convergence of numerical solutions and hybrid difference schemes. Synthe- sis of the hybrid schemes and special meshes constructed using a priori information about the behaviour of solutions in layers and eliminating layers in transformed coordinates allows achieving uniform convergence to the second order and to increase the real accuracy of approximate solutions practically without complicating the algorithms typical of homogeneous schemes of the up-wind type. The developed technologies of generation of special meshes and construction of difference schemes can be applied for numerical solution of problems with boundary and internal layers of various types, including modelling of viscous fluid and gas flows.