Все выпуски

В развитии науки важную роль играют научные школы — объединения исследователей, связанные общей проблемой, идеями и методами, используемыми для решения проблемы. Научные школы формируются вокруг лидера и объединяющей идеи.

За время научной деятельности академика А. С. Холодова вокруг него сформировалось несколько научных школ. В обзоре делается попытка представить основные научные направления, вокруг которых сформировались яркие коллективы с общими системами взглядов и подходами к исследованиям. В обзоре отмечается эта общая основа. Во-первых, это развитие группы численных методов для решения систем дифференциальных уравнений в частных производных гиперболического типа — сеточно-характеристические методы. Во-вторых, описание численных методов в пространствах неопределенных коэф- фициентов. Этот подход развивался как для всех типов уравнений в частных производных, так и для обыкновенных дифференциальных уравнений.

На основе предложенных А. С. Холодовым численных подходов сложились научные коллективы, работающие в разных предметных областях. Это математическое моделирование динамики плазмы, динамики деформируемого твердого тела, некоторых задач биологии, биофизики, медицинской физики и биомеханики. Сравнительно новые направления — решение задач на графах (процессы транспортировки электроэнергии, моделирование транспортных потоков на дорожной сети и т. д.).

В обзоре делается попытка отследить деятельность научных школ от момента их зарождения до настоящего времени, проследить связь работ А. С. Холодова с работами его учеников и коллег. Полный обзор деятельности всех научных школ, сформировавшихся вокруг Александра Сергеевча, невозможен ввиду огромного количества и разнообразия научных результатов.

Делается также попытка связать деятельность научных школ с появлением научно-образовательной школы в Московском физико-техническом институте.

In the science development an important role the scientific schools are played. This schools are the associations of researchers connected by the common problem, the ideas and the methods used for problems solution. Usually Scientific schools are formed around the leader and the uniting idea.

The several sciences schools were created around academician A. S. Kholodov during his scientific and pedagogical activity.

This review tries to present the main scientific directions in which the bright science collectives with the common frames of reference and approaches to researches were created. In the review this common base is marked out. First, this is development of the group of numerical methods for hyperbolic type systems of partial derivatives differential equations solution — grid and characteristic methods. Secondly, the description of different numerical methods in the undetermined coefficients spaces. This approach developed for all types of partial equations and for ordinary differential equations.

On the basis of A. S. Kholodov’s numerical approaches the research teams working in different subject domains are formed. The fields of interests are including mathematical modeling of the plasma dynamics, deformable solid body dynamics, some problems of biology, biophysics, medical physics and biomechanics. The new field of interest includes solving problem on graphs (such as processes of the electric power transportation, modeling of the traffic flows on a road network etc).

There is the attempt in the present review analyzed the activity of scientific schools from the moment of their origin so far, to trace the connection of A. S. Kholodov’s works with his colleagues and followers works. The complete overview of all the scientific schools created around A. S. Kholodov is impossible due to the huge amount and a variety of the scientific results.

The attempt to connect scientific schools activity with the advent of scientific and educational school in Moscow Institute of Physics and Technology also becomes.

В работе рассматриваются результаты применения неявного итерационного полинейного рекуррентного метода решения систем разностных эллиптических уравнений, возникающих при численном моделировании динамики несжимаемой вязкой жидкости. Исследование проводится на примере решения задачи о стационарном течении в плоской каверне с подвижной крышкой, сформулированной в естественных переменных ($u, \,v, \,p$) при больших значениях чисел Re (до 20 000) и сеточных разрешений (до 2049×2049). Демонстрируется высокая эффективность метода при расчете полей поправки давления. Анализируются проблемы решения задачи при больших числах Re.

In the paper the results of applying the implicit line-by-line recurrence method for solving of systems of elliptic difference equations, arising, in particular, at numerical simulation of dynamics of incompressible viscous fluid are considered. Research is conducted on the example of the problem about a steady-state two-dimensional lid-driven cavity flow formulated in primitive variables ($u,\, v,\, p$) for large Re (up to 20 000) and grids (up to 2049×2049). High efficiency of the method at calculation of a pressure correction fields is demonstrated. The difficulties of constructing a solution of the problem for large Rе are analyzed.

Предложено обобщение блочного клеточного автомата Марголуса на гексагональную сетку. Проведена статистическая обработка результатов вероятностных клеточно-автоматных вычислений для ряда модификаций схемы, решающей тестовую задачу диффузии вещества. Показано, что выбор блоков в виде гексагонов на 25% эффективнее, чем в виде Y-блоков. Показано, что алгоритмы имеют полиномиальную сложность, причем степень полинома для параллельных вычислителей лежит в пределах 0.6÷0.8, а для последовательных — в пределах 1.5÷1.7. Исследовалось влияние внедренных в поле клеточного автомата дефектных ячеек на скорость сходимости.

The generalization of Margolus’s block cellular automaton on a hexagonal grid is formulated. Statistical analysis of the results of probabilistic cellular automation for vast variety of this scheme solving the test task of diffusion is done. It is shown that the choice of the hexagon blocks is 25% more efficient than Y-blocks. It is shown that the algorithms have polynomial complexity, and the polynom degree lies within 0.6÷0.8 for parallel computer, and in the range 1.5÷1.7 for serial computer. The effects of embedded into automaton’s field defective cells on the rate of convergence are studied also.

В работе представлена новая версия программного комплекса FlowVision, предназначенного для автоматизации инженерных расчетов в области вычислительной гидродинамики: FlowVision 3.09.05. Программный комплекс (ПК) FlowVision используется для решения различных прикладных задач в различных областях промышленности. Его популярность основана на том, что он позволяет решать сложные нетрадиционные задачи, находящиеся на стыке различных дисциплин, с одной стороны, и, с другой стороны, на парадигме полной автоматизации таких трудоемких для инженера процессов, как построение расчетной сетки. FlowVision — это программный комплекс, полностью отчуждаемый от разработчиков. Он имеет развитый графический интерфейс, систему задания расчетного проекта и систему визуализации течений различными методами — от построения контуров (для скалярных переменных) и векторов (для векторных переменных) на плоскостях и поверхностях до объемной визуализации расчетных данных. Кроме этого, ПК FlowVision предоставляет пользователю возможность вычислять интегральные характеристики на поверхностях и в ограниченных объемах.

ПК основан на конечно-объемном подходе к аппроксимации основных уравнений движения жидкости. В нем реализованы явный и неявный методы решения этих уравнений. ПК имеет автоматический построитель неструктурированной сетки с возможностью ее локальной динамической адаптации. В ПК реализован двухуровневый параллелизм, позволяющий эффективно проводить расчеты на компьютерах, имеющих распределенную и общую память одновременно. FlowVision обладает широким спектром физико-математических моделей: турбулентности (URANS, LES, ILES), горения, массопереноса с учетом химических превращений и радиоактивного распада, электрогидродинамики.

FlowVision позволяет решать задачи движения жидкостей со скоростями, соответствующими несжимаемому или гиперзвуковому режимам за счет использования все-скоростного метода расщепления по физическим переменным для решения уравнений Навье–Стокса. FlowVision позволяет решать междисциплинарные задачи с использованием различных средств моделирования, например: моделировать многофазные течения методом VOF, обтекание подвижных тел с помощью эйлерова подхода при неподвижной расчетной сетке, моделировать вращающиеся машины с использованием метода скользящей сетки, решать задачи взаимодействия жидкости и конструкций методом двухстороннего сопряжения FlowVision с конечно-элементными кодами. В данной работе показаны примеры решения задач-вызовов: a) посадка космического корабля на воду при торможении ракетными двигателями, где есть граница раздела «воздух–вода», подвижные тела и взаимодействие сверхзвуковой струи газа с границей раздела «вода–воздух»; б) моделирование работы человеческого сердца с искусственными и живыми клапанами, спроектированными на базе томографических исследований, с использованием двухстороннего сопряжения «жидкостной» расчетной области с конечно-элементной моделью мышц сердца.

The work submits new release of the FlowVision software designed for automation of engineering calculations in computational fluid dynamics: FlowVision 3.09.05. The FlowVision software is used for solving different industrial problems. Its popularity is based on the capability to solve complex non-tradition problems involving different physical processes. The paradigm of complete automation of labor-intensive and time-taking processes like grid generation makes FlowVision attractive for many engineers. FlowVision is completely developer-independent software. It includes an advanced graphical interface, the system for specifying a computational project as well as the system for flow visualization on planes, on curvilinear surfaces and in volume by means of different methods: plots, color contours, iso-lines, iso-surfaces, vector fields. Besides that, FlowVision provides tools for calculation of integral characteristics on surfaces and in volumetric regions.

The software is based on the finite-volume approach to approximation of the partial differential equations describing fluid motion and accompanying physical processes. It provides explicit and implicit methods for time integration of these equations. The software includes automated generator of unstructured grid with capability of its local dynamic adaptation. The solver involves two-level parallelism which allows calculations on computers with distributed and shared memory (coexisting in the same hardware). FlowVision incorporates a wide spectrum of physical models: different turbulence models, models for mass transfer accounting for chemical reactions and radioactive decay, several combustion models, a dispersed phase model, an electro-hydrodynamic model, an original VOF model for tracking moving interfaces. It should be noted that turbulence can be simulated within URANS, LES, and ILES approaches. FlowVision simulates fluid motion with velocities corresponding to all possible flow regimes: from incompressible to hypersonic. This is achieved by using an original all-speed velocity-pressure split algorithm for integration of the Navier-Stokes equations.

FlowVision enables solving multi-physic problems with use of different modeling tools. For instance, one can simulate multi-phase flows with use of the VOF method, flows past bodies moving across a stationary grid (within Euler approach), flows in rotary machines with use of the technology of sliding grid. Besides that, the software solves fluid-structure interaction problems using the technology of two-way coupling of FlowVision with finite-element codes. Two examples of solving challenging problems in the FlowVision software are demonstrated in the given article. The first one is splashdown of a spacecraft after deceleration by means of jet engines. This problem is characterized by presence of moving bodies and contact surface between the air and the water in the computational domain. The supersonic jets interact with the air-water interphase. The second problem is simulation of the work of a human heart with artificial and natural valves designed on the basis of tomographic investigations with use of a finite-element model of the heart. This problem is characterized by two-way coupling between the “liquid” computational domain and the finite-element model of the hart muscles.

Во второй части статьи, носящей более прикладной характер, завершается рассмотрение трех классических уравнений математической физики (Лапласа, диффузии и волнового) простейшими численными схемами в формулировке клеточных автоматов (КА). На нескольких примерах, относящихся к гексагональной сетке, показана специфика такого решения и подтверждаются выводы первой части, в частности о выполнении свойства консервативности и эффекте избыточной гексагональной симметрии (ИГС).

При решении задачи Неймана для колебаний круглой мембраны показана критичность требований к дискретизации условий для граничных КА-ячеек. Для квазиодномерной задачи «диффузия в полупространство» сравниваются КА-расчеты, проводимые по простой схеме и с использованием обобщенного блочно-поворотного механизма Марголуса. При решении смешанной задачи для классического случая колебания круглой мембраны с закрепленными концами показано, что одновременное применение метода Кранка–Николсон и учет членов второго порядка позволяет избежать ИГС-эффекта, наблюдаемого нами для более простой схемы. С точки зрения КА центральное место занимает уравнение диффузии, на пути решения которого на бесконечных временах находится решение краевой задачи для уравнения Лапласа, а путем введения вектор-переменной становится разрешимо волновое уравнение (по крайней мере скалярное).

На примере центрально-симметричной задачи Неймана продемонстрирован новый способ введения пространственных производных в postfix-процедуру КА, отражающую временные производные (основанием является уравнение непрерывности). Для случая центральной симметрии эмпирически найдено значение константы, связывающее эти производные. Показано, что препятствием к применению КА-методов для таких задач являются низкая скорость сходимости и точность, лимитируемая точностью дискретизации границ, а не формальной точностью метода (4-й порядок); наша рекомендация состоит в использовании техники multigrid. При решении квазиодномерного уравнения диффузии (двумерным КА) показано, что блочно-поворотный КА (по механизму Марголуса) более эффективен, чем простой КА.

The second part of paper is devoted to final study of three classic partial differential equations (Laplace, Diffusion and Wave) solution using simple numerical methods in terms of Cellular Automata. Specificity of this solution has been shown by different examples, which are related to the hexagonal grid. Also the next statements that are mentioned in the first part have been proved: the matter conservation law and the offensive effect of excessive hexagonal symmetry.

From the point of CA view diffusion equation is the most important. While solving of diffusion equation at the infinite time interval we can find solution of boundary value problem of Laplace equation and if we introduce vector-variable we will solve wave equation (at least, for scalar). The critical requirement for the sampling of the boundary conditions for CA-cells has been shown during the solving of problem of circular membrane vibrations with Neumann boundary conditions. CA-calculations using the simple scheme and Margolus rotary-block mechanism were compared for the quasione-dimensional problem “diffusion in the half-space”. During the solving of mixed task of circular membrane vibration with the fixed ends in a classical case it has been shown that the simultaneous application of the Crank–Nicholson method and taking into account of the second-order terms is allowed to avoid the effect of excessive hexagonal symmetry that was studied for a simple scheme.

By the example of the centrally symmetric Neumann problem a new method of spatial derivatives introducing into the postfix CA procedure, which is reflecting the time derivatives (on the base of the continuity equation) was demonstrated. The value of the constant that is related to these derivatives has been empirically found in the case of central symmetry. The low rate of convergence and accuracy that limited within the boundaries of the sample, in contrary to the formal precision of the method (4-th order), prevents the using of the CAmethods for such problems. We recommend using multigrid method. During the solving of the quasi-diffusion equations (two-dimensional CA) it was showing that the rotary-block mechanism of CA (Margolus mechanism) is more effective than simple CA.

Работа посвящена теоретическому обоснованию неявного итерационного полинейного рекуррентного метода решения систем разностных уравнений, которые возникают при аппроксимации двумерных эллиптических дифференциальных уравнений на регулярной сетке. Высокая эффективность этого метода практически подтверждена при решении сложных тестовых задач, а также задач течения и теплообмена вязкой несжимаемой жидкости. Однако теоретические положения, объясняющие высокую скорость сходимости и устойчивость метода, до сих пор оставались за кадром внимания, что и послужило причиной проведения настоящего исследования. В работе подробно излагается процедура эквивалентных и приближенных преобразований исходной системы линейных алгебраических уравнений (СЛАУ) как в матрично-векторной форме, так и виде расчетных формул метода. При этом для наглядности изложения материала ключевые моменты преобразований иллюстрируются схемами изменения разностных шаблонов, отвечающих преобразованным уравнениям. Конечная цель процедуры преобразований — получение канонической формы записи метода, из которого следует его корректность в случае сходимости решения. На основе анализа структур и элементных составов матричных операторов проводится оценка их норм и, соответственно, доказывается сходимость метода для произвольных начальных векторов.

В специальном случае слабых ограничений на искомое решение производится оценка нормы оператора перехода. Показывается, что с ростом размерности матрицы этого оператора величина его нормы уменьшается пропорционально квадрату (или кубу, в зависимости от версии метода) шага сеточного разбиения области решения задачи. С помощью простых оценок получено необходимое условие устойчивости метода. Также даются рекомендации относительно выбора по порядку величины оптимального итерационного параметра компенсации. Теоретические выводы проиллюстрированы результатами решения тестовых задач. Показано, что при увеличении размерности сеточного разбиения области решения количество итераций, необходимых для достижения заданной точности решения, при прочих равных условиях уменьшается. Также продемонстрировано, что если слабые ограничения на решение нарушены при выборе его начального приближения, то в полном соответствии с полученными теоретическими результатами скорость сходимости метода существенно уменьшается.

In the article a theory of the implicit iterative line-by-line recurrence method for solving the systems of finite-difference equations which arise as a result of approximation of the two-dimensional elliptic differential equations on a regular grid is stated. On the one hand, the high effectiveness of the method has confirmed in practice. Some complex test problems, as well as several problems of fluid flow and heat transfer of a viscous incompressible liquid, have solved with its use. On the other hand, the theoretical provisions that explain the high convergence rate of the method and its stability are not yet presented in the literature. This fact is the reason for the present investigation. In the paper, the procedure of equivalent and approximate transformations of the initial system of linear algebraic equations (SLAE) is described in detail. The transformations are presented in a matrix-vector form, as well as in the form of the computational formulas of the method. The key points of the transformations are illustrated by schemes of changing of the difference stencils that correspond to the transformed equations. The canonical form of the method is the goal of the transformation procedure. The correctness of the method follows from the canonical form in the case of the solution convergence. The estimation of norms of the matrix operators is carried out on the basis of analysis of structures and element sets of the corresponding matrices. As a result, the convergence of the method is proved for arbitrary initial vectors of the solution of the problem.

The norm of the transition matrix operator is estimated in the special case of weak restrictions on a desired solution. It is shown, that the value of this norm decreases proportionally to the second power (or third degree, it depends on the version of the method) of the grid step of the problem solution area in the case of transition matrix order increases. The necessary condition of the method stability is obtained by means of simple estimates of the vector of an approximate solution. Also, the estimate in order of magnitude of the optimum iterative compensation parameter is given. Theoretical conclusions are illustrated by using the solutions of the test problems. It is shown, that the number of the iterations required to achieve a given accuracy of the solution decreases if a grid size of the solution area increases. It is also demonstrated that if the weak restrictions on solution are violated in the choice of the initial approximation of the solution, then the rate of convergence of the method decreases essentially in full accordance with the deduced theoretical results.

В работе обсуждается возможность моделирования турбулентных сжимаемых течений газа с использованием моделей турбулентности $k-\varepsilon$ стандартная (KES), $k-\varepsilon$ FlowVision (KEFV) и SST $k-\omega$. Представлена новая версия модели турбулентности KEFV. Показаны результаты ее тестирования. Проведено численное исследование истечения сверхзвуковой перерасширенной струи из конического сопла в безграничное пространство. Результаты сравниваются с экспериментальными данными. Демонстрируется зависимость результатов от сетки. Демонстрируется зависимость результатов от турбулентности, задаваемой на входе в сопло. Делается вывод о том, что в двухпараметрических моделях турбулентности необходимо учитывать сжимаемость. Для этого подходит простой способ, предложенный Вилкоксом в 1994 г. В результате область применимости трех указанных двухпараметрических моделей заметно расширяется. Предлагаются конкретные значения констант, управляющих учетом сжимаемости в подходе Вилкокса. Эти значения рекомендуется задавать в моделях KES, KEFV и SST при моделировании сжимаемых течений.

Дополнительно рассмотрен вопрос о том, как получать правильные характеристики сверхзвукового турбулентного течения с использованием двухпараметрических моделей турбулентности. Расчеты на разных сетках показали, что при задании ламинарного потока на входе в сопло и пристеночных функций на его поверхностях ядро потока остается ламинарным вплоть до 5-й бочки. Для получения правильных характеристик нужно либо на входе в расчетную область задавать два параметра, характеризующие турбулентность втекающего потока, либо задавать «затравочную» турбулентность в ограниченной области на выходе из сопла, охватывающей зону предполагаемого ламинарно-турбулентного перехода. Последняя возможность реализована в модели KEFV.

Simulation of turbulent compressible gas flows using turbulence models $k-\varepsilon$ standard (KES), $k-\varepsilon$ FlowVision (KEFV) and SST $k-\omega$ is discussed in the given article. A new version of turbulence model KEFV is presented. The results of its testing are shown. Numerical investigation of the discharge of an over-expanded jet from a conic nozzle into unlimited space is performed. The results are compared against experimental data. The dependence of the results on computational mesh is demonstrated. The dependence of the results on turbulence specified at the nozzle inlet is demonstrated. The conclusion is drawn about necessity to allow for compressibility in two-parametric turbulence models. The simple method proposed by Wilcox in 1994 suits well for this purpose. As a result, the range of applicability of the three aforementioned two-parametric turbulence models is essentially extended. Particular values of the constants responsible for the account of compressibility in the Wilcox approach are proposed. It is recommended to specify these values in simulations of compressible flows with use of models KES, KEFV, and SST.

In addition, the question how to obtain correct characteristics of supersonic turbulent flows using two-parametric turbulence models is considered. The calculations on different grids have shown that specifying a laminar flow at the inlet to the nozzle and wall functions at its surfaces, one obtains the laminar core of the flow up to the fifth Mach disk. In order to obtain correct flow characteristics, it is necessary either to specify two parameters characterizing turbulence of the inflowing gas, or to set a “starting” turbulence in a limited volume enveloping the region of presumable laminar-turbulent transition next to the exit from the nozzle. The latter possibility is implemented in model KEFV.

Статья носит методический характер и посвящена решению трех классических уравнений математической физики (Лапласа, диффузии и волнового) простейшими численными схемами в формулировке клеточных автоматов (КА). Особое внимание уделяется законам сохранения вещества и неприятному эффекту избыточной гексагональной симметрии (ИГС).

Делается вывод о том, что по сравнению с классическими конечно-разностными методами, хотя локальная функция перехода (ЛФП) КА терминологически эквивалентна шаблону вычислительной двухслоевой явной схемы, различие состоит в замене матричных (direct) методов (например, метода прогонки для трехдиагональной матрицы) итерационными. Из этого следуют более жесткие требования к дискретизации условий для граничных КА-ячеек.

Для гексагональной сетки и консервативных граничных условий записана корректная ЛФП для граничных ячеек, справедливая, по крайней мере, для границ прямоугольной и круговой формы. Предложена идея разделения ЛФП на internal, boundary и postfix. На примере этой задачи заново осмыслено значение числа Куранта–Леви как соотношения скорости сходимости КА к решению задачи, данному на фиксированный момент времени, и скорости изменения самого решения в динамике.

The paper has methodical character; it is devoted to three classic partial differential equations (Laplace, Diffusion and Wave) solution using simple numerical methods in terms of Cellular Automata. Special attention was payed to the matter conservation law and the offensive effect of excessive hexagonal symmetry.

It has been shown that in contrary to finite-difference approach, in spite of terminological equivalence of CA local transition function to the pattern of computing double layer explicit method, CA approach contains the replacement of matrix technique by iterative ones (for instance, sweep method for three diagonal matrixes). This suggests that discretization of boundary conditions for CA-cells needs more rigid conditions.

The correct local transition function (LTF) of the boundary cells, which is valid at least for the boundaries of the rectangular and circular shapes have been firstly proposed and empirically given for the hexagonal grid and the conservative boundary conditions. The idea of LTF separation into «internal», «boundary» and «postfix» have been proposed. By the example of this problem the value of the Courant-Levy constant was re-evaluated as the CA convergence speed ratio to the solution, which is given at a fixed time, and to the rate of the solution change over time.

В настоящей работе приводятся результаты численного моделирования свободного сдвигового течения с помощью схемы «КАБАРЕ», реализованной в приближении слабой сжимаемости. Анализ схемы проводится на основе изучения свойств неустойчивости Кельвина–Гельмгольца и порождаемой ею двумерной турбулентности, с использованием интегральных кривых кинетической энергии и энстрофии, картин временной эволюции завихренности, спектров энстрофии и энергии, а также дисперсионного соотношения для инкремента неустойчивости. Расчеты проводились для числа Рейнольдса $\text{Re} = 4 \times 10^5$, на квадратных последовательно сгущаемых сетках в диапазоне $128^2-2048^2$ ячеек. Внимание уделено проблеме «недоразрешенности слоев», проявляющейся в возникновении лишнего вихря при свертывании двух вихревых листов (слоев вихревой пелены). Данное явление существует только на грубых сетках $(128^2)$, однако, полностью симметричная картина эволюции завихренности начинает наблюдаться только при переходе к сетке $1024^2$ ячеек. Размерные оценки отношения вихрей на границах инерционного интервала показывают, что наиболее подробная сетка $2048^2$ ячеек оказывается достаточной для качественного отображения мелкомасштабных сгустков завихренности. Тем не менее можно говорить о достижении хорошей сходимости при отображении крупномасштабных структур. Эволюция турбулентности, в полном соответствии с теоретическими представлениями, приводит к появлению крупных вихрей, в которых сосредотачивается вся кинетическая энергия движения, и уединенных мелкомасштабных образований. Последние обладают свойствами когерентных структур, выживая в процессе нитеобразования (филаментации), и практически не взаимодействуют с вихрями других масштабов. Обсуждение диссипативных характеристик схемы ведется на основе анализа графиков скорости диссипации кинетической энергии, вычисляемой непосредственно, а также на основе теоретических соотношений для моделей несжимаемой жидкости (по кривым энстрофии) и сжимаемого газа (по влиянию тензора скоростей деформации и эффектов дилатации). Асимптотическое поведение каскадов кинетической энергии и энстрофии подчиняется реализующимся в двумерной турбулентности соотношениям $E(k) \propto k^{−3}$, $\omega^2(k) \propto k^{−1}$. Исследование зависимости инкремента неустойчивости от безразмерного волнового числа показывает хорошее согласие с данными других исследователей, вместе с тем часто используемый способ расчета инкремента неустойчивости не всегда оказывается достаточно точным, вследствие чего была предложена его модификация.

Таким образом, реализованная схема, отличаясь малой диссипативностью и хорошим вихреразрешением, оказывается вполне конкурентоспособной в сравнении с методами высокого порядка точности.

In present paper we reexamine the properties of CABARET numerical scheme formulated for a weakly compressible fluid flow basing the results of free shear layer modeling. Kelvin–Helmholtz instability and successive generation of two-dimensional turbulence provide a wide field for a scheme analysis including temporal evolution of the integral energy and enstrophy curves, the vorticity patterns and energy spectra, as well as the dispersion relation for the instability increment. The most part of calculations is performed for Reynolds number $\text{Re} = 4 \times 10^5$ for square grids sequentially refined in the range of $128^2-2048^2$ nodes. An attention is paid to the problem of underresolved layers generating a spurious vortex during the vorticity layers roll-up. This phenomenon takes place only on a coarse grid with $128^2$ nodes, while the fully regularized evolution pattern of vorticity appears only when approaching $1024^2$-node grid. We also discuss the vorticity resolution properties of grids used with respect to dimensional estimates for the eddies at the borders of the inertial interval, showing that the available range of grids appears to be sufficient for a good resolution of small–scale vorticity patches. Nevertheless, we claim for the convergence achieved for the domains occupied by large-scale structures.

The generated turbulence evolution is consistent with theoretical concepts imposing the emergence of large vortices, which collect all the kinetic energy of motion, and solitary small-scale eddies. The latter resemble the coherent structures surviving in the filamentation process and almost noninteracting with other scales. The dissipative characteristics of numerical method employed are discussed in terms of kinetic energy dissipation rate calculated directly and basing theoretical laws for incompressible (via enstrophy curves) and compressible (with respect to the strain rate tensor and dilatation) fluid models. The asymptotic behavior of the kinetic energy and enstrophy cascades comply with two-dimensional turbulence laws $E(k) \propto k^{−3}, \omega^2(k) \propto k^{−1}$. Considering the instability increment as a function of dimensionless wave number shows a good agreement with other papers, however, commonly used method of instability growth rate calculation is not always accurate, so some modification is proposed. Thus, the implemented CABARET scheme possessing remarkably small numerical dissipation and good vorticity resolution is quite competitive approach compared to other high-order accuracy methods

Cтроится семейство явных разностных схем на пятиточечном шаблоне для численного решения линейного уравнения переноса. Анализ свойств разностных схем проводится в пространстве неопределенных коэффициентов. Такие пространства впервые были введены в рассмотрение А. С. Холодовым. Для исследования свойств разностных схем ставилась задача линейного программирования. В качестве целевой функции обычно рассматривался коэффициент при главном члене невязки. Для построения монотонных разностных схем ставилась задача оптимизации с ограничениями типа неравенств. Ограниченность такого подхода становится ясной с учетом того, что аппроксимация разностной схемы определяется лишь на классических (гладких) решениях дифференциальной задачи.

В соответствие разностной схеме ставится некоторый функционал, определяющий свойства разностной схемы. Функционал должен быть линейным по коэффициентам схемы. Возможно, что функционал зависит от сеточной функции — решения разностной задачи или проекции на сетку решения дифференциальной задачи. Если первые члены разложения в ряд Тейлора этого функционала по сеточным параметрам совпадут с условиями классической аппроксимации, такой функционал будем называть обобщенным условием аппроксимации. В статье показано, что такие функционалы существуют. Для линейного уравнения с постоянными коэффициентами построение такого функционала возможно и для обобщенного (негладкого) решения дифференциальной задачи.

Построение разностной схемы с заданными свойствами тогда опирается на решение задачи поиска минимума функционала.

Построены семейства функционалов как для гладких решений исходной дифференциальной задачи, так и для обобщенных решений. Построены новые разностные схемы, основанные на анализе функционалов методами линейного программирования. При этом использован аппарат исследования пары самодвойственных задач линейного программирования. Найдена оптимальная монотонная разностная схема, обладающая первым порядком аппроксимации на гладком решении. Обсуждается возможность применения построенных новых схем для построения гибридных разностных схем повышенного порядка аппроксимации на гладких решениях.

Приводится пример численной реализации простейшей разностной схемы с обобщенной аппроксимацией.

A set of implicit difference schemes on the five-pointwise stensil is under construction. The analysis of properties of difference schemes is carried out in a space of undetermined coefficients. The spaces were introduced for the first time by A. S. Kholodov. Usually for properties of difference schemes investigation the problem of the linear programming was constructed. The coefficient at the main term of a discrepancy was considered as the target function. The optimization task with inequalities type restrictions was considered for construction of the monotonic difference schemes. The limitation of such an approach becomes clear taking into account that approximation of the difference scheme is defined only on the classical (smooth) solutions of partial differential equations.

The functional which minimum will be found put in compliance to the difference scheme. The functional must be the linear on the difference schemes coefficients. It is possible that the functional depends on net function – the solution of a difference task or a grid projection of the differential problem solution. If the initial terms of the functional expansion in a Taylor series on grid parameters are equal to conditions of classical approximation, we will call that the functional will be the generalized condition of approximation. It is shown that such functionals exist. For the simple linear partial differential equation with constant coefficients construction of the functional is possible also for the generalized (non-smooth) solution of a differential problem.

Families of functionals both for smooth solutions of an initial differential problem and for the generalized solution are constructed. The new difference schemes based on the analysis of the functionals by linear programming methods are constructed. At the same time the research of couple of self-dual problems of the linear programming is used. The optimum monotonic difference scheme possessing the first order of approximation on the smooth solution of differential problem is found. The possibility of application of the new schemes for creation of hybrid difference methods of the raised approximation order on smooth solutions is discussed.

The example of numerical implementation of the simplest difference scheme with the generalized approximation is given.