Уравнение теплопроводности: закон Фурье и решение

Уравнение теплопроводности описывает, как температура распределяется в теле и меняется со временем, а закон Фурье связывает это распределение с потоком тепла через сечение. Вместе эти два соотношения - основа расчёта нагрева и остывания стержней, стен, обмоток и любых других тел, где тепло передаётся не излучением и не конвекцией, а обычной теплопроводностью материала. Ниже разберём, как из закона Фурье получается дифференциальное уравнение в частных производных, что такое коэффициент температуропроводности и почему он определяет скорость выравнивания температуры, и как решить типовую задачу методом разделения переменных. Чтобы сразу увидеть, как решение выглядит на живом примере, покрути калькулятор ниже: он решает уравнение для стержня с горячим пятном посередине и показывает, как это пятно расплывается и остывает.

Закон Фурье: поток тепла и температурный градиент

Закон Фурье - эмпирическое соотношение, установленное ещё в начале XIX века: плотность теплового потока через сечение пропорциональна градиенту температуры и направлена в сторону убывания температуры:

$q = -k\,\frac{\partial T}{\partial x},$

где $q$ - плотность потока тепла (Вт/м²), $k$ - коэффициент теплопроводности материала (Вт/(м·К)), а знак минус отражает физический смысл: тепло течёт от горячего к холодному, то есть в сторону, противоположную направлению роста температуры. Чем круче температурный градиент, тем интенсивнее поток - это прямая аналогия с диффузией вещества, только вместо концентрации переносится тепловая энергия.

Рассчитать для своей задачи

На схеме видно, как перепад температуры $\Delta T$ на коротком участке $\Delta x$ задаёт наклон касательной - именно он и определяет величину потока в этой точке. Там, где линия температуры круче, поток больше; там, где профиль почти горизонтален, поток близок к нулю.

Уравнение теплопроводности: вывод из закона Фурье и баланса энергии

Чтобы получить уравнение, описывающее температуру во времени, применим закон Фурье к тонкому слою стержня толщиной $dx$ и запишем баланс энергии: разница входящего и выходящего потоков идёт на изменение внутренней энергии слоя. После стандартных преобразований (запись потока в сечениях $x$ и $x+dx$, разложение в ряд и деление на $c\rho\,dx$) получается одномерное уравнение теплопроводности:

$\frac{\partial T}{\partial t} = a\,\frac{\partial^2 T}{\partial x^2}, \qquad a = \frac{k}{c\rho},$

где $c$ - удельная теплоёмкость материала, $\rho$ - его плотность, а $a$ называется коэффициентом температуропроводности. Именно $a$, а не сам коэффициент теплопроводности $k$, определяет, насколько быстро выравнивается температура: металл с высокой теплопроводностью, но и высокой теплоёмкостью может прогреваться медленнее, чем кажется на первый взгляд, - важна именно комбинация всех трёх свойств.

Решение методом разделения переменных: ряд Фурье

Для стержня длиной $L$ с концами, закреплёнными на нулевой температуре ($T(0,t)=T(L,t)=0$), решение ищут в виде суммы гармоник, каждая из которых удовлетворяет уравнению и граничным условиям по отдельности:

$T(x,t) = \sum_{n=1}^{\infty} b_n \sin\frac{n\pi x}{L}\, e^{-a(n\pi/L)^2 t},$

где коэффициенты $b_n$ находятся разложением начального профиля $T(x,0)$ в ряд Фурье по синусам:

$b_n = \frac{2}{L}\int_0^L T(x,0)\sin\frac{n\pi x}{L}\,dx.$

Если начальный профиль уже задан одной гармоникой, $T(x,0) = T_0\sin(\pi x/L)$, всё упрощается: единственный ненулевой коэффициент $b_1 = T_0$, и решение принимает компактный вид $T(x,t) = T_0\sin(\pi x/L)\,e^{-t/\tau}$, где $\tau = L^2/(a\pi^2)$ - характерное время затухания. Именно этот частный случай удобен для быстрой оценки: амплитуда убывает по чистой экспоненте, а форма профиля не меняется.

Рассчитать для своей задачи

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

Почему высшие гармоники Фурье гасятся быстрее и сглаживают резкие перепады

Ключевое наблюдение решения - показатель экспоненты $-a(n\pi/L)^2 t$ растёт как квадрат номера гармоники $n$. Это значит, что вторая гармоника гаснет вчетверо быстрее первой, третья - в девять раз быстрее, и так далее. Именно поэтому резкий прямоугольный перепад температуры, который математически раскладывается в ряд из бесконечного числа гармоник, теряет свои острые углы почти мгновенно: высокочастотные составляющие, отвечающие за резкость краёв, исчезают первыми, а низкочастотная плавная «шапка» остаётся дольше всех.

Рассчитать для своей задачи

Этот же эффект объясняет, почему в калькуляторе выше карта распределения температуры быстро теряет чёткие границы горячей зоны: за первые доли характерного времени $\tau$ пропадают именно резкие грани, а дальше остаётся плавное, всё более широкое и низкое тепловое пятно.

Характерное время выравнивания температуры

Величина $\tau = L^2/(a\pi^2)$ - это не просто вспомогательное число: она показывает, за какое время амплитуда основной (медленно затухающей) гармоники уменьшается в $e \approx 2{,}718$ раза, и удобна для сравнения материалов между собой. Для стержня длиной 20 см характерное время сильно различается:

• медь ($a \approx 1{,}162$ см²/с): $\tau \approx 35$ с;
• сталь ($a \approx 0{,}139$ см²/с): $\tau \approx 4{,}9$ мин;
• кирпич ($a \approx 0{,}0056$ см²/с): $\tau \approx 2{,}0$ ч;
• дерево ($a \approx 0{,}0018$ см²/с): $\tau \approx 6{,}3$ ч.

Разница на четыре порядка объясняет, почему металлическую ручку сковороды нельзя держать голой рукой уже через пару секунд, а деревянная ручка остаётся тёплой на ощупь ещё долго: у металлов коэффициент температуропроводности на порядки выше, чем у дерева и кирпича, поэтому тепло по ним «разбегается» несравнимо быстрее.

Стационарное распределение температуры и закон Фурье в задачах

Если ждать достаточно долго, распределение температуры перестаёт меняться во времени - наступает стационарный (установившийся) режим, и уравнение теплопроводности упрощается до $d^2T/dx^2 = 0$. Его решение - линейная функция координаты. Для стержня, соединяющего два тела с температурами $T_1$ и $T_2$ на концах, стационарный профиль:

$T(x) = T_1 + \frac{T_2 - T_1}{L}\,x,$

а поток тепла через любое сечение по закону Фурье одинаков и равен $q = k\,(T_1 - T_2)/L$. Это самая частая формулировка задач на закон Фурье: градиент постоянен, стержень греется или остывает равномерно вдоль всей длины, и достаточно одной формулы без всякого ряда Фурье.

Пример решения типовой задачи

Возьмём стальной стержень с начальным профилем в виде первой гармоники $T(x,0) = 100\sin(\pi x/L)$ и найдём, через сколько времени максимальная температура упадёт вдвое. Для стали коэффициент температуропроводности $a \approx 0{,}139$ см²/с. При длине стержня $L = 20$ см характерное время:

$\tau = \frac{L^2}{a\pi^2} = \frac{400}{0{,}139 \cdot \pi^2} \approx 292\ \text{с}.$

Поскольку в этом частном случае амплитуда убывает точно по экспоненте $e^{-t/\tau}$, момент, когда она уменьшится вдвое, находится из условия $e^{-t/\tau} = 0{,}5$:

$t = \tau\ln 2 \approx 292 \cdot 0{,}693 \approx 202\ \text{с} \approx 3\ \text{мин}\ 22\ \text{с}.$

Обратите внимание: ответ не зависит от начальной амплитуды 100 °C - удваивание или уменьшение начальной температуры не меняет время полураспада, потому что уравнение линейно, а экспонента затухания зависит только от $a$, $L$ и номера гармоники.

Частые ошибки

• Путают коэффициент теплопроводности $k$ и температуропроводности $a$. Это разные величины с разными единицами измерения: $k$ - Вт/(м·К), $a$ - м²/с. Скорость выравнивания температуры задаёт именно $a = k/(c\rho)$.
• Забывают перевести сантиметры в метры при расчёте потока по закону Фурье. Градиент, посчитанный в °C/см, нужно умножить на 100, чтобы получить °C/м перед подстановкой в формулу с $k$ в Вт/(м·К).
• Считают, что при $t=0$ сумма нескольких гармоник ряда Фурье точно равна разрывной ступеньке. Вблизи самого разрыва сумма конечного числа гармоник даёт характерный «выброс» (эффект Гиббса) - это не ошибка вычислений, а особенность разложения разрывной функции в ряд.
• Путают знак в законе Фурье. Поток направлен против градиента температуры, поэтому в формуле обязателен минус: $q = -k\,\partial T/\partial x$, а не $q = k\,\partial T/\partial x$.
• Применяют формулу стационарного линейного профиля к нестационарной задаче. Линейное распределение температуры справедливо только после того, как процесс установился; в начальный переходный период нужно решать полное уравнение с временной производной.

FAQ

Чем отличается уравнение теплопроводности от закона Фурье? Закон Фурье - это формула для потока тепла через сечение, $q = -k\,\partial T/\partial x$. Уравнение теплопроводности - дифференциальное уравнение для самой температуры, полученное из закона Фурье и баланса энергии; оно описывает, как температура меняется во времени и в пространстве.

Что такое коэффициент температуропроводности и чем он отличается от теплопроводности? Коэффициент теплопроводности $k$ показывает, насколько хорошо материал проводит тепло. Коэффициент температуропроводности $a = k/(c\rho)$ учитывает ещё и теплоёмкость с плотностью, поэтому именно он определяет скорость выравнивания температуры во времени.

Почему решение уравнения теплопроводности ищут в виде ряда Фурье? Потому что синусы автоматически удовлетворяют граничным условиям на закреплённых концах, а уравнение линейно - значит, сумма решений тоже решение. Раскладывая произвольное начальное условие в ряд по этим синусам, получают точное решение задачи.

Коротко

Закон Фурье $q = -k\,\partial T/\partial x$ связывает поток тепла с температурным градиентом, а из него и баланса энергии выводится уравнение теплопроводности $\partial T/\partial t = a\,\partial^2 T/\partial x^2$ с коэффициентом температуропроводности $a = k/(c\rho)$. Для стержня с закреплёнными на нуле концами решение раскладывается в ряд Фурье $T(x,t) = \sum b_n\sin(n\pi x/L)e^{-a(n\pi/L)^2 t}$, где высшие гармоники затухают быстрее и первыми сглаживают резкие перепады температуры. Характерное время выравнивания $\tau = L^2/(a\pi^2)$ различается на порядки между металлами и изоляторами, а в установившемся режиме распределение температуры становится линейным.