Лапласиан

Ротор

По определению проекция (rot a) n ротора векторного поля a в некоторой точке P на направление фиксированного вектора n представляется формулой:

.

Следовательно, мы можем получить проекцию rot a на направление вектора e 1, вычислить циркуляцию вектора a вдоль контура MM 2 N 1 M 3 M (см. рис. 9.1) и разделив ее на d s1. Представим эту циркуляцию в виде суммы четырех слагаемых, отвечающих отрезкам MM 2, M 2 N 1, N 1 M 3 и M 3 M, и вычислим каждое слагаемое отдельно. Начнем с первого из них. Проекция вектора a на направление равна a 2, следовательно, циркуляция вектора a вдоль (с точностью до бесконечно малых величин первого порядка относительно d s1) равна

. (9.14)

Циркуляция вдоль N 1 M 3 отличаются от только что полученного выражения тем, что на N 1 M 3 третья координата равна q 3+ dq 3, а не q 3, как на MM 2, и, кроме того, направление отрезка N 1 M 3 противоположно направлению e 2. Поэтому циркуляция вдоль N 1 M 3 равна

. (9.15)

Аналогично получаем для циркуляции вдоль M 2 N 1 и N 3 M выражения

(9.16)

и

. (9.17)

Сложив величины (9.14)-(9.17), получим, что циркуляция вектора a вдоль контура MM 2 N 1 N 3 M равна

.

Деля полученное выражение на H 2 H 3 dq 2 dq 3, т.е. на площадь грани MM 2 N 1 M 3, получим, что компонента (rot a)1 вектора rot a в направлении базисного вектора e 1 равна

. (9.18)

Аналогично вычисляются две другие компоненты:

, (9.19)

. (9.20)

Символически это можно записать в виде

(9.21)

Оператор Лапласа DU функции U определяется формулой

.

Пользуясь формулами для градиента и дивергенции, находим

. (9.22)


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: