ФестнщФТУММ/Proceedings of VSUET ISSN 2226-910X E-ISSN 2310-1202
DOI: http://doi.org/1Q.2Q914/231Q-12Q2-202Q-1-23Q-236_Оригинальная статья/Research article_
УДК 661.935, 519.633.2_Open Access Available online at vestnik-vsuet.ru
Численный анализ циклического процесса адсорбционного _разделения синтез-газа_
Евгений И. Акулинин 1 akulinin-2QQ6@yandex.ru QQQQ-QQQ2-3816-62Q4 Олег О. Голубятников 1 golubyatnikov_ol@mail.ru 0000-0002-0530-0964 Дмитрий С. Дворецкий 1 dvoretsky@tambov.ru 0000-0002-4352-810Х
_Станислав И. Дворецкий 1 sdvoretsky@mail.tstu.ru 0000-0002-7571-3748
1 Тамбовский государственный технический университет, ул. Советская, 1Q6, г. Тамбов, 392QQQ, Россия_
Аннотация. Проведено численное исследование динамики циклического процесса короткоцикловой безнагревной адсорбции при разделении синтез-газа и концентрировании водорода. Изучено влияние режимных (давления на выходе компрессора и вакуум-насоса, длительности цикла «адсорбция-десорбция», коэффициента обратного потока) и конструктивных (высоты насыпного слоя адсорбента и внутреннего диаметра адсорбера) параметров установки на степень извлечения и чистоту концентрируемого водорода, производительность установки в заданном диапазоне изменения температуры, состав и давление исходной газовой смеси. Установлено, что при уменьшении содержания водорода от 68 до 48 об. % и одновременном увеличении содержания диоксида углерода в исходной газовой смеси от 27 до 47 об. % требуемое значение чистоты водорода на уровне 99.99 об. % может быть достигнуто за счет уменьшения давления на выходе вакуум-насоса от Q,75^1Q5 до Q,5^1Q5 Па. При этом увеличивается длительность стадии адсорбции от 12Q до 15Q с, а степень извлечения водорода снижается от 55 до 52% вследствие увеличения доли потока, отбираемого на регенерацию адсорбента. Повышение температуры исходной газовой смеси от 293 до 323 К приводит к необходимости уменьшения длительности стадии адсорбции от 148 до 42 с из-за уменьшения равновесных концентраций диоксида и
моноксида углерода в адсорбенте._
Ключевые слова безнагревная адсорбция, цеолитовый адсорбент NaX, синтез-газ, водород, динамика, математическое моделирование, вычислительный эксперимент.
Numerical study of the dynamics of synthesis gas adsorption _separation_
Evgeny I. Akulinin 1 akulinin-2QQ6@yandex.ru QQQQ-QQQ2-3816-62Q4 Oleg O. Golubyatnikov 1 golubyatnikov_ol@mail.ru 0000-0002-0530-0964 Dmitry S. Dvoretsky 1 dvoretsky@tambov.ru 0000-0002-4352-810X
_Stanislav I. Dvoretsky 1 sdvoretsky@mail.tstu.ru 0000-0002-7571-3748
1 Tambov State Technical University, Sovetskaya str., 1Q6, Tambov, 392QQQ, Russia_
Abstract. A numerical study of the dynamics of pressure swing adsorption process for synthesis gas separation and hydrogen concentration was carried out. The influence of regime parameters (compressor and vacuum pump outlet pressures, duration of the "adsorption-desorption" cycle, backflow coefficient) and design parameters (height of the adsorbent bulk layer and the internal bed diameter) on the recovery and purity of hydrogen, plant performance in a given range of changes in the temperature, composition and pressure of the initial gas mixture was studied. It is found that when the hydrogen content decreases from 68 to 48 vol. % and a simultaneous increase in the carbon dioxide content in the initial gas mixture from 27 to 47 vol. % required hydrogen purity value at 99.99 vol. % can be achieved by reducing the pressure at the outlet of the vacuum pump from Q.75^1Q5 to Q.5^1Q5 Pa. At the same time, the duration of the adsorption stage increases from 12Q to 15Q seconds, and the degree of hydrogen extraction decreases from 55 to 52% due to an increase in the proportion of the flow selected for the regeneration of the adsorbent. Increasing the temperature of the initial gas mixture from 293 to 323 K leads to the need to reduce the duration of the adsorption stage from 148 to 42 s due to a decrease
in the equilibrium concentrations of carbon dioxide and monoxide in the adsorbent.Click to enter "Abstract"_
Keywords: heatless adsorption, zeolite adsorbent 13X, hydrogen, synthesis, dynamic, mathematical modeling, simulation.
Введение
Циклические адсорбционные процессы и, в частности, короткоцикловая безнагревная адсорбция (КБА) являются наиболее распространенным способом разделения газовых смесей и концентрирования в них целевых продуктов. Процессы КБА не предполагают наличия
Для цитирования Акулинин Е.И., Голубятников О.О., Дворецкий Д.С., Дворецкий С.И. Численный анализ циклического процесса адсорбционного разделения синтез-газа // Вестник ВГУИТ. 2020. Т. 82. № 1. С. 230-236. doi:10.20914/2310-1202-2020-1-230-236
внешнего источника теплоты, а теплота, выделяемая на стадии адсорбции, используется для осуществления процесса десорбции адсорбтива из адсорбента (регенерации адсорбента).
Циклические адсорбционные процессы широко применяют в промышленности для разделения газовых смесей и концентрирования
For citation
Akulinin E.I., Golubyatnikov O.O., Dvoretsky D.S., Dvoretsky S.I.
Numerical study of the dynamics of synthesis gas adsorption separation.
Vestnik VGUIT [Proceedings of VSUET]. 2Q2Q. vol. 82. no. 1. pp. 23Q-236.
(in Russian). doi:1Q.2Q914/231Q-12Q2-2Q2Q-1-23Q-236_
© 2020, Акулинин Е.И. и др. / Akulinin E.I. et al.
This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License
водорода, кислорода, азота, диоксида углерода и других газов. Одной из актуальных задач в области адсорбционного разделения является извлечение водорода с чистотой 99.99 об. % из водородсодержащих технологических потоков и синтез-газа [1-4].
Цель работы - численное исследование системных связей и закономерностей циклических адсорбционных процессов разделения синтез-газа и концентрирования водорода, влияния конструктивных и режимных параметров установки КБА на степень извлечения и чистоту водорода в заданном диапазоне изменения возмущающих воздействий (температуры, состава и давления исходной водородсодержащей смеси).
Описание циклического процесса адсорбционного разделения синтез-газа и концентрирования водорода
Технологический процесс концентрирования водорода методом адсорбционного разделения синтез-газа с циклически изменяющимся давлением осуществляется в четырехадсорберной установке КБА с гранулированным адсорбентом NaX (рисунок 1) [5].
Установка предназначена для концентрирования водорода до 99.99 об. %. Исходная газовая смесь подается в установку после предварительной осушки (стадия осушки в работе не рассматривается). Подъем давления в адсорберах 2а, Ь, с, d производится открытием управляемых клапанов 3а, Ь, с, d, через которые газовая смесь подводится к слою адсорбента в адсорберах на стадии адсорбции. Продукционный водород отводится из адсорберов на стадии адсорбции через клапаны 5а, Ь, с, d и направляется потребителю через ресивер (на рисунке 1 не показан). Противоточная регенерация слоя адсорбента осуществляется путем пропускания продувочной (обогащенной водородом) газовой смеси через клапаны 6а, Ь, с, d, а через клапаны 4а, Ь, с, d осуществляется отвод продувочного газа на стадии десорбции (регенерации адсорбента).
Согласованная работа клапанов установи КБА (рисунок 1) осуществляет последовательно в каждом адсорбере цикл «адсорбция-десорбция». Циклограмма изменения давления в одном адсорбере в течение одного цикла продолжительностью 400 с приведена на рисунке 2.
В качестве адсорбентов в установках КБА при осуществлении циклических адсорбционных процессов широко используют гранулированные активные угли и цеолиты СаА, NaX, обладающие наибольшей емкостью и селективностью по СО2 и СО [3, 4].
Продуктовый поток
Обратный поток
Рисунок 1. Технологическая схема установки КБА разделения синтез-газа и концентрирования водорода
Figure 1. Technological scheme of the PSA unit for synthesis gas separation and hydrogen production
'ads
pn
'ads
18
16
14
12
10
8
6
4
2
10s Па
105 РЭ
1 fei 1 'des 1 'eq i
-.--- -■-1-— f—-1
0 50 100 150 200 250 300
350 f, f,
Рисунок 2. Циклограмма изменения давления в адсорбере установки КБА в течение одного цикла «адсорбция-десорбция»
Figure 2. Cyclogram of dynamics pressure in the adsorber of the PSA unit during one cycle "adsorption-desorption"
При адсорбции (десорбции) H2, CO2 и CO цеолитовым адсорбентом в адсорберах 2а, б, в, г установки КБА (рисунок 1) протекают следующие массо- и теплообменные процессы: а) диффузия H2, CO2 и CO в потоке газовой смеси; б) массообмен H2, CO2, CO и теплообмен между газовой фазой и адсорбентом; в) адсорбция H2, CO2, CO на поверхности и в микропорах гранул цеолитового адсорбента с выделением теплоты и десорбция H2, CO2, CO из микропор и с поверхности гранул с поглощением теплоты.
Анализ результатов физического моделирования показал, что диффузия H2, CO2, CO и распространение теплоты в газовом потоке и гранулированном адсорбенте осуществляется в основном в продольном направлении относительно движения потока газовой смеси в адсорбере (по высоте слоя адсорбента). При этом процесс адсорбции преимущественно CO2 и CO
адсорбентом осуществляется в смешанно-диффузионной области и определяется коэффициентами внешней массоотдачи и внутренней диффузии, скоростью газовой смеси в слое адсорбента, а также равновесными соотношениями концентраций Н2, СО2 и СО в фазах.
При численном анализе процесса концентрирования водорода принимали следующие допущения: 1) исходная газовая смесь - синтез-газ содержит 1 - Н2 с концентрацией 48-68 об. %, 2 - СО2 с концентрацией 27-47 об. %, 3 - СО с концентрацией 5 об. % является 3-компонент-ной газовой смесью и рассматривается как идеальный газ, что вполне допустимо при давлении в адсорбере до 200-105 Па [6, 7]; 2) диффузия Н2, СО2, СО и распространение теплоты в газовой и твердой фазах осуществляется только в аксиальном направлении движения потока газовой смеси в адсорбере (по высоте адсорбента) [7-9]; 3) в качестве адсорбента используется гранулированный цеолит №Х с диаметром гранул 1,5 мм [4]; 4) адсорбционное равновесие описывается уравнением Дубинина-Радушкевича [10]; 5) десорбционные ветви изотерм адсорбции Н2, СО2, СО на цеолите совпадают с адсорбционными [10]; 6) температура газа в ресивере равна температуре газа на выходе из адсорбера, тепловые потери в окружающую среду пренебрежимо малы.
Математическое описание адсорбционного процесса разделения синтез-газа и концентрирования водорода включает дифференциальные
уравнения с частными производными покомпонентного материального баланса компонентов Н2, СО2, СО в потоке газовой фазы и распространения теплоты в газовой и твердой фазах; обыкновенные дифференциальные уравнения кинетики адсорбции-десорбции, изменения скорости потока газовоздушной смеси и давления газовой смеси по высоте насыпного слоя адсорбента[11-14].
Для решения системы дифференциальных уравнений с соответствующими начальными и граничными условиями использовали метод прямых в программной среде МаЙаЬ [15, 16]. Анализ точности математической модели проводили с использованием фактической средне-квадратической ошибки между рассчитанными по модели и экспериментальными значениями чистоты водорода на выходе из адсорбера. Фактическая среднеквадратическая ошибка составила 0,15 об. %, что позволяет использовать математическую модель для исследования и оптимизации циклических режимов функционирования установки КБА для разделения синтез-газа и концентрирования водорода [1-4].
Численный анализ циклического процесса адсорбционного разделения синтез-газа
Для исследования влияния режимных и конструктивных параметров установки КБА на показатели эффективности функционирования циклического процесса адсорбционного концентрирования водорода проводили вычислительные эксперименты (таблица 1).
Исходные данные для вычислительного эксперимента Source data for computational experiment
Таблица 1.
Table 1.
Варьируемые переменные Variables Диапазон изменения The range of variation Номинальное значение Nominal value
Длительность стадии адсорбции /ads, с Duration of the adsorption stage tads, s 10-180 40
Давление на выходе компрессора Pln, х 10 Па The outlet pressure of the compressor Pln, x10 Pa 5-30 17,50
Давление на выходе вакуум-насоса , х 10 Па The pressure at the outlet of vacuum-pump , x 10 Pa 0,5-1,00 0,75
Коэффициент обратного потока 9, отн. ед. | Backflow coefficient 9, [-] 1,25-15,00 5,83
Концентрация компонентов в исходной смеси у111, об.% Concentration of components in the initial mixture y™, vol.% H2 CO2 CO 48-68 27-47 5 68 47 5
Температура исходной смеси Tgm, К The temperature of the initial mixture Tgln , K 273-323 298
Отношение длины адсорбера к его диаметру L/Da The ratio of the length of the adsorber to its diameter L/Da 4-8 6
На рисунке 3 представлены зависимости чистоты продукционного водорода при различных значениях давления на выходе компрессора и длительности стадии адсорбции ^=100 с.
Рисунок 3. Динамика чистоты продукционного водорода в ресивере установки при длительности стадии адсорбции 4ds=100 с и различных значениях давления на выходе компрессора Pm: 1 - 5-105 Па; 2 - 17,5405 Па; 3 - 30-105 Па
Figure 3. Dynamics of the purity of production hydrogen in the unit receiver at the duration of the adsorption stage 4ds =100 s and different values of pressure at the compressor outlet Pm: 1 - 5-105 Pa; 2 - 17.5-105 Pa; 3 - 30-105 Pa
Анализ рисунка 3 показывает, что увеличение давления на выходе компрессора от 5-105 до 30-105 Па способствует ускоренному выходу установки на стационарный периодический режим от ~30 до ~10 циклов за счет повышения перепада давлений на впускном клапане данного адсорбера и увеличения расхода исходной газовой смеси на входе в адсорбер установки на стадии адсорбции.
об.% vol,%o 99,99^
99,88 99,77 99,66 99,55 99,44 99,33
G0Llt, нл/мин Gou\ l/min STP
г 60
Рисунок 4. Зависимость чистоты продуктового водорода >>1out , степени его извлечения п и производительности установки Gout при номинальном режиме работы установки
Figure 4. The dependence of purity of product hydrogen y1out , recovery П and capacity Gout at a nominal regime mode
Анализ графиков зависимостей чистоты продуктового водорода y1out , степени его извлечения п и производительности установки Gout от продолжительности стадии адсорбции tads при номинальном режиме ее работы (рисунок 4, 5) показывает, что увеличение давления на входе в установку в 1,7 раза приводит к снижению степени извлечения водорода на ~10% и увеличению производительности в 2,5 раза. Этот факт объясняется увеличением равновесной концентрации извлекаемых компонентов в адсорбенте при увеличении давления Pm от 5-105 до 30-105 Па.
y^ut, об.% у™1, vol.% 99,99
99,88
99,77
99,66
99,55
99,44
99,33
Gm\ нл/мин Gm", l/min STP 1000
a b
Рисунок 5. Зависимость чистоты продуктового водорода у"* , степени его извлечения п и производительности установки Gout при давлении на выходе компрессора Pm: a - 5-105 Па; b - 30-105 Па
Figure 5. Dependence of the purity of the product hydrogen y1out , recovery п and capacity Gout at the compressor outlet pressure Pm: a - 5-105 Pa; b - 30-105 Pa
Повышение температуры исходной газовой смеси от 293 до 323 К приводит к необходимости уменьшения длительности стадии адсорбции
от 148 до 42 с (рисунок 6, а) из-за уменьшения равновесных концентраций диоксида и монок-сида углерода в адсорбенте.
При уменьшении содержания водорода от 68 до 48 об. % и одновременном увеличении содержания диоксида углерода в исходной газовой смеси от 27 до 47 об. % адсорбционная емкость исчерпывается быстрее, а требуемое значение чистоты водорода на уровне 99.99 об. % не достигается (рисунок 6, Ь). При этом производительность установки снижается на ~25%. Расчеты показывают, что увеличение давления на выходе компрессора от 17,5-Ш5 до 23,5-105 Па позволяет повысить производительность на ~25%, однако степень извлечения водорода при этом снижается относительно номинального режима от 55 до 50%.
Увеличение давления на выходе вакуум-насоса от 0,75 • 105 до 1 • 105 Па приводит к уменьшению: длительности стадии адсорбции
уТ\ 0б.%
99,99- !li 60
99,88 50
99,77 40-
99,66- 30
99,55 20
99,44- 10
99,33 0
%
40
80
120 160
нл/мин l/min STP 450
400 350 300 250 200 150 100 50 0
.. с
„, S
от 120 до 50 с, степени извлечения водорода от 55 до 48% (рисунок 7, Ь). Уменьшение давления на выходе вакуум-насоса от 0,75-105 до 0,5^ 105 Па приводит к повышению длительности стадии адсорбции от 120 до 150 с, при которой обеспечивается требуемая чистота продуктового водорода у™' =99.99 об. %. Степень извлечения водорода при этом снижается от 55 до 52% вследствие увеличения доли потока, отбираемого на регенерацию адсорбента (рисунок 7, а). Расчеты показывают, что для увеличения степени извлечения до 55% (при сохранении чистоты водорода 99.99 об. %) в данном режиме необходимо уменьшить коэффициент обратного потока на ~5% относительно его номинального значения.
Got", нл/мин G^', i/mrn STP 450
vol.%
99,99-, Tb
60-
99,00 50-
98,01 ■ 40
30-
97,02
20-
96,03 10-
95,04 0-
Рисунок 6. Зависимость чистоты продуктового водорода установки Gout при: а - тш =323 К; Ь - у™ =48 об.%
, степени его извлечения п и производительности
Figure 6. Dependence of the purity of the product hydrogen
y10Ut , recovery п and capacity Gout at: a - Tgin=323 K;
=48 vol.%
y^ut, об.% vol.%
99,99 n. 60-
99,98 50
99,97 40
99,96 30-
99,95 20-
99,94 10
99,93 0
%
G°ut, нл/мин У1и1, 06.% GDL", l/min STP Kul, vol.%
40
80 120 160
-450 99,99- '1. 60-
400
-350 50-
-300 40-
250 200 99,00 30-
150 20
100
-50 10-
■0 98,01 oJ
^ads ' c ^ads ■ s
Goul, нл/мин G0l,t, l/min STP -450
Рисунок 7. Зависимость чистоты продуктового водорода установки Gout при р£8 : а - 0,5 • 105 Па; Ь - 1,0 • 105 Па
y° , степени его извлечения п и производительности
Figure 7. Dependence of the purity of the product hydrogen b - 1,0 • 105 Pa
y1out , recovery п and capacity Gout at PdieI1s: a - 0,5 • 105 Pa;
b
a
b
a
Увеличение длины насыпного слоя адсорбента от 0,9 до 1,2 м приводит к снижению степени извлечения водорода в среднем на 5% вследствие увеличения перепада давлений в слое адсорбента (рисунок 8, Ь). Уменьшение длины адсорбера от 0,9 до 0,6 м, напротив,
нл/мин l/min STP 450
приводит к повышению степени извлечения водорода в среднем на 2%, при этом чистота водорода не достигается на требуемом уровне 99.99 об. % вследствие быстрого исчерпания адсорбционной емкости адсорбента (рисунок 8, а).
y°ut, vol.%
99,76- 'It бО-
99,18
98,60- бО-
98.02- 40-
97,44- 30-
96,86
96,28 20-
95,70 10-
95,12 0-
vol, %
99,99- n,
60
99,88
50-
99,77
40-
99,66-
99,55- 30
99,44 20-
99,33- 10
99,22 0
%
40
80
G°ut, нл/мин G°ul, l/min STP 450 400 1-350 300 250 1-200 150 100 50 0
I г с
, s
120 160
b
Рисунок 8. Зависимость чистоты продуктового водорода у"- , степени его извлечения п и производительности установки Gout при: а - L/DA=4; b - L/DA=8
Figure 8. Dependence of the purity of the product hydrogen у"- , recovery п and capacity Gout at: a - L/Da=4; b - L/Da=8
Заключение
На основе современных методов системного анализа, математического моделирования и вычислительного эксперимента получены новые результаты для теории и практики создания технологических процессов адсорбционного разделения многокомпонентных газовых смесей с циклически изменяющимся давлением.
Определены наиболее опасные возмущающие воздействия (состав и температура исходной водородсодержащей смеси), наиболее эффективные управляющие воздействия (длительность
стадии адсорбции, давление на выходе компрессора и вакуум-насоса, коэффициент обратного потока) и области допустимых режимов, необходимые для постановки и решения задачи оптимизации стационарного периодического режима функционирования установки КБА при разделении синтез-газа и концентрировании водорода.
Благодарности
Работа выполнена при финансовой поддержке Министерства науки и высшего образования России в рамках гранта Президента РФ МК-1604.2020.8.
Литература
1 Ohs B., Lohaus J., Marten D., Hannemann-Tamas R. et al. Optimized hollow fiber sorbents and pressure swing adsorption
process for H-2 recovery //Ind. Eng. Chem. Res. 2018. V. 57. № 14. P. 5093-5105.
2 Agarwal A., Biegler L., Zitiiey S. Superstructure-based optimal synthesis of pressure swing adsorption cycles for precombustion C02 Capture // Ind. Eng. Chem. Res. 2010. V. 49. № 11. P. 5066-5079.
3 Arora A., Bajaj I., Iyer S., Hasan M. Optimal synthesis of periodic sorption enhanced reaction processes with application to hydrogen production//Сотр. Chem. Eng. 2018. V. 115. P. 89-111.
4 Ben-Mansour R., Habib M.A., Bamidele O.E., Basha M. et al. Carbon capture by physical adsorption: Materials, experimental investigations and numerical modeling and simulations - A review // Applied Energy. 2016. V. 161. P. 225-255.
5 Ruthven D.M., Farooq S., Knaebel K.S. Pressure swing adsorption. New York, 1993. 376 p.
6 Кириллин B.A., Сычев B.B., Шейндпин A.E. Техническая термодинамика. Москва: Издательский дом МЭИ, 2016. 496 с.
7 Shafeeyan M.S., Wan Daud W.M.A., Shamiri A. A review of mathematical modeling of fixed-bed columns for carbon dioxide adsorption // Chem. Eng. Res. Des. 2014. V. 92. № 5. P. 961-988.
8 Tao W., Ma S., Xiao J., Benard P. et al. Simulation and optimization for hydrogen purification performance of vacuum pressure swing adsorption // Energy Procedia. 2019. V. 158. P. 1917-1923.
9 Han Z.-Y., Xing R., Zhang D.-H., Shen Y.-H. et al. Vacuum pressure swing adsorption system for N2/CH4 separation under uncertainty // Chem. Eng. Res. and Des. 2019. V. 142. P. 245-256.
10 Дубинин M.M. Адсорбция и пористость. Москва: Изд. ВАХЗ, 1972. 123 с.
11 Tavan Y., Hosseini S.H., Olazar М. A note on an integrated process of methane steam reforming injunction with pressure-swing adsorption to produce pure hydrogen: Mathematical modeling // Ind. Eng. Chem. Res. 2015. V. 54. №№ 51. P. 12937-12947.
12 Akulinin E.I., Golubyatnikov O.O., Dvoretsky D.S., Dvoretsky S.I. Numerical study of cyclic adsorption processes of air oxygen enrichment in dynamics // Journal of Physics: Conference Series. 2019. V. 1278. № 1. 012005.
13 Akulinin E.I., Golubyatnikov O.O., Dvoretsky D.S., Dvoretsky S.I. Optimizing pressure-swing adsorption processes and installations for gas mixture purification and separation // Chemical Engineering Transactions. 2019. V. 74. P. 883-888.
a
14 Grossmann I.E., Apap R.M., Calfa B.A., Garcia-Herreros P. et al. Mathematical Programming Techniques for Optimization under Uncertainty and Their Application in Process Systems Engineering // Theoretical Foundations of Chemical Engineering. 2017.
V. 51. №6. P. 893-909.
15 Rice R.G., Do D.D. Applied Mathematics and Modeling for Chemical Engineers; 2 ed. New Jersey, 2012.
16 Гольдштейн А. Л. Оптимизация в среде MATLAB. Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2015 1 92 с.
References
1 Ohs B., Lohaus J., Marten D., Hannemann-Tamas R. et al. Optimized hollow fiber sorbents and pressure swing adsorption process for H 2 recovery. Ind. Eng. Chem. Res. 2018. vol. 57. no. 14. pp. 5093-5105.
2 Agarwal A., Biegler L., Zitney S. Superstructure-based optimal synthesis of pressure swing adsorption cycles for precombustion CO2 Capture. Ind. Eng. Chem. Res. 2010. vol. 49. no. 11. pp. 5066-5079.
3 Arora A., Bajaj I., Iyer S., Hasan M. Optimal synthesis of periodic sorption enhanced reaction processes with application to hydrogen production. Comp. Chem. Eng. 2018. vol. 115. pp. 89-111.
4 Ben-Mansour R., Habib M.A., Bamidele O.E., Basha M. et al. Carbon capture by physical adsorption: Materials, experimental investigations and numerical modeling and simulations - A review. Applied Energy. 2016. vol. 161. pp. 225-255.
5 Ruthven D.M., Farooq S., Knaebel K.S. Pressure swing adsorption. New York, 1993. 376 p.
6 Kirillin V.A., Sychev V.V., Sheindlin А.Е. Technical thermodynamics. Moscow, Publishing house MEI, 2016. 496 p. (in Russian).
7 Shafeeyan M.S., Wan Daud W.M.A., Shamiri A. A review of mathematical modeling of fixed-bed columns for carbon dioxide adsorption. Chem. Eng. Res. Des. 2014. vol. 92. no. 5. pp. 961-988.
8 Tao W., Ma S., Xiao J., Benard P. et al. Simulation and optimization for hydrogen purification performance of vacuum pressure swing adsorption. Energy Procedia. 2019. vol. 158. pp. 1917-1923.
9 Han Z.-Y., Xing R., Zhang D.-H., Shen Y.-H. et al. Vacuum pressure swing adsorption system for N2/CH4 separation under uncertainty. Chem. Eng. Res. and Des. 2019. vol. 142. pp. 245-256.
10 Dubinin M.M. Adsorption and porosity. Moscow, 1972. (in Russian).
11 Tavan Y., Hosseini S.H., Olazar M A note on an integrated process of methane steam reforming in junction with pressure-swing adsorption to produce pure hydrogen: Mathematical modeling. Ind. Eng. Chem. Res. 2015. vol. 54. no. 51. pp. 12937-12947.
12 Akulinin E.I., Golubyatnikov O.O., Dvoretsky D.S., Dvoretsky S.I. Numerical study of cyclic adsorption processes of air oxygen enrichment in dynamics. Journal of Physics: Conference Series. 2019. vol. 1278. no. 1. 012005.
13 Akulinin E.I., Golubyatnikov O.O., Dvoretsky D.S., Dvoretsky S.I. Optimizing pressure-swing adsorption processes and installations for gas mixture purification and separation. Chemical Engineering Transactions. 2019. vol. 74. pp. 883-888.
14 Grossmann I.E., Apap R.M., Calfa B.A., Garcia-Herreros P. et al. Mathematical Programming Techniques for Optimization under Uncertainty and Their Application in Process Systems Engineering. Theoretical Foundations of Chemical Engineering. 2017. vol. 51. no. 6. pp. 893-909.
15 Rice R.G., Do D.D. Applied Mathematics and Modeling for Chemical Engineers; 2 ed. New Jersey, 2012.
16 Goldstein A.L. Optimization in MATLAB. Perm, Publishing house Perm. nat. researched Polytechnic University, 2015. 192 p. (in Russian).
Сведения об авторах Information about authors
Евгений И. Акулинин к.т.н., доцент, кафедра технологий и Evgeny I Akulinin Cand. Sci. (Engin.), associate professor, оборудования пищевых и химических производств, Тамбовский technologies and equipment of food and chemical industries государственный университет инженерных технологий, ул. Со- department, Tambov State Technical University, Sovetskaya str., ветская, 106, г. Тамбов, 393200, Россия, akiiinin-2006(o!yandex.ru 106, Tambov, 393200, Russia, akulinin-2006®!yandex.ru
https://orcid.org/0000-0002-3816-6204 https://orcid.org/0000-0002-3816-6204
Олег О. Голубятников к.т.н., ст. пргподаватель. кафедра техно- Qleg Q ^^^ Cand Scl (Engm senlor lecturer_
логий и оборудования пищажк и химических пр°изводств, technologies and equipment of food and chemical industries Тамбовский государственный университет инженерных техно-
r £ ^ department, Tambov State Technical University, Sovetskaya str.,
логий, ул. Советская, 106, г. Тамбов, 393200, Россия, .Л. „ ' т^ , , ,
, , 106, Tambov, 393200, Russia, golubyatnikov ol@mail.ru goTubyatnikov olfffimail.ru ° - _
https://orcKforg/0000-0002-0530-0964 https://ordd.org/0000-0002-0530-0964
Дмитрий С. Дворецкий д.т.н., профессор, зав. кафедрой техно- Dmitry S. Dvoretsky Dr. Sci. (Engin.), professor , technologies
логий и оборудования пищевых и химических производств, and equipment of food and chemical industries department,
Тамбовский государственный технический университет, ул. Со- Tambov State Technical University, Sovetskaya str., 106,
ветская, 106, г. Тамбов, 393200, Россия, dvoretsky®!tambov.ru Tambov, 393200, Russia, dvoretsky®!tambov.ru
https://orcid.Org/0000-0002-4352-810X https://orcid.org/0000-0002-4352-810X
Станислав И. Дворецкий д.т.н., профессор, кафедра технологий Stanislav I Dvoretsky Dr. Sci. (Engin.), professor, technologies
и оборудования пищевых и химических производств», Тамбов- and equipment of food and chemical industries department,
ский государственный технический университет, ул. Совет- Tambov State Technical University, Sovetskaya str., 106, Tambov,
екая, 106, г. Тамбов, 393200, Россия, sdvoretsky®!mail.tstu.ru 393200, Russia, sdvoretsky®?mail.tstu.ru
https://orcid.org/0000-0002-7571-3748 https://orcid.org/0000-0002-7571-3748
Вклад авторов Contribution
Все авторы в равной степени принимали участие в написании All authors are equally involved in the writing of the manuscript and are
рукописи и несут ответственность за плагиат responsible for plagiarism
Конфликт интересов Conflict of interest
Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.
Поступила 09/12/2019 После редакции 18/12/2019 Принята в печать 16/01/2020
Received 09/12/2019 Accepted in revised 18/12/2019 Accepted 16/01/2020