Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,749

РЕШЕНИЕ КОНТАКТНОЙ ЗАДАЧИ ТЕПЛОПРОВОДНОСТИ ДЛЯ ПОЛИГОНА ТВЕРДЫХ БЫТОВЫХ ОТХОДОВ

Середа Т.Г. 1
1 ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет»
К физическим процессам, протекающим в массиве твердых бытовых отходов (ТБО), относятся тепло и массоперенос, растворение и диффузия. Движущей силой теплопереноса является разность температур, при этом температура оказывает большое влияние на весь спектр процессов, протекающих на полигоне ТБО. Воздействие температуры в значительной степени определяет характер протекания химико-биологических процессов и влияет на состояние системы в целом. В связи с этим разработана задача теплопроводности в массиве ТБО с использованием метода конечных элементов, включающего в себя дифференциальную и вариационную постановки задачи теплопроводности на полигоне ТБО. Решена задача теплопроводности с использованием метода конечных элементов. Для расчета температурных полей разработан пакет программ «work», предназначенный для определения стационарных и нестационарных температурных полей в плоских областях. В пакете реализован стандартный метод конечных элементов (МКЭ), применяемый для решения различных задач теплопроводности.
твёрдые бытовые отходы (ТБО)
АСУТП
полигон ТБО
1. Артемов Н.И. [и др.] Технологии автоматизированного управления полигоном твердых бытовых отходов / Н.И. Артемов, Т.Г. Середа, С.Н. Костарев, О.Б. Низамутдинов // Международный журнал экспериментального образования. – 2010. – № 11. – С. 43.
2. Костарев С.Н., Середа Т.Г., Михайлова М.А. Системный анализ управления отходами. – LAP LAMBERT Academic Publishing GmbH & Co. KG, 2012.
3. Костарев С.Н. Математическая модель управления состоянием полигона твердых бытовых отходов : дис. … канд. техн. наук : 05.13.18. – Пермь, 2003. – 199 с.
4. Костарев С.Н., Середа Т.Г., Михайлова М.А. Программно-аппаратный комплекс управления качеством фильтрационных стоков // Экологические системы и приборы. – 2014. – № 3. – С. 39–46.
5. Костарев С.Н., Середа Т.Г., Михайлова М.А. Разработка автоматизированной системы мониторинга и управления природно-техническими системами утилизации отходов // Фундаментальные исследования. – 2013. – № 6–2. – С. 273–277.
6. Костарев С.Н., Середа Т.Г. Комплексное внедрение автоматизированных систем управления технологическими процессами на санитарных полигонах твердых бытовых отходов // Экологические системы и приборы. – 2014. – № 4. – С. 20–28.
7. Свидетельство о регистрации программы для ЭВМ 2009612494.– М.: ФИПС, 2009.
8. Середа Т.Г. [и др.]. Снижение пожаровзрывоопасности объектов депонирования отходов // Пожарная безопасность. – 2008. – № 3. – С. 84–89.
9. Середа Т.Г. Обоснование технологических режимов функционирования искусственных экосистем хранения отходов: дис. … д-ра техн. наук : 03.00.16; МГУП. – М., 2006.
10. Середа Т.Г., Костарев С.Н. Разработка методов проектирования автоматизированных систем обработки информации и управления искусственными экосистемами хранения отходов // Экологические системы и приборы. – 2006. – № 11. – С. 21–24.
11. Середа Т.Г. Источники загрязнения и методы защиты человека в техносфере: учеб. пособие для студентов вузов. – Пермь, Изд-во «ПНИПУ», 2012.
12. Сегерлинд Л. Применение метода конечных элементов: пер. с англ. – М.: Мир, 1979. – 390 с.
13. Таблицы физических величин: справочник / под ред. И.К. Кикоина. – М.: Атомиздат, 1976. – 1008 с.
14. Kostarev S.N., Sereda T.G. Аutomated process control of sanitary municipal solid waste landfill // World Applied Sciences Journal. – 2013. – Т. 22. – № SPL.ISSUE2. – С. 64–69.
15. Sereda T.G., Kostarev S.N., Elancheva E.N. Study safety environmental protection landfills using models anaerobic digesters // Applied Mechanics and Materials. – 2014. – Vol. 682. – Р. 339–345.

Математическое моделирование в экологии приобретает все большее значение, что может принести существенные результаты в области моделирования макродлительных экспериментов, прогнозировании и синтезе свойств исследуемой системы. В работе рассматривается моделирование тепловых полей на полигоне захоронения ТБО [5]. Температура является определяющим фактором, влияющим на химические и биологические процессы, протекающие в теле полигона [1, 2].

Целью представленной работы являлось моделирование процессов теплопереноса на полигонах ТБО.

Материал и методы исследования

При обосновании методов теплопереноса использовался метод конечных элементов, методы математического и имитационного моделирования. Объектом исследований являлись тепловые потоки, протекающие на объектах утилизации отходов, рассматриваемые как объекты моделирования и управления. Материалом исследования являлись твёрдые бытовые отходы.

Результаты исследования
и их обсуждение

Дифференциальная постановка задачи. Процесс теплопереноса на полигоне ТБО задан в области W в R3 c границей S, sereda01.wmf sereda02.wmf на части границы Si заданы граничные условия 1-го рода (изолированный экран), на части Sk заданы граничные условия 3-го рода (конвективный теплообмен) (рис. 1) [8, 9].

Уравнение теплопроводности в сплошной среде имеет следующий вид [2]:

sereda03.wmf (1)

где t – температура; КХХ, КУУ, КZZ – коэффициенты теплопроводности в направлениях х, у, z; Q – внутренний источник тепла, при граничном условии

sereda04.wmf (2)

где h – коэффициент теплообмена, кВт/м2; t – температура на границе; tос – температура окружающей среды; q – точечный источник тепла.

pic_7.tif

Рис. 1. Сечение полигона с граничными условиями 1-го и 3-го рода

Для нахождения приближенного решения дифференциальной задачи (1), (2) с помощью таких численных методов, как метод Ритца, Бубнова – Галеркина, конечных элементов (МКЭ), требуется построение вариационных аналогов исходных дифференциальных задач.

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

sereda05.wmf (3)

Для минимизации функционала на множестве узловых значений с различными характеристиками материала объемный интеграл должен быть представлен в виде суммы интегралов, каждый из которых вычисляется по отдельному элементу. Матрица теплопроводности элемента имеет следующий вид [12, 3]:

sereda06.wmf (4)

В качестве дискретного элемента для разбивки используем треугольный элемент. Функция формы для линейного треуголь­ного элемента имеет вид

sereda07.wmf β = i, j, k. (5)

sereda08.wmf (6)

Матрица градиентов [В]:

sereda09.wmf (7)

Матрица свойств материала [D]:

sereda10.wmf (8)

Теперь объемный интеграл примет следующий вид:

sereda11.wmf (9)

Предполагая толщину элемента единичной, заменим dV на dA. Подынтегральное выражение в (9) постоянно и может быть вынесено за знак интеграла:

sereda12.wmf (10)

Вычисляя произведение матриц, имеем

sereda13.wmf (11)

Второй интеграл должен быть вычислен по поверхности:

sereda14.wmf

Подставляя в матрицу [N] функции формы и выполняя мат­ричное умножение, получаем

sereda15.wmf (12)

Программная реализация

Для расчета температурных полей разработан пакет программ «work», предназначенный для определения стационарных и нестационарных температурных полей в плоских областях [4, 6, 7]. В пакете реализован стандартный метод конечных элементов (МКЭ), применяемый для решения различных задач теплопроводности. Функционально пакет разделен на две отдельные части: Grid – для разбивки рабочей области на конечные элементы; Vega – для решения плоской задачи теплопроводности. Рассмотрим основные этапы работы с программой «Work».

расчет теплофизических свойств фракций отходов

Для определения коэффициентов теплоемкости и теплопроводности использованы источники [11, 13] (таблица).

Теплофизические свойства фракций отходов

Компонент

Теплопроводность, кал/(см∙с∙град)

Теплоемкость, ккал/ кг

Бумага

144–308.10–6

3390

Пищевые отходы

130.10–6

890

Дерево

120.10–6

3410

Металл

0,136

49

Стекло

400–470.10–6

23

Резина

500.10–6

5670

Камни

2629.10–6

23

Текстиль

95–120.10–6

3460

Разбивка рабочей области
на конечные элементы

Для разбивки рабочей области на конечные элементы (КЭ) необходимо предварительно разделить ее на зоны, каждая из которых образована восемью опорными узлами. Данную операцию можно провести либо вручную, либо с использованием специальных утилит для автоматического задания координат и номеров опорных узлов. До разбивки следует выполнить эскиз исследуемой области, выбрать необходимое число зон, опорных узлов, занумеровать зоны и узлы, нумерация сторон в каждой зоне ведется против хода часовой стрелки. Информация по разбивке формируется в файлах «Work» с расширениями in1, in2, in3. Далее запускается модуль «Maingrid». Происходит формирование файла граничных условий. Для формирования файла запускается программа «Bc_trm». Исходные данные для «Bc_trm» собирают в текстовом файле «Work.ggg». Расчет температурных полей создается в файле «Work.reg».

Для построения диаграммы температурных полей рис. 2 использована программа AXUM 5.0.

pic_8.tif

Рис. 2. Распределение температурных полей по сечению полигона

Применение методов системного анализа и кибернетики в управлении системами ПТО открывает возможность исследования и организации безопасного функционирования объекта депонирования ТБО в условиях, когда информация последовательно накапливается и реализуется в виде алгоритмов для ЭВМ. В этом случае управляющие действия будут направлены на минимизацию жизненного цикла природно-технических систем утилизации отходов, что имеет важное значение [14]. Внедрение автоматизированного рабочего места конструктора-проектировщика полигона ТБО уже на этапе проектирования позволит ускорить выполнение проектной документации и повысить безопасность эксплуатации природно-технической системы депонирования отходов на этапах её эксплуатации и рекультивации [8, 10, 15].

Рецензенты:

Козлов В.В., д.т.н., профессор кафедры ВМКСиС, Пермский военный институт внутренних войск МВД России, г. Пермь;

Костарев С.Н., д.т.н., доцент, профессор кафедры «Строительное производство и материаловедение», Пермская государственная сельскохозяйственная академия,
г. Пермь.

Работа поступила в редакцию 19.12.2014.


Библиографическая ссылка

Середа Т.Г. РЕШЕНИЕ КОНТАКТНОЙ ЗАДАЧИ ТЕПЛОПРОВОДНОСТИ ДЛЯ ПОЛИГОНА ТВЕРДЫХ БЫТОВЫХ ОТХОДОВ // Фундаментальные исследования. – 2014. – № 12-5. – С. 936-940;
URL: https://fundamental-research.ru/ru/article/view?id=36252 (дата обращения: 24.09.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074