Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лекции математика / 04 Лекции 16 Теория поля Эктов

.pdf
Скачиваний:
75
Добавлен:
02.04.2015
Размер:
296.7 Кб
Скачать

1. Основные понятия теории поля

Теория поля лежит в основе многих представлений современной физики, механики, математики. Основными ее понятиями являются градиент, поток, потенциал, ротор, дивергенция, циркуляция и др. Эти понятия важны и в усвоении основных идей математического анализа функций многих переменных.

Полем называется область G пространства, в каждой точке которой определено значение некоторой величины.

В физических задачах обычно встречаются величины двух типов: скаляры и векторы. В соответствии с этим рассматривают два вида полей.

Если каждой точке M этой области поставлено в соответствие некоторое число U (M ) , говорят, что в

области задано (определено) скалярное поле. Примерами скалярных полей являются поле температур внутри некоторого нагретого тела (в каждой точке M этого тела задана соответствующая температура U (M ) ), поле

освещённости, создаваемое каким-либо источником света. Пусть в пространстве фиксирована система

координат. Тогда задание поля U (M )

равносильно

заданию одной функции U = U (x, y, z) ,

где x, y,z --

координаты точки M в этой системе координат. Значения функции U (x, y, z) совпадают со значениями поля U (M ) ,

поэтому для неё сохранён тот же символ.

Если каждой точке M этой области поставлен в соответствие определённый вектор a(M ) , говорят, что

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

некоторой скоростью v , не зависящей от времени (но

различной, вообще говоря, в разных точках); поставив в соответствие каждой точке M из G вектор v(M ) , мы получим векторное поле, называемое полем скоростей.

Если a(M ) - некоторое векторное поле в

пространстве, то взяв в этом пространстве фиксированную прямоугольную декартову систему координат, можно

представить a(M ) как упорядоченную тройку скалярных

функций: a(M ) = (P(x, y, z),Q(x, y, z), R(x, y, z)) . Эти

функции будем считать непрерывными и имеющими непрерывные частные производные первого и второго порядка. Кроме того, область G полагаем односвязной.

Если функция U (M ) (или a(M ) ) не зависит от

времени, то скалярное (векторное) поле называется стационарным; поле, которое меняется с течением времени, называется нестационарным. Далее будем рассматривать только стационарные поля.

2. Основные характеристики скалярного и векторного полей

Вектор, координатами которого являются значения частных производных функции U (x, y, z) в точке

M (x, y, z) , называют градиентом функции и обозначают

gradU (x, y, z) , т.е.

 

 

 

 

 

 

 

 

U (M )

,

U (M )

,

U (M )

 

(1)

gradU (x, y, z) =

 

 

x

 

y

 

z

 

 

Известно, что градиент задаёт в точке M направление наибыстрейшего возрастания функции U (x, y, z) . Говорят, что скалярное поле U порождает

векторное поле градиента U .

Линией градиента скалярного поля U (M ) называют

2

всякую кривую, касательная к которой в каждой точке направлена по gradU в этой точке.

Таким образом, линии градиента поля - это те линии, вдоль которых поле меняется быстрее всего.

Чтобы сформулировать ещё одно свойство градиента, напомним определение поверхности уровня.

Поверхностью уровня функции (поля) U = U (x, y, z)

называется поверхность, на которой функция (поле) сохраняет постоянное значение. Уравнение поверхности уровня имеет вид U (x, y, z) = C .

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

Потоком Π векторного поля a = (P,Q, R) через

поверхность σ называется поверхностный интеграл

∫∫(P cosα +Q cos β + R cosγ) dS

σ

или, короче, ∫∫a n dS, где через n = (cosα, cos β, cosγ)

σ

обозначили единичный вектор нормали к поверхности σ , определяющий её сторону.

Дивергенцией векторного поля a(M ) в

точке

M

 

lim

∫∫

a ns

dS

 

 

 

 

называется предел

σ

 

 

,

где

v

- объём

 

v

 

 

v0

 

 

 

 

 

 

области Ω G , содержащей

точку M , а σ

-

граница

области Ω, который обозначается diva(M ) .

 

 

 

Если частные

производные

P ,

Q ,

R

 

 

 

 

 

 

x

y

 

z

непрерывны, то

3

 

 

 

 

s

 

 

 

 

P +

Q +

R .

 

 

 

 

 

 

 

 

(2)

 

 

 

div a(M ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

z

 

 

 

 

 

 

 

 

 

Ротором (или вихрем) векторного поля a = (P,Q, R)

называется следующий вектор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

R

Q

s

P

R

s

 

Q

P

s

 

rot a

=

y

z

i +

z

x

j +

 

x

y

k.

 

 

 

 

 

 

 

 

 

 

 

 

 

Ротор векторного поля удобно записывать в виде

символического детермината

 

s

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot a =

 

x

 

 

 

y

 

z

,

 

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

P

 

 

 

Q

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

,

 

 

где под умножением одного из символов

 

на

 

x

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

некоторую

функцию

 

 

 

 

 

понимается

 

 

выполнение

соответствующей

 

операции

 

 

 

дифференцирования

(например,

Q означает

Q

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть L -- замкнутая кривая в области Ω. Интеграл

 

 

 

 

 

P dx +Q dy + R dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

называется циркуляцией поля a = (P, Q, R)

 

по кривой L и

обозначается

a d r,

 

где

 

d r = (dx, dy, dz) .

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Формулы Стокса и Остроградского-Гаусса

Обозначим через L -- некоторый замкнутый контур, а σ -- поверхность, натянутая на этот контур.

4

Предполагается, что выбор направления на контуре согласован с выбором стороны поверхности (при обходе контура в выбранном направлении выбранная сторона находится слева).

В этих условиях имеет место

формула Стокса

 

adrs

= ∫∫rot as

ns

dS.

(4)

L

σ

 

 

 

Формула Стокса говорит о том, что циркуляция векторного поля вдоль некоторого контура равна потоку ротора векторного поля через поверхность, натянутую на этот контур.

Пусть теперь Ω - некоторая замкнутая ограниченная область, а σ - граница этой области. Тогда справедлива

формула Остроградского-Гаусса

 

∫∫ a ns dS = ∫∫∫div as dx dy dz.

(5)

σΩ

Напомним, что поверхностный интеграл в левой части формулы (5) берется по внешней стороне поверхности σ .

Формула Остроградского-Гаусса означает, что тройной интеграл по области от дивергенции векторного поля равен потоку этого поля через поверхность, ограничивающую данную область.

4. Оператор Гамильтона. Некоторые виды скалярных и векторных полей

Английским математиком и механиком Гамильтоном был введён векторный дифференциальный оператор

=

i +

j +

k,

(6)

x

y

z

 

 

 

 

 

называемый оператором набла.

Следует сразу отметить, что аналогия между символическим вектором и "настоящими" векторами не

5

полная. Именно, формулы, содержащие символический вектор , аналогичны обычным формулам векторной алгебры в том случае, если они не содержат произведений переменных величин (скалярных и векторных), то есть до тех пор, пока не приходится применять входящие в операции дифференцирования к произведению переменных величин.

Используя набла-вектор, градиент скалярного поля

U (M ) , дивергенцию и ротор векторного поля

a(M )

можно записать в виде

 

gradU = sU ,

(7)

div a = as,

(8)

rot a = ×a.

(9)

Целесообразность введения символического векторасостоит в том, что с его помощью удобно получать и записывать различные формулы векторного анализа.

Продемонстрируем сказанное на примерах.

Задача 1. Доказать, что ротор градиента скалярного поля U (M ) равен 0 , то есть rot(gradU ) = 0.

Решение.

Докажем сначала это равенство, не используя оператор Гамильтона. Таким образом,

rot(gradU ) = rot

U (M )

, U (M ) ,

U (M ) =

(1)

 

 

(3)

 

x

y

z

6

i j

=x y U U

x y

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

U

U

 

U

U

 

 

 

 

 

 

 

 

 

U

=

zy

 

 

i +

xz

 

j +

 

 

 

yz

 

zx

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

yx

xy

k = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поскольку по теореме Шварца непрерывные смешанные производные равны.

Теперь, используя форму записи градиента (7) и ротора (9) через , имеем rot(gradU ) = × U.

Так как вектор U (произведение вектора на скаляр U ) коллинеарен вектору , то их векторное

произведение равно 0 .

Задача 2. Записать дивергенцию градиента скалярного поля div(gradU ) , используя .

Решение.

Образуя дивергенцию от gradU , получим

div(gradU ) = div

U si + U sj +

U ks =

(1)

 

 

(2)

x

y

z

=2U + 2U + 2U . x2 y2 z2

Оператор

2

+

2

+

2

называют оператором

x2

y2

z2

 

 

 

 

Лапласа и обозначают символом :

7

= 2 + 2 + 2 . x2 y2 z2

Поскольку скалярный квадрат вектора равен квадрату его модуля, то = 2 . Таким образом, div(gradU ) = 2U.

Векторное поле a(M ) называется потенциальным,

если его можно представить как градиент некоторого скалярного поля U (M ) :

a = gradU.

Само скалярное поле U называется при этом потенциалом векторного поля a .

Для того, чтобы векторное поле a(M ) было

потенциальным, необходимо и достаточно, чтобы

выполнялось условие

s

 

 

(10)

 

rot a = 0.

Необходимость выполнения равенства (10) доказана (см. задачу 1, рассмотренную выше).

Потенциал векторного поля можно найти по формуле

x

y

z

U (M ) = P(x, y, z)dx + Q(x0 , y, z)dy + R(x0 , y0 , z)dz +C,

x0

y0

z0

 

 

(11)

где (x0 , y0 , z0 ) - произвольная точка области G .

Векторное поле a(M ) , дивергенция которого

тождественно равна нулю, называется соленоидальным (трубчатым).

Для того чтобы сформулировать одно из важнейших свойств соленоидального поля, введём понятия векторной линии и векторной трубки.

Линия L , лежащая в G называется векторной

8

линией, если в каждой точке этой линии направление касательной к ней совпадает с направлением векторного поля в этой точке.

Известно, что векторная линия является интегральной кривой системы дифференциальных уравнений

dx

=

dy

=

dz

.

(12)

P

Q

 

 

 

R

 

В частности, если векторное поле есть поле скоростей стационарного потока жидкости, то его векторные линии -- это траектории частиц жидкости.

Векторная трубка - это замкнутое множество Φ точек области G , в котором задано векторное поле a(M ) , такое, что всюду на его граничной поверхности вектор нормали n ортогонален a(M ) .

Векторная трубка состоит из векторных линий поля a(M ) . Векторная линия целиком содержится в Φ, если

одна точка линии содержится в Φ.

Интенсивностью трубки Φ в сечении называется поток поля a(M ) через это сечение.

Если поле соленоидальное, то выполняется закон сохранения интенсивности векторной трубки.

Для поля скоростей v(M ) несжимаемой жидкости приsотсутствии стоков и источников (то есть при условии div v(M ) = 0 ) закон сохранения интенсивности векторной

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

Ниже приводятся некоторые типичные задачи с решениями.

9

Задача 3. Найти поверхности уровня скалярного поля

U (M ) = x2 + y2 z .

Решение.

 

 

 

Поверхности

уровня

определяются

уравнениями

x2 + y2 z = C или

x2 + y2

= z +C . Отсюда

видно, что

поверхности уровня представляют собой семейство эллиптических параболоидов, осью симметрии которых является ось Oz .

Задача 4.

В скалярном поле U (M ) = xy2 + z2 найти

градиент в точке M0 (2,1,1) .

 

 

 

 

 

 

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем значения

частных производных

функции

U (M ) в точке M0 :

 

 

 

 

 

 

 

 

 

 

 

 

U

|M0 = y2 |M0 = 12 = 1,

U

|M0 = 2xy |M0 = 2 2 1 = 4,

 

x

y

 

 

 

 

 

U

 

|

= 2z |

= 2 (1) = 2.

 

 

 

 

 

 

z

 

 

 

 

 

 

 

M0

 

M0

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gradU (M0 ) = si +4sj 2ks.

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

 

 

Задача 5.

Вычислить дивергенцию векторного поля

a(M ) = 2xy2 i yz j +3z2 k

в точке M0 (1,2,1) .

 

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем

P = 2xy2 , Q = yz , R = 3z2 . Найдём значение

соответствующих частных производных в точке M0 :

 

P |

= 2 y2 |

 

= 2 (2)2 = 8,

Q

|

= z |

= 1,

 

 

 

 

x

M0

 

 

M

0

 

 

 

y

M0

M0

10