Массообмен в абсорбере с трехфазным псевдоожиженным слоем
Расчет локальных скоростей, кинетических характеристик, процессов массопереноса и коэффициентов диффузии трехфазного псевдоожиженного слоя и поверхности контакта фаз. Блок-схема алгоритма инженерного расчета процессов физической абсорбции и хемосорбции.
Рубрика | Химия |
Вид | автореферат |
Язык | русский |
Дата добавления | 14.12.2017 |
Размер файла | 362,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Размещено на http://www.allbest.ru/
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата технических наук
Массообмен в абсорбере с трехфазным псевдоожиженным слоем
05.17.08 - Процессы и аппараты химических технологий
Миронов Евгений Викторович
Иваново 2007
1. Общая характеристика работы
Актуальность проблемы. Задача развития, совершенствования и повышения эффективности технологического оборудования с двух и трехфазными системами и процессов массопереноса в них связана с глубокой проработкой физико-химических, гидродинамических, энергетических, кинетических характеристик и тепло-массообменных процессов в условиях нестационарного движения фаз с учетом конструктивных особенностей аппарата. Решение этой задачи возможно при комплексном рассмотрении нестационарного движения газовой и жидкой фаз и подвижной дисперсной фазы (элементов насадки или катализатора) в аппаратах интенсивного действия. Важен также единый методологический подход к расчету энергетических и кинетических характеристик, движущей силы процесса, тепло-массообменных характеристик и поверхности контакта фаз, позволяющий сократить трудозатраты по оценке качества конструкторских и проектных решений. Разработка методики расчета аппаратов с трехфазным слоем с учетом влияния их конструктивного оформления и реальной гидродинамики фаз является актуальной задачей.
Цель работы. Создание инженерной методики расчета процесса массопереноса в аппарате с трехфазным псевдоожиженным слоем, включающей разработку физических и математических моделей, алгоритмов и программ расчета: кинетических характеристик, движущей силы и поверхности контакта фаз, процессов массопереноса в аппаратах с трехфазными псевдоожиженными слоями с учетом нестационарности движения фаз на единичном элементе насадки и в капле.
Для достижения поставленной цели необходимо решить следующие задачи:
- выполнить анализ исследований по динамике движения газовой и жидкой фаз с энергетическими затратами на реализацию их движений;
- создать физическую и математическую модели движения фаз в секции с трехфазным псевдоожиженным слоем, алгоритм и программу для расчета кинетических характеристик и нестационарной поверхности контакта фаз;
- создать математическую модель и программу для расчета скоростей движения пленки жидкости по поверхностям насадки, стенок аппарата и в отверстиях газораспределительной решетки;
- создать математическую модель и программу для расчета движущей силы процесса массопереноса при нестационарной изотермической абсорбции на микроуровне в газовой ячейке и в секции аппарата с трехфазным псевдоожиженным слоем с учетом конвективного и диффузионного массопереноса.
Научная новизна работы:
- разработаны физическая и математическая модели движения газовой и жидкой фаз в аппарате с трехфазным псевдоожиженным слоем, алгоритм и программа для расчета скоростей движения фаз с использованием полуэмпирической (К-е) модели теории турбулентности;
- установлена взаимосвязь между масштабами энергоемких пульсаций газовой фазы и амплитудами пульсации жидкой и твердой дисперсной фазами, затратами потенциальной и кинетической энергии;
- разработана математическая модель и программа для расчета скоростей движения пленки жидкости по элементу поверхности: насадки, стенки аппарата и в отверстии газораспределительной решетки с учетом турбулентного характера движения газовой и жидкой фаз;
- согласно физической и математической моделей движения фаз в секции аппарата с трехфазным псевдоожиженным слоем и в газовых ячейках усовершенствована двухпараметрическая модель нестационарной абсорбции на микроуровне с учетом конвективного и диффузионного потоков как на границе раздела фаз в движущихся пленках и каплях жидкости, так и в целом по газовой ячейке.
Практическая ценность и реализация работы.
На основе разработанных математических моделей созданы алгоритмы и методика инженерного расчета по проектированию массообменных и пылеулавливающих аппаратов с трехфазным псевдоожиженным слоем.
Разработанная методика расчета кинетических и массообменных характеристик использована при расчете процессов: процесса физической абсорбции и хемосорбции для очистки газовых выбросов в скрубберах производства нитроаммофоски для очистки газовых выбросов от NH3, NOx, HF, SiF4; SO2, и пыли в производстве синтетических моющих средств.
Материалы диссертационной работы используются в учебном процессе Ивановского государственного химико-технологического университета в курсах ПиАХТ и МАХП.
Автор защищает:
1. Энергетический подход к построению физической и математической модели движения газовой и жидкой фаз в аппарате с трехфазным псевдоожиженным слоем.
2. Математическую модель, учитывающую амплитудно-частотные параметры трехфазной системы, алгоритм и программу расчета кинетических характеристик и мгновенной поверхности контакта фаз.
3. Математическую модель, алгоритм и программу расчета движущей силы процесса массопереноса и массового потока в газовой и жидкой фазах, учитывающих амплитудно-частотные характеристики движения фаз, равновесие и химическое взаимодействие фаз.
4. Инженерную методику расчета нестационарного изотермического процесса абсорбции и хемосорбции, учитывающую одновременный конвективный и диффузионный массоперенос на границе раздела фаз в пленках и каплях жидкости.
Апробация работы. Основные положения диссертации докладывались на международных научных конференциях: международная научная конференция «Теоретические основы создания, оптимизации и управления энерго- и ресурсосберегающими процессами и оборудованием» (г. Иваново, 2007); седьмая международная научная конференция «Теоретические и экспериментальные основы создания новых высокоэффективных процессов и оборудования» (г.Иваново, 2005г.); XVIII международная конференция молодых ученых по химии и химической технологии «МКХТ - 2004» (г. Москва, 2004г.); международных студенческих научных конференциях «Фундаментальные науки - специалисту нового века» (г. Иваново, 2002 - 2004 гг.); второй китайско-российско-корейский международный симпозиум по химии и технологии новых материалов (г.Иваново, 2003г.); третья международная научно-практическая конференция «Проблемы строительства, инженерного обеспечения и экологии городов» (г. Пенза, 2001г.); пятая международная конференция «Теоретические и экспериментальные основы создания новых высокоэффективных химико-технологических процессов и оборудования» (г. Иваново, 2001г.); Third international students' conference «environment, development, engineering» (Poland, Cracow, 2001г.); международная конференция и пятый международный симпозиум молодых ученых, аспирантов и студентов (г. Москва, 2001г.); международная студенческая конференция «Развитие, окружающая среда, химическая инженерия» (г.Иваново, 2000г.); международная научная конференция «Жидкофазные системы и нелинейные процессы в химии и химической технологии» (г. Иваново, 1999 г.).
2. Основное содержание работы
В первой главе представлен литературный обзор, посвященный проблемам физического и математического моделирования движения газовой, жидкой и твердой дисперсной фаз, расчета кинетических характеристик, движущей силы и массопереноса в трехфазном псевдоожиженном слое.
Решению проблем физического и математического моделирования аппаратов с трех- и двухфазными системами посвящено большое количество исследований различных школ. Особенно в этом направлении необходимо отметить работы Гельперина Н.И., Айнштейна В.Г., Кваши В.Б., Захарова М.К. (Московская государственная академия тонкой химической технологии им. М.В. Ломоносова); Мухленова И.П., Тарата Э.Я., Островского Г.М., Абиева Р.Ш. (Санкт-Петербургский государственный технологический институт); Балабекова О.С., Сабырханова Д., Шарыгина М.П. (Технологический университет г. Чимкент); Кутепова А.М., Тимонина А.С., Покусаева Б.Г., Соломаха Г.П., Бляхера И.Г., Живайкина Л.Я. (Московский государственный университет инженерной экологии) по исследованию гидродинамики, тепло- и массопереноса в двух- и трехфазных системах. Проблемам математического моделирования аппаратов с двух- и трехфазными системами с учетом иерархической структуры процессов посвящены работы Кафарова В.В., Дорохова И.Н., Мешалкина В.П., Гордеева Л.С., Вальдберга А.Ю. (Российский государственный химико-технологический университет им. Д.И. Менделеева). Решению фундаментальных проблем моделирования двухфазных химических реакторов посвящены работы Слинько М.Г., Вязьмина А.В., Каминского В.А. (Научно-исследовательский физико-химический институт им. Л.Я. Карпова), Малюсова В.А., Шкадова В.Я., Дильмана В.В., Кулова В.Я. (Институт общей и неорганической химии им. Н.С. Курнакова); Холпанова Л.П. (Институт проблем химической физики РАН); Мошинского А.И., Юдаева В.Ф. (ФГУП РНЦ «Прикладная химия», г. Санкт-Петербург), Полянина А.Д. (Институт проблем механики РАН, г. Москва), Накорякова В.Е. (Институт теплофизики им С.С. Кутателадзе СО РАН г. Новосибирск), С.П. Рудобашта (Московский государственный агроинженерный университет им. В.П. Горячкина).
Отмечено, что в настоящее время значительное место в теории и в практике занимают эмпирические и полуэмпирические, а также балансовые модели массообмена. Отсутствуют работы, в которых в рамках единого системного подхода при математическом описании движения фаз были бы использованы энергетические и амплитудно-частотные характеристики, в том числе и с движущейся границей фазового перехода.
Анализ процесса массопереноса в двухфазной системе газ-жидкость в условиях быстрообновляемой поверхности раздела фаз показывает необходимость учета конвективного переноса в жидкой и газовой фазах. Поверхностная конвекция оказывает большое влияние на интенсификацию процессов массопереноса в жидкой и газовой фазах, как при турбулентных, так и при переходных режимах.
Исследования по массопереносу, выполненные для двухфазной системы, не затрагивают сложную гидродинамическую обстановку в аппаратах интенсивного действия, генерируемую турбулентностью с вихрями, составленными из волновых пакетов с широкими амплитудно-частотными спектрами, характерными для промышленных аппаратов.
Выполненный анализ литературных источников показал недостаточно глубокую проработку физических и математических моделей в трехфазном псевдоожиженном слое и двухфазных системах на микро и макроуровнях, наличие большого количества допущений, ставящих под сомнение адекватность физических и математических моделей реальным физическим и химическим процессам, отсутствие в физических и математических моделях взаимосвязанного характера движения газовой, жидкой и твердой фаз (элементов насадки) и влияния на это движение конструктивных особенностей аппарата.
Вторая глава посвящена расчету основных энергетических характеристик системы (турбулентная кинетическая энергия, скорость диссипации кинетической энергии, потенциальная энергия и потенциальная мощность газового потока), которые определяют амплитудно-частотные характеристики, величины коэффициентов турбулентной и эффективной диффузии и поверхность контакта фаз в трехфазной псевдоожиженной системе.
В качестве объекта для определения величин значений кинетических характеристик и поверхности контакта фаз выбран аппарат - абсорбер с трехфазным псевдоожиженным слоем (псевдоожиженные газовым потоком подвижные элементы шаровой насадки с движущейся по ним и стенке аппарата пленке жидкости и подвижные капли жидкости, сорванные с движущихся пленок жидкости или образованные в результате дробления струй).
Абсорбер представляет собой цилиндро-конический аппарат, включающий: циклонную (вихревую) газораспределительную камеру; 2ч3 секции с подвижной шаровой насадкой и с крупнодырчатыми (отверстиями диаметром 6ч20 мм), щелевыми или трубными решетками; сепаратор; расширитель, включающий устройство с тангенциально расположенными лопатками для отделения капель от газового потока.
Формализация конструктивного оформления аппарата с псевдоожиженным слоем насадки проводилась следующим образом, аппарат разбивался на ступени (газораспределения, контактная, сепарационная), каждая ступень, в свою очередь, разделялась на зоны. Газораспределительная ступень представлялась как: ядро, основной тангенциальный поток, пристенная зона. Контактная ступень подразделялась на прирешеточную зону, пристенную зону, зону псевдоожижения и зону орошения.
В каждой зоне выделялись элементы, обтекаемые газовым потоком (капли жидкости первичного и вторичного дробления, элементы шаровой насадки с пленкой жидкости, перемычки между отверстиями с пленкой жидкости в решетке) и каналы (между элементами насадки, каплями, стенкой и элементами насадки, в отверстиях решетки) (рис.1).
Для зон орошения и псевдоожижения определялись диаметры газовых ячеек и каналов между каплями и элементами насадки. Газовая ячейка (рис.1 в) представляет собой единичный элемент насадки с движущейся пленкой жидкости, обтекаемый газожидкостным (капельным) потоком или каплю жидкости, окруженную газовым потоком. Зоны (орошения, псевдоожижения, ядро и основной тангенциальный поток в газораспределителе) дополнительно по высоте зоны разбивались на «псевдотарелки». Исходя из расхода жидкости, удерживающей способности слоя и скоростей газового потока в данных локальных областях по кривой дробления капель в потоке газа с учетом режимов критической деформации и распыления и числа We, определялись диаметры капель первичного дробления, их число и поверхность.
а б в
Рис. 1. Физическая модель трехфазного псевдоожиженного слоя и газовой ячейки: 1 - газораспределительная решетка, 2 - контуры циркуляции, 3 - элементы подвижной насадки, 4 - стенка аппарата, 5 - агрегаты с нисходящим движением, 6 - агрегаты с восходящим движением, 7 и 8 - «псевдотарелки»
Структура трехфазного псевдоожиженного слоя была формализована, как ряд «псевдотарелок». Высота тарелки принималась равной диаметру газовой ячейки в соответствии с порозностью локальной зоны и диаметром насадки (dШ). Количество элементов насадки (NШ) в аппарате, амплитуда пульсаций стохастически образующихся агрегатов из элементов насадки по высоте слоя и порозность псевдоожиженного слоя определяют расстояние между «псевдотарелками» по вертикали (hТ), между элементами насадки по горизонтали(dГ) и количество элементов насадки на одной «псевдотарелке» (lЭ). Число «псевдотарелок» (mТ) зависит от динамической высоты слоя и связано с порозностью. Каждая тарелка в зависимости от режима псевдоожижения и диаметра агрегатов разбита по горизонтали на n агрегатов и k их элементов, совершающих на данной высоте слоя периодически расширение и сжатие с определенными частотами и амплитудами.
Анализ экспериментальных данных исследования динамики слоя стереометрическим методом показал, что с ростом высоты псевдоожиженного слоя частота пульсаций уменьшается, а амплитуда увеличивается. Частота и амплитуда в поперечном сечении слоя постоянны для агрегатов при их восходящем или нисходящем движении на данной высоте. Рост амплитуды пульсации, по-видимому, происходит как за счет увеличения объема газового пузыря, так и за счет роста диаметра каналов, в которых движутся агрегаты, что приводит к уменьшению частоты пульсаций. В трехфазном псевдоожиженном слое больше высоты слоя от решетки 200 мм появляются низкочастотные пульсации (рис. 2). С увеличением плотности орошения с 25м3/м2ч до 90м3/м2ч пульсационный режим движения в псевдоожиженном слое часто переходит в циркуляционный по всей высоте слоя.
Рис. 2. Изменение частоты пульсаций в трехфазном псевдоожиженном слое по динамической высоте слоя
Движение газовой и жидкой фаз представляют собой взаимосвязанные колебательные течения, направленные навстречу друг друга. Механизм пульсационного движения фаз обусловлен нелинейным характером локальных инерционных сил в трехфазной псевдоожиженной системе, а также сил, действующих на деформирующейся границе газ-жидкость.
Размер газовой ячейки определяется в зависимости от объема псевдоожиженного слоя, количества элементов насадки, общей и локальной порозности в слое при размещении их по равностороннему треугольнику(i=1) или квадрату(i=2).
. (1)
Диаметр канала между элементами равен:
. (2)
Количество «псевдотарелок» определяются:
. (3)
В третьей главе выполнен расчет локальных скоростей и коэффициентов диффузии трехфазного псевдоожиженного слоя и поверхности контакта фаз.
Изменение частоты пульсации трехфазного псевдоожиженного слоя по высоте происходит согласно универсальной структуре, называемой обратным (перевернутым) «генеалогическим древом».
Так как механизм турбулентного переноса импульса, массы и тепла один и тот же, то коэффициент турбулентной диффузии DT и коэффициент турбулентной вязкостиТ должны быть близки по величине и определяться по одной методике:
. (4)
Эта формула позволяет наиболее полно учесть энергетические, частотные и конструктивные параметры аппарата с трехфазным псевдоожиженным слоем. Величина коэффициента турбулентной диффузии на границе раздела фаз в газовой ячейке может изменяться по экспоненциальному (или линейному) закону в зависимости от количества волн в волновом пакете. Уравнение для расчета коэффициента эффективной вязкости TijГ в случае экспоненциального закона распределения имеет вид:
, при i=1,2,3 (каналы);(5)
при j=1,2,3 (элемент насадки).(6)
Ввиду того, что перенос импульса и массы в трехфазном псевдоожиженном слое определяется по средней величине из i, j-тых волновых движений, суммарный коэффициент турбулентной диффузии в газовой фазе в зависимости от числа волновых пакетов равен:
. (7)
Коэффициент эффективной диффузии в жидкой фазе принят изменяющимся по экспоненциальному (или линейному) закону от коэффициента турбулентной диффузии на границе раздела фаз до коэффициента молекулярной диффузии на поверхности насадки или аппарата:
.(8)
На границе раздела фаз коэффициенты турбулентной диффузии определяются для жидкости - волновым движением при обтекании пленки жидкости газовым потоком, обтеканием элементов насадки, протеканием газовой фазы по каналу между элементами насадки и обтеканием капель жидкости, дроблением возникающих турбулентных вихрей, бомбардировкой пленки каплями, срывом капель с поверхности пленки, соударением и вращением элементов насадки.
Определение поверхности контакта фаз. Для расчета поверхности контакта фаз принимаем известное положение, согласно которому поверхность пленки покрыта регулярными конусообразными трехмерными волнами с соответствующими амплитудами и частотами.
Кинетическая энергия расходуется на процесс волнообразования и на преодоление сил поверхностного натяжения. Таким образом, поверхность образующейся капли равна:
, (9)
(S=QFN+Skont)
здесь - энергия единичной волны.
Поверхность капель, генерируемая в результате взаимодействия турбулентного газового потока с пленкой жидкости, в слое насадки определится по следующей формуле:
. (10)
Рис.3. Изменение поверхности контакта фаз (S) в трехфазном псевдоожиженном слое (Dа=0,19м, dнас=27мм, Но=0,2м, wг=3,5м/с, П=50м3/м2ч). QFN - поверхность пленок жидкости, Skont - поверхность капель.
Общая волновая поверхность межфазного контакта в аппарате с трехфазным псевдоожиженным слоем, согласно рассчитанным спектрам пленок и капель, может быть выражена в виде зависимости:
. (11)
Расчет толщины пленки жидкости д производится по формулам:
дн = Uг * lнас / нгт * (нІгт/g)1/3 на элементе насадки; (12)
да-р = Uж * lапп / нгт * (нІгт/g)1/3 на стенке аппарата и в отверстиях решетки. (13)
Определение средней скорости течения пленки жидкости по поверхности элемента насадки, в отверстии газораспределительной решетки и на стенке аппарата выполнено для условий динамического равновесия, когда сила тяжести пленки и вышележащего газожидкостного слоя уравновешиваются силой внутреннего трения и трением на границе фаз. На поверхности насадки, решетке и стенке аппарата принято условие прилипания, на границе раздела фаз - равенство касательных напряжений при турбулентном режиме. Скорость движения пленки жидкости на поверхности насадки и по стенке аппарата рассчитывается по формуле:
V=g/(kІнTmax ) · ekд[(kr-1)e-kr +1]+ (гд-фrs)/(kсжнФмах) · ekд(1-e-кr). (14)
Скорость движения пленки жидкости в отверстии решетки определяется по формулам для движения жидкости по стенке аппарата с учетом действия вышерасположенного столба жидкости:
Vжпл = Vж + (2g·Hслссл/сж)Ѕ. (15)
Разработаны алгоритмы расчета скорости стекания пленки жидкости по поверхности насадки, стенке аппарата и в отверстии решетки.
Расчет по трехслойной схеме течения пленки жидкости, выполненной Даклером и Бергелином, дает аналогичные результаты с расчетом по нашим алгоритмам. Сравнение результатов расчета профилей скоростей в турбулентной пленке, свободно стекающей по вертикальной поверхности, проведенных по уравнению данной работы и уравнению Кутателадзе и Стырковича показывает: способ определения скоростей течения турбулентной пленки по данным математическим моделям и основные существующие способы расчета параметров пленки дают сопоставимые результаты.
В четвертой главе для расчета движущей силы на микроуровне качестве элементарного объема аппарата принята газовая ячейка, ограниченная поверхностью нулевого градиента концентрации абсорбируемого компонента. В центре газовой ячейки расположен элемент шаровой насадки с движущейся по нему пленкой жидкой фазы, обтекаемый газожидкостным потоком.
Основой математической модели процессов физической абсорбции и хемосорбции в газовой ячейке аппарата с трехфазным псевдоожиженным слоем принята следующая система балансовых уравнений:
- газовая фаза
; (16)
- жидкая фаза
. (17)
где t - время средне- или высокочастотного пульсационного цикла в потоке, соизмеримого с временем прохождения газовым потоком поверхности контакта фаз (пленки жидкости или капли).
Для замыкания системы дифференциальных уравнений согласно рассмотренной физической модели, принимаем следующие начальные, граничные условия и условия сопряжения:
- по газовой фазе.
;;;; (18)
- по жидкой фазе
;;;; (19)
- на межфазной поверхности
;
|r=(R+д)+0
|r=(R+д)-0 + . (20)
Из условий (на межфазной поверхности) - первое условие дает связь концентрации абсорбируемого компонента в газовой и жидкой фазах через коэффициент распределения б по линейной зависимости (система подчиняется закону Генри), а второе - является условием равенства диффузионных и конвективных потоков на границе раздела фаз газ-жидкость.
Кроме вышеуказанных, при разработке модели приняты следующие допущения: тепловой эффект хемосорбции незначителен, процесс испарения и конденсации пленки жидкости отсутствует.
Для решения задачи используем пространственно временную сетку. Заменяя преобразованную задачу ее конечно-разностным аналогом и используя в качестве итерационной продольно-поперечную экономичную схему Писмана-Рекфорда, которая аппроксимирует дифференциальные уравнения на каждом полушаге, получим систему алгебраических уравнений. По разработанному алгоритму составлена программа расчета. В результате расчета получены профили концентраций в газовой и жидкой (пленке и капле) фазах, значения диффузионного и конвективного потоков, эффективность «псевдотарелок» по высоте.
На рис. 4 и 5 показано изменение безразмерной концентрации в газовой и жидкой фазе для газовой ячейки с каплей и элементом насадки. Время процесса соответствует времени прохождения газовым потоком единичного элемента. C=Cт/С0 - безразмерная концентрация.
газ капля жидкости
Рис. 4. Изменение безразмерной концентрации (С) в газовой и жидкой фазах при физической абсорбции SO2 в газовой ячейке, t=0,0077с. (Dа=0,19м, dнас=27мм, Но=0,2м, wг=3,5м/с, П=50м3/м2ч, б=0,00074)
газа пленка жидкости
газ пленка жидкости
Рис. 5. Изменение безразмерной концентрации в газовой и жидкой фазах при физической абсорбции SO2 в газовой ячейке, а - t=0,0077с., б - t=0,0761с. (Dа=0,19м, dнас=27мм, Но=0,2м, wг=3,5м/с, П=50м3/м2ч, б=0,00074)
Кривые соответствуют шагам по угловой координате и от 0 до р. 20-ая позиция дг в газовой фазе и 1-ая позиция дж в жидкой фазе соответствуют границе раздела фаз, 1-ая позиция дг в газовой фазе - радиусу газовой ячейки, 10-ая позиция дж в жидкой фазе - центру жидкой капли или поверхности элемента насадки.
Для определения гидродинамической структуры жидкостного потока в качестве трассера использовали однонормальный раствор KCl. Изменение концентрации определялось при помощи кондуктометрического датчика. Гидродинамическую структуру слоя и степень продольного и поперечного перемешивания газового потока определяли методом газового трассера. В качестве трассера использовали газ гелий. Кривые отклика обрабатывались по методу моментов с получением значений среднего времени пребывания, дисперсии, скоростей движения фаз, коэффициентов продольной и радиальной диффузии, удерживающей способности по фазам и числа ячеек идеального перемешивания.
Движение элементов меченой насадки фиксировали методом высокоскоростной киносъемки с последующей обработкой кадров киносъемки. При этом определялись траектория элемента насадки, средняя восходящая и нисходящая скорость ее движения, скорость вращения элементов меченой насадки, а также амплитуды и частоты пульсаций поверхности слоя.
Измерение толщины пленки жидкости, на движущимся элементе насадки, осуществлялось методом емкостного зонда. Спектральный анализ осуществляли по осциллограммам мгновенных пульсаций толщины пленки жидкости на стенке аппарата и подвижном элементе шаровой или кольцевой насадки. Обработка данных выполнялась на базе быстрого преобразования Фурье. В результате получали спектры: модуля и плотности мощности, характеризующего энергетический вклад от конструктивных элементов аппарата и агрегатов из элементов насадки.
Величина статической поверхности контакта фаз определялась химическим способом посредством абсорбции CO2 12% раствором NaOH. Величина динамической поверхности контакта фаз определялась стереометрическим методом. Мгновенная удельная межфазная поверхность в псевдоожиженном слое орошаемой насадки растет с увеличением скорости ожижающего агента и плотности орошения. При малых плотностях орошения мгновенная удельная поверхность контакта фаз, измеренная в верхней части слоя, имеет небольшие значения (88 - 260 м2/м3), для мгновенной удельной поверхности контакта фаз над решеткой (в целом по слою) характерны значительно большие величины (880 - 1100 м2/м3).
Определены гидродинамические параметры и коэффициенты математических моделей при Re газа от 2000 до 70000 , Re жидкости от 100 до 3000 и We от 1,7 до 821 для аппаратов объемного контактирования диаметром 0,19 и 1,8 м с псевдоожиженным слоем орошаемой насадки. Наибольший энергетический вклад обеспечивают частоты, определяющие конструктивные параметры аппарата, размеры каналов в агрегатах насадок и каналов между насадками. Анализ амплитудно-частотных характеристик трехфазного псевдоожиженного слоя выявил наличие следующих классов частот: низких (от 0.5 (1.5) до 15 1/с), средних (от 15 до 75 1/с) и высоких (от 75 до 350 1/с и выше при обтекании капель). Соответственно в газовой фазе для низких частот характерны следующие масштабы пульсаций от 9.5*10-3 до 1.5*10-3 м, для средних - от 1.5*10-3 до 3*10-4 м, для высоких - от 3*10-4 до 3*10-5 и меньше. В жидкой фазе амплитуды пульсаций пленки жидкости имеют следующие значения: для низких частот от 3.9*10-3 до 2.1*10-4 м, для средних от 2.1*10-4 до 1*10-5 м, для высоких от 1*10-5 м. и менее.
Пульсационная скорость газа изменялась - от 0,1 до 1,1 м/с; скорость жидкой фазы - от 0,022 до 0,5 м/с; коэффициенты радиальной диффузии газа от 10-4 до 0,8х10-3 м2/с; продольной - от 0,005 до 0,035 м2/с; жидкой фазы - от 0,0005 до 0,027 м2/с; диаметры капель - от 0,00004 до 0,0045 м при их числе, соответственно, от 758320 до 1162 шт. Частота пульсаций поверхности трехфазного псевдоожиженного слоя изменяется от 1 до 4 1/с, амплитуда от 0,4 до 0,03 м; скорости восходящего и нисходящего движения насадки от 0,2 до 0,5 м/с.
Экспериментальные исследования по двухфазным барботажным и двух- и трехфазным псевдоожиженным слоям показывают, что порозность в двух и трехфазных слоях по высоте изменяется. Изменение порозности в соответствии с полученными расчетными и экспериментальными данными принято по линейному закону.
Исследования процесса образования капель при обтекании единичного элемента насадки газовым потоком, совершающим возвратно-поступательные движения, выполняли на модельной установке. В результате обработки эксперимента определяли диаметр капель, который изменялся от 9 мм до 0,2 мм. Средний диаметр капель 1,5-2 мм., причем часть капель представляют из себя дисперсную систему из газо-жидкостных пузырьков с размерами на 1-2 порядка меньше чем диаметры капель.
Рис.6. Сопоставление экспериментальных данных и расчетных значений выходных концентраций SO2 в аппарате с псевдоожиженным слоем орошаемой насадки для Dа=0,19м, dнас=27мм, Но=0,2м, wг=3,0м/с, П=70м3/м2ч, Не=0,00074 и начальных концентрациях 200, 500, 1000, 5000, 10000 мг/м3
Относительная погрешность расчетных и экспериментальных данных не превышает 10 - 15%.
В шестой главе представлена блок-схема алгоритма инженерного расчета процессов физической абсорбции и хемосорбции.
Основные выводы
1. С позиции системного анализа выполнена формализация аппарата с трехфазной псевдоожиженной системой. В соответствии с классификацией уровней иерархии в аппарате выделены: ступени, зоны, газовые ячейки, каналы и элементы, обтекаемые газовом потоком.
2. Разработаны физические и математические модели движения фаз в трехфазной псевдоожиженной системе. В основе физических моделей положено обтекание движущейся пленки жидкости пульсирующим газожидкостным потоком в слое подвижной шаровой насадки. Математическое описание газожидкостного и твердого дисперсного потока учитывает их энергетические характеристики и выполнено на основе формул полуэмпирической теории турбулентности и вынужденных колебаний трехфазной системы.
3. Разработаны методики расчета с учетом конструктивных особенностей аппарата и его внутренних устройств энергетических и кинетических характеристик движения фаз и генерируемой при этом мгновенной (быстро обновляемой) поверхности контакта фаз на основе амплитудно-частотных характеристик волновых пакетов.
4. Разработана математическая модель движения пленки жидкости по элементу насадки, вертикальной стенке аппарата и в отверстии решетки и получены расчетные уравнения для определения толщин и скоростей движения пленок.
5. На базе фундаментального двухпараметрического уравнения конвективной диффузии с учетом амплитудно-частотных характеристик волновых пакетов разработана математическая модель, алгоритм и программа расчета движущей силы процесса массопереноса и массовых потоков газовой и жидкой фаз.
6. На экспериментальной и модельной установках получены необходимые кинетические и массообменные характеристики для проверки адекватности расчетных параметров математических моделей. Подтверждена адекватность математических моделей массопереноса по экспериментальным данным и исследованиям других авторов.
7. Разработана методика инженерного расчета процесса абсорбции и хемосорбции. Методика инженерного расчета использована при проектной проработки вариантов реконструкции промышленных абсорберов для очистки газовых выбросов от SO2, NH3, HF, NOх.
Обозначения: Dа - диаметр аппарата, dн - диаметр насадки, Hd - динамическая высота слоя, С - коэффициент сопротивления, k i,j - кинетическая энергия, еi,j - энергия диссипации, U*,U' - динамическая и пульсационная скорости, у - поверхностное натяжение, fij - частота пульсаций, mвij - масса единичной волны, Aжij - амплитуда пульсаций, nвij - количество волн на поверхности, Nш - количество элементов насадки, л - длина волны, д - толщина пленки жидкости, lнас, lапп - масштабы энергоемких пульсаций для насадки и аппарата, k - коэффициент затухания, ф -касательные напряжения, г - удельный вес, Rш - радиус элемента насадки, Rя - радиус газовой ячейки, W(Ci,T) - скорость химической реакции псевдопервого порядка, б - константа Генри.
массоперенос псевдоожиженный абсорбция
Результаты диссертационной работы изложены в 18 публикациях, основными из которых являются следующие:
1. Миронов, Е.В. Моделирование процесса массопереноса в трехфазном псевдоожиженном слое / Миронов Е.В., Блиничев В.Н., Миронов В.П. // Международная научная конференция «Теоретические основы создания, оптимизации и управления энерго- и ресурсосберегающими процессами и оборудованием». - Иваново, 2007. - т. 1, С. 259 - 265.
2. Миронов, Е.В. Исследование кинетических характеристик для химико-технологических аппаратов с трехфазной системой (ТФС) / Миронов Е.В., Носков С.А., Миронов В.П. // Известия ВУЗов. Химия и химическая технология. - Иваново, 2000. - т. 43, Вып. 6, С. 92-96.
3. Миронов, Е.В. Определение движущей силы и кинетики процессов массопереноса в трехфазной системе на микроуровне / Миронов Е.В., Блиничев В.Н., Миронов В.П. // Успехи в химии и химической технологии. - М., 2004. - Т. XVIII, № 5 (45), С. 107 - 111.
4. Миронов, Е.В. Исследование процесса массопереноса в трехфазной системе с учетом динамических и энергетических характеристик / Миронов Е.В, Блиничев В.Н., Миронов В.П. // Седьмая международная научная конференция «Теоретические и экспериментальные основы создания новых высокоэффективных процессов и оборудования». - Иваново, 2005г. - С. 228 - 230.
5. Mironov, Е. Modeling the Mass- and Heat Transfer in the Apparatus with Two- and Three Phase Systems / Е. Mironov, V. Mironov // Материалы Второго китайско-российско-корейского международного симпозиума по химии и технологии новых материалов. - Иваново, 2003. - С. 33 - 35.
6. Миронов, Е.В. Исследование амплитудно-частотных и кинетических характеристик в аппаратах с трехфазной системой / Миронов Е.В., Миронов В.П. // Межвузовский сборник научных трудов «Процессы в дисперсных средах». - Иваново, 2002. - С. 131 - 135.
7. Миронов, Е.В. Анализ эффективности экологических систем / Миронов В.П., Миронов Е.В. // Третья международная научно-практическая конференция «Проблемы строительства, инженерного обеспечения и экологии городов. - Пенза, 2001г. - С. 24 - 26.
8. Миронов, В.П. Энергетический подход при расчете кинетических характеристик и поверхности контакта фаз в двух и трехфазных системах / Миронов В.П., Кокина Н.Р., Миронов Е.В. // Пятая международная конференция «Теоретические и экспериментальные основы создания новых высокоэффективных химико-технологических процессов и оборудования». - Иваново, 2001. - С. 319 - 324.
9. Mironov, Y. The power approach to calculation of kinetic characteristics and surfaces of contact of phases in two and three-phase systems / Yegeniy Mironov // Third international students' conference «Environment, Development, Engineering». - Poland, Cracow, 2001. - P. 119 - 122.
10. Миронов, Е.В. Теоретическое определение кинетических характеристик для химико-технологических аппаратов с трёхфазной системой (ТФС) / Миронов Е.В., Носков С.А., Прусов С.С. // Международная научная конференция «Жидкофазные системы и нелинейные процессы в химии и химической технологии». - Иваново, 1999 г., - С. 73 - 75.
11. Миронов, Е.В. Экспериментальное определение кинетических характеристик для химико-технологических аппаратов с трёхфазной системой (ТФС) / Миронов Е.В. , Носков С.А., Прусов С.С. // Международная научная конференция «Жидкофазные системы и нелинейные процессы в химии и химической технологии». - Иваново, 1999 г., - С.78 - 79.
Размещено на Allbest.ru
Подобные документы
Физическая сущность абсорбционных процессов. Принципиальная схема циркуляции абсорбента на установках масляной и низкотемпературной абсорбции. Технологические схемы процесса низкотемпературной абсорбции. Основной недостаток низкомолекулярных абсорбентов.
реферат [1,4 M], добавлен 04.04.2017Стадии взаимодействия газообразных реагентов на поверхности твердого катализатора. Соотношение скоростей химической реакции и диффузии на примере необратимой реакции. Расчет адиабатических реакторов для реакций, протекающих в кинетической области.
презентация [428,6 K], добавлен 17.03.2014Диализ - процесс, основанный на различии скоростей диффузии веществ через полупроницаемую мембрану, разделяющую концентрированный и разбавленный растворы. Промышленные аппараты для мембранных процессов. Схема устройства и распределения потоков в аппарате.
курсовая работа [3,8 M], добавлен 02.05.2013Массообменные процессы. Основное уравнение массопередачи. Кинетика диффузионных процессов. Равновесие при абсорбции, дистилляция и ректификация. Простая перегонка. Схема непрерывно действующей ректификационной установки. Экстракция и кристаллизация.
лекция [612,4 K], добавлен 26.02.2014Описание технологической схемы очистки фторсодержащих газов экстракции. Материальный баланс процесса абсорбции в полом абсорбере. Тепловой и механический расчет. Выбор конструкционного материала. Диаметр абсорбера и скорость газа. Расчет вентилятора.
курсовая работа [226,9 K], добавлен 23.04.2015Сущность понятия диффузия. Классификация методов экспериментального исследования. Феноменологическая теория диффузии. Описание безградиентных методов. Основа молекулярно-кинетической теории, ее возможности. Термодинамическая сторона диффузионных явлений.
реферат [22,7 K], добавлен 20.01.2010Физико-химические основы процессов окисления SO2 в системе двойного контактирования и абсорбции. Расчет значения констант равновесия и выхода продукции. Материальный и тепловой балансы процессов. Разработка технологической схемы получения серной кислоты.
дипломная работа [207,8 K], добавлен 23.06.2014Сочетание абсорбции с десорбцией. Поверхностные, барботажные абсорберы. Тарельчатая колона со сливными устройствами. Области применения абсорбционных процессов. Очистка газа от примесей вредных компонентов. Материальный баланс и расход абсорбента.
реферат [165,8 K], добавлен 30.05.2013Зависимость высоты слоя и расхода воздуха от скорости газа. Графическое определение критической скорости газа, определение веса слоя. Теоретическая величина скорости начала взвешивания частиц и сравнение ее со значением, полученным экспериментально.
лабораторная работа [436,3 K], добавлен 18.12.2013Метод капиллярного электрофореза: история появления, основной принцип. Двойной электрический слой. Схема процессов, происходящих на поверхности кварца. Формирование двойного электрического слоя и ход потенциала на границе раздела кварц-электролит.
реферат [217,2 K], добавлен 08.01.2012