Все выпуски

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

Conjugate natural convection in a porous medium saturated with a gas surrounded by the finite thickness solid walls at presence of a local heat source has been numerically analyzed. Boundary problem has been formulated in dimensionless variables such as "stream function - vorticity vector - temperature" and it has been solved by finite difference method. The effect levels of the heat source, the medium permeability, the transient factor and the heat conductivity of the solid walls on flow patterns and heat transfer modes have been determined.

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

The paper deals with the statement of a problem of dynamic interaction of arbitrary solid bodies and its test solutions in the context of discrete element modeling. For discretization we use description of bodies with arbitrary shapes, composed of rigidly bound spheres. The clumps were built with different characteristics, which allowed to estimate their influence on the process of clump construction and the smoothness of obtained surface. A system of equations of motion relative to global axes for a clump of spheres is presented. The forces of interaction between the spheres are determined based on the Hertz-Mindlin contact model with due account for viscous damping. A problem of interaction of two spheres was chosen as a test case. Spheres' trajectories composed of clumps of spheres were calculated. The results were compared with the results for the case of motion and interaction of spheres in one-particle approximation.

Рассмотрена математическая модель дозвуковых нестационарных турбулентных течений несжимаемого газа, основанная на методе крупных вихрей. Приводятся описания модели подсеточной турбулентности и вычислительного алгоритма, представлены результаты параметрических расчетов турбулентных течений несжимаемого газа в прямоугольном канале при различных числах Рейнольдса.

The Large Eddy Simulation mathematical model for the subsonic unsteady turbulent flow of incompressible gas has been considered. The subgrid turbulence model, computational algorithm and numerical results of parametrical study of incompressible gas flow of in square channel under various Reynolds numbers has been presented.

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

Mathematical simulation of convective-conductive-radiative heat transfer in a cubical cavity filled with diathermanous medium has been carried out. The domain of interest is a closed volume having heat-conducting solid walls of finite thickness with diffuse grey inner surfaces. The outer surfaces of two vertical walls are isothermal while the other walls are adiabatic. The boundary-value problem has been formulated in dimensionless variables such as “vector potential-vorticity vector-temperature’’ in the Boussinesq approximation and taking into account the diathermancy of the continuous medium. An analysis of surface thermal radiation has been conducted on the basis of the net-radiation method in the form of Poljak. The formulated transient boundary-value problem has been solved by finite difference method in a wide range of the Rayleigh number, thermal conductivity ratio and surface emissivity. Correlations for the average convective and radiative Nusselt numbers at the characteristic internal solid-fluid interface have been obtained. The comparison between the obtained three-dimensional results and the two-dimensional data has been conducted. It has been found, that on the basis of a three-dimensional model it is possible to analyze the formation of intensive transverse flows from two vertical surfaces which are absent in a two-dimensional model. It has been also shown, that the solution of convective-radiative heat transfer problems in the conjugate statement leads to essential changes in distributions of local and integral parameters in comparison with the non-conjugate model, which first of all is related to a more correct description of the thermal radiation in diathermanous media due to taking into account the thermal conduction of the solid walls.

Численно исследуется явление возникновения скоростного бафтинга при обтекании профиля 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.

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

We consider the infinite non-Hermitian finite-difference Kitaev model simulating a one-dimensional superconducting wire. Non-Hermitianity is introduced into the model using delta-shaped imaginary potentials that simulate the gains and losses of the amplitudes of Majorana bound states (MBS). In a rigorous mathematical approach, the conditions for the existence of eigenfunctions describing the MBSs are found, as well as the dependence of the eigenfunctions on the model parameters and the effect of non-Hermiticity on the MBSs. Two regimes are considered, near the topological interphase boundary and at zero chemical potential.

В работе представлены результаты расчетного исследования локальной структуры восходящего газожидкостного потока в вертикальной трубе. Математическая модель основана на использовании двухжидкостного эйлерова подхода с учетом обратного влияния пузырьков на осредненные характеристики и турбулентность несущей фазы. Турбулентная кинетическая энергия жидкости рассчитывается с применением двухпараметрической изотропной модели турбулентности $k - \varepsilon$, модифицированной для двухфазных сред. Для описания динамики распределения пузырьков по размерам используются уравнения сохранения количества частиц для отдельных групп пузырьков с различными диаметрами для каждой фракции с учетом процессов дробления и коалесценции. Численно исследовано влияние изменения степени дисперсности газовой фазы, объемного расходного газосодержания, скорости дисперсной фазы на локальную структуру и поверхностное трение в двухфазном потоке. Сравнение результатов моделирования с экспериментальными данными показало, что разработанный подход позволяет адекватно описывать турбулентные газожидкостные течения в широком диапазоне изменения газосодержания и начальных размеров пузырьков.

The results of numerical simulation of the structure of a two-phase flow of a gas-liquid bubble mixture in a vertical ascending flow in a pipe are presented. The mathematical model is based on the use of the two-fluid Eulerian approach taking into account the inverse influence of bubbles on averaged characteristics and turbulence of the carrying phase. The turbulent kinetic energy of a liquid is calculated using equations for the transfer of Reynolds stresses. To describe the dynamics of bubble size distribution, the equations of particle number conservation for individual groups of bubbles with different constant diameters for each fraction are used taking into account the processes of breakup and coalescence. The influence of changes in the degree of dispersion of the gas phase, volume flow gas content and the velocity of the dispersed phase on the local structure and surface friction in the two-phase flow is numerically investigated. Comparison of simulation results with experimental data has shown that the developed approach allows an adequate description of turbulent gas-liquid flows in a wide range of changes in gas content and initial bubble sizes.

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

Mathematical models were developed and the nonlinear boundary value problem of dynamics thinwalled shells of the arbitrary form under action shock pulse is formulated. Dependence of processes of deformation on speed loading, compressibility of a material, finite deformations and large displacements of a shell middle surface, formation and kinetic of plasticity zones of a material during action of a shock wave are considered. Parameterization of a shell surface is carried out by bi-cubic splines. For the description of nonlinear, time and speed dependents of a shell material behavior with anisotropic hardening the generalized model of microplasticity is developed on the account of viscosity of deformation, hysteresis losses and Baushinger's effect. The solution of boundary value problems on the basis of difference schemes is constructed. Results of modeling of nonlinear wave processes in a assemble shell under action of explosion also are presented.

Численно исследуются газодинамические процессы, протекающие в начальный момент работы сверхзвукового сопла с высокой степенью геометрического расширения. Основное внимание уделяется изучению механизмов потери течением осевой симметрии за счет неустойчивости образующихся в сверхзвуковой части сопла зон отрывного течения. Модель нестационарного течения вязкого теплопроводного сжимаемого газа по соплу основана на системе уравнений сохранения в форме Навье-Стокса. Турбулентность исследуемого течения моделируется методом отсоединенных вихрей 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.

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

The paper presents simulation of heat and mass transfer gas (air) bubble in a liquid (water) cooling the heated metal surface.