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

книги / Физика и философия подобия от преонов до метагалактик

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
41.03 Mб
Скачать

422

§48.1. Уравнения гравитации в пространстве Минковского

 

 

diviУг

н и 00

(547)

 

------------ J 'G t

 

г

dt

где £/°° — скалярная компонента в тензоре Ulk (533) и одновременно плотность

 

энергии гравитационного поля (511),

 

/ — плотность тока массы (плотность импульса вещества),

 

G — гравитационное ускорение,

 

 

SJX^ CU01,

SrY=cU°2, Srz=cU°*, здесь 5 ^ ,5 ^ , Srzкомпоненты вектора Sr ;

UQl, U02,

U03 — компоненты тензора (533),

 

с — скорость света.

 

 

При выводе (547) следует использовать векторное соотношение:

 

 

div[A х В] = В-rot А - А-тоХВ.

 

Смысл равенства (547) заключается в том, что поток энергии в некоторый обьем через его поверхность приводит к увеличению гравитационной энергии в этом объе­ ме и совершению работы по ускорению вещества.

В тензоре tlk (537) содержится вектор плотности потока механической энергии

Sv :

 

 

s u = f °_С у 1^1 у > или

s S = ct°p,

где индексу) = 1,2 ,3 ,

= Sw ,

S V2 = S ^ ,

top — компоненты тензора (537).

Этот вектор несколько отличается от обычно упоминаемого вектора Умова.

Дивергенция Sv с учетом (540) равна:

 

 

 

 

 

 

 

-

Рос

+

р0 dh.r |

d l p

y +РЛ

 

э / р 0сг + Р 0\

 

 

1 - v 2/ c 2

<*(i -

v 2/ c 2 )

 

a r(i

-

V 2/ C 2 j

 

= - * 1

+ & f l + й. divV + ± (P.Qc1 + A )

d(L - Pa)

 

 

dt

i - V 2/c 2

 

d t \ l - V

2/c 2)

dt

здесь/ 00 — скалярная компонента в тензоре tik (537).

 

 

Сложим дивергенции векторов

Sr и

Sv

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

U00 и /°° по времени и используем (542), (545), (547):

 

 

 

 

 

 

 

 

ЛГ/00

 

 

 

а/00

 

 

 

 

 

 

divSr + —— + divSu + —

=

 

 

 

 

 

г

dt

 

v

 

dt

 

 

-----/• G +

PQC2 + PQ diwV +

d (poc

+

A

 

d/

 

 

 

1

- V2/c2

 

d/^1 -

V2/с 2

 

_

P QC

- divF

+ ~

 

 

Pog

 

 

1 "

V*/с 2

 

d t [ l - V 2/c 2)

V T - K2/c 2 d t[ J — V2/ c 2

 

_

Ppg

■divF +

 

 

4

 

Poc*

1 = 0 .

 

 

 

 

V1

 

 

 

 

i - ^ / c 2

 

 

 

—V 2/c 2;

§48.1. Уравнения гравитации в пространстве Минковского

423

Вспоминая, что е

РоС2

=

pc2 — полная энергия единицы объема

- V2/c 2

Vl

 

 

 

вещества без учета давления, для дивергенции скорости находим уравнение неразрывности, эквивалентное (539):

divr = - I *

= - ! £ .

г dt

р dt

Кроме этого, дивергенция суммарного потока механической и гравитационной энергии из элемента объема равняется скорости суммарных потерь энергии в этом элементе объема:

div(5f + S„) = - —(U00 + /°°).

(548)

dt

 

Если рассматривать еще и электромагнитное поле, то в левую часть (548) под знак дивергенции следует добавить вектор Пойнтинга Sp = х Я], где Е%И — соответственно электрическая и магнитная напряженности поля, а в правую часть — плотность электромагнитной энергии или скалярную компоненту тензора плотности

энергии-импульса электромагнитного поля, равную fV00 = E D3 + В Н ), где

D9— электрическое смещение, В — магнитная индукция.

Учитывая (547) и аналогичное соотношение для электромагнитного поля, полу­

чим:

 

 

 

 

 

 

 

dt^

+ J*G + j ‘E ,

 

 

 

 

div S v

 

 

 

 

dt

 

 

 

гдеj — плотность электрического тока,

 

 

 

Е — напряженность электрического поля.

 

 

Рассмотрим

теперь

скалярные

компоненты

тензоров

плотности

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

 

£

= # 00

+ f00 + W 00 =

(549)

= -

(G2 + с2а 2) +

V /с

+ L - Р0 + U E -D9 + В Н) = const,

 

ъяу

1 -

2

 

здесь В — полная плотность энергии вещества и поля, у — гравитационная постоянная,

G — гравитационное ускорение, с — скорость света,

Q — гравитационное кручение,

р 0 — плотность покоящегося вещества, PQ— изотропное давление,

L — величина согласно (546),

V—средняя скорость движения вещества в данном обьеме. Используя (548) и учитывая электромагнитное поле, можно записать:

-— = div(SV + S„ + Sr) = divS. dt

Используем теперь следующие выражения:

424

§48.1. Уравнения гравитации в пространстве Минковского

 

= — + F*gradS

согласно (540),

 

dt

dt

 

 

 

 

div(£F) = £divF + F*grad£,

 

divF

_1_ ф

 

P о

 

где p =

т/l - KJ/c J

 

 

p dt

 

 

В результате получим;

 

 

 

 

- — + div(£F) +

p d t

= div5.

 

 

Л

 

Предположим теперь, что рассматривается небольшой кусок вещества, изолиро­ ванный от окружающей среды. Тогда изменением внешних полей можно пренебречь

— они не проникают в такую систему благодаря изоляции. Полная внутренняя энер­ гия системы при ее неизменном объеме будет также постоянна и можно положить

Л _

 

 

 

 

 

 

= 0. В этом случае имеем:

 

 

 

dt

 

 

£ do

 

 

 

 

div(£F -

S)

£ и S\

 

 

 

+ —— = 0, или подставляя

 

 

 

 

 

р

dt

 

 

 

 

div[(I

-

?0)F

+ U*V + IV™V - Sr -

5,1 +

U<*-+ t ^ + W ^ dp = Q

 

L

 

 

 

1

p

dt

To, что находится под знаком дивергенции, следует по-видимому отождествить с вектором типа Умова или вектором плотности потока энергии, возникающим вслед­ ствие колебаний плотности вещества р. В замкнутой системе колебания приводят к изменениям локальных плотности вещества, давления, плотности потенциальной энергии и излучения. При этом надо учесть, что в небольшом куске вещества можно пренебречь обычной гравитацией, но необходимо учитывать ядерную гравитацию между атомами вещества, то есть в U00 и Sr должна входить постоянная ядерной 1равитации Г из (422).

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

J = Роэ

с Роэ

у р о

*) = (СРэ . Л

= (

’ -Jl - К 7 с

 

л/l - К 2/С 2

 

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

j — плотность электрического тока, V— скорость движения зарядов.

Четырехмерный вектор электромагнитного потенциала имеет вид:

 

 

 

 

Л ' = ( £ , Л ) ,

(550)

 

 

 

С

 

где — электрический скалярный потенциал,

 

А — векторный потенциал электромагнитного поля.

 

Векторы

/

и

А1 записываются практически также,

как векторы

1равитационного поля J 1 и /У, а волновое уравнение электромагнитного поля с точностью до знака и коэффициентов подобно (528):

§48.2. Три избранные задачи

425

□ V = J - f , с2е0

здесь г0 — электрическая постоянная.

Мы видим, что без всяких противоречий электромагнитное и гравитационное поля описываются одними и теми же формулами, следовательно, они образуют еди­ ный объект — электрогравитационное поле. Лишь одновременный учет всех компо­ нент такого поля позволяет правильно записывать законы сохранения энергии, импульса и момента импульса и находить силы, действующие между объектами.

§48.2. Три избранные задачи

1.Рассмотрим задачу о гравитацион­

ном взаимодействии малой пробной массы Л/, и большой массы М2 в двух

системах отсчета — в системе К \ в кото­ рой обе массы покоятся, а масса М2на­ ходится в начале координат, и в системе К, которая движется относительно К' влево вдоль оси X ' (рисунок 85) со ско­

ростью V0. Компоненты силы тяготения

Рис. 85. Массы М х и М г покоятся в системе

между массами в системе отсчета К'

отсчета К', но двигаются вправо относительно

имеют вид:

системы отсчета К со скоростью V0 (или система

 

отсчета К движется влево относительно К' со

 

скоростью V0).

уМ, М2х'

уМ. М7х'

 

 

(551)

F

(,'’ + ;-+ zT

 

 

 

 

Г. _

у-Л/, Мг у'

 

_

уМ хМ2т!

 

F y ~ ~

if'1"

z

к ,ъ >

 

 

 

здесь у — гравитационная постоянная,

х\ у \ z' — координаты точки, где находится масса Мх, R! —расстояние между массами М2иМ х.

Перейдем в систему отсчета К>в которой обе массы вместе с системой отсчета К' движутся вправо вдоль оси X со скоростью К0. Согласно (524) и (525) потенциалы от движущейся массы М2 в системе К будут таковы:

уМ г

г м г

Ф = -

 

 

1( х - К 0 2 + у 2 + г2

R j l - к 1/ С1 ’

К22 i 1 -

у 2/с 1

 

Dx = ^

, Dr =D Z = 0 ,

 

здесь предполагается, что при /= 0 масса М2находилась в начале координат сис­ темы отсчета К, а расстояние R до массы ЛГ, в системе отсчета К учитывает эффект сокращения длины вдоль оси X за счет движения массы Мх. Определим теперь ком­ поненты гравитационного ускорения G и кручения Q с помощью (506):

 

Эф

дРх _

уМ 2(х - V0t)

уМ 2(х

- VQt)V02 _

*

Эх

д/

Л3(1 - у 2/с2)'-5

Л1 с2(1 -

К2/с2)15

426 §48.2. Три избранные задачи

_ уМ 2(х - Г,0 R ^ l - V i/с г '

Так как Dr и Dz равны нулю, то компоненты GYи Gz получаются проще:

 

 

 

G

-

дхР -

Vм г У

 

 

 

 

 

 

Г

 

Эу

Л371 - УЦс1 '

 

 

 

 

 

G

=

_

= _

 

УМ г г

 

 

 

 

 

 

z

 

sz

л 3Vi - у,2/с 1

 

 

 

Компоненты кручения по (506) равны:

 

 

 

 

 

 

 

 

 

 

 

Э^

dz

 

 

 

 

 

 

32).

 

dDz

dDx

YM2zV t

 

 

 

 

Qr =

 

 

dX

bz

 

R’C2J 1 - К,2 /с 2 ’

 

 

 

 

 

dz

 

 

 

 

 

 

 

-

aZ)r

 

dDy

dD* -

 

 

 

 

 

 

aZ)x _

R 'c 'y j T ^ v f f ?

 

 

 

 

z

Эх

 

a?

aj;

 

 

 

 

С помощью (505) находим компоненты силы, действующей на массу ЖГ,,

двигающейся со скоростью V= (Vx , VY, Vz), причем Vx = V0, VY= VZ = Q:

 

 

 

Fx = M XGX + Mx(VrQz -

Vz Qy) = M XGX = -

 

(552)

2V

=

+ MX(VZQX -

VXQZ) = M xGr -

M, VXQZ = -

~

 

^

 

+ З Д 2 2 , -

Vr Qx) = M XGZ + M xVx Qr = - yM ' Ml Zj l

V° /C\

 

Согласно (393) поперечные по отношению к скорости VQприращения импульсов

dp

и dp, создаваемые силами в системах отсчета К' и К, должны быть равны:

 

 

 

 

 

 

 

dpY ~ dpу,

dpz

—dp2 *

 

 

 

Так как по (395) F = — , то можно записать:

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

F; dt' = Fr dt,

/ 5 d f = Fz dt.

 

(553)

 

Система отсчета

К

движется

влево относительно системы

отсчета

К ' со

скоростью VQ, поэтому время в К замедлено по отношению к К' согласно (364):

 

 

 

 

 

 

dt

=

dt'

 

 

(554)

 

 

 

 

 

 

 

 

 

Vi - v t / c 1'

Подставим значения сил (551), (552) и интервалы времени (554) в (553), после сокращения одинаковых величин получим:

У* У

§48.2. Три избранные задачи

427

Поперечные к скорости К0 координаты массы М] в обеих системах отсчета совпадают: у' = у, т! = z, следовательно, R' = R. Раскрывая значения R' и R через координаты, получим:

что совпадает с преобразованиями Лоренца (386). Кроме этого, при R = R' сила Fj[ в (551) равняется силе Fx в (552), если координату х' выразить через х и вре­ мя t и учесть равенства у 1= у, т! = z. Следовательно, использование гравитацион­ ных потенциалов хр к D, гравитационного ускорения G и кручения Q в данной задаче находится в согласии с теорией относительности Эйнштейна.

При получении компонент силы Fy и Fz в системе отсчета К было видно, что вследствие движения параллельных масс Мх и М2 возникает сила поперечного от­ талкивания (в членах, содержащих кручение Qy и Qz). Данная сила отталкивания вполне эквивалентна силе отталкивания F{ в (503). Благодаря этой силе при боль­ ших скоростях движения параллельно движущиеся массы притягиваются слабее друг к другу — суммарная сила стремится к нулю по мере приближения скорости к скоро­ сти света. Как следствие этого, две быстровращающиеся в одну сторону вокруг одной оси нейтронные звезды должны притягиваться друг к другу слабее, чем если бы они не вращались. В § 45 обсуждалась связь нуклонов в дейтроне, для которой известно, что нуклоны вращаются в одну сторону. Полагая, что все вышесказанное справедли­ во и для ядерной гравитации (с заменой постоянной у на постоянную ядерной гра­ витации Г из (422) ), приходим к тому, что вращающиеся нуклоны в дейтроне притягиваются ядерной гравитацией слабее, чем если бы они были неподвижны.

2. В качестве одного из применений уравнений (508) рассмотрим стационарное решение для гравитационного скалярного потенциала хр галактики, имеющей массивное ядро и разреженное гало. Так как хр не зависит от времени, уравнение для потенциала примет вид:

Ахр = 4яу р ,

(555)

где у — гравитационная постоянная, р — средняя пространственная плотность вещества звезд в галактике.

Предположим, что галактика имеет форму диска, тогда хр в основном зависит

только от одной координаты — радиального расстояния г. Разделим

хр и р

на две

компоненты: хр{г) = хр0(г) + хр\г), р(г) = р0(г) + р\г)у где хр0(г) и

р 0(г)

задают

основное распределение гравитационного потенциала и плотности звездной массы, а величины хр\г) и р\г) являются некоторыми небольшими стационарными отклонениями от основного распределения. Зная хр'(г), можно найти потенциальную энергию звезды с массой Ms , отсчитываемую от положения равновесия при потенциале xpQ(r): U - Ms хр\г). Звезды, движущиеся вокруг ядра галактики, образуют звездный газ, который можно охарактеризовать некоторой звездной температурой Т. С помощью величин U и Т можно, как это следует из статистической механики, найти распределение плотности р'(г) в пространстве:

(556)

здесь р '— амплитуда отклонения плотности вещества от среднего значения р 0(г), М5— средняя масса звезды в галактике,

Ks — звездная постоянная Больцмана (183),

Т — звездная температура кинетического движения звезд.

428

§48.2. Три избранные задачи

Ввиду линейности (555) можно разделить на два уравнения для основного потенциала тр0(г) и возмущения гр'(г):

dr1

d r

Мы будем решать только второе уравнение, для которого известно распределение плотности в виде (556). Пользуясь малостью величины М г)> разложим экспоненту в (556) в ряд Тейлора и возьмем самые большие члены:

р'(г) ~ р' -

p'M s tp\r)

(557)

 

KS T

С учетом (557) уравнение для t//(r) примет ввд:

d2rl>\r) + 4n yp 'M s Ар'(г) = 4п у р \

dr2

KS T

 

Общее решение этого уравнения таково:

 

KS T

+ A sin(r

\А лур'М 5 + §).

M r) = M s

KS T

где Аур некоторые константы.

Подставляя М г) в (557), найдем добавку к основному распределению плотности звезд:

Р’(г) =

- A p’M s sin(г

\4 я ур 'М 5

+ P ) .

 

KST

1 KS T

 

Колебания плотности

p'(r) на фоне

некоторой

средней плотности р 0(г)

выглядит как кольца в диске вокруг ядра галактики. Учитывая периодичность

синуса, можно найти расстояние

между двумя соседними максимумами

плотности, один из которых находится в точке г, а другой в точке г+дг:

 

 

sin(r

l4nyp'M s

Р)

= sin

14л уp'Mr

1

 

 

 

+

+ dr)

+ P

 

 

 

J

Kt T

 

 

 

У~~кГт

j =

 

 

 

 

= sin (г

!4 xyp 'M s

+ P + 2л),

 

 

 

 

 

 

1

KS T

 

 

 

 

 

 

 

dr

=

л К 5Т

 

(558)

 

 

 

УР'М3

 

 

 

 

 

 

 

 

 

Оценим dr в нашей Галактике, полагая, что Ks ~ 1033 Дж/К

по (182),

10б К

по (164) — (174),

р’ ~ 5-10‘ 21

кг/м3

— средняя плотность

вещества звезд

по

пространству вокруг Солнца, средняя масса звезды М5 ~ Ю30 кг (смотри

§

21),

у = 6,67 • 10“ *1м3• кг" 1• с-2 — гравитационная постоянная. В результате находим:

 

 

dr ~ Ю20 метра ~ 3 кпк.

Спиральная структура нашей Галактики действительно имеет характерный шаг порядка 3 кпк. Выбросы вещества и вращение галактических ядер превращает кольца в спирали, но в некоторых галактиках кольца плотности наблюдаются очень хорошо.

По-видимому, можно назвать по крайней мере еще два примера явлений, в которых кольца плотности должны играть существенную роль. В частности, хорошо

 

 

§48.2. Три избранные задачи

429

 

 

 

 

известны кольца вокруг Сатурна и дру­

 

 

гих больших планет, состоящие из пыли

 

 

и обледенелых обломков различных раз­

 

Р(ху у, Z)

меров вплоть до нескольких метров в по­

 

перечнике. В качестве другого примера

 

 

возьмем

упорядоченное

расположение

 

 

планет Солнечной системы, которое мо­

 

 

жет быть следствием квазипериодиче-

 

 

ского

распределения

плотности

 

 

вещества в допланетном диске. В резуль­

 

 

тате планетезимали вначале оформля­

 

 

лись в кольцевые структуры, а затем в

 

 

отдельные сгущения вещества. Ход из­

 

 

менения наблюдаемого расстояния меж­

 

 

ду соседними планетами от текущего

 

 

радиуса скорее всего был задан еще на

Рис.86.0пределение фавитационного излуче­

стадии существования сплошного газо­

ния в точке Рот планеты, вращающейся по кру­

во-пылевого диска, при

этом важную

говой орбите вокруг звезды в плоскости XOY.

роль играло падение температуры веще­ ства диска по мере удаления от Солнца.

3.Решим задачу о гравитационном излучении от планеты, обращающейся вокруг

звезды по круговой орбите с радиусом (смотри рисунок 8 6 ). Звезда находится в начале координат в точке 0, планета в момент наблюдения — в точке L(xQ, у0 , z0), потенциалы гравитационного поля ищутся в произвольной удаленной точке Р(х, у , z). Для кругового движения планеты в плоскости XOY с угловой частотой о) можно записать:

Орбитальная скорость:

V0 = Rgw,

 

Координаты планеты:

xQ= R$ cos <о t, yQ= R0sin <o/, 2^ = 0.

Скорости планеты: Vx

=

= - (oi^sin^,

VY = — = (oRgCoscot,

 

 

at

dt

Vx + Vy2 = VI

Расстояние от планеты до точки наблюдения Р(х, у, г) с учетом запаздывания передачи гравитационного возмущения равно:

 

 

 

г'

= л](х

~ x 0(t'))2 + ( у - у , ( П ) 2 + Z2,

здесь t

,

= t

г'

R

запаздывающее время,

 

---------t

------ —

сс

с— скорость распространения гравитационного возмущения,

R = -yjx2 + у 2 + z2 — среднее расстояние от планеты до точки Р.

Если скорость VQ невелика по сравнению со скоростью света и скоростью гравитонов, то соотношения для потенциалов отдельной движущейся массы (524) и (525) можно упростить:

V = -

уМУ

 

Dx = у Ma)RQsin cot' А

Y M W R Q COScat'

430

§48.2. Три избранные задачи

где Ц— скалярный потенциал, определяемый в точке Р в момент времени / через положение планеты в более ранний момент времени

у — гравитационная постоянная, М — масса планеты,

D векторный потенциал, V— вектор скорости планеты.

Поскольку планета движется только в плоскости XOY, то Vz = 0, значит и Dz = 0, и используя выражение Q = rot D из (506), можно найти вначале компоненту кручения Qx в точке Р :

0 _ dDz

BDY _

BDy

_ Y M W R Q d / coso)(t

R /c )s

dy

dZ

dZ

CJ dZ

r '

Преобразуем частную производную и вычислим ее отдельные части:

dz

 

 

 

г dz

 

 

 

 

 

+ cos[a>(? -

Л /с )] 1 (1 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

г'

 

 

 

 

d_

 

^

 

_ ^ , + / + г. А )] =

(DZ s i n (Dt'

 

- < / 0 ] ш у

 

 

 

T c

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

]

dr'

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 ) = _______.

 

 

 

 

 

 

 

 

 

 

 

 

dZ

Г'

 

 

r '2

dz

 

 

 

 

 

 

 

Собирая все вместе, находим для Q:

 

 

 

 

 

 

 

 

 

 

 

 

^

_

yM(oR0/Q)zsino)t'

 

cosa>t'dr\

 

 

 

 

 

 

 

 

 

с

2

'

п

 

I

 

г

ft

Г"'*

 

 

 

 

 

 

 

 

 

 

 

Rcr

 

 

 

dz

 

 

 

 

Подобным же образом находим другие компоненты кручения:

 

 

 

о =

bD*

_

dDz

-

Лгг

-

У^й>Л0

'

о)z cos*»' _

sin cpt' dr'

Y

fa

 

dx

 

 

 

гЛ

R cr'

 

 

~'2

 

 

 

 

dz

 

 

 

 

 

 

 

 

dz

 

Q

-

dA

_

3-A =

__ /Л/й>Л0 .COXSinft)/' _

cosw/'dr'

_

 

г

dx

 

dy

 

 

 

c 2

 

 

Л е г '

 

r '2 dx

 

 

 

 

 

 

 

wycostyf'

 

sin w/'dr'

 

 

 

 

 

 

 

 

 

 

 

 

Лег'

 

 

 

 

).

 

 

 

 

 

 

 

 

 

 

 

 

 

r '2 dy

 

 

 

 

 

 

Поскольку все отношения z/Л , х / Л , у/Л меньше единицы, все частные произ-

d r' d r' d r'

водные — ,— ,— также меньше единицы, то при достаточно больших расстояниях dz dx dy

г' члены, содержащие в знаменателе г'2, будут малы и ими можно будет пренебречь. Выпишем оставшиеся члены в компонентах кручения, которые будут в основном от­ ветственны за 1равитационное излучение:

_ уЛ/а>2Л0г smart'

Л/ _

уМо)2R^ z coscot'

 

Л ^ ?

 

Щ

Л ? ?

Qz = - -

A

(xsinaj/'

- ycosatf')*

 

z

R r ’c*

 

 

 

Вычислим теперь х-компоненту ускорения гравитационного поля по (506):

 

 

 

§48.2. Три избранные задачи

 

431

q _

__

dtp __

dDx

___ д .

уМ .

a yAfa>J?0 sin6>/' _

*

~

дх

dt

дх{

г' }

дГ

г’с2

 

}

=

_

УМ дг' __

УМ(о2RpCosoPt'

yMajR0smcot' dr'

 

 

 

г'2 дх

 

г'с2

 

 

г'2 с2

dt

*

Вновь пренебрегая всеми членами, содержащими г'1в знаменателе, для главных членов компонент Gx и Gy гравитационного ускорения получаем:

п*

-

yMo)2R0 coscot'

 

ь х

гrсi

9

 

__

y^/VJ^sineitf'

 

С /у

---------------- п

«

Поскольку Dz = 0, то Gz = - * t

= _ Y M d S =

_ о .

az

Поток гравитационной энергии в единицу времени через единичную площадь, перпендикулярную потоку, определяется вектором (512), подобным вектору Пойнтинга для электромагнитного поля:

 

 

 

S r = —— [С х О].

 

 

 

 

 

 

г

4 л /

1

 

 

 

Найдем компоненты вектора Sr для полей G' и Q ':

 

 

SJX =

-

 

= “ УМ ы

До.. ^

^ ^ > -

ycos<y/')sinotf',

 

4 лгу

 

 

4яс г

R

 

 

 

Sly

------— (Gz i3y

- G'XQ'Z) =

У—. -^ .J^0 (x sin wt‘

-

cosM/') coscu/',

 

4 я у

*

 

4 л с У 2Л

 

 

 

 

*^12

//1» n ;

Л '

_ _

У-Л/ Cl> R Q Z

 

(<?x i2r

GyQji) -

-4jtcy

2R

 

 

 

4 л у

 

 

 

 

 

На больших расстояниях, когда Л много больше радиуса орбиты i?0, можно положить г' ~ Л. Модуль вектора плотности потока энергии будет равен:

S r = <jslx + S$y + S]z =

ypf g ^(xsineDt* - у coscot1)2 + z2.

 

4 л с л

После усреднения по времени приблизительно получаем:

«

yM 2co4R20

 

г

47ic3R2

 

Из всех компонент вектора плотности потока энергии только

постоянна, так

что вектор Sr в каждой точке осциллирует вокруг компоненты Sn . Зададим в сфери­ ческих координатах не зависящую от времени составляющую радиальной компонен­ ты вектора Sr с тем, чтобы проинтегрировать поток энергии по поверхности сферы с радиусом R, для чего спроектируем s * на радиус R :

4Ж С К

здесь Q — угол, отсчитываемый от оси Z; Q является одной из сферических координат и изменяется от 0 до ж.