Publication text
(PDF):
Read
Download
Введение В настоящее время концентраты фруктовых и овощных соков, будучи основой производства как восстановленных соков, так и многих продуктов питания и безалкогольных напитков, пользуются широким спросом. Среди существующих способов получения концентрированных соков самым распространенным является выпаривание, которое всегда приводит к негативным сложным физико-химическим изменениям исходного сока, что является причиной непрерывного поиска принципиально новых технологий и оборудования, обеспечивающих высокие вкусовые качества и биологическую активность конечного продукта. Одним из видов такого перспективного и высокоэффективного технологического оборудования для концентрирования жидкостей рассматривается конвекционный аппарат с тангенциально-лопаточным завихрителем потока газа. Обзор классических и современных работ, посвященных исследованию динамического взаимодействия частиц дисперсной фазы с несущей средой в двухфазных системах, представленные в работах [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 не оказывает существенного влияния на геометрию камеры.