Все выпуски

Работа посвящена использованию бикомпактных схем для численного решения эволюционных уравнений гиперболического типа. Основным преимуществом схем этого класса является сочетание двух положительных свойств: пространственной аппроксимации высокого четного порядка на шаблоне, всегда занимающем одну ячейку сетки, и спектрального разрешения, лучшего по сравнению с классическими компактными конечно-разностными схемами того же порядка пространственной аппроксимации. Рассматривается одна особенность бикомпактных схем — жесткая привязка их пространственной аппроксимации к декартовым сеткам (с ячейками-параллелепипедами в трехмерном случае). Она делает затруднительным применение бикомпактных схем к решению задач в сложных расчетных областях в рамках подхода неструктурированных сеток. Предлагается решать эту проблему путем применения известных методов аппроксимации границ сложной формы и соответствующих им краевых условий на декартовых сетках. Обобщение бикомпактных схем на задачи в геометрически сложных областях проводится на примере задач газовой динамики и уравнений Эйлера. В качестве конкретного метода, позволяющего учесть на декартовых сетках влияние твердых границ произвольной формы на течение газа, выбирается метод свободной границы. Приводится краткое описание этого метода, выписываются его уравнения. Для них строятся бикомпактные схемы четвертого порядка аппроксимации по пространству с локально-одномерным расщеплением. Компенсационный поток метода свободной границы дискретизируется со вторым порядком точности. Для интегрирования по времени в получаемых схемах применяются неявный метод Эйлера и $L$-устойчивый жестко-точный трехстадийный однократно диагонально-неявный метод Рунге–Кутты третьего порядка точности. Разработанные бикомпактные схемы тестируются на трех двумерных задачах: о стационарном сверхзвуковом обтекании с числом Маха, равным трем, одного круглого цилиндра и группы изт рех круглых цилиндров, а также о нестационарном взаимодействии плоской ударной волны и круглого цилиндра в канале с плоскопараллельными стенками. Полученные результаты хорошо согласуются с результатами других работ: твердые тела физически корректно влияют на поток газа, давление в контрольных точках на поверхностях тел рассчитывается с точностью, в целом отвечающей выбранному разрешению сетки и уровню численной диссипации.

This work is dedicated to application of bicompact schemes to numerical solution of evolutionary hyperbolic equations. The main advantage of this class of schemes lies in combination of two beneficial properties: the first one is spatial approximation of high even order on a stencil that always occupies only one mesh cell; the second one is spectral resolution which is better in comparison to classic compact finite-difference schemes of the same order of spatial approximation. One feature of bicompact schemes is considered: their spatial approximation is rigidly tied to Cartesian meshes (with parallelepiped-shaped cells in three-dimensional case). This feature makes rather challenging any application of bicompact schemes to problems with complex computational domains as treated in the framework of unstructured meshes. This problem is proposed to be solved using well-known methods for treating complex-shaped boundaries and their corresponding boundary conditions on Cartesian meshes. The generalization of bicompact schemes on problems in geometrically complex domains is made in case of gas dynamics problems and Euler equations. The free boundary method is chosen as a particular tool to introduce the influence of arbitrary-shaped solid boundaries on gas flows on Cartesian meshes. A brief description of this method is given, its governing equations are written down. Bicompact schemes of fourth order of approximation in space with locally one-dimensional splitting are constructed for equations of the free boundary method. Its compensation flux is discretized with second order of accuracy. Time stepping in the obtained schemes is done with the implicit Euler method and the third order accurate $L$-stable stiffly accurate three-stage singly diagonally implicit Runge–Kutta method. The designed bicompact schemes are tested on three two-dimensional problems: stationary supersonic flows with Mach number three past one circular cylinder and past three circular cylinders; the non-stationary interaction of planar shock wave with a circular cylinder in a channel with planar parallel walls. The obtained results are in a good agreement with other works: influence of solid bodies on gas flows is physically correct, pressure in control points on solid surfaces is calculated with the accuracy appropriate to the chosen mesh resolution and level of numerical dissipation.