Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,593

МОДЕЛЬ ТУРБУЛЕНТНОСТИ В ЖИДКОЙ ФАЗЕ БАРБОТАЖНОГО СЛОЯ

Лаптев А.Г. 1 Лаптева Е.А. 1
1 ФГБОУ ВПО «Казанский государственный энергетический университет»
Рассмотрено применение модели локальной изотропной турбулентности для определения коэффициентов турбулентной вязкости в газожидкостном слое при пузырьковом и пенном режимах. Получены выражения для расчета коэффициентов турбулентной вязкости, связанные с динамической скоростью трения и средней скоростью диссипации энергии. Представлены результаты расчета данных коэффициентов при различных скоростях газа в барботажном слое.
изотропная турбулентность
газожидкостный слой
турбулентная вязкость
динамическая скорость
диссипация энергии
1. Кутателадзе С.С. Консервативные свойства пристенной турбулентности / С.С. Кутателадзе // ТОХТ. – 1971. – Т. 5, № 1. – С. 3–12.
2. Таунсенд А.А. Структура турбулентного потока с поперечным сдвигом. – М.: Наука, 1959. – 400 с.
3. Барабаш В.М. Процессы переноса в турбулентных потоках с интенсивным внешним источником турбулизации / В.М. Барабаш // Теор. основы хим. технол. – 1994. – Т. 28, № 2. – С. 110–117.
4. Kawase Y. Mathematical models for design of biorectors applications of Kolmogoroff’s theory of isotropic turbulence / Y. Kawase, M. Moo-Young // Chem. Eng. J. – 1990. – V. 43, № 5. – P. 1319–1341.
5. Лаптев А.Г. Математическое моделирование теплоотдачи при турбулентном обтекании пучков труб / А.Г. Лаптев, С.Г. Дьяконов, В.И. Елизаров // Теплоэнергетика. – 1992. – № 12. – С. 34–38.
6. Брагинский Л.Н. Перемешивание в жидких средах: физические основы и инженерные методы расчета / Л.Н. Брагинский, В.И. Бегачев, В.М. Барабаш. – Л.: Химия, 1984.
7. Доманский И.В. Обобщение различных случаев конвективного теплообмена с помощью полуэмпирической теории турбулентного переноса / И.В. Доманский, В.Н. Соколов // Теор. основы хим. технол. – 1968. – Т. 2, № 5. – С. 761–767.
8. Соколов В.Н. Газожидкостные реакторы / В.Н.Соколов, И.В.Доманский. – Л.: Машиностроение, 1976.
9. Лаптева Е.А., Фарахов Т.М. Модель турбулентной вязкости в барботажных аппаратах // Вестник КГЭУ, 2012, № 4. С. 26-33
10. Лаптев А.Г. Модели пограничного слоя и расчет тепломассообменных процессов. Казань: Изд-во Казан. ун-та, 2007. – 500 с.
11. Нигматуллин Р.И. Динамика многофазных сред / Р.И. Нигматуллин. – М.: Наука, 1987.
12. Шлихтинг Г. Теория пограничного слоя. / Г. Шлихтинг. – М.: Наука, 1974.
13. Броунштейн Б.И. Гидродинамика, массо- и теплообмен в колонных аппаратах / Б.И. Броунштейн, В.В. Щеголев. – Л.: Химия, 1988.

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

Первые фундаментальные исследования турбулентных течений были выполнены еще в 19 веке Г. Гагеном, О. Рейнольдсом, Ж. Буссинеском, Г. Лоренцем и в начале 20 века развиты многими учеными (Л. Прандтль, Т. Кармен, В. Экман, И. Бюргерс, А. Колмогоров, Х. Драйден, Г. Клаузер, Дж. Тейлор, Г. Шлихтинг, Л. Лойцянский и др). Разработаны различные полуэмпирические теории пристенной турбулентности и отрывных течений. Осреднение по Рейнольдсу уравнений динамики вязкой жидкости приводит к системе уравнений Рейнольдса, называемые уравнениями для первых моментов. В общем случае проблема замыкания этих уравнений сводится к установлению связи тензора турбулентных напряжений с тензором осредненных скоростей деформаций.

Одна из первых попыток упростить расчет турбулентных течений принадлежит Буссинеску, который ввел понятие турбулентной вязкости. Далее эта гипотеза была развита Прандтлем и Карманом. При этом делались различные допущения о характере изменения коэффициента lapt1.wmf. Опытами было установлено, что nт изменяется по поперечному сечению потока, возрастая по мере удаления от стенки до некоторого максимума и затем убывая на оси потока. Часто с достаточной для практических расчетов принимают среднее значение lapt2.wmf в ядре потока, т.е. за пределами пристенного слоя.

Целью данной работы является получить приближенные расчетные выражения для коэффициента средней турбулентной вязкости в барботажном слое.

Коэффициент турбулентной вязкости

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

Коэффициент турбулентной вязкости по определению записывается в виде:

lapt3.wmf (1)

где u′ – пульсационная скорость, м/с; l – масштаб турбулентного движения, м.

В вязком подслое lapt4.wmf [1]

lapt5.wmf (2)

где u′, u′ – продольная и поперечная составляющая пульсационной скорости, м/с; u* – динамическая скорость трения, м/с; u – продольная скорость в пограничном слое, м/с; y –поперечная координата, м; δ1 – толщина вязкого подслоя, м; n – коэффициент кинематической вязкости, м2/с.

а в основной зоне генерации турбулентности lapt6.wmf

lapt7.wmf (3)

где d – толщина гидродинамического пограничного слоя, м.

Во внешней части пограничного слоя (окрестности оси канала)

lapt8.wmf (4)

Для основной зоны генерации турбулентности коэффициент турбулентной вязкости (1) с учетом (3) получит вид

lapt9.wmf (5)

Принимая модель локальной изотропной турбулентности запишем скорость диссипации энергии [2]

lapt10.wmf (6)

Из выражений (5) и (6) получим среднее значение nт при lapt11.wmf

lapt12.wmf (7)

Как отмечено, выражение (6) имеет более широкую область применения [3-9], а не только для условий локальной однородности и изотропности потока.

Турбулентность в жидкой фазе

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

Рассмотрим два режима барботажа в газожидкостных аппаратах [8]:

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

При скорости газа на полное сечение барботажного аппарата (без учета жидкости) lapt14.wmf м/с выражение для расчета динамической скорости получено в виде [8].

lapt15.wmf (8)

а при lapt16.wmf м/с

lapt17.wmf (9)

где ϕг – среднее объемное газосодержание.

Газосодержание вычисляется по известным эмпирическим выражениям [8], например по формуле

lapt18.wmf, (10)

где ρг, ρж – плотности газа и жидкости, кг/м3; s –поверхностное натяжение, Н/м; ∆ρ – разность плотностей жидкости и пузырька, кг/м3.

При выводе выражений (8), (9) было принято, что газ барботирует в условно неподвижную жидкость.

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

lapt19.wmf (11)

τгр – среднее касательное напряжение на поверхности пузырей, Па; lapt20.wmf – скорость движения пузырей, м/с; uгр – скорость движения межфазной поверхности, м/с; Vж – объем жидкости, где происходит диссипация энергии, м3; F – площадь поверхности пузырька, м2.

Учитывая, что толщина гидродинамического слоя на пузырьке на порядок меньше его диаметра запишем объем в виде суммы

lapt21.wmf (12)

где Vсл – объем вихревого следа за пузырьком, м3; lapt22.wmf – объем пограничного слоя на пузырьке, м3; d – средняя толщина пограничного слоя, м.

Представляя объем вихревого следа в виде полуэллипсоида вращения, в центр которого вписана дисперсная частица (пузырек), запишем [10]

lapt23.wmf (13)

где R – радиус пузырька, м; S – длина вихревого следа, м.

При числах lapt24.wmf длина вихревого следа примерно равна [11]

lapt25.wmf (14)

Тогда lapt26.wmf (15)

Тогда из (11) – (15) при F = 4πR2 получим

lapt27.wmf (16)

и выражение (7) запишется в виде

lapt28.wmf (17)

Учитывая, что lapt29.wmf среднюю толщину пограничного слоя найдем из выражения для пластины с учетом движения межфазной поверхности [12]

lapt30.wmf (18)

где L = πR полупериметр длины обтекания пузырька, м.

Известно также выражение lapt31.wmf [13]; lapt32.wmfчисло Рейнольдса для пузырька.

Среднюю скорость на границе раздела можно вычислить, используя функцию Адамара – Рыбчинского

lapt33.wmf (19)

где lapt34.wmf отношение коэффициентов динамической вязкости дисперсной и сплошной фаз; j –угол.

В выражениях (17), (18) диаметр dп =2R и скорость всплытия пузыря вычисляются по известным зависимостям

lapt35.wmf (20)

lapt36.wmf (21)

где do – диаметр отверстия барботера, м.

Из (7), (8) и (17) в итоге получим

lapt37.wmf (22)

Режим стесненного движения пузырей.

Переход к такому режиму можно оценить по условию lapt38.wmf где jж – удельный объем жидкости, находящийся в пограничных слоях всех пузырей и их вихревых следах. Из геометрических выводов имеем lapt39.wmf Как показывают расчеты, для системы воздух – вода, такой режим начинается при lapt40.wmf м/с.

Если lapt41.wmf то имеем lapt42.wmf и при таком режиме энергия диссипирует во всем объеме жидкой фазы и в газовых пузырях. Тогда выражение (11) запишем, используя удельную поверхность газожидкостного слоя

lapt43.wmf (23)

lapt44.wmf (24)

где Vг.ж – объем газожидкостного слоя, м3; Vж – объем жидкости в слое, м3.

Если скорость газа в слое превышает скорость свободного всплытия пузыря, то lapt45.wmf

Учитывая связь

lapt46.wmf (25)

уравнение (24) получит форму

lapt47.wmf (26)

Из (7) и (26) получим

lapt48.wmf (27)

где диаметр пузырьков

lapt49.wmf (28)

где Vo – расход газа через отверстие, м3/с; Sк – площадь поперечного сечения аппарата, м2; n – число отверстий в барботере.

Окончательно с u* (9) имеем

lapt50.wmf (29)

На рисунке дана зависимость nтж от скорости газа в барботажном слое.

lapt51.wmf м2/с

lap1.tiff

Зависимость коэффициента nтж от скорости газа: 1 – расчет по (22); 2 – расчет по (29). Система воздух – вода

Выводы

В статье на основе использования модели локальной изотропной турбулентности получены выражения для приближенного расчета средних значений коэффициента турбулентной вязкости в ядре жидкой фазы барботажного слоя. В результате расчетов получено, что при небольших скоростях газа lapt52.wmf и как следствие малых значениях газосодержания lapt53.wmf коэффициент турбулентной вязкости увеличивается с повышением скорости газа, примерно lapt54.wmf Этот режим характеризуется нестесненным движением пузырей и турбулентность вызвана, главным образом вихревым следом. При lapt55.wmf м/с начинается стесненное движение пузырей и с повышением скорости газа коэффициент турбулентной вязкости значительно снижается. Это можно объяснить уменьшением масштаба турбулентности l в выражении (5), так как расстояние между пузырями уменьшается. При скоростях газа lapt56.wmf м/с в отверстиях барботера образуются газовые струи (факелы), которые являются источниками генерации турбулентности. Для такого режима расчет ведется с учетом газовых струй [9].

Полученные выражения (22) и (29) пригодны для lapt57.wmf м/с и могут использоваться для расчетов высокослойных барботажных реакторов, аэротенков, флотаторов и других аппаратов с газожидкостными средами, когда концевыми (входными) эффектами можно пренебречь.


Библиографическая ссылка

Лаптев А.Г., Лаптева Е.А. МОДЕЛЬ ТУРБУЛЕНТНОСТИ В ЖИДКОЙ ФАЗЕ БАРБОТАЖНОГО СЛОЯ // Международный журнал прикладных и фундаментальных исследований. – 2013. – № 12. – С. 18-22;
URL: https://applied-research.ru/ru/article/view?id=4489 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674