Уравнение теплопроводности¶

Уравнение теплопроводности — это дифференциальное уравнение в частных производных параболического типа, описывающее процесс распространения тепла в сплошной среде (твёрдом теле, жидкости или газе) во времени и пространстве. В математической физике оно является базовой моделью для процессов диффузии и теплопереноса, а также для широкого класса явлений, подчиняющихся аналогичным законам (диффузия вещества, фильтрация жидкости, распространение нейтронов в ядерных реакторах). Уравнение устанавливает связь между изменением температуры в данной точке и потоком тепла через окружающую область.

¶Математическая формулировка

¶Общий вид

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

\[ \frac{\partial u}{\partial t} = a \nabla^2 u + f(\mathbf{x}, t) \]

где:

• \( u(\mathbf{x}, t) \) — температура в точке с координатами \(\mathbf{x}\) в момент времени \(t\);
• \( t \) — время;
• \( a = \frac{k}{\rho c_p} \) — коэффициент температуропроводности (м²/с), где \(k\) — теплопроводность, \(\rho\) — плотность, \(c_p\) — удельная теплоёмкость при постоянном давлении;
• \( \nabla^2 \) — оператор Лапласа (в декартовых координатах: \(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}\));
• \( f(\mathbf{x}, t) \) — объёмная плотность источников тепла (например, от химических реакций или электрического тока).

¶Одномерный случай

Для стержня или тонкой пластины, где температура зависит только от одной пространственной координаты \(x\), уравнение упрощается:

\[ \frac{\partial u}{\partial t} = a \frac{\partial^2 u}{\partial x^2} + f(x, t) \]

¶Стационарное уравнение

Если процесс не зависит от времени (установившийся режим), \(\frac{\partial u}{\partial t} = 0\), и уравнение переходит в уравнение Лапласа (при отсутствии источников) или Пуассона (при наличии источников):

\[ \nabla^2 u = -\frac{f}{a} \]

¶Физический смысл

Уравнение теплопроводности является прямым следствием закона сохранения энергии (первого начала термодинамики) и закона теплопроводности Фурье. Согласно закону Фурье, плотность теплового потока \(\mathbf{q}\) пропорциональна градиенту температуры:

\[ \mathbf{q} = -k \nabla u \]

Знак минус указывает, что тепло переносится от более нагретых областей к менее нагретым. Подстановка этого выражения в уравнение баланса тепла (изменение внутренней энергии равно сумме притока тепла через границы и выделения тепла источниками) приводит к дифференциальному уравнению второго порядка.

¶Классификация и виды

¶По размерности пространства

• Одномерное — описывает нагрев стержня, пластины или проволоки.
• Двумерное — для тонких пластин, мембран или плоских слоёв.
• Трёхмерное — для объёмных тел произвольной формы.

¶По наличию источников

• Однородное (\(f = 0\)) — тепло распространяется только за счёт теплопроводности, без внутренних источников.
• Неоднородное (\(f \neq 0\)) — учитывает внутреннее тепловыделение (например, при радиоактивном распаде, химической реакции или прохождении электрического тока).

¶По нелинейности

• Линейное — коэффициенты \(a\) и \(k\) не зависят от температуры.
• Нелинейное — теплофизические свойства среды (теплопроводность, теплоёмкость) зависят от температуры, что приводит к более сложным уравнениям, часто решаемым численно.

¶Граничные и начальные условия

Для однозначного решения уравнения теплопроводности необходимо задать:

¶Начальное условие (условие Коши)

Распределение температуры в начальный момент времени \(t = 0\):

\[ u(\mathbf{x}, 0) = u_0(\mathbf{x}) \]

¶Граничные условия

Определяют тепловое взаимодействие тела с окружающей средой. Различают три основных типа:

1. Граничное условие первого рода (Дирихле) — задана температура на поверхности тела: \(u|_{\partial \Omega} = g(\mathbf{x}, t)\).
2. Граничное условие второго рода (Неймана) — задан тепловой поток через поверхность: \(-k \frac{\partial u}{\partial n}|_{\partial \Omega} = q(\mathbf{x}, t)\).
3. Граничное условие третьего рода (Робина) — задан закон теплообмена с окружающей средой (закон Ньютона — Рихмана): \(-k \frac{\partial u}{\partial n}|_{\partial \Omega} = h (u|_{\partial \Omega} - u_{\text{окр}})\), где \(h\) — коэффициент теплоотдачи.

¶Методы решения

¶Аналитические методы

• Метод разделения переменных (Фурье) — позволяет представить решение в виде ряда по собственным функциям оператора Лапласа. Применим для тел простой геометрии (прямоугольник, круг, шар, цилиндр) и линейных граничных условий.
• Метод функций Грина — использует фундаментальное решение уравнения (тепловое ядро) для построения решения при произвольных начальных и граничных условиях.
• Интегральные преобразования — преобразование Фурье по пространственным переменным или преобразование Лапласа по времени сводят уравнение к алгебраическому или обыкновенному дифференциальному уравнению.

¶Численные методы

• Метод конечных разностей — дискретизация производных по времени и пространству. Наиболее распространённые схемы: явная (Эйлера), неявная (Кранка — Николсона) и схемы переменных направлений.
• Метод конечных элементов — разбиение расчётной области на треугольные или четырёхугольные элементы. Широко применяется в инженерных расчётах (ANSYS, COMSOL, ABAQUS).
• Метод граничных элементов — сводит задачу к интегральным уравнениям на границе области, что эффективно для задач с бесконечными или полубесконечными областями.

¶Применение

Уравнение теплопроводности является фундаментальным инструментом в теплофизике, материаловедении, геофизике, астрофизике и технике. Оно используется для:

• расчёта температурных полей в деталях машин, двигателях, электронных компонентах;
• моделирования тепловых режимов зданий, теплоизоляции и систем отопления;
• анализа процессов закалки, отжига и сварки металлов;
• прогнозирования распространения тепла в грунте, горных породах и ледниках;
• расчёта тепловых потоков в ядерных реакторах и тепловыделяющих элементах;
• моделирования климатических процессов и теплового баланса Земли.

¶История

Первые теоретические основы теплопроводности были заложены французским математиком и физиком Жаном Батистом Жозефом Фурье (1768–1830). В 1822 году он опубликовал монографию «Аналитическая теория тепла», в которой впервые вывел дифференциальное уравнение теплопроводности и разработал метод разделения переменных для его решения. Работы Фурье стали основой для развития теории рядов и интегралов Фурье, а также для всей математической физики XIX века. Впоследствии уравнение обобщалось на случай анизотропных сред, нелинейных процессов и сред с фазовыми переходами (задача Стефана).

¶Интересные факты

• Уравнение теплопроводности является частным случаем уравнения диффузии (второго закона Фика) — математически они идентичны, если заменить температуру на концентрацию, а температуропроводность — на коэффициент диффузии.
• В квантовой механике уравнение Шрёдингера по форме напоминает уравнение теплопроводности, но с мнимой единицей перед производной по времени, что приводит к волновым свойствам решения.
• При очень низких температурах (ниже 10 К) классическое уравнение теплопроводности перестаёт работать из-за квантовых эффектов (баллистический теплоперенос фононов).
• В задачах с фазовыми переходами (плавление, затвердевание) граница раздела фаз движется, и её положение определяется дополнительным условием Стефана — это так называемая задача Стефана, которая является нелинейной и решается, как правило, численно.

¶Источники

• Фурье Ж. Б. Ж. Аналитическая теория тепла. — М.: Наука, 1968.
• Тихонов А. Н., Самарский А. А. Уравнения математической физики. — М.: Изд-во МГУ, 1999.
• Карслоу Г., Егер Д. Теплопроводность твёрдых тел. — М.: Наука, 1964.
• Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Том 6: Гидродинамика. — М.: Физматлит, 2001 (глава о теплопроводности).
• Патанкар С. В. Численные методы решения задач теплообмена и динамики жидкости. — М.: Энергоатомиздат, 1984.