Науки о Земле
УДК 51-7, 533,62-1, 550
DOI: 10.21209/2227-9245-2020-26-3-6-14
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ ПЛОТНОСТИ МАТЕРИАЛА ДИСПЕРСНЫХ ВКЛЮЧЕНИЙ НА ПРОЦЕСС РАЗЛЕТА ГАЗОВЗВЕСИ В ВАКУУМ
NUMERICAL STUDY OF THE INFLUENCE OF MATERIAL DENSITY OF DISPERSED INCLUSIONS ON THE PROCESS OF GAS-SUSPENSION SPLITTING IN VACUUM
Математически моделируется течение однородного газа и неоднородной среды, представляющей собой взвесь твердых частиц в газе - газовзвесь. Целью работы является изучение влияния плотности материала частиц твердой компоненты смеси на процесс истечения смеси в вакуум и выявление отличий от процесса истечения в вакуум однородного газа. При моделировании процесса истечения учитывались вязкость, сжимаемость и теплопроводность газа. Математическая модель, примененная в исследовании, реализует континуальную методологию моделирования течения неоднородной среды.
Особенностью этой методики является то, что при описании движения смеси записывается полная гидродинамическая система уравнений движения для каждой из компонент смеси. При этом системы уравнений движения компонент смеси связаны слагаемыми, отвечающими за межфазное силовое и тепловое взаимодействие.
Важность учета межкомпонентного взаимодействия в смеси связана с тем, что динамика гетерогенных сред - смесей, в которых компоненты имеют различное агрегатное состояние, во многом определяется эффектами, связанными с межкомпонентным взаимодействием. Система уравнений математической модели включает уравнения непрерывности для плотности несущей среды и «средней плотности» дисперсной компоненты смеси.
Для описания сохранения импульса несущей среды решалось уравнение Навье - Стокса, для дисперсной компоненты смеси также записывалось уравнение сохранения импульса с учетом слагаемых, отвечающих за межкомпонентное взаимодействие. Уравнения сохранения энергии компонент смеси решались с учетом межкомпонентного теплообмена. Система уравнений математической модели, дополненная краевыми условиями, решалась явным конечно-разностным методом второго порядка точности.
В результате моделирования выявлены отличия в распределении параметров сплошной среды при распространении в вакуум чистого газа и газовой взвеси частиц. Выявлено влияние плотности материала частиц на процесс истечения несущей среды и дисперсной компоненты газовзвеси в вакуум
Ключевые слова: газовзвесь; уравнение Навье - Стокса; численное моделирование; уравнения в частных производных; краевая задача; моделирование технологических процессов; межкомпонентное взаимодействие; газовая динамика; скоростные потоки газа; истечение в вакуум
In the article, the flow of a homogeneous gas and inhomogeneous medium representing a suspension of solid particles in a gas - gas suspension is mathematically modeled. The aim of the work is to study the effect of material density of the solid component particles of the mixture on the process of the mixture outflow into vacuum and to identify differences from the process of a homogeneous gas outflow into vacuum. When simulating the outflow process, the viscosity, compressibility and thermal conductivity of the gas were taken into account. The
А. А. Ахунов,
Казанский национальный исследовательский технический университет, г. Казань Ahunov.Adel@yandex.ru
Д. А. Тукмаков,
Казанский научный центр РАН, г. Казань
tukmakovDA@imm.knc.ru
А. Akhunov, Kazan National Research Technical University, Kazan
D. Tukmakov,
Kazan Scientific Center RAS
© А.А. Ахунов, Д. А. Тукмаков, 2020
6
mathematical model described in this article implements a continuous methodology for modeling the inhomo-geneous medium flow.
A feature of this technique is that when describing the motion of a mixture, a complete hydrodynamic system of motion equations is recorded for each of the components of the mixture. In this case, the systems of motion equations of the mixture components are connected with the components responsible for interphase force and thermal interaction.
The importance of taking into account the intercomponent interaction in a mixture is due to the fact that the dynamics of heterogeneous media, mixtures in which the components have different aggregate states, are largely determined by the effects associated with intercomponent interaction. The system of the mathematical model equations includes continuity equations for the density of the carrier medium and the "average density" of the dispersed component of the mixture.
To describe the momentum conservation of the carrier medium, the Navier - Stokes equation has been solved; for the dispersed component of the mixture; the equation of momentum conservation has also been written, taking into account the terms responsible for the intercomponent interaction. The energy conservation equations for the mixture components have been solved, taking into account inter-component heat transfer. The system of the mathematical model equations, supplemented by boundary conditions, has been solved by an explicit finite-difference method of the second order of accuracy.
As a result of the simulation, differences in the distribution of the parameters of a continuous medium during the propagation of pure gas and gas suspension of particles into a vacuum are revealed as well as the influence of the material density of the particles on the process of the carrier medium and dispersed component expiration of the gas suspension in vacuum
Key words: gas suspension; Navier - Stokes equation; numerical simulation; partial differential equations; boundary value problem; modeling of technological processes; intercomponent interaction; gas dynamics; high-speed gas flows; outflow into vacuum
Б ведение. Течения неоднородных сред встречаются как в естественной природе, так и в технике. В ряде монографий описана методология моделирования течений неоднородных сред [9; 11; 16]. Теоретическому и экспериментальному исследованию динамики неоднородных сред посвящены публикации в периодических изданиях [1-8; 10; 12; 14-16; 18-20].
Во многом интерес к исследованию таких сред связан с оптимизацией работы реактивных двигателей и экранированием промышленных взрывов взвесями жидких капель либо твердых частиц [1; 3]. Изучение динамических процессов в неоднородных средах связано с проблемами сварочной техники, экологии и здравоохранения [4; 6-8; 10]. Течения неоднородных сред встречаются также в технологиях обогащения полезных ископаемых [2].
Истечение неоднородной среды в разряженный газ происходит при работе газозаборных аппаратов. В статье «Аспирация аэрозоля в цилиндрический пробоотборник из низкоскоростного нисходящего потока и из неподвижной среды» аналитически смоделирован процесс аспирации аэрозоля в двумерной области, без учета взаимообратного влияния несущей и дисперсной среды.
Предполагалось, что динамика дисперсной компоненты не оказывает влияния на движение несущей среды [5]. Анализ теоретических и экспериментальных исследований динамики аэрозолей указывает, что процессы межфазного взаимодействия приводят к существенным отличиям характеристик течений однородных сред и дисперсных смесей [9].
Объектом исследования являются газовзвеси, или аэрозоли - неоднородные среды, состоящие из взвешенных в газе частиц. В случае, если дисперсная компонента газовзвеси является жидкостью, такие среды называют газокапельными, если дисперсная компонента газовзвеси состоит из твердых частиц - запыленными.
Предметом исследования является один из видов течений газовзвеси, так называемое «истечение в вакуум». Для невязкого газа существуют представленные в научной литературе аналитические решения. Однако реальные газы являются вязкими средами, в связи с чем при моделировании возникает необходимость учесть вязкость среды. В данной работе вязкость несущей среды учитывается.
Целью исследования является определение того, каким образом плотность дис-
персной компоненты среды влияет на истечение химически и термодинамически (без фазовых переходов) инертной газовзвеси, при этом, как правило, инертные газовзвеси являются запыленными средами, т. е. взвесями твердых частиц. Необходимость подобных теоретических исследований для промышленности связана с тем, что наличие дисперсных включений оказывает значительное влияние на характер процессов и рабочие характеристики устройств, в которых используются потоки неоднородных сред [9; 16].
Задача, которая решалась в процессе исследования, представляла собой краевую задачу с начальными значениями для нелинейной системы уравнений в частных производных. Решение поставленной задачи осуществлено методом конечных разностей.
Методология и методы исследования. Движение среды описывалось двумя системами уравнений - системой уравнений движения несущей фазы и системой уравнения движения дисперсной компоненты. Система уравнений движения несущей фазы включала три уравнения: уравнение непрерывности плотности газа (1); уравнение Навье - Сток-са [13; 17] в одномерном приближении (2); уравнение непрерывности полной энергии газа (3). В совокупности уравнения (1)-(3) образуют полную гидродинамическую систему уравнений движения сплошной среды в нестационарной постановке для течения с одномерной геометрией. Уравнение (4) описывает функцию сдвиговых вязких напряжений газа в одномерном течении.
Ф+<зИ=0 8( & '
д(р») , 8
ь
81 дх
гкъ
&
ас /
4 8и т = —и —. 3 8х
(1)
(2)
(3)
(4)
Функция вязких напряжений является выражением, описывающим вязкие потери несущей среды, выражающиеся через пространственные производные от координатных составляющих векторы скорости, в одномерном случае скорость газа является скалярной величиной. Функция (4) входит как в уравнение сохранения импульса (2) (Навье - Стокса), так и в уравнение сохранения энергии газовой фазы смеси (3). В общем
случае функция является нелинейной, однако в данной работе рассматривается течение газа, для которого названная функция является линейной.
Сжимаемость несущей среды учитывалась представлением плотности газа как функции зависящей от пространственной и временной координат. Динамика ударно-волновых, детонационных, сверхзвуковых и, в целом, высокоскоростных течений газа, включая и истечения газа в вакуум, требует математических моделей, учитывающих сжимаемость среды. Отсутствие сжимаемости течения среды в одномерном случае привело бы к сведению уравнения нераз-
и\ ды
рывности массы газа (1) к уравнению — = 0,
дх
из чего следовало бы, что функция скорости является стационарной по пространственной координате функцией временной переменной п=п(,)
Динамика дисперсной фазы аэрозольной смеси описывается системой уравнений подобной системе уравнений (1)-(3), за рядом исключений. Прежде всего, для дисперсной компоненты решается уравнение сохранения «средней плотности» (5), а не физической плотности материала дисперсных включений [9]. «Средняя плотность» - произведение истинной (физической) плотности материала частиц и объемного содержания дисперсной фазы, являющегося функцией пространственной и временной переменных. Таким образом, при математическом моделировании дисперсная компонента также имеет некий аналог «сжимаемости», соответствующей сжимаемости несущей среды - газа. При этом сами частицы являются несжимаемыми.
Также система уравнений динамики дисперсной компоненты включает уравнение сохранения импульса (6), отличающееся от уравнения Навье - Стокса (2) отсутствием функции вязких напряжений. Для дисперсной компоненты смеси записывается уравнение сохранения энергии (7). Уравнения математической модели (2), (3), (6), (7) в правых частях имеют слагаемые, отвечающие за обмен импульсом (в уравнениях сохранения импульса) и энергией (в уравнениях сохранения энергии) между компонентами смеси.
Ф, , 5(Р1^1> _ 0
дх
d(Pi"i) , 8
8t 8xK 1 i; 8x
5(ei) 8 i \ ^ аГ+ъЬ*^-
А =аРы. ei =PicPTv
„ J a „ , I / (8u du
Fs=-4щСМи~щ\{и-щ), FA=ap^- + u-
, (8u 8u 8ul 8ul .
F„„ =0.5ap.\ — + u---L-u,—L ,
a™ I 8t 8x dt 1 8x 1
(6) (7)
ei=PiCPTP
(9)
F = Fs + FA+Fim,Cd
24 4
Re, Re,
+ 0.4,
М1 =|и -щ\/с ,
Re1=p\u-u1\2r/\1, Рг = уСр\1/Х .
В уравнениях (1), (5) правых частей нет, что связано с отсутствием массообменных процессов между компонентами, так называемой инертностью газовзвеси (отсутствием химических превращений и фазовых переходов). Здесь а = а(Г,х) - функция, описывающая объемное содержание дисперсной фазы. Коэффициент сопротивления Сб - соответствует коэффициенту аэродинамического сопротивления частиц сферической формы.
В уравнениях использованы следующие физические величины: р - давление, р -плотность, и - скорость; Т - температура, е -энергия. Индекс «1» обозначает физическую величину дисперсной компоненты смеси. Физическая величина без индекса описывает изменение параметров газовой фазы.
Температура газа находится из уравнения
T=(y-l)(e/p-0,5u)/R,
(8)
где Я - газовая постоянная.
Силовое взаимодействие несущей и дисперсной фаз учитывает несколько разнородных сил, действующих на дисперсные частицы со стороны газа: силу Стокса FS, динамическую силу Архимеда FA и силу присоединенных масс Fаm. При этом силовое воздействие имеет взаимообратный характер [9]. Математическая модель предполагает монодисперсный состав твердой компоненты газовзвеси - все частицы имеют одинаковый размер и одинаковые физические свойства - плотность и теплоемкость материала.
Тепловая энергия взвешенной в газе дисперсной фазы определяется как
где Ср - удельная теплоемкость единицы массы вещества, из которого состоят частицы.
Тепловой поток между компонентами смеси описывается следующим выражением:
^=6раМи1Ъ(Т-Т1)/(2г)2. (10)
В математической модели использованы относительные числа Нуссельта, Маха, Рейнольдса, а также число Прандтля [9]
0.55 п 033
Щ = 2ехр(-М1) + 0.459 Re, Рг
(11)
Для искомых функций задавались следующие граничные условия [5; 17; 20]:
м(л1) = 0,м1(^1) = 0, е(г,\) = е(г,2), е1(^1) = е1(л2), = р^,2),р1 (М) = А (^.2), и (/, ы) = и (г, N -1), щ {г, ы) = щ N -1), е (*, ЛГ) = е (*, ЛГ -1), е, (Л = е, N -1), = 1), =
В начальный момент времени компоненты смеси покоились:
и(0,0 = 0, и1(0,г) = 0.
Для всех искомых функций в моделируемой области течения задавались начальные значения в левой половине канала (х^/2)
е(= ¿о> = ею'
Р ( 0 г ) = А, . А ( 0 г ) = А о
и в правой половине канала (х^/2) е(0,/) = 0, ^(0,/) = 0,
р(0,/) = 0,р1(0,/) = 0.
Система уравнений течения двухком-понентной среды решалась явным конечно-разностным методом второго порядка точности по пространственным переменным [17]. Численное решение проводилось на равномерной сетке с количеством узлов вдоль оси х - N = 1000. Методология моделирования тестировалась сопоставлением с численными решениями, известными из ли-
тературы [18], и результатами физических экспериментов [12].
Основные уравнения динамики неоднородных сред разработаны в монографии «Основы механики гетерогенных сред» [11]. Тематика высокоскоростных течений газовзвесей (как правило, химически инертных запыленных сред) подробно рассмотрена в ряде статей и монографий [3; 9; 15; 16]. Вместе с тем истечение неоднородных сред в вакуум исследовано в меньшей мере, чем ударно-волновые и скоростные течения дисперсных сред. Также существенный интерес представляет рассмотрение взаимообратных эффектов, возникающих при истечении газовзвеси в вакуум, и то, каким образом параметры дисперсной среды влияют на динамику газа при истечении в вакуум двухкомпо-нентной смеси.
Результаты исследования. На рис. 1 представлено схематичное изображение моделируемой физической области, в которой одна половина канала заполнена газовзвесью, а в другой - вакуум.
Значения скорости при разлете в вакуум в невязком газе - 725 м/с - больше, чем в вязком газе (557 м/с), при этом в вязком газе с дисперсной компонентой скорость спутно-го следа (315 м/с) существенно ниже, чем в
Рис. 2. Пространственное распределение продольной скорости: 1 - в аналитическом решении;
2 - в численном решении для однородного вязкого газа;
3 - в численном решении для газовзвеси с плотностью
материала частицр10=2500 кг/м3 / Fig. 2. Spatial distribution of the longitudinal velocity:1 - in the analytical solution; 2 - in the numerical solution for a homogeneous viscous gas; 3 - in a numerical solution for a gas suspension with a particle material density of р10 = 2500 kg /m3
однородной среде (рис. 2). Численные расчеты для разлета чистого газа и газовзвеси в вакуум имеют отличия. Учет межкомпонентного взаимодействия между несущей средой и дисперсной компонентой приводит к тому, что результаты расчетов численной модели газовзвеси выявляют более низкие скорости спутного потока газа.
Рассмотрим влияние плотности материала частиц твердой компоненты смеси на процесс разлета газовзвеси в вакуум. При последовательном увеличении плотности материала частиц происходит уменьшение скорости, достигаемой спутным потоком газа (рис. 3). Также замедляется процесс истечения газовой компоненты смеси в вакуум.
. » « о *e о
о • 0 * " e
4 V» «
Рис. 1. Схематичное изображение канала с запыленной средой / Fig. 1. Schematic representation of a channel with a dusty environment
Рис. 3. Численно рассчитанные пространственные распределения скорости газа при разлете в вакуум: 1 - вязкого газа с дисперсной компонентой с плотностью материала р10=750 кг/м3; 2 - вязкого газа с дисперсной компонентой с плотностью материала р10=2500 кг/м3; 3 - вязкого газа с плотностью материала р10=19050 кг/м3/ Fig. 3 Numerically calculated spatial distributions of the gas velocity during the expansion of a: 1 - viscous gas with a dispersed component with a material densityp10 = 750 kg / m3 into vacuum; 2 - viscous gas with a dispersed component with a material density p10 = 2500 kg / m3; 3 - viscous gas with a material densityp10 = 19050 kg / m3
Для частиц из менее плотных материалов распределение давления близко к распределению давления в течении однородного газа. С увеличением физической плотности дисперсной компоненты смеси происходит замедление течения газа (рис. 4). Интенсивность скоростного скольжения компонент смеси в процессе разлета газовзвеси в вакуум возрастает с увеличением плотности частиц (рис. 5). Это связано с наличием у более
плотных частиц большей массы и, следовательно, большей инерционности. Таким образом, интенсивность скоростного скольжения прямо пропорциональна массе частицы.
Результаты расчетов демонстрируют, что для более плотных материалов дисперсной фазы замедляется массоперенос дисперсной компоненты газовзвеси в вакуум (рис. 6).
\и-и I, м/с 1 250-
0.8 х, М
Рис. 4. Пространственное распределение давления
газа при разлете в вакуум: 1 - чистого газа; 2 - газовзвеси с плотностью материала дисперсной фазы р10=750 кг/м3; 3 - газовзвеси с плотностью материала дисперсной фазы р10=2500 кг/м3; 4 - газовзвеси с плотностью материала дисперсной фазы р10=19050 кг/м3 / Fig. 4. Spatial distribution of gas pressure during expansion into a vacuum: 1 - of pure gas; 2 - gas suspensions with the density of the dispersed phase
materialр10 = 750 kg / m3; 3 - gas suspensions with the density of the dispersed phase material р10 = 2500 kg / m3; 4 - gas suspensions with the density of the dispersed phase material р10 = 19050 kg / m3
Рис. 5. Пространственное распределение величины межфазного скоростного скольжения при разлете в вакуум: 1 - газовзвеси с плотностью материала дисперсной фазы р10=750 кг/м3; 2 - газовзвеси с плотностью материала дисперсной фазы р10=2500 кг/м3; 3 - газовзвеси с плотностью материала дисперсной фазы р10=19050 кг/м3/ Fig. 5. Spatial distribution of the magnitude of interfacial velocity slip during expansion into vacuum of a: 1 - gas suspension with a material density of the dispersed phase р10 = 750 kg / m3; 2 - gas suspensions with the density of the dispersed phase material
р10 = 2500 kg / m3; 3 - gas suspension with the density of the dispersed phase material р10 = 19050 kg / m3
Рис. 6. Пространственное распределение объемного содержания дисперсной фазы газовзвеси при разлете в вакуум: 1 - газовзвеси с плотностью материала дисперсной фазы р10=750 кг/м3; 2 - газовзвеси с плотностью материала дисперсной фазы р10=2500 кг/м3; 3 - газовзвеси с плотностью материала дисперсной фазы р10=19050 кг/м3 / Fig. 6. The spatial distribution of the volumetric content of
the dispersed phase of the gas suspension during expansion into vacuum of a: 1 - gas suspension with the density of the material of the dispersed phase р10 = 750 kg / m3; 2 - gas suspensions with the density of the dispersed phase materialр10 = 2500 kg / m3; 3 - gas suspension with the density of the dispersed phase material р10 = 19050 kg / m3
Выявленные закономерности обусловлены тем, что более плотные частицы за счет большей инерционности имеют меньшую скорость в потоке двухкомпонентной среды. При этом интенсивность межкомпонентного взаимодействия определяется суммарной силой аэродинамического сопротивления частиц, а значит, возрастает одновременно с возрастанием интенсивности межкомпонентного скоростного скольжения, которое прямо пропорционально массе частицы. Более массивные частицы разгоняются до меньшей скорости, при этом среда имеет большее суммарное межкомпонентное силовое взаимодействие. Таким образом, частицы более плотных материалов при одинаковом размере оказывают большее влияние на динамику смеси и, в конечном счете, сильнее замедляют течение несущей среды - газа. Выявленные закономерности показывают, что в истечении газовзвеси в вакуум имеют место взаимообратные эффекты, вызванные взаимодействием между несущей и дисперсной фазами смеси.
Заключение. В результате численного моделирования разлета газовзвеси в вакуум обнаружено, что процесс истечения неоднородной среды имеет существенные различия
с истечением в вакуум чистого газа. Численное моделирование выявило следующие закономерности:
1) увеличение плотности материала частиц дисперсной компоненты газовзвеси приводит к уменьшению скорости спутного потока газа. Процесс истечения газовой и дисперсной компонент смеси в газовзвесях с более массивными частицами имеет большую длительность;
2) для газовзвесей, дисперсная компонента которых состоит из легких частиц, истечение неоднородной среды в вакуум имеет существенное сходство с процессом истечения в вакуум чистого однородного газа;
3) интенсивность межфазного скоростного скольжения в процессе разлета газовзвеси в вакуум прямо пропорциональна плотности материала дисперсных включений.
Анализ численных расчетов показывает, что замедление скорости истечения газовзвеси в вакуум, в частности газовой компоненты смеси, происходит при увеличении физической плотности твердой компоненты дисперсной среды. Результаты исследований возможно применить при разработке технологических процессов, связанных с течениями дисперсных сред.
Список литературы _
1. Арефьев К. Ю., Воронецкий А. В., Сучков С. А. Расчетное исследование особенностей дробления и испарения капель в газодинамических течениях с циклическими ударными волнами // Известия высших учебных заведений. Машиностроение. 2015. № 10. С. 17-30.
2. Баландин О. А., Верхотуров А. Р. Влияние интенсивности электромагнитных волн на движение твердой частицы // Вестник Забайкальского государственного университета. 2013. № 11. С. 17-21.
3. Бедарев И. А., Федоров А. В. Структура и устойчивость ударной волны в газовзвеси с двумя давлениями // Вычислительные технологии. 2015. Т. 20, № 2. С. 3-19.
4. Борисов И. С., Романов С. А. Особенности моделирования процесса образования сварочного аэрозоля при сварке горно-шахтного оборудования // Экология и безопасность в техносфере: современные проблемы и пути решения: сб. ст.: в 2 т. Т. 1. Томск: ТПУ, 2015. С. 186-190.
5. Ванюнина М. В., Галеев Р. С., Зарипов Ш. X., Скворцов Э. В. Аспирация аэрозоля в цилиндрический пробоотборник из низкоскоростного нисходящего потока и из неподвижной среды // Прикладная механика и техническая физика. 2005. Т. 46, № 2. С. 122-129.
6. Кириченко К. Ю., Дрозд В. А., Чайка В. В., Гридасов А. В., Голохваст К. С. Сварочный аэрозоль как источник опасных для здоровья техногенных нано- и микрочастиц: гранулометрический анализ // Известия Самарского научного центра Российской академии наук. 2015. Т. 17, № 5-2. С. 662-665.
7. Кондратьев И. И. Оценка скорости осаждения аэрозоля по материалам геохимических исследований и по данным мониторинга радиоактивного загрязнения атмосферы // Метеорология и гидрология. 2016. № 5. С. 78-87.
8. Кузнецова А. Н. Анализ сварочных аэрозолей // Фундаментальные проблемы науки: сб. ст. Тюмень, 2016. С. 41-46.
9. Кутушев А. Г. Математическое моделирование волновых процессов в аэродисперсных и порошкообразных средах. СПб.: Недра, 2003. 284 с.
10. Ляпунов H. А., Бовтенко В. А., Безуглая E. П., Столпер Ю. М. Обоснование нового подхода к оценке качества дозированных аэрозолей для ингаляций на этапе их разработки // Научные ведомости Белгородского государственного университета. Серия: Медицина. Фармация. 2016. Т. 33, № 5. С. 170-179.
11. Нигматулин Р. И. Основы механики гетерогенных сред. Москва: Наука, 1978. 336 с.
12. Нигматулин Р. И., Губайдуллин Д. А., Тукмаков Д. А. Ударно-волновой разлет газовзвесей // Доклады Академии наук. 2016. Т. 466, № 4. С. 418-421.
13. Тукмаков А. Л., Мубаракшин Б. Р., Тонконог В. Г. Моделирование процесса одоризации природного газа // Инженерно-физический журнал. 2016. Т. 89, № 1. С. 127-132.
14. Тукмаков Д. А. Численное моделирование ударно-волновых течений в газовзвеси с неоднородной концентрацией дисперсной фазы // Авиационная техника. 2019. № 1. С. 54-59.
15. Усманова Р. Р., Заиков Г. E. Постановка разностной краевой задачи к расчету параметров вихревого течения газовзвеси // Вестник Казанского технологического университета. 2015. Т. 18, № 1. С. 350-354.
16. Федоров А. В., Фомин В. М., Хмель Т. А. Волновые процессы в газовзвесях частиц металлов. Новосибирск: Параллель, 2015. 305 с.
17. Fletcher C. A. Computational techniques for fluid dynamics. Berlin: Springer-Verlang, 1988. 409 p.
18. Gubaidullin D. A., Tukmakov D. A. Numerical research of evolution of the shock wave in gas-particles suspension with account uneven distribution of particles // Mathematical Models and Computer Simulations. 2015. Vol. 7, No. 3. P. 246-253.
19. Tukmakov D. A. Numerical study of polydisperse aerosol dynamics with the drops destruction // Lobachevskii Journal of Mathematics. 2019. Vol. 40, No. 6. P. 824-827.
20. Zhuoqing A., Jesse Z. Correlating the apparent viscosity with gas-solid suspension flow in straight pipelines // Powder Technology. 2019. Vol. 345. P. 346-351.
References _
1. Arefiev K. Y, Voronetsky A. B., Suchkov S. A. Izvestiya vysshih uchebnyh zavedeniy. Mashinostroenie (News of higher educational institutions. Engineering), 2015, no. 10, pp. 17-30.
2. Balandin O. A., Verkhoturov A. R. VestnikZabaykalskogo gosudarstvennogo universiteta (Transbaikal State University Journal), 2013, no. 11, pp. 17-21.
3. Bedarev I. A., Fedorov A. V. Vychislitelnye tekhnologii (Computing technology), 2015, vol. 20, no. 2, pp. 3-19.
4. Borisov I. S., Romanov S. A. Ekologiya i bezopasnost v tehnosfere: sovremennye problemy iputi resh-eniya: sb. st.: v 2 t. T. 1. (Ecology and safety in the technosphere: modern problems and solutions: collected articles). Tomsk: TPU, 2015, pp. 186-190.
5. Vanyunina M. V., Galeev R. S., Zaripov Sh. Kh., Skvortsov E. V. Prikladnaya mekhanika i tehnicheskaya fizika (Journal of Applied Mechanics and Technical Physics), 2005, vol. 46, no. 2, pp. 122-129.
6. Kirichenko K. Yu., Drozd V. A., Chayka V. V., Gridasov A. V., Golokhvast K. S. Izvestiya Samarskogo nauch-nogo tsentra Rossiyskoy akademii nauk (Bulletin of the Samara Scientific Center of the Russian Academy of Sciences), 2015, vol. 17, no. 5-2, pp. 662-665.
7. Kondratiev I. I. Meteorologiya igidrologiya (Russian Meteorology and Hydrology), 2016, no. 5, pp. 78-87.
8. Kuznetsova A. N. Fundamentalnye problemynauki: sb. st. (Fundamental problems of science: collected articles), Tyumen, 2016, pp. 41-46.
9. Kutushev A. G. Matematicheskoe modelirovanie volnovyh protsessov v aerodispersnyh i poroshkoo-braznyh sredah (Mathematical modeling of wave processes in aerodispersed and powdery media). St. Peterburg: Nedra, 2003. 284 p.
10. Lyapunov N. A., Bovtenko V. A., Bezuglaya E. P., Stolper Yu. M. Nauchnye vedomosti Belgorodskogo gosudarstvennogo universiteta. Seriya: Meditsina. Farmaciya (Scientific reports of the Belgorod State University. Series: Medicine. Pharmacy), 2016, vol. 33, no. 5, pp. 170-179.
11. Nigmatulin R. I. Osnovy mehaniki geterogennyh sred (Fundamentals of the mechanics of heterogeneous media). Moscow: Nauka, 1978. 336 p.
12. Nigmatulin R. I., Gubaidullin D. A., Tukmakov D. A. Doklady Akademii nauk (Reports of the Academy of Sciences), 2016, vol. 466, no. 2, pp. 418-421.
13. Tukmakov A. L., Mubarakshin B. R., Tonkonog V. G. Inzhenerno-fizicheskiy zhurnal (Journal of Engineering Physics and Thermophysics), 2016, vol. 89, no. 1, pp. 127-132.
14. Tukmakov D. A. Aviatsionnaya tehnika (Aeronautical engineering), 2019, no.1, pp. 54-59.
15. Usmanova R. R., Zaikov G. E. VestnikKazanskogo tehnologicheskogo universiteta (Bulletin of the Kazan Technological University), 2015, vol. 18, no. 1, pp. 350-354.
16. Fedorov A. V., Fomin V. M., Khmel T. A. Volnovye protsessy vgazovzvesyah chastic metallov (Wave processes in gas-suspended particles of metals). Novosibirsk: Parallel, 2015. 301 p.
17. Fletcher C. A. Computation techniques for fluid dynacmics (Computation techniques for fluid dynacmics). Berlin: Springer-Verlang, 1988. 409 p.
18. Gubaidullin D. A., Tukmakov D. A. Mathematical Models and Computer Simulations (Mathematical Models and Computer Simulations), 2015, vol. 7, no. 3, pp. 246-253.
19. Tukmakov D. A. Lobachevskii Journal of Mathematics (Lobachevskii Journal of Mathematics), 2019, no. 6, pp. 824-827.
20. Zhuoqing A. Jesse Z. Powder Technology (Powder Technology), 2019, vol. 345, pp. 346-351.
Работа выполнена при поддержке гранта РФФИ № 19-01-00442
Коротко об авторах_
АхуновАдельАйратович, аспирант, кафедра теплотехники и энергетического машиностроения, Казанский национальный исследовательский технический университет, г. Казань, Россия. Область научных интересов: вакуумная и компрессорная техника Ahunov.Adel@yandex.ru
Тукмаков Дмитрий Алексеевич, канд. физ.-мат. наук, научный сотрудник, Институт механики и машиностроения, Казанский научный центр РАН, г. Казань, Россия. Область научных интересов: математическое моделирование, механика жидкости и газа tukmakovDA@imm.knc.ru
Briefly about the authors _
Akhunov Adel, postgraduate, Heat Engineering and Power Engineering department, Kazan National Research Technical University, Kazan, Russia. Sphere of scientific interests: vacuum and compressor technology
Tukmakov Dmitry, candidate of physical and mathematical sciences, researcher, Kazan Scientific Center of the Russian Academy of Sciences, Kazan, Russian Federation. Sphere of scientific interests: mathematical modeling, fluid and gas mechanics
Образец цитирования
Ахунов А. А., Тукмаков Д. А. Численное исследование влияния плотности материала дисперсных включений на процесс разлета газовзвеси в вакуум // Вестник Забайкальского государственного университета 2020. Т. 26, № 3. С. 6-14. DOI: 10.21209/2227-9245-2020-26-3-6-14.
Akhunov А., Tukmakov D. Numerical study of the influence of material density of dispersed inclusions on the process of gassuspension splitting in vacuum // Transbaikal State University Journal, 2020, vol. 26, no. 3, pp. 6-14. DOI: 10.21209/22279245-2020-26-3-6-14.
Статья поступила в редакцию: 13.03.2020 г. Статья принята к публикации: 23.03.2020 г.