CC BY
17
УДК 620.93
Фрактальная модель слоя твердых бытовых отходов при его переработке
в термическом реакторе
О.Б. Колибаба, Р.Н. Габитов, И.А. Скворцов, Т.В. Муратова ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
г. Иваново, Российская Федерация E-mail: tevp@tvp.ispu.ru
Авторское резюме
Состояние вопроса: При моделировании процесса тепло- и массообмена в термическом реакторе слой твердых бытовых отходов представляют в виде упорядоченной структуры - шаров, цилиндров. В реальности слой является неупорядоченной полидисперсной структурой с различной порозностью и теплофизическими свойствами. Существующие модели недостаточно точно учитывают реальную структуру слоя, что приводит к погрешности 10-15 % и необходимости организации и осуществления большого количества опытов и получения большого количества зависимостей для различных вариантов пористых тел. В связи с этим возникает необходимость разработки модели, позволяющей корректно учесть реальную структуру слоя твердых бытовых отходов.
Материалы и методы: В качестве модели слоя твердых бытовых отходов для исследования тепломассообменных процессов использована структура, представляющая собой фрактальный элемент - пространственный аналог квадратного ковра Серпинского, известный как губка Менгера. При моделировании использованы численные методы с применением программного комплекса Ansys и экспериментальные методы с использованием лабораторных стендов.
Результаты: Предложена геометрическая модель слоя твердых бытовых отходов в виде губки Менгера, реализованная для расчета процесса конвективной сушки слоя твердых бытовых отходов с погрешностью 3,5 % в зоне подсушки термического реактора шахтного типа, позволяющая исследовать процессы тепло- и массообмена в слое твердых бытовых отходов в программном комплексе ANSYS.
Выводы: Предложенная модель может применяться в расчетах при проектировании термических реакторов для переработки многокомпонентных органических отходов, в том числе и твердых бытовых отходов.
Ключевые слова: твердые бытовые отходы, модель слоя, конвективная сушка, фрактал, тепло- и массообмен, губка Менгера.
Fractal model of solid waste layer under recycling in a thermal reactor
O.B. Kolibaba, R.N. Gabitov, I.A. Skvortsov, T.V. Muratova Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: tevp@tvp.ispu.ru
Abstract
Background: Models of heat and mass transfer in a thermal reactor represent household solid waste bed as an ordered structure consisting of balls or cylinders. In fact, the bed is a disordered polydisperse structure with different porosity values and thermal properties. The models in use cannot properly account for the real bed structure, which leads to a 10-15 % error and the need to conduct numerous experiments and to find dependences for a variety of porous bodies. This fact makes it urgent to develop a model that would take a proper account of the real household solid waste bed structure. Materials and methods: A fractal element structure - a three-dimensional analogue of the Sierpinski square carpet known as the Menger sponge - was used as a solid waste bed model to study heat and mass transfer processes. The modeling was conducted by numerical methods in the Ansys software package and by experimental methods on laboratory stands. Results: A geometric model of solid waste bed has been suggested. It represents the bed as a Menger sponge realized for calculating convective drying of the solid waste bed with an accuracy of 3,5 % in the thermal shaft reactor drying area. The model allows investigating the processes of heat and mass transfer in a solid waste bed in the ANSYS software package. Conclusions: The proposed model can be used in calculations in thermal reactor design to process multicomponent organic waste, including household solid waste.
Key words: household solid waste, bed layer, convective drying, fractal, heat and mass transfer, Menger sponge.
DOI: 10.17588/2072-2672.2015.5.032-037
Твердые бытовые отходы (ТБО), подвергаемые переработке в термическом реакторе, представляют собой двухфазную систему «газ -твердое тело» и являются, по существу, пористой структурой с порозностью, лежащей в диапазоне 0,4-0,7. Осуществление термических
процессов характеризуется переносом теплоты и массы вещества как в пределах каждой из фаз, так и превращениями (фазовыми и химическими) на границе раздела. Моделирование процессов нагрева и термического разложения невозможно без постановки и решения краевых
задач нестационарного переноса теплоты и массы вещества.
Существуют два подхода к описанию закономерностей теплообмена и массообмена в пористых телах.
Первый подход предполагает замену реальной пористой структуры ее моделью - условным квазиоднородным изотропным пористым телом с усредненными величинами пороз-ности, температуры, давления и т. д. Такому условному пористому телу приписывают непрерывное температурное поле и условные эффективные теплофизические коэффициенты: эффективный коэффициент удельной теплоемкости сэ и эффективный коэффициент теплопроводности Хэ. Полагают, что температуры газа в порах и материала мало отличаются друг от друга и в расчетах принимаются одинаковыми. Для описания температурного поля пористого тела используют уравнение
% х дt „ д^ СэТ- = рс— + (1)
дх дх дх2
где w - скорость фильтрации теплоносителя, м/с; с - средняя удельная теплоемкость теплоносителя, кДж/кг°С; р - плотность теплоносителя, кг/м3; f - порозность тела.
Во втором случае тепломассообмен описывают системой дифференциальных уравнений для твердого элемента пористого тела (внутренняя задача) и теплоносителя (внешняя задача).
В общем случае решение внешней задачи тепломассообмена основано на совместном рассмотрении уравнений:
• движения Навье - Стокса
1
— + идгас1и = д—дгадР + \>Ч2 и; (2)
дх р
• неразрывности вязкого несжимаемого
потока
др
дх
+ div ри = 0;
(3)
конвективно-диффузионного переноса
Я/'"4
—+ G- gradC = DV2C;
дх
• и уравнения энергии
— + vgradt = aW2t±^,
дх ер
(4)
(5)
где д - ускорение свободного падения, м/с ; р -плотность сушильного агента, кг/м3; Р - давление, Па; V - коэффициент кинематической вязкости, м2/с; С - концентрация, м3/м3; й - коэффициент молекулярной диффузии, м2/с; t - температура, °С; а - коэффициент температуропроводности теплоносителя, м/с ; - мощность внутреннего источника (стока) теплоты.
Решение внутренней задачи тепломассообмена основано на рассмотрении уравнения теплопроводности Фурье (уравнение энергии)
^ + (И, дгаМм) = а дх смрм
и уравнения диффузии дС
С + (и, gradСм) = ОЧ2См.
дх
(7)
Система (2)-(7) дополняется уравнениями граничных условий с учетом известных начальных полей температур, давлений и концентраций.
Использование при расчете температурных полей пористых тел как первого, так и второго подходов вызывает ряд затруднений. Прежде всего, это моделирование самой структуры пористого тела. Обычно при моделировании исходят из того, что твердый «скелет» состоит из тел простейшей формы (шаров, цилиндров и т.п.). Отклонения от известных форм учитывают введением всевозможных поправок, которые бывает трудно оценить [1, 2].
Сложности также вызывает определение эффективного коэффициента теплопроводности в уравнении (1), который, как известно, не является аддитивной величиной, и он должен учесть, кроме кондуктивной составляющей, радиационный и конвективный теплоперенос в порах. Известные расчетные зависимости [3] не только сложны и громоздки, но в своей основе базируются на той же модели пористой структуры из тел простейшей формы.
Широко применяемые на практике экспериментальные и расчетно-экспериментальные методы обратных задач теплопроводности приводят к необходимости организации и осуществления большого количества опытов и получения большого количества зависимостей для различных вариантов пористых структур.
Дополнительные трудности возникают также при определении коэффициента теплоотдачи а, для расчета которого обычно используют критериальные уравнения, полученные с учетом того, что элементами слоя являются тела правильной формы (чаще всего это шары одинакового размера) [4,5].
За последнее время для исследования процессов гидродинамики и тепломассоперено-са активно используется 3D-моделирование, которое позволяет более точно учесть неоднородность структуры тела, его теплофизические характеристики и т.д. Основной трудностью при 3й-моделировании пористых тел является построение в CAD-пакете его структуры.
Задача настоящего исследования заключается в формировании подхода к моделированию хаотической полидисперсной многокомпонентной пористой структуры, позволяющего, с одной стороны, унифицировать структуру, а с другой - учесть ее характерные особенности с последующим распространением подхода на аналогичные пористые тела. Модель должна адекватно отражать структуру реального тела (устойчивость, изотропность, сохранение формы) и корректность математического описания.
Предлагается на базе современных пакетов прикладных программ, с точки зрения удобства математического моделирования пористого тела и протекающих в нем теплообменных процессов, использовать фрактальный подход, который позволяет унифицировать внутреннюю структуру геометрии пористого тела.
Примем в качестве модели реального пористого тела (слоя ТБО) фрактальную структуру -пространственный аналог квадратного ковра Серпинского (triadic Sierpinski carpet), известную как губка Менгера (Menger sponge) [6] (рис. 1).
Рис. 1. Фрактальная модель пористого тела - губка Менгера
Механизм ее построения состоит в следующем. Каждая грань куба, имеющая единичную длину, делится на 9 равных квадратиков так же, как и при построении квадратного ковра Серпинского. В результате исходный куб разбивается на 27 одинаковых кубиков с длиной ребра, равной 1/3. Затем, удаляя 7 кубиков (один центральный и 6 из центра каждой из граней), противоположные грани исходного куба соединяются сквозным центральным отверстием квадратной формы. В результате из 27 остается 20 маленьких кубиков. Такая итерационная процедура с вырезанием сквозных отверстий и последующего превращения каждого оставшегося кубика в 20 еще более мелких (размера в три раза меньше исходного) продолжается до бесконечности. В результате этих операций образуется идеально самоподобный объект, называемый губкой Менгера. Фрактальная размерность губки Менгера равна й = 1п 20 / 1п 3 = 2,7268. Поскольку 2 < й < 3, то это говорит о том, что губка имеет нулевой объем, но обладает бесконечной площадью поверхности своих пор.
Предлагаемая модель отвечает необходимым требованиям. Наличие непрерывной последовательности контактирующих элементов в любом рассматриваемом направлении обеспечивает устойчивость «скелета». Требование изотропности обеспечивает кубическая симметрия пространственной решетки.
Твердый «скелет», пронизанный сеткой пустот, образует структуру второго порядка. Сквозные поры заполнены газовым телом.
В общем случае твердый «скелет» может являться многокомпонентной смесью (структура третьего порядка). В качестве модели твердого «скелета», состоящего из нескольких компонентов, можно использовать одну из известных моделей Г.Н. Дульнева [1].
Предлагаемая фрактальная модель позволяет, с одной стороны, систематизировать внутреннюю геометрию пористого тела, а с другой стороны, внести определенный элемент неупорядоченности, что максимально приближает ее к структуре реального пористого тела (слоя ТБО).
Появляется возможность заменить реальную пористую структуру с определенным диапазоном порозности эквивалентной фрактальной моделью при условии соответствия температурных полей модели и образца, что позволит сократить количество расчетов при моделировании.
Для замены реальной пористой структуры эквивалентной фрактальной моделью необходимо определить требуемый порядок итерации модели. Очевидно, что порядок итерации т будет зависеть от порозности структуры Для построения зависимости была рассчитана пороз-ность модели при нескольких порядках итерации. Искомая зависимость представлена на рис. 2.
0.7
0 6
0 5
0.4
т
дат; fcv1
0.2
Рис. 2. Зависимость порядка итерации модели от порозности структуры
Анализ полученной зависимости (рис. 2) показывает, что, например, при порозности тела, равной 0,5, не получаем целого значения порядка итерации (т = 3,4). Это говорит о том, что при замене реальной пористой структуры фрактальной моделью можно принять модель как 3-го, так и 4-го итерационного порядка. Иными словами, одной этой зависимости недостаточно при выборе эквивалентной модели реальной структуры. Необходимо ввести параметр, который учитывал бы степень развития поверхности пор, приближая, таким образом, модель к реальной структуре. В качестве такого параметра примем эффективный размер который представляет собой отношение объема тела к общей поверхности пор, заключенной в этот объем:
R =
V
(8)
пор
В тепловом отношении параметр будет определять эффективность пористой структуры
(максимальная поверхность теплообмена заключена в минимальном объеме).
Графическая связь между эффективным размером и числом разбиений исходного элемента модели, т е. коэффициентом модификации кт, приведена на рис. 3. Зависимости представлены для трех порядков итерации.
0 5 10 15
кт
Рис. 3. Зависимость эффективного размера от коэффициента модификации кт при различных порядках итерации т
Совместное использование графиков рис. 2 и рис. 3 дает возможность корректно подобрать фрактальную модель реальной пористой структуры. А именно, зная порозность реального пористого тела, по графику рис. 2 находим порядок итерации фрактальной модели. Затем по графику рис. 3, предварительно рассчитав величину для реального тела и зная порядок итерации, определяем коэффициент модификации, т.е. число разбиений (слоев) в модели.
Предлагаемая модель обладает рядом достоинств:
1. Наличие пор различных размеров, распределенных в теле в соответствии с описанной выше итерационной процедурой, позволяет учесть полидисперсность структуры.
2. Геометрия структуры пористого тела (наличие каналов квадратного сечения в губке Менгера) облегчает процедуру нахождения эквивалентного размера пор.
3. Наличие развитой поверхности теплообмена, заключенной в малый объем, характеризует эффективность структуры.
4. Построение геометрии тела с порами прямоугольного сечения (уход от классической губки Менгера) [6] позволяет учесть изменение геометрии структуры пористого тела в результате его термической обработки (например, усадка слоя в процессе сушки в зоне сушки термического реактора).
Фрактальная модель пористого тела на основе губки Менгера использовалась для расчета полей температур и скоростей в слое твердых бытовых отходов, расположенных в
зоне сушки термического реактора шахтного типа [7].
Процесс прогрева и сушки слоя ТБО осуществляется газом, поступающим противотоком из нижерасположенной зоны пиролиза. В результате фильтрации газовый поток полностью заполняет пространство между элементами слоя. Содержащаяся в материале влага отводится из зоны сушки вместе с сушильным агентом.
Постановка задачи формулируется системой (2), (5), (6) в совокупности с уравнениями начальных и граничных условий.
Поскольку перед началом процесса слой ТБО находится в состоянии термодинамического равновесия с окружающей средой, то принимаем начальную температуру слоя равной температуре окружающей среды (20°С). Скорость движения газа в начальный момент времени равна нулю.
Задание граничных условий в данной постановке задачи в явном виде затруднительно, в связи с тем что процесс сушки слоя ТБО является нестационарным. Теплофизические коэффициенты и коэффициенты теплооотдачи зависят от порозности тела и температуры, которая является искомой величиной. По сути, мы имеем дело с нелинейной сопряженной задачей теплообмена и гидродинамики, решение которой практически невозможно ввиду сложности геометрии пористой структуры и трудности нахождения коэффициентов теплообмена.
Выйти из положения позволяют возможности современного математического моделирования. В ПК ANSYS можно корректно решить поставленную задачу, не задавая граничных условий, а задавая лишь входные параметры потоков. Программа, используя внутренний математический аппарат, который заключается в совместном численном решении вышеописанных уравнений, сама определяет необходимые параметры. При этом возникает дополнительная сложность в выборе модели течения флюида (модели турбулентности), поскольку в настоящее время до сих пор не получено аналитическое решение уравнения Навье-Стокса, а каждая отдельно взятая модель имеет определенный диапазон применения.
При моделировании приняты следующие допущения:
1. Исходя из процентной составляющей морфологического состава ТБО, в качестве модели твердого «скелета» пористого тела принята модель трехкомпонентной смеси «с вкраплениями» [1].
2. Теплофизические свойства сушильного агента и материала приняты постоянными, усредненными в интервале рабочих температур.
3. Влияние массообмена на теплообмен учитывалось введением функции внутреннего источника теплоты в уравнении энергии - теплоты, затраченной на испарение влаги в материале. Функция внутреннего источника (стока)
теплоты получена на основе экспериментальных данных и имеет вид qv = 101383 - 5,43 т.
4. Турбулентность движения греющей среды учтена с помощью k-e модели турбулентности.
Реализация модели осуществлена в конечно-элементном пакете ANSYS [8, 9]. Выбор данного продукта обусловлен тем, что он позволяет совместно решать задачи теплообмена, гидравлики и массоперноса при их взаимном влиянии друг на друга.
Построение геометрии модели производилось в среде SolidWorks с последующим импортированием в ПК ANSYS-Fluent. В соответствии с расчетным кодом ANSYS-Workbench создана расчетная сетка конечных элементов с количеством узлов 109524. Как показывает практика расчетов, данной точности достаточно.
В результате решения задачи получены температурные поля газового тела, твердого «скелета» и поле скоростей газа внутри пор, графическая интерпретация которых представлена на рис. 4.
VU^Îwarne 1 а з isfe-OT? 31530*002 3 l-We'UOi1 3 14Ze*002 ЕШ а) щ 3168е«00? зиэгюог з мь.-ео2 э.142е*«в
® з пгЯю N Р X L_ IM
Marty в 7 0Й»Ч1ГС) - у! w
1 1 7ГЭе"000 L «••ч
б)
Л
Рис. 4. Температурное поле скелета в конце всего процесса (а) и газового тела в конце процесса (б) и поле скоростей движения газа (в)
Для проверки адекватности модели использовались результаты лабораторного эксперимента.
На рис. 5 приведена установка ЭЛБ-С-1 для исследования процессов тепломассообмена между слоем ТБО и сушильным агентом. Установка состояла из участка трубы (1) с металлической сеткой для удержания слоя ТБО (рабочая зона) и электрического калорифера (2) для нагрева сушильного агента. Подача сушильного агента осуществлялась вентилятором
(4), расположенным перед калорифером. Скорость потока и его температура регулировались при помощи блока управления (3) с экраном для вывода информации. Температура и влажность сушильного агента фиксировалась при помощи датчиков RHP-2S11 (5) с выводом показаний на экран блока управления. В процессе эксперимента температура образца измерялась пятью термопарами типа ТХА, установленными в слое ТБО. Преобразователь аналого-цифрового сигнала ОВЕН МВА8 передавал сигнал на персональный компьютер. Сушку образца выполняли до установления постоянства его массы.
Рис. 5. Сушильная установка ЭЛБ-С-1
Результаты сравнения расчетных и экспериментальных данных приведены на рис. 6. Относительная погрешность расчета не превышает 3,5 %.
315
т, к
310 305 300 295
2 ЛГ Ш
\
-A V
r 1
320
т, К
315
305 300
?
-
1
О 5000 10000 15000 20000
т, с
Рис. 6. Графики зависимости: а - средней по объему температуры скелета от времени; б - средней по площади температуры газа на выходе из канала от времени; 1 - эксперимент; 2 - расчет на модели
Заключение
Математическая реализация модели осуществлена в программном пакете Ansys Fluent.
Адекватность предложенной модели структуры слоя ТБО на основе губки Менгера подтверждена экспериментальными данными.
Относительная погрешность расчета не превышает 3,5 %.
Список литературы
1. Дульнев Г.Н., Заричняк Ю.П. Теплопроводность смесей и композиционных материалов. - Л.: Энергия, 1974. -264 с.
2. Попов И.А. Гидродинамика и теплообмен в пористых теплообменных элементах и аппаратах. Интенсификация теплообмена. - Казань: Центр инновационных технологий, 2007. - 240 с.
3. Васильев Л.Л., Танаева С.А. Теплофизические свойства пористых материалов. - Минск: Наука и техника, 1971. - 268 с.
4. Кирсанов Ю.А., Назипов Р.А., Данилов В.А. Геометрические и теплофизические характеристики высокопористой структуры // Известия вузов. Авиационная техника. -2010. - № 2.
5. Кирсанов Ю.А., Назипов Р.А., Данилов В.А., Башкирцев Г.В. Математическая модель тепловых процессов и методика исследования теплоотдачи в пористом цилиндре // Известия Самарского научного центра Российской академии наук. - 2010. - Т. 12, № 4. - С. 90-96.
6. Смирнов Б.М. Физика фрактальных кластеров. -М.: Наука, 1991.
7. Пат. №2433344 Российская Федерация МПК F23G. Установка для термического разложения несортированных твердых бытовых отходов / Д.А. Долинин, Р.Н. Габитов, Е.С. Семин, О.В. Самышина, О.Б. Колибаба, О.И. Горинов, В.А. Горбунов. Опубл. 10.11.2011.
8. Бруяка В.А., Фокин В.Г., Кураева Я.В. Инженерный анализ в ANSYS Workbench: учеб. пособие. - Самара: Самар. гос. техн. ун-т, 2013. - 149 с.
9. Басов К.А. Графический интерфейс комплекса ANSYS. - М.: ДМК Пресс, 2006. - 248 с.
References
1. Dul'nev, G.N., Zarichnyak, Yu.P. Teploprovodnost' smesey i kompozitsionnykh materialov [Thermal conductivity of composite materials and mixtures]. Leningrad, Energiya, 1974. 264 p.
2. Popov, I.A. Gidrodinamika i teploobmen v poristykh teploobmennykh elementakh i apparatakh. Intensifikatsiya tep-loobmena [Hydrodynamics and heat transfer in porous heat exchange elements and devices]. Kazan', Tsentr innovatsion-nykh tekhnologiy, 2007. 240 p.
3. Vasil'ev, L.L., Tanaeva, S.A. Teplofizicheskie svoystva poristykh materialov [Thermal properties of porous materials]. Minsk, Nauka i tekhnika, 1971. 268 p.
4. Kirsanov, Yu.A., Nazipov, R.A., Danilov, V.A. Geome-tricheskie i teplofizicheskie kharakteristiki vysokoporistoy struk-tury [Geometric and thermal characteristics of a highly porous structure]. Izvestiya vuzov. Aviatsionnaya tekhnika, 2010, no. 2.
5. Kirsanov, Yu.A., Nazipov, R.A., Danilov, V.A., Bash-kirtsev, G.V. Matematicheskaya model' teplovykh protsessov i metodika issledovaniya teplootdachi v poristom tsilindre [A mathematical model of thermal processes and a technique of heat exchange research in porous cylinders]. Izvestiya Samarskogo nauchnogo tsentra Rossiyskoy akademii nauk, 2010, vol. 12, no. 4, pp. 90-96.
6. Smirnov, B.M. Fizika fraktal'nykh klasterov [Fractal cluster physics]. Moscow, Nauka, 1991.
7. Dolinin, D.A., Gabitov, R.N., Semin, E.S., Samyshina, O.V., Kolibaba, O.B., Gorinov, O.I., Gorbunov, V.A. Ustanovka dlya termicheskogo razlozheniya nesortirovannykh tverdykh bytovykh otkhodov [Unsorted domestic solid waste thermal decomposition unit]. Patent RF, no. 2433344, 2011.
8. Bruyaka, V.A., Fokin, V.G., Kuraeva, Ya.V. Inzhener-nyy analiz v ANSYS Workbench [Engineering analysis in ANSYS Workbench]. Samara, 2013. 149 p.
9. Basov, K.A. Graficheskiy interfeys kompleksa ANSYS [Graphical interface of the ANSYS complex]. Moscow, DMK Press, 2006. 248 p.
Колибаба Ольга Борисовна,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры энергетики теплотехнологий и газоснабжения, e-mail: tevp@tvp.ispu.ru
Габитов Рамиль Наилевич,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», инженер кафедры энергетики теплотехнологий и газоснабжения, e-mail: tevp@tvp.ispu.ru
Скворцов Илья Александрович,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», магистрант 1 -го курса, e-mail: tevp@tvp.ispu.ru
Муратова Татьяна Викторовна,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», магистрант 1 -го курса, e-mail: tevp@tvp.ispu.ru
