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

ИССЛЕДОВАНИЕ ТЕМПЕРАТУРНОГО РЕЖИМА В ПРОЦЕССЕ НАГРЕВА НЕОДНОРОДНОЙ СРЕДЫ «ТЭН-ПЕСОК-ВОЗДУХ»

Телков М.Г. 1 Тимошенко П.О. 1 Суханов И.А. 1 Наимов А.Н. 1 Синицын А.А. 1
1 Вологодский государственный технический университет
Построена и исследована математическая модель процесса нагрева неоднородной среды «ТЭН-песок-воздух» в виде начально-краевой задачи для уравнения теплопроводности в полярной системе координат с краевыми условиями, учитывающими особенности тепловых процессов на границах неоднородной среды. Приведено решение начально-краевой задачи и предложен алгоритм расчета тепловых характеристик процесса нагрева и регулярного температурного режима. Исследование актуально для разработки компьютерной виртуальной модели процесса нагрева неоднородной среды с целью наглядной демонстрации процесса нагрева за небольшой промежуток времени, а также для расчета соответствующих тепловых характеристик на основе реальных экспериментальных данных. Исследования проведены сотрудниками Вологодского государственного технического университета при поддержке Министерства образования и науки РФ в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009–2013 годы.
начально-краевая задача
собственные значения и собственные функции
стационарный режим
регулярный режим
расчет тепловых характеристик
1. Самарский А.А. Избранные труды. – М.: МАКС Пресс, 2003. – С. 531.
2. Самарский А.А., Михайлов А.П. Математическое моделирование. – М.: Физматлит, 2002. – 320 с.
3. Синицын А.А., Тимошенко П.О. К разработке виртуального лабораторного стенда «Изучение процессов теплопередачи в условиях свободной конвекции» // Вузовская наука – региону: материалы седьмой всерос. науч.-техн. конф., 27 февр. 2009 г. – Вологда: ВоГТУ, 2009. – Т. 1. – С. 162–163.
4. Телков М.Г., Наимов А.Н. Существование и полнота собственных векторов задачи Штурма-Лиувилля в форме Кнезера // Вузовская наука – региону: материалы десятой всероссийской научно-технической конференции. – Вологда: ВоГТУ, 2012. – С. 120–125.
5. Тимошенко П.О., Синицын А.А. Основные результаты разработки виртуального лабораторного стенда кафедры теплогазоснабжения и вентиляции // Молодые исследователи – регионам: материалы Всеросс. науч. конф. студентов и аспирантов. – Вологда: ВоГТУ, 2009. – Т. 1. – С. 151–153.

Статья посвящена исследованию температурного режима в процессе нагрева неоднородной среды, состоящей из термоэлектрического нагревателя (ТЭН), помещенного в стальную цилиндрическую трубу, песка, которым заполнено пространство между внутренней стальной трубой и внешней медной трубой, и воздуха, которым окружена внешняя труба. Данное исследование актуально для исследований, описанных в статьях
[3, 2] по разработке компьютерной виртуальной модели процесса нагрева неоднородной среды с целью наглядной демонстрации процесса нагрева за небольшой промежуток времени, а также для расчета соответствующих тепловых характеристик на основе реальных экспериментальных данных.

С целью расчета температурного режима в процессе нагрева неоднородной среды «ТЭН-песок-воздух» были проведены три эксперимента с различными электрическими мощностями ТЭН и были произведены замеры температуры с помощью термопары в четырех точках (на поверхности внешней и внутренней труб и в двух внутренних точках) в последовательные моменты времени (через каждые 5 миллисекунд) от начала процесса до момента стабилизации температуры. Показания термопары фиксировались аналого-цифровым преобразователем и обрабатывались компьютерной программой Table Curve 2D. Программа для каждой точки замера выводила свой вариант формулы зависимости температуры от времени (рисунок) с конкретными значениями коэффициентов.

При этом возникают следующие вопросы:

1. Как связаны между собой коэффициенты получаемых формул и тепловые характеристики процесса – коэффициенты теплопроводности, теплоемкости, теплопередачи?

2. Можно ли находить коэффициенты теплопередачи от ТЭНа к песку и теплоотдачи песка воздуху?

3. Какова динамика температуры в пространстве и во времени?

Для ответа на эти вопросы используется методология работы [1], состоящая из следующих этапов:

I. Построение математической модели процесса нагрева в виде начально-краевой задачи для уравнения теплопроводности с краевыми условиями, учитывающими неоднородность среды.

II. Решение краевой задачи методом разделения переменных, где необходимо построить базис собственных функций, удовлетворяющих краевым ус-
ловиям.

III. Расчет тепловых характеристик процесса на основе экспериментальных данных посредством формул стационарного и регулярного режимов.

рис_138.tif

Зависимость температуры от времени и вариант формулы этой зависимости
в программе Table Curve 2D

Исследуемый процесс нагрева «ТЭН-песок-воздух» качественно отличается от процесса, рассмотренного в работе [1]. Вследствие этого получается другая начально-краевая задача, где базис собственных функций явно находить невозможно. Существование и полнота собственных функций соответствующей задачи Штурма-Лиувилля доказаны в работе [4]. С их помощью выводится формула решения начально-краевой задачи, и, в частности, формула регулярного режима. Регулярным режимом определяется «кривая разгона» – кривая изменения температуры от начала процесса нагрева до момента стабилизации.

Построение математической модели

Для построения математической модели процесса нагрева неоднородной среды «ТЭН-песок-воздух» приведем следующие предположения и обозначения.

1. Температура есть функция U(t, x) переменных t и x, где t – время, t ≥ 0, x –
расстояние от центра внутренней трубы до точки, где измеряется температура, r ≤ x ≤ R, r > 0, r – радиус внутренней трубы, R – радиус внешней трубы; в проведенных экспериментах r = 10 мм, R = 50 мм. Следовательно, температура зависит лишь от расстояния до центра внутренней трубы и она постоянна вдоль любой цилиндрической поверхности, концентрично расположенной относительно внутренней трубы.

2. Обозначим через U(t, r) температуру на поверхности внутренней трубы в момент времени t, U(t, r + 0) – температуру песка на границе с внутренней трубой, U(t, x) – температуру песка во внутренней точке x, r < x< R, U(t, R – 0) – температуру песка на границе с внешней трубой, U(t, R) – температуру воздуха вблизи внешней трубы. При этом температурой медной трубы пренебрегаем из-за большого значения коэффициента теплопроводности меди.

3. Неоднородная среда характеризуется следующими константами: cr – удельная теплоемкость внутренней (стальной) трубы, cR – удельная теплоемкость внешней (медной) трубы, Q – количество тепла, выделяющееся из единичной площади внутренней трубы за единичный промежуток времени, hr – коэффициент теплопередачи от внутренней трубы к песку, h – коэффициент теплоотдачи песка воздуху, αR – коэффициент теплопотери воздуха вблизи внешней трубы, k, ρ, c – коэффициент теплопроводности, плотность и удельная теплоемкость песка, Eqn150.wmf – коэффициент температуропроводности песка.

Математическая модель процесса нагрева будет состоять из уравнений баланса во внутреннем пространстве песка, на внутренней поверхности трубы, на границах «ТЭН-песок», «песок-воздух», в воздухе вблизи внешней трубы.

Уравнение теплового баланса во внутреннем пространстве песка, согласно закону Фурье, [1], имеет вид:

Eqn151.wmf t > 0, r < x < R, (1)

где правая часть есть оператор Лапласа в полярной системе координат с радиус-вектором х.

На поверхности внутренней трубы тепловой баланс зададим условием

Eqn152.wmf t > 0, (2)

где Eqn153.wmf – количество тепла на нагревание внутренней трубы; Eqn154.wmf – количество тепла, отдаваемое песку.

На границе «ТЭН-песок» условие теплового баланса имеет вид:

Eqn155.wmf t > 0. (3)

А на границе «песок-воздух» имеем:

Eqn156.wmf t > 0. (4)

Тепловой баланс в воздухе вблизи внешней трубы зададим условием

Eqn157.wmf (5)

где Eqn158.wmf – количество тепла на нагревание внешней трубы; αRU(t, R) – количество тепла, теряемое воздухом вблизи внешней трубы.

В начальный момент неоднородная среда имеет температуру окружающей среды:

U(0, x) = u0, r ≤ x ≤ R. (6)

Таким образом, математическая модель процесса нагрева представляет собой начально-краевую задачу для дифференциального уравнения теплопроводности (1) в полярной системе координат с краевыми условиями (2), (3) – на внутренней границе, (4), (5) – на внешней границе, и начальным условием (6). Решением U(t, x), t ≥ 0, r ≤ x ≤ R начально-краевой задачи (1)–(6) определяется температурный режим в процессе нагрева неоднородной среды «ТЭН-песок-воздух». Функция U(t, x) при каждом t > 0 терпит разрыв в точках x = r и x = R .

Стационарный режим

Стационарным режимом, как в работе [1], называем стабилизацию температуры в процессе нагрева:

Eqn159.wmf r ≤ x ≤ R,

где Eqn160.wmf – температура на расстоянии x от центра ТЭНа при стабилизации. В формулах (1)‒(5), переходя к пределу при t → +∞, выводим, что функция Eqn160.wmf должна удовлетворять следующим условиям:

Eqn161.wmf r < x < R,

Eqn162.wmf

Eqn163.wmf

Eqn164.wmf

Eqn165.wmf

Отсюда находим Eqn160.wmf:

Eqn166.wmf (7)

По формуле (7) на основе экспериментальных данных можно находить значения тепловых констант αR, k, hr, hR. А именно, пусть известны значения температуры в четырех точках замера при стабилизации:

Eqn167.wmf, Eqn168.wmf Eqn169.wmf, Eqn170.wmf

где r < r2 < r3 < R. Полагая

Eqn171.wmf,Eqn172.wmf

Eqn173.wmf Eqn174.wmf

относительно неизвестных αR, k, hr, hR, получим следующую систему линейных урав-
нений:

Eqn175.wmf

Находим неизвестные:

Eqn176.wmf (8)

Решение начально-краевой задачи

Решение задачи (1)–(6) находим в следующем виде:

Eqn177.wmf (9)

где числа 0 < λ1 < λ2 < ... < λn < ... и функции определяются как собственные значения и собственные функции задачи Штурма-Лиувилля в форме Кнезера [1]:

Eqn178.wmf (10)

В работе [4] доказаны существование и полнота собственных функций задачи (10) в пространстве функций

Eqn179.wmf

со скалярным произведением

Eqn180.wmf

Другими словами, любая функция f(x) ∈ H единственным образом представима в виде ряда по собственным функциям задачи (10):

Eqn181.wmf

где Bi = ⟨f, Vi⟩, ⟨Vi, Vi⟩ = 1, i = 1, 2, ...
и ⟨Vi, Vj⟩ = 0 при i ≠ j. При этом для первого собственного значения λ1 имеет место неравенство:

Eqn182.wmf

Полагая Eqn183.wmf, находим
коэффициенты Ai, i = 1, 2, ... представ-
ления (9):

Eqn184.wmf i = 1, 2, ... . (11)

В этом случае ряд в правой части формулы (9) при каждом t > 0 сходится в обычном смысле.

Регулярный режим

Согласно терминологии работы [5], регулярный температурный режим определяется приближенной формулой

Eqn185.wmf

где λ1 – первое собственное значение, а V1(x) – соответствующая нормированная собственная функция задачи Штурма-Лиувилля (10); их можно находить приближенно численно-аналитическими методами. Коэффициент A1 определяется формулой

Eqn186.wmf

Таким образом, можно задать следующий алгоритм расчета регулярного температурного режима в процессе нагрева неоднородной среды «ТЭН-песок-воздух»:

1. По экспериментальным данным Eqn167.wmf, Eqn168.wmf, Eqn169.wmf, Eqn170.wmf, полученным в ходе замеров температуры в четырех точках при стабилизации, находим приближенные значения тепловых констант с помощью формул (8).

2. Численно-аналитическими методами находим приближенные значения первого собственного значения λ1 задачи Штурма-Лиувилля (10) и соответствующей нормированной собственной функции V1(x).

3. Находим стационарный режим Eqn160.wmf формулой (7).

4. Вычисляем коэффициент A1 формулой

Eqn187.wmf

5. Регулярный режим рассчитывается формулой

Eqn188.wmf r ≤ x ≤ R , t ≥ 0. (12)

Рецензенты:

Горбунов В.А., д.ф.-м.н., профессор, главный специалист ООО НПФ «ЭнергоКИТ», г. Вологда;

Калягин Ю.А., д.т.н., профессор, главный конструктор Общества с ограниченной ответственностью Научно-производственный центр «Информационные и энергетические технологии» (ООО НПЦ «Инэнтех»), г. Вологда.

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


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

Телков М.Г., Тимошенко П.О., Суханов И.А., Наимов А.Н., Синицын А.А. ИССЛЕДОВАНИЕ ТЕМПЕРАТУРНОГО РЕЖИМА В ПРОЦЕССЕ НАГРЕВА НЕОДНОРОДНОЙ СРЕДЫ «ТЭН-ПЕСОК-ВОЗДУХ» // Фундаментальные исследования. – 2012. – № 11-2. – С. 458-462;
URL: http://www.fundamental-research.ru/ru/article/view?id=30558 (дата обращения: 14.11.2019).

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

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