Все выпуски

Численно исследуется явление возникновения скоростного бафтинга при обтекании профиля NACA0012 трансзвуковым потоком. Формулируется математическая модель, основанная на алгоритмах высокого порядка аппроксимации, позволяющая рассчитывать нестационарные отрывные течения. Модель базируется на интегрировании квазигидродинамических уравнений. Проводится параметрическое исследование обтекания профиля высокоскоростным потоком вязкого газа в зависимости от угла атаки. Анализируются как мгновенные, так и осредненные картины течения. Получены распределения пульсационных характеристик течений при различных углах атаки. Выявляются закономерности возникновения отрыва пограничного слоя, определено влияние скачков уплотнения на характер течения вблизи поверхности профиля. Определяется критический угол атаки, при котором начинает иметь место скоростной бафтинг.

Numerically, the phenomenon of the appearance of high-speed bufting is investigated for the case of a transonic flow past the NACA0012 airfoil. A mathematical model based on high-order approximation algorithms is formulated, which makes it possible to calculate nonstationary separated flows. The model is based on the integration of quasi-hydrodynamic equations. A parametric investigation of high-velocity viscous gas flow past an airfoil as a function of the angle of attack is carried out. Both instantaneous and averaged flow patterns are analyzed. The distributions of the pulsation characteristics of flows are obtained at different angles of attack. Regularities in the onset of detachment of the boundary layer are revealed, and the effect of shock waves on the nature of the flow near the surface of the airfoil is determined. The critical angle of attack at which high-speed bufting begins is determined.

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

In this paper, we numerically examine the separated flow in an overexpanded nozzle developed during the starting phase of a rocket engine featuring the side-load effect. We use detached eddy simulation (DES), namely, its version referred to as delayed DES, and compare our results with the existing experimental and numerical data. The calculations are based on finite volume approximation of the Navier-Stokes equations for perfect gas and pressure-based formulation with PISO coupling. By investigating the unsteadiness associated with the shock-induced separation, a better insight can be obtained as to how these affect the nozzle start-up. Three-dimensional, transient, turbulent computational fluid dynamics methodology has been demonstrated to capture major side-load physics with supersonic nozzles. The main flow properties (wall pressure levels, shock waves pattern, side load locus) are rather well reproduced. We suggest using a delayed detached eddy simulation for more accurate predictions of the flow separation.

В работе построен алгоритм повышенного порядка точности на основе WENO схем для моделирования динамики многокомпонентного реагирующего газа с учетом процессов диффузии, теплопроводности и химических реакций. Проведены расчеты для течения газа в проточном реакторе для термического пиролиза этана с внешним обогревом реакционной зоны. В рассматриваемых течениях скорость движения газа много меньше скорости распространения звука в газовой смеси, что обуславливает использование уравнений Навье-Стокса в приближении малых чисел Маха для описания исследуемых процессов. Расчет уравнений химических реакций выделяется в отдельный шаг, где скорость реакции определяется на основе выражений Аррениуса. Для построения модели химической кинетики принята кинетическая схема пиролиза этана, представляющая собой разветвленный радикальный механизм. Проведены расчеты дозвукового течения газа с учетом процессов диффузии, химических реакций и их тепловых эффектов для различных температур нагревательных элементов. Сравнение с экспериментальными данными показало, что $1.97\,\%$-ная конверсия этана в расчетах достигается для $648\,^{\circ}$C на выходе металлического реактора, что близко к экспериментальным значениям, составляющим $2.1\,\%$. Сравнение данных экспериментов по термическому пиролизу этана с данными, полученными в ходе вычислительного эксперимента, показало высокую степень достоверности полученных результатов.

The article considers a high-order accuracy algorithm for modelling the dynamics of multicomponent reactive gas taking into account the processes of diffusion, thermal conductivity and chemical reactions, based on WENO schemes. Computations for gas flow in a flowing reactor for thermal ethane pyrolysis with external heating of the reaction zone are carried out. The velocity of gas motion in explored flows is much less then sound velocity in gas mixture, which motivates using the Navier-Stokes equations in approximation of low Mach numbers for describing the processes under study. Computation of chemical kinetics equations is singled out as a separate step. The velocity of chemical reactions is defined by Arrhenius expressions. The ethane pyrolysis kinetic scheme is used for constructing the model, which is a branched radical mechanism. Computations of subsonic gas flow taking into account the processes of diffusion, chemical reactions and their thermal effects for different temperature of heating elements are carried out. Comparison with experimental data shows that $1.97\,\%$ conversion of ethane is reached at $648\,^{\circ}$C at the outflow of metal reactor. This result is close to $2.1\,\%$, which is obtained by experiment. Comparison of experimental data of thermal ethane pyrolysis with numerical experimental data shows a high level of reliability of the results obtained.

Описаны результаты линейного анализа устойчивости плоскопараллельного течения несжимаемой жидкости над слоем насыщенной пористой среды при различных значениях ее пористости. Рассматривается ограниченная двухслойная система, состоящая из слоя однородной недеформируемой пористой среды конечной толщины и слоя несжимаемой однородной жидкости над ним. Пористый слой ограничен снизу твердой стенкой, верхняя граница жидкости рассматривается как свободная, но недеформируемая. Выполнен анализ линейной устойчивости стационарного течения в такой системе в условиях существования бимодальной нейтральной кривой и варьировании пористости нижнего слоя. Продемонстрирован переход между двумя основными модами неустойчивости: длинноволновой, связанной с точками перегиба в профиле течения, и коротковолновой, обусловленной большим поперечным градиентом скорости течения вблизи границы раздела жидкости и пористой среды. Уменьшение пористости влечет стабилизацию длинноволновых возмущений без существенного изменения критического волнового числа. Коротковолновые возмущения при этом дестабилизируются, а их критическое волновое меняется в широких пределах. При значении пористости меньше 0.7 инерционные слагаемые в уравнении фильтрации и величина механических напряжений на границе раздела возрастают настолько, что доминирующим механизмом развития неустойчивости становится аналог неустойчивости Кельвина-Гельмгольца. В узком интервале пористости реализуется полоса устойчивости течения, разделяющая ветви нейтральной кривой.

The stability of incompressible fluid plane-parallel flow over a layer of a saturated porous medium is studied. The results of a linear stability analysis are described at different porosity values. The considered system is bounded by solid wall from the porous layer bottom. Top fluid surface is free and rigid. A linear stability analysis of plane-parallel stationary flow is presented. It is realized for parameter area where the neutral stability curves are bimodal. The porosity variation effect on flow stability is considered. It is shown that there is a transition between two main instability modes: long-wave and short-wave. The long-wave instability mechanism is determined by inflection points within the velocity profile. The short-wave instability is due to the large transverse gradient of flow velocity near the interface between liquid and porous medium. Porosity decrease stabilizes the long wave perturbations without significant shift of the critical wavenumber. Simultaneously, the short-wave perturbations destabilize, and their critical wavenumber changes in wide range. When the porosity is less than 0.7, the inertial terms in filtration equation and magnitude of the viscous stress near the interface increase to such an extent that the Kelvin-Helmholtz analogue of instability becomes the dominant mechanism for instability development. The stability band realizes in narrow porosity area. It separates the two branches of the neutral curve.

Рассмотрены закрученные ламинарные осесимметричные течения вязких несжимаемых жидкостей в потенциальном поле массовых сил. Исследования течений осуществляются в цилиндрической системе координат. В течениях отдельно рассматриваются области, в которых осевая производная окружной скорости не может принимать нулевое значение в какой-нибудь открытой окрестности (существенно закрученные течения), и области, в которых эта производная равна нулю (область со слоистой закруткой). Показано, что для областей со слоистой закруткой можно применять известный метод (метод вязких вихревых доменов), разработанный для незакрученных течений. Для существенно закрученных течений получена формула для вычисления радиально-осевой скорости воображаемой жидкости через окружную компоненту завихренности, окружную циркуляцию реальной жидкости и частные производные этих функций. Частицы этой воображаемой жидкости «переносят» вихревые трубки радиально-осевой составляющей завихренности с сохранением интенсивности этих трубок, а также «переносят» величину окружной циркуляции и произведение окружной составляющей завихренности на некоторую функцию расстояния до оси симметрии. Предложен неинтегральный способ восстановления поля скорости по полю завихренности. Он сводится к решению системы линейных алгебраических уравнений с двумя переменными. Полученный результат предлагается использовать для распространения метода вязких вихревых доменов на закрученные осесимметричные течения.

Swirling laminar axisymmetric flows of viscous incompressible fluids in a potential field of body forces are considered. The study of flows is carried out in a cylindrical coordinate system. In the flows, the regions in which the axial derivative of the circumferential velocity cannot take on zero value in some open neighborhood (essentially swirling flows) and the regions in which this derivative is equal to zero (the region with layered swirl) are considered separately. It is shown that a well-known method (the method of viscous vortex domains) developed for non-swirling flows can be used for regions with layered swirling. For substantially swirling flows, a formula is obtained for calculating the radial-axial velocity of an imaginary fluid through the circumferential vorticity component, the circumferential circulation of a real fluid, and the partial derivatives of these functions. The particles of this imaginary fluid “transfer” vortex tubes of the radial-axial vorticity component while maintaining the intensity of these tubes, and also “transfer” the circumferential circulation and the product of the circular vorticity component by some function of the distance to the axis of symmetry. A non-integral method for reconstructing the velocity field from the vorticity field is proposed. It is reduced to solving a system of linear algebraic equations in two variables. The obtained result is proposed to be used to extend the method of viscous vortex domains to swirling axisymmetric flows.

В работе рассматриваются результаты решения задачи стационарного течения вязкой несжимаемой жидкости в плоском канале с обратным уступом и прогреваемой нижней стенкой в широком диапазоне числа Рейнольдса $100\leqslant \text{Re}\leqslant 1000$ и параметра расширения потока $1.11 \leqslant ER \leqslant 10$. Исследование выполнено путем численного интегрирования системы двумерных уравнений Навье-Стокса в переменных «скорость-давление» на равномерных сетках с шагом 1/300. Достоверность полученных результатов подтверждается их сравнением с литературными данными. Приводятся подробные картины течения и перегрева жидкости в зависимости от двух основных параметров задачи: $\text{Re}$ и $ER$. Показывается, что с одновременным ростом параметров $\text{Re}$ и $ER$ существенно усложняется структура течения - увеличиваются количество вихрей и их размеры вплоть до образования вихря за уступом с двумя центрами вращения. Также показывается, что характерная высота зоны прогрева течения слабо зависит от $\text{Re}$ и $ER$ и в конечном счете ближе к выходу из канала составляет приблизительно половину его высоты. Для всех центров вихрей определяются их основные характеристики: координаты, экстремумы функции тока, завихренности. Анализируется сложное немонотонное поведение профилей коэффициентов трения, сопротивления и теплоотдачи (числа Нуссельта) по длине канала. Показывается, что эти коэффициенты в одинаковой степени сильно зависят как от числа Рейнольдса, так и от параметра расширения канала, достигая своих максимальных значений при максимальных значениях $\text{Re}$ и $ER$.

The paper deals with the results of solving the problem of steady-state flow of a viscous incompressible fluid in a plane channel with a backward-facing step and a heated bottom wall for the Reynolds number in the range $100\leqslant \text{Re}\leqslant1000$ and the expansion ratio of a plane channel in the range $1.11 \leqslant ER \leqslant 10$. The study was carried out by numerical integration of the 2-D Navier-Stokes equations in velocity-pressure formulation on uniform grids with a step which equals to 1/300. Correction of the results is confirmed by comparing them with the literature data. Detailed flow patterns and fields of stream overheating depending on two basic parameters of the problem $\text{Re}$ and $ER$ are demonstrated. It is shown that with the increase of parameters $\text{Re}$ and $ER$ the structure of flow becomes much more complicated, that is, there is an increase of the number of vortices and their sizes up to the formation of a vortex behind the backward-facing step with two centers of rotation. It is also stated that the typical height of the heating zone of the flow depends weakly on $\text{Re}$ and $ER$ and eventually, near the exit of the channel, equals approximately half of the channel height. For all centers of vortices their main characteristics are defined: location, extremums of stream function, vorticity. Complex nonmonotonic behaviors of the coefficients of friction, hydrodynamic resistance and heat transfer (Nusselt number) along the channel are analyzed. It is shown that these coefficients strongly depend both on Reynolds number and on expansion ratio, reaching the maximum values at the maximum values of $\text{Re}$ and $ER$.

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

The problem of process automation for the development of information systems based on the application of the original neural network architecture is considered. An analysis of existing approaches to the automation of information systems design is carried out. Recommendations for the information systems architecture, aimed at reducing the negative impact of human factor, are formulated. A general concept of neural network architecture in the form of a structural model is presented. Definitions of the main entities and components are given. The key differences of the neural network architecture are: the independence of the key entities of information systems and the possibility of automation of their design and interaction based on the use of neural networks; isolation of the mathematical software of architecture; separation of models of information processes and functional elements from control systems and information representation systems; taking into account the influence of the environment on the processes of movement of information flows, the elements of control and system visualization; the possibility of adapting structural units of information systems to the characteristics of the subject area, the parameters of user equipment without the need to make significant changes to the architecture. The concept of a neural network data channel, its structure and generalized mathematical software are considered. The decomposition of the structural mode is implemented. The structural diagrams of each entity of the neural network architecture of information systems, the description of the main components, the neural network data channels used to connect the entities and their components are presented. The scope of application of the neural network architecture is analyzed.