В случае трех измерений, контурные линии становятся поверхностями. Малое смещение Dr (= iDx + jDy + kDz) приводит к изменению f, которое записывается как
где точками обозначены члены более высоких порядков. Это выражение можно записать в виде скалярного произведения
Разделим правую и левую части этого равенства на Ds, и пусть Ds стремится к нулю; тогда
где dr/ds – единичный вектор в выбранном направлении. Выражение в круглых скобках – вектор, зависящий от выбранной точки. Таким образом, df/ds имеет максимальное значение, когда dr/ds указывает в том же направлении, выражение, стоящее в скобках, является градиентом. Таким образом,
– вектор, равный по величине и совпадающий по направлению с максимальной скоростью изменения f относительно координат. Градиент f часто записывается в виде
Это означает, что оператор С существует сам по себе. Во многих случаях он ведет себя, как вектор, и фактически является «векторным дифференциальным оператором» – одним из наиболее важных дифференциальных операторов в физике. Несмотря на то, что С содержит единичные векторы i, j и k, его физический смысл не зависит от выбранной системы координат.
Какова связь между Сf и f? Прежде всего предположим, что f определяет потенциал в любой точке. При любом малом смещении Dr величина f изменится на
Если q – величина (например масса, заряд), перемещенная на Dr, то работа, выполненная при перемещении q на Dr равна
Так как Dr – перемещение, то qСf – сила; –Сf – напряженность (сила на единицу количества), связанная с f. Например, пусть U – электростатический потенциал; тогда E – напряженность электрического поля, задается формулой
E = –СU.
Допустим, что U создается точечным электрическим зарядом в q кулонов, помещенным в начало координат. Значение U в точке P (x, y, z) с радиус-вектором r задается формулой
где e0 – диэлектрическая постоянная свободного пространства. Поэтому
откуда следует, что E действует в направлении r и его величина равна q/(4pe0r3).
Зная скалярное поле, можно определить связанное с ним векторное поле. Также возможно и обратное. С точки зрения математической обработки скалярными полями оперировать легче, чем векторными, так как они задаются одной функцией координат, в то время как векторное поле требует три функции, соответствующие компонентам вектора в трех направлениях. Таким образом, возникает вопрос: дано векторное поле, может ли мы записать связанное с ним скалярное поле?
Дивергенция и ротор. Мы видели результат действия С на скалярную функцию. Что произойдет, если С применить к вектору? Имеются две возможности: пусть U (x, y, z) – вектор; тогда мы можем образовать векторное и скалярное произведения следующим образом:
Первое из этих выражений – скаляр, называемый дивергенцией U (обозначается divU); второе – вектор, названный ротор U (обозначается rotU).
Эти дифференциальные функции, дивергенция и ротор, широко используются в математической физике.
Представьте, что U – некоторый вектор и что он и его первые производные непрерывны в некоторой области. Пусть P – точка в этой области, окруженная малой замкнутой поверхностью S, ограничивающей объем DV. Пусть n – единичный вектор, перпендикулярный к этой поверхности в каждой точке (n меняет направление при движении вокруг поверхности, но всегда имеет единичную длину); пусть n направлен наружу. Покажем, что
и
Здесь S указывает, что эти интегралы берутся по всей поверхности, da – элемент поверхности S.
Для простоты мы выберем удобную для нас форму S в виде небольшого параллелепипеда (как показано на рис. 12) со сторонами Dx, Dy и Dz; точка P – центр параллелепипеда. Вычислим интеграл из уравнения (4) сначала по одной грани параллелепипеда. Для передней грани n = i (единичный вектор параллелен оси x); Da = DyDz. Вклад в интеграл от передней грани равен
Рис. 12. ОБЪЕМ ДЛЯ ВЫЧИСЛЕНИЯ ДИВЕРГЕНЦИИ. Стороны параллелепипеда Dx, Dy и Dz; P – центр; n – единичные векторы, перпендикулярные к поверхностям.
На противоположной грани n = –i; эта грань дает вклад в интеграл
Используя теорему Тейлора, получим общий вклад от двух граней
Заметим, что DxDyDz = DV. Аналогичным образом можно вычислить вклад от двух других пар граней. Полный интеграл равен
и если мы положим DV ® 0, то члены более высокого порядка исчезнут. По формуле (2) выражение в скобках – это divU, что доказывает равенство (4).
Равенство (5) можно доказать таким же образом. Воспользуемся снова рис. 12; тогда вклад от передней грани в интеграл будет равен
и, используя теорему Тейлора, получим, что суммарный вклад в интеграл от двух граней имеет вид
т.е. это два члена из выражения для rotU в уравнении (3). Другие четыре члена получатся после учета вкладов от других четырех граней.
Что, в сущности, означают эти соотношения? Рассмотрим равенство (4). Предположим, что U – скорость (жидкости, например). Тогда nЧU da = Un da, где Un является нормальной компонентой вектора U к поверхности. Поэтому, Un da – это объем жидкости, протекающей через da в единицу времени, а – это объем жидкости, вытекающей через S в единицу времени. Следовательно,
– скорость расширения единицы объема вокруг точки P. Отсюда дивергенция получила свое название; она показывает скорость, с которой жидкость расширяется из (т.е. расходится от) P.
Чтобы объяснить физическое значение ротора U, рассмотрим другой поверхностный интеграл по маленькому цилиндрическому объему высотой h, окружающему точку P; плоско-параллельные поверхности могут быть ориентированы в любом направлении, которое мы выбираем. Пусть k –единичный вектор перпендикулярный к каждой поверхности, и пусть площадь каждой поверхности DA; тогда полный объем DV = hDA (рис. 13). Рассмотрим теперь интеграл
Подынтегральное выражение – уже упоминавшееся ранее тройное скалярное произведение. Это произведение будет равно нулю на плоских поверхностях, где k и n параллельны. На кривой поверхности
где ds – элемент кривой как показано на рис. 13. Сравнивая эти равенства с соотношением (5), получаем, что
Мы по-прежнему предполагаем, что U – скорость. Чему в таком случае будет равна средняя угловая скорость жидкости вокруг k? Очевидно, что
если DA ® 0. Это выражение максимально, когда k и rotU указывают в одном и том же направлении; это означает, что rotU – вектор, равный удвоенной угловой скорости жидкости в точке P. Если жидкость вращается относительно P, то rotU № 0, и векторы U будут вращаться вокруг P. Отсюда и возникло название ротора.
Теорема дивергенции (теорема Остроградского – Гаусса) является обобщением формулы (4) для конечных объемов. Она утверждает, что для некоторого объема V, ограниченного замкнутой поверхностью S,
и справедлива для всех непрерывных векторных функций U, имеющих непрерывные первые производные всюду в V и на S. Мы не будем приводить здесь доказательство этой теоремы, но ее справедливость можно понять интуитивно, представляя объем V разделенным на ячейки. Поток U через поверхность, общую для двух ячеек обращается в нуль, и только ячейки, находящиеся на границе S внесут вклад в поверхностный интеграл.
Теорема Стокса является обобщением уравнения (6) для конечных поверхностей. Она утверждает, что
где C – замкнутая кривая и S – любая поверхность, ограниченная этой кривой. U и ее первые производные должны быть непрерывны всюду на S и C.