Дифференциальное уравнение распространения тепла было получено в гл. 2. В соответствии с уравнением (2.44),
(3.25)
(Помните, что это соотношение приближенное, но для некоторых веществ вроде металлов выдерживается неплохо.) Применимо оно, конечно, только в тех частях тела, где нет ни выделения, ни поглощения тепла. Выше мы вывели другое соотношение (3.21), которое выполняется тогда, когда количество тепла сохраняется. Если мы это уравнение скомбинируем с (3.25), то получим
или
(3.26)
если c — величина постоянная. Напоминаю, что q — это количество тепла в единичном объеме, а С·С = С2 — лапласиан, т. е. оператор
Если мы теперь сделаем еще одно допущение, сразу возникнет одно очень интересное уравнение. Допустим, что температура материала пропорциональна содержанию тепла в единице объема, т. е. что у материала есть определенная удельная теплоемкость. Когда это допущение верно (а так бывает часто), мы можем писать
или
(3.27)
Скорость изменения количества тепла пропорциональна скорости изменения температуры. Коэффициент пропорциональности cvздесь — удельная теплоемкость на единицу объема материала. Подставляя (3.27) в (3.26), получаем
(3.28)
Мы обнаружили, что быстрота изменения со временем температуры Т в каждой точке пропорциональна лапласиану от Т, т. е. вторым производным от пространственного распределения температур. Мы имеем дифференциальное уравнение — в переменных х, у, z и t — для температуры Т.
Дифференциальное уравнение (3.28) называется уравнением диффузии тепла, или уравнением теплопроводности. Часто его пишут в виде
(3.29)
где D — постоянная. Она равна x/cv.
Уравнение диффузии появляется во многих физических задачах: о диффузии газов, диффузии нейтронов и других. Мы уже обсуждали физику некоторых таких явлений в вып. 4, гл. 43. Теперь перед вами полное уравнение, описывающее диффузию в самом общем виде. Немного позже мы займемся решением уравнения диффузии, чтобы посмотреть, как распределяется температура в некоторых случаях. А сейчас вернемся к рассмотрению других теорем о векторных полях.
§ 5. Циркуляция векторного поля
Мы хотим теперь рассмотреть ротор поля примерно так же, как рассматривали дивергенцию. Мы вывели теорему Гаусса, вычисляя интеграл по поверхности, хотя с самого начала отнюдь не было ясно, что мы будем иметь дело с дивергенцией. Откуда же можно было знать, что для ее получения надо интегрировать по поверхности? Этот результат вовсе не был очевиден. И столь же неоправданно мы сейчас вычислим другую характеристику поля и покажем, что она связана с ротором. На этот раз мы подсчитаем так называемую циркуляцию векторного поля. Если С — произвольное векторное поле, мы возьмем его составляющую вдоль кривой линии и проинтегрируем эту составляющую по замкнутому контуру. Интеграл называется циркуляцией векторного поля по контуру. Мы уже раньше в этой главе рассматривали криволинейный интеграл от Сy. Сейчас мы то же самое проделываем с произвольным векторным полем С.
Пусть Г — произвольный замкнутый контур в пространстве (воображаемый, разумеется). Пример мы видим на фиг. 3.7. Криволинейный интеграл от касательной составляющей С по контуру записывается в виде
(3.30)
Фиг. 3.7. Циркуляция вектора С но кривой Г есть криволинейный интеграл от С t (касательной составляющей С).
Фиг. 3.8. Циркуляция по всему контуру есть сумма циркуляции по двум контурам: Г 1 =Г a +Г ab и Г 2 =Г ь +Г a Ь .
Заметьте, что интеграл берется по всему замкнутому пути, а не от одной точки до другой, как это делалось раньше. Кружочек на знаке интеграла должен нам напоминать об этом. Такой интеграл называется циркуляцией векторного поля по кривой Г. Название связано с тем, что первоначально так рассчитывали циркуляцию жидкости. Но название это, как и поток, было распространено на любые поля, даже такие, в которых «циркулировать» нечему.
Забавляясь той же игрой, как с потоком, мы можем показать, что циркуляция вдоль контура есть сумма циркуляции вдоль двух меньших контуров. Положим, что, соединив две точки (1) и (2) первоначальной кривой с помощью некоторой линии, мы разбили кривую на два контура Г1 и Г2 (фиг. 3.8). Контур Г1 состоит из Гa — части первоначальной кривой слева от (1) и (2) и «соединения» Гab. Контур Г2 состоит из остатка первоначальной кривой плюс то же соединение.
Циркуляция вдоль Г1 есть сумма интеграла вдоль Га и вдоль ГаЬ. Точно так же и циркуляция вдоль Г2 есть сумма двух частей, одной вдоль Гb, другой — вдоль Гab. Интеграл вдоль Гab для кривой Г2 имеет знак, противоположный тому знаку, который он имел для кривой Г1, потому что направления обхода противоположны (в обоих криволинейных интегралах направления поворота нужно брать одни и те же).
Повторяя прежние аргументы, мы можем убедиться, что сумма двух циркуляции даст как раз криволинейный интеграл вдоль первоначальной кривой Г. Интегралы по Гab сократятся. Циркуляция по одной части плюс циркуляция вдоль другой равняется циркуляции вдоль внешней линии. Этот процесс разрезания большого контура на меньшие можно продолжить. При сложении циркуляции по меньшим контурам смежные части будут сокращаться, так что сумма их сведется к циркуляции вдоль единственного первоначального контура.
Теперь предположим, что первоначальный контур — это граница некоторой поверхности. Существует бесконечное множество поверхностей, границей которых служит все тот же первоначальный замкнутый контур. Наши результаты не зависят, однако, от выбора этих поверхностей. Сперва мы разобьем наш первоначальный контур на множество малых контуров, лежащих на выбранной поверхности (фиг. 3.9).
Фиг. 3.9. Некоторая поверхность, ограниченная контуром Г.
Поверхность разделена на множество маленьких участков, каждый примерно в форме квадрата. Циркуляция по Г есть сумма циркуляции по всем маленьким контурам.
Какой бы ни была форма поверхности, но если малые контуры сделать достаточно малыми, всегда можно будет считать каждый из них замыкающим достаточно плоскую поверхность. Кроме того, каждый из них можно сделать очень похожим на квадрат. И циркуляцию вокруг большого контура Г можно найти, подсчитав циркуляции по всем квадратикам и сложив их.