Изменить размер шрифта - +
Для работы с ним ваш браузер должен поддерживать Java. Апплет занимает около 11 Кбайт, так что, если у вас медленный Интернет, загрузка отнимет какое-то время<sup>128</sup>.

Следует заметить, что программа использует алгоритмы вычислений, основанные на общей теории относительности в ее современном виде. Эффекты физики вселенной Диаспоры апплет не учитывает.

Период вращения и расстояние между компонентами системы связаны законом Кеплера. Квадрат периода T пропорционален кубу расстояния а:

Здесь M= Yn<sub>x</sub><sup>j</sup>Tm<sub>2</sub>- совокупная масса компонентов, a G - гравитационная постоянная. Общая энергия излучения гравитационных волн обратно пропорциональна пятой степени расстояния между звездами:

Здесь ц = (т<sub>х</sub>т^)1{т<sub>х</sub> + т<sub>2</sub>) - приведенная масса системы. Краткое обоснование соотношения таково: амплитуда гравитационного излучения каждой звезды пропорциональна ее массе т. и центробежному ускорению ш<sup>2</sup>а., где а - расстояние от центра масс системы, а со = 2л/Т, а также обратно пропорциональна расстоянию до наблюдателя г. По определению Jn<sub>l</sub>O<sub>l</sub> = m<sub>2</sub>a<sub>v</sub> так что звезды генерируют гравитационные импульсы в точности одинаковой амплитуды, а разность фаз между ними составляет 180 градусов. Временная задержка при прохождении звезд не дает волнам скомпенсировать друг друга и поставляет дополнительную разность фаз, пропорциональную (а<sub>х</sub> + а<sub>2</sub>)со. Итак, амплитуда волны, регистрируемой удаленным наблюдателем, составляет

A ~ pa<sup>2</sup>co<sup>3</sup>/r

Мощность излучения пропорциональна квадрату амплитуды. Подставляя вместо со полученное из закона Кеплера соотношение со<sup>2</sup> ~ Mla<sup>2</sup>, получаем

L- А<sup>2</sup>г<sup>2</sup>

L - M<sup>2</sup>H<sup>2</sup>Ia<sup>2</sup>

Точные численные коэффициенты при Ghc можно найти из анализа размерностей, но все же множитель 32/5 требует учета общей теории относительности. Точное описание движения пробных частиц также выходит за рамки настоящего очерка, но оценить искривление орбиты можно, умножив соотношение для А на 772. Получаем

dx/x ~ ра<sup>2</sup>со<sup>2</sup>/г

Воспользовавшись законом Кеплера и подставляя численные коэффициенты при Ghc, имеем

dx/x ~ (G<sup>2</sup>Mp)/(c<sup>4</sup>ar)

Вычисляя время до столкновения звезд, учтем, что суммарная энергия системы (кинетическая плюс потенциальная) составит

E = -GMp/a + Hi<sub>1</sub>W<sup>2</sup>O<sub>1</sub><sup>2</sup>^ + т<sub>2</sub>(л<sup>2</sup>а<sup>2</sup>/2

E = -GM\x/a + In<sub>1</sub>CO<sup>2</sup> а<sup>2</sup>12

E = -GMyJa + GM\x!2a E = -GM\il2a

Дифференциальное уравнение для ее изменения во времени получается из условия равенства этой энергии и энергии гравитационных волн<sup></sup>:

GMu

dE/dt= - da/dt dE/dt = — L

<sup>C</sup>-^-da/dt = —(32C<sup>4</sup>MV)/5c<sup>s</sup>a<sup>5 </sup>da/dt = — (64С<sup>3</sup>М<sup>2</sup>д/5с<sup>5</sup>а<sup>3</sup>) dt/da = -(5c<sup>5</sup> a<sup>3</sup> /MG<sup>3</sup>M<sup>2</sup> ц)

Интегрирование элементарно и привносит дополнительный множитель 1/4:

t = -(5c<sup>5</sup>cd/256G<sup>3</sup>MPn)

Здесь предполагается, что звезды столкнутся при 1=0.

Быстрый переход