автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.16, диссертация на тему:Применение математических методов для интегрирования нелинейных уравнений теории гравитации и анализ их решений

ВВЕДЕНИЕ.

Глава 1. Математические модели теории гравитации . 9 5 1. Уравнения гравитационного поля.

1.1 Общая Теория Относительности.

1.2 Скалярно-тензорные теории. 2. Пространство Вайдья. 3. Метрика статического сферически симметричного скалярного источника.

5 4. Математические методы исследования движения частиц в гравитационных полях.

Глава 2. Применение математических методов для анализа пространств со скалярным полем.

5. Интегрирование уравнений для изотропных геодезических.

6. Интегрирование уравнений для времениподобных геодезических.

7. Система уравнений, определяющих видимое положение звезд.

8. Алгоритм численных расчетов и полученные результаты. 9. Основные наблюдательные проявления предсказываемых гравитационных эффектов.

Глава 3. Анализ математической модели гравитационного поля массивной скалярной звезды . . 50 10. Интегрирование уравнений движения безмассовых частиц в гравитационном поле массивной скалярной звезды.

§ 11. Интегрирование уравнений движения массивных частиц в гравитационном поле массивной скалярной звезды.

§ 12. Условие отсутствия у фотонов круговых орбит в гравитационном поле массивной скалярной звезды.

§ 13. Математическое моделирование искривления лучей света в гравитационном поле массивной скалярной звезды.

§ 14. Разработка алгоритма численных расчетов и полученные результаты.

Глава 4. Применение математических методов для анализа гравитационных эффектов в пространствах специального вида.

§ 15. Применение метода Боголюбова-Митропольского для исследования движения частиц в пространстве Вайдья.

§ 16. Применение параметрического метода для исследования радиальных геодезических в метрике Вайдья.

§17. Применение методов тензорной алгебры для расчета взаимодействия гравитационного поля с электромагнитной волной в кольцевом лазере.

С помощью математического моделирования [1] изучают разнообразные явления природы и приближенно описывают их на языке математических уравнений. Однако при разработке математической модели изучаемых явлений невозможно учесть все физические факторы процесса, поэтому приходится выделять главные и отбрасывать второстепенные. Затем, разработав аналитический или численный метод решения полученных уравнений, составляют программы и, проводя вычислительный эксперимент [2], изучают явление в целом. Для этого используются компьютеры, которые достигли производительности больших вычислительных машин.

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

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

Следует отметить, что применение численных методов решения уравнений гравитации и проведение вычислительного эксперимента [3] в области гравитации ограничено случаями слабых полей и медленных движений источников, создающих эти поля. Основной причиной этого является большое число слагаемых в неупрощенных уравнениях теории гравитационного поля. Как показано в нашей работе [4], в общем случае, когда отличны от нуля все десять компонент метрического тензора и этот тензор зависит от времени и всех трех пространственных координат, каждое из десяти уравнений теории Эйнштейна в координатном базисе содержит более десяти тысяч слагаемых. Вполне понятно, что не только решать эти уравнения, но даже просто выписать их не в символическом виде, а в явном виде (относительно компонент метрического тензора), без использования системы аналитических вычислений на компьютере не имеет смысла. А так как для анализа гравитационных явлений обычно требуется знание метрического тензора в очень большой пространственно-временной области, то применение вычислительного эксперимента для изучения релятивистских гравитационных явлений в настоящее время не получило широкого распространения, несмотря на резко возросшую мощь вычислительной техники.

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

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

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

Во второй главе диссертации проведен математический анализ пространства - времени, создаваемого сферически симметричным безмассовым скалярным полем. Гравитационные свойства такого источника ранее не исследовались и впервые это было сделано в наших работах [5-8]. В этом случае нелинейные уравнения геодезических допускают аналитическое решение и позволяют найти законы движения и уравнения траекторий для массивных и безмассовых частиц.

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

В § 6 нами проинтегрированы системы нелинейных дифференциальных уравнений, описывающих движение массивных частиц в гравитационном поле безмассового скалярного источника, и найдены их законы движения и уравнения траекторий. Показано, что на нерадиаль-но движущиеся массивные частицы в этом гравитационном поле также действуют силы гравитационного отталкивания. Кроме того впервые обнаружено, что на покоящиеся частицы это гравитационное поле вообще не оказывает никакого действия и они остаются в покое. Отмечено, что ещё одним точным решением уравнений Эйнштейна, обладающим таким свойством, является метрика Брдички [9].

В § 7 проведены анализ и выбор системы уравнений, которые описывают распространение электромагнитных лучей в гравитационном поле сферически симметричного источника. На основе этих уравнений в § 8 разработаны алгоритмы для численных расчетов искажений, вносимых гравитационным полем скалярного источника в угловые координаты звёзд при астрометрических измерениях. Проведенные по этому алгоритму численные вычисления позволили нам [8] получить величину угловых искажений 6(р в зависимости от параметров задачи и тем самым выработать способ исключения этих искажений из данных астрометрических измерений.

В § 9 на основе проведенного нами анализа [10, 11] перечислены основные наблюдательные проявления предсказываемых эффектов в гравитационном поле сферически симметричного источника и сформулированы рекомендации по их поиску в астрофизических условиях.

В третьей главе диссертации проведен анализ геодезических в гравитационном поле массивной скалярной звезды. В § 10 нами [11] решена система нелинейных уравнений геодезических для безмассовой частицы, движущейся в гравитационном поле звезды, имеющей массу и скалярный заряд, и найдены уравнения ее траектории в квадратурах.

В § 11 аналогичная задача нами решена для массивной частицы. На основе полученных решений в § 12 проведен анализ и найдено условие д2 < 3при выполнении которого у безмассовых частиц возможны круговые орбиты в гравитационном поле массивной скалярной звезды.

В § 13 построена математическая модель искривления лучей света в гравитационном поле массивной скалярной звезды, в § 14 нами разработан алгоритм и проведены численные расчеты гравитационного искривления световых лучей, проходящих мимо этой звезды [11].

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

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

В § 16 для поиска не приближенного, а точного решения этой задачи в случае радиального движения частиц нами был применен [12] параметрический метод, впервые разработанный в работе [13].

В последнем § 17 проведено применение тензорных соотношений, доказанных в работах И.П.Денисовой [14], к расчету взаимодействия гравитационного поля с электромагнитной волной в кольцевом лазере. На основе полученных результатов нами [15] предложен эксперимент по проверке принципа эквивалентности Эйнштейна для фотонов с относительной точностью до Ю-21, что на семь порядков превышает точность планирующихся в США космических экспериментов по проверке принципа эквивалентности в рамках программы Gravity-Probe-B [16].

В заключении приведены основные результаты, полученные в диссертации и выносимые на защиту.

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

XV Международной конференции по Общей теории относительности и теории гравитации, Декабрь 16-21, 1997, (Pune, Индия),

IX научной конференции "Современные проблемы вычислительной математики и математической физики", посвященной восьмидесятилетию академика РАН Александра Андреевича Самарского, Москва, МГУ, 22-25 февраля 1999 года,

X Российской гравитационной конференции "Теоретические и экспериментальные проблемы Общей теории относительности и гравитации", Секция "Релятивистская астрофизика и космология", Владимир, 20-27 июня 1999 года,

Боголюбовской конференции по проблемам математической и теоретической физики, (Москва - Дубна - Киев), 27 сентября - 6 октября 1999 года и опубликованы в работах [5, 6, 7, 8, 10, 11, 12, 15].

ЗАКЛЮЧЕНИЕ.

Перечислим основные результаты, составляющие содержание настоящей диссертации и выносимые на защиту.

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

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

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

4. Построена математическая модель распространения лучей света в гравитационном поле массивной скалярной звезды и на основе вычислительного эксперимента впервые получена и проанализирована зависимость угла отклонения от прицельного расстояния и соотношения между массой звезды и ее скалярным зарядом.

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

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

7. На основе методов тензорной алгебры проведен расчет взаимодействия гравитационного поля с электромагнитной волной в кольцевом лазере и получена оценка точности измерений в эксперименте по проверке принципа эквивалентности Эйнштейна для фотонов.

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

1. Самарский. А.А. Введение в численные методы. М.: Наука, 1997. - 239 с.

2. Самарский А.А. Численные методы. М.: Наука, 1982. - 272 с.

3. Самарский А.А., Попов Ю.П. "Вычислительный эксперимент в физике", в книге "Наука и человечество". М.: АН СССР, 1975. с. 281-291.

4. Денисова И.П., Зубрило А.А. Применение компьютерной алгебры Reduce для интегрирования уравнений теории гравитации методом Ньюмена Пенроуза. // Математическое моделирование, 2000, № 2, с. 59-67.

5. Зубрило А.А. Об одном частном решении в общей теории относительности. // Теоретическая и математическая физика, 1999, т. 118, № 2, с. 317-320.

6. Denisova I.P., Mehta B.V., Zubrilo A. A. The investigation of the model of gravitational repulsion in Einstein's general theory of relativity. // General Relativity and Gravitation, 1999, v. 31, № 6, p. 821-837.

7. Денисова И.П., Зубрило А.А. Математическая модель рассеивающей гравитационной линзы. // Математическое моделирование, 2000, № 2, с. 68-74

8. Brdicka М. On gravitational waves. // Proc. Roy. Irish. Acad. 1951, v. A54, p. 137.

9. Денисова И.П., Зубрило А.А. Возможные астрофизические проявления существования скалярных зарядов. Труды X Российской гравитационной конференции, Владимир, 20-27 июня 1999 года, с. 145.

10. Денисов В.И., Зубрило А.А. Искажение хода световых лучей гравитационным полем массивной скалярной звезды. Препринт НИИЯФ МГУ № 2000 10/614, 8 с.

11. Denisova I.P., Zubrilo A.A. Application of a parametric method for investigation radial geodesies in Vaidya's metric. // Gravitation and Cosmology, 2000. Vol. 22, No. 2.

12. Григорьев В.И., Денисова И.П. Параметрический метод интегрирования уравнений геодезического движения в пространстве Вайдья. Препринт НИИЯФ МГУ № 2000 11/615. 8 с.

13. Denisova I.P., Dalai М. Development of the method of potentials for the problems of gravitation-electromagnetic conversion. // Journ. of Math. Phys., 1997, 38, p. 5820-5831.

14. Денисов В.И., Зубрило А.А., Кравцов Н.В., Пинчук В.Б. Об использовании кольцевых лазеров для измерения релятивистских эффектов. // Квантовая электроника, 1999, 26, № 2 с. 171-174.

15. Everitt C.W.F. Gravity Probe В. In Proceedings of the 6-th Marcel Grossmann Meeting on General Relativity, ed. Humitaka Sato, Singapore, World Scientific Publishing Co., 1992, p. 1632.

16. Спасский Б.И. История физики. М.: Высшая школа, 1977, ч. 1, с. 135, с. 177.

17. Альберт Эйнштейн и теория гравитации. М.: Мир, 1979. - 529 с.

18. Физический энциклопедический словарь. М.: Советская энциклопедия, 1983. - 928 с.

19. Зельдович Я.Б., Новиков И.Д. Теория тяготения и эволюция звезд. М.: Наука, 1971. - 484 с.

20. Блохинцев Д.И. Пространство-время в микромире. М.: Наука, 1970. - 359 с.

21. Экспериментальные тесты теории гравитации, под ред. В.Б. Брагинского, В.И. Денисова, М.: МГУ, 1989 г. - 255 с.

22. Петров А.З. Пространства Эйнштейна. М.: Физматгиз, 1961. -363 с.

23. Петров А.З. Новые методы в общей теории относительности. М.: Наука, 1966. - 495 с.

24. Freund P.G.O., Maheshwari A., Schömberg Е. Finite-Range gravitation. // Astrophysical Journal, 1969, vol. 157, p. 857-867.

25. M.Visser. Mass for graviton. // General Relativity and Gravitation, 1998, v. 30, № 12, p. 1717-1728.

26. Эйнштейн А. Собрание научных трудов, т.1. M.: Наука, 1965. -700 с.

27. П.A.M. Дирак. Общая теория относительности. М.: Атомиздат, 1978. - 65 с.

28. Ландау Л.Д., Лифшиц Е.М. Теория поля. М.: Наука, 1988, - 510 с.

29. Тихонов А.Н., Самарский A.A. Уравнения математической физики.- М.: Наука, 1972. 735 с.

30. Кошляков Н.С. Уравнения в частных производных математической физики. М.: Высшая школа, 1970, с. 488-534.

31. Рихтмайер Р. Принципы современной математической физики, т.1.- М.: Мир, 1982. 488 с.

32. Захаров В.Д. Гравитационные волны в теории тяготения Эйнштейна. М.: Наука, 1972. - 198 с.

33. Денисов В.И., Мехта Б.В. Космологическая модель скалярно тензорной теории гравитации. // Вестник Московского Университета. Сер. 3. Физ., астр. 1996, № 3, с. 17-22.

34. Bronnikov R.A. Scalar-tensor theory and scalar charge. // Acta Physica Polon, 1973, B4, p. 251-266.

35. Саакян Г.С. Пространство-время и гравитация. Ереван, Из-во Ереванского университета, 1985. - 334 с.

36. Turner E.L. Statistics of the Hubble diagnostic. // Astrophys. J. 1979, 230, p. 291.

37. Vaidya P.C. Non-static gravitation field with spherical symmetry. Proceedings of the Indian Academy of Science, 1948, vol. 14, № 1, p. 53-54.

38. Фишер И.З. Скалярное мезостатическое поле с учетом гравитационных эффектов. // Журнал экспериментальной и теоретической физики, 1948, № 18, р. 636-640.

39. Бочарова Н.М., Бронников К.А., Мельников В.Н. Об одном точном решении уравнений Эйнштейна и скалярного поля. // Вестник МГУ, Физ., Астрономия, 1970, № 6, с. 706-709.

40. Фок В.А. Теория пространства, времени и тяготения. М.: Физ-матгиз, 1961. - 504 с.

41. Арнольд В.И. Математические методы классической механики. -М.: Наука, 1979 . 431 с.

42. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1971. - 1108 с.

43. Денисов В.И. Исследование свойств решения релятивистской теории гравитации в окрестности сингулярной сферы. // Теоретическая и математическая физика, 1997, т. 111, № 1, с. 144-148.

44. Денисов В.И. Изотропные геодезические в окрестности сингулярного статического сферически симметричного решения. // Теоретическая и математическая физика, 1997, т.111, № 2, с.312-320.

45. Блиох П.В., Минаков А.А. Гравитационные линзы. Киев.: Науко-ва думка, 1989, - 240 с.

46. Наумов Н.Д. Об операторном методе расчета движения частиц. // Космические исследования, 1993, т.31, № 4, с. 117-119.

47. Наумов Н.Д. Автомодельные движения заряженных частиц. // Физика плазмы, 1993, т.19, № 11, с. 1406-1408.

48. Bondi Н. Negative mass in General Theory of Relativity. // Rev. Mod. Phys., 1957, Vol. 29, p. 423.

49. Piran T. On Gravitational Repulsion. // General Relativity and Gravitation, 1997, v.29, 11, p.1363-1370.

50. Чандрасекар С. Математическая теория черных дыр. М.: Мир, 1986. 276 с.

51. Alcock С et al. Possible Gravitational Microlensing of a Star in the Large Magellanic Cloud. // Nature, 1993, Vol. 265 p. 621-622.

52. Aubourg E et al. Evidence for Gravitational Microlensing by Dark Objects in the Galactic Halo. // Nature, 1993, Vol. 265, p. 623-624.

53. Рихтмайер P. Принципы современной математической физики, т.2. М.: Мир, 1984. - 383 с.

54. Крамер Д. и другие. Точные решения уравнений Эйнштейна, М.: Энергоиздат, 1982, 416 с.

55. Григорьев В.И., Денисова И.П. Некоторые новые соотношения тензорной алгебры. // Вестник Московского университета, сер.З, физика, астрономия, 1996, №6, с.3-8.

56. Григорьев В.И., Денисова И.П., Некоторые новые соотношения для тензора электромагнитного поля в псевдоримановом пространствевремени. // Вестник Московского Университета, сер.З, физика, астрономия 1997, № 5, с. 15-17.

57. Денисова И.П. Нахождение частного решения уравнений релятивистской теории гравитации методом неопределенных координат. // Теоретическая и математическая физика, 1997, т.112, № 3, с. 501-512.

58. Денисова И.П. Метрика плоской эллиптически поляризованной электромагнитной волны. // ДАН, 1998, 360, с. 335-336.

59. Irene P. Denisova and Binita V. Mehta, Tensor Expressions for Solving Einstein's Equations by the Method of Sequential Approximation. // General Relativity and Gravitation, 1997, vol. 29, c. 583-590.

60. Боголюбов H.H. Избранные труды в трех томах. К.: Наукова думка, 1969-1971. т.1 - 644 е., т.2 - 522 е., т.З - 487 с.

61. Боголюбов H.H. Колебания. Механика в СССР за 30 лет. М.: Гостехиздат, 1950. с. 99-114.

62. Боголюбов H.H., Митропольский Ю.А., Самойленко A.M. Метод ускоренной сходимости в нелинейной механике. К.: Наукова думка, 1969. - 247 с.

63. Боголюбов H.H., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М.: Наука, 1974. - 504 с.

64. Lindquist R.W., Schwartz R.A., Misner C.W. Vaidya's Radiating Schwarzschild Metric. // Phys. Rev., 1965, vol. 137B, p. 1364-1368.

65. Денисов В.И. Введение в электродинамику материальных сред М.: МГУ, 1989. - 168 с.

66. Einstein A. Die Grungladen der allgemeinen Relativitätstheorie. Leipzig. // Annalen der physilc, 1916, 49, s. 769.

67. Румянцев C.B. Соотношение электродинамики Борна Инфельда и электродинамики Максвелла в Общей теории относительности. Всб.: "Проблемы теории гравитации, элементарных частиц", вып. 9, М.: Атомиздат, 1978, с. 66-68.

68. Брагинский В.Б., Панов В.И. Проверка эквивалентности инертной и гравитационой масс. //Журнал экспериментальной и теоретической физики, 1971, № 61, с. 873-879.

69. Дубровин Б.А., Новиков С.П., Фоменко А.Г. Современная геометрия. М.: Наука, 1979. - 759 с.

70. D. Bardas et al. Gravity Probe В. Hardware development progress towards the flight instrument. In Proceedings of the 6-th Marcel Grossmann Meeting on General Relativity, ed. Humitaka Sato, Singapore, World Scientific Publishing Co., 1992, p. 382-393.

71. Y.M. Xiao, et al. Gravity Probe B. The precision gyroscope. In Proceedings of the 6-th Marcel Grossmann Meeting on General Relativity, ed. Humitaka Sato, Singapore, World Scientific Publishing Co., 1992, p. 394-396.

72. Stedman G.E., Li Z., Bilger H.R. Sideband analysis and seismic detection in large ring lasers. // Appl. Opt., 1995, 34, p. 7390-7396.

73. Stedman G.E., Johnsson M.T., Li Z., Rowe C.H., Bilger H.R. Tviolation and microhertz resolution in a ring laser. // Opt. Lett., 1995, 20, p. 324-326.

74. Stedman G.E., Li Z., Rowe C.H., McGregor A.D., Bilger H.R. Harmonic analysis in a large ring laser with backscatter-induced pulling. // Phys. Rev., 1995, A51, p. 4944-4958.

75. Stedman G.E., Bilger H.R., Li Ziyuan, Poulton M.P., Rowe C.H., Vetharaniam I. and Wells P.V. Canterbury ring laser and tests for nonrecip-rocal phenomena. // Austr.J. Phys., 1993, 46, p. 87-101.

76. Denisov V.I. The usin'g of laser gyroscope for the measurement of NUT's parameter. Preprint of Institute of Nuclear Physics MSU, № 98-32/533, 8 c.

77. Денисов В.И., Денисов M.И., Кравцов H.В., Пинчук В.Б. Проверка основных принципов теории гравитации с помощью лазерного гироскопа в Земных условиях. Препринт НИИЯФ МГУ № 96 29 /436, 8 с.