CC BY
21
DOI 10.23859/1994-0637-2017-1-76-11 УДК 004.891
© Синицын Н.Н., Телин Н.В., Плашенков В.В., Андреев А.С., 2017
Синицын Николай Николаевич
Доктор технических наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail: sinitsyn@chsu.ru
Телин Николай Владимирович
Доктор технических наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail: telin_nv@mail.ru
Плашенков Валерий Владимирович
Доктор военных наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail: sinitsyn@chsu.ru
Андреев Александр Сергеевич
Кандидат технических наук, доцент, Череповецкое высшее инженерное военное училище радиоэлектроники (Череповец, Россия) E-mail: sinitsyn@chsu.ru
РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ РАСЧЕТА ТЕМПЕРАТУРНОГО ПОЛЯ ПРИ ПРОГРЕВЕ ГРАНУЛ ИЗ БИОТОПЛИВА*
Sinitsyn Nikolay Nikolaevich
Doctor of Technical Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail: sinitsyn@chsu.ru
Telin Nikolay Vladimirovich
Doctor of Technical Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail: telin_nv@mail.ru
Plashenkov Valeriy Vladimirovich
Doctor of Military Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail: sinitsyn@chsu.ru
Andreev Aleksandr Sergeyevich
PhD (Technology), Associate professor, Cherepovets higher engineering military school of radio electronics (Cherepovets, Russia) E-mail: sinitsyn@chsu.ru
DEVELOPMENT OF MATHEMATICAL MODEL CALCULATION OF TEMPERATURE FIELD ON HEATING GRANULES OF BIOFUELS
Аннотация. Представлена математическая модель одномерного симметричного процесса прогрева гранул из биотоплива в форме бесконечного цилиндра. Математическая модель включает в себя сквозное уравнение теплопроводности общее для сухой зоны, влажной зоны, двухфазной зоны, содержащей пар и воду, при граничных условиях третьего рода. Используется метод конечных разностей, явная схема аппроксимации производных. Представлена методика расчета поля температур и координаты фронта испарения влаги. Представлено тестирование модели реального процесса методом оценки баланса теплоты в теле.
Summary. In the article a mathematical
model of a one-dimensional symmetrical heating process the granules of the biofuel in the form of an infinite cylinder is proposed. Mathematical model comprises a through heat equation , general for the dry zone, wet zone and two-phase zone, that contains the vapor and water, in the boundary conditions of the third kind. It uses the finite difference method, the explicit scheme of approximation of derivatives. The technique for calculating the temperature field and coordinates of the front vapor are proposed. Testing the model of the actual process by evaluating the heat balance in the subject is proposed.
*
В работе принимали участие М.С. Виноградова, Э.В. Гаркавченко, Т.В. Гневашева, В.П. Кузнецова, А.И. Павлова.
Ключевые слова. Прогрев, фазовый пе- Keywords. Warming, the phase transition, реход, гранулы, биотопливо. pellets, biofuels.
Введение
Энергетическое пользование древесного топлива открывает большие возможности в плане охраны окружающей среды и создания экологически чистых производств. Древесные топливные гранулы - облагороженное топливо, полученное путем прессования сухой размолотой первичной или вторичной древесины, коры или их смеси в различных пропорциях. По стандартам различных стран диаметр гранул составляет от 6 до 12 мм, а длина от 5 до 30 мм. Однородность гранулометрического состава гранул позволяет обеспечить полную локализацию и автоматизацию всех циклов технологического процесса котлоагрегата с меньшими затратами. Оценка температурного поля гранулы при ее прогреве необходима при расчетах выгорания топлива, оптимизации аэродинамики топочной камеры с целью изменения выноса недогоревших частиц при проектировании способа сжигания топлива, горелочного устройства. Прогрев частиц твердого топлива в камере сгорания происходит за счет конвективного и радиационного теплообмена между частицей и потоком горячих газов, которые с помощью тех или иных конструктивных мероприятий подводятся к частице.
Основная часть
Расчет нестационарной температуры по сечению гранулы является актуальной задачей. Для изучения закономерностей охлаждения гранул разработано математическое описание процесса. При разработке описания принято во внимание, что гранула имеет форму бесконечного цилиндра. Охлаждение осуществляется принудительной конвекцией холодного воздуха. Математическая модель одномерного симметричного процесса охлаждения включает сквозное уравнение теплопроводности, общее для зоны содержащей пар влажной зоны, и зоны конденсации влаги
дТ д „ дТ 1 дТ
Сэф -р--= — (1—) +---
дг дг дг r дг
(1)
интегрируемая в области: 0 < г < Б,
- начальное условие:
- граничные условия: при г = 0
при г = s
о < т < Тк; Т (г,0) = Т °
>Л = о
дг
= а(Т " ТсР )
(2)
(3)
(4)
где а - коэффициент теплоотдачи, Тср - температура среды; Т - текущая температура, г - текущий радиус; 1 - коэффициент теплопроводности; р - плотность; сэф - эффективная удельная теплоемкость; 5 - радиус цилиндра; т и тк - текущее и конечное время.
В зоне конденсации влаги выделяется теплота, которую можно учесть в уравнении (1) с помощью эффективной теплоемкости сэф:
Сэф _
С (Т), Т < Гс; С2 (Тл)(1 -у) + с2 (Тс)у + |Ь-,Тс < Т < Тл; С2 (Т), Т > Тл
где Тл = Тф + АТ /2, Тс = Тф - АТ /2 - фиктивные температуры начала и окончания
конденсации влаги в элементарном объеме; ДТ - фиктивный интервал температуры конденсации воды; с(Т) - теплоемкость материала; с1 и с2 - теплоемкости влажного и содержащего пар материала; g - доля влаги в элементе объема; Ь - удельная теплота конденсации влаги; у - доля влажного материала; Тф - температура фазового перехода.
Коэффициент теплопроводности и плотность определим по формулам:
Х =
Р =
А.1, Т < Тс
X (1 -у) + Х2^, Тс < Т <Тл
X2, Т >Тл
Р1, Т < Тс;
Р1 (1 -^) + р2Тс <Т <Тл;
Р2, Т > Тл
где р1 и р2 - плотности влажного и содержащего пар материала; и Х2 - коэффициенты теплопроводности влажного и содержащего пар материала. Величина у определяется по формуле:
у =
1, Т < Тс; Т - Т
Т„ - Тс
Т < Т < Т
0, Т > Тл
На рис. 1 показана схема расчетной области конденсации влаги в грануле. Система уравнений (1)-(4) решена численным методом. При использовании метода конечных разностей значения температуры находят в узлах расчетной области, координаты которых рассчитываются по формуле: г =(/ - 0,5)-Аг, для дискретных
моментов времени т = Ах - п , где I = 0, N +1, N - количество узлов внутри расчетной области, 0 и N + 1 - номера фиктивных узлов, находящихся за пределами об-
Аг А л ~
ласти на расстоянии —; Аг = ^ ~ расстояние между узлами; Дт - расчетный шаг по
времени. Для краткости температуры Т(г, тп) обозначают Т" .
Рис. 1. Схема расчетной области 1 - сухая зона; 2 - влажная зона; 3 - двухфазная зона, содержащая пар и воду; ес, 8ф и ел - коэффициенты изотерм начала испарения влаги, фазового перехода влаги и окончания испарения влаги
При использовании явной схемы аппроксимации производных по координате температуру в N внутренних узлах в момент времени п + 1 определяют по зависимости:
Тп+1 — Тп
Дг
с ( ).р(т;п ).дг2
х. 1 (+1 - Т )-х. 1 (( - т--1)
гн— 2
г— 2
п
где г — 1, N , X 1 — X
гн— 2
7^п . грп Т+1 + Тг
X 1 — X
I— 2
тПг + Тп
Температуру в начальный момент времени задают по формуле: Ti — Т , для
г — о, N +1.
Температуру в фиктивных узлах в момент времени п+1 определяют по выражениям:
То — Т1; ^+1 —
(1 -х) TN +2ХТС1
аДг
1+ Х ' " 2Х
Координату фазового перехода воды определяют по координате изотермы конденсации влаги в поле температуры в цикле по г — 2..N из условия:
Т-1 - Тф
Если Тг-1 < Тф < Т ,то 8ф — Дг I г — I + Дг
Т-1 - Тг
Численное решение при явной схеме аппроксимации является условно устойчи-
дг2
вым. В этом случае расчетный шаг определяем по формуле: Дт —-, где ку > 4 ;
куа
а - коэффициент температуропроводности.
Погрешность численного решения в данном случае будет зависеть от настроечных параметров алгоритма N, ку и ДТ. Необходимо эти параметры выбрать таким
образом, чтобы погрешность результатов моделирования не превосходила данную величину.
Тестирование модели реального процесса производим методом оценки погрешности баланса тепла в теле в процессе моделирования реального объема [1]. Для этого составляется уравнение баланса в момент времени т: для цилиндра длиной 1 м Qщ = Qт ,где Qпр - количество теплоты, пришедшее в тело за время т; QТ - количество теплоты, оставшееся в теле к моменту времени т.
Далее конкретизируем приведенные величины для случая исследования тепловых процессов в бесконечном цилиндре:
Qпр = 1^2™^ и Qт = с-р}(г— Та )).
о
При использовании метода конечных разностей интегралы вычисляются численно. В связи с этим левая и правая части баланса теплоты будут отличаться. Для
Юпр — &
оценки погрешности моделирования определим: е =J-• 100% .
Qпр
При варьировании количества узлов погрешность будет изменяться. Допустим, что е = / (Ж), можно изучить эту зависимость и определить оптимальное количество узлов, при котором погрешность баланса теплоты не превзойдет допустимую погрешность.
Формулы для определения Qпр и QТ в конечно-разностной форме при Х = Х(Т) имеют вид:
Qпр = | >-(Т )Т^——ГТж • 2»Дт = ДТ £ I () ( „ — Тк ) 22
N N 2
Qт =с•р^Я£(Т — Т0) — 1)2 Дг2 =с•р^я£(Т( — Тср)[(/ — 1)Дг]2,
I=1 I=1
где N - количество внутренних узлов в расчетной области; N + 1 - фиктивный узел; Дт, Дг - расчетный шаг по времени и координате; 1 - индекс узла; п - индекс временного шага; п* - индекс, соответствующий текущему моменту времени т* = Дт • п * .
При тестировании исследовали влияние настроечных параметров конечно-разностного решения задачи Стефана на результаты и погрешность моделирования. Исследование влияния количества узлов N на погрешность моделирования температурного поля приведено на рис. 2 при настроечном коэффициенте Ку = 4.
%
20
15
10
ЛТ=2
N
0 20 40 60
100 120 140 160
Рис. 2. Зависимость погрешности баланса теплоты от количества узлов при ДТ=2°С
5
0
Как видно из рис. 2, зависимость е = ДЫ) имеет экстремальный характер. С увеличением N погрешность сначала уменьшается, затем увеличивается. При настроечном параметре ДТ = 2 °С оптимальное число узлов N = 60.
Исследование влияния фиктивного интервала АТ (рис. 3) показало, что имеется оптимальное значение ДТ = 2 °С, при котором численное решение по динамике фронта испарения влаги имеет наименьшее значение средней относительной погрешности баланса теплоты.
Рис. 3. Зависимость относительной погрешности баланса теплоты от ДТ (параметра модели: N = 60, Ку = 4)
Уменьшение А Т приводит к тому, что количество узлов, попадающих в фиктивный интервал (Д1 на рис. 1), становится недостаточным для аппроксимации источника испарения влаги. При увеличении ДТ>2°С погрешность повышается, так как увеличивается различие характера температурного поля в окрестности коэффициента фазового перехода.
Вывод
Предложена математическая модель процесса прогрева гранулы и методика тестирования конечно-разностного решения задачи прогрева гранулы с учетом фазового перехода воды (испарение влаги).
Литература
1. Грибкова Ю.В., Габелая Д.И., Кабаков З.К. Разработка методики оценки статей баланса теплоты в заготовке, отливаемой на МНЛЗ // Международный журнал прикладных и фундаментальных исследований. 2016. №1. С. 132-134.
2.Синицын Н.Н., Кабаков З.К., Степанова А.В., Малинов А.Г. Модель замораживания железорудного концентрата // Вестник Череповецкого государственного университета. 2013. №2. Т.1. С. 19-22.
3.Синицын Н.Н., Кабаков З.К., Домрачев Д.А. Математическая модель сушки коры деревьев при высокоинтенсивном нагреве // Вестник Череповецкого государственного университета. 2013. №2. Т.2. С. 24-28.
References
1. Gribkova Yu.V., Gabelaya D.I., Kabakov Z.K. Razrabotka metodiki otsenki statey balansa teploty v zagotovke, otlivayemoy na MNLZ [The development of methodology for assessing articles of the heat balance in the workpiece, the cast on the caster]. Mezhdunarodnyy zhurnal prikladnykh i fundamental'nykh issledovaniy [International Journal of Applied and Basic Research], 2016, no.1, pp. 132-134.
2. Sinitsyn N.N., Kabakov Z.K., Stepanova A.V., Malinov A.G. Model' zamorazhivaniia zhele-zorudnogo kontsentrata. [Model of freezing of iron ore concentrate]. Vestnik Cherepoveckogo gosu-darstvennogo universiteta [Bulletin of the Cherepovets State University], 2013, no. 2, pp. 19-22.
3.Sinitsyn N.N., Kabakov Z.K., Domrachev D.A. Matematicheskaia model' sushki kory derev'ev pri vysokointensivnom nagreve. [A mathematical model of drying bark at high intensity heat]. Vestnik Cherepoveckogo gosudarstvennogo universiteta [Bulletin of the Cherepovets State University], 2013, no. 2, pp. 24-28.
Синицын Н.Н., Телин Н.В., Плашенков В.В., Андреев А.С. Разработка математической модели расчета температурного поля при прогреве гранул из биотоплива // Вестник Череповецкого государственного университета. 2017. №1. С. 81-86.
Sinitsyn N.N., Telin N.V., Plashenkov V.V., Andreev A.S Development of mathematical model calculation of temperature field on heating granules of biofuels. Bulletin of the Cherepovets State University, 2017, no. 1, pp. 81-86.
