Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,749

SOLUTION OF THE CONTACT HEAT CONDUCTION PROBLEMS FOR SOLID WASTE LANDFILLS

Sereda T.G. 1
1 State National Research Politechnical University of Perm
The physical processes occurring in the array of municipal solid waste (MSW) are heat and mass transfer, dissolution and diffusion. The driving force is the temperature difference between the heat transfer, the temperature has a great influence on the whole spectrum of processes occurring at the landfill. Influence of temperature largely determines the character of the chemical and biological processes and affect the state of the VET system as a whole. In this regard, developed heat conduction problem in an array of MSW using the finite element method, includes a differential and variation formulation of the problem of heat conduction at the landfill. The problem of thermal conductivity using the finite element method. For the calculation of the temperature fields developed software package «Work», designed to determine the stationary and non-stationary temperature fields in flat areas. The package with the standard finite element method (FEM) is used to solve various problems of heat conduction.
Municipal solid waste (MSW)
ASUTP
sanitary landfill municipal solid waste
1. Artemov N.I., Sereda T.G., Kostarev S.N., Nizamutdinov O.B. Tehnologii avtomatizirovannogo upravlenija poligonom tverdyh bytovyh othodov – Mezhdunarodnyj zhurnal jeksperimentalnogo obrazovanija, 2010, no. 11, pp. 43.
2. Kostarev S.N., Sereda T.G., Mihajlova M.A. Sistemnyj analiz upravlenija othodami. LAP LAMBERT Academic Publishing GmbH & Co. KG, 2012. ISBN 978-3-8465-8298-5.
3. Kostarev S.N. Matematicheskaja model’ upravlenija sostojaniem poligona tverdyh bytovyh jekosistem hranenija othodov: dis. … dokt. tehn. nauk: 03.00.16; MGUP. M., 2006.
4. Kostarev S.N., Sereda T.G., Mihajlova M.A. Programmno-apparatnyj kompleks upravlenija kachestvom fil’tracionnyh stokov // Jekologicheskie sistemy i pribory. 2014. no. 3. pp. 39–46.
5. Kostarev S.N., Sereda T.G., Mihajlova M.A. Razrabotka avtomatizirovannoj sistemy monitoringa i upravlenija prirodno-tehnicheskimi sistemami utilizacii othodov // Fundamental’nye issledovanija. 2013. no. 6–2. pp. 273–277.
6. Kostarev S.N., Sereda T.G. Kompleksnoe vnedrenie avtomatizirovannyh sistem upravlenija tehnologicheskimi processami na sanitarnyh poligonah tverdyh bytovyh othodov // Jekologicheskie sistemy i pribory. 2014. no. 4. pp. 20–28.
7. Svidetelstvo o registracii programmy dlja JeVM 2009612494. M.: FIPS, 2009.
8. Sereda T.G. [i dr.]. Snizhenie pozharovzryvoopasnosti obektov deponirovanija othodov // Pozharnaja bezopasnost’. 2008. no. 3. pp. 84–89.
9. Sereda T.G. Obosnovanie tehnologicheskih rezhimov funkcionirovanija iskusstvennyh Sereda T.G. [i dr.]. Snizhenie pozharovzryvoopasnosti obektov deponirovanija othodov // Pozharnaja bezopasnost’. 2008. no. 3. pp. 84–89.
10. Sereda T.G., Kostarev S.N. Razrabotka metodov proektirovanija avtomatizirovannyh sistem obrabotki informacii i upravlenija iskusstvennymi jekosistemami hranenija othodov – Jekologicheskie sistemy i pribory. 2006, no. 4, pp. 21–24.
11. Sereda T.G. Istochniki zagrjaznenija i metody zashhity cheloveka v tehnosfere: Ucheb. posobie dlja studentov vuzov. – Perm’, Izd-vo «PNIPU», 2012.
12. Segerlind L. Primenenie metoda konechnyh jelementov: Per. s angl. M.: Mir, 1979.
13. Tablicy fizicheskih velichin: Spravochnik / Pod red. I.K. Kikoina. M.: Atomizdat, 1976.
14. Kostarev S.N., Sereda T.G. Аutomated process control of sanitary municipal solid waste landfill // World Applied Sciences Journal. 2013. Т. 22. no. SPL.ISSUE2. pp. 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 Vol. 682 (2014). pp. 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.