Я сделал наивный подсчет высоты экваториальной выпуклости Земли и обнаружил, что она должна быть около 10 км. Реальная высота около 20 км. У меня вопрос: почему такое несоответствие?
Расчет, который я сделал, состоял в том, чтобы представить, как я помещаю мяч на вращающуюся Землю. Куда бы я его ни поместил, он не должен двигаться.
Гравитационный потенциал на единицу массы шара равен , с высота над расстоянием от полюса до центра Земли (назовем это ) и гравитационное ускорение.
Гравитация хочет притянуть мяч к полюсам, подальше от выпуклости. Он уравновешивается центробежной силой, которая имеет потенциал на единицу массы, с угловая скорость Земли и угол от северного полюса. Это происходит из того, что в инерциальной системе отсчета было бы кинетической энергией мяча, и превращения ее в потенциал в ускоряющей системе отсчета.
Если мяч не движется, этот потенциал должен быть постоянным, поэтому
мы могли бы также положить константу равной нулю и написать
Для Земли,
Это дает 10,8 км, когда , поэтому экваториальная выпуклость должна быть примерно такой большой.
Согласно Википедии , диаметр Земли на экваторе на 42,72 км шире, чем между полюсами, а это означает, что выпуклость примерно в два раза больше, чем я ожидал. (Википедия указывает диаметр; я оценил радиус.)
Откуда берется лишняя выпуклость? Мой простой расчет использует и как константы, но ни один из них не изменяется более чем на процент или около того. Это правда, что Земля не имеет однородной плотности, но мне непонятно, как это должно повлиять на расчет, если распределение плотности по-прежнему сферически-симметрично (или почти так).
(В Википедии также есть выражение без производных, которое согласуется с моим.)
Ошибка заключается в том, что вы предполагаете, что распределение плотности «почти сферически симметрично». Это достаточно далеко от сферической симметрии, если вы хотите рассчитать субведущие эффекты первого порядка, такие как экваториальная выпуклость. Если ваша цель состоит в том, чтобы вычислить отклонения уровня моря от сферической симметрии (до первого порядка), непоследовательно пренебрегать такими же большими поправками первого порядка к сферической симметрии с другой стороны — источника гравитации. Другими словами, термин в вашем потенциале неправильно.
Только представьте, что Земля представляет собой эллипсоид с экваториальной выпуклостью, она не вращается и на поверхности нет воды. Каким был бы потенциал на поверхности или потенциал на фиксированном расстоянии от центра эллипсоида? Вы де-факто предполагали, что в этом случае будет где - фиксированный радиус Земли (сферического распределения вещества) и - фактическое расстояние зонда от начала координат (центра Земли). Однако этим анзацем вы только признали переменное расстояние зонда от сферически симметричного источника гравитации: вы все еще пренебрегли вкладом выпуклости в несферичность гравитационного поля.
Если включить несферически-симметричную поправку к гравитационному полю Земли, примерно изменится на , и, соответственно, искомая выпуклость придется удвоить, чтобы компенсировать вращательный потенциал. Эвристическое объяснение фактора заключается в том, что истинный потенциал над эллипсоидом зависит от «нечто промежуточного» между расстоянием от центра масс и расстоянием от поверхности. Другими словами, «постоянная потенциальная поверхность» вокруг эллипсоидального источника материи находится «точно посередине» между фактической поверхностью эллипсоида и сферической поверхностью. поверхность.
Я постараюсь добавить более точные формулы для гравитационного поля эллипсоида в обновленной версии этого ответа.
Обновление: гравитационное поле эллипсоида
Я численно проверил, что гравитационное поле эллипсоида имеет именно тот эффект деления пополам, который я набросал выше, используя код Монте-Карло Mathematica, чтобы избежать двойных интегралов, которые можно было бы вычислить аналитически, но пока это меня раздражает.
Я взял миллионы случайных точек внутри вытянутого эллипсоида с «радиусами». ; обратите внимание, что разница между двумя радиусами . Среднее значение , обратное расстояние между случайной точкой эллипсоида и выбранной точкой над эллипсоидом равно меньше, если выбранная точка находится над экватором, чем если она находится над полюсом, при условии, что расстояние от начала координат одинаково для обеих выбранных точек.
Код:
{xt, yt, zt} = {1.1, 0, 0};
runs = 200000;
totalRinverse = 0;
total = 0;
For[i = 1, i < runs, i++,
x = RandomReal[]*2 - 1;
y = RandomReal[]*2 - 1;
z = RandomReal[]*2 - 1;
inside = x^2/0.81 + y^2/0.81 + z^2 < 1;
total = If[inside, total + 1, total];
totalRinverse =
totalRinverse +
If[inside, 1/Sqrt[(x - xt)^2 + (y - yt)^2 + (z - zt)^2], 0];
]
res1 = N[total/runs / (4 Pi/3/8)]
res2 = N[totalRinverse/runs / (4 Pi/3/8)]
res2/res1
Описание
Используйте приведенный выше код Mathematica: его цель — вычислить одну чисто числовую константу из-за пропорциональности несферичности гравитационного поля выпуклости; масса; Постоянная Ньютона очевидна. Последнее число, которое печатает код, является средним значением . Если в начале выбрано {1.1, 0, 0} вместо {0, 0, 1.1}, программа генерирует 0,89 вместо 0,94. Это доказывает, что гравитационный потенциал эллипсоида ведет себя как на расстоянии от происхождения, где — локальная высота поверхности относительно идеализированной сферической поверхности.
В приведенном выше коде я выбрал эллипсоид с радиусами (0,9, 0,9, 1), который представляет собой вытянутый сфероид (длинный, палообразный), в отличие от Земли, которая близка к сплюснутому сфероиду (плоскому, дискообразному). Так что пусть вас не смущают некоторые знаки — они работают нормально.
Бонус от Исаака
Мариано С. указал на следующее решение довольно известного автора:
И) Плоскостность. Здесь мы хотели бы аналитически вычислить решение Любоша Мотля в первом порядке по параметру плоскостности ,
где и - экваториальный и полярный радиусы Земли соответственно, и . ( отныне будет означать равенство до членов более высокого порядка в .) Мы предполагаем, что Земля представляет собой массивный сплюснутый эллипсоид
с равномерной плотностью и объем
Эксцентриситет
II) квадруполь. Будем считать, что гравитация ньютоновская. Что нам нужно вычислить, так это вклад квадрупольного момента в гравитационный потенциал
По соображениям симметрии одно из главных направлений квадрупольного момента должно быть вдоль полярного -ось, а два других главных направления должны иметь одинаковые собственные значения и лежать в экваториальной плоскости. самолет. Таким образом должен быть в форме
где и – экваториальное и полярное собственные значения соответственно; где мы ввели сокращенное обозначение и ; и где является полярным углом. Поскольку квадрупольный момент не может внести вклад в закон Гаусса, мы должны потребовать, чтобы лапласиан обращается в нуль, что приводит к тому, что полярное собственное значение в два раза меньше экваториального собственного значения, . Другими словами,
Поэтому достаточно вычислить гравитационный потенциал в точке на полярном -ось, где у нас есть азимутальная симметрия .
III) Кольцо. Вычислим вклад в гравитационную потенциальную энергию кольца, параллельного плоскость и с центром на -ось. Из-за азимутальной симметрии, мы можем сфокусироваться на точке на кольце с и , и которые лежат в самолет. Пусть эта точка имеет двумерные полярные координаты . Другими словами, точка имеет трехмерные координаты. . Из уравнения эллипса
после некоторых элементарных вычислений получаем
Далее, расстояние с кольца в точку дается косинусным соотношением
Площадь «внешней» поверхности кольца
с бесконечно малой «радиальной» толщиной
Таким образом, объем кольца
(В конце концов мы хотим интегрировать по полярному углу от к . Это соответствует интегрированию по от к в отрицательном направлении. Поэтому является отрицательным.)
IV) Потенциал. Вклад кольца в гравитационную потенциальную энергию в точку является
Интегрирование по косинусу полярного угла для формирования тонкой оболочки выходов (с помощью программы MAPLE)
Далее мы интегрируем по чтобы получить полярный потенциал массивного эллипсоида
V) Монополь и квадруполь. Первый срок есть, конечно, именно монопольный потенциал
и второй срок - полярный квадрупольный потенциал
Таким образом, мы знаем, что гравитационный квадрупольный потенциал в произвольной точке (не обязательно на -ось) есть
а полный гравитационный потенциал равен
VI) Поверхность. С этого момента будем рассматривать только точки с
на поверхности эллипсовидной Земли. Тогда монопольный потенциал в становится
а квадрупольный потенциал в становится
так что полный гравитационный потенциал в становится
VII) Обсуждение. Возьмем за точку отсчета северный полюс, т. е. вычтем гравитационный потенциал на Северном полюсе. Тогда гравитационный потенциал в становится
где
- разница экваториального и полярного радиусов, а
Теперь мы должны добавить центробежный потенциал
К заказу в что мы работаем, мы видим, что общий потенциал постоянна (не зависит от точки поверхности ), если
Вывод: находим фактор в отличие от исходного монопольного аргумента Марка Эйхенлауба.
Обновление: через полгода после того, как этот ответ был опубликован на Phys.SE, 18 ноября 2011 г., страница Википедии изменила указанное математическое выражение для параметра плоскостности на , и теперь полностью согласен с этим ответом.
Были некоторые сомнения по поводу ответа Любоша (который я принял), так что это просто проверка.
Я скопировал метод, описанный Любошем, и нашел разность потенциалов для эллипсоида с разными эксцентриситетами. Конечно же, для сплюснутого сфероида, если вы сделаете расстояние от центра до экватора долей больше, чем расстояние между центром и полюсом, потенциал составляет примерно долю на экваторе меньше.
Чтобы решить всю проблему, мы должны были бы принять во внимание различную плотность Земли, но в качестве грубой оценки это, похоже, работает.
Например, для сплюснутого сфероида
среднее значение в составляет около 0,996, а в среднем на составляет около 0,991.
Код Python ниже (извините за дилетантство)
import random
import math
points = 10000000
e = .01
rad = 1+e
diam = 2*rad
pot= 0
count = 0
for i in range(1,points):
x = diam*random.random()-rad
y = diam*random.random()-rad
z = diam*random.random()-rad
r = math.sqrt((x-rad)*(x-rad)+y*y+z*z)
if x*x/(rad*rad)+y*y/(rad*rad) + z*z < 1:
pot = pot + 1/r
count = count + 1
print pot/count
pot2 = 0
count = 0
for j in range(1,points):
x = diam*random.random()-rad
y = diam*random.random()-rad
z = diam*random.random()-rad
r = math.sqrt(x*x+y*y+(z-1.0)*(z-1.0))
if x*x/(rad*rad)+y*y/(rad*rad) + z*z < 1:
pot2 = pot2 + 1/r
count = count + 1
print pot2/count
Начните с невозмущенного гравитационного потенциала для однородной сферы массы M и радиуса R, внутренней и внешней:
Добавьте квадрупольное возмущение, вы получите
Масштабные коэффициенты M и R нужны только для того, чтобы безразмерный, спад просто так, что внешнее решение решает уравнение Лапласа, а согласование решений должно гарантировать, что на любом эллипсоиде около сферы радиуса R два решения равны по порядку . Причина, по которой это работает, заключается в том, что решения совпадают как по значению, так и по первой производной при x = R, поэтому они остаются согласованными по значению в ведущем порядке, даже если они отклонены от сферы. Приказ квадрупольные члены равны на сфере и, следовательно, соответствуют ведущему порядку.
Эллипсоид, который я выберу, решает уравнение:
Z-диаметр увеличивается на долю , а диаметр x уменьшился на . Так что отношение полярного радиуса к экваториальному равно . К ведущему заказу
Мы уже сопоставили значения внутреннего и внешнего решений, но нам нужно сопоставить производные. взяв букву «д»:
На первый заказ в , только первый член второго уравнения изменяется из-за того, что r не является постоянным на эллипсоиде. Специализируясь на поверхности эллипсоида:
Приравнивая входную и выходную производные, части, пропорциональные отменить (как и должно быть — тангенциальные производные равны, потому что две функции равны на эллипсоиде). Остальные тоже должны отменить, так что
Таким образом, вы находите связь между и . Решение для дает
Это означает, глядя на уравнение в скобках, что эквипотенциалы на 60% сжаты, как эллипсоид.
Теперь есть условие, что это уравновешивается вращением, а это означает, что эллипсоид является эквипотенциальным, если добавить центробежный потенциал:
Чтобы сделать эквипотенциальность эллипсоида требует, чтобы равно оставшемуся , так что, позвонив (ускорение свободного падения) под названием «g», и по наименованию «С» (центробежный)
Фактическая разница экваториального и полярного диаметров находится путем умножения на 3/2 (см. выше):
вместо наивной оценки . Таким образом, наивная оценка умножается на два с половиной для вращающейся сферы с одинаковой плотностью.
Предыдущее решение является как внутренним, так и внешним для вращающегося однородного эллипсоида, оно точно по r и имеет только главный порядок отклонения от сферической симметрии. Таким образом, он сразу расширяется, чтобы придать форму Земли неравномерному распределению внутренней массы. Оценка с однородной плотностью на удивление хороша, и это потому, что существуют конкурирующие эффекты, в значительной степени нейтрализующие поправку на неоднородную плотность.
Два конкурирующих эффекта: 1. внутреннее распределение более эллиптическое, чем поверхностное, потому что внутреннее решение чувствует, что вся окружающая эллиптическая Земля деформирует его, а дополнительная плотность деформирует его больше. 2. Эллиптичность внутренности подавляется спад квадрупольного решения уравнения Лапласа, который быстрее обычного потенциала. Таким образом, хотя внутренняя часть несколько более деформирована, спад более чем компенсирует эффект внутренней дополнительной плотности, который делает Землю более сферической, хотя и ненамного.
Эти конкурирующие эффекты смещают поправочный коэффициент с 2,5 до 2, что на самом деле довольно мало, учитывая, что внутренняя часть Земли чрезвычайно неоднородна, а центральная часть более чем в три раза плотнее внешних частей.
Точное решение немного сложное, поэтому я начну с дурацкой модели. Это предполагает, что Земля представляет собой однородный эллипсоид с массой M и параметром эллиптичности , плюс точечный источник посередине (или сфера, это не имеет значения), что составляет дополнительную массу внутри, массу M'. Внутренний потенциал задается суперпозицией. С центробежным потенциалом:
Он имеет схематическую форму сферического плюс квадруполь (включая центробежную силу внутри F и G).
Условие, что эллипсоид является эквипотенциалом, находится заменой с внутри F(r) и установив D-часть в ноль:
В этом случае вы получите приведенное ниже уравнение, которое сводится к предыдущему случаю, когда :
где центробежная сила, а сила тяжести на поверхности. Следует отметить, что сферическая часть центробежного потенциала всегда вносит подчиненный член, пропорциональный к уравнению и должны быть опущены. Результат
Таким образом, если вы выберете M' равным 0,2 M, вы получите правильный ответ, так что дополнительный экваториальный радиус будет в два раза больше наивной величины .
Это говорит о том, что потенциал на поверхности Земли изменен только по сравнению с оценкой однородного эллипсоида путем добавления сферы с 20% общей массы в центре. Это несколько мало, учитывая, что неоднородная плотность внутри содержит около 25% массы Земли (возмущающая масса вдвое превышает плотность на половине радиуса, то есть около 25% от общего количества). Небольшая разница связана с эллиптичностью ядра.
Главное, что упущено из вышеизложенного, это то, что центр также несферичен и, таким образом, добавляет к несферичности D часть потенциала на поверхности. Этот эффект в основном противодействует общей тенденции дополнительной массы в центре делать поверхность более сферической, хотя и несовершенной, так что остается коррекция.
Вы можете рассматривать его как суперпозицию однородных эллипсоидов среднего радиуса s с параметром эллиптичности для увеличивается по мере продвижения к центру. Каждый из них однороден внутри, с массовой плотностью где - дополнительная плотность Земли на расстоянии s от центра, так что . Эти эллипсоиды накладываются на эллипсоид однородной плотности с плотностью равна поверхностной плотности земной коры:
Я учту и известно, так что я также знаю , это (отрицательная) производная по s, которая является плотностью эллипсоида, которую вы добавляете к s, и я также знаю:
Количество является раз больше дополнительной массы внутри по сравнению с однородной Землей при плотности земной коры. Обратите внимание, что не зависит от эллипсоидальной формы в ведущем порядке, потому что все вложенные эллипсоиды представляют собой квадрупольные возмущения и, следовательно, содержат тот же объем, что и сферы.
Каждый из этих концентрических эллипсоидов сам по себе является эквипотенциальной поверхностью для центробежного потенциала плюс потенциал внутреннего и внешнего эллипсоидов. Итак, как только вы знаете форму потенциала всех этих наложенных эллипсоидов, которая имеет форму сферического + квадрупольного + центробежного квадруполя (центробежная сферическая часть всегда дает второстепенную поправку, поэтому я ее опускаю):
Вы знаете, что каждый из этих вложенных эллипсоидов является эквипотенциальным
так что уравнение, требующее, чтобы это был эквипотенциал при любом s, имеет вид
Чтобы найти форму F и G, вы сначала выражаете внутреннее/внешнее решение для однородного эллипсоида с точки зрения плотности и радиус R:
Вы можете проверить знак и численное значение коэффициентов, используя правило 3/5 для внутренних эквипотенциальных эллипсоидов, раздельное согласование сферических и D возмущений при r=R и размерный анализ. я положил фактор на дне так что правая часть решает уравнение Лапласа в свободной форме с постоянными.
Теперь вы можете наложить все эллипсоиды, установив на каждом эллипсоиде быть , параметр на каждом эллипсоиде быть , и быть . Я только собираюсь дать внутреннее решение в r (выполняя интегрирование по частям на сферической части, где, как вы знаете, окажется ответ, и отбрасывая некоторую аддитивную константу C):
Первые два члена являются внутренним решением для постоянной плотности. . Третий член — это полный сферический вклад, такой же, как и в сферически-симметричном случае. Четвертый член представляет собой наложенный внешний потенциал эллипсоидов внутри r, а последний член представляет собой наложенный внутренний потенциал эллипсоидов вне r.
Отсюда вы можете прочитать сферическую и квадрупольную части:
Так что интегральное уравнение для утверждает, что форма является эквипотенциальной на любой глубине.
Это уравнение можно решить численно для любого профиля массы внутри, чтобы найти . Это сложно сделать вручную, но вы можете получить качественное представление.
Рассмотрим эллипсоидальное возмущение внутри однородного эллипсоида плотности. Если вы позволите этой массе осесть вдоль эквипотенциала, она примет ту же эллипсоидальную форму, что и поверхность, потому что внутреннее решение для однородного эллипсоида является квадратичным и, следовательно, имеет точно вложенные эллипсоиды той же формы, что и эквипотенциалы. Но эта дополнительная плотность будет давать меньше, чем ее доля эллиптического потенциала поверхности, уменьшаясь как третья степень отношения радиуса Земли к радиусу возмущения. Но это создаст более сильные эллипсы внутри, так что внутренняя часть всегда будет более эллиптической, чем поверхность.
Точное решение слишком сложно для бумажных и карандашных расчетов, но ищите [здесь]( http://www.google.com/imgres?hl=en&client=ubuntu&hs=dhf&sa=X&channel=fs&tbm=isch&prmd=imvns&tbnid=hjMCgNhAjHnRiM:&imgrefurl= http://www.springerimages.com/Images/Geosciences/1-10.1007_978-90-481-8702-7_100-1&docid=ijMBfCAOC1GhEM&imgurl=http://img.springerimages.com/Images/SpringerBooks/BSE%253D5898/BOK %253D978-90-481-8702-7/PRT%253D5/MediaObjects/WATER_978-90-481-8702-7_5_Part_Fig1-100_HTML.jpg&w=300&h=228&ei=ZccgUJCTK8iH6QHEuoHICQ&zoom=1&iact=hc&vpx=210&vpy=153&dur=4872&hovh=182&hovw=240&tx =134&ty=82&sig=108672344460589538944&page=1&tbnh=129&tbnw=170&start=0&ndsp=8&ved=1t:429,r:1,s:0,i:79&biw=729&bih=483 ) ,вы видите, что разумно смоделировать Землю как две концентрические сферы радиуса и с общей массой и и и .
я возьму
и
то есть внутренняя сфера имеет диаметр 3000 км и плотность в два раза больше, что примерно соответствует действительности. Наложение потенциалов и нахождение уравнения для s (двухточечное усечение интегрального уравнения), вы найдете
Где
- гравитационная сила и центробежная сила на единицу массы, как обычно. Использование параметров и определение и , находят:
(это точные десятичные дроби, есть знаменатели 100 и 25). Вычитание двух уравнений дает:
(все еще точные дроби) Что дает уравнение
Так что фактор впереди , вместо наивного 2. Это дает экваториальный диаметр 44,3 км, а не 42,73, что достаточно близко, чтобы модель по существу объяснила все, что вы хотели знать.
Значение также интересно, он говорит вам, что ядро Земли на 9% более эксцентрично, чем внешний эллипсоид самой Земли. Учитывая, что точность модели находится на уровне 3%, это должно быть очень точно.
В этом ответе я представлю основу для использования, а затем включу в нее предыдущие ответы. Позвольте мне суммировать ценности, которые мы здесь имеем. Я буду использовать те же обозначения (насколько это возможно), что и все остальные, и Википедию для сплюснутого сфероида, где - большой, экваториальный, радиус.
Вот мой подход:
Массу Земли можно принять как комбинацию двух фигур: внутренней сферы с радиусом и краевой объем, который представляет собой сплюснутый сфероид за вычетом внутренней сферы. Земля имеет некоторую среднюю плотность, определяемую , но его можно разделить на два разных типа материалов: ядро и кору. Требуемая общая масса будет тогда требовать, чтобы , в то время как средняя плотность Земли требует , ограничение и на одну степень свободы. При написании кода мы можем сказать, что всегда подразумевает в соответствии с , плотность земной коры также может означать плотность ядра. Затем ограничение эквипотенциальной поверхности фиксирует другую степень свободы, которую можно использовать для итеративного нахождения значения . Иллюстрация:
Потенциал из внутренней сферы прост. Я напишу это для точки на экваторе и на полюсе в сочетании с другими потенциалами.
Очевидно, что самое сложное — вычислить потенциал по смехотворно неправильной форме края. Однако перед этим важно подумать о физических последствиях такого взгляда на проблему. Начнем с того, какие плотности относятся к Земле? Здесь указана средняя плотность, а затем плотность в пределах от 0 до 200 км от поверхности.
Когда мы на самом деле решим задачу, мы укажем плотность коры, а затем это будет означать плотность ядра сферы. Является ли этот метод точным? Нет. Главное, что он упускает, это то, что центральная сфера не является сферически-симметричным распределением материи. Теоретически плотности будут расслаиваться в соответствии с постоянными потенциальными линиями. Другими словами, если бы у Земли было ядро высокой плотности, оно тоже было бы сплюснутым сфероидом. По этой причине введение внутренней сферы упускает некоторые детали, но эта модель все еще может быть довольно хорошей.
Реализация этого немного сложна, как указывали другие, но сосредоточение вычислений на краевом объеме очень помогает. Я также использовал метод Монте-Карло, но вычел объем ядра. Другими словами, я выстреливал точками случайным образом в крайний объем настолько эффективно, насколько это было разумно. Для этого я использовал взвешенный метод, и, похоже, он сработал. С 5 миллионами итераций я обнаружил случайные вариации в расчете потенциала при иметь стандартное отклонение около эквивалента гравитационного потенциала и около эквивалента гравитационного потенциала на . Причина более высокого отклонения в связано с тем, что вблизи него масса больше, а выборка была несмещенной в плоскости yz для расчетов потенциал и несмещенный в плоскости xy для расчетов потенциал. В любом случае, когда я повторял значение, я установил допуск для , потому что это должно быть значительно выше, чем статистическая вариация, и этого достаточно для расчета выпуклости. Напомним, это мой метод:
Для расчета потенциала при
Я так и сделал, и для разных значений плотности корки получил следующее.
Я подумал, что это хорошие результаты, потому что первый и последний из них близки к предыдущим ответам в пределах числовой погрешности, а разумное значение плотности земной коры приближается к реальному значению .
Если кому-то интересно, я могу посмотреть код для этого на github или что-то подобное. В противном случае он немного длиннее, чем другие, размещенные здесь, поэтому пока я не буду загромождать это пространство.
Я смотрю на две модели "жирной земли":
Для обоих анализов я использую следующий результат: рассмотрим тонкую оболочку из материала с плотностью , средний радиус и толщина
(Это полярные координаты с , второй полином Лежандра.) Заметим, что , так является уплощение. Легко показать, что средний радиус действительно а чистая масса раковины равна 0. (Подумайте о поверхностной плотности массы , что отрицательно, где . Не обращайте внимания на афизичность отрицательной массы; на практике эта оболочка будет наложена на сферу.)
Тогда в первом порядке по f поле, создаваемое этой оболочкой, равно
(Здесь это объем сферы с радиусом .) Я включаю схему расчета в конце.)
1. сферически-симметричный внутренний слой + поверхностный слой
Пусть средняя плотность внутренней части и плотность на поверхности . Тогда для описанного выше асферического поверхностного слоя потенциал на поверхности равен (в первом порядке по f):
где первый член включает изменение сферического потенциала (« "член) второй - квадрупольный потенциал, а последний - псевдопотенциал вращения (записанный в терминах полиномов Лежандра). Грубо говоря, у полюсов поверхность ближе к центру, а значит, и глубже в потенциале хорошо, но уменьшение близлежащей массы частично нейтрализует этот эффект.
Для гидростатического равновесия этот потенциал не должен зависеть от , так:
Решение:
или
Числа: я использую км, кг/м , и
Случаи.
2. Две наложенные друг на друга сплюснутые сферы, каждая с постоянной плотностью.
Следуя анализу Рона Маймона, теперь я добавляю основную массу. Обозначение:
Предполагается, что обе поверхности сферы находятся в гидростатическом равновесии (и, следовательно, сплющены). Два коэффициента уплощения связаны через свои квадрупольные поля.
Накладывая поля для двух масс, потенциал на поверхности ядра равен:
Потенциал на внешней поверхности равен:
Требование для гидростатического равновесия, чтобы не было -зависимость дает пару уравнений:
Решение:
Обратите внимание, что , в соответствии с предыдущим анализом.
Числа. Используя значения Рона Маймона: , один вычисляет и , в результате чего значение для рассчитано в части 1 случае 4, , увеличивается до 2,002 (+4%), если учитывать сжатие сердцевины. Этот эффект меньше, чем эффект массы ядра в снижении , что уменьшило с 2,5 до 1,923 (-23%).
Наконец, вот схема расчета потенциала, связанного с «квадрупольной» оболочкой: аппроксимируйте оболочку как поверхность переменной плотности массы на среднем радиусе . Эту поверхность можно разложить на кольца постоянных с массой
Потенциал, создаваемый этим кольцом, рассчитан в Jackson, Classical EM , Section 3.3:
Подставляя массовую плотность (которая пропорциональна ) и интегрирование по обнуляет все члены, кроме член, в силу ортогональности полиномов Лежандра, дающий указанный результат. (Разработка для похож.)
Здесь я хотел бы численно проверить теоретическое предсказание фактора в отличие от исходного монопольного аргумента Марка Эйхенлауба. На практике это означает вычисление разницы гравитационных потенциалов между Северным полюсом и экватором и деление на соответствующую разницу потенциалов монополей. По численным причинам на практике лучше вычислить обратную (= обратную) дробь, которую затем следует сравнить с . Поскольку мои навыки программирования ограничены, я только что написал медленный код MAPLE, чтобы выполнить эту работу.
b:=100; f:=.10; a:=b*(1+f); V1:=evalf(4*Pi*a^2*b/3);
xa:=a; ya:= 0; za:=0; xb:=0; yb:=0; zb:=b;
U1a := evalf(V1/sqrt(xa^2 + ya^2 + za^2));
U1b := evalf(V1/sqrt(xb^2 + yb^2 + zb^2));
Ua:=0;Ub:=0;V:=0;
for x from -a-.5 by 1 to a+.5 do
for y from -a-.5 by 1 to a+.5 do
for z from -b-.5 by 1 to b+.5 do
if (x/a)^2 + (y/a)^2 + (z/b)^2 < 1 then
Ua:=Ua + 1/sqrt((x-xa)^2 + (y-ya)^2 + (z-za)^2);
Ub:=Ub + 1/sqrt((x-xb)^2 + (y-yb)^2 + (z-zb)^2);
V:=V+1;
end if;od;od;od;
b;f;Ua;U1a;Ub;U1b;V;V1;Ub-Ua;U1b-U1a;(U1b-U1a)/(Ub-Ua);
Результат с полярным радиусом был
Тот факт, что оценка не улучшается от плоскостности к является артефактом решетки, потому что шаг решетки того же порядка, что и разница в радиусе между северным полюсом и экватором.
Мартин Гейлз
ОХОТНИК НА ТРОЛЛЕЙ
Марк Эйхенлауб
Марк Эйхенлауб
Марк С
Qмеханик
Рон Маймон