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

МЕХАНИКА

УДК 536.423:532.52

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

Б. А. Снигерев

Снигерев Борис Александрович, доктор технических наук, ведущий научный сотрудник, Институт механики и машиностроения — обособленное структурное подразделение ФИЦ КазНЦ РАН, Россия, 420111, Казань, ул. Лобачевского, д. 2/31, snigerev@imm.knc.ru

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

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

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

DOI: Г^^ш^Л 0.18500/1816-9791-2019-19-2-182-195

ВВЕДЕНИЕ

Двухфазные пузырьковые турбулентные течения широко распространены в различных областях промышленности, таких как химическая, пищевая, фармацевтическая, атомная и теплоэнергетика. Сложность моделирования этих потоков связана с большим количеством явлений различной природы, поскольку такие течения характеризуются сильным взаимным влиянием несущей и дисперсной фазы, сопровождаются процессами теплообмена, фазовых переходов, коалесценции, дробления и т. д. Поэтому моделирование турбулентных пузырьковых потоков представляет немалый практический интерес и ему посвящено большое количество публикаций [1-3]. Тем не менее, недостаточно изученными остаются вопросы о детальной турбулентной структуре для двухфазных потоков, полидисперности распределения пузырьков, закономерностях динамики двухфазного потока. Это связано с тем, что на структуру течения и гидравлическое сопротивление в газожидкостных потоках влияет большое число гидродинамических и геометрических параметров, что особенно важно в инженерных приложениях. В настоящее время в литературе имеется большое количество экспериментальных данных и расчетных рекомендаций по гидродинамическому сопротивлению и теплоотдаче при течении газожидкостных потоков в трубах разного диаметра (см., например, монографию [3]). Экспериментальному изучению турбулентных пузырьковых течений посвящены работы [4-8]. Важная информация содержится в [5], где изложены результаты экспериментального исследования гидродинамического сопротивления при течении пароводяной смеси в обогреваемой трубе. Определено влияние обогрева на сопротивление трения смеси и дана эмпирическая зависимость для расчета потери давления трения в области небольших паросодержа-ний. Также для сравнения приводятся результаты для необогреваемых труб. Показано, что на относительный перепад давления в трубе при течении газожидкостного потока оказывает влияние не только давление, но и массовая скорость. С увеличением массовой скорости происходит уменьшение относительной потери давления на трение, а также уменьшение значения паросодержания, при котором начинается «кризис сопротивления». В работах [5, 6] в результате измерений гидравлического сопротивления и трения в вертикальных каналах была обнаружена область течения с аномальными значениями этих величин. Измеренные величины сопротивления на порядок превышают значения, получающиеся по соответствующим расчетным методикам. Эта область соответствует режиму пузырькового течения при малых приведенных скоростях жидкой фазы. Наряду с экспериментальными исследованиями газожидкостных потоков проводятся численные моделирования двухфазных пузырьковых течений. На основе эйлерова континуального подхода в работах [9-11] предложены эффективные методы, получены аппроксимации уравнения, описывающие

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

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

1. ПОСТАНОВКА ЗАДАЧИ

Для математического описания пространственного движения двухфазной среды применяется эйлеров подход, основанный на модели взаимопроникающих и взаимодействующих континуумов, который сводится к описанию условий разделенного движения фаз и определению величин, описывающих межфазное взаимодействие. Система уравнений Навье - Стокса, осредненных по Рейнольдсу, для описания динамики газожидкостного потока имеет следующий вид [2,11]:

д д

д1 (6т ат) + — (т 6т Утз) = 0, ^ + ад = 1, т = 1,д, г,з = 1, 2, 3, (1) д д дР

— (ат 6т Утз) + (ат 6т Ут, Утз ) = -ат + ат 6т д% + Мт, +

д

дхз

(дУт* дУтз

ат Мей, т --г -

\ дхз дх,

2 д

3 дХ, К ^-т-ёхг3

п

В системе уравнений (1)-(2) Ь — время, 61,6д — плотности сплошной и газовой фаз соответственно, а\,ад — объемные концентрации несущей и дисперсной фазы, д, — сила тяжести, Ут, — компоненты скорости т фазы, Мт, — силы межфазного взаимодействия, Р — давление, ме«,т — эффективная динамическая вязкость т фазы. Для описания турбулентных характеристик газожидкостного потока применяется расширенная двухпараметрическая модель турбулентности, определяемая через эффективную вязкость , состоящую из динамической вязкости несущей среды М\ат,1, турбулентной вязкости Мь,1, вязкости Мв1,1 [12], учитывающая влияние наличия пузырьков в несущей среде на турбулентные характеристики газожидкостного потока в виде

МвЦ,1 = М1ат, I + Мь,1 + МВ1,1

Для вычисления му применяется формула Колмогорова, а для мшд — соотношение, предложенное в работе [12]:

С 61 к2 6

Мг, I = —-1-, Мв1,1 = С^ъ рг аг ds \Уд, — У,| , Мей, д = — Мей, I, С^ъ = 0.6.

£1 61

Для определения кинетической энергии турбулентности к и скорости диссипации жидкой фазы б используется модель турбулентности к — б, модифицированная для двухфазных сред [11,12]:

д

д

дЬ дх, дхз \ Щ дхз) к

Ог = ^ №, I

тп 2

1 дУтг + ( ЗУтз \

дху \ дх,

Значения постоянных в модели турбулентности приняты равными

С^ = 0.09, С 1 = 1.44, С£2 = 1.92, а; = 1.0.

Для полидисперсного потока вклад межфазного взаимодействия в эйлеровом приближении складывается из силы сопротивления ЫВ, подъемной силы Ек присоеди-

-> у у

ненной силы ""ум

Ыи = -Ыдг = -

Эти силы рассчитываются по следующим формулам:

3 С

Ы£ = 4 адв1 у |Уд, - V;,|

М™ = Сум ад й (дУт1 - , Сум = 0.5. (11) Здесь Св — коэффициент сопротивления частицы:

_ (24(1+ 0.15Яв0-687)/Яв, Яв < 1000,

Св = л

[0.44, Яв > 1000,

где Яв = 2д;\Уд, - Уи— число Рейнольдса относительного движения фаз, — средний диаметр пузырьков.

Коэффициент подъемной силы Ск определяется как

!шт[0.2881апЬ

-0.29, Ео^ > 10,

/(Ео^) = 0.00105Ео^ - 0.0159Ео;^ - 0.0204Ео^ + 0.474, Ео^ = - дд)^/а.

Здесь /(Ео^) — функция, зависящая от числа Этвоша Ео^, а — коэффициент поверхностного натяжения.

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

В данной работе принимаются во внимание только процессы коалесценции и дробления, которые являются основными механизмами, приводящими к полидисперсному распределению пузырьков по размерам. Предполагается, что пузырьки при этом сохраняют сферическую форму. Эволюция спектра капель описывается кинетическим уравнением для удельной концентрации пузырьков п(^,£) [11].

Для расчета полидисперсности пузырьков применяется метод эйлеровых фракций, основанный на разбиении всего спектра частиц на отдельные дискретные классы с фиксированными размерами ^(г = 1,Ы/, Ыf — число фракций). Количество частиц в каждом классе, объем которых находится в пределах между V, и ,

определяется через N = . Метод предполагает возможность обмена

пузырьками между классами частиц с разными диаметрами в результате процессов коалесценции и разрушения.

Введем функцию /%, определяемую как доля концентрации пузырьков ¿-го размера от общего числа пузырьков, для которой выполняются соотношения Л% = а%/ад, ад/% = Nу%. Запишем уравнение для эволюции функции /% в виде

д д ....

т

Ва = Ба^,В*а = Па¿V, В? = ВЬ ¿V, = ВЬ ¿V.

«/ У «/ У «/ V «/ У

Слагаемое Ва определяет прирост концентрации пузырьков за счет слияния частиц с суммарным объемом V, а Па — убыль за счет слияния частиц V с остальными. Изменение концентрации за счет процессов дробления описывают следующие два слагаемых: Ва определяет прирост за счет дробления частиц большего объема, Па — убывание в результате дробления частиц с объемом V. Для завершения постановки задачи требуется задание начального распределения частиц на входе.

В данной работе для определения функций Ва, ВЬ, Па, ПЬ используются эмпирические выражения, полученные в [13,14].

Для замыкания системы уравнений (1)—( 14) требуется найти силы межфазного взаимодействия, определяемые соотношениями (9)-(13). Неизвестный параметр выражающий средний объемно-поверхностный диаметр пузырьков в объеме V, определяется через /%:

М /М^ \

^ = Е л /

Движение газожидкостного потока описывается системой уравнений сохранения для несущей и газовой фаз (1)-(2), дополненной модифицированной двух-параметрической моделью турбулентности для двухфазной среды (5)-(7). Силы межфазного взаимодействия (9)-(13), действующие на элементарный объем среды, определяются через характеристики потока, доли объемного газосодержания ад, среднего диаметра пузырьков Зная характеристики потока в текущий момент времени, решаем систему уравнений (14) и находим распределение /% для каждой группы пузырьков, что позволяет определить в каждом элементарном конечном объеме сетки. Распределение /% дает дополнительную информацию о распределении частиц по размерам в объеме потока двухфазной среды. В рамках данного подхода движение полидисперсного двухфазного газожидкостного потока описывается замкнутой системой уравнений (1)-(7), (9)-(15).

Рассчитанные в данной работе потери давления на трение сравниваются с известными расчетными соотношениями, полученными в результате обработки большого числа экспериментальных данных для различных режимов течения двухфазного потока. Для расчета потери давления на трение в двухфазном потоке АР^р различают два типа моделей: гомогенные и модели с разделением движения фаз [2,15]. Потери давления на трение определяются через значения градиента давления двухфазного потока (ф/^я)^р. Для гомогенной модели двухфазная среда рассматривается в виде смеси с переменной плотностью дТР, вязкостью ц,Тр. Градиент давления вычисляется

по аналогичной формуле для однофазной жидкости в виде

G2

2 D Qtp

где GTP = Ql(Jl) + Qg(Jg) — полный расход смеси, а для плотности gTP и вязкости смеси ßTP применяются соотношения

Qtp =

\Qg Qi J

В (17) x — массовое содержание газовой фазы x = Qg(Jg)/GTP. Для определения коэффициента трения fTP применяется соотношение, аналогичное для однофазного потока в виде

= i16/RTP ReTP < 2000, = GtpD .

[0.079 х Re-p ReTp > 2000, »tp

Для сравнения потери давления рассчитываются с помощью эмпирического соотношения, полученного в работе [15], где авторы исходят из известной формулы для перепада давления при течении однородной смеси. Градиент потери давления на трение для двухфазного потока (dp/dz)9P выражается в виде некоторой функции, зависящей от градиентов давления для однофазных потоков dp/dz)9, (dp/dz)9, предполагая, что однородные среды с жидкостью или газом при течении по трубе имеют тот же расход GTP. В работе [15] аппроксимация для определения (dp/dz)9P предложена в виде

F = (dp/dz)l0 + 2[(dp/dz)To - (dp/dz)lß] x.

В соотношении (19) градиенты давления для однофазных потоков dp/dz )T, (dp/dz)T вычисляются по формулам

G2

(dp/dz)T = Gr: f'o, (dp/dz)T = Gr fgo, 2 D Qi 2 D Qg

где flo, fgo — коэффициенты трения для несущей и дисперсной фаз, вычисляемые аналогично (18).

2. ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ

Для численного решения систем уравнений используется метод конечных объемов с использованием библиотеки OpenFOAM, предоставляющей пользователю инструменты для аппроксимации дифференциальных операторов на структурированных и неструктурированных сетках [16]. При помощи метода Гаусса интегралы по контрольному объему сводятся к поверхностным, а значение искомой функции на поверхности ячейки получается путем интерполяции ее значений в центроидах соседних ячеек. Для конвективных слагаемых дифференциальных уравнений применяется процедура QUICK. Для дискретизации производной по времени использовалась неявная схема Эйлера первого порядка точности. Коррекции поля давления осуществляются по конечно-объемной согласованной процедуре SIMPLEC. В ходе реализации

метода конечных объемов задача сводится к решению систем линейных алгебраических уравнений. Решение системы выполняется методом сопряженных градиентов с предобусловливанием. Расчетная область представляет собой канал с входной и выходной частями и твердыми стенками. Во входном сечении 51 задаются следующие условия:

51 : V1 = V?, Узд = V?, V/ = V1 = V? = У2д = 0, р = р?, ад = ао, к; = к?, е; = е?, И; = 6, / = 1, ¿з = 0.0025, /% = 0 (г = 1, 6, г = 3), = 0.001, ^ах = 0.012,

на выходной границе —

52 : /дх3 = д/%/дх3 = др/дх3 = дк;/дх3 = /дх3 = 0.

На твердой стенке 53 принимаются условия прилипания и граничные условия для турбулентных переменных в виде

Уф = V; = V1 = V1 = 0, /% = 0, к; = 0,е; = ц(д22) = 0,

где (^,т) — компоненты единичного вектора нормали к поверхности. Условия на границе симметрии имеют вид 54: V; = V? = 0.

3. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ

Проверка представленной модели была проведена путем сравнения результатов расчетов с экспериментальными данными. Был выбран эксперимент, проведенный в работе [17], в котором изучалось полидисперсное пузырьковое течение в вертикальной трубе. Эксперимент проводился с водовоздушной смесью при атмосферном давлении и постоянной комнатной температуре.

В таблице представлены основные параметры для выбранного эксперимента, где ), (¿д) — приведенные скорости жидкости и газа соответственно.

Характеристика экспериментов [17] Characteristics of the experiments [17]

Jl) м/с (Jg) М/С a d мм D мм

L1 0.986 0.0473 0.051 2.5 50.8

L2 0.986 0.321 0.231 2.5 50.8

На входе подается монодисперсный пузырьковый поток. С учетом приведенных значений сред вычислялись начальные значения скорости текущей среды. Двухфазное течение в вертикальной трубе предполагается осесимметричным, поэтому для численного моделирования выбирается расчетная область, состоящая из кругового сектора с радиусом г0 = П/2 = 2.54 х 10 3 м, длиной Н = 3.06 м и углом раствора 50. Проведены численные расчеты на сетках конечных объемов, состоящих из Ие = 32600, 84200, 168600 узлов. В сечении плоскости ох?х2 число разбиений по координатам вдоль оси трубы и по длине для различных сеток составляет = 20 х 200, Ь2 = 40 х 400, Ь3 = 60 х 600. Численное решение на сетке Ь2 отличается от Ь3 менее чем на 2%, поэтому все приведенные численные результаты выполнены с помощью сетки конечных объемов с разбиением Ь2.

Задаются следующие параметры компонентов двухфазной среды вода-воздух: р; = 998.2 кг/м3, рд = 0.2 кг/м3, = 1.1 х 10-3 Па-с, ц? = 1.7 х 10-5 Па-с. Распределение частиц по размерам представлено в виде шести групп пузырьков (И; = 6) с минимальным диаметром = 10-3 м и максимальным ¿¡тах = 1.2 х 10-2 м.

На рис. 1 показаны картины течения газожидкостного потока при ) = 0.986 м/с, {<1д) = 0.321 м/с в виде изолинии осевой скорости несущей фазы У^, изолинии объемного газосодержания ад, среднего диаметра пузырьков .

б / b

Рис. 1. Картина течения газожидкостного потока (Х1) в виде изолиний осевой скорости Уд, м/с (а), объемного газосодержания ад (б) и среднего

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

Fig. 1. Flow pattern of the gas-liquid flow

Видно, что на входном участке наблюдается треугольная область перестройки поля скоростей с равномерного профиля скорости пузырьков на входе к параболическому по мере продвижения вверх по трубе, а также увеличение скорости пузырьков в центральной части канала и замедление скоростей пузырьков, находящихся ближе к стенкам. Две зоны равномерного распределения скорости разделяет полоса, характеризующаяся наиболее резким изменением как направления, так и величины скорости. Наиболее выраженно проявляется наличие данной зоны в распределении объемного газосодержания во входной части трубы, что иллюстрирует рис. 1, б. Отчетливо видна зона с равномерным распределением объемного газосодержания, что является следствием потока пузырьков одинакового диаметра во входной части трубы. Поскольку на входе в эксперименте [17] подаются пузырьки постоянного размера с диаметром = 2.5 мм, то отчетливо наблюдается сужающаяся по ширине трубы зона с равномерным распределением ад.

Представлены расчетные и экспериментальные данные по распределению характеристик несущего и полидисперсного потока в развитом турбулентном течении в сечении х2/В = 53.6 вдоль радиуса трубы (г — расстояние от оси трубы). На рис. 2, а приведены сравнительные данные осредненного профиля осевой составляющей скорости несущей фазы газожидкостного потока для (^) = 0.98 м/с

и {¿д) = 0.0473,0.321 м/с (режимы , Ь2). Введение небольшой концентрации пузырьков (ад = 0.051) приводит к тому, что профиль скорости жидкой фазы становится практически плоским в ядре потока (рис. 2, а, кривая /). Влияние пузырьков проявляется в том, что вблизи стенки происходит увеличение градиента скорости. Повышение скорости и концентрации пузырьков (ад = 0.231) приводит к уменьшению влияния дисперсной фазы на профиль скорости. Градиенты профиля скорости в пристеночной области уменьшаются, а сам профиль становится более параболическим, соответствующим однофазному потоку.

а / a б / b

Рис. 2. Сравнение расчетных результатов

При небольших скоростях течения смеси наблюдаются пиковая концентрация пузырьков в пристеночной области и равномерное пологое распределение концентрации пузырьков в ядре потока при различных концентрациях газа (рис. 2, б). Аналогичные результаты в эпюре распределения ag вдоль радиуса трубы подтверждаются экспериментальными данными многих исследователей. Так, седлообразные профили газосодержания вдоль сечения трубы обнаружены экспериментально при течении восходящего газожидкостного потока в [8-10]. В этих работах пиковое распределение объемного газосодержания в пристеночной области и равномерное в ядре потока (называемого «скин-эффектом») получено при течении пузырьковой смеси в трубе при Re, = 6.0 х 103 - 6.0 х 104.

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

На рис. 3 в сечении x3/D = 53.6 показано распределение частиц по размерам в виде доли частиц N соответствующего диаметра d от общего числа пузырьков в области потока Ns. Переход от режима L1 (рис. 3, а) к режиму L2 (рис. 3, б) приводит к увеличению Ns = 0.62 х 106 до 1.22 х 106. Увеличение скорости потока ведет к более интенсивному вовлечению пузырьков в турбулентные пульсационные движения, которые сопровождаются процессами коалесценции пузырьков.

0.0 0.005 0.01 d 0.0 0.005 0.01 d

а / a б / b

Рис. 3. Распределение частиц по размерам для режимов течения L1, L2:

Ni/Ns для L1

and L2

На рис. 4, а показано распределение среднего объемно-поверхностного диаметра пузырьков ds вдоль радиуса трубы для двух режимов. При малых скоростях течения жидкости (кривая 1) наблюдается равномерное распределение среднего диаметра пузырьков по сечению канала. В этом случае в единице объема находится небольшое число пузырьков и их влияние мало, также интенсивность процессов коалесцен-ции и дробления при этом незначительна, заметно некоторое равновесное состояние. При этом число пузырьков, образующихся в результате слияния, примерно равно количеству частиц, которые распались на более мелкие в результате столкновения с другими частицами. Увеличение скорости и концентрации газовой фазы приводит к росту ds до 4 мм, что свидетельствует о преобладании интенсивности процессов коалесценции пузырьков над их скоростью распада.

а / a б / b

Рис. 4. Сравнение расчетных (линии 1, 2) и экспериментальных данных L1, L2 (символы) для среднего диаметра пузырьков ds (а); результатов

APTP/APf (б): 1 - расчет, 2 - [15], 3 - [2] Fig. 4. Comparison of calculated (lines 1, 2) and experimental data L1, L2 (symbols) for the mean diameter of bubbles ds (a); the results of APTP/APf (b): 1 - calculation, 2 - [15], 3 - [2]

На рис. 4, б представлены результаты расчета потери давления на трение в двухфазном потоке APjP по отношению к потерям давления APF в однофазном потоке

при идентичных условиях. Здесь расчеты данной работы помечены линией 1, по формуле (16) —линией 3, (19) — линией 2. Результаты расчетов потерь давления по гомогенной модели (рис. 4, б, линия 3) дают заниженные результаты, а расчеты по соотношениям [15] (рис. 4, б, линия 2) близки с результатами данной работы. Увеличение трения для двухфазного потока по сравнению с однофазным подтверждаются в экспериментах [4-7]. Это объясняется возрастанием градиента скорости в жидкости в пристенной зоне (рис. 2, a) за счет более пологой формы в центральной части для двухфазного по сравнению с однофазным. Увеличение расхода газовой фазы приводит к росту гидравлического сопротивления трения, однако, как отмечено в работах [4-7], с увеличением общего расхода смеси относительное давление АРрР/АРР уменьшается. Добавление газовой фазы в турбулентный поток приводит к возрастанию перепада давлений на трение. Это влияние возрастает с увеличением объемного газосодержания. Увеличение скорости двухфазного потока приводит к уменьшению относительного давления на трение за счет более сильного роста трения в однофазном потоке по сравнению с пузырьковым потоком. Результаты расчетов удовлетворительно согласуются с зависимостями из работы [15].

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

Библиографический список

1. Броунштейн Б. И., Щеголев В. В. Гидродинамика, массо-теплообмен в колонных аппаратах. Л. : Химия, 1988. 336 с.

2. Баттерворс Д., Хьюитт С. Теплопередача в двухфазном потоке. М. : Энергия, 1980.

3. Нигматулин Р. И. Основы механики гетерогенных сред. М. : Наука, 1978. 336 с.

4. Бурдуков А. П., Валукина Н. В., Накоряков В. Е. Особенности течения газожидкостной пузырьковой смеси при малых числах Рейнольдса // ПМТФ. 1975. № 4. С. 137-141.

ЗАКЛЮЧЕНИЕ

328 с.

5. Тарасова Н. В., Леонтьев А. И. Гидравлическое сопротивление при течении пароводяной смеси в обогреваемой вертикальной трубе // ТВТ. 1965. Т. 3, № 1. С. 115-123.

6. Ганчев Б. Г., Пересадько В. Г. Процессы гидродинамики и теплообмена в опускных пузырьковых потоках // ИФЖ. 1985. Т. 49, № 2. С. 181-189.

7. Горелик Р. С., Кашинский О. Н, Накоряков В. Е. Исследование опускного пузырькового течения в трубе // ПМТФ. 1987. № 1. С. 69-74.

8. Kashinsky O. N.., Lobanov P. D., Pakhomov M. A. Experimental and numerical study of downward bubbly flow in a pipe // Int. J. Heat Mass Transfer. 2006. Vol. 49. P. 3717-3722. DOI: https://doi.org/10.1016/j-.ijheatmasstransfer.2006.02.004

9. Зайчик Л. И, Мукин Р. В., Мукина Л. С., Стрижов В. Ф. Развитие диффузионно-инерционной модели для расчета пузырьковых турбулентных течений. Изотермическое полидисперсное турбулентное течение в вертикальной трубе // ТВТ. 2012. Т. 50, № 5. С. 665-675. DOI: https://doi.org/10.1134/s0018151x12040220

10. Пахомов М. А., Терехов В. И. Моделирование турбулентной структуры течения в восходящем полидисперсном газожидкостном потоке // Изв. РАН. МЖГ. 2015. № 2. С. 8389.

11. Pfleger D., Gomes S., Wagner G. H., Gilbert N. Hydrodynamics simulations of laboratory scale bubble columns: fundamentals studies on the eulerian-eulerian modeling approach // Chem. Eng. Sci. 1999. Vol. 54, № 4. P. 5091-5095. DOI: https://doi.org/10.1016/s0009-2509(99)00261-4

12. Sato Y., Sadatomi M., Sekoguchi K. Momentum and heat transfer in two-phase bubble flow-I. Theory // Int. J. Multiphase Flow. 1981. Vol. 7. P. 167-177. DOI: https://doi.org/10.1016/0301-9322(81)90003-3

13. Luo H., Svendsen H. Theoretical model for drop and bubble breakup in turbulent dispersions // AIChE. J. 1996. Vol. 46. P. 1225-1231. DOI: https://doi.org/15.1024/aic.712452004

14. Prince M. J., Blanch H. W. Bubble coalescence and break-up in air sparged bubble columns // AIChE. J. 1990. Vol. 36. P. 1485-1489. DOI: https://doi.org/10.1002/aic.690361004

15. Muller-Steinhagen H., Heck K. A simple friction pressure drop correlation for two-phase flow in pipes // Chem. Eng. Prog. 1986. Vol. 26. P. 297-308. DOI: https://doi.org/10.1016/0255-2701(86)80008-3

16. Weller N. G., Tabor G., Jasak H. A tensorial approach to computational continuum mechanics using object oriented tecniques // Computers in Physics. 1998. Vol. 12. P. 620624. DOI: https://doi.org/10.1063/1.168744

17. Hibiki T., Ishii M., Xiao Z. Axial interfacial area transport of vertical bubble flows // Int. J. Heat Mass Transfer. 2001. Vol. 44. P. 1869-1871. DOI: https://doi.org/10.1016/S0017-9310(00)00232-5

Образец для цитирования:

Снигерев Б. А. Влияние пузырьков на структуру течения и трение в восходящем турбулентном газожидкостном потоке // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2019. Т. 19, вып. 2. С. 182-195. 001: https://doi.org/10.18500/1816-9791-2019-19-2-182-195

The Effect of Bubbles on the Structure of Flow and the Friction in Upward Turbulent Gas-Liquid Flow

B. A. Snigerev

Boris A. Snigerev, https://orcid.org/0000-0003-2134-3355, Institute of Mechanics and Engineering, FRC Kazan Scientific Center of RAS, 2/31 Lobachevsky St., Kazan 420111, Russia, snigerev@imm.knc.ru

This paper presents the computational study results of the ascending gas-liquid flow local structure in a vertical pipe. The mathematical model is based on the use of two-fluid Eulerian approach taking into account the inverse influence of bubbles on averaged characteristics and turbulence of the carrying phase. The equations conservation of mass and momentum quantity of motion in the form of Navier-Stokes equations averaged over Reynolds for each phase are written down. For turbulent stresses the relations under the assumption of the Boussinesq hypothesis are written. Turbulent viscosity for the carrier liquid phase is determined using a two-parameter turbulence model modified for two-phase media. In the equations for transferring the kinetic energy of turbulence and its dissipation, additional terms for the kinetic energy, caused by pulsations of bubbles are introduced. As the main forces, the following components are considered: Archimedes force, resistance force, attachment force, Magnus rotational force, near-wall frictional force. To describe the bubble size distribution in a two-phase flow, an equation for preservation of the number of particles, taking into account the processes of coalescence and breakup is written. To solve the number of bubbles conservation equation, an approach based on the method of fractions is applied. The particle size distribution spectrum is divided into a number of fractions with fixed boundaries. It is assumed that bubbles can be exchanged between different fractions as a result of coalescence and breakups. In the framework of this method, the bubble size distribution is approximated by a piecewise equiprobable. Thus, the problem of describing the spectrum of drops by size reduces to the solution of equations for the volume concentrations of individual fractions. The effect of changes in the degree of dispersion of the gas phase volumetric gas flow rate, the velocity of the dispersed phase to the local structure and surface friction in a two-phase flow are numerically investigated. Comparison of the simulation results with experimental data showed that the developed approach allows to adequately describe turbulent gas-liquid flows in a wide range of gas content and initial bubble sizes.

Keywords: multi-phase flow, turbulence, Euler method, pressure drop. Received: 07.07.2018 / Accepted: 11.11.2018 / Published online: 28.05.2019

References

1. Brounchtein B. I., Chegolev B. B. Gidrodinamika, masso-teploobmen v kolonnykh appa-ratax [Gydrogynamics, mass and heat transfer in bubble columns]. Leningrag, Khimia, 1988. 336 p. (in Russian).

2. Butterworth D., Hewitt G. F. Two-phase flow and heat transfer. Oxford Press, 1977. 321 p. (Russ. ed.: Moscow, Energiya, 1980. 328 p.).

3. Nigmatulin R. I. Osnovy mekhaniki geterogennykh sred [Fundamentals of multiphase media]. Moscow, Nauka, 1978. 336 p. (in Russian).

4. Burdukov A. P., Valukina N. V., Nakoryakov V. E. Special characteristics of the flow of a gas-liquid bubble-type mixture with small Reynolds number. J. Appl. Mech. Tech. Phys., 1975, vol. 16, iss. 4, pp. 592-597. DOI: https://doi.org/10.1007/BF00858302

5. Tarasova N. V., Leontiev A. I. Gidravlicheskoe soprotivlenie pri techenii parovodianoi smesi v obogrevaemoi trube [Hydraulic resistance when the flow of steam-water mixture in a heated tube]. Teploviz. Vys. Temp., 1965, vol. 3, no. 1, pp. 115-123 (in Russian).

6. Ganchev B. G., Peresadko V. G. Processes of hydrodynamics and heat exchange in descending bubble flows. Journal of Engineering Physics, 1985, vol. 49, iss. 2, pp. 879-885. DOI: https://doi.org/10.1007/BF00872635

7. Gorelik P. C., Kashinsky O. N., Nakoryakov V. E. Study of downward bubbly flow in a vertical pipe. J. Appl. Mech. Tech. Phys., 1987, vol. 28, iss. 1, pp. 64-67. DOI: https://doi.org/10.1007/BF00918774

8. Kashinsky O. N., Lobanov P. D., Pakhomov M. A. Experimental and numerical study of downward bubbly flow in a pipe. Int. J. Heat Mass Transfer, 2006, vol. 49, pp. 3717-3722. DOI: https://doi.org/10.1016/j-.ijheatmasstransfer.2006.02.004

9. Zaichik L. I., Mukin R. V., Strizhov V. F. Development of a diffusion-inertia model for calculating bubble turbulent flows: isothermal polydispersed flow in a vertical pipe. High Temperature, 2012, vol. 50, iss. 5, pp. 665-675. DOI: https://doi.org/10.1134/s0018151x12040220

10. Paxomov M. A., Terexov V. I. Modeling of the turbulent flow structure of an upward a polydisperse gas-liquid flow. Fluid Dynamics, 2015, vol. 50, iss. 2, pp. 229-239. DOI: https://doi.org/10.1134/S0015462815020076

11. Pfleger D., Gomes S., Wagner G. H., Gilbert N. Hydrodynamics simulations of laboratory scale bubble columns: fundamentals studies on the eulerian-eulerian modeling approach. Chem. Eng. Sci., 1999, vol. 54, no. 4, pp. 5091-5095. DOI: https://doi.org/10.1016/s0009-2509(99)00261-4

12. Sato Y., Sadatomi M., Sekoguchi K. Momentum and heat transfer in two-phase bubble flow-I. Theory. Int. J. Multiphase Flow, 1981, vol. 7, pp. 167-177. DOI: https://doi.org/10.1016/0301-9322(81)90003-3

13. Luo H., Svendsen H. Theoretical model for drop and bubble breakup in turbulent dispersions. AIChE. J., 1996, vol. 46, pp. 1225-1231. DOI: https://doi.org/15.1024/aic.712452004

14. Prince M. J., Blanch H. W. Bubble coalescence and break-up in air sparged bubble Columns. AIChE. J., 1990, vol. 36, pp. 1485-1489. DOI: https://doi.org/10.1002/aic.690361004

15. Muller-Steinhagen H., Heck K. A simple friction pressure drop correlation for two-phase flow in pipes. Chem. Eng. Prog., 1986, vol. 26, pp. 297-308. DOI: https://doi.org/10.1016/0255-2701(86)80008-3

16. Weller N. G., Tabor G., Jasak H. A tensorial approach to computational continuum mechanics using object oriented tecniques. Computers in Physics, 1998, vol. 12, pp. 620-624. DOI: https://doi.org/10.1063/L168744

17. Hibiki T., Ishii M., Xiao Z. Axial interfacial area transport of vertical bubble flows. Int. J. Heat Mass Transfer, 2001, vol. 44, pp. 1869-1871. DOI: https://doi.org/10.1016/S0017-9310(00)00232-5

Cite this article as:

Snigerev B. A. The Effect of Bubbles on the Structure of Flow and the Friction in Upward Turbulent Gas-Liquid Flow. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2019, vol. 19, iss. 2, pp. 182-195 (in Russian). DOI: https://doi.org/10.18500/1816-9791-2019-19-2-182-195