NUMERICAL STUDY OF JUICE CONCENTRATION IN THE CONVECTION APPARATUS WITH SWIRLING HEAT-TRANSFER AGENT FLOW
Authors:
1.
Kazan National Research Technological University
2.
Kazan National Research Technological University
Type:
Article
Pages:
from 119
to 124
Status:
Published
Received:
30.06.2016
Accepted:
30.06.2016
Published:
30.06.2016
Subject area:
UDK 66.048.5:664.857.3
UDK 66 Химическая технология. Химическая промышленность. Родственные отрасли
UDK 60 Прикладные науки. Общие вопросы
UDK 66 Химическая технология. Химическая промышленность. Родственные отрасли
UDK 60 Прикладные науки. Общие вопросы
Language:
Russian
Keywords:
Concentrate, juice, vortex camera, évaporation, drop, trajectory, air tangent vortex generator
Abstract (English):
The paper presents a highly potent vortex apparatus for fruit and vegetable juice concentration in the “gentle” regime ensuring the preservation of biologically active substances and high flavor quality of the product. Calculation of the convection apparatus with air tangent vortex generator with the concentrate’s known properties is based on the material and heat balances, together with condition that necessary evaporation time, induced by the process kinetics, should be provided by the chamber and vortex generator constructions, determined by hydro- aerodynamic laws and heat- and mass transfer distinctive features in the apparatus. The main requirement for the concentration termination is to achieve necessary mass content of soluble solids in the product. The designed mathematical model allows us to describe the behavior of an evaporating drop in a swirling gas flow over the wide range of calculated parameters and operating regimes for the vortex camera. Numerical calculation shows that average evaporation time for various drop diameters decreases 8.8 times with circulation ratio increasing from 0 to 5 in the concentrate liquid. The fact of the relative apparatus height reduction under the rise of the average consumed gas velocity in a vortex generator with decreasing blades’ inclination of the vortex generator and a liquid-gas ratio in various cases of the initial liquid drop diameter has been established. The numerical calculation indicates that the narrowing of annular rotating drop layer of liquid on increasing the circulation ratio for “smooth” entering of a drop onto its trajectory needs observance of the principle of equality between the initial drop tangential velocity and gas tangential velocity. According to the calculation, the gas temperature within 105 ± 25°C with liquid circulation ratio over two does not influence significantly the drop trajectory.
The paper presents a highly potent vortex apparatus for fruit and vegetable juice concentration in the “gentle” regime ensuring the preservation of biologically active substances and high flavor quality of the product. Calculation of the convection apparatus with air tangent vortex generator with the concentrate’s known properties is based on the material and heat balances, together with condition that necessary evaporation time, induced by the process kinetics, should be provided by the chamber and vortex generator constructions, determined by hydro- aerodynamic laws and heat- and mass transfer distinctive features in the apparatus. The main requirement for the concentration termination is to achieve necessary mass content of soluble solids in the product. The designed mathematical model allows us to describe the behavior of an evaporating drop in a swirling gas flow over the wide range of calculated parameters and operating regimes for the vortex camera. Numerical calculation shows that average evaporation time for various drop diameters decreases 8.8 times with circulation ratio increasing from 0 to 5 in the concentrate liquid. The fact of the relative apparatus height reduction under the rise of the average consumed gas velocity in a vortex generator with decreasing blades’ inclination of the vortex generator and a liquid-gas ratio in various cases of the initial liquid drop diameter has been established. The numerical calculation indicates that the narrowing of annular rotating drop layer of liquid on increasing the circulation ratio for “smooth” entering of a drop onto its trajectory needs observance of the principle of equality between the initial drop tangential velocity and gas tangential velocity. According to the calculation, the gas temperature within 105 ± 25°C with liquid circulation ratio over two does not influence significantly the drop trajectory.
Keywords:
Concentrate, juice, vortex camera, évaporation, drop, trajectory, air tangent vortex generator
Concentrate, juice, vortex camera, évaporation, drop, trajectory, air tangent vortex generator
Введение В настоящее время концентраты фруктовых и овощных соков, будучи основой производства как восстановленных соков, так и многих продуктов питания и безалкогольных напитков, пользуются широким спросом. Среди существующих способов получения концентрированных соков самым распространенным является выпаривание, которое всегда приводит к негативным сложным физико-химическим изменениям исходного сока, что является причиной непрерывного поиска принципиально новых технологий и оборудования, обеспечивающих высокие вкусовые качества и биологическую активность конечного продукта. Одним из видов такого перспективного и высокоэффективного технологического оборудования для концентрирования жидкостей рассматривается конвекционный аппарат с тангенциально-лопаточным завихрителем потока газа. Обзор классических и современных работ, посвященных исследованию динамического взаимодействия частиц дисперсной фазы с несущей средой в двухфазных системах, представленные в работах [1-4], показал, что использование закрученных потоков в технических устройствах ограничено отсутствием надежных экспериментальных и расчетных данных о характеристиках течения и тепломассообмене в таких потоках. Аппарат для концентрирования фруктовых и овощных соков в закрученном газовом потоке работает следующим образом (рис. 1). Греющий газ, представляющий собой паровоздушную смесь, предварительно нагретый в калорифере, подается через осевой патрубок газораспределительного устройства в вихревую камеру. Проходя через тангенциально-лопаточный завихритель, газ закручивается, и в рабочей зоне камеры образуется вращающийся газовый поток. Подвод сока осуществляется через патрубки в кольцевой трубчатый ороситель, расположенный вблизи крышки внутри рабочей камеры. Жидкий продукт дробится высокоскоростным потоком газа на капли, образуя вращающийся газожидкостный слой. Такое взаимодействие обеспечивает высокую тепломассообменную эффективность процесса испарения воды с поверхности капель, что способствует интенсификации процесса концентрирования. Насыщенный влагой газ удаляется через центральный выходной патрубок. В выходном патрубке капли жидкости отделяются от газового потока, жидкость через щелевой отсекатель попадает в сливной стакан, откуда самотеком выводится из камеры. Чтобы предотвратить разложение витаминов и других полезных веществ, а также пригорание продукта и карамелизацию сахаров, стенка выходного патрубка постоянно ох-лаждается. При работе вихревой камеры без рециркуляции жидкости вся жидкая фаза с начальной концентрацией сухих веществ однократно взаимодействует с газом и выводится из камеры с требуемой конечной концентрацией. При работе камеры с рециркуляцией жидкости прошедшая через рабочую зону камеры жидкость попадает в циркуляционный контур и, смешиваясь со свежей исходной жидкостью, вновь поступает в аппарат. Согласно рециркуляционной схеме содержание сухих веществ в поступающей жидкости выше, чем без рециркуляции. При помощи рециркуляции можно интенсифицировать процесс концентрирования в камере и увеличить конечную концентрацию сухих веществ. Рис. 1. Конвекционный аппарат для концентрирования соков: С - сок; Г - газ; ОА - охлаждающий агент; К - концентрированный сок; 1 - корпус; 2 - крышка; 3 - днище; 4 - тангенциальный завихритель; 5 - лопасти; 6 - газораспределительное устройство; 7 - патрубок ввода газа; 8 - патрубок ввода жидкости; 9 - трубчатый ороситель; 10 - патрубок выхода газа; 11 - щелевой отсекатель; 12 - патрубок выхода концентрата; 13 - охлаждающая рубашка; 14, 15 - технологические патрубки При проектировании вихревого оборудования необходим расчет параметров аппарата и эффективности, основанные на реальной гидроаэродинамической обстановке в рабочей зоне и особенностей тепло- и массообмена. В работе [5] приведена математическая модель движения испаряющейся капли жидкости в закрученном потоке теплоносителя в конвекционном аппарате. Целью работы является численное исследование процесса с использованием разработанной математической модели и выбор рабочих конструктивных и режимных параметров аппарата. Объекты и методы исследований Моделирование проводилось при следующих допущениях: - активные силы, в решающей степени определяющие движение капли в несущем газовом потоке, - сила аэродинамического сопротивления и сила тяжести; - не учитываются циркуляционные течения внутри капель, их взаимодействие между собой, дробление и коалесценция; - тепломассоперенос осуществляется на поверхности капли; - теплофизические параметры паровоздушной смеси в пределах пограничного слоя рассчитываются при среднеарифметической температуре и концентрации водяного пара; - принимается полное смешение газа по толщине капельного слоя; - все сухие вещества сока рассматриваются как один компонент с заданными теплофизическими свойствами; - жидкость в начальный момент времени находится при температуре фазового перехода. Движение капли жидкости в неподвижной цилиндрической системе координат в вихревой камере описывалось системой уравнений (1) где , , - цилиндрические координаты; , , - компоненты вектора абсолютной скорости капли; , , - компоненты вектора скорости газа; - плотность газа и жидкости; - коэффициент лобового сопротивления капли; - число Рейнольдса капли; - время; - диаметр капли; - скорость капли относительно газа; - коэффициент массоотдачи; - концентрация пара воды в паровоздушной смеси на поверхности капли и вне пограничного слоя соответственно. Численное решение системы уравнений (1) получено четырехшаговым методом Рунге-Кутта с начальными условиями = 0, = , = 0, = , = , = , = 0, = , (2) где - начальный радиус ввода жидкости (радиус оросителя); - расстояние от крышки камеры (высота ввода жидкости); - скорость истечения жидкости из оросителя; - начальный диаметр капли. Характерными элементами алгоритма сформированного расчета являются: 1) определение температуры на поверхности испаряющихся капель методом последовательных приближений из условия термодинамического равновесия на межфазной поверхности (3) где - удельная теплота парообразования на поверхности капли; - температура газа на удалении от капли; 2) влагосодержание и температура паровоздушной среды корректируются в ходе расчета по выражениям (4) (5) где , - влагосодержание и температура газа на входе; - отношение массовых расходов жидкости и газа на входе в камеру; - отношение конечной и начальной плотности капли; - отношение конечного и начального диаметра капель; - средняя изобарная теплоемкость паровоздушной смеси, сухого воздуха и водяного пара соответственно. Расчеты проводились для газожидкостной системы паровоздушная смесь-яблочный сок (Г-Ж). Необходимым условием окончания процесса концентрирования является достижение требуемого конечного содержания сухих веществ в продукте, которому соответствует конечный диаметр капли . Результаты и их обсуждение В ходе расчетов определялось поведение испаряющихся капель по времени в зависимости от диаметра капель, кратности циркуляции жидкости, начальной скорости дисперсной фазы и температуры воздушной среды. Расчет траекторий капель в конвекционном аппарате радиусом , работающем по принципу принудительной циркуляции концентрируемой жидкости, показал, что достаточная высота цилиндрической части камеры и требуемый диаметр выходного патрубка определяются кратностью циркуляции жидкой фазы . Капли с увеличением кратности циркуляции имеют бо́льшую массовую долю сухих веществ на входе и для достижения нужной степени концентрирования проходят более короткий путь, сокращая время пребывания продукта в камере. Результаты расчета времени испарения капель до конечного диаметра в камере при разных значениях кратности циркуляции жидкости представлены на рис. 2. При отсутствии циркуляционного контура ( 0) время пребывания капель диаметром 0,25 м в аппарате с радиусом 0,35 м составляет 3,56 с, а с циркуляцией от 0 до 5 время снижается в 7,91 раза; капель диаметром 0,30 м в аппарате с радиусом 0,46 м - в 9,07 раза и капель диаметром 0,37 м в аппарате с радиусом 0,65 м - в 9,37 раза. Рис. 2. Время достижения каплей заданного конечного содержания сухих веществ в зависимости от кратности циркуляции. Расчетные параметры: 35 ; 12 шт.; 13 м/с; 0,2; , мм: 1 - 0,25; 2 - 0,30; 3 - 0,37 Во время расчета движения испаряющихся капель были изучены следующие параметры, варьирование которыми оказывает непосредственное влияние на траекторию капель: среднерасходная скорость газа в живом сечении завихрителя изменялась от 10 до 20 м/с (рис. 3), угол наклона профилированных лопастей тангенциального завихрителя - от 25 до 45 градусов (рис. 4) и соотношение массовых расходов жидкости и газа - от 0,1 до 0,5 (рис. 5) для начальных диаметров капель = 0,25; 0,30 и 0,37 мм и при 3 при степени концентрирования 70 % значения = 0,225, 0,27 и 0,333 мм соответственно. Рис. 3. Зависимость высоты камеры от скорости газа в живом сечении завихрителя: 3; 35 ; 12 шт.; 0,2; , мм: 1 - 0,25; 2 - 0,30; 3 - 0,37 Рис. 4. Зависимость высоты камеры от угла наклона лопастей тангенциального завихрителя: 3; 12 шт.; 0,2; 15 м/с; , мм: 1 - 0,25; 2 - 0,30; 3 - 0,37 Рис. 5. Зависимость высоты камеры от отношения массовых расходов жидкой и газовой фаз: 3; 35 ; 12 шт.; 15 м/с; , мм: 1 - 0,25; 2 - 0,30; 3 - 0,37 Анализ рассчитанных траекторий капель дисперсной жидкой фазы на входе в аппараты с тангенциальным завихрителем разного диаметра показал, что не соблюдается подобие траекторий капель при равных значениях режимных параметров. В общем это связано с различной степенью влияния следующих определяющих факторов: начальных условий движения, источника и положения потока дисперсной фазы, направления истечения дисперсной фазы. Расчеты при разных кратностях циркуляционного контура показали, что величина начальной тангенциальной компоненты скорости капель определяется толщиной кольцевого слоя , который сужается при увеличении и требует значение , близкое к тангенциальной составляющей скорости газового потока на входе в завихритель. Высокий уровень тангенциальной скорости практически предопределяет мгновенный выход капель на равновесную траекторию, что крайне необходимо при высоких значениях кратности циркуляции (рис. 6). Рис. 6. Зависимость начальной тангенциальной скорости капли от ширины относительного кольцевого слоя: 0,48 м; 35 ; 12 шт.; 13 м/с; 0,2; 0,25 мм На рис. 7 представлены результаты расчета относительной высоты вихревой камеры при изменении температуры греющего газа для различных случаев циркуляции жидкой фазы. Для капель диаметром 0,25 мм с увеличением температуры газа в пределах от 80 до 130 °С величина уменьшается. Вместе с тем при 0 и 1 наблюдается более значительное влияние температуры среды на относительную высоту камеры, чем при кратности циркуляции 2, 3, 4 и 5. Рис. 7. Зависимость геометрии камеры от температуры газа: 0,35 м; 35 ; 12 шт.; 15 м/с; 0,2; 0,25 мм; : ● - 0; ▼ - 1; × - 2; ▲ - 3; ♦ - 4; ■ - 5 В целом по результатам проведенного численного эксперимента движения капель сока при испарении в конвекционном аппарате с тангенциально-лопаточным завихрителем можно сделать следующие выводы: 1) с увеличением кратности циркуляции по концентрируемой жидкости от 0 до 5 время испарения капель разного диаметра снижается в среднем в 8,8 раза; 2) уменьшение относительной высоты рабочей зоны аппарата наблюдается при росте среднерасходной скорости газа в живом сечении завихрителя, при уменьшении угла наклона профилированных лопастей завихрителя и соотношения массовых расходов жидкости и газа для разных начальных диаметров капель жидкости; 3) сужение кольцевого вращающегося капельного слоя жидкости при увеличении кратности циркуляции требует обеспечения начальной тангенциальной составляющей скорости капель, близкой к тангенциальной скорости газового потока, для более «плавного» выхода на равновесную траекторию; 4) температура паровоздушной среды в пределах (105 25) °С при кратности циркуляции 2 не оказывает существенного влияния на геометрию камеры.
1. Busroyd, R. Techenie gaza so vzveshennymi chasticami / R. Busroyd; perevod s angl. pod red. Z.R. Gorbisa. - M.: Mir, 1975. - 380 s.
2. Chepkasov, V.M. Vliyanie struktury gazovogo potoka na dvizhenie dispersnoy fazy v vihrevom separatore / V.M. Chepkasov, A.A. Ovchinnikov, N.A. Nikolaev // Izv. vuzov. Himiya i him. tehnol. - 1981. - T. 24. - № 5. - S. 639-642.
3. Sikalo, S. Hydrodynamics and heat transfer investigation of air-water dispersed flow / S. Sikalo, N. Delalic, E. N. Ganic // Intern. J. Exp. Thermal Fluid Sci. - 2002. - V. 25. - № 7. - P. 511-521.
4. Laptev, S.A. Zakonomernosti povedeniya dvuhfaznogo gazozhidkostnogo potoka v vihrevyh apparatah / S.A. Laptev // Vestnik tehnologicheskogo universiteta. - 2015. - T. 18. - № 23. - S. 116-118.
5. Har'kov, V.V. Modelirovanie teplo- i massoobmena pri koncentrirovanii sokov v vihrevoy kamere / V.V. Har'kov // Nauchno-tehnicheskiy vestnik Povolzh'ya. - 2016. - № 1. - S. 37-44.
