Теорема Гаусса. Дивергенция

Поверхностный интеграл можно преобразовать в объёмный; в этом заключается содержание одной из важнейших теорем векторного анализа—теоремы Гаусса.

Рассмотрим сначала поток dN произвольного, но дифференцируемого вектора а через поверхность бесконечно малого параллелепипеда и выберем для удобства вычислений направление осей координат ординат х, у, z так, чтобы они совпадали с ребрами этого параллелепипеда dх, dу, dz (рис. 1 5).

Рис. 1.5. Бесконечно малый параллепипед.

Интеграл сводится в этом случае к сумме шести интегралов по каждой из граней параллелепипеда. Воспользовавшись известной из интегрального исчисления теоремой о среднем, можно каждый из этих шести интегралов представить как произведение площади грани на некоторое среднее значение нормальной слагающей вектора а на данной грани.

Рассмотрим сначала поток вектора а через две параллельные грани 1 и 2, перпендикулярные к оси х. Поток через переднюю грань 2 равен:

где у и z — некоторые средние значения координат y и z на грани 2 и а ₂.—значение вектора а на грани 2; поток через заднюю грань 1 равен:

где a₁ значение вектора а на грани 1, ибо внешняя нормаль к этой грани направлена прямо противоположно оси х. Стало быть, общий поток через грани 1 и 2 равен:

(а _2x- а _{1 х} ) dуdz.

Разность а _{2 х} — а_1х есть приращение слагающей вектора а_х при изменении координаты х на расстояние dх между гранями 1 и 2. С точностью до бесконечно малых второго порядка приращение это равно:

a_2x (х+ dх, у, z) — а_1х (х, у, z) = dх,

где ввиду бесконечной малости параллелепипеда под можно понимать значение этой производной в любой точке параллелепипеда. Таким образом, общий поток через обе грани, перпендикулярные к оси х, равен:

dх dу dz.

Для потоков через пары граней, перпендикулярных к осям y и z получим аналогично:

dх dу dz и dх dу dz

Складывая полученные выражения, получим общий поток вектора а через все шесть граней элементарного параллелепипеда:

d = (1.10)

Стоящую в скобках сумму производных вектора а по осям координат принято для краткости обозначать символом div а:

(1.11)

Если, кроме того, ввести для бесконечно малого элемента объема обозначение dV:

dV=dxdydz

то выражение потока dN примет вид

dN =div a dV. (1.12)

Эту формулу, выражающую поток вектора а через поверхность бесконечно малого параллелепипеда, нетрудно обобщить для поверхности произвольной формы и размеров. Рассмотрим произвольную замкнутую поверхность S. Разобьем ограниченный ею объем V системой взаимно перпендикулярных плоскостей на совокупность бесконечно малых кубических элементов. Конечно, крайние_; смежные с поверхностью S, элементы объема, вообще говоря, не будут иметь кубической формы; однако путем дальнейшего дробления их можно достигнуть того, чтобы грани крайних кубиков с любой степенью точности совпадали с заданной поверхностью S. Вычислим с помощью уравнения (16*) поток вектора а через поверхность каждого кубика, лежащего внутри S, и сложим полученные выражения:

В этом уравнении тройной интеграл означает, что суммирование подынтегрального выражения должно быть произведено по всем элементам трехмерного объема V, заключенного внутри поверхности S. Однако на протяжении всей этой книги мы обозначали интегралы любой кратности одним-единственным знаком ; различение же интегралов разной кратности достигалось различным обозначением элементов интегрирования: элемент объема (трехкратного интеграла) обозначался через dV, элемент поверхности (двукратного интеграла) через dS, элемент линии (одинарного интеграла) через dl.

Грани всех элементарных кубиков, составляющих в совокупности объем V, могут быть разделены на два класса — грани внешние, совпадающие с элементами поверхности S, и грани внутренние, отграничивающие смежные кубики друг от друга. Очевидно, что в сумму поток вектора а через каждую внутреннюю грань войдет дважды: при подсчете потока через поверхность кубика, лежащего по одну сторону от этой грани, и при подсчете потока через поверхность кубика, лежащего по другую сторону от нее. Так как нормаль к грани, внешняя по отношению к первому кубику, противоположна нормали к той же грани, внешней по отношению ко второму кубику, то оба потока через эту грань будут иметь противоположные знаки. Следовательно, все члены суммы , относящиеся ко внутренним граням, сократятся, и сумма эта сведется к сумме потоков вектора а через одни лишь внешние грани кубиков, совпадающие с элементами поверхности S. Таким образом, оказывается равной потоку N вектора а через заданную поверхность S, и, стало быть,

N= (1.13)

Это выражение. представляет собой теорему Гаусса: поток вектора а, являющегося непрерывной функцией точки, через произвольную замкнутую поверхность S равен интегралу дивергенции этого вектора по объему V, ограниченному этой поверхностью.

Если поверхность S столь мала, что во всех лежащих внутри нее точках div а можно считать величиной постоянной, то в уравнении (1.13) div а можно вынести за знак интеграла. Стало быть, поток dN через бесконечно малую замкнутую поверхность S произвольной формы выражается той же формулой (1.12):

dN=

как и поток через поверхность элементарного параллелепипеда. Так как эта формула справедлива лишь в предельном случае бесконечно малой поверхности, то ее правильнее записать в следующей форме:

. (1.14)

Правильнее всего считать эту формулу определением понятия дивергенции: дивергенция вектора а в данной точке поля есть предел, к которому стремится отношение потока вектора а через произвольную, окружающую эту точку, поверхность к величине ограниченного этой поверхностью объема (при 0). Из этого определения дивергенции следует, что значение ее вовсе не зависит от выбора системы координат, т. е. что дивергенция вектора есть истинный скаляр. Исходя из (1.14) и воспользовавшись (1.12), мы в частном случае декартовых координат, очевидно, вновь придем к (1.10).

Отметим в заключение, что в гидродинамике дивергенция скорости жидкости v имеет непосредственное физическое значение. Действительно, в каждой точке жидкости

равна рассчитанному на единицу объема количеству жидкости, вытекающей из элемента объема dV, окружающего рассматриваемую точку. Название «дивергенция», что значит по-латыни расхождение или расходимость, было избрано для этой величины именно потому, что жидкость растекается или расходится из тех и только из тех точек или участков занимаемого ею пространства, в которых div v >0. Очевидно, что в этих точках должны быть расположены источники жидкости. По аналогии, те точки поля произвольного вектора а, в которых diva 0, принято называть истоками этого поля. Численная же величина div a называется силой, или обильностью истоков поля; в зависимости от знака дивергенции сила истоков может быть как положительной, так и отрицательной. Иногда отрицательным истокам поля дают название стоков поля. Векторные поля, у которых div a =0, называются свободными от источников, или соленоидальными, или вихревыми.