CC BY
22
УДК 66.096.5
МЕХАНИКА И ТЕПЛООБМЕН АНСАМБЛЯ ЧАСТИЦ В ПСЕВДООЖИЖЕННОМ СЛОЕ
В.Е. МИЗОНОВ, Л.Н. ОВЧИННИКОВ, доктора техн. наук, А.В. ОГУРЦОВ, канд. техн. наук,
А.В. МИТРОФАНОВ, асп., K. TANNOUS, Dr.-Eng.
Предложена математическая модель и выполнено экспериментальное исследование механики и теплообмена ансамбля частиц в псевдоожиженном слое. Модель основана на теории цепей Маркова, для которой подобраны эмпирические зависимости для расчета коэффициентов переноса. Во всех случаях достигается удовлетворительное совпадение расчетных и экспериментальных данных, это позволяет сделать вывод, что предложенная математическая модель может быть использована как метод расчета параметров псевдоожижения.
Ключевые слова: ансамбль частиц, псевдоожижение, вектор состояния, переходная матрица, теплообмен, экспериментальное распределение температуры.
MECHANICS AND HEAT EXCHANGE OF PARTICLE ENSEMBLE IN FLUIDIZED BED
V.E. MIZONOV, Doctor of Engineering, L.N. OVCHINNIKOV, Doctor of Engineering,
A.V. OGURTSOV, Candidate of Engineering, A.V. MITROFANOV, Post Graduate Student,
K. TANNOUS, Doctor of Engineering,
The article presents a new mathematical model and experimental study of mechanics heat transfer in a particle ensemble in fluidized bed. The model is based on the theory of Markov's chains, for which appropriate correlations to calculate transfer coefficients are selected. In all cases comparison of model predictions with experimental results is made, it allows us to conclude that the model can be used as a computational method to calculate fluidization.
Key words: particle ensemble, fluidization, state vector, transition matrix, heat exchange, experimental temperature distribution.
Технологические процессы, реализуемые в псевдоожиженном слое, часто сопровождаются теплообменом между твердыми частицами и ожи-жающим агентом. И хотя только в ряде случаев назначение самого процесса сводится к теплообмену (например, охлаждение твердого материала после выхода из реакционной камеры), проблема отвода или подвода теплоты к твердым частицам имеет большое значение для осуществления всех типов процесса [1].
Феноменологически теплообмен между ожи-жающим газом и частицами твердого материала описывается уравнением [1]
q = а^Д^ (1)
где а - коэффициент теплоотдачи от ожижающе-го агента к твердой частице (или наоборот); Р -поверхность теплообмена (поверхность твердых частиц); ДТ - разность температур ожижающего агента и твердых частиц; д - количество теплоты, переданное в единицу времени от ожижаю-щего агента к твердым частицам (или наоборот).
На практике же использовать уравнение (1) оказывается затруднительным. Если при расчете простейших поверхностных теплообменников трудность представляет лишь определение коэффициента теплообмена, то в псевдоожиженном слое из-за неравномерности распределения потока ожижающей среды между частицами возникают затруднения и с определением величины температурного напора, а в ряде случаев и поверхности
теплообмена. Это связано с тем, что поверхность теплообмена, величина коэффициента теплоотдачи и разность температур ожижающего агента и твердых частиц определяются локальными характеристиками движения твердой фазы.
Псевдоожиженный слой зачастую является структурно неоднородным, поэтому в каждом малом его объеме в общем случае содержится разное количество частиц и создается свой режим обтекания этих частиц потоком ожижающего агента (причем эти явления взаимно определяют друг друга). Коэффициент теплоотдачи тоже нельзя считать константой для всего слоя, так как он определяется режимом движения жидкости в каждом микрообъеме слоя. В общем случае с уменьшением объемной концентрации частиц твердой фазы коэффициент теплообмена снижается до величины, характерной для чистого потока ожижающей среды.
Однако, в связи с недостаточной изученностью гидромеханических параметров псевдоожижения, при инженерных расчетах делается далеко идущее допущение, что по высоте слоя объемная концентрация твердой фазы постоянна, а значит, постоянна и интенсивность теплоотдачи.
Ранее нами была предложена математическая модель, описывающая поведение твердых частиц в восходящем потоке газа с учетом стесненности обтекания частиц в слое и их взаимодействия друг с другом [2]. Эта математическая модель позволила рассчитывать распределение концентраций частиц и локальных скоростей газа
по высоте аппарата в установившемся режиме и показала удовлетворительную точность прогноза этих характеристик. Результаты численного и натурного экспериментов показали неравномерность распределения концентраций и скоростей частиц по высоте аппарата в установившемся режиме работы. Учет этой неравномерности необходим для более высокой степени достоверности расчетов тепломассообменных процессов в псевдоожиженном слое.
Рис. 1. Расчетная схема процесса
Предлагаемая схема расчета процесса теплообмена между потоком ожижающего агента и частицами твердой фазы показана на рис. 1.
Она состоит из двух параллельных цепей ячеек: одна - для частиц, а другая - для газа. Перемещение частиц рассматривается через малые конечные промежутки времени Д1
Текущее состояние процесса представлено вектором-столбцом Бр, состоящим из содержания материала в каждой ячейке. В течение Ы-го перехода состояние Брм переходит в состояние Бр +1, причем эволюция состояния описывается рекуррентным матричным равенством
БрМ+1=РрМ Бр1'1, (2)
где N - номер перехода; Ррм - матрица переходных вероятностей для частиц материала, механизм формирования которой подробно изложен в работе [2].
При наличии теплообмена между газом и частицами происходит поперечная теплопередача между сходственными ячейками цепей и продольный перенос теплоты вместе с частицами и газом вдоль соответствующих цепей. Распределение температуры вдоль цепей может быть описано векторами температур Тр и Тд, а векторы текущих запасов теплоты в ячейках рассчитаны как Ор=СрБр .*Тр , Од=СдБд .*Тд, где Ср и Сд - теплоемкости частиц и газа; Бр и Бд - содержания частиц и газа в ячейках, в общем случае переменные по высоте слоя; символ «.*» означает поэлементное умножение векторов. Эволюция теплового состояния цепей описывается рекуррентными матричными равенствами, связывающими последовательные состояния цепей через промежуток времени Д1:
ОрМ+1=РрМ(ОрМ+а.*Р.*(ТдМ- Трм)), (3)
ОдМ+1=РдМ(ОдМ - а .*Р.*(Тд'1- Тр")+ О*), (4)
где Од - вектор подачи теплоты с входящим газом при его температуре Тд0, имеющий один ненулевой элемент ОдТд0Сд Д1 для первой ячейки. Матрицы Ррм и Рдм описывают продольное перемещение частиц и газа и их теплоты вдоль цепей. В матрице Ррм содержатся вероятности переходов в течение Д1 из данной ячейки во все
~ .-ч N
другие ячейки цепи, в матрице Рд - доли газа, переносимого в течение этого времени в соседние по ходу ячейки.
Поскольку равенства (3), (4) отражают очевидные уравнения баланса, то получаемые по ним решения всегда будут правдоподобными. Количественная же достоверность расчетных прогнозов напрямую зависит от достоверности определения коэффициентов теплоотдачи в ячейках, которые меняются по высоте слоя.
Для выбора подходящей расчетной зависимости необходимы экспериментальные данные, представляющие именно локальную кинетику прогрева частиц в псевдоожиженном слое. В этой связи нами были выполнены экспериментальные исследования на лабораторной установке (рис. 2), состоящей из воздуходувки (1), подававшей воздух через калорифер 2 в цилиндроконический аппарат с псевдоожиженным слоем 3. Предварительно аппарат путем подачи нагретого воздуха прогревался до установившейся температуры. Затем на решетку мгновенно засыпалась навеска карбамида 300 г с фракциями 23 мм и начиналось измерение температуры газа с интервалом 20 с термопарами, расположенными на высоте 10, 50 и 75 мм над решеткой.
Высота слоя материала при отсутствии по-
Рис. 2. Схема и основные размеры экспериментальной установки
дачи воздуха составляла 45 мм, а во взвешенном состоянии - 75 мм, то есть последняя термопара находилась фактически на верхней границе слоя в его рабочем состоянии. Замеры температуры продолжались до полной стабилизации температуры слоя, равной температуре подаваемого подогретого до 95оС воздуха.
Результаты расчетного и экспериментального исследований переходного процесса прогрева слоя при различных зависимостях для расчета коэффициента теплообмена показаны на рис. 3-5, где линии соответствуют расчетной кинетике прогрева, а точки - экспериментальным данным.
1 1 1 1 1 1 1 1 1111
• ?•••* г-1 Г -1
г / ■
/У' ^ / V .* *
к 1 со ‘ ^ -■ -
1 1 1 1 г***. 1
г .*1 1 / 1 1 *• 1 1 Л- 1 1
• 1 1 • 1 1 1 1
30>1
т, с
0 40 60 80 100 120 140 160 180 200 220 240
Рис. 3. Изменение температуры воздуха в различных сечениях слоя: 1 - 10 мм; 2 - 50 мм; 3 - 75 мм (Ыи = 0,4(Ре/е)2/3Рг1/3 [3])
ta,0C
т, с
При этом критерий Рейнольдса Ре = рассчитывался не для всего слоя, как это принято при традиционных расчетах [3], а для каждой ячейки отдельно (ш - локальная скорость обтекания частиц ожижающим газом, разная в разных ячейках слоя; V - коэффициент кинематической вязкости газа).
Рис. 4. Изменение температуры воздуха в различных сечениях слоя: 1 - 10 мм; 2 - 50 мм; 3 - 75 мм (формула А.В. Чечеткина Ыи=0,123Ре0,83 [3])
После засыпки навески в нагретый газ его температура из-за теплоотдачи к материалу резко снижается, причем снижение выше в верхней части слоя, так как туда поступает уже охлажденный в нижней части слоя газ. Постепенно все температуры повышаются, стремясь к температуре подаваемого нагретого газа, что естественно при имевшейся хорошей теплоизоляции аппарата.
Ключевым вопросом согласования расчетных и экспериментальных данных является выбор эмпирической зависимости для вычисления критерия Нуссельта (Ыи = а5/Х; 5 - средний диаметр ожижаемой фракции; X - коэффициент теплопроводности), который необходим для расчета коэффициента теплоотдачи от газа к частицам а. Принципиально важно, что эти зависимости были взяты из независимых от модели экспериментальных результатов.
Рис. 5. Изменение температуры воздуха в различных
сечениях слоя: 1 - 10 мм; 2 - 50 мм; 3 - 75 мм (Ыи = 0,018Ре0,991 [4])
Согласно полученным результатам (рис. 3-5), все используемые эмпирические зависимости для расчета критерия Нуссельта позволили в количественном отношении довольно точно смоделировать кинетику теплообмена частиц с ожи-жающим агентом.
Однако наибольшее соответствие расчетных и экспериментальных данных достигается при использовании эмпирической зависимости, полученной в работе [4].
Таким образом, предложенная нами ранее математическая модель [2] и эмпирическая зависимость для расчета коэффициента теплоотдачи [4], с одной стороны, являются независимыми, а с другой - обеспечивают хорошее соответствие расчетных и экспериментальных данных, а поэтому может быть рассмотрена как достоверная научная основа компьютерного метода расчета псевдоожиженного слоя.
Список литературы
1. Гельперин Н.И., Айнштейн В.Г., Кваша В.Б. Основы техники псевдоожижения. - М.: Химия, 1967. - 664 с.
2. Митрофанов, А.В., Мизонов В.Е., Огурцов А.В., Таппоив К. Расчетно-экспериментальное исследование аппарата с двумерным псевдоожиженным слоем // Вестник ИГЭУ. - 2009. - Вып. 3. - С. 20-22.
3.Расчеты аппаратов кипящего слоя: Справочник / под ред. И.П. Мухленова, Б.С. Сажина, В.Ф. Фролова. - Л.: Химия, 1986. - 352 с.
4.Овчинников Л.Н. Моделирование процесса сушки минеральных удобрений во взвешенном слое // Изв. вузов «Химия и химич. технология». - 2009. - Т. 52. - Вып. 7. - С. 122-124.
Мизонов Вадим Евгеньевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор, зав. кафедрой прикладной математики, телефон (4932) 26-97-06.
Овчинников Лев Николаевич,
ГОУВПО «Ивановский государственный химико-технологический университет»,
доктор технических наук, профессор кафедры процессов и аппаратов химической технологии,
телефон (4932) 38-22-18.
Огурцов Антон Валерьевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры прикладной математики, телефон (4932) 26-97-06.
Митрофанов Андрей Васильевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», аспирант кафедры прикладной математики, телефон (4932) 26-97-06.
КаИа Таппоиэ,
Университет г. Кампинас,
профессор кафедры химической инженерии, доктор технических наук, каИа@Тед.ип1сатр.Ьг
