Ранее мы уже обсуждали (гл. 3, § 4) поток тепла. Вообразите кусок какого-то материала, необязательно однородного (в разных местах может быть разное вещество), в котором температура меняется от точки к точке. Как следствие этих температурных изменений возникает поток тепла, который можно обозначить вектором h. Он представляет собой количество тепловой энергии, которое проходит в единицу времени через единичную площадку, перпендикулярную потоку. Дивергенция h есть скорость ухода тепла из данного места в расчете на единицу объема:
Ñ•h = Скорость ухода тепла на единицу объема.
(Мы могли, конечно, записать уравнение в интегральном виде, как мы поступали в электродинамике с законом Гаусса, тогда оно выражало бы тот факт, что поток через поверхность равен скорости изменения тепловой энергии внутри материала. Мы не будем больше переводить уравнения из дифференциальной формы в интегральную и обратно, это делается точно так же, как в электростатике.)
Скорость, с которой тепло поглощается или рождается в разных местах, конечно, зависит от условий задачи. Предположим, например, что источник тепла находится внутри материала (возможно, радиоактивный источник или сопротивление, через которое пропускают ток). Обозначим через s тепловую энергию, производимую этим источником в единице объема за 1 сек. Кроме того, могут возникнуть потери (или, наоборот, дополнительное рождение) тепловой энергии за счет перехода в другие виды внутренней энергии в данном объеме. Если и — внутренняя энергия в единице объема, то —du/dt будет тоже играть роль «источника» тепловой энергии. Итак, имеем
(12.5)
Мы не собираемся здесь обсуждать полное уравнение, величины в котором изменяются со временем, потому что мы проводим аналогию с электростатикой, где ничто не зависит от времени. Мы рассмотрим только задачи с постоянным потоком тепла, в которых постоянные источники создают состояние равновесия. В таких случаях
(12.6)
Нужно иметь, конечно, еще одно уравнение, которое описывает, как поток течет в разных местах. Во многих веществах поток тепла примерно пропорционален скорости изменения температуры с положением: чем больше разность температур, тем больше поток тепла. Мы знаем, что вектор потока тепла пропорционален градиенту температуры. Константа пропорциональности К, зависящая от свойств материала, называется коэффициентом теплопроводности
(12.7)
Если свойства материала меняются от точки к точке, то К=К (х, у, z) и есть функция положения. [Уравнение (12.7) не столь фундаментально, как (12.5), выражающее сохранение тепловой энергии, потому что оно зависит от характерных свойств вещества.] Подставляя теперь уравнение (12.7) в (12.6), получаем
(12.8)
что в точности совпадает по форме с (12.4). Задачи с постоянным потоком тепла и задачи электростатики одинаковы. Вектор потока тепла h соответствует Е, а температура Т соответствует j.
Фиг. 12.1. Поток тепла в случае цилиндрической симметрии (а) и соответствующая задача из электричества (б).
Мы уже отмечали, что точечный тепловой источник создает поле температур, меняющееся, как 1/r, и поток тепла, меняющийся, как 1/r2. Это есть не более чем простой перенос утверждений электростатики, что точечный заряд дает потенциал, меняющийся, как 1/r, и электрическое поле, меняющееся, как 1/r2. Вообще мы можем решать статистические тепловые задачи с той же степенью легкости, как и задачи электростатики.
Рассмотрим простой пример. Пусть имеется цилиндр с радиусом а при температуре T1? поддерживающейся за счет генерации тепла в цилиндре. (Это может быть, скажем, проволока, по которой течет ток, или трубка с конденсацией пара внутри цилиндра.) Цилиндр покрыт концентрической обшивкой из изолирующего материала с теплопроводностью К. Пусть внешний радиус изоляции равен b, а в наружном пространстве поддерживается температура T2 (фиг. 12. 1, а). Нам нужно определить скорость потери тепла проволокой или паропроводом (все равно чем), проходящим по центру цилиндра. Пусть полное количество тепла, теряемого на длине трубы L, равно G, его-то мы и хотим найти.
Как надо решать такую задачу? У нас есть дифференциальные уравнения, но поскольку они такие же, как в электростатике, то математическое решение их нам уже известно. Аналогичная задача электростатики относится к проводнику радиусом а при потенциале j1, отделенном от другого проводника радиусом b при потенциале j2, с концентрическим слоем диэлектрика между ними (фиг. 12.1, б). Далее, поскольку поток тепла h соответствует электрическому полю Е, то наша искомая величина G соответствует потоку электрического поля от единичной длины (другими словами, электрическому заряду на единице длины, деленному на e0). Мы решали электростатическую задачу с помощью закона Гаусса. Нашу задачу о потоке тепла будем решать таким же способом.
Из симметрии задачи мы видим, что h зависит только от расстояния до центра. Поэтому мы окружим трубку гауссовой поверхностью — цилиндром длиной L и радиусом r. С помощью закона Гаусса мы выводим, что поток тепла h, умноженный на площадь поверхности 2prL, должен быть равен полному количеству тепла, рождаемому внутри, т. е. тому, что мы назвали G:
(12.9)
Поток тепла пропорционален градиенту температуры
или в данном случае величина h равна
Вместе с (12.9) это дает
(12.10)
Интегрируя от r=а до r=b, получаем
(12.11)
Разрешая отнсительно G, находим
(12.12)
Этот результат в точности соответствует формуле для заряда цилиндрического конденсатора:
Задачи одинаковые и имеют одинаковые решения. Зная электростатику, мы тем самым знаем, сколько тепла теряет изолированная труба.
Рассмотрим еще один пример. Пусть мы хотим узнать поток тепла в окрестности точечного источника, расположенного неглубоко под поверхностью земли или же вблизи поверхности большого металлического предмета. В качестве локализованного источника тепла может быть и атомная бомба, которая взорвалась под землей и представляет собой мощный источник тепла, или же небольшой источник радиоактивности внутри железного блока — возможностей очень много.
Рассмотрим идеализированную задачу о точечном источнике тепла, мощность которого G, на расстоянии а под поверхностью бесконечной однородной среды с коэффициентом теплопроводности К. Теплопроводностью воздуха над поверхностью среды мы пренебрежем. Мы хотим определить распределение температуры на поверхности среды. Насколько горячо будет прямо над источником и в разных местах на поверхности?
Как же решить эту задачу? Она похожа на задачу по электростатике, в которой имеются два материала с разной диэлектрической проницаемостью x по обе стороны от разделяющей их границы. Здесь что-то есть! Возможно, это похоже на точечный заряд вблизи границы между диэлектриком и проводником или что-нибудь вроде этого. Посмотрим, что происходит вблизи границы. Физическое условие состоит в том, что нормальная составляющая h на поверхности равна нулю, поскольку мы предположили, что потока из блока нет. Мы должны задать вопрос: в какой электростатической задаче возникает условие, что нормальная компонента электрического поля Е (представляющая собой аналог h) равна нулю у поверхности? Нет такой!
Это один из тех случаев, к которым следует относиться с осторожностью. По физическим причинам могут быть определенные ограничения тех математических условий, которые возникают в каком-либо случае. Поэтому если мы проанализировали дифференциальное уравнение только для некоторых ограниченных примеров, то вполне можем упустить ряд решений, возникающих в других физических условиях. Например, нет материала, обладающего диэлектрической проницаемостью, равной нулю, а теплопроводность вакуума равна нулю. Поэтому нет электростатического аналога идеального теплоизолятора. Мы можем, однако, попытаться использовать те же методы. Попробуем вообразить, что произошло бы, если бы диэлектрическая проницаемость была равна нулю. (Разумеется, в реальных условиях диэлектрическая проницаемость никогда не обращается в нуль. Но может представиться случай, когда вещество имеет очень большую диэлектрическую проницаемость, так что диэлектрической проницаемостью воздуха вне среды можно пренебречь.)
Как же найти электрическое поле, у которого нет составляющей, перпендикулярной к поверхности? Иначе говоря, такое поле, которое всюду касательно к поверхности? Вы заметите, что эта задача обратна задаче о точечном заряде вблизи проводящей плоскости. Там нам нужно было поле, перпендикулярное
к поверхности, потому что проводник всюду находился при одном и том же значении потенциала.
В задаче об электрическом поле мы придумали решение, вообразив за проводящей плоскостью точечный заряд. Можно воспользоваться снова этой же идеей. Попытаемся выбрать такое «изображение» источника, которое автоматически обращало бы в нуль нормальную компоненту поля вблизи поверхности. Решение показано на фиг. 12.2. Электрическое изображение источника с тем же знаком и той же величины, находящееся на расстоянии а над поверхностью, дает поле, горизонтальное повсюду у поверхности. Нормальные компоненты от обоих источников взаимно уничтожаются.
Итак, наша задача о потоке тепла решена. Температура во всем пространстве одинакова по непосредственной аналогии с потенциалом от двух одинаковых точечных зарядов. Температура Т на расстоянии r от одного точечного источника G в бесконечной среде равна
(12.13)
(Это, конечно, полностью аналогично j= q/4pe0r.) Температура точечного источника и, кроме того, его изображения равна
(12.14)
Эта формула дает нам температуру всюду внутри блока. Несколько изотермических поверхностей приведено на фиг. 12.2.
Показаны также линии h, которые можно получить из выражения h =-КÑТ.
В самом начале мы интересовались распределением температуры на поверхности. Для точки на поверхности находящейся на расстоянии р от оси, r1=r2=Ö (р2 + а2),
Фиг. 12.2. Поток тепла и изотерма у точечного источника тепла, расположенного на расстоянии а под поверхностью тела с хорошей теплопроводностью. Вне тела показано мнимое изображение источника.
следовательно,
(12.15)
Эта функция также изображена на фиг. 12.2. Естественно, что температура прямо над источником выше, чем вдали от него. Такого рода задачи часто приходится решать геофизикам. Теперь мы видим, что это те же самые задачи, которые мы решали в электричестве.