УДК. 664.002.05 ББК 65.9(2)304.25 М-52
Меретуков Заур Айдамирович, кандидат технических наук, доцент кафедры технологии, машин и оборудования пищевых производств ФГБОУ ВПО «Майкопский государственный технологический университет», e-mail: zamer@radnet. ru;
Кошевой Евгений Пантелеевич, доктор технических наук, профессор, «Заслуженный деятель науки Российской Федерации», заведующий кафедрой машины и аппараты пищевых производств Кубанского государственного технологического университета, e-mail: Koshevoi@kubstu. ru;
Косачев Вячеслав Степанович, доктор технических наук, профессор кафедры машин и аппаратов пищевых производств факультета машиностроения и автосервиса Кубанского государственного технологического университета.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МАССОПЕРЕНОСА ПРИ СВЕРХКРИТИЧЕСКОМ ЭКСТРАГИРОВАНИИ МАСЛИЧНЫХ МАТЕРИАЛОВ ДВУОКИСЬЮ УГЛЕРОДА
(рецензирована)
Целью работы является получение математической модели экстрагирования слоя растительного масличного материала (подсолнечника) с учетом особенностей растворения масла в двуокиси углерода при сверхкритических условиях.
Ключевые слова: растительные масла, экстракция, сверхкритический СО2, растворимость, масличность.
Meretukov Zaur Aydamirovich, Candidate of Technical Sciences, assistant professor of the Department of Technology, Machinery and Equipment of Food Production of FSBEIHPE "Maikop State Technological University», e-mail: zamer@radnet.ru;
Koshevoi Eugene Panteleevich, Doctor of Technical Sciences, professor, "Honored Scientist of the Russian Federation", head of the Department of Machines and Equipment for Food Production, Kuban State Technological University, e-mail: Koshevoi@kubstu.ru;
Kosachev Vyacheslav Stepanovich, Doctor of Technical Sciences, professor of the Department of Machines and Equipment for Food Production of the Faculty of Machine-building and Service of the Kuban State Technological University.
MATHEMATICAL MODEL OF MASS TRANSFER IN SUPERCRITICAL EXTRACTION OF OIL MATERIALS WITH CARBON DIOXIDE
(reviewed)
The aim of the research has been to obtain a mathematical model of extraction of vegetable oil seed layer mate-rial (sunflower) taking into consideration features of dissolution of oil in carbon dioxide under supercritical conditions.
Keywords: vegetable oils, extraction, supercritical CO2, solubility, oil percentage.
Известен целый ряд работ, в которых исследовалась экстракция масличного материала семян подсолнечника, подробно исследовалась кинетика экстракции семян подсолнечника двуокисью углерода при сверхкритических условиях [1]. Были получены экспериментальные данные для экстракции, проводимой при давлениях 20, 30, 40, 50 и 60 МПа; температурах 313, 333 и 353 K, расходах двуокиси углерода 1-4, и 6 cм3 CO2 минута-1; средние диаметры частиц 0.23, 0.55, 1.09, 2.18 мм. В каждом эксперименте известная масса сухих молотых семян была приблизительно 4 г была помещена в колонку для экстракции неподвижного слоя объемом 10 мл (ячейка из нержавеющей стали высотой 57 мм и внутренним диаметром 20 мм), сверху и снизу ячейки были помещены фильтры с размером отверстий 0.5 мм для предотвращения выноса частиц.
Как отмечено во всех работах по экстракции высокомасличных материалов и, в том числе семян подсолнечника [1-8], начальный период процесса лимитируется растворимостью экстрагируемого масла в сверхкритическом СО2 и характеризуется постоянной скоростью. Представленные в работе [1] опыты по кинетике экстракции в зависимости от давления (20, 30, 40, 50 и 60 MQa) и температуры (313, 333 и 353) К могут быть обобщены от безразмерного времени N*t (рисунок 1), где N постоянная скорость экстракции в начальном периоде.
• Еехр Д Emod
Рис. 1 - Обобщенные зависимости по кинетике экстракции при давлениях 20, 30, 40, 50 и 60 МПа и температурах 313 К, 333 К и 353К. Все точки (эксперимент) внутри треугольников (модель)
Каждая экспериментальная кривая экстракции Е = q(t)/q0 - являющаяся зависимостью отношения текущего содержания масла в экстрагируемом материале от времени к его исходному содержанию, может быть представлена двумя участками - периодом постоянной скорости и
периодом падающей скорости. Таким образом, для описания кинетики экстракции можно
предложить двухзонную модель:
Е^Ш) = 1-А-(Ш) (1)
Е2(ш) = B•eхр [- х ^ К)] (2)
В этих уравнениях Е относительное (к исходному) содержание масла; N скорость в период постоянной скорости; t - время процесса; ^ - обобщенное безразмерное время, х - приведенный коэффициент экстракции.
На рисунке 1 каждая экспериментальная точка сопоставлена с расчетом по двухзонной модели (представлена треугольником). Видно как с высокой точностью точки совмещаются с треугольниками - ошибка не превышает 2...3 %. Таким образом, в каждой экспериментальной кривой в соответствии с моделью выделяются периоды постоянной и падающей скорости.
Для проверки механизма экстракции, который как принято в работах [1-8] лимитируется растворимостью масла в двуокиси углерода, сопоставим скорости экстракции в первом периоде кинетика со скоростью процесса растворения Крастворимости:
N,
растворимости
= ys Uf Pf /qo
где ys - растворимость масла, кг/кгСО2, которая определяется по уравнению [9]:
exp (40,361 -18,708/ T + 2,186840/ T2 )-(0,001 • pf )°’724
у =-
(3)
(4)
Pf
где Т - температура, К; pf - плотность двуокиси углерода в сверхкритических условиях (в зависимости от Р и Т при условиях опыта); и объемный расход двуокиси углерода; q0 - количество исходного масла в экстрагируемом материале.
На рисунке 2 представлен корреляционный график между N
растворимости
и N,
кинетика
. Степень
связи между сопоставляемыми скоростями достаточно высокая, что подтверждается механизмом растворимости в начальном периоде экстракции.
tN
1Чимкинетика, сек'1
Рис. 2 Корреляционный график между N растворимости и N кинетика
В двухзонной модели важным является граница между зонами, т.е. точка, где происходит стыковка уравнений (3) и (4). Эту границу можно отмечать по шкале относительного извлечения (координата у) и по шкале обобщенного времени N t (координата х). Между этими значениями по обеим координатам существует практически функциональная связь (рисунок 3). Для определения граничного значения можно использовать зависимость от количества растворенного экстракта (рисунок 4).
Представленные обработанные данные (рисунок 1) [1] характеризуются высоким рассеиванием, что связано с невысокой точностью экспериментов на установке объемом 10 мл. В работе [3] использовалась самая крупная установка емкостью 4 литра с исходной навеской 500 г и соответственно высотой слоя 0,45 м измельченных семян подсолнечника с исходной масличностью 45,1%. Расход двуокиси углерода во всех опытах был постоянный 5,1 кг/ч. Исследова- ния кинетики проводилось при трех режимах (Б1 Р = 30МПа, t = 50°С; Б2 Р = 34 МПа, 50°С; Б3 Р = 34 МПа,
70°С). Результаты расчетов, обработанны по представленной выше методике получения двухзонной модели, сведены в таблицу 1.
^Ммор
Рис. 3 - Связь граничных значений двухзонной модели по обеим координатам кинетической
экстракции
1,000
о о
о/
у = -С ,0262х + 0,9 ^ = 0,8697 238 О \ О
О О о^э! о э о о
1 01 5 2 0 2 5 3 0 3 5 4
*МОйУ®
Рис. 4 - Зависимость граничного значения от количества растворенного масла в первом периоде
экстракции
Таблица 1 - Результаты расчетов параметров двухзонной модели
Режим (^МПаДД) tмoD, сек EмOD, кг/кг +.між t N мOD -\tLIN - N , сек 1 -\tLIN Х
S 1(30;323) 7648 0,642 0,362 4,729E-05 8,12Ш-05 1,71722569
S2(34,323) 8487 0,459 0,537 6,332E-05 8,960E-05 1,415067852
S3(34,343) 10398 0,423 0,569 5,468E-05 6,612E-05 1,209263298
В таблице и EMOD граничные значения, отделяющие первый период с постоянной
скоростью на кривой экстракции в исходных переменных. Скорость экстракции в период постоянной скорость На границе периодов обобщенное время ^№^ш. Приведенный коэффициент
экстракции во второй период обобщенной кинетической кривой падающей скорости по экспоненциальной зависимости %.
Обобщенная кинетическая кривая с уравнениями двухзонной модели представлена на рисунке
5. Точность уравнений высокая.
E1(t•N) = 1-1,0115-(Ши№) R2 = 0.99872 (5)
Б2(^ = 0,96234^хр [- 1,38621 ^ ^ш)] R2 = 0.97484 (6)
• Е, кд/кд ОЕМООиЫ ДЕМОйЕХР
Рис. 5 - Обобщенная кинетическая кривая с уравнениями двухзонной модели по результатам опытов Б1,
82 и 83
Уравнение для расчета времени, необходимого для достижения требуемой масличности в конце процесса экстракции, имеет вид:
fK nLIN
1 EMOD __
ln 0,96234 + ln(Emod / E)'
1,0115 1,38621
или в упрощенном виде с незначительном потерей точности:
t =-
1
N
LIN
1 - E +
1 emod +
ln (EMOD / Ek )
(7)
(8)
1,38621
Таким образом, всю кинетику определяет величина постоянной скорости экстракции в первом начальном периоде NLIN. По таблице 1 видно, что наибольшая = 6,332E-05 сек-1. При этом EMOD = 0,459 кг/кг. Эта величина больше остаточного содержания масла после прессования с введением СО2 10,5 и 12,9%, это соответствует EPRESS = 0,1482 и EPRESS = 0,1803. Таким образом, эти значения меньше EMOD и процесс не лимитируется растворимостью и происходит в период падающей скорости. Для расчета времени процесса в этом случае может быть использована формула.
t = ■
1
NL
ln (Epress / Ek )
1,38621
(9)
Для определения эффективности предварительного прессования перед экстракцией сравним длительность процесса экстракции после прессования по формуле (7) с длительностью процесса экстракции от исходной масличности по формуле (6). При этом приняты конечные масличности 1 или 2%, соответствующие экстракции с применением гексана.
Таблица 2 - Результаты расчета обобщенной длительности процесса экстракции t*N без предварительного прессования и с предварительным прессованием
Исходная относительная масличность
Epress = 0,1482 Epress = 0,1803 E = 1
Ek = 0.0101 t*N = 1,93765 t*N = 2,07911 t*N = 3,29422
Ek = 0.0204 t*N = 1,43052 t*N = 1,57198 t*N = 2,78709
Основной полученный результат применения предварительного прессования это возможность примерно в 2 раза сокращение длительности процесса экстракция. При проектировании технологического процесса режимы процесса и величина конечных масличностей как на стадии прессования, так и экстракции должны дополнительно обосновываться.
Литература:
1. Salgin U., Doker O., Calimli A. Extraction of sunflower oil with supercritical CO2: Experiments and modeling // Journal of Supercritical Fluids. 2006. Vol. 38. Р. 326-331.
2. Mathematical Modeling of Sunflower Seed Extraction by Supercritical CO2 / Perrut M. and others // Ind. Eng. Chem. Res. 1997. Vol. 36. Р. 430-435.
3. Experimental and modelling of supercritical oil extraction from rapeseeds and sunflower seeds / Boutin O. and others // Chemical engineering research and design. 2011. Vol. 89.
4. Extraction of sunflower (Heliantus annuus L.) oil with supercritical CO2 and subcritical propane: Experimental and modeling / Nimet G. and others // Chemical Engineering Journal. 2011. Vol. 168. Р. 262268.
5. Roy B.C., Sasaki M., Goto M. Effect of Temperature and Pressure on the Extraction Yield of Oil from Sunflower Seed with Supercritical Carbon Dioxide // Journal of Applied Sciences. 2006. Vol. 6, №1. Р. 71-75.
6. A mathematical model of supercritical carbon dioxide extraction of sunflower oil / Esmaeili A. and others // 18th International Congress of Chemical and Process Engineering (24-28 August 2008). Prague, 2008.
7. Process optimisation in sunflower oil extraction by supercritical CO2 / Bravi M. and others // Chemical Engineering Science. 2002. Vol. 57. Р. 2753-2764.
8. Fiori L. Supercritical extraction of sunflower seed oil: Experimental data and model validation // J. of Supercritical Fluids. 2009. Vol. 50. Р. 218-224.
9. del Valle J.M., Aguilera J.M. An improved equation for predicting the solubility of vegetable oils in supercritical CO2 // Ind. Eng. Chem. Res. 1988. Vol. 27. Р. 1551.
References:
1. Salgin U., Doker O., Calimli A. Extraction of sunflower oil with supercritical CO2: Experiments and modeling // Journal of Supercritical Fluids. 2006. Vol. 38. Р. 326-331.
2. Mathematical Modeling of Sunflower Seed Extraction by Supercritical CO2 / Perrut M. and others // Ind. Eng. Chem. Res. 1997. Vol. 36. P. 430-435.
3. Experimental and modelling of supercritical oil extraction from rapeseeds and sunflower seeds / Boutin O. and others // Chemical engineering research and design. 2011. Vol. 89.
4. Extraction of sunflower (Heliantus annuus L.) oil with supercritical CO2 and subcritical propane: Experimental and modeling / Nimet G. and others // Chemical Engineering Journal. 2011. Vol. 168. P. 262268.
5. Roy B.C., Sasaki M., Goto M. Effect of Temperature and Pressure on the Extraction Yield of Oil from Sunflower Seed with Supercritical Carbon Dioxide // Journal of Applied Sciences. 2006. Vol.6, №1. P. 71-75.
6. A mathematical model of supercritical carbon dioxide extraction of sunflower oil / Esmaeili A. and others // 18th International Congress of Chemical and Process Engineering (24-28 August 2008). Prague, 2008.
7. Process optimisation in sunflower oil extraction by supercritical CO2 / Bravi M. and others // Chemical Engineering Science. 2002. Vol. 57. P. 2753-2764.
8. Fiori L. Supercritical extraction of sunflower seed oil: Experimental data and model validation // J. of Supercritical Fluids. 2009. Vol. 50. P. 218-224.
9. del Valle J.M., Aguilera J.M. An improved equation for predicting the solubility of vegetable oils in supercritical CO2 // Ind. Eng. Chem. Res. 1988. Vol. 27. P. 1551.