МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛОМАССОПЕРЕНОСА ПРИ РАСПЫЛИТЕЛЬНОЙ СУШКЕ РАСТИТЕЛЬНЫХ ЭКСТРАКТОВ
Аннотация и ключевые слова
Аннотация (русский):
Представлены результаты экспериментально-аналитических исследований процесса сушки растительных экстрактов в диспергированном состоянии. Разработана математическая модель тепломассопереноса при распылительной сушке растительных экстрактов и рассчитаны частные поля температур для промышленной организации производства сухих растительных экстрактов. Последующие теоретические и экспериментальные исследования комплекса физико-химических характеристик объектов переработки позволят дать рекомендации по рациональному проектированию конструкции сушильного аппарата и оптимальной организации технологических процессов. Предложенная математическая модель тепломассопереноса необходима для оперативного управления процессом обезвоживания. Дальнейшее развитие получают экспериментальные исследования продуктов растительного происхождения. Полученные результаты используются с целью создания сушильных комплексов для широкого перечня продуктов, с возможностью оперативного варьирования режимов работы в зависимости от вида сырья.

Ключевые слова:
моделирование технологических процессов, автоматизация техноло- гических процессов, тепломассообмен, сушка, распылительная сушка, растительные экстракты, тепломассоперенос, м
Текст
Введение В настоящее время для целей пищевой, косметической и фармацевтической отраслей промышленности активно развивается производство сухих растительных экстрактов, которые являются источниками биологически активных веществ и витаминов. Перспективно производство порошковых форм сухих растительных экстрактов. Порошки обладают рядом бесспорных преимуществ по сравнению с кусковыми и пастообразными продуктами и технологически рациональны в производственной практике. Для организации промышленного производства порошков в широком ассортименте перспективен способ, основанный на распылительной конвективной сушке специально подготовленных растительных экстрактов [1]. Для совершенствования способа распылительной сушки растительных экстрактов необходимы не только рациональные режимные параметры, но и сведения о кинетике процесса обезвоживания частиц, их гигроскопической и равновесной влажности, времени пребывания частиц в сушильной камере и допустимой температуре нагрева. При сушке распыленных частиц, движущихся в сушильной камере, их траектория имеет преимущественно вероятностный характер, энергия и масса распределяются в турбулентных потоках взаимодействующих фаз. Информации о свойствах, размере, форме, структуре частиц дисперсной фазы недостаточно для анализа и, тем более, описания их движения и процессов тепломассообмена при нестационарных краевых условиях, поэтому именно использование системы дифференциальных уравнений гидродинамики и тепломассопереноса для расчета процессов обезвоживания ограничено и невозможно без квазистационарных допущений в механизме протекания процесса и изменениях характеристик и свойств объекта обработки. Для переноса результатов экспериментов в опытных установках необходимо соблюдение условий масштабного перехода. Приближенное моделирование сушки следует осуществлять на основе комплексного применения расчетных и экспериментальных методов исследования. Температурные режимы при сушке термолабильных продуктов растительного происхождения чрезвычайно важны. Необходимость расчета значений температуры в течение процесса сушки и в каждой пространственной точке объекта обезвоживания очевидна из-за сложности экспериментального определения полей температур в хаотично движущихся тонкодисперсных частицах. С учетом высокой интенсивности процесса распылительной сушки реализация математической модели с целью определения температурных полей не вызывает сомнений. При моделировании сушки взвешенных капель (частиц) растительных экстрактов малого диаметра их можно представить как «тонкий слой», который предполагает малые значения теплового и массообменного чисел Био (, ) при однородности полей температур и влагосодержания в материале. Если принять частицу шарообразной, то энергоподвод можно считать равномерным по всей поверхности частицы, поэтому за абсциссу, от которой зависит температурная функция, принят радиус частицы. При этом можно не использовать полярные координаты, т. к. эволюция значений температуры не зависит от угла и направления энергоподвода. В отличие от моделей процесса сушки [2], где пренебрегают усадкой, при распылительной сушке диспергированных частиц усадка неизбежна. При численном решении уравнения переноса тепла усадка учитывалась дискретно (ступенчато) с шагом, обеспечивающим необходимую точность, что и предполагает метод конечных разностей, которыми заменяются дифференциалы. Конечные условия (распределение температур по радиусу частицы) предыдущего участка являются начальными для последующего. За вторую координату при составлении разностной сетки берется время сушки или связанная с ней влажность продукта (экспериментальные данные), также разбитые на характерные участки. Окончательная совместная компоновка эволюций температур на каждом участке позволяет получить общее поле температур в процессе сушки и динамику изменения влагообменных коэффициентов с заданной нами точностью. Ввиду того, что точную зависимость размера частицы от времени или влажности определить невозможно, принимаем линейное изменение диаметра от начального до конечного, который определен экспериментально на основе дисперсного анализа. Разработка методики расчета проводилась с учетом допущений и положений, обоснованных в [2]. Нахождение нестационарных полей температур связано с решением системы дифференциальных уравнений влаго- и теплопереноса. Целесообразной является реализация математических численных методов, с помощью которых можно с высокой точностью рассчитать параметрические поля, т. к. аналитические решения достаточно громоздки и некорректны. Рассчитаны частные поля температур при конвективной распылительной сушке растительных экстрактов. Уравнение переноса тепла при одномерной задаче в случае объемного энергоподвода имеет вид [3]: , где ρ(x, Т, w), сρ(x, Т, w), λ(w, x, Т) – функции плотности, теплоемкости и температуропроводности; W(w, x, Т) – функция распределения объемной плотности поглощенной энергии излучения по глубине оптически тонкого слоя; r(w, x, t) – теплота испарения, определенная при термодинамическом анализе механизма внутреннего массопереноса; x – координата; ε – коэффициент фазовых превращений. Здесь r – количество тепла для испарения 1 кг влаги с учетом энтропийной составляющей: В нашем случае внутренний источник тепла (объемный энергоподвод при инфракрасной и сверхчастотной сушке) отсутствует: . Кроме того, можно подставить вместо дифференциальное изменение средней по слою влажности , а также с достаточной точностью принять коэффициент фазовых превращений ε = 1, т. к. перенос влаги в процессе высокоинтенсивной сушки происходит в основном в виде пара при изотропности (равномерности w по слою) структуры [2]. Поскольку структура при cyшке изотропна, теплофизические характеристики практически не зависят от x, и если вынести за знак дифференциала, разделить уравнение на получим , где – коэффициент температуропроводности. Имея функции и, соответственно, , после преобразования (далее знак среднего и варьируемые параметры будем опускать) получим . (1) Начальные условия для первой зоны зададим, учитывая равномерное распределение температуры в начальный момент времени, соответствующий , в виде: при , т. е. . При этом конечные условия (распределение температур по радиусу частицы) предыдущего участка являются начальными для последующего. Так как фактически отражателем является поверхность продукта на соседних частицах, имеющего такую же , то лучистым теплообменом между поверхностями частиц можно пренебречь. Теплообмен на границе задается граничными условиями 2-го рода. Ввиду того, что основной целью облучения является создание оптимального распределения плотности теплового потока на поверхности продукта при сушке, имеют место граничные условия , где α – коэффициент теплоотдачи, Вт/(м2 · К). Обозначим: x – координата глубины частицы, м; – температура окружающей среды, К; – температура на поверхности частицы, К. Решением уравнения (1) при краевых условиях является функция , подставив в которую , получим искомую функцию . Расчеты проводились по зонам сушки. Чтобы оба аргумента, по которым ведется дифференцирование, возрастали в процессе обезвоживания, заменим влажность w концентрацией сухих веществ c. Для численного решения дифференциального уравнения параболического типа в частных производных при разнородных граничных условиях удобен метод конечных разностей. При решении одномерного уравнения параболического типа отдают предпочтение неявным схемам, которые, безусловно, устойчивы при любом соотношении шагов по x и c. На рисунке в качестве примера приведены результаты расчетов полей температур для распылительной сушки пектинового экстракта. Температурное поле при Wн = 0,15 кг/кг; T = 423 К: с, х – номера шагов (Δс, Δх) по осям с (концентрация) и х (глубина частицы). Размер шага – отношение диапазона варьирования (по с 0,15–0,9 кг/кг, по х: при с = 0,15–0,3, х = 0,000015 м; с = 0,3–0,45, х = 0,000012 м; с = 0,45–0,6, х = 0,000009 м; с = 0,6–0,75, х = 0,000006 м; с = 0,75–0,9, х = 0,000003 м) Анализ температурных полей свидетельствует о наличии малых температурных градиентов, что предопределяет равномерность прогрева и обезвоживания частиц, «мягкие» режимы сушки. Это подтверждает анализ расчетных данных, т. к. температура продукта в течение конвективной сушки не превышает T = 353 К. Отмечены периодические скачки температуры, которые соответствуют пикам скорости сушки. На участках с неизменной температурой удаляется свободная влага [3, 4]. Рост температуры наблюдается при расширении структурной влаги в замкнутых объемах (закрытых жидкостными менисками капиллярах и клетках) – парниковый эффект при снижении скорости сушки. Заключение Таким образом, нами проведено моделирование тепломассопереноса при распылительной сушке растительных экстрактов и рассчитаны частные поля температур для последующей промышленной организации производства сухих растительных экстрактов. Дальнейшие теоретические и экспериментальные исследования комплекса физико-химических характеристик объектов исследования позволят выработать рекомендации по рациональному проектированию конструкции сушильного аппарата, а также разработать физико-математические модели тепломассопереноса для оперативного управления процессом обезвоживания. Дальнейшее развитие получают экспериментальные исследования продуктов растительного происхождения. Результаты могут быть использованы для создания сушильных комплексов для широкого перечня продуктов, с возможностью оперативного варьирования режимов работы в зависимости от вида сырья.
Список литературы

1. Максименко Ю. А. Термодинамика внутреннего массопереноса при взаимодействии плодоовощных продуктов с водой // Вестн. Астрахан. гос. техн. ун-та. – 2012. – № 1 (53). – С. 41–45.

2. Алексанян И. Ю., Буйнов А. А. Высокоинтенсивная сушка пищевых продуктов. Пеносушка. Теория. Практика. Моделирование. – Астрахань: Изд-во АГТУ, 2004. – 380 с.

3. Лыков А. В. Теория сушки. – М.: Энергия, 1968. – 471 с.

4. Гинзбург А. С. Основы теории и техники сушки пищевых продуктов. – М.: Пищ. пром-сть, 1975. – 527 с. REFERENCES


Войти или Создать
* Забыли пароль?