Фиг. 3.9. Некоторая поверхность, ограниченная контуром Г.
Поверхность разделена на множество маленьких участков, каждый примерно в форме квадрата. Циркуляция по Г есть сумма циркуляции по всем маленьким контурам.
Какой бы ни была форма поверхности, но если малые контуры сделать достаточно малыми, всегда можно будет считать каждый из них замыкающим достаточно плоскую поверхность. Кроме того, каждый из них можно сделать очень похожим на квадрат. И циркуляцию вокруг большого контура Г можно найти, подсчитав циркуляции по всем квадратикам и сложив их.
§ 6. Циркуляция по квадрату; теорема Стокса
Как нам найти циркуляцию по каждому квадратику? Все зависит от того, как квадрат ориентирован в пространстве. Если ориентация его подобрана удачно (к примеру, он расположен в одной из координатных плоскостей), то расчет сделать легко. Так как пока мы не делали никаких предположений об ориентации осей координат, мы вправе выбрать их так, чтобы тот квадратик, на котором мы сосредоточили свое внимание, оказался в плоскости ху (фиг. 3.10). Если результат расчета будет выражен в векторной записи, то можно говорить, что он не зависит от специальной ориентации плоскости.
Фиг. 3.10. Вычисление циркуляции вектора С по маленькому квадратику.
Мы хотим теперь найти циркуляцию поля С по нашему квадратику. Криволинейное интегрирование легко проделать, если квадратик сделать таким маленьким, чтобы вектор С на протяжении одной стороны квадрата менялся очень мало. (Это предположение выполняется тем лучше, чем меньше квадратик, так что на самом деле речь идет о бесконечно малых квадратиках.) Отправившись от точки (х, у) — в левом нижнем углу фигуры,— мы обойдем весь квадрат в направлении, указанном стрелками. Вдоль первой стороны, отмеченной цифрой 1, касательная составляющая равна Сх(1), а расстояние равно Dх. Первая часть интеграла равна Cx(1) Dх, Вдоль второй стороны получится Су(2) Dy. Вдоль третьей мы получим -Сx(3) Dх, а вдоль четвертой -Cy(4) Dy. Знаки минус стоят потому, что нас интересует касательная составляющая в направлении обхода. Весь криволинейный интеграл тогда равен
(3.31) Посмотрим теперь на первый и третий члены. В сумме они дают
(3.32)
Вам может показаться, что в принятом приближении эта разность равна нулю. Но это только в первом приближении. Мы можем быть более точными и учесть скорость изменения Сх , тогда можно написать
(3.33)
В следующем приближении пойдут члены с (Dy)2, но ввиду того, что нас интересует в конечном счете только предел при Dy®0, то этими членами можно пренебречь. Подставляя (3.33) в (3.32), мы получаем
(3.34)
Производную при нашей точности можно брать в точке (х, у). Подобным же образом оставшиеся два члена можно написать в виде
(3.35)
и циркуляция по квадрату тогда равна
(3.36)
Интересно, что в скобках получилась как раз z-компонента ротора С. Множитель DxDy— это площадь нашего квадрата. Так что циркуляцию (3.36) можно записать как
(СXС)zDа.
Но z-компонента это на самом деле компонента, нормальная к элементу поверхности.
Фиг. 3.11. Циркуляция вектора С по Г равна поверхностному интегралу от нормальной компоненты вектора СXС.
Поэтому циркуляцию вокруг квадратика можно задать и в инвариантной векторной записи:
(3.37)
В результате имеем: циркуляция произвольного вектора С по бесконечно малому квадрату равна произведению составляющей ротора С, нормальной к поверхности, на площадь квадрата.
Циркуляция по произвольному контуру Г легко теперь может быть увязана с ротором векторного поля. Натянем на контур любую подходящую поверхность S (как на фиг. 3.11) и сложим между собой циркуляции по всем бесконечно малым квадратикам на этой поверхности. Сумма может быть записана в виде интеграла. В итоге получится очень полезная теорема, называемая теоремой Стокса [по имени физика Стокса].
ТЕОРЕМА СТОКСА
(3.38)
где S — произвольная поверхность, ограниченная контуром Г. Теперь мы должны ввести соглашение о знаках. На приведенной ранее фиг. 3.10 ось z показывает на вас, если система координат «обычная», т. е. «правая». Когда в криволинейном интеграле мы брали «положительное» направление обхода, то циркуляция получилась равной z-компоненте вектора СXC. Обойди мы контур в другую сторону, мы бы получили противоположный знак. Как вообще узнавать, какое направление надо выбирать для положительного направления «нормальной» компоненты вектора СXC? «Положительную» нормаль надо всегда связывать с направлением так, как это сделано было на фиг. 3.10. Общий случай показан на фиг. 3.11.