CC BY
22
УДК 621.184.64 ББК 31.31 О-92
Схаляхов Анзаур Адамович, кандидат технических наук, доцент, профессор кафедры технологий, машин и оборудования пищевых производств, декан технологического факультета Майкопского государственного технологического университета, 385000, Республика Адыгея, г. Майкоп, ул. Первомайская, 191, т.: (8772) 570412;
Блягоз Хазрет Рамазанович, доктор технических наук, профессор, профессор кафедры технологий, машин и оборудования пищевых производств, ректор Майкопского государственного технологического университета, 385000, Республика Адыгея, г. Майкоп, ул. Первомайская, 191, т.: (8772) 570011;
Кошевой Евгений Пантелеевич, доктор технических наук, профессор, заведующий кафедрой машин и аппаратов пищевых производств Кубанского государственного технологического университета, 350072, г. Краснодар, ул. Московская, 2, т.: (861) 2752279;
Косачев Вячеслав Степанович, доктор технических наук, профессор, профессор кафедры машин и аппаратов пищевых производств Кубанского государственного технологического университета, 350072, г. Краснодар, ул. Московская, 2, т.: (861) 2752279.
ОПИСАНИЕ ТЕПЛООБМЕНА В ПОЛОВОЛОКОННЫХ МЕМБРАНАХ ДЛЯ ВНУТРЕННЕЙ
СОПРЯЖЕННОЙ ЗАДАЧИ
(рецензирована)
Получение в результате численного моделирования обобщенной кинетической зависимости теплообмена позволяет обобщить результаты численного решения сопряженной задачи теплообмена от внутреннего потока в половолоконных мембранах и может быть применена в практических расчетах.
Ключевые слова: половолоконные мембраны, сопряженная задача, теплообмен,
моделирование, кинетическая зависимость.
Skhalyakhov Anzaur Adamovich, Cand. of Technology, associate professor, professor of the chair of technologies, machines and equipment of food products, dean of technological faculty of Maikop State Technological University, 191,Pervomayskaya St. Maikop, AdighRepublic, zip: 385000, tel.: (8772)570412;
Blyagoz Khazret Ramazanovich, Doctor of Technology, professor, professor of the chair of technologies, machines and equipment of food products, dean of technological faculty of Maikop State Technological University, 191,Pervomayskaya St. Maikop, Adigh Republic, zip: 385000, tel.: (8772) 570011;
Koshevoy Eugeniy Panteleevich, Doctor of Technology, professor, professor, head of the chair of technologies, machines and equipment of food products, of Kuban State Technological University, 2, Moscovskaya St., Krasnodar, zip: 350072, tel.: (861) 2752279;
Kosachev Vyacheslav Stepanovich, Doctor of Technology, professor, professor of the chair of technologies, machines and equipment of food products, of Kuban State Technological University, 2, Moscovskaya St., Krasnodar, zip: 350072, tel.: (861) 2752279.
DESCRIPTION OF HEAT EXCHANGE IN EMPTY-FIBER MEMBRANES FOR INTERNAL
ADJOINT PROBLEM
As a result of computational modeling integrated kinetic dependency of heat-exchange was received. It allows generalizing the results of computation solution of adjoint problem of heat-exchange from internal stream in empty-fiber membranes and can be applied practically.
Keywords: empty-fiber membranes, adjoint problem, heat-exchane, modeling, kinetic dependency.
Получение обобщенной кинетической зависимости теплообмена на основе расчетной схемы [1] выполнялось методом численного моделирования путем варьирования параметров, входящих в конечно-разностные уравнения: конструктивные - Ятр, Rct, Ьтр, ^тр, Стр, ртр (соответственно, радиусы трубы и стенки, длина трубы, теплофизические свойства трубы - теплопроводность, теплоемкость, плотность) и эксплуатационные - Рнач, Ркон, ^ж, Сж, рж, дж (соответственно, давления в начале и конце
трубы, теплофизические свойства жидкости - теплопроводность, теплоемкость, плотность и вязкость). Эти параметры полностью определяют поле температур в рассматриваемом численном решении. При этом ряд этих параметров может быть сгруппирован в виде комплексов, характеризующих температуропроводность потока жидкости аж = Аж /(Сж ■ рж),
температуропроводность стенки атр = Хтр /(Стр ■ ртр), а также число Прандтля Рг = Сж ■ иж / Яж,
учитывающего влияние теплофизических свойств теплоносителя на теплоотдачу потока жидкости к стенке трубы.
Существенно сократить число исследуемых параметров возможно, используя безразмерные критерии подобия, принятые при описании процессов переноса тепла [2] и учитывающие специфику
Т - Г
гут х од
а а а , частности использовали безразмерную температуру Т =
т - т
іа+ од
( Р - Р
нач кон о2 Ятр
•2-Я
тр
ч 4 ■ /иж ■ Ь .
осевое число Рейнольдса потока жидкости Яе = ------------------------------------------------, число Пекле
Иж / Рж
8 • а • Ах-і- иж • Ьт
Ре = РгЯе^/2, число Фурье Рож =------------ж—--------ж—т, где Ах - шаг интегрирования по оси
Ятр \Рнач Ркон )
потока жидкости; і - номер узла сетки по оси потока жидкости, а также связанное с ним число Греца Ог = 4/ Ро^ , симплекс температуропроводности КА = аж / а , симплекс геометрии
С п Л
•п-я 2
Я
тр
Р - Р
нач кон Т)2
4-и • Ь тр
V г^ж тр J
тр
ТГ тр ^ тт-г V > ж тр / ^
К = п~г——, водяной эквивалент потока жидкости 1фх = -------------------------------------- --- --------------рж • Сж .
1-^-----—г , оидлпип і ішіила /і\ кімі гі ' ' х —
уЦЯСт - Ятр ) 2
Фрагмент температурного поля, полученный при моделировании теплообмена для полипропиленовых мембран, представлен на рисунке 1.
Рисунок 1. Двумерное температурное поле сопряженной задачи теплообмена в
половолоконной мембране
Проведенное численное моделирование процесса теплообмена в половолоконной мембране показало, что интенсивность теплообмена быстро снижается, температурное поле быстро выравнивается по оси потока и разность температур стремится к нулю. Установлено, что интенсивный теплообмен в половолоконных мембранах наблюдается на активном начальном участке (Ьуаг). Изменение температуры носит нелинейный характер и это показывает, что допущение [3] линейной зависимости изменения температуры на внутренней стенке трубы в условиях конвективного теплообмена является неточным.
Из представленных данных видно, что температурное поле на внутренней стенке имеет излом, характерный для задач сопряженного типа [4], следовательно, расчетная схема достаточно подробно моделирует этот тип процесса теплообмена.
Для оценки этой динамики температурного поля было построено изменение температуры жидкости в пограничном слое (Я=0.09 мм) и на внутренней стенке трубы (Я=0.1 мм) (рисунок 2).
0,35
0,3
0,25
0,2
0,15
0,1
0,05
О о
9
о
а
<& ) ® © @ @ (§ —@—@— @@
8
Ро
10
12
14
16
(Ж=0,09 мм <^=0,1 мм
Рисунок 2. Изменение относительной температуры жидкости в пограничном слое и на внутренней
стенке трубы
Отмечаем, что температуры жидкости в пограничном слое и на внутренней стенке трубы практически совпадают и эффект сопряженности проявляется на коротких расстояниях. Для определения влияния на эти расстояния интенсивного теплообмена и исходных параметров процесса были проведены расчеты при изменении исходных параметров, связанных с гидродинамикой процесса. Для оценки интенсивности теплообмена использовали кинетические параметры процесса, определяемые уравнением [5], связывающего среднюю избыточную температуру с длиной участка трубы с активным теплообменом:
фх =
N
1----- е
N..
1)
0
0
2
4
6
где ^=Ь/Ьтаг - относительная длина активного участка теплообмена, определяемая как отношение текущей длины к общей длине активного участка.
Учитывая, что в представленном выше уравнении два определяемых параметра N и Кх, а также собственно длина активного участка Ьуаг, что значительно меньше числа узлов сетки по оси потока, можно с достаточной точностью определить эти величины по кривой изменения избыточной температуры по длине трубы.
В этом случае для определения указанных выше параметров минимизировали следующий функционал:
п2
(©у -1) -ф
Ах
V ^уаг J
У
2)
где © - средняя избыточная температура.
Типичный вид кривой изменения средней избыточной температуры по длине канала после минимизации функционала представлен на графике (рисунок 3).
О 0-1 ^
Рисунок 3. Изменение средней избыточной температуры по длине канала после минимизации
функционала (2)
Представленный график (рисунок 3) показывает высокую степень согласия результатов численных расчетов и аппроксимации зависимостью (1). Отметим, что полученное значение N близко к нулю. Это соответствует принятым граничным условиям сопряженной задачи [1]. Установленный параметр Ьуаг позволяет рассчитать площадь активного участка теплообмена:
Ааг = 2 'П' Яст ' Ауаг
3)
Установленный в результате минимизации функционала (1) параметр (Ых) с учетом водяного эквивалента ^х) позволяет определить кажущийся коэффициент теплоотдачи:
X X
К =
то 7-г
^ 4)
уаг
Коэффициент теплоотдачи использовали в дальнейшем для расчета критерия Нуссельта по соотношению:
Л
ср
где
Л ср
Л Тр •
П • R 2 • Т — П • R 2 • Т I • О
1чст Ттр 1хТр Ттр) !"^тр_
+ Л
ж
П • ^тр • Ттр) • рж
П • ^ст • Ттр П • ^тр • Ттр^) • Ртр + (П • ^тр • Ттр)^ рж
Учитывая значительные вычислительные трудности при реализации численных расчетов сопряженной задачи, дальнейшее моделирование было направлено на установление критериальных зависимостей между рассмотренными выше параметрами. Интервалы варьирования представлены в таблице 1.
Таблица 1. - Интервалы варьирования критериев сопряженной задачи теплообмена
Параметр Минимум Максимум
Pнач, кПа 172 222
Кь 0,436 0,894
Кд 0,679 5,659
Р-ешах 17,909 103,615
РОж 9,559 17,602
Ыи 1,164 2,078
Ьуаг, м 0,046 0,298
Ых 12,398 20,497
^^тр, мм 0,1 0,15
^^ст, мм 0,15 0,25
Лр, Вт/(м • град) 0,03 0,25
Ре 69 400
02 0,9090 1,6738
На первом этапе моделирования исследовали виляние начального давления и соответственно расхода потока на теплообмен. Обобщая полученные данные, получили зависимости активного участка теплообмена Lvaг в зависимости от критерия Пекле при изменении параметров Кь и Кд (рисунок 4).
Ре
Рисунок 4. Влияние критерия Пекле на длину активного участка теплообмена при изменении
параметров Кь и КА
Анализ данных (рисунок 4) показал, что для различных параметров Кь и Кд зависимости активной длины Lvar представляют собой расходящийся пучок прямых из начала координат, наклон которых ть обусловлен параметрами Кь и Кд.
Данные численного моделирования могут быть обобщены в виде следующей зависимости:
тЬ (КЬ > КА )= (Ь0 + Ь1 'КЬ + Ь2 -Ка + Ь12 • КЬ •КА + ^1 'КЬ + Ь22 ' КА ) (6)
где Ь0= 0,000129; Ь1= 0,0011433; Ь2= 0,00013191; Ь12= -0,000011;
Ьц= -0,00064587 ; Ь22= -0,00001412.
Реальная длина половолоконной мембраны должна быть не больше определяемой длине активного участка теплообмена.
С учетом обработки данных моделирования получена зависимость критерия Нуссельта от критериев Пекле, Греца и симплексов, учитывающих геометрию и теплофизические свойства стенки мембраны:
Мы = 3.66 - 8.08 • Ре~°-2334 - С?-0 1474 . К°Л787 • Ка-°Д33°
^ 7)
Пределы изменения переменных в представленной зависимости: Ре=69...10105;
0г=0,828.. .3,443; Кь=0,436.1,066; Кд=0,679.5,659.
Особенностью уравнения (7) является то, что оно в пределе при Кд——0 в случае высокой проводимости тепла стенкой принимает значение №=3,66, соответствующее известному решению задачи Нуссельта - Гретца [2]. Уравнение (7) адекватно с отклонением 10% описывает различные варианты теплообмена от внутреннего потока в половолоконных мембранах.
В результате проведенной работы можно сделать вывод: полученная в результате
моделирования зависимость (7) обобщает результаты численного решения сопряженной задачи
теплообмена от внутреннего потока в половолоконных мембранах и может быть применена в практических расчетах.
Литература:
1. Схаляхов А.А. Постановка сопряженной задачи переноса тепла через трубчатую мембрану от протекающего в ней потока // Новые технологии. 2009. Вып. 4.
2. Справочник по теплообменникам. В 2 т. Т.1: пер. с англ. М.: Энергоатомиздат, 1987. 560 с.
3. Лыков А.В. Теория теплопроводности. М.: Высшая школа, 1967. 600 с.
4. Протодьяконов И.О., Марцулевич Н.А., Марков А.В. Явления переноса в процессах химической технологии. Л.: Химия, 1981. 264 с.
5. Схаляхов А.А. Мембранная экстракция двуокисью углерода в пищевой технологии. Краснодар, 2007. 159 с. Деп. в ВИНИТИ 12.02.07, № 128-В2007.
