Обоснование кинетического метода расчета многокомпонентной изотермической абсорбции

Тип работы:
Реферат
Предмет:
Физика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 66. 081. 2
ОБОСНОВАНИЕ КИНЕТИЧЕСКОГО МЕТОДА РАСЧЕТА МНОГОКОМПОНЕНТНОЙ ИЗОТЕРМИЧЕСКОЙ АБСОРБЦИИ
Неумоина Н. Г., Белов А. В.
Камышинский технологический институт, филиал ФГБОУ ВПО «Волгоградский государственный технический университет», Камышин, Волгоградская область, e-mail: fpt@kti. ru
Предложен кинетический метод расчета многокомпонентной изотермической абсорбции. Основой для разработки метода послужила нелокальная версия термодинамики, предложенная Майковым В. П. В рамках разрабатываемого подхода установлен принцип дискретизации термодинамических параметров, который позволяет ввести в термодинамику время как существенный параметр и определить минимальный макроскопический масштаб для описания процессов переноса массы, тепла, импульса. Для системы «газ-жидкость» получено уравнение удельного мольного потока компонента через границу раздела фаз (уравнение массопередачи). При выводе этого уравнения применено приближение теории пограничного слоя.
В пределах пограничного слоя перепад концентраций происходит по линейному закону. В коэффициентах массоотдачи, отвечающих за транспорт компонента в газовой и жидкой фазе, удалось выделить гидродинамическую и физическую составляющие. Гидродинамическая составляющая зависит от конкретной гидродинамической обстановки в зоне раздела фаз и фактически определяет толщину пограничного слоя. Физическая составляющая показывает, от каких теплофизических свойств контактирующих фаз газа и жидкости зависит коэффициент массоотдачи.
Ключевые слова: абсорбция, массопередача, массоотдача, коэффициент переноса
GROUND OF KINETIC METHOD OF CALCULATION TO MULTICOMPONENT ISOTHERMAL ABSORPTION
Neumoina N.G., Belov A.V.
A Kamyshin technological institute of the «Volgograd state technical university», Kamyshin,
Volgograd area, e-mail: fpt@kti. ru
_The kinetic method of calculation of multicomponent isothermal absorption is proposed. The method is based on the nonlocal version of thermodynamics developed by Maikov V.P. Within a framework of the developed approach, a principle of discretization of thermodynamics parameters is set. This principal allows to introduce time to thermodynamics as an essential parameter and to define a minimum macroscopic scale for description of transfer processes of mass, heat and impulse. For the «gas-liquid» system, an equation of specific molar stream of the component through the boundary of phases' separation (equation of the masstransfer) was derived. For the conclusion of this equation approaching of theory of frontier layer is applied. Within the limits of frontier layer the overfall of concentrations takes place on a linear law. In the coefficients of massreturn responsible for the transport of the component in gas and liquid phase, it was succeeded to hydrodynamic and physical constituents. Hydrodynamic constituent depends on a certain hydrodynamic situation on the border of division of phases and actually determines the thickness of frontier layer. Physical constituent shows: what thermophysical properties of contacting phases of gas and liquid the coefficient of massreturn depends on.
Keywords: absorption, masstransfer, massreturn, coefficient of transfer
В промышленности процесс абсорбции применяется для разделения и осушки углеводородных газов. Увеличение масштабов производства при возросшем в настоящее время внимании к вопросам экологии и качеству готовой продукции требует дальнейшей разработки научно обоснованных и надежных методов расчета абсорбционной аппаратуры. Наиболее полное описание процессов, происходящих при массопередаче в системах «газ-жидкость», возможно на основе кинетических методов расчета. Но применение этих методов затруднено в связи с необходимостью определения коэффициентов массоотдачи и массопередачи [1−3]. Предлагаемый метод расчета кинетики изотермической абсорбции для многокомпонентных систем основан на уравнении массопередачи, которое получено в рамках нелокальной версии термодинамики (НВТ) [4]. Исполь-
зование этого подхода позволяет частично устранить основную трудность, присущую кинетическим методам расчета.
Описание кинетики массообменного процесса на основе НВТ стало возможным в связи с тем, что установленный В. П. Майковым принцип дискретизации термодинамических параметров [4] позволяет внести в термодинамику время как существенный параметр. Рассмотрим кратко основные закономерности дискретизации термодинамических параметров. Приводимые ниже теоретические выводы не опираются на какие-либо модельно-механизменные представления и сохраняют в значительной степени термодинамический характер.
Принцип квантования термодинамических параметров явился следствием того, что естественная граница точности измерения количества теплоты (среднее значение теплового шума кТ где к — константа
¦ FUNDAMENTAL RESEARCH № 10, 2013 ¦
Больцмана, к = 1,38−10~23 Дж/К- Т — термодинамическая температура, К) была принята в качестве минимального значения количества теплоты. Записав второй закон термодинамики в конечных приращениях
AS-AQ/T,
(1)
где ДQ — количество теплоты, Дж- ДS -приращение энтропии, Дж/К- и приравняв ДQ = кТ получим минимальное приращение энтропии
AS = k.
(2)
Здесь и далее минимальные приращения параметров будут обозначаться значком Д.
Используя величину ДQ как минимальную энергию в соотношении неопределенностей энергия-время квантовой физики, можно получить характерный масштаб времени
Ат = П/(2кТ), (3)
где И — постоянная Планка, И = 1,05−10−34 Джх.
Величина Дт определяет размеры пространства, в котором устанавливается локальное термодинамическое равновесие и формируется температура как термодинамический параметр (для нормальных условий Дт имеет порядок 10−14 е). Объем этого пространства для однородной изотропной среды составит
4 з п (ПсЛ3
Гм = ~пг3=- -, (4)
3 6
АГ=6 Гг4.
п сРрН3с3' Am — AVpM,
(6)
(7)
где kS — адиабатический модуль сжатия, н/м2- ср — мольная изобарная теплоемкость. Дж/(кмоль К) — р — молярная плотность, кмоль/м3- M — молярная масса, кг/кмоль.
Выбор условия P = const вызван тем, что в состоянии динамического (флуктуа-ционного) равновесия каждая отдельная макроячейка выступает лишь как область пространства, охваченная электромагнитным взаимодействием за время Ат, т. е. макроячейка не является объемом, который был бы физически фиксирован в определенных границах.
Изменение массы макроячейки — элементарная масса Ат содержит меньше одной частицы. Это значит, что элементарную массу можно рассматривать только как квазичастицу, которая в данном случае является акустическим фононом. Этот акустический фонон и отвечает за перенос массы, тепла и импульса в данной среде.
Так как минимальная скорость распространения акустического фонона (квазичастицы) равна скорости распространения звука в данной среде (это можно доказать), то можно ввести еще один характерный радиус взаимодействия — это радиус упругих взаимодействий между макроячейками
кТ) '
где т — Ате — he /(2 кТ) — характерный радиус взаимодействия, здесь с — скорость света в данной среде, м/с.
Далее этот объект будем называть макроячейкой. Так как макроячейка задает размеры пространства, в котором устанавливается термодинамическое равновесие, то для неё применимы соотношения равновесной термодинамики. Поэтому состояние макроячейки может быть охарактеризовано с помощью температуры, энтропии, массы количества частиц и других термодинамических параметров. Но эти параметры можно рассчитать лишь с конечной долей определенности. Возникающая неопределенность носит объективный характер и связана с тем, что макроячейки постоянно обмениваются между собой элементарными порциями количества теплоты AQ = kT.
Если макроячейка получает элементарную порцию количества теплоты AQ при P = const, то объем, температура и масса макроячейки изменяются на величину
А? = с~Ат,
(8)
где с — скорость звука в данной среде, м/с.
Известно, что уравнения массообмена, с помощью которых вводятся коэффициенты массопереноса, не имеют фундаментальной теоретической основы. Аппарат НВТ позволяет сформулировать основания такой теории.
Определим удельный равновесный мольный поток /-го компонента следующим образом
АГРсм
Jm =
F Ат
Уі
(9)
где р — молярная плотность смеси, кмоль/м, у — мольная доля компонента- Е-поверхность, через которую проходит поток вещества, м2.
Поверхность Е определяется как отношение характерного для упругих взаимодействий объема, А V к характерному линейному размеру А1. Тогда выражение для удельного равновесного потока для 7-го компонента запишется следующим образом
AV = кТ/к
'S'
(5)
(1G)
Если подставить значения скорости звука и плотности смеси для идеальногазовых условий, то получим следующее выражение:
yj (cP /Cy)RTM '-
(11)
Сравнение этого выражения с удельным потоком частиц, полученным в рамках статистической физики [5, с. 110]
N =
JlnRTM'-
(12)
Ju. = AJU. =
Ml Mi
dJ,
Mi
dy
AK,.
(13)
Линейное приближение предполагает, что перепад концентраций Ду мал. Оценим величину этого перепада. Для этого в некотором приграничном объеме, охватывающем как газовую, так и жидкую фазы, выделим ламинарный газовый слой толщиной 5 (рисунок). Пусть в этом слое имеет место
ЭУ,
перепад концентраций (у, — - у), (где У1 -равновесная концентрация на границе раздела фаз). С позиций НВТ нельзя считать, что концентрация /-го компонента изменяется непрерывно. В силу квантовых эффектов профиль концентраций размытый. Величину этой размытости проще всего оценить, предполагая линейный профиль концентраций в пределах малой толщины слоя 5
приводит к выводу о полной их идентичности с точностью до постоянного множителя порядка 2п, что подтверждает правильность выбора характерного линейного размера.
Перейдем непосредственно к описанию процесса переноса вещества из одной фазы в другую. В литературе по массообмену [1−3] принято считать, что при массоотдаче к турбулентному потоку или от него сопротивление переносу массы в значительной степени сосредоточено в тонком слое, прилегающем к границе раздела фаз. Примем эту гипотезу в качестве рабочей. Будем также считать, что для описания неравновесных процессов в приграничном слое можно применить линейное приближение. Тогда результирующая неравновесная плотность потока массы в линейном приближении для изотермических и изобарических условий составит
Ay,. =lP^
'- Ъу! М'-
(14)
где П-Ьу! М — число элементарных ступенек, каждой из которых соответствует перепад Ду/.
Условия на межфазной границе
Подставляя выражение равновесного мольного потока (10) в (13) и учитывая (14), получим уравнение для неравновесного мольного потока
дс
CSyPcMy РсМуУі
Sy dpi

CMy
дУі
S, Ac
(15)
Представим мольный поток изведения плотности потока
гок в виде про-^& gt-уі(Уі-Уі) на
Вместо скорости переноса субстанции, н удобно ввести безразмерную характери-некоторую скорость переноса^*с'-у, которую стику процесса переноса массы — коэффициент переноса V, который определяется выражением
будем называть скоростью переноса субстанции
JMyi = ^уі(Уі-Уі)^
ССу ~ С5уАТ / 8^.
Фи — РсМу
РемуУі Эс
''-Sy
dyt
Sy, dpCMy
+ Уі
dyt
(16)
(17)
(18)
Vy=CCylcSy
(19)
Тогда величина неравновесного мольного потока со стороны газовой фазы окончательно запишется следующим образом:
1ц*=УуСя& amp-уі{Уі-Уі) —
(2G)
В FUNDAMENTAL RESEARCH № 1G, 2G13 В
На границе раздела фаз допускается существование термодинамического равновесия. Это условие определяется выражением
у*=к, хп (21)
где к. — константа фазового равновесия для 1-го компонента.
Учитывая условие фазового равновесия (21) и записав выражение для удельного мольного потока в газовой фазе (20) и аналогичное ему выражение для жидкой фазы, получим уравнение для удельного мольного потока через границу раздела фаз
мУ1 = & amp- л (у" - к’Х& gt-)'
А-
ка Р* + Р*- '
Р* = УуС8у§/'
Р- =
где кр. — коэффициент массопередачи /-го компонента- в ., в. — коэффициенты мас-соотдачи /-го компонента соответственно в газовой и жидкой фазах. Величина ф. (ф.), входящая в уравнение для коэффициента массоотдачи, может быть рассчитана теоретически (см. уравнение (18)). Она зависит от физических свойств смеси газа или жидкости (скорости звука, плотности, молярной массы). Безразмерный коэффициент переноса субстанции V (у) характеризует гидродинамические условия, в которых протекает процесс массопередачи и теоретически в рамках поставленной задачи определен быть не может. В общем случае коэффициент переноса V меньше единицы. В пределе он может стать равным единице, и тогда °су = %, 5у = Д1, то есть процесс переноса массы будет происходить с максимально возможной скоростью, равной скорости звука в данной среде, а толщина пограничного ламинарного слоя 5 будет равна характерному линейному размеру Д1.
При использовании уравнения массо-передачи (22) для расчета массообменного процесса коэффициентами переноса V и V необходимо задаваться и корректировать их в зависимости от того, какое распределение концентраций по высоте аппарата необходимо получить.
Таким образом, в рамках нелокальной версии термодинамики удалось раскрыть смысл коэффициентов мас-соотдачи, выделить в них физическую и гидродинамическую составляющие. Это
(22)
(23)
(24)
(25)
позволит в дальнейшем по предложенной модели сформировать методику расчета массообменного аппарата, в котором протекает процесс однокомпонентной или многокомпонентной абсорбции. Как видно из предложенной модели, для каждого ком -понента будет своя плотность потока в газовой (у ~yt) или жидкой ф^. (*, — л: *)
фазе и один для всех компонентов коэффициент переноса v или v который отвечает за гидродинамическую обстановку, существующую в данных условиях массопере-носа. Алгоритм расчета по предложенной модели и результаты расчета однокомпонентной и многокомпонентной изотермической абсорбции будут изложены в следующем сообщении.
Список литературы
1. Александров И. А. Массопередача при ректификации и абсорбции могокомпонентных смесей. — Л.: Химия, 1974. — 320 с.
2. Кафаров В. В. Основы массопередачи. — М.: Высшая школа, 1979. — 439 с.
3. Майков В. П. Расширенная версия классической термодинамики — физика дискретного пространства-времени. -М.: МГУИЭ. 1997 — 160 с., ил.
4. Трайбус М. Термостатика и термодинамика. — М.: Энергия, 1974. — 504 с.
5. Шервуд Т., Пигфорд Р., Уилки Ч. Массопередача. -М.: Химия, 1982. — 696 с.
References
1. Aleksandrov I.A. Masstransfer during rectification and absorption of multicomponent mixtures. Leningrad: Chemistry, 1974. 320 p.
2. Kafarov V.V. Bases of masstransfer. Moscow.: Higher school, 1979. 439 p.
3. Maikov V.P. The extended version of classic thermodynamics is physics of discrete space-time. Moscow: MGUIE. 1997. 160 p.
4. Traibus M. Thermostatics and thermodynamics. Moscow: Energy, 1974. 504 p.
5. Sherwood T., Pigford R., Wilke Ch. Masstransfer, Moscow: Chemistry, 1982. 696 p.
Рецензенты:
Фокин В. М., д.т.н., профессор, заведующий кафедрой «Энергоснабжение и теплотехника», ФГБОУ ВПО «Волгоградский государственный архитектурно-строительный университет», г. Волгоград-
Семенов Б. А., д.т.н., профессор, заведующий кафедрой «Промышленная теплотехника», ФГБОУ ВПО «Саратовский государственный технический университет им. Ю.А. Гагарина», г. Саратов.
Работа поступила в редакцию 10. 09. 2013.

ПоказатьСвернуть
Заполнить форму текущей работой