автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.16, диссертация на тему:Математического моделирование процесса массопереноса в барботажных устройствах
Автореферат диссертации по теме "Математического моделирование процесса массопереноса в барботажных устройствах"
В БАРБ0ТА1ННХ УСТРОЙСТВАХ
Специальность 05.13.16 - применение вычислительной техники,
математического моделирования и математических методов- в научных
исследованиях
АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук
Волгоград 1237
Работа выполнена на кафедре "Процессы и аппараты химических м производств" Волгоградского государственного технического университета.
Научный руководитель - доктор технических наук, заслуженный деятель науки и техники РФ, профессор Тябин Н.В.
Научный консультант - кандидат технических наук, доцент
Официальные оппоненты - доктор технических наук, заслуженный
Ведущее предприятие: АО "Химпром" им.С.М.Кирова.
Защита состоится июня 1997г. в 40 час на заседании диссертационного совета К.063.76.05 в Волгоградском государственном техническом университете в аудитории 209 по адресу: 40X56, г.Волгоград, пр.Ленина, 28.
С диссертацией можно ознакомится в библиотеке Волгоградского государственного технического университета.
Автореферат разослан мая 1997г.
Ученый секретарь диссертационного совета
Трусов С. А.
деятель науки и техники РФ, профессор Бойков Г.П.,
кандидат технических наук, доцент Герасименко В.А.
кандитат технических наук, доцент
В.И.Водопьянов
АКТУАЛЬНОСТЬ РАБОТЫ. В развитии теоретических основ процессов массоперетгоса большую роль играпт фундаментальные исследования в области гидродинамики и массообмена. За последние годы в области технологии, где используются среды с различного рода включениями: пузырями, каплями, частицами достигнут значительный прогресс. Тем не менее, до сих пор ощущается недостаток в достоверных физических и математических моделях для описания разнообразных процессов контакта фаз и массопере-носа ме.жду ними., Определяют интенсивность мзссопереноса процессом колонных аппаратов является барботирование в распылительных и барбо-тажных колоннах, в таких технологических процессах, как нагрев жидкой фазы водяным "острым" паром, очистка сточных вод, в том числе от газообразных и твердых примесей, аэрирование культиваторов, экстрагирование и т.п. Дисперсионные среды, которые применяются в этих процессах, обладают как малой, так и высокой еязкостью. Теоретические и экспериментальные исследования процессов барботирования направлены на оптимизацию работы технологического оборудования и выбор оптимальных конструкций барботажных устройств, а разработка достоверных моделей процесса и методик инженерных расчетов является весьма актуальной задачей и представляет теоретический и прикладной интерес.
Продуктивным методом исследования сложных технологических процессов является математическое моделирование, которое при незначительных временных и материальных затратах позволяет получить достоверные результаты, определяет пути совершенствования оборудования и перспективы его развития. Математическая модель должна быть корректно обоснована с . точки зрения принятых допущении и определения области применимости порученных результатов. Для достаточно сложных математических моделей применяются современные методы дифференциально-интегрального исчисления, векторной алгебры и методы численного решения с использованием ЭЕМ.
ЦЕЛЬ РАБОТЫ. Целью настоящей работы является: .
1. Разработка методик инженерных расчетоз контактных разделительных аппаратов: колонн, кубов и других устройств, в которых происходит барботирование газовых или паровых пузырей с целью осуществления процесса массопереноса.
2. Нахождение оптимального коэффициента массоотдачи в широком диапазоне изменения параметров процесса путем аналитического преобразования и численного интегрирования исходных уравнений математической модели, сопоставления , апроксимации опытных данных и полученных результатов математической модели.
• 3. Определение геометрических параметров конструкций барботажных устройств в зависимости от нагрузки по дисперсной фазе.
НАУЧНАЯ НОВИЗНА:
1.Получена достоверная зависимость среднестатистического диаметра пузыря, образующегося в жидкости при. барботировании, от диаметра отверстия барботажного устройства в широком диапазоне изменения параметров процесса барботирования.
2.Из анализа опытных данных и аналитических результатов, предложенных различными авторами, получена зависимость среднестатистической скорости подъема пузыря от его диаметра, вязкости и плотности жидкости при различных режимах движения.
3. Разработана математическая модель для пузырькового режима образования дисперсной фазы, позволяющая описать процессы массопереноса с учетом диаметра отверстия барботажного устройства, расхода газа,' плотности и вязкости сплошной фазы и коэффициента диффузии.
4.Разработан алгоритм и программа численного решения уравнений Навье-Стокса и массопереноса при движении совокупности частиц в высоковязкой жидкости для определения объемного коэффициента массоотдачи.
Полученные результаты теоретических исследований используются для
описания и оптимизации технологических и конструктивных характеристик барботажных устройств с целью интенсификации процесса массопереноса при разделении многофазных систем.
ПРАКТИЧЕСКАЯ ЦЕННОСТЬ. На основе проведенных теоретических исследований разработаны методики инженерных расчетов абсорбционных, ректификационных и экстракционных колонных аппаратов. Полученные теоретические зависимости для определения геометрических и технологических параметров процесса барботирования могут быть использованы при разработке процессов, реализуемых в химической, пищевой, микробиологической и других отраслях промышленности. Результаты диссертационной работы используются в учебном процессе Волгоградского государственного технического университета.
РЕАЛИЗАЦИЯ В ПРОМЫШЛЕННОСТИ: Полученные в диссертационной работе методики расчета колонных аппаратов предложены Волгоградскому акционерному обществу "Химпром" в качестве пособий при конструировании ректификационных и абсорбционных колонных аппаратов.
АПРОБАЦИЯ РАБОТЫ: Основные результаты работы докладывались на международном симпозиуме, проводимом в рамках международного конгресса "Экология, жизнь, здоровье"(Волгоград, 1996г.); 1 межвузовской конференции молодых ученых Волгоградской области (Волгоград, 1994г.); на научно-технических конференциях в Волгоградском.государственном техническом университете (Волгоград, 1994-1997гг.).
ПУБЛИКАЦИИ. По материалам выполненных исследований опубликовано 5 печатных работ.
СТРУКТУРА И ОБЪЕМ РАБОТЫ. Диссертационная работа состоит из
- б -
.введения, пяти глав, общих выводов и приложения. Диссертация содержит 1?5* страниц машинописного текста, 25 рисунков, 2 программы. Список литературы включает 134 наименования. Приложение содержит 2О страниц, 2 программы.
НА ЗАЩИТУ ВЫНОСЯТСЯ СЛЕДУЮЩИЕ РЕЗУЛЬТАТЫ:
- математическая модель процесса образования дисперсной фазы при барботировании газа через жидкость;
- определение зависимости диаметра пузыря при барботировании в широком диапазоне изменения параметров процесса;
- зависимость скорости движения частиц дисперсной фазы в жидкости в широком диапазоне изменения чисел подобия, включающих диаметр частиц, плотность, вязкость жидкости и коэффициент поверхностного натяже-
ф
ния на границе раздела фаз;
- математическая модель, описывающая структуру дисперсной фазы при барботировании;
- зависимости относительных величин газосодержания и поверхности контакта фаз от определяющих чисел подобия;
- методика расчета коэффициента массоотдачи при барботировании газа через жидкость при пузырьковом режиме в маловязких системах и больших расходах газа;
- методика расчета коэффициента массоотдачи при барботировании газа через высоковяз^ую жидкость;
- методика инженерных расчетов интенсивности массопереноса в кон-• тактных аппаратах при разделении .¡родуктов органического и нефтехимического синтеза на основании полученных зависимостей, на примере абсорбционных", ректификационных и экстракционных аппаратов.
СОДЕРЖАНИЕ РАБОТЫ ' '
Во введении обоснована актуальность теш, дана общая
характеристика работы. _. ........... ..... ------
В первой главе приведен обзор работ по экспериментальному и теоретическому исследованию образования и движения твердых, жидких и газообразных частиц в жидкости и массопереносу в дисперсных системах.
Во второй главе предлагаются универсальные зависимости для определения диаметра газовых пузырей (капель) и скорости их всплытия от расхода газа и определяющих чисел подобия, включающих вязкость, плот-еость жидкости и размеры отверстий, и показана корректность выбора чисел подобия, а также физическая и математическая модель пузырькового режима газожидкостного слоя и массопереноса в слое при числах Рейноль-дса больше 10.
В третьей главе представлен алгоритм решения и содержание прог-, раммы для определения поля скоростей при стесненном всплытии пузырей при. ламинарном, режиме -
В четвертой главе приведен алгоритм численного решения для расчета коэффициента массоотдачи в высоковязкой жидкости.
В пятой главе приводится методика инженерного ¡.хчета процессов, протекающих в колонных аппаратах, на примере экстракции.
ФИЗИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПУЗУРЬКОВОГО РЕЖМА ' ГА30ВДК0СТН0Г0 СЛОЯ И МАССОПЕРЕНОС ПРИ ЧИСЛАХ РЕЙНОЛЬЛСА БОЛЕЕ 10
Предлагается-модель, связанная с определением сил, действующих на пузырь от момента его образования до отрыва. Форма пузыря рассматривается как шаровидная и размер пузыря характеризуется диаметром с1.
Скорость газа в отверстии определяется в виде
40
. . С»
где 0 - расход газа в отверстии и а0 - диаметр отверстия.
Движущей силой, определяющей отрыв пузыря, является архимедова сила при известных плотностях среды р и дисперсной фазы р0
А = pg —- . (2)
К консервативным силам относятся:
ltd3
сила тяжести J = p0g —— ; (3)
6
сила поверхностного натяжения R = хдф , (4)
где б - коэффициент поверхностного натяжения. ■
В нашей работе предлагается ует того, что рост пузыря связан с перемещением его поверхности и созданием движения жидкости в прилегающих областях, которое будет оказывать сопротивление на рост пузыря с силой
Fe = — *Pd2*2 • ® '
о
Входящая в уравнение (5) величина й определяется максимальной скоростью перемещения пузыря в момент его роста и является характерной скоростью роста пузыря равной скорости изменения текущего радиуса газового пузыря Rn
®n „ Q
1Г = 2 id2
Коэффициент сопротивления определяется как сумма коэффициентов сопротивлении при ламинарном и турбулентном течениях
in
Число Рейнольдса Reí определяется через характерную скорость роста пузыря и вязкость жидкости д
Reí = — , (8)
и
Отрыв пузыря происходит при достижении равенства архимедовой силы суше трех других сил
А = J + R + Fc . (9)
' Сопоставление обобщенной зависимости, вытекающей из уравнения
0=2-^ = 2-^. (6)
(9), с экспериментальными данными, приведенное в диссертационной работе, позволяет определить 'коэффициент сопротивления, входящий' в выражение динамического воздействия на рост пузыря как
24
X =
Reí
+ 11
(Ю)
Данное выражение коэффициента сопротивления можно использовать для всех режимов образования газового пузыря.
Из анализа выражения сил и уравнении (9)-(10), определены необходимые и достаточные числа подобия, описывающие процесс образования пузыря: число Рейнольдса Re0, число Архимеда Аг0, число Вебера We0:
Мор . d0\>o2P
Т5Г Р^Р
Re0 =
Aro =
WeQ =
(И)
ц - " .6
Использованаэ выражений этих чисел подобия позволяет получить уравнение для определения безразмерного диаметра
с!0
d =
(12)
в следующем виде d5 - 6
Reo2 ¿2
WeoAro
Re
- 20,25 —2 d - 2,15 Ar0
Rec
0
(13)
Ar0
Сопоставление уравнения (13) с уравнениями, описанными в литературе, было проведено для различных чисел Архимеда Аг0 в диапазоне от 1
до 100000 при
Re0' Wec
2
от 0,1 до 10. Указанные диапазоны охватывают прак-
тически все области образования пузырей в дисперсных системах.. На рис. 1 изображены зависимости безразмерного диаметра пузыря от числа Re0, построенные в соответствии с известными уравнениями, указанными на рис.1, и полученного в данной работе уравнения (13) для чисел Аг0
равных 100 и 10000 при
Rec
Wec
1. Максимальное расхождение этих расче-
тов не превышало 30 процентов,, что. является вполне допустимым для практических расчетов. На этом же рис.1 в качестве сопоставления приведена кривая, построенная в соответствии с уравнением Гаддиса Е.С.
г
- и -
(линия 1), которое рекомендуется для расчета диаметра образующихся пу-; аырей во всем диапазоне изменения параметров. Как видно из рис.1, рас-хсщение значении, полученных по уравнению Гаддиса Е.С. со значениями уравнений, ёййеываемых в литературе, достигает 100. и более процентов. Приведенное в Настоящей работе уравнение (13) является трансцейдент-ным, однако дбстаточно точно описывает зависимость диаметра пузырей от определяющих параметров для любых режимов образования.- пузырей .. и используется в данной работе для построения модели газожидкостного слоя в широком диапазоне изменения параметров.
Сопоставление зависимостей описывающих движение пузырей и жидких капель для различных режимов течения} показывает их недостаточную корректность для различных режимов движения жидкотекучих частиц. Наиболее полно уравнения ~ля определения скорости движение пузырей газа приведены в монографии Уоллиса Г., которые в безразмерной форме представлены на рис.2 в виде зависимостей от чисел Рейнольдса Re, Архимеда Аг и модифицированного числа Вебера WeM:
где и - скорость всплытия газового пузыря.
В результате преобразований, варьирования и подбора коэффициентов, в. соответствии с у равнениями приведенными в монографии - Уоллиса Г., получается
+ 0,^Г
Сопоставление уравнения (15) с уравнениями, приведенными в монографии, было проведено для различных значений Меи в диапазоне изменения числа Архимеда Аг от 1 до 100000. На рис.2 изображены зависимости Ре(Аг,Ием), построенные, .в соотвествии с известными и уравнением (15). Графики были построены для Ием= 0,01. Такое значение Уем принимает для системы газ-вода при диаметре пузырей порядка 6*10мм. Указанный
gd3
(14)
размер, как показывает анализ работы барботажных устройств, является оптимальным, так как с увеличением диаметра пузырей существенно уменьшается поверхность контакта фаз и их время пребывания в жидкости, а уменьшение размера пузырей требует изготовления дисперсных устройств с очень малыми отверстиями. Как видно из рис.2, уравнение (15) удовлетворительно согласуется с известными зависимостями в широком диапазоне чисел Рейнольдса, Архимеда и Бебера и может быть рекомендовано для расчета скорости всплытия газовых пузырей. Полученное выражение (15) используется в настоящей работе для построения модели газожидкостного слоя.
Для определения гидравлических параметров пузырькового режима в работе применяется ячеистая модель (рис.3), согласно которой над каждым отверстием в барботажном устройстве движется ряд пузырей газа с постоянным расстоянием друг от друга; диаметр пузырей и расстояние между ними определяются как среднестатистические величины. Тогда весь объем барботажного слоя может быть разбит на прямоугольные ячейки, причем, внутри объема каждой ячейки содержится единственный пузырь шаровидной формы.
Геометрические параметры барботажного слоя при пузырьковом pern« псевдоожижения определяются размерами пузырей и их количеством в объеме слоя. Количество пузырей по высоте Нслоя определяется соотношением количества отверстий п0 и расстояния между ними h
ПоНсЛОЯ с ПрНсдОяО п . _- = 6 -Шэ~, (16)
Поверхность контакта фаз находится через диаметр пузыря
F = nitd2. (17)
Удельная поверхность межфазного контакта между газом и жидкостью определяется отношением суммарной поверхности газовых пузырей к объему барботажного слоя при его поперечном сечении S
а = . (18)
^Полоя
Безразмерная удельная поверхность контакта фаз с учетом чисел Рейнольдса, определяющих процессы образования пузыря Re0 и всплыть пузыря Re
actödo „ JReo
Re '
а =
f,
= 6
(19)
св
где «б - коэффициент эффективного сечения барботера и fCB - Доля свободного сечения барботажного устройства.
Для составления математической модели, в соответствии с уравнением (15), число Рейнольдса Re выражается через Re0, Аг0, We0
Re = d2Re,
(20)
Другим ванжым геометрическим параметром является газосодержание барботажного слоя, то есть отношение суммарного объема пузырей к объему слоя
9 = - ; (21)
ООНслоя
Для описания модели пузырькового слоя следует ввести безразмерную величину газосодержания, которая будет равна
Ф = ф
«б
a d
(22)
fea 6 Re Таким образом, математическая модель для определен поверхности контакта фаз и газосодержания основана на системе нелинейных трансцендентных алгебраических уравнений: (13), (19), (20), (22).
В результате решения уравнений математической модели были определены зависимости безразмерной удельной поверхности контакта фаз и безразмерного газосодержалия от чисел Re0, Аг0 и We0. На рис.4 показаны
зависимости a=f(Re0,We0), которые были построены для Аго=100. При малом поверхностном натяжении б, высоких We0, удельная поверхность
контакта фаз с увеличением Re0 уменьшается за счет увеличения диаметра пузыря. При небольших We0 и Re0 порядка 102 значение безразмерной поверхности контакта фаз падает до минимального 1,41 (линия 1) что обусловлено значительным увеличением диаметра пузыря и его скорости всплывтия - число Re. Аналогичный характер имеют зависимости и для других значений Аг0.
На рис.5 приведены зависимости газосодержания от диаметра отверс- -тия при различных вязкостях жидкости, скорости газа в отверстиях i>0=ím/c и заданном свободном сечении тарелки ( fCB=0.1; fCB=0.05; ?св=0,01). Первое значение выбрано в соответствии с рекомендуемыми значениями свободного сечения для существующих барботажных устройств. Последние два значения выбраны в соответствии с рекомендациями, которые вытекают из'анализа результатов полученных в данной работе. Из графиков видно, что увеличение диаметра отверстия для жидкостей с вязкостью менее 10~5м2/с приводит к увеличению газосодержания за счет увеличения размеров пузырей, однако расстояния между ними по высоте увеличивается, вследствии увеличения их скорости подъема и времени образования пузырей. При этом, как следует из графика (рис.4), удельная поверхность контакта фаз уменьшается. Для жидкостей с вязкостью более 10~5м2/с (кривые 1, 2) наблюдаются минимумы на кривой зависимости относительного газосодержания от диаметра отверстия, что связано с изменением режима движения пузырей (при ламинарном режиме движения U пропорциональна d2, при турбулентном режиме движения U пропорциональна d0'5). Увеличение вязкости приводит к смещению минимума на кривой зависимости газосодержания от диаметра отверстий ?(d0) в область больших значений диаметров. Газосодержание порядка 0,5 и более практически получить не удается и приведено на рис.5 для наглядности и подтверждения, полученного в работе результата о необходимости уменьшения свободного сечения барботажного устройства.
Интенсивность массопереноса определяется скоростью движения
<р
0,8 О,Г 0,6
Ч5 ' 0,4 0,3
аг си
-(
fc¿
Г
t
W-fo4^
О Í 2 3 4 5 6 Рис.5. Зайисимосги прн
— fe.fi «0.1 -—fc6»0,05
•— fei =001
cío.
лоч.
Y Weo-Í
2 Weo-ÍO
3 Veo =00
4 Veo=íOOO
ÍOO /ООО Reo
Рис. б. Зависимости Shyvt>=f("Reo,Wfeo) при Аг? <00 и Scu = ЮО
(fcÄ°o,oi)
Ol
I
пузыря и его разменами. Если пренебречь влиянием сопротивления массо-переносу внутри газового пузыря, то интенсивность массообменного процесса характеризуется коэффициентом массоотдачи ßyv.
Для определения коэффициента массоотдачи в жидкой фазе (с поверхности пузыря в жидкость или наоборот) наиболее достоверным является критериальное уравнение, которое имеет вид
Shyvo = 6 pS ^ (2 + O.öRe^Sc0'33) , (23)
Re ctg
Syvdo2 v где ShyVo = —г--число Шервуда и Scy = —— число Шмидта.
Dy Dy
Таким образом, математическая модель массопереноса при пузырьковом режиме включает уравнения (13), (20), (23). Данная модель не учитывает искажение '^ормы пузыря при всплытии, их коалесценцию, взаимное влияние пузырей друг на друга и достаточно справедливо при газосодер-кании не более 0,3 и размерах пузырей не более 6мм:
В результате расчетов по уравнениям математической модели построены зависимости числа Шервуда ShyVo от чисел Re0, А~ь и We0. В качестве примера на рис.6 показаны графики Shyvo^f(Re0,We0) при Аго=100. Как видно из графика, при больших Re0 процесс массопереноса идет более интенсивно за счет увеличения диаметра пузырей. Вслучае низкого поверхностного натяжения, Weo=1000, . на кривой 4 „зависимости ShyV0=f(Re0) наблюдается минимум, появление которого вызвано изменением режима движения пузырей. Таким образом, при небольших Re0 рост пузыря влечет за собой резкое увеличение его скорости всплытия и уменьшение времени пребывания в жидкости, Shyvo падает•(кривая 4). При Reo>100 рост пузыря вызывает незначительное увеличение его скорости, Shyvo возрастает.
Зависимости объемного коэффициента массотдачи в сплошной фазе от диаметра отверстия и вязкости жидкости приведены на рис.7 при скорости газа в отверстиях й0=1ы/с. На графике показаны зависимости коэффициен-■ та массоотдачи для свободных сечений барботера fCB=0.05 и fCB=0,01. Как видно из рис.7, при увеличении fCB и уменьшении do, происходит
Рис. J . 3émmocTu Jfypfído) m VÖ=1m/C.
интенсификация процесса массопереноса, при этом диаметр пузырей не превышает 5+бмм. Таким образом, при конструировании барботажных устройств следует, уменьшать- диаметры-отверстий и увеличивать свободное сечение барботажного устройства. Причем, необходимо соблюдение условия .1 < с1, которое позволяет определить максимально допустимое значение свободного сечения не более 0,5, а невыполнение его может привести к коалесценции пузырей газа. .
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МАССОПЕРЕНОСА В ПУЗЫРЬКОВОМ СЛОЕ ПРИ ЧИСЛАХ РЕЙНОЛЬДСА МЕНЬШЕ 10 С УЧЕТОМ СТЕСНЕННОГО ДВИЖЕНИЯ ЧАСТИЦ Для ламинарных течений вязкой жидкости без учета образования вихрей, в настоящее время применяются апробированные уравнения Нявье-Стокса; Если течение рассматривать как обратимое, "-то есть пузыри газа неподвижны, а жидкость натекает на пузырь со скоростью равной скорости движения пузыря, то вокруг каждой частицы удается выделить равноценный объем, в центре которого находится пузырь газа. Следует отметить, что рассматриваемый объем жидкости является произвольно выбранным из всех объемов. Такая картина течения корректна при достаточном удалении рассматриваемого объема жидкости с пузырем газа от.. твердой поверхности барботера и верхней границы слоя. Если рассматривать объем, выделенный вокруг частицы, как цилиндрический и не учитывать влияние малых объемов на границах, то течение будет осесимметричным и уравнения Навье-Стокса, с учетом безразмерных переменных:
(24)
примут вид:
+ _1 _0г + 9^0Г ЭР
Зг2 г Эг г2 Зг2 Зг
& + + . (26)
Зг2 г Зг 3z2 32
Для решения данной системы с учетом уравнения неразрывности строится система, которая включает в себя уравнения (25), (26) и третье уравнение, содержащее в качестве неизвестного давление
+ с*)
Зг2 г Зг 3Z2
Для ламинарного режима течение внутри объема будет осесимметрич-ным, поэтому рассматривается только часть объема вокруг пузыря газа (рис.8) и граничные условия вытекают из естественных условий прилипания к поверхности пузыря и симметрии течения относительно оси и поверхности раздела объемов:
1) при г = L ; при г = 0 и г > 1 :
0Г = 0 ; Ц* = 0 ; Щ = 0 ; (22)
Зг Зг
2) при z = h ; при z = 0 и г > 1 :
• иг = 0 ; —2 = 0 ; Р = const ; (29)
. 3i
3) при г2 + z2 = 1 :
Ur = 0 ; Uz = 0 . (30)
Система уравнений (25)-(27), будучи представленной в разностной схеме, решается методом Гаусса-Зейделя, применяемым для системы алгебраических уравнений. Для аналогичных задач метод Гаусса-Зейделя абсолютно сходится, устойчив, и погрешность расчета составляет величин-)' Д2. Причем, не рекомендуется брать малый шаг сетки, так как при зтом резко возрастает время счета. Расчет полей безразмерного давления и безразмерных компонент скорости осуществляется до достижения заданной точности
6Р + 5U < 5. (31)
Исходным уравнением для расчета концентраций внутри движущейся жидкости (рис.9) является уравнение конвективного ■массопереноса, которое с учетом безразмерных переменных:
C-C2 Ur U?. Ud
с =-; Ur - — Ре ; Uz = Ре, Ре = — - число Пекле, (32)
С0-С2 U U Dy
для рассматриваемой задачи имеет вид
иг Щ + и2 Щ = ^ + 11 Щ + . ' (33)
Эг 9г Зг2 г 9г 9z2 Граничные условия вытекают из условий непрерывности концентраций жидкости на поверхности пузыря, материального баланса с учетом потока количества вещества и симметрии концентраций относительно оси и границ рассматриваемой области:
при 2 = +Н с = С2 = 0 ; (34)
- - Cl"C2
при z = -Н с = ci --- ; (35)
С0-С2
- - Эс
при г = 0 . и г = L -- 0 ; (36)
Эг
при г2 + 22 = 1 Со = 1 . (37)
Расчет поля безразмерных концентраций осуществляется методом Га-усса-Зевделя до достижения заданной точности.
С целью нахождения плотности потока диффуг'и рассчитывается безразмерный градиент концентрации к поверхности пузыря через частные производные безразмерных концентраций
9с = CnCi,]] - Ср 9z Дп
где Cnti.Jl - значения безразмерной концентрации в узлах сетки близ лежащих, на шаг Дп, от поверхности пузыря относительно оси г, и
JL = - °° , " " (39)
Эг • Дп
где СпП.ЗЗ -значения безразмерной концентрации в узлах сетки близ ледащих, на шаг Дщ, от поверхности пузыря относительно оси г.
Безразмерный градиент концентрации на поверхности пузыря определяется в виде
<1с , - 8с — = егая с * —
1-г2 + I
<1п Эг 9г
Плотность потока диффузии в безразмерной форме
Ом
Ом^п
dc
Оу(С0-С2) с!п
и безразмерный поток диффузии принимает ш
+1
* = вд^У = ^
(40)
(41)
(42)
Из уравнения "материального баланса находится значение безразмерной концентрации С1 на нижней границе рассматриваемой области, при 2=-Н
Ом
С1-
(43)
лРе / и2г<3г 0
которое сопоставляется с предыдущими значениями и сравнивается с ошибками расчета.
Среднш по объему концентрацию для численного счета удобнее представить в виде двух определенных интегралов:
Сер (г)
-1
№ -Н
-1___
| с{г,г)дг +
-Н
+Н +1
+н___
| С(Г,2)С12 +1
>0.5 +
Г< 1?п
! +Н____
-3- |_с (г, 2) бг
2Н -Н
! к
Сер --I ССр(г)гс1г . (45)
С2 О
Коэффициент массоотдачи характеризуетсся числом Шервуда ИЬуу
вуус!2 20м
= = _ ^ _ . (46)
0у ЯЬ2Н ЛсСр
Таким образом, алгоритм нахождения коэффициента массоотдачи в высоковязкой жидкости построен на последовательном решении системы уравнений (25)-(27), (31), (33), (38)-(46) и имеет целью определение числа Шервуда, а также позволяет построить зависимости числа Шервуда от определяющих чисел подобия 5Ьуу=Ше0,Аго,Иео,5с). Результаты расчета числа Шервуда при {?ео<10 изображены на рис.6 пунктирными линиями. Кая видно из графиков, значения полученные с учетом стесненного движения пузырей несколько вше значений найденных по методике приведенной £ предыдущем параграфе, что обуславливается возникновением циркуляции жидкости между пузырями и интенсификацией процесса массопереноса.
ОБЩЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ
1.На основании динамической модели образования пузырей газа и капель жидкости на бароотажных устройствах получено уравнение, позволяющее определить размер пузырей (капель) в зависимости от расхода газа, диаметра отверстия и реодинамичеалх свойств жидкости.
2.Полученная зависимость, представленная в безразмерном виде, позволила определить основные числа подобия, описывающие процесс образования пузырей, и установить критические значения для определения диапазонов ламинарного, переходного и турбулентного режимов.
3. Обобщены известные экспериментальные и расчетные зависимости по
скорости всплытия газовых пузырей и составлено обобщенное критериальное уравнение для процесса движения пузырей без учета их коалесценции. Уточнены и составлены, критические значения числа Рейнольдса для определения режимов всплытия.
• _ 4. Разработана ячеистая модель для пузырькового режима барботиро-вания, которая позволяет связать его гидродинамические параметры с факторами, определяющими процесс образования, отрыва и движения пузырей в жидкости, и получены исходные уравнения математической модели пузырькового слоя.
5.Составлен алгоритм расчета и программа для определения связи между безразмерными параметрами и числами подобия, которые были получены при обработки математической модели пузырькового слоя.
6.Установлено, что для создания устойчивого пузырькового режима-барботировакия для маяовязких жидкостей (у=10~6*10~7м2/с) следует уменьшать диаметр отверстий и свободное сечение барботажного устройства вплоть до величин 6=0,2мм и Гсв=1%, а для высоковязких жидкостей (V > 10~4м2/с) размеры отверстий и свободное сечение можно увеличивать до с1=6мм и ГСв=6£.
7.Установлены зависимости безразмерного газосодержания и безразмерной удельной поверхности контакта фаз от чисел Рейнольдса, Архимеда, Вебера, которые позволили определить влияние вязкости и плотности жидкости, расхода газа, диаметра отверстий и свободного сечения на условия существования (Кф и" гидродинамические параметры газожидкостного слоя.
8.В соответствии с известными критериальными уравнениями, определяющими массоперенос в системе газ-жидкость гидродинамическая модель пузырькового слоя дополнена моделью массопереноса для пузырькового режима барботирования при переходном и турбулентном режимах образования и движения пузырей в жидкости.
9.На основании приведенной модели получена зависимость числа Шер-
: вуда и объемного коэффициента массопередачи от чисел подобия Рейноль-дса, Архимеда, Вебера и основных параметров процесса (расход газа, ' вязкость и плотность жидкости, диаметр отверстия, свободное сечение барботажного устройства).
10.Расчеты, проведенные на основании математической модели пузырькового слоя, показывают, что для маловязких жидкостей объемный коэффициент массоотдачи при пузырьковом режиме' барботирования имеет практически теже значения, что и дль более интенсивных режимов: струйного и инверсионного барботирования, и пузырьковый режим является вполне приемлимым для этих жидкостей. Однако, при этом отверстие барботажного устройства следует уменьшать до 1+2®, что в настоящее время является трудно выполнимым технологически. Для высоковязких жидкостей (v > 10~4м2/с) пузырьковый режим обеспечивает наиболее высог/ю интенсивность массообменного процесса вследствии развитой поверхности контакта фаз.
11. Разработана математическая модель движения конглемерата пузырей в вязкой жидкости при малых числах Рейнольдса (Re<10которая позволяет определить поля скоростей и давления в жидкости, на основании которых получено численное решение задачи диффузии от поверхности конгламерата пузырей в поток жидкости.
12.Построены зависимости объемного коэффициента массоотдачи от размера отверстий, вязкости жидкости при заданных свободном сечении барботажного устройства и газосодержания барботажного слоя, позволяющая прогнозировать интенсивность массопереноса при ламинарном режиме образования и движения пузырей.
Основное содержание диссертации опубликовано в следующих работах:
1.Жильцова O.A., Трусов С.А., Тябин Н.В. Гидродинамика капельного диспергирования // 1 Межвузовская научно-практическая конференция студентов и молодых ученых Волгоградской области: Тезисы докладов.-.Вол-
гоград, 1994г. - СМЛ5
2.Жильцова O.A., Трусов С.А., Тябин Н.В. Структура газсэдкостной системы при барботировзнии сточных вод // Международный симпозиум, проводимый в рамках международного конгресса "Экология, жизнь, здоровье": Тезисы докладов.- Волгоград, 1395г.- С. 17.
3.Жильцова O.A., Трусов С.А., Тябин Н.В. Определение размеров частиц дисперсной фазы при барботировании сточных вод // Международный симпозиум, проводимый в рамках международного конгресса "Экология, жизнь, здоровье": Тезисы докладов.- Волгоград, 1996г.- С.16.
4.Жильцова O.A., Трусов С.А., Тябин Н.В. Обобщенное критериальное уравнение для всплытия пузырей в вязкой жидкости // Реология, процессы и аппараты химической технологии: Межвуз. сб. науч. тр./ ВолгГТУ.-Волгоград, 1996г.- С.27-30.
•5.Жильцова O.A., Трусов С.А., Тябин Н.В. Определение диаметра газовых пузырей, всплывающих в жидкости. //' Реология, процессы и аппараты химической технологии: Межвуз. сб. науч. тр./ ВолгГТУ.- Волгоград, 1995г.- С.74-79.
Жильцова Ольга Александровна
■—-
-
Похожие работы
- Повышение эффективности атмосферных деаэрационных установок с барботажными устройствами
- Повышение эффективности работы малых биореакторов для анаэробной переработки органических отходов животноводства
- Совершенствование технологии десорбции кислорода в струйно-барботажных деаэраторах атмосферного давления
- Исследование и разработка вихревых аппаратов с вращающимся многофазным слоем
- Получение биогаза в биореакторе с барботажным перемешиванием
-
- Системный анализ, управление и обработка информации (по отраслям)
- Теория систем, теория автоматического регулирования и управления, системный анализ
- Элементы и устройства вычислительной техники и систем управления
- Автоматизация и управление технологическими процессами и производствами (по отраслям)
- Автоматизация технологических процессов и производств (в том числе по отраслям)
- Управление в биологических и медицинских системах (включая применения вычислительной техники)
- Управление в социальных и экономических системах
- Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей
- Системы автоматизации проектирования (по отраслям)
- Телекоммуникационные системы и компьютерные сети
- Системы обработки информации и управления
- Вычислительные машины и системы
- Применение вычислительной техники, математического моделирования и математических методов в научных исследованиях (по отраслям наук)
- Теоретические основы информатики
- Математическое моделирование, численные методы и комплексы программ
- Методы и системы защиты информации, информационная безопасность