Вестник СПбГУ. Прикладная математика. Информатика... 2021. Т. 17. Вып. 2 УДК (517.97+514.822): 531.51 MSC 49S05; 83С40
О плотности потока импульса гравитационного поля
О. И. Дривотин
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
Для цитирования: Дривотин О. И. О плотности потока импульса гравитационного поля // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2021. Т. 17. Вып. 2. С. 137-147. https://doi.org/10.21638/11701/spbul0.2021.204
Импульс рассматривается на основе подхода, широко используемого в вариационном исчислении и теории оптимального управления, в рамках которого исследуется вариация некоторого целевого функционала. В физической теории таким функционалом является действие. При этом вариация действия при сдвиге вдоль некоторого векторного поля может быть выражена в виде поверхностного интеграла некоторой дифференциальной формы третьей степени, через которую и определяется плотность потока импульса. В работе получено уравнение баланса импульса, показывающее, что импульс поля переходит в импульс массы. Рассмотрен пример, иллюстрирующий структуру потока импульса для распределения массы в виде однородного тонкого слоя.
Ключевые слова: вариация действия гравитационного поля, плотность потока импульса гравитационного поля, уравнения баланса импульса, тонкий слой с равномерным распределением массы.
1. Введение. В задачах вариационного исчисления и оптимального управления импульсы широко используются и представляют собой линейные формы, выражающие вариации функционалов через приращения искомых величин [1, 2]. В настоящей работе такой подход применяется для вариационной задачи, где функционалом является действие гравитационного поля. Показано, что вариация рассматриваемого функционала может быть сведена к поверхностному интегралу некоторой дифференциальной формы третьей степени. Впервые эта форма была введена в работе [3]. Записано уравнение баланса импульса и введен вектор, который может считаться вектором плотности потока импульса. При этом плотность потока импульса всегда ассоциируется с некоторым векторным полем. В частности, если рассматриваемое векторное поле образовано касательными векторами к координатным линиям временной координаты, то данный вектор представляет собой вектор плотности потока энергии.
В качестве простейшего примера приведен тонкий слой с равномерным распределением плотности массы, для гравитационного поля которого найдены потоки импульса.
Плотность распределения массы рассматривается на основе ковариантного подхода, разработанного в работах [4-6].
Задача построения тензорного аппарата для описания потока импульса гравитационного поля представляет большой интерес, в частности, в связи с экспериментальным обнаружением гравитационных волн [7]. В работе [8] для характеристики
© Санкт-Петербургский государственный университет, 2021
плотности потока энергии и импульса использован так называемый псевдотензор энергии-импульса. Но предложенный в этой работе псевдотензор не является тензором, и таким образом, его применение в рамках современной теории пространства-времени как нековариантного объекта некорректно с математической точки зрения. Насколько нам известно, ковариантный анализ плотностей потоков импульса ранее не проводился. Отсутствие в физической теории ковариантного объекта, описывающего потоки энергии и импульса гравитационного поля, отмечено, например, в работе [9].
2. Вариация действия для гравитационного поля. Пусть В — некоторая область в четырехмерном пространстве-времени. Рассмотрим гравитационное поле, создаваемое некоторым распределением массы в области В. Гравитационное поле описывается метрическим тензором д, имеющим компоненты д^, г, к = 0, 3.
Запишем действие для изучаемой системы в виде суммы двух членов, первый из которых есть действие собственно гравитационного поля, а второй — действие, характеризующее взаимодействие гравитационного поля и массы:
5 = Бд + Бт. (1)
Действие собственно поля представляется известным классическим выражением [3, 8]
Sg = -aJ ВП, (2)
в
а выражение для Бт будет записано далее. Здесь
с3
\6пС
С — гравитационная постоянная, с — скорость света в вакууме, Н — скалярная кривизна:
Я = д'3 Яц. (3)
При записи (3) и всех остальных формул пользуемся правилом суммирования Эйнштейна, согласно которому по всем совпадающим верхним и нижним индексам производится суммирование. Далее, д3 — контравариантные компоненты метрического тензора, определяемые соотношением
д'3 дзк = Ъ, (4)
в котором 61к — символы Кронекера; Я^* — компоненты тензора Риччи Я, равные Я* = дГ*/дхк - дГкк/дх3 + ГткГт -Г^зПь Г* — символы Кристоффеля 2-го рода; П = А (¿ж1 А 3,х2 А 3,х3 — четырехмерная форма объема; д — определитель
матрицы компонент метрического тензора; х' — координаты в области В.
Варьируя выражение (1) и применяя принцип стационарного действия, можно получить уравнение гравитационного поля. В рассматриваемом случае принцип стационарного действия формулируется следующим образом. Пусть заданы распределение массы в области В и метрический тензор д'к на границе области дВ. Тогда гравитационное поле должно быть таким, чтобы вариация действия при произвольной допустимой вариации метрического тензора внутри области была бы равна нулю.
Рассмотрим сначала вариацию Бд. При варьировании метрического тензора вариации его ковариантных и контравариантных компонент связаны линейным образом в соответствии с равенством (4):
5дг к дк 1 + дт 5дт1 = 0. (5)
а =
Будем варьировать контравариантные компоненты. При этом
бяд = -а ! ((^í7ifcДífcv/^ + í;ifc<^Дífcv/^ + в
+ д'1кИ1к5л/^д)(1х0 Л(1х1 Л(1х2 Л(1х3. (6)
Поскольку ¿д = Д^к¿дгк, где Дгк = ддгк — алгебраическое дополнение элемента матрицы коптравариантпых компонент метрического тензора, имеем равенство
Объединяя в (6) первый и третий члены в подынтегральном выражении, запишем вариацию Б в виде суммы двух членов:
йвд = -а I (ък - 6д"п - а У д'к6Н,кП. (7)
в в
Покажем, что второй член в (7) можно выразить в виде некоторого интеграла по поверхности дБ, которую считаем достаточно гладкой. Имеем равенство
/ д Г д Г \
? о _ л гк И I рш ртр! \
О-Щк — 0 I -Т—1 -7Г~Г ■+" 1 гк1 1т ~ 1 И 1 кт ] ■
Хотя символы Кристоффеля Гк,- те являются тензорами, их вариации ¿Гк представляют собой тензоры, что можно увидеть, анализируя закон преобразования вариаций символов Кристоффеля при переходе от одних координат к другим. Ковариантную производную вариации символа Кристоффеля по координате представим следующим образом:
_ <9(5Гг
Отсюда имеем выражения
I д£Гг- ! ! !
= ^ ; + Гг„(5Г"д, - Т^5Тпк — Т™к5Ты
(ЛГг ^ = _— 4- Гг ЛГ™ — Г™ ЛГг — Г™ лгг
v гк) -т 0хт тп ^ тг пк тк пг'
I д^Гг- 1 1 1
= $хк + ^кп^и ~ Гь<5Г„г - Гк1ЗГ^.
Разность между ними равна
1 1 д5Г1-к д6Г1ч 1 1
- = О^Г + Г1пЗГ7к ~ гизгпк~
- ГП5Г1Ы - Гкм5ГЪ + Г1,5Г1п1 + Г(Г1Ы = = (5Г\к).1 - (¿Г^).к + 5(ГПкГ\п) - ¿(ГП^Пк)
и совпадает с ¿Е^к. Тогда получим, что
д*к5Ък = дгк (¿Г\к)., - д*к (ЗГ1и);к = дгк (¿Гк)., - / (5Г™т);к =
= (д'к ¿Г'к); г - (д1 ¿ГТт); к = ^ г. Здесь ™ 1 = -д1к ¿Г\к + д11 ¿Гт.
Рассмотрим дифференциальную форму третьей степени где * обозначает оператор дуализации (Ходжа) [3], действующий по такому правилу:
( = Ю1.
При этом, как нетрудно убедиться,
где 0гзк1 = 1/0^к1 — контравариантные компоненты формы объема. Вычисляя кова-риантную производную вектора 'Ш, находим, что
= -(№Ы ( );1 = -&3'к1 ( * Ю^к );1 = 1 = 0123.
уМ
Отсюда вытекает, что (Ю1);10 = 3Таким образом,
1(ю1).1 0 = ! з
Б Б дБ
Действие точечной частицы массы т в гравитационном поле, движущейся по траектории х(А), А € [Аа, Ав], равно [3]
л
в
= -тсу у д^кп^ ¿А. (8)
лл
Здесь П(Л) = ¿х/ЗА. Поэтому интеграл не зависит от параметризации. Считаем, что ¿А/Зв > 0, где в — собственное время частицы, Аа < Ав.
Точечная частица характеризуется вырожденным распределением тока массы, носитель которого совпадает с ее траекторией. Будем рассматривать распределение массы, описываемое дифференциальной формой третьей, т. е. максимально возможной, степени. При этом под плотностью тока массы понимаем такую дифференциальную форму третьей степени интегрирование которой по некоторой достаточно гладкой ориентированной поверхности V дает массу, проходящую через эту поверхность в направлении ее ориентации:
т = J 1. (9)
у
Кроме того, будем считать, что в каждой точке некоторой подобласти Б области О задано непрерывно дифференцируемое векторное поле п, определяющее скорость движения массы в рассматриваемой точке. Интегральные кривые этого поля можно считать мировыми линиями, вдоль которых движется масса.
Также примем, что в Б можно построить такую пространственно-подобную поверхность V, которая однократно пересекает все интегральные кривые.
Для вычисления массы по формуле (9) разбиваем поверхность, через которую проходит масса, на ячейки, суммируем по всем ячейкам и производим предельный переход, уменьшая ячейки. Массу в малой ячейке можно представить как некоторую
частицу, которой соответствует некоторая интегральная кривая, проходящая через эту ячейку, для которой действие дается выражением (8). Чтобы получить действие для распределения, описываемого формой 1, просуммируем вместо массы для каждой ячейки действие для частицы, представляемой данной ячейкой. А тогда действие для рассматриваемого распределения принимает вид [3]
-с J | У I (10)
УпБ л /
где V — некоторая поверхность, пересекающая по одному разу все траектории, а А = х(Аа) и В = х(Ав) — точки пересечения траектории с границей области Б.
Заметим, что так как при движении не происходит рождения и уничтожения массы, £их 1 = 0, где £их обозначает производную Ли вдоль векторного поля пл. Это означает, что при интегрировании по любой поверхности, пересекающей однократно все траектории, интеграл (10) имеет одно и то же значение.
Вариация действия Бт при изменении гравитационного поля, в котором движется частица, равна
с Г ( Г и\х)и(х)69гк ^ \
2упб \лА V№кп(л)пкл)
Учитывая, что ^9гк^х)и\\) = Лэ/йХ, имеем выражение
= J (|u^SgгkdX\J-
v пБ \ла
с ¡Л? Хк \ 3
\ / / 53 ¿тпр ¿Хт Л ¿Хп Л ¿ХР
v пБ \ла / 2<р=Р
С I : (1хк Х 3
Лв к
2 1 \ I ~ ¿Х69^ ЗА 53 £1тпр ¿Х"1 Л ¿ХП Л ¿ХР.
Vпб Vа ) 2<п<р
Здесь и = ¿х/Зв, ] — вектор, определяемый из соотношения
<7тпр - -1'т п л/ Ц ~1уп п ¡>.1 ■>
который будем называть вектором плотности тока массы, а £1тпр — символы Леви— Чивита.
А
динату х1, которая вместе с координатами на поверхности V образует систему координат в области Б. Тогда
— ~~ 2
/ хВ \ 3
■I
J ! и13да аX1 ^дз]' £1тпр Лхт А <1хп А д,хр.
упБ \хА ) 2«
&т
Поскольку = 0, \J-gj1 можно занести во внутренний интеграл. Пусть
ориентация поверхности задается вектором и. Отсюда имеем, что
~Т2 J(Зш)1и\*>дил/^дг1тпрс],х1 Л <1хт Л <1хп Л ¿хр.
в
Учитывая соотношение (5), получим, что
_ с /■-.„„ „•,„ 1
-'т
5Бт = ^ / ШдиП = — / Тц6диП,
2 7 " и 2с ^
в в
где ТгI = С ЩI щ — компоненты тензора Т, представляющего собой вариационную производную действия массы по метрике, а и и щ — ковариантные векторы, соответствующие векторам и и ] : jl = ди]г, иг = дгкик.
Объединяя все члены, запишем вариацию действия в виде
5Б = J -а - -Ндгк^ +
в 'дБ
6дгкП + а \ *W. (И)
При заданном условии 6дгк = 0 та границе области В, второй интеграл в (11) равен нулю. А тогда для стационарности действия необходимо обращение в нуль подынтегрального выражения в первом интеграле в (11):
(12)
Уравнение (12) известно как уравнение Эйнштейна.
3. О законах сохранения для гравитационного поля. Пусть носителем распределения массы является некоторая область Б, Б С О, и заданы координаты, которые назовем исходными.
Считаем, что можно выбрать такое векторное поле Vг = йхг /¿1, вдоль которого
и
и,
если перенести, например, точки пересечения интегральных кривых с поверхностью V, провести стационарные кривые, проходящие через эти точки, а затем взять векторы смещения точек построенных кривых относительно соответствующих точек исходных кривых.
Произведем перенос Ли поверхности V и заданного на ней тензора 7 вдоль данного поля на приращение параметра 6¡. При этом перенесенные интегральные кривые, проходящие через точки Ьух € Ьу V, являются стационарными по построению.
Рассмотрим вариацию действия Бт при описываемом переносе. Учтем, что вариация действия отдельной частицы Б\ (8) при изменении траектории при условии, что траектория остается стационарной:
6Б\ = (¿Б, 6ц = -тс(и, v)\A6¡,
поскольку импульс частицы как дифференциал действия [8] при рассматриваемом условии р = йБ\ = -тси.
Тогда изменение действия для массы, проходящей через некоторую трехмерную ячейку на поверхности V, определяемую векторами 5х(1), 5х(2), 5х(3), задающими ее ребра, можно представить в виде
,<5Ж(2),<5Ж(З) = (ГР,^\Л ,^х(2) ,^х(3))
с точностью до величин более высокого порядка малости по 5х^к), к = 1, 2, 3. Суммируя по ячейкам и переходя к пределу, получим выражение
5Бт = !(р,У)\Л 3. (13)
у
Если V есть базисное векторное поле некоторой координаты, то (13) можно рассматривать как изменение соответствующей компоненты импульса для всего распределения массы в целом при переходе от поверхности, на которой расположены начальные точки траекторий, к поверхности, где лежат конечные точки траекторий.
Вариацию Бд при описываемом переносе можно записать, согласно выражению (2), следующим образом:
5Бд =-а ! ^(Ш). (14)
в
Вместе с тем при интегрировании по перенесенным областям в выражениях (2), (4) можно взять в качестве координат точек координаты их прообразов. Соответствующие компоненты скоростей и формы плотности потока массы при переносе не меняются. Интегралы меняются лишь за счет изменения компонент метрического тензора, и вариацию функционала можно представить в виде (11).
Считая также, что распределение массы переносится вместе с гравитационным полем, которое им создается, т. е. уравнения гравитационного поля продолжают выполняться, получим, что в (11) отличен от нуля только поверхностный интеграл. Тогда, учитывая, что в соответствии с (12) К = —Т/(2ас), где Т = дг->Тц = ,
находим, что
/(р,у) + ^ / = I *мг, (15)
у в дв
где
Ш(у) = а(-9гкГ\к + г?т).
Уравнение (15) можно считать уравнением баланса импульса, а вектор Ш — вектором плотности потока импульса. При этом импульс, втекающий в область Б, частично уходит на изменение импульса массы и частично передается гравитационному полю.
4. Пример. Однородный слой в однородном гравитационном поле. Рассмотрим слабое однородное гравитационное поле, описываемое метрическим тензором
\\дгк || = <0.^(1+ 2а(х - х0),-1,-1,-1), х > хо, (16)
в координатах х0 = аЬ, х1 = х, х2 = у, х3 = г. Под слабым полем понимаем поле с координатами, в которых матрица компонент метрического тензора мало отличается от матрицы diag(1, -1, -1, -1). В рассматриваемом случае это означает, что а(х - хо) ^ 1.
Нетрудно убедиться, что тензор (16) удовлетворяет уравнению (12) с точностью до членов первого порядка малости по а(х — х0). Такое поле может создаваться покоящейся плоскостью с равномерным распределением массы, ортогональной оси х и расположенной при х = хо. Из уравнения динамики легко получить, что ускорение движущейся в таком поле массы имеет компоненту вдоль осп х, равную и> = —ас2.
Пусть х0 < 0. Рассмотрим топкий слой массы, ортогональный оси х, находящийся в данном поле при х = 0 и создающий слабое поле, метрический тензор для которого запишем в виде
Jdiag(1 — 2Ьх, — 1, — 1, —1), х< —й, " = jdiag(1 + 2Ьх, —1, —1, —1), х>й.
Здесь 2й — толщина слоя. Постоянная Ь определяется поверхностной плотностью массы а: Ь = 2пОа/с2.
Нетрудно понять, что поле, являющееся результатом наложения внешнего поля (16) и собственного поля рассматриваемого слоя, можно охарактеризовать метрическим тензором
\\9ъ,
{diag(1 — 2(—а + Ь)х, —1, —1, —1), х < —й, diag(1 + 2(а + Ь)х, —1, —1, —1), й> 0.
При этом отличные от нуля символы Кристоффеля 2-го рода равны
I —2(а — Ь)2х, х < —й,
рх _ р0
г 00 = г 0х
—2(а + Ь)2х, х > й.
Вычисляя компоненты вектора плотности потока х-компопепты импульса W(х)
(V = е(ж)), получим W0х■) = 0,
К)
2(а — Ь)2, х < —й, 2(а + Ь)2, х > й.
Удерживая при вычислении только главные члены, имеем выражение
I —2(а — Ь)2йх0 Л 3,уАд,г, х < —й, *W(х) = < у ' п
(х) [— 2(а + Ь) йх Л йу Л ¿г, х > й.
Вычислим поток импульса Ф^, соответствующий вектору W(х), в направлении —е(х) через ячейку, находящуюся при некотором х < —й и образованную векторами смещения е(х0)с5Ь, е(у)5у, е(г)5г,тд$ е(х0), е(х), е(у), е(г) — базисные векторы координат х0, х, у, г. Если задать ориентацию элемента трехмерной поверхности, представляю-
—е(х),
*W(е(х0),е(х),е(у)), поскольку векторы е(х0), е(х), еу образуют на этой поверхности правоориентированный набор. Итак, имеем равенство
Ф( = * W(е(х0)сЖ, е(у)5у, е(г)с6г) = —2(а — Ь)2сЖ5у5г.
Вычислим теперь поток импульса Фг в направлении е(х) через ячейку, находящуюся при некотором х > й, образованную такими же векторами смещения, как
и ранее. Векторы е(х0), е(х), еу образуют на соответствующей ориентированной поверхности левоориентированный набор. Тогда
Здесь 5т — масса ячейки.
Предположим, что описываемый слой неподвижен. Тогда его импульс не меняется и первый член в левой части уравнения (15) равен нулю. Вычислим теперь второй член. Имеем
ТП = (7оос2м°^°л/—д<1х° Л <1х Л ¿у Л ¿г = (1 + 2ах)(? 1хуг(1х0 Л <1х Л ¿у Л ¿г.
Интегрируя по ячейке, образованной векторами смещения с51е(0), 5хе(х), с5уе(у), с5ге(г), получаем
5х
900
(15) выполняется.
Если же рассматриваемый слой ускоряется в приложенном внешнем поле, то второй член в левой части уравнения (15) обращается в нуль, поскольку его можно вычислять в пространственно-временных координатах ускоренной системы отсчета [10], ассоциированной с конгруэнцией интегральных кривых векторного поля и, которые можно характеризовать как мировые линии наблюдателей, формирующих систему отсчета. В этой системе отсчета компоненты метрического тензора не зависят от ж, и = 0. Первый же член в таком случае равен изменению импульса рассматриваемой ячейки слоя. Тогда уравнение (15) означает, что изменение импульса массы ячейки равно потоку импульса гравитационного поля, втекающего в данную ячейку.
5. Заключение. В настоящей работе на основе вариационных методов разработан подход, позволяющий описывать плотности потоков импульса гравитационного поля. При этом импульс рассматривается так же, как в вариационном исчислении и теории оптимального управления, а именно, как вариационная производная функционала действия. В представленном подходе импульс всегда ассоциируется с некоторым векторным полем. Введен тензор, интегрируя который по некоторой замкнутой поверхности можно выразить вариацию действия гравитационного поля внутри этой поверхности. Такой тензор является фактически вариационной производной действия при вариации области.
Рассматривается распределение массы, движение которой происходит с определенной скоростью в каждой точке пространства-времени (пылевидная материя). На основе анализа вариации действия при сдвиге вдоль некоторого векторного поля получено уравнение, содержащее объемный и поверхностные интегралы. Его можно интерпретировать как уравнения баланса импульса массы и гравитационного поля, поскольку одним из членов, дающих вклад в баланс со стороны массы, является изменение импульса этой массы. Потому упомянутый тензор, который интегрируется по поверхности, можно назвать плотностью потока импульса гравитационного поля.
Фг = — (е(х0)с5г, е(у)5у, е(г)5г) = 2(а — Ь)2с5г5у5г.
4=
Суммарный поток, вытекающий из изучаемой ячейки за время 5Ь, равен Ф = Ф; + Фг = 8аЬс5Ь5у5г = —аи)515у5г = —5т->л51.
В качестве примера описан однородный слой массы во внешнем однородном гравитационном поле, для которого уравнение баланса импульса выполняется, что можно рассматривать как подтверждение представленного подхода.
Литература
1. Понтрягин Л. С., Болтянский В. Г., Гамкрелидзе Р. В., Мичценко Е. Ф. Математическая теория оптимальных процессов. М.: Наука, 1976. 392 с.
2. Дривотин О. И. О численном решении задачи оптимального управления на основе метода, использующего вторую вариацию траектории // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2019. Т. 15. Вып. 2. С. 283-295. https://doi.org/10.21638/11702/spbul0.2019.211
3. Дривотин О. И. Математические основы теории поля. СПб.: Изд-во С.-Петерб. ун-та, 2010. 168 с.
4. Drivotin О. I. Covariant formulation of the Vlasov equation // Proceedings of International Particle Accelerators Conference (IPAC'2011). San Sebastian: Kursaal, 2011. P. 2277-2279.
5. Drivotin О. I. Degenerate solutions of the Vlasov equation // Proceedings of Russian Accelerator Conference (RuPAC'2012). St. Petersburg: St. Petersburg State University, 2012. P. 376-378.
6. Дривотин О. И. Ковариантное описание распределений в фазовом пространстве // Вестник Санкт-Петербургского университета. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2016. Вып. 3. С. 39-52. https://doi.org/10.21638/11701/spbul0.2016.304
7. Abbot В. P., Abbot R., Abbot Т. D. et al. LIGO scientific collaboration and Virgo collaboration // Physical Review Letters. 2016. Vol. 116. Iss. 6. P. 061102.
8. Ландау Л. Д., Лифшиц Е. М. Теория поля. М.: Наука, 1973. 504 с.
9. Carrol S. Spacetime and geometry. An introduction to general relativity. San Fransisco: Addison Wesley, 2004. 513 p.
10. Drivotin О. I. Rigorous definition of the reference frame // Вестник Санкт-Петербургского университета. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2014. Вып. 4. С. 25-36.
Статья поступила в редакцию 11 апреля 2020 г. Статья принята к печати 5 апреля 2021 г.
Контактная информация:
Дривотин Олег Игоревич — д-р. физ.-мат. наук, проф.; o.drivotin@spbu.ru
On momentum flow density of the gravitational field
О. I. Drivotin
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Drivotin О. I. On momentum flow density of the gravitational field. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2021, vol. 17, iss. 2, pp. 137-147. https://doi.org/10.21638/11701/spbul0.2021.204 (In Russian)
Momentum is considered on the basis of the approach widely used in the calculus of variations and in the optimal control theory, where variation of a cost functional is investigated. In physical theory, it is the action functional. Action variation under Lie dragging can be expressed as a surface integral of some differential form. The momentum density flow is defined using this form. In this work, the momentum balance equation is obtained. This equation shows that the momentum field transforms into a momentum of a mass. Examples showing the momentum flow structure for a mass distribution representing a uniform thin layer are provided.
Keywords: action variation of the gravitational field, momentum flow density of the gravitational field, momentum balance equation, thin layer with uniform mass distribution.
References
1. Pontryagin L. S., Boltyanskiy V. G., Gamkrelidze R. V., Mischenko E. F. Matematicheskaya teoriya optimal'nyh protsessov [Mathematical theory of optimal processes]. Moscow, Nauka Publ., 1976, 392 p. (In Russian)
2. Drivotin O. I. O chislennom reshenii zadachi optimal'nogo upravleniya na osnove metoda, ispol'zuyushchego vtoruyu variatsiyu trayektorii [On numerical solution of the optimal control problem based on a method using the second variation of a trajectory]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2019, vol. 15, iss. 2, pp. 283-295.
https://doi.org/10.21638/11702/spbul0.2019.211 (In Russian)
[]
St. Petersburg, St. Petersburg University Press, 2010, 168 p. (In Russian)
4. Drivotin O. I. Covariant formulation of the Vlasov equation. Proceedings of International Particle Accelerators Conference (IPAC'2011). San Sebastian, Kursaal Publ., 2011, pp. 2277-2279.
5. Drivotin O.I. Degenerate solutions of the Vlasov equation. Proceedings of Russian Accelerators Conference (RuPAC'2012). St. Petersburg, St. Petersburg State University Publ., 2012, pp. 376-378.
6. Drivotin O. I. Kovariantnoe opisanie raspredeleniy v fazovom prostranstve [Covariant description of phase space distributions]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2016, vol. 12, iss. 3, pp. 39-52. https://doi.org/10.21638/11701/spbul0.2016.304 (In Russian)
7. Abbot B.P., Abbot R., Abbot T.D. et al. LIGO scientific collaboration and Virgo collaboration.
Physical Review Letters, 2016, vol. 116, iss. 6, p. 061102.
[]
(In Russian)
9. Carrol S. Spacetime and geometry. An introduction to general relativity. San Fransisco, Addison Wesley Publ., 2004, 513 p.
10. Drivotin O. I. Rigorous definition of the reference frame. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2014, iss. 4, pp. 25-36.
Received: April 11, 2020. Accepted: April 05, 2021.
Author's information:
Oleg I. Drivotin — Dr. Sci. in Physics and Mathematics, Professor; o.drivotin@spbu.ru