Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,087

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

Игнатов В.Е. 1 Марков А.А. 1 Шахов С.В. 1 Нестеров Д.А. 1
1 Воронежский государственный университет инженерных технологий
В статье для характеристики в струйном аппарате двухфазного потока, состоящего из жидкой и газовой фаз, определены показатели и основные уравнения, которым подчиняется работа: законы сохранения массы, импульса и энергии. При этом двухфазный поток в камере смешения струйного аппарата разделен на два участка: начальный и зону молочно-белой эмульсии. Закономерности течения двухфазной среды на начальном участке в камере смешения описаны системой дифференциальных уравнений. Входными данными для расчета парамеров в зоне гомогенного потока стали данные последнего шага расчета начального участка. После чего получена система дифференциальных уравнений, описывающих процессы гидродинамики и массопередачи в области молочно-белой эмульсии. Полученная модель позволяет адекватно описать массообменные процессы в струйных аппаратах.
струйные аппараты
напитки
насыщение
кислород
инжекция
1. Разработка способа приготовления бутилированной кислородонасыщенной воды и комплекса для его осуществления / А.А. Марков, С.В. Шахов, В.Е. Игнатов, Р.В. Шевляков // Современная российская наука глазами молодых исследователей: Сб. статей. (по итогам II Всероссийской научно-практической конференции молодых ученых и специалистов) Т.2. – Красноярск: Изд. Научно-инновационный центр, 2012. – С. 201–205.
2. Разработка установки приготовления кислородонасыщенной питьевой воды / А.А. Марков, Е.Д. Бурхан, В.Е. Игнатов, С.В. Шахов // Успехи современного естествознания. – 2012. – № 6 – С. 142–142.
3. Пат. 2311849 Россия, МПК A 23 L 2/54 (2006.01), A 23 L 2/38 (2006.01). Способ приготовления бутилированной кислородонасыщенной воды и комплекс для его осуществления / Кленчищева С.И. № 2005139146/13; Заявл. 15.12.2005; Опубл.10.12.2007.
4. Рид Р., Шервуд Б. Свойства газов и жидкостей (Определение и корреляция): пер. с англ., ред. В.Б. Когона,. – Л.: Химия, 1971. – 704 с.
5. Циклаури Г.В., Данилин В.С., Селезнев Л.И. Адиабатные двухфазные течения– М.: Атомиздат, 1973. – 446 с.

Струйные аппараты очень широко применяются практически во всех отраслях промышленности: в энергетике, в теплофикационных установках, холодильных установках, в водопроводных установках, для транспортировки жидких и твердых материалов. Струйные аппараты успешно применяются для насыщения воды и безалкогольных напитков различными газами, например, диоксидом углерода или кислородом [1, 2, 3].

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

Для характеристики в струйном аппарате двухфазного потока, состоящего из жидкой и газовой фаз, применяются следующие показатели [5]:

φ – объемная концентрация фазы, численно равная отношению объема, занимаемого фазой, к общему объему двухфазной среды.

Очевидно, что

φр + φг = 1, (1)

где φр и φг – соответственно объемные концентрации жидкой и газовой фаз; Х – массовая концентрация фазы, численно равная отношению массы фазы в единице объема к массе двухфазного потока.

xр + xг = 1. (2)

Плотность двухфазного потока определяется по зависимости

ρ = ρрφр + ρгφг = 1. (3)

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

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

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

1. Статическое давление постоянно по сечению камеры смешения.

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

pic_39.tif

Принципиальная схема струйного аппарата:1 – входной штуцер; 2 – рабочее сопло; 3 – приемная камера; 4 – патрубок отвода газа; 5 – игла шприца; 6 – газовая камера сепаратора; 7 – жидкосная камера сепаратора; 8 – капиляр для отвода газа; 9 – резиновая прокладка; 10 – корпус струйного аппарата; 11 – диффузор; 12 – штуцер отвода насыщенной рабочей жидкости

Запишем эти уравнения в дифференциальной форме.

2. Уравнения неразрывности [5]

для газовой фазы:

Eqn41.wmf (4)

для жидкой фазы:

Eqn42.wmf (5)

3. Уравнение сохранения импульса [5]:

для газовой фазы

Eqn43.wmf (6)

для жидкой фазы:

Eqn44.wmf

Eqn44.wmf (7)

где F – сила трения на поверхности раздела жидкостной струи и газовой фазы.

Скорость газовой фазы на выходе в камеру смешения определяется по уравнению

Eqn45.wmf (8)

Значение величины φг определяется по заданному значению основного геометрического параметра

Eqn46.wmf (9)

Сила трения на поверхности раздела фаз, отнесенная к единице длины,

Eqn47.wmf (10)

где сf – коэффициент сопротивления.

Сила трения, отнесенная к единице объема [5]

Eqn48.wmf (11)

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

Eqn49.wmf (12)

где K – коэффициент массопередачи при спутном движении газовой и жидкой фаз.

Коэффициент сопротивления сf при спутном движении жидкостной и газовой струи определяется при числах Рейнольдса Re < 5∙105 по формуле Блазиуса [5]

Eqn50.wmf (13)

Число Рейнольдса подсчитывается по формуле

Eqn51.wmf (14)

где Z – расстояние от выходного сечения сопла до рассчитываемого. При числах Re больше 5∙105

Eqn52.wmf (15)

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

Eqn53.wmf (16)

где C – постоянная Сатерленда. Для умеренных температур можно принять

C = 184,5.

Значение динамической вязкости при умеренных давлениях и температуре 273,2 К

Eqn54.wmf

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

Eqn55.wmf (17)

Закономерности течения двухфазной среды на начальном участке в камере смешения описываются системой дифференциальных уравнений (4), (5), (6), (7). Параметры потоков рассчитываются в пределах каждого шага Δl. Граничными условиями для расчета параметров потоков являются: Eqn95.wmf ‒ начальная концентрация газа в рабочей жидкости; Pн – давление газа в приемной камере; tр – температура рабочей жидкости; tн – температура газа в приемной камере; u0 – коэффициент объемной подачи; m – значение основного геометрического параметра; Pр –давление подаваемой рабочей жидкости.

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

Для расчета двухфазного потока в области молочно-белой эмульсии примем следующие упрощения:

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

2. Изменение структуры потока не сопровождается заметным растворением углекислого газа в рабочей жидкости.

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

На основании допущения об отсутствии массопередачи в сечении распада струи объемное газосодержание потока остается постоянным. Из уравнения сохранения импульса определим скорость двухфазного потока [5]:

Eqn56.wmf (18)

где ωII–III – скорость двухфазного потока после добавления струи; Eqn57.wmf – скорость соответственно потоков газа и рабочей жидкости в конце начального участка.

Уравнение неразрывности для двухфазного потока запишется в виде [5]:

Eqn58.wmf (19)

Отсюда

Eqn59.wmf (20)

Плотность рабочей жидкости ρр значительно больше плотности углекислого газа и можно записать [5]:

Eqn60.wmf (21)

Отсюда

Eqn61.wmf (22)

Массовая концентрация газовой фазы за сечением распада [5]:

Eqn62.wmf (23)

По уравнениям (21), (22), (23) определяются плотность, объемное и массовое газосодержание в гомогенном потоке за сечением распада II–II.

Закон сохранения импульса напишется в виде [5]:

Eqn63.wmf

Дифференцируя это уравнение по dZ, получим [5]:

Eqn64.wmf (24)

где Π = π⋅dк – параметр канала камеры смешения; mр, mг, – соответственно массовые расходы жидкой и газовой фаз; τω – касательное напряжение между двухфазным потоком и стенкой камеры вследствие сил трения.

Уравнение сохранения массы для гомогенного двухфазного потока в дифференциальной форме имеет вид [5]:

Eqn65.wmf (25)

Дифференцируя уравнение (21) и подставляя полученное выражение в уравнения сохранения массы, получаем [5]:

Eqn66.wmf (26)

Подставляя формулу (26) в (24), получаем уравнение изменения давления

Eqn67.wmf (27)

Запишем соотношения масс газовой фаз [5]:

Eqn68.wmf (28)

Eqn69.wmf (29)

Поделив эти уравнения, получаем [5]:

Eqn70.wmf (30)

Дифференцируя уравнения (30), получаем [5]:

Eqn71.wmf (31)

Количество тепла, выделяющегося в результате растворения газа в воде, крайне незначительно, и процесс можно считать изотермическим, то есть подчиняющимся закону изотермического сжатия

Eqn72.wmf (32)

В дифференциальном виде уравнение изотермического сжатия имеет вид:

Eqn73.wmf (33)

Подставляя это выражение в формулу (31), получаем:

Eqn74.wmf (34)

Скорость изменения массовой концентрации газа Eqn75.wmf может быть определена из уравнения массопередачи

Eqn76.wmf (35)

Eqn77.wmf (36)

Учитывая, что, Eqn78.wmf получаем:

Eqn79.wmf (37)

Выражаем плотность газа из уравнения изотермы:

Eqn80.wmf (38)

Подставляем это выражение в формулы (33) и (31)

Eqn81.wmf (39)

Комбинируя это уравнение с уравнением (27), получаем зависимость давления от газосодержания в данном сечении:

Eqn82.wmf (40)

Рассчитывая, что Eqn83.wmf

Eqn84.wmf (41)

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

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

Eqn86.wmf (42)

Дополнительные уравнения:

Eqn87.wmf

Касательное напряжение трения двухфазного потока о стенки камеры определятся по формуле [5]:

Eqn88.wmf (43)

где

Eqn89.wmf (44)

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

Рецензенты:

Глотова И.А., д.т.н., доцент, заведующая кафедрой технологии переработки животноводческой продукции Воронежского государственного агроуниверситета имени Императора Петра I, г. Воронеж;

Шашкин А.И., д.ф.-м.н., доцент, заведующий кафедрой дифференциальных уравнений, декан факультета прикладной математики, информатики и механики ФГБОУ ВПО «Воронежский государственный университет», г. Воронеж.

Работа поступила в редакцию 21.12.2012.


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

Игнатов В.Е., Марков А.А., Шахов С.В., Нестеров Д.А. РАСЧЕТ ПАРАМЕТРОВ МАССООБМЕННЫХ ПРОЦЕССОВ В КАМЕРЕ СМЕШЕНИЯ СТРУЙНЫХ АППАРАТОВ // Фундаментальные исследования. – 2013. – № 1-2. – С. 392-396;
URL: http://www.fundamental-research.ru/ru/article/view?id=30958 (дата обращения: 31.03.2020).

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

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