В геофизике модель геопотенциала представляет собой теоретический анализ измерения и расчета эффекты Земли гравитационного поля.
Содержание
- 1 Закон Ньютона
- 1.1 Случай однородной сферы
- 2 Отклонения гравитационного поля Земли от гравитационного поля однородная сфера
- 3 Рекурсивные алгоритмы, используемые для численного распространения орбит космических аппаратов
- 4 Сферические гармоники
- 5 Ссылки
- 6 Внешние ссылки
Закон Ньютона
Диаграмма двух массы, притягивающие друг друга
Закон всемирного тяготения Ньютона утверждает, что гравитационная сила F, действующая между двумя точечными массами m1и m 2 с центром масс разделение r определяется выражением
где G - гравитационная постоянная и r̂ - радиальный единичный вектор. Для объекта с непрерывным распределением массы каждый элемент массы dm можно рассматривать как точечную массу, поэтому интеграл объема по протяженности объекта дает:
| | (1) |
с соответствующим гравитационным потенциалом
| | (2) |
где ρ = ρ (x, y, z) - массовая плотность в элементе объема и направление от элемента объема к точечной массе.
Случай однородной сферы
В частном случае сферы со сферически симметричной плотностью массы, тогда ρ = ρ (s), т.е. плотность зависит только от радиального расстояния
Эти интегралы можно вычислить аналитически. Это теорема оболочек, согласно которой в данном случае:
| | (3) |
с соответствующим потенциалом
| | (4) |
где M = ∫ V ρ (s) dxdydz - полная масса сферы.
Отклонения гравитационного поля Земли от поля однородной сферы
На самом деле Земля не совсем сферическая, в основном из-за ее вращения вокруг полярной оси, что делает ее форму слегка сжатой. Если бы эта форма была точно известна вместе с точной плотностью массы ρ = ρ (x, y, z), интегралы (1) и (2) можно было бы оценить численными методами, чтобы найти более точную модель для Земли. гравитационное поле. Однако на самом деле ситуация обратная. Наблюдая за орбитами космического корабля и Луны, гравитационное поле Земли может быть определено довольно точно, и наилучшая оценка массы Земли получается путем деления произведения GM, определенного в результате анализа орбиты космического корабля, на значение G, определенное на более низкую относительную величину. точность при использовании других физических методов.
Из определяющих уравнений (1) и (2) ясно (взяв частные производные от подынтегрального выражения), что вне тела в пустом пространстве следующие дифференциальные уравнения действительны для поля, вызванного по телу:
| | (5) |
| | (6) |
Функции формы , где (r, θ, φ) - сферические координаты, которые удовлетворяют уравнению в частных производных (6) (уравнение Лапласа ), называются сферической гармоникой. функции.
Они принимают вид:
| | (7) |
где сферические координаты (r, θ, φ), указанные здесь в декартовых (x, y, z) для справок:
| | (8) |
также P n - это многочлены Лежандра, а P n для 1 ≤ m ≤ n - связанные функции Лежандра.
Первые сферические гармоники с n = 0,1,2,3 представлены в таблице ниже.
n | Сферические гармоники |
---|
0 | |
1 | |
|
|
2 | |
|
|
|
|
3 | |
|
|
|
|
|
|
Чт Модель гравитационного потенциала Земли представляет собой сумму
| | (9) |
где , а координаты (8) являются относительно стандартной геодезической системы отсчета расширен в пространстве с началом в центре эллипсоида и оси в направлении полярная ось.
зональные термины относятся к терминам формы:
и тессеральные термины термины относятся к терминам формы:
Зональные и тессеральные термины для n = 1 опущены в (9). Коэффициенты при n = 1 с членами как m = 0, так и m = 1 соответствуют произвольно ориентированному дипольному члену в многополюсном разложении. Гравитация физически не проявляет дипольного характера, поэтому интеграл, характеризующий n = 1, должен быть равен нулю.
различным коэффициентам J n, C n, S n затем присваиваются значения, для которых наилучшее возможное совпадение между вычисленными и получены наблюдаемые орбиты космических аппаратов.
As P n (x) = −P n (−x) ненулевые коэффициенты J n для нечетного n соответствуют отсутствию симметрии «север-юг» относительно экваториальной плоскости для распределения массы Земли. Ненулевые коэффициенты C n, S n соответствуют отсутствию вращательной симметрии вокруг полярной оси для распределения массы Земли, то есть «трехосности» Земли.
Для больших значений n приведенные выше коэффициенты (которые делятся на r в (9)) принимают очень большие значения, когда, например, в качестве единиц используются километры и секунды. В литературе принято вводить произвольный «опорный радиус» R, близкий к радиусу Земли, и работать с безразмерными коэффициентами
и записать потенциал как
| | (10) |
Главный член (после члена −μ / r) in (9) - это «член J 2 »:
Относительно системы координат
| | (11) |
Рисунок 1: Единичные векторы. Это не правильно. Должна быть тета, а не лямбда
проиллюстрированные на рисунке 1 компоненты силы, вызванной «членом J 2 », равны
| | (12) |
В прямоугольной системе координат (x, y, z) с единичными векторами (x̂ ŷ ẑ) компоненты силы:
| | (13) |
Составляющие силы, соответствующие «члену J 3 »
равны
| | (14) |
и
| | (15) |
Точные численные значения коэффициентов различаются (несколько) между разными моделями Земли, но для самых низких коэффициентов все они почти полностью совпадают.
Для JGM-3 значения следующие:
- μ = 398600,440 км⋅с
- J2= 1,75553 × 10 км⋅с
- J3= −2,61913 × 10 км⋅с
. Например, на радиусе 6600 км (около 200 км над поверхностью Земли) J 3 / (J 2 r) составляет около 0,002, то есть поправка к «J 2 сила »из« члена J 3 »составляет порядка 2 промилле. Отрицательное значение J 3 подразумевает, что для точечной массы в экваториальной плоскости Земли гравитационная сила слегка наклонена к югу из-за отсутствия симметрии массового распределения Земли «север-юг».
Рекурсивные алгоритмы, используемые для численного распространения орбит космических аппаратов
Орбиты космических аппаратов вычисляются посредством численного интегрирования уравнения движения. Для этого необходимо вычислить гравитационную силу, т.е. градиент потенциала. Эффективные рекурсивные алгоритмы были разработаны для вычисления силы тяжести для любых и (максимальная степень зональных и тессеральных условий) и такие алгоритмы используются в стандартном программном обеспечении для распространения орбиты.
Доступные модели
Самыми ранними моделями Земли, обычно используемыми NASA и ESRO / ESA, были «Земля Годдарда. Модели », разработанные Центром космических полетов им. Годдарда, обозначались« GEM-1 »,« GEM-2 »,« GEM-3 »и так далее. Позже стали доступны «Совместные модели земной гравитации», обозначенные «JGM-1», «JGM-2», «JGM-3», разработанные Центром космических полетов Годдарда в сотрудничестве с университетами и частными компаниями. В более новых моделях обычно использовались термины более высокого порядка, чем в их предшественниках. EGM96 использует N z = N t = 360, что дает 130317 коэффициентов. Также доступна модель EGM2008.
Для обычного спутника Земли, требующего точности определения / предсказания орбиты в несколько метров, "JGM-3" усечено до N z = N t = 36 ( 1365 коэффициентов) обычно достаточно. Неточности моделирования сопротивления воздуха и, в меньшей степени, давления солнечной радиации будут превышать неточности, вызванные ошибками моделирования гравитации.
Безразмерные коэффициенты , , для первых зональных и тессеральных терминов (с использованием = 6378,1363 км и = 398600,4415 км / с) модели JGM-3 равны
Зональные коэффициентыn |
---|
2 | -0,1082635854D-02 |
---|
3 | 0,2532435346D-05 |
---|
4 | 0,1619331205D-05 |
---|
5 | 0,2277161016D-06 |
---|
6 | -0,5396484906D-06 |
---|
7 | 0,3513684422D-06 |
---|
8 | 0.2025187152D-06 |
---|
Тессеральные коэффициентыn | m | C | S |
---|
2 | 1 | - 0.3504890360D-09 | 0.1635406077D-08 |
---|
2 | 2 | 0.1574536043D-05 | -0.9038680729D-06 |
---|
3 | 1 | 0.2192798802D-05 | 0.2680118938D-06 |
---|
3 | 2 | 0,3090160446D-06 | -0,2114023978D-06 |
---|
3 | 3 | 0.1005588574D-06 | 0,1972013239D-06 |
---|
4 | 1 | -0,5087253036D-06 | - 0.4494599352D-06 |
---|
4 | 2 | 0.7841223074D-07 | 0.1481554569D-06 |
---|
4 | 3 | 0.5921574319D-07 | -0.1201129183D-07 |
---|
4 | 4 | -0.3982395740D-08 | 0,6525605810D-08 |
---|
Следовательно, согласно JGM-3, км / с = км / с и км / с = км / с
Сферические гармоники
Ниже приводится краткое описание сферических гармоник, используемых для моделирования гравитационного поля Земли. Сферические гармоники являются производными от поиска гармонических функций вида
| | (16) |
где (r, θ, φ) - это сферические координаты, определенные уравнениями (8). Непосредственными вычислениями получаем, что для любой функции f
| | (17) |
Представляем выражение (16) в (17) получаем, что
| | (18) |
В качестве члена
зависит только от переменной и сумма
зависит только от переменных θ и φ. Получается, что φ является гармоническим тогда и только тогда, когда
| | (19) |
и
| | (20) |
для некоторой константы
From ( 20) тогда следует, что
Первые два члена зависят только от переменной , а третий - только от переменной .
Из определения φ как сферической координаты ясно, что Φ (φ) должна быть периодической с периодом 2π, и поэтому должно быть, что
| | (21) |
и
| | (22) |
для некоторого целого числа m в качестве семейства решений (21), тогда равны
| | (23) |
При подстановке переменной
уравнение (22) принимает форму
| | (24) |
Из (19) следует, что в чтобы получить решение с
необходимо иметь это
Если P n (x) является решением дифференциального уравнения
| | (25) |
следовательно, потенциал, соответствующий m = 0
, который является вращательно-симметричным относительно оси z, является гармонической функцией
Если - решение дифференциального уравнения
| | (26) |
с m ≥ 1 имеет потенциал
| | (27) |
, где a и b - произвольные константы, является гармонической функцией, которая зависит от φ и поэтому не осесимметрична относительно оси z
Дифференциальное уравнение (25) имеет вид дифференциальное уравнение Лежандра, для которого полиномы Лежандра определяют
| | (28) |
- решения.
Произвольный коэффициент 1 / (2n!) Выбирается так, чтобы P n (−1) = - 1 и P n (1) = 1 для нечетных n и P n (−1) = P n (1) = 1 для четного n.
Первые шесть полиномов Лежандра:
| | (29) |
Решения дифференциального уравнения (26) являются ассоциированными функциями Лежандра
| | (30) |
Следовательно,
Ссылки
- Эль Ясберг Теория полета искусственных спутников Земли, Израильская программа научных переводов (1967)
- Лерх, Ф.Дж., Вагнер, Калифорния, Смит, Д.Э., Сэндсон, М.Л., Браунд, Дж. Э., Ричардсон, Дж. А., "Гравитационное поле" Модели Земли (GEM1 и 2) », Отчет X55372146, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1972
- Лерх, Ф.Дж., Вагнер, Калифорния, Патни, М.Л., Сандсон, М.Л., Браунд, Дж. Э., Ричардсон, Дж. А., Тейлор, Вашингтон, "Модели гравитационного поля GEM3 и 4", Отчет X59272476, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1972
- Лерх, Ф.Дж., Вагнер, Калифорния, Ричардсон, Дж. А., Браунд, Дж. Э., «Модели Земли Годдарда (5 и 6)», Отчет X92174145, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1974
- Лерх, Ф.Дж., Вагнер, Калифорния, Клоско, С.М., Белотт, Р.П., Лаубшер, Р. Рейлор, WA, «Улучшение гравитационной модели с использованием Geos3 Altimetry (GEM10A и 10B)», Весеннее ежегодное собрание Американского геофизического союза 1978 года, Майами, 1978
- Лерх, Ф.Дж., Клоско, С.М., Ла Убшер, Р.Е., Вагнер, К.А., "Улучшение гравитационной модели с использованием Geos3 (GEM9 и 10)", Журнал геофизических исследований, том. 84, В8, с. 3897-3916, 1979
- Lerch, F.J., Putney, B.H., Wagner, C.A., Klosko, S.M., "Модели Земли Годдарда для океанографических приложений (GEM 10B и 10C)", Морская геодезия, 5 (2), с. 145-187, 1981
- Лерх, Ф.Дж., Клоско, С.М., Патель, Великобритания, «Уточненная гравитационная модель из Лагеоса (GEML2)», 'Технический меморандум НАСА 84986, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1983
- Лерх, Ф.Дж., Нерем, Р.С., Патни, Б.Х., Фелсентрегер, Т.Л., Санчес, Б.В., Клоско, С.М., Патель, Великобритания, Уильямсон, Р.Г., Чинн, Д.С., Чан, Дж.К., Рахлин, К.Э., Чендлер, Н.Л., Маккарти, Дж. Дж., Маршалл, Дж. А., Латке, С.Б., Павлис, Д.В., Роббинс, Дж. У., Капур, С., Павлис, Е.К., "Геопотенциальные модели Земли по данным спутникового слежения, высотомера и наблюдений за поверхностной гравитацией: GEMT3 и GEMT3S », Технический меморандум НАСА 104555, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1992
- Лерх, Ф.Дж., Нерем, Р.С., Патни, Б.Н., Фелсентрегер, Т.Л., Санчес, Б.В., Маршалл, Дж. Клоско, С.М., Патель, ГБ, Уильямсон, Р.Г., Чинн, Д.С., Чан, Дж. К., Рахлин, К. Э., Чендлер, Н. Л., Маккарти, Дж. Дж., Лутке, С. Б., Павлис, Н. С., Павлис, ЭЦ, "А" Геопотенциальная модель на основе данных спутникового слежения, высотомера и данных поверхностной силы тяжести: GEMT3 ", Journal of Geophysical Research, Vol. 99, № B2, стр. 2815-2839, 1994
- Нерем, Р.С., Лерх, Ф.Дж., Маршалл, Дж. А., Павлис, Е.К., Патни, Б.Х., Тэпли, Б.Д., Эанс, Р.Дж., Райс, Дж., Уоткинс, М.М., Клоско, С.М., Чан, Дж. К., Лутке, С.Б., Патель, Г.Б., Павлис, Н.К., Уильямсон, Р.Г., Рапп, Р.Х., Бианкал, Р., Ноуэль, Ф., "Разработка модели гравитации для Topex / Посейдон: Совместные гравитационные модели 1 и 2 », Журнал геофизических исследований, Vol. 99, No. C12, p. 24421-24447, 1994a
Внешние ссылки