ИССЛЕДОВАНИЕ ВЛИЯНИЯ ИСХОДНЫХ ЗНАЧЕНИЙ СКОРОСТИ НА ПРОЦЕССЫ ПЕРЕМЕШИВАНИЯ И РАСПРОСТРАНЕНИЯ ТРЕХМЕРНЫХ ТУРБУЛЕНТНЫХ СТРУЙ Текст научной статьи по специальности «Физика»

¿kuiSliVERSUM:

№Щ80)_ТЕХНИЧЕСКИЕ НАУКИ_ноябрь. 2020 г.

ИССЛЕДОВАНИЕ ВЛИЯНИЯ ИСХОДНЫХ ЗНАЧЕНИЙ СКОРОСТИ НА ПРОЦЕССЫ ПЕРЕМЕШИВАНИЯ И РАСПРОСТРАНЕНИЯ ТРЕХМЕРНЫХ ТУРБУЛЕНТНЫХ СТРУЙ

Махмудов Содикжон Ахмаджонович

канд. техн. наук, доцент, Андижанский машиностроительный институт, Республика Узбекистан, г. Андижан

Эшонхужаев Дилмурод Одилович

инженер-исследователь, Андижанский машиностроительный институт, Республика Узбекистан, г. Андижан E-mail: eshonxojayevdilmurod34@,gmail.com

RESEARCH OF THE INFLUENCE OF THE INITIAL VALUES OF THE VELOCITY ON THE PROCESSES OF MIXING AND PROPAGATION OF THREE-DIMENSIONAL TURBULENT JETS

Sodikzhon Makhmudov

Cand. tech. Sciences, Assoc. prof., Andijan Machine-Building Institute, Uzbekistan, Andijan

Dilmurod Eshonkhuzhaev

Research Engineer, Andijan Machine-Building Institute, Uzbekistan, Andijan

АННОТАЦИЯ

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

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

ABSTRACT

This scientific article investigates the influence of initial values on three-dimensional isothermal and non-isothermal jet flows flowing from a rectangular device on the right side.

In the case of turbulent exchange in the boundary layer of three-dimensional jet flows, a generalization of L. Prandtl's model taking into account temperature variability is proposed and the values of the constants in it are determined.

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

Keywords: spatial coordinates turbulent Prandtl number, dynamic coefficient of viscosity, total enthalpy of a gas mixture, temperature, density of a gas mixture, kinematic coefficient of turbulent viscosity.

Введение. В целях исполнения В данном работе рассматривается задача истечения струи из квадратного сопла и распространения ее в спутном потоке воздуха.

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

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

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

Библиографическое описание: Махмудов С.А., Эшонхужаев Д.О. Исследование влияния исходных значений скорости на процессы перемешивания и распространения трехмерных турбулентных струй // Universum: технические науки : электрон. научн. журн. 2020. 11(80). URL: https://7universum.com/ru/tech/archive/item/10940 (дата обращения: 25.11.2020).

№ 11(80)

UNIVERSUM:

ТЕХНИЧЕСКИЕ НАУКИ

ноябрь, 2020 г,

Относительно коэффициента турбулентной вязкости использована эффективная вязкость с учетом трёх мерности задачи и сжимаемости объема:

Уэфф = Vл + Ра(Ь2(Ыу) + ^)ф2 + ф2 + &фа (1)

Вышеизложенная задача описывается системой нелинейных уравнений.

Уравнения неразрывности:

дуй ^ I dpv ^ dpw _ g дх L ду dz

Уравнения движения по координате Х

(2)

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

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

I. х = Ъ1) 0<у<а,0<г<Ь

и = и2,у = 0 , ш = 0 ,р = р2 ,Н = Н2,

2) а<у < у+т , Ь <г <г+т

ди

ди

ди

pu — + pv—~ + pw— = dx Ldy dz

дх + Ь2 ду ду) + дг дг) Уравнения движения по координате У ди ди ди дР 1 д / ди

(3)

ÖP I о / ои\ о / ои\ 2 д / dv\ д ( dv\

YLdy{ßd-z) + YAßdy) (4)

dy\ dz/ dz\ dy, Уравнения движения по координате Z

dw

dw

dw

dP 4 д /dw\

+

---1---

dz 3 dz I d fdw\

l?dy[ß~dy) +

dw dz

I d idv\ 2 d idv\

Уравнения состояние газовой смеси Р=рТ

Связь полной энтальпии с температурой

?. ?. ?

H = cvT +

(5)

(6)

(7)

соответственно.

Метод и результаты исследования. Для решения системы уравнений (1)-(2) с граничными условиями их необходимо привести к безразмерному виду. Использования безразмерной формы записи исходных уравнений и граничных условий позволяет существенно снизить степен конкретизации решаемой задачи. С другой стороны, важной с точки

u = u1tv = 0 ,w = 0 ,р=р1 ,Н = Н1 , (8) II. х>0: 1) z = 0,0 <y < y+m

du dv dH

— = 0,— = 0,w = 0 ,— = 0 dz dz dz

2) y = 0,0 < z < z+m

du dw dH

— =0 ,v = 0 ,— =0 ,— = 0 dy dy dy

3) z-z+OT , y — y+t

и = u1 ,v = 0 ,w = 0 ,H = H1

(9)

В настоящей работе приводится численное исследование, имеющее целью, во-первых: проверить разработанные алгоритм и методику решения, приведенные в этой работе реализованные в виде программы для ПЕВМ на языке Фортран, которые позволяют получить численные результаты в широком диапазоне исходных параметров, во вторых выбрать модели первых моментов относительно коэффициента турбулентного обмена с учетом трёхмерности объемной сжимаемости и определения эмпирических констант, входящих в эти модели.

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

ментальных работ и2 = 38т/с иг 300 к

Рл = Р? =Р„+п, = 1 ат

= 0 ТЛ = Т2 =

Заключение. Основные результаты расчетов приведены в виде графиков на рисунках (1) - (7). На рис. 1а, б приведено поперечное распределение продольной скорости по осям оу и о2 в зависимости от

9

№ 11(80)

UNIVERSUM:

ТЕХНИЧЕСКИЕ НАУКИ

ноябрь, 2020 г,

скорости струйного потока. Пунктирная линия соответствует затопленной струе, а оплошная линия профилю продольной скорости при иг=5 т/с . Присутствие скорости струйного потока приводит к медленному затуханию скорости струи, т.е. при иг=5 т/с дальнобойность струи увеличивается примерно на 12 -13 %, а при скорости ступного потока и1=10 т/с дальнобойность струи увеличивается примерно на 16 -18 % по сравнению щ = 0 (рис 2 и 3).

На рис. 4 приведены поперечные распределения продольной скорости по оси оу вдоль струи в зависимости от исходного значения скорости основного потока. Здесь также можно отметить, что увеличение исходного значения скорости основного потока приводить к дальнобойности струи. Например, при исходном значении скорости основного потока и2 = 50т/с против и2 = 38т/с дальнобойность струи увеличивается примерно на 7-7,51% , что соответствует физике процесса. На рис. 5 приведены границы зоны смещения в плоскости оу2 в разных сечениях вдоль струи. Здесь можно отметить, что в начальных сечениях расширение границы зоны смешения по оси 02 отстает по сравнению с осью оу примерно до 5 калибров от входного сечения при истечении струи из квадратного сопла и границы зоны смешения принимают при этом эллипсоидную формула при х=15 калибров и больше струя переходить к круглой форме и к осесимметричным струям.

На рис. 6,7 приведены границы зоны смешения струи по осям оу и о2. Из рисунков видно, что при распространении струи в затопленном потоке воздуха граница зоны смешения шире, чем при скорости спутного потока иг=5 т/с . Здесь также можно отметить, что увеличение исходного значения скорости спутного потока приводит к сужению границы зоны смешения.

Рисунок 2. Поперечное распределение продольной скорости в разных сечениях по осям OY и OZ

Рисунок 3. Поперечное распределение продольной скорости при разных спутностях

Рисунок 4. Поперечное распределение продольной скорости по осям OY в разных сечениях

Рисунок 1. Поперечное распределение продольной скорости в зависимости от спутности по осям OYи OZ

Рисунок 5. Границы зоны смешения струй в разных сечениях

Л.

UNIVERSUM:

ТЕХНИЧЕСКИЕ НАУКИ

ноябрь, 2020 г,

Рисунок 6. Границы зоны смешения

Рисунок 7. Границы зоны смешения при разных ступностях

Список литературы.

1. Абрамович Г.Н. Теория турбулентных струй. Москва. Наука, 1984 г.

2. Алиев Ф., Жумаев З.Ш., Струйные течения реагирующих газов. Ташкент. Фан 1987 г.

3. Бай-Ши-И. Теория струй. Москва. Физматгиз. 1960 г.

4. Мак - Гирик Дж. Дж., Роди В. Расчет трехмерных турбулентных струй. Москва. Машиностроения 1984 г.

5. Палатник И.Б., Темирбаев Д.Ж. О распространении свободных турбулентных струй вытекающих из насадка прямоугольной формы. «Проблемы теплоэнергетики и прикладной теплофизики». Алма ата. 1964 г.