УДК 66.061.351
МОДЕЛИРОВАНИЕ КИНЕТИКИ МАССОПЕРЕДАЧИ В УСЛОВИЯХ САМОПРОИЗВОЛЬНОЙ МЕЖФАЗНОЙ КОНВЕКЦИИ НА ПЛОСКОЙ ГРАНИЦЕ РАЗДЕЛА ФАЗ
Л.А. Стенин, С.А. Ермаков
Кафедра «Процессы и аппараты химической технологии»,
ФГАОУ ВПО «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина», г. Екатеринбург; leonidd123@e1.ru
Представлена членом редколлегии профессором В.И. Коноваловым
Ключевые слова и фразы: массопередача; метод Монте-Карло; самопроизвольная межфазная конвекция.
Аннотация: Представлен способ аппроксимации кинетических кривых массопередачи в условиях самопроизвольний межфазной конвекции с помощью метода Монте-Карло. С помощью численных экспериментов найдены границы его применимости.
Обозначения
с - концентрация, кмоль/м3;
Кт - коэффициент массопередачи, м/с;
L - концентрация переносимого вещества в режиме самопроизвольний межфазной конвекции в нулевой момент времени, кмоль/м3;
M - концентрация переносимого вещества в «диффузионном» режиме в нулевой момент времени, кмоль/м3;
Sуд - удельная межфазная поверхность, м-1; t - время, с;
вд - подэкспоненциальный множитель для «диффузионного» режима;
Рсмк - подэкспоненциальный множитель для режима самопроизвольний межфазной конвекции.
В данной работе показано моделирование переноса масляной кислоты из бензола в водный раствор гидроксида натрия в условиях самопроизвольной межфазной конвекции. При описании процессов массопередачи вещества через плоскую границу раздела фаз можно применять [1] двухэкспоненциальную функцию
ф) = Ье-Рсмк^ + Ме-рд^ . (1)
Предполагается, что массопередача состоит из двух процессов; каждое слагаемое в уравнении (1) описывает отдельный процесс. Первый процесс - это мас-соперенос в режиме самопроизвольной межфазной конвекции (СМК). Второй процесс - это массоперенос по диффузионно-конвективному механизму («диффузионный» режим). Оба процесса идут одновременно.
Поиск решения осуществлялся методом Монте-Карло. Для каждого эксперимента было проведено 700 000 бросков случайных чисел Ь, M, рсмк и рд, рассчитывалась сумма квадратов отклонения расчетной кривой от экспериментальной. Результатом работы программы является минимальная сумма квадратов отклонений расчетной кривой от экспериментальной и соответствующие ей коэффициенты L, M, Рсмк и Рд.
В настоящей работе была поставлена задача: оценить надежность аппроксимации экспериментальных данных и установить границы применимости данного метода. Для этого использовались коэффициенты Ь, М, рсмк и рд, экспериментально полученные при массопереносе масляной кислоты концентрацией 1,58 кмоль/м3 из бензола в водный раствор гидроксида натрия концентрацией 1,56 кмоль/м3, так как при этой концентрации был проведен эксперимент. Коэффициенты Ь = 1,04 кмоль/м3, М = 0,54 кмоль/м3, рсмк = 1,05-10 4 м/с, Рд = 1,01-10-5 м/с были приняты за «эталон». С их помощью была построена эталонная кривая изменения концентрации в отдающей фазе от времени. Кривая была построена по тринадцати точкам, что согласовано с экспериментальными данными. Затем к каждой точке эталонной кривой случайным образом добавлялась абсолютная погрешность в диапазонах ±2, ±4, ±6, ±8 и ±10 %. Таким способом проводилось «зашумление» эталонных данных (рис. 1).
Затем с помощью программы были найдены коэффициенты Ь, М, рсмк и Рд и суммы квадратов отклонений. Данные представлены в табл. 1.
Из табл. 1 можно видеть, что при абсолютной погрешности экспериментальных данных ± 6 % погрешность расчета коэффициентов Ь, М, рсмк и рд не превышает ± 5 %. Показателем достоверности аппроксимации может служить сумма квадратов отклонений расчетной и экспериментальной кривых. При удовлетворительной достоверности аппроксимации, как видно из табл. 1, сумма может быть
Рис. 1. Моделирование зависимости концентрации масляной кислоты в отдающей фазе (бензоле) от времени:
---- эталонная кривая;----- кривая, полученная добавлением
к эталонной кривой абсолютной погрешности в диапазоне ± 10 %
Таблица 1
Степень отклонения расчетных коэффициентов от эталонных при различной степени «зашумления» эталонных данных
Степень «зашумления» (±), % AL, % AM, % АРсмк, % Авд, % Сумма квадратов отклонений
2 0,96 1,85 1,49 1,76 9,24-10-4
4 1,44 2,78 3,90 3,41 1,0810-3
6 2,40 4,63 4,17 4,98 5,73-10-3
8 18,3 35,2 24,3 59,8 1,03-10-2
10 6,73 13,0 8,08 24,4 7,99-10-3
Таблица 2
Отклонение массового потока первого и второго процессов в зависимости от погрешности, прибавляемой к эталонной кривой
Степень «зашумления» эталонной кривой (±), % Отклонение массового потока, %
1 процесс 2 процесс
2 0,55 2,96
4 0,83 5,40
6 2,28 3,66
8 13,6 99,5
10 3,7 32,0
порядка 10-3, но не выше 5,73-10-3. При суммарном влиянии различных факторов (абсолютной погрешности) на процесс переноса вещества более 8-10 % отклонения коэффициентов значительны и сумма квадратов отклонений велика - от 8-10-3 до значений порядка 10-2.
В присутствии быстрой химической реакции массопередача осуществляется из лимитирующей фазы, и коэффициент массопередачи рассчитывают по формуле (2)
km (2)
m dlS„ c(t) ' '
Интересно было оценить отклонение массового потока (площадь под кривой в координатах коэффициент массопередачи - движущая сила) в зависимости от степени «зашумления» эталонной кривой. Площадь под кривой во всех случаях была найдена численным интегрированием по рассчитанным коэффициентам L, M, рсмк и рд. Площадь под эталонной кривой для первого процесса равна 1,362-10-4 кмоль/(м2-с). Для второго процесса - 6,77-10-5 кмоль/(м2-с).
Таким образом, из табл. 2 можно сделать выводы о границах применимости аппроксимации двухэкспоненциальной кривой с помощью метода Монте-Карло. Данный метод применим при абсолютной погрешности экспериментальных данных не более ± 6 %. При абсолютной погрешности до ± 6 % массовые потоки в режиме СМК и в диффузионно-конвективном режиме можно рассчитать с погрешностью не более ± 4 %. Основным ориентиром, свидетельствующим о величине погрешности, может служить сумма квадратов отклонения расчетных данных от эталонной кривой. Полученная информация может быть полезна при исследованиях массопередачи, сопровождаемой быстрой химической реакцией в условиях СМК.
Список литературы
1. Вайсов, Д.В. Аппроксимация двухэкспоненциальной кривой экспериментальных данных процесса массопередачи с химической реакцией в условиях межфазной конвекции / Д.В. Вайсов, С.А. Ермаков // Тр. СвердНИИхиммаш. Сер. Оборудование для оснащения технолог. производств. - 2000. - № 7. - С. 134-141.
Modeling of Kinetics of the Mass Transfer under Conditions of Spontaneous Interfacial Turbulence on the Plane Interfacial Boundary
L.A. Stenin, S.A. Ermakov
Department “Processes and Devices of Chemical Engineering",
Ural Federal University named after First President of Russia B.N. Eltsin, Ekaterinburg; leoniddl23@e1 .ru
Key words and phrases: mass transfer; Monte-Carlo method; spontaneous interfacial turbulence.
Abstract: The paper presents the approximation technique of kinetic curves of mass transfer using Monte-Carlo method in conditions of spontaneous interfacial turbulence. Its applicability boundaries are determined via numerical experiments.
Modellierung der Kinetik des Massenübergang in den Bedingungen der spontanen zwischenphasischen Konvektion auf der flachen Grenze
der Phasenteilung
Zusammenfassung: Es ist das Verfahren der Approximation der kinetischen Kurven des Massenübergang bei den Bedingungen der spontanen zwischenphasischen Konvektion mit Hilfe der Methode von Monte-Karlo dargelegt. Mit Hilfe der Zahlexperimenten sind die Grenze seiner Benutzung gefunden.
Modélage de la cinétique du transfert de masse dans les conditions de la convection spontannée interfaciale sur une interface plate des phases
Résumé: Est présenté le moyen de l’approximation des courbes cinétiques du transfert de masse dans les conditions de la convection spontannée interfaciale à l’aide de la méthode Monte-Carlo. A l’aide des expériments numériques sont trouvées les limites de son application.
Авторы: Стенин Леонид Алексеевич - ведущий инженер кафедры «Процессы и аппараты химической технологии»; Ермаков Сергей Анатольевич -доктор технических наук, профессор кафедры «Процессы и аппараты химической технологии», ФГАОУ ВПО «УрФУ им. первого Президента России Б.Н. Ельцина», г. Екатеринбург.
Рецензент: Ермаков Анатолий Александрович - доктор технических наук, профессор, заведующий кафедрой «Процессы и аппараты химической технологии», ФГАОУ ВПО «УрФУ им. первого Президента России Б.Н. Ельцина», г. Екатеринбург.