Расчет распада произвольного разрыва в двухскоростном потоке с несжимаемыми компонентами
Б.Л. Кащырев <Boris.Kantsyrev@mail.ru> Институт Океанологии им. П.П.Ширшова, Россия, Москва, Нахимовский просп. 36
Аннотация. Представлен анализ влияния вида пульсационных слагаемых в системе осредненных уравнений гидродинамики пузырькового потока с несжимаемыми фазами на волновые свойства модели. Получены соотношения, определяющие решение задачи о распаде произвольного разрыва в газожидкостном пузырьковом потоке с несжимаемыми компонентами.
Ключевые слова: пульсационные слагаемые; характеристика; автомодельное решение; решение Римана; распад произвольного разрыва.
1. Введение
В настоящее время в литературе уже устоялись представления о виде системы уравнений двухскоростного движения гетерогенных сред, о структуре слагаемых, входящих в неё. Однако установление соотношений, которые определяют вид коэффициентов при дифференциальных слагаемых, до сих пор представляет собой еще не решенную и актуальную задачу. Действительно, указанными коэффициентами определяется тип и волновые свойства указанной системы. Недостаток информации о коэффициентах часто приводит к ухудшению расчетной модели, несмотря на более детальный учёт межфазных сил и совершаемой ими работы. Поэтому в данной работе рассмотрено получение соотношений, уточняющих коэффициенты при дифференциальных слагаемых в уравнениях импульса и энергии гетерогенного пузырькового потока. Это дает возможность конкретизировать волновые свойства системы осредненных уравнений, т.е. характеристические скорости и автомодельные решения типа центрированной волны и получить решение задачи Римана, необходимое для построения численной схемы, аналогичной схеме С.Годунова [1] . Рассмотрим режим движения с достаточно большими характерными временами, для которых можно пренебречь сжимаемостью фаз, а скорость изменения размера пузырька по порядку величины не превышает
относительной скорости фаз [2, гл.6, §6]. Прежде, чем записать систему уравнений двухскоростного движения газожидкостного потока учтём, что при выводе осреднённых уравнений многофазной гидродинамики получаются дополнительные слагаемые, соответствующие вкладу пульсационных напряжений, аналогичных рейнольдсовым напряжениям в турбулентных потоках. Действительно, в соответствии с [3], (гл 1, §1), систему уравнений гидродинамики для фазы сплошной среды , обозначенной индексом «¡» можно представить в виде :
0, ад = д_ ^ +у.к о= у^'к + (1)
& дГ
где
4к=0; с;к; а>;-д;к
/' = о ; д, ; д, и
При осреднении каждое из слагаемых в уравнениях заменяется осреднённым. Вводя обозначение для отклонения от среднего значения :
АЕ[=^-Е{, (2)
можно показать, как это было сделано в ([1],гл 1, §2), что среднее значение
(Р = (РГ )(е] )(и1к')+(Р:'ае; ди;к). (3)
Второе слагаемое в правой части ( 3 ) соответствует пульсационным добавкам.
Таким образом, в осреднённой системе уравнений под знаком дифференцирования оказываются не только осреднённые значения переменных, но и пульсационные. Их расчёт в общем случае представляет собой весьма сложную задачу, поэтому в данной работе рассмотрен специальный «ламинарный режим» движения дисперсной смеси. В таком режиме отсутствует хаотичное движение дисперсной фазы, а пульсации скоростей несущей (жидкой) фазы определяются разностью осреднённых скоростей фаз. Как было показано в [ 4, гл 16], этот режим действительно реализуется при не слишком больших скоростях потоков. В данной работе рассмотрено влияние указанных пульсационных слагаемых кинетической энергии и тензора напряжений несущей фазы на волновые свойства системы уравнений пузырькового потока. Кинетическая энергия мелкомасштабных движений и пульсационная составляющая тензора поверхностных напряжений, соответствующая пространственно одномерному осредненному 174
уравнению движения жидкой фазы в указанном ламинарном режиме могут быть представлены [3] соответственно в виде;
к\ = 0.5а2х(а2)°2
Г^ = -0.5а21|/(а2) о2
где а2 - объемное газосодержание, и - относительная скорость фаз и = I \ -I \ V 2 5 V х- соответственно макроскопические (осредненные) скорости дисперсной и несущей фазы, % и у -т.н. пульсационные коэффициенты, которые являются функциями объемного газосодержания При а2 > 0 они рассматриваются, как искомые величины.
2. Система уравнений
Запишем уравнения баланса импульса и энергии для двухфазного потока в целом с учетом сделанных выше предположений.
д^до^У, + р2а2У2) З^ра, + р2а2 + р^У,2 + р2а2У22-рЛ,) ^ Ы дг л
аг ' (6)
а ^«х р^ад ^ + -VII- с,
дг
Индекс 1 относится к жидкой (несущей) фазе, а 2- к дисперсной фазе. В дальнейшем приняты обозначения: с1 К/(Й =с !с\ + I', с /дг, к=1,2. а! + а2 =1,
р = р1а1+р2а2-давление в потоке , р1-р2 = Лд°м2 где Л = — 1- —
4 ^ а
- соответственно, разность между осредненными величинами давления в первой фазе и на поверхности пузырька постоянного радиуса [1, гл.З, стр 130], 1Г=а 1Г1 + а2Г2,
С71 =~р + р°П1 Сг = -р1ау1 - р2а2У2 + У2р° П1
-напряжения поверхностных сил и работа поверхностных сил с учетом пульсационных составляющих. В дальнейшем предположим, что фазы
несжимаемы, фазовый переход не учитывается (газо-жидкостная среда). Уравнения баланса тепловой энергии имеют следующий вид:
(1 € 2
Ае 2 (7)
Р2°Ч = Х2^а1а2 (У2 ~ У + Оп
Где %1 , %2 - коэффициенты, учитывающие доли работы вязких сил, сосредоточенных на межфазной границе, приходящиеся на жидкую и газовую фазы. Уравнения неразрывности для фаз имеют вид:
<3а, <3а, V, —1 + —— = 0,
^ (8) 5а. да^К, —- +—— = 0 дг
В соответствии с рассуждениями, изложенными в [4,§ 139], уравнения двухскоростного движения не должны противоречить законам сохранения импульса и полной энергии, а условие непротиворечивости позволяет конкретизировать вид уравнений движения.. Условием непротиворечивости уравнений движения и законов сохранения в данном случае оказывается соотношение, отражающее тот факт, что притоки тепла от фаз к межфазной границе равны по величине и противоположны:
0и + 0я = о
При этом (5-6) с учётом (7,8) можно преобразовать к следующему виду [6.7];
I ---- = ~V- - Yp.OC ОС , ,
dt dz F dz AFl 7 \ dt dt У
+
0A 2 5a 0 du „ (9)
+ p Ф и —- + p, Ф и— + F, dz dz
d(p°2a2V2) 8(p°2a2v:-a2a1)_ да 2 dt dz dz
2 да2 0 du - p Ф и —- - p Ф и— + F2
~ 1 a ~ 1 и ^ 2
dz dz
dy, d F;
di di
где =(1(а1а2(\|/-Л)) а,
ёа, 2
(Х~ ^а2)+а1а2
¿7. скх„
скх2у_,
Первое слагаемое в правых частях (3.2.5), (3.2.6) представляет межфазную силу Рахматулина. Учет кинетической энергии мелкомасштабных движений в уравнениях (5,6) обеспечивает в (9,10) учет силы присоединенных масс;
хккЧаГ^-^'
причем полученное в [1] значение % =0.5 для а2=0 соответствует известному из литературы предельному значению коэффициента присоединенных масс для уединённой сферы. Слагаемые в правых частях (9) и (10), пропорциональные Фа и Фи соответствуют межфазным силам, обусловленным мелкомасштабными движениями и коллективными взаимодействиями дисперсных частиц с несущим потоком], -Рь-Рг в правой части (9),(10) представляют внешние объемные силы и межфазное взаимодействие, обусловленное вязкостью:
= Р¡а1д2 - Ри, ,Р2 = р°2а2дг + Ри
Р12 =
(П)
.Представляя скорость каждой из компонент в зависимости от полного объёмного расхода И7 и проскальзывания и после исключения из (9), (10) слагаемых, содержащих дР/дг, получим уравнение для проскальзывания фаз:
т у
тг \ди иКи)--
-и2Кс
да2
а ^
дг
(12)
где
К„ =
ОС1-СХ2 — + х(а1-а2)— -р2 р2
фц Р1
а1а2 Р2
J_
P,
l + x
OCj 0U
P,
Pi
P2J Л
.Pi P2J F = -[K)lU + g'(p1-p2)l
P,
PiP:
< рУ
1 + X —
p 2JJ
, Я =9z
dW d t
Рассмотрим уравнение (12) совместно с уравнениями неразрывности для жидкой и газовой фазы.
Sotj i SotjFj q
dt
dz
да0 да^К, ——"- = 0
dt
dz
(13)
(14)
При отсутствии сжимаемости фаз из (13), (14) следует условие независимости полного объемного расхода IV от координаты (что позволяет рассматривать IV как граничное условие - заданную функцию времени) и уравнение для расчёта объемного газосодержания:
da2/dt + д( остГ 2) /9 z = 0.
(15)
Учитывая, что Г 2= W+ щ U, получим из (14)
da2/dt+aia2 dU/d z + [W+{ ai - a2)U]da2/dz = 0.
(16)
Система уравнений (12) и (16) замкнута. В дальнейшем для определённости будем считать, что И'-сош! Поскольку конкретизация любых моделей тесно связана с обобщением экспериментальных данных, в настоящей работе был рассмотрен один из важных режимов течения двухскоростного потока с несжимаемыми фазами- режим непрерывных волн, соответствующий известной модели «потока-дрейфа» [5]. В основу этой модели положены эмпирические соотношения, связывающие относительную скорость фаз и газосодержание, вида г/=г/0(а2). Подставляя эмпирическое соотношение в уравнение неразрывности (14), можно привести его к виду уравнения переноса
да
Я
дг
(17)
где скорость распространения волны
Как известно, в [3,5] на основании ряда сопоставлений результатов расчётов и экспериментальных данных было показано, что указанная модель правильно описывает квазистационарные волновые процессы в двухфазных потоках. Существенно, что эмпирические соотношения для проскальзывания позволяют моделировать не только движение гладких волн, но и скачков газо-и паросодержания, т.н. «непрерывные ударные волны» ( Г. Уоллис). Более того, в [3, Т2 ] с позиций модели дрейфа был успешно проведён анализ такого существенно нестационарного процесса, как распад произвольного разрыва в потоке с несжимаемыми фазами. Таким образом, в рамках указанной модели правильный результат расчёта может быть получен сравнительно просто. С другой стороны, более детальное двухскоростное моделирование, когда каяедой фазе соответствует отдельное дифференциальное уравнение движения, не всегда приводит к правильному результату даже в том случае, когда система уравнений газожидкостного потока гиперболична. Например, первоначально разрывное решение может «размываться» со временем и модель не описывает движение скачка. С учётом этих обстоятельств, представляется целесообразным определить параметры системы уравнений двухскоростного движения исходя из соответствия с волновыми свойствами модели дрейфа в таких режимах течений, когда указанные эмпирические соотношения выполняются. Оказывается, что зависимости %(а2) и \|/(а2) можно единственным образом определить так, чтобы эмпирические соотношения для проскальзывания были частными решениями системы уравнений двухскоростного движения.
3. Волновые свойства системы уравнений гидродинамики пузырькового потока с несжимаемыми фазами
Существенно, что при х(аг) и Ч^гХ определённых таким образом уравнения (12), (16) можно привести к дивергентному виду:
(18)
-i* -i T = const a, --—-
dt dz l^w(«2)J
Данная форма записи компактна и удобна для реализации расчетного метода С.К. Годунова. Следует отметить, что проведённый в [21] анализ также привёл авторов к выводу о совпадении одной из конвективных характеристик системы уравнений многоскоростной гидродинамики со скоростью Кн. и к целесообразности такой конкретизации коэффициентов при дифференциальных слагаемых, моделирующих межфазные силы,
которые обеспечивали бы дивергентный вид уравнений импульса.
4. Расчет распада произвольного разрыва в двухскоростном потоке с несжимаемыми компонентами
Решение задачи о распаде произвольного разрыва (задача Римана) в среде с двумя скоростями создает возможность построения точных численных схем типа схемы Годунова [2] для интегрирования системы уравнений гидродинамики газо- и парожидкостного потока. На практике точное решение нелинейной задачи Римана актуально при моделировании течений с большими градиентами параметров потока, но с другой стороны -весьма сложно. Поэтому в настоящее время разработаны численные методы с линейной аппроксимацией этого решения. Например, в [23-25] представлен основанный на линейной аппроксимации, численный метод P. Roe для двухфазных потоков. Другими словами, решение задачи Римана, используемое на промежуточном этапе вычислений, получается в линеаризованной форме. При этом не всегда гарантируется успешное решение расчетных проблем, обусловленных наличием больших градиентов газосодержания , скоростей, и других параметров течения. Один из перспективных путей решения этих проблем - построение решения задачи Римана, учитывающего движение волн Римана, и разрывов. Например, в работе [22] рассмотрена не гиперболическая система двух уравнений газо-жидкостного потока с несжимаемыми компонентами и получено решение нелинейной задачи Римана с учетом движения контактных разрывов и волн разрежения - сжатия. Основные трудности на этом пути, по-видимому, заключаются в в том, чтобы правильно учесть влияние скоростной неравновесности на решение задачи о распаде произвольного разрыва гидродинамических параметров (давления, скорости, энергии , объемной доли) для каждой из компонент газожидкостного потока. Конкретизация соотношений, определяющих волновые свойства системы уравнений (12),(16), позволяет перейти к построению решения задачи о распаде произвольного разрыва. Поскольку в отсутствие сжимаемости компонент можно отделить задачу нахождения поля скоростей от задачи расчета распределения давления, последнее определяется из решения для 180
скоростей. Полный объемный расход предполагается известной постоянной величиной. Таким образом, также как и в [16] , исходными параметрами считаются проскальзывание и и газосодержание а^. заданные по обе стороны от разрыва. В дальнейшем терминами 'волна сжатия' и 'волна разрежения' обозначаются центрированные волны, в которых изменение плотности среды происходит за счет изменения газосодержания. Кроме того: 1.-в соответствии с классической постановкой рассматривается движение среды в отсутствие объемных и вязких сил. Поэтому система уравнений примет вид:
а(р;а/К + Р;а2у;) | д{р + Р;а^2 + Р;<хгугг - Р;д) = 0 (19)
дЬ дг
дТ ЭТУ
— +-^ = 0 (20)
дг
а(р;а2) ^ д(£а2У2)_ (21)
(3f дг
2,- в соответствии с данным в [1] обоснованием структуры решения предполагается, что первоначальный разрыв распадается на комбинацию вторичных разрывов и центрированных волн, для которых при /=г/[ все зависимые переменные считаются функциями аргумента X. Рассмотрим автомодельные решения системы (19)-(21) типа центрированной волны. Характеристикам ^ и ^ соответствует решение:
ёШа2 = и Фш(а2), (22)
которое получается из (20)-(21) , а также соотношение, определяющее давление, следующее из (19):
с1Р/с1а2 = Р1 и2 Фр1(а2), ( 23)
Фр1 (¡=1,2 - соответствует номеру характеристики) -функция газосодержания, определяемые из (19)- (21). Соотношение, определяющее Фш(а2), для "к= >., =>.(|| имеет вид:
Фи1(а2)=(1/£4)с1ада2 (24)
а для Х= >.9 =12
Фи2( <Х2)=1/<Х1 (25)
По терминологии, принятой в [5 ], данные волновые решения соответствуют динамическим волнам, поскольку они обусловлены силами, пропорциональными градиентам параметров потока (и и а2). Зная зависимости (22, 23), можно определить также изменение давления в гладких волнах, соответствующих характеристикам. На рис 1 они представлены в виде зависимостей давления от плотности среды. Плотность первой фазы 500 кг /м3, второй -10 кг /м3 . По своему смыслу эти зависимости аналогичны адиабатам Пуассона. Для расчета соотношений на скачке, исходя (19)-(21) можно получить также соотношения на разрыве в газожидкостном пузырьковом потоке с несжимаемыми компонентами. В теории одномерных двухфазных течений [5] подобные разрывы называются «непрерывными ударными волнами». Рассматривая всевозможные конфигурации вторичных волн и разрывов, можно показать [7] ,что задача о распаде произвольного разрыва имеет единственное решение.
Рис 1.
На рис 2 представлено распределение газосодержания по координате для различных моментов времени, полученное в соответствии с представленным алгоритмом расчета. Кривая '0' представляет собой распределение газосодержания в начальный момент времени. Слева от разрыва газосодержание а2 равно 0.05 , справа - 0.15. Кривые, соответствующие расчетным вариантам а и Ь соответствуют распаду разрыва с образованием двух гладких волн, причем на кривой а имеется область с 182
нулевым газосодержанием между фронтами волн . Кривой с соответствует движение волны и разрыва. Кривая с1 соответствует решению, построенному в виде комбинации двух разрывов. Следует отметить, что логика представленной схемы распада разрыва соответствует задаче Римана для совершенного газа. Действительно, в классическом решении задачи Римана, как известно, возможно образование двух волн разрежения, волны разрежения и разрыва, а также образование двух разрывов. На основе полученного алгоритма расчета распада произвольного разрыва в потоке с несжимаемыми фазами был создан расчетный код, использующий решение задачи Римана, полученное в [4].
0,25
0,2
0,15
0,1
0,05
Рис 2.
Распределения газосодержания по координате для двух фиксированных моментов времени. В момент времени 1=0 (распределение '0') слева от разрыва 02=0.05 и 11=1 м/с, -справа 02=0.15. Кривые а-(1 соответствуют моменту времени 1=1 с. и различным начальным значениям проскальзывания и справа от разрыва. Соответственно для кривых а- 1.5 м/с , Ь-1.2 м/с, с- 0.8 м/с, (1 - 0.6 м/с.
Результат следующего тестирования показан на рис 3,. где представлены распределения газосодержания для нестационарного противоточного
'2
и
0.5
г[т] 1
течения в вертикальном канале с закрытым дном. Верхнее поперечное сечение канала моделируется, как граница с ёмкостью бесконечно большого объема, заполненной жидкостью Продольная координата увеличивается снизу вверх Видно, что граница растущего уровня жидкости не размывается в течение всего процесса. Кроме того, видно, что стекающий вниз поток постепенно сжимается вследствие ускорения на участке соответствующем продольной координате, большей 0.7 м
1.2
1
С
о 0,8
<5
ч— 0 R
TJ
3
ту 04
—I
0,2
0
Г IV 1)
4 V
Length, m
Рис 3. Г (сек)= 1.-0; 2.-2; 3.-4; 4.-6; 5.-8;
5. Заключение
Получена система уравнений гидродинамики для газожидкостного пузырькового потока, с волновыми свойствами, определяемыми вкладом пульсационных слагаемых. Представлена гидродинамическая модель пузырькового потока, в рамках которой одной из характеристических скоростей газожидкостной дисперсной смеси следует считать скорость переноса возмущений газосодержания, соответствующую модели потока-дрейфа, а другой- скорость пузырьков в потоке.. Осредненная система уравнений имеет правильные волновые свойства только при учете различия давлений компонент. Получены автомодельные решения системы уравнений, соответствующие ее характеристикам. Как следствие данной системы уравнений, записанных в форме законов сохранения, получены соотношения, определяющие ее разрывные решения. Созданная расчётная 184
программа позволяет решать представленную систему уравнений гидродинамики двухскоростного пузырькового потока с помощью метода С.К. Годунова, адаптированного для численного моделирования двухфазной среды.
Список литературы
[1]. Годунов С.К, Численное решение многомерных задач газовой динамики, М.Наука, 1976, 400 стр
[2]. Нигматулин Р.И. Механика сплошной среды, М. Геотар, 2014, 639 стр
[3]. Р.И. Нигматулин. Основы механики гетерогенных сред. М. Наука, 1978, 336 стр
[4]. JI.Д. Ландау, Е.М. Лифшиц, Гидродинамика Т6, М. Наука, 711стр
[5]. Г. Уоллис, Одномерные двухфазные течения, М. Мир, 1972, 440 стр
[6]. В.Г. Kantsyrev. Riemann Problem When Modeling Dual-Speed Bubble Flow.//Proceedings of ICAPP'06, Reno, NV USA, June 4-8, 2006, Paper 6016, pp 1526-1529
[7]. В.Г. Kantsyrev. Modification of Godunov Computation Method for Modeling Non-Stationary Gas-Liquid Flow.// Proceedings of ICAPP'07.Paper 7061, Nice Acropolis, France,May 13-18, 2007
Calculation of the disintegration of any break in the flow of a two-speed two-phase incompressible
B.L. Kantsvrev <Boris.kantsyrev(a)mail.ru>
Institute of Oceanology Russia,Moscow, Nahimovskv broad street 36
Abstract. It is proposed to modify the Godunov computation method for modeling the gas-liquid bubble flow. To have a numerical scheme, effect of the coefficients - at the differential constituents in the hydrodynamic averaged equations set - on the hydrodynamic wave properties is discussed. Specific solution of Riemann problem has been obtined.
Keywords: pulsation terms; characteristic; self-determination; Riemann solution; decay of an arbitrary discontinuity.
References
[1]. Godunov, S. K., A Difference Scheme for Numerical Solution of Discontinuous Solution of Hydrodynamic Equations, Math. Sbomik, 47, 271-306, translated US Joint Publ. Res. Service, JPRS 7225 Nov. 29, 1960
[2]. R.I. Nigmatulin Continuum Mechanics, M. Geotar, 2014, 639 pages, (in Russian).
[3]. R.I. Nigmatulin. Dynamics of multiphase media . Vol .1-2, Hemisphere, New-York, 1991.
[4]. L.D. Landau, E.M. Lifshitz Fluid Mechanics ( Volume 6 of A Course of Theoretical Physics ) Pergamon Press 1959
[5]. G. Wallis. One dimensiona Two-Phase Flow. McGraw-Hill, New York, 1969.
[6]. B.L. Kantsyrev. Riemann Problem When Modeling Dual-Speed Bubble Flow.//Proceedings of ICAPP'06, Reno, NV USA, June 4-8,2006, Paper 6016, pp 1526-1529
[7]. B.L. Kantsyrev. Modification of Godunov Computation Method for Modeling Non-Stationary Gas- Liquid Flow.// Proceedings of ICAPP'07.Paper 7061, Nice Acropolis, France,May 13-18, 2007