Моделирование процесса сушки структуры потоков теплоносителей в водонагревательной установке Текст научной статьи по специальности «Химические технологии»

МОДЕЛИРОВАНИЕ ПРОЦЕССА СУШКИ СТРУКТУРЫ ПОТОКОВ ТЕПЛОНОСИТЕЛЕЙ

В ВОДОНАГРЕВАТЕЛЬНОЙ УСТАНОВКЕ

Султанова Шахноза Абдувахитовна

соискатель, Ташкентский государственный технический университет, 100095, Узбекистан, г. Ташкент, ул. Университетская, дом № 2

E-mail: sh. sultanova@yahoo. com

MODELING OF THE DRYING PROCESS OF THE STRUCTURE OF HEAT FLOWS IN THE WATER HEATER INSTALLATION

Shaxnoza Sultanova

aspirant, Tashkent State Technical University, 100095, Uzbekistan, Tashkent, University st, 2

АННОТАЦИЯ

В статье изложены результаты экспериментальных исследований проведённые как на лаборатории и теории. В результате проведения процесса сушки в установке были подобраны соответственные температуры, а также направление хода теплоносителя Осуществлено математическое моделирование температурного поля на основе трехфазной структуры потоков. В широком диапазоне варьирования конструктивных и режимных параметров, идентифицированы их номинальные значения.

ABSTRACT

The article describes the results of experimental research conducted both at the laboratory and theory. As a result of carrying out the drying process, the corresponding temperatures were selected in the installation, as well as the direction of the coolant flow. A mathematical simulation of the temperature field based on the three-phase flow structure was carried out. In a wide range of variation of design and mode parameters, their nominal values are identified.

Ключевые слова: продукт, установка, температура, теплоноситель, сушка, процесс, модель, агент, тепловой поток.

Keywords: product, installation, temperature, coolant, drying, process, model, agent, heat flow.

Процессы обезвоживания (сушки), то - есть процессы удаления влаги из продукта, имеют большое значение в различных отраслях пищевой промышленности и неразрывно связаны друг с другом. По мнению доктора технических наук Слуцкой Т.И., доктора технических наук Мезенова О.Я., доктора технических наук Ким Э.Н. под сушкой понимается процесс удаления из материалов любой жидкости, в результате чего увеличивается относительное содержание сухой части. При этом, в процессе сушки пищевых продуктов удаляется главным образом вода, поэтому под сушкой понимают процесс обезвоживания материалов [2, с. 36], [4, с. 203].

Основное отличие сушки от механических способов отделения влаги или ее химического поглощения - повышение температуры и переход растворителя из жидкой фазы в газообразную. Эффективное применение процесса сушки невозможно без использования современных высокоинтенсивных сушильных установок [3, с.36-39].

Кинетика сушки влажного материала включает в себя определение оптимальных параметров сушиль-

ного агента (температура, расход, влажность и давление), а также конструкцию и габариты сушильной установки. [1, с .102], [5, с.68].

В установившимся режиме уравнения теплового баланса полученных на основе теоритических исследований по секциям и фазам можно записать следующим образом:

по фазе теплоносителя:

!170 - ¿71 - - = 0 - для 1 - ого поддона ¿71 - 1Т2 - ^(Г2 - ¿2) = 0 - для 2 - ого поддона ^

1ТП - 1ТП - - £п) = 0 - для п - ого поддона по паровой фазе:

С£0 - а1 + - £^=0 - для 1-ого поддона С£0 - С£:2 + - £2)=0 - для 2-ого поддона (2) С£0 - + - £„)=0 - для п -ого поддона,

Здесь: ^-температура на поверхности соответствующего змеевика (принимаем равной средней температуре воды); ^-температура воздуха в соответствующей секции; 70-входная температура воды; t0-

Библиографическое описание: Сафаров Ж.Э., Дадаев Г.Т. Исследование процесса нагревания и охлаждения при использовании различных материалов для аккумулияции тепловой энергии // Universum: Технические науки: электрон. научн. журн. 2017. № 11(44) . URL: http://7universum.com/ru/tech/archive/item/5234

входная температура воздуха; ¿-произведение расхода воды на теплоемкость (некоторый аналог потока энтропии); О-аналогично для воздуха; К-отношение теплообменного потока к разности температур фаз.

Для перехода к безразмерным переменным введем обозначения отношений интенсивностей поток в

ков энтрапии к Ь: к = — ,9=~ •

Тогда уравнения для 1-го поддона можно выразить так:

Г, =

Тр + М1 1+к д10+кТ1 д + к

(3)

Такое представление удобно для решения уравнений теплового баланса по отношению к неизвестным значениям температур методом итерации. Оно особенно удобно при рассмотрении противоточных и сложных организаций потоков и, в дальнейшем, с учетом нелинейностией.

Расчет производится последовательно по поддонам снизу вверх. Вычисления показывают, что достаточно 10-кратных измерений для достижения требуемой точности.

Рисунок 1. Схема прямоточного движения тепловых потоков

ч =

Прежде чем начать исследование точности идентификации модели, необходимо получить квазистационарные оценки g и к для исследуемого материала. На рис. 1. представлены экспериментальные значения профиля температуры по высоте, где в качестве исследуемого материала были использованы лекарственные растения.

Идентификация модели на экспериментальных данных показывает, что:

1. Модель верно отражает понижение температуры воды по высоте и повышение температуры воздуха приблизительно по экспоненциальному закону;

2. Перекос расчётных и экспериментальных кривых даже для последних стадий процесса свидетельствуют о необходимости включения в модель параметра оттока тепла из секции за счет испарения.

Представленное на рис. 2. сравнение приближений расчетных и экспериментальных профилей свидетельствует о том, что допущение о равенстве температур газовой фазы и материала не позволяет получить требуемый уровень адекватности.

Рисунок 2. Расчетные и экспериментальные профили температур с учетом потока на испарение

На этом этапе также проверяется противоточная схема организации потоков теплоносителя и воздуха (рис. 3 и рис. 4)

Уравнение баланса для первого поддона: по жидкой фазе

ЬТ2-ЬТ1-К(Т1-11)-2=0,

по газовой фазе

О0-Ои+К(Т1-Г1)=0,

(4)

(5)

Рисунок 4. Противоточная схема организации движения тепловых потоков

в итерационной форме

Т, =

и =

Т2-г+Ы1 1+к '

д1р+кТ1 д + к

(6) (7)

Рисунок 3. Противоточная схема организации движения тепловых потоков

На рис. 5. представлены профили противоточной схемы организации потоков. Результаты красноречиво свидетельствуют о недопустимых неравномер-ностях сушки по поддонам.

Реализация процедуры идентификации показала, что оценки параметров модели лежат в диапазоне к около ~0,36, g около ~0,56, а г принимает минимальное значение.

Следует отметить, что решение задачи идентификации сильно зависит от начальных приближений параметров модели. Для проверки результатов также получено решение системы линейных уравнений методом обратной матрицы.

Матрица коэффициентов системы уравнений (45) представляет собой четырёх- диагональную матрицу:

А,+и=к Аи+2=1 А+3, i+1i=g

1=1,4 . . . 2П i=1,3 . . . 2П-3 i=1,3 . . . 2П-3

Aii=-1-k Aii=-g-k

Аи+1=к

1=Х,3 . . . 2П-1 !=2,4 . . . 2П i=1,3 . . . 2П-1

Время расчета по методу обратной матрицы приблизительно в два раза меньше, чем время расчета по итерационному методу.

Рисунок 5. Профили расчетных и экспериментальных значений температур при противоточной схеме

организации потоков

Сама матрица выглядит следующим образом:

Т1 T2 t2 Tз tз Тп tn

-1-к k 1 Г

k -g-k -ЯЪ

-1-k k 1 г

g k -Я-к 0

-1-k X

Кроме того, как в предыдущих случаях, применяем раздельный учет температур воды, материала и воздуха. Соответствующая схема организации потоков представлена на рис. 6.

Уравнения для стационарного режима для 3-х фазной модели:

для 1 -го поддона

^„о - - ^пт(Г„1 - Гт1) + КПД(Т'П1 - 7^1) = 0 - Т'т1) - ^тя(Гт1 - Тд1) = 0

- - (^П1 - + (^Ш1 - = 0

(8)

Рисунок 6. Трехфазная модель движения тепловых потоков

Введем обозначения:

кпт _ т. . кпд _ , _ ктд _

ктд; - — г (9)

Приведем рассматриваемые соотношения к итерационному виду:

—

Тп0 + кптГт1 + кпдТд1

91

_ кптГп1 + ктдТд1-г

кпт+ктд _ 9Тдо + кПдТП1 + ктдТт1

З+кпд+ктд

(10)

На рис. 7. представлены результаты идентификации трехфазных моделей.

Получены оценки параметров модели (при условии Тю=95 0С; Tgo=12 0С): я=0,5752; ^=0,0261; kmg=0,0174; kng=0,0936; г=0. Последнее значение отражает конечную стадию процесса.

Для оценки точности идентификации проведено статистическое исследование с имитацией погрешности измерения температуры по нормальному закону с дисперсией 2 0С. Средние ошибки оценок стремятся к нулю, что показывает состоятельность оценок. Средние же абсолютные ошибки: (%=7%; дкпт=13%; &1т=15%>; дк.щ,=13%о.

-о- средняя температура в трубопроводе поддона; —* - средняя температура материала; - средняя температура газовой фазы.

Рисунок 7. Идентификация трехфазной модели

Результаты показывают определленую целесо-бразность отдельного (автономного) расстмотрения теплообменных аспектов процесса с целью

выработки конструктивных рекомендаций относительно структуры организации контактирующих потоков теплоносителя.

Список литературы:

1. Лебедев П. Д. Расчет и проектирование сушильных установок. М.: Госэнергоиздат, 1963. -320 с.

2. Султанова Ш.А. Разработка конвективной сушильной установки контейнерного типа предназначенной для лекарственных растений. Международный НТЖ «Химическая технология. Контроль и управление». Ташкент, 2017 №1 (73). С. 36-40.

3. Сафаров Ж.Э., Султанова Ш.А. Изучение теплофизических характеристик растительного материала для качественной сушки. Продовольча iндустрiя АПК науково_практичний журнал. 2017 №1-2. С. 36-39.

4. Технология продуктов из гидробионтов/ С.А. Артюхова и др.; под ред.Т.М. Сафроновой и В.И. Шендерюка.-М: Колос 2001. -496 с.

5. Чернобыльский И.И., Тананайко Ю.М. Сушильные установки в химической промышленности. Киев: Техника, 1969. -279 с.

1 + кпт+кпд

Т.

V