CC BY
14
2017 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 4 (62). Вып. 3
МАТЕМАТИКА. МЕХАНИКА. АСТРОНОМИЯ
АСТРОНОМИЯ
УДК 521.14 МЯО 31В05
КОЭФФИЦИЕНТЫ СТОКСА СЖАТОГО ЭЛЛИПСОИДА ВРАЩЕНИЯ,
ЭКВИДЕНСИТЫ КОТОРОГО ПОДОБНЫ ЕГО ПОВЕРХНОСТИ*
К. В. Холшевников1'2, Д. В. Миланов1, В. Ш. Шайдулин1
1 Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
2 Институт прикладной астрономии РАН,
Российская Федерация, 191187, Санкт-Петербург, наб. Кутузова, 10
Теория фигур равновесия активно развивалась в XIX столетии, когда выяснились причины, по которым наблюдаемые массивные небесные тела (Солнце, планеты, спутники) обладают близкой к эллипсоидальной формой. Было установлено, что существуют и в точности эллипсоидальные фигуры. Гравитационный потенциал таких фигур представляется рядом Лапласа, коэффициенты которого (постоянные Стокса 1„) определяются некоторым интегральным оператором. В случае однородного эллипсоида вращения был найден общий член ряда, а для некоторых других распределений масс найдены первые члены ряда. Здесь мы получили общий член ряда для произвольного распределения масс при условии, что эквиденситы (поверхности равной плотности) подобны внешней поверхности эллипсоида вращения. Получены также простые оценки и асимптотика 1п. Библиогр. 13 назв.
Ключевые слова: гравитационный потенциал, ряд Лапласа, эллипсоид.
Введение. Представление гравитационного потенциала эллипсоидов со времен Ньютона являлось важной областью математического естествознания. Сотни работ посвящены этой теме. Результаты собраны в многочисленных монографиях (см., например, [1-8]). В частности, еще Лежандр [9] сумел вычислить коэффициенты Стокса ряда Лапласа однородного эллипсоида вращения, и этот результат приводится в учебниках и монографиях. В работах [7, 10] он использован для доказательства точности оценок скорости убывания общего члена ряда Лаплапса в случае тела аналитической структуры.
* Работа выполнена при финансовой поддержке Программы проведения фундаментальных исследований СПбГУ по приоритетным направлениям (грант 6.37.341.2015). (¡5 Санкт-Петербургский государственный университет, 2017
В настоящей статье мы обобщаем результат Лежандра на случай эллипсоида вращения, эквиденситы (поверхности равной плотности) которого подобны внешней поверхности эллипсоида, а в остальном распределение масс произвольно.
Коэффициенты Стокса. Пусть тело Т ограничено поверхностью сжатого эллипсоида вращения 5 с экваториальной полуосью а и полярной полуосью с < а. Обозначим через е, в первый и второй эксцентриситеты меридионального сечения эллипсоида:
л/а2 — с2 л/а2 — с2 в
Р =-, е = —, р :
угт^2' VT^'
- = (1 -е2)(1 +/32) = 1. (1)
а V71 + /5
Плотность д тела Т предполагается интегрируемой и постоянной на эллипсоидах, подобных поверхности 5.
По симметрии ряд Лапласа гравитационного потенциала V тела Т содержит лишь зональные гармоники четного порядка [7]
п=0
Здесь 0 — гравитационная постоянная, Рп — многочлен Лежандра со стандартной
нормировкой Рп (1) = 1, М — масса Т,
1п = / вгПРп(^)<1т> м=^ (з)
r
T T
где q — плотность, dr — элемент объема. Поскольку /о = 1, а при нечетном n коэффициенты Стокса /п равны нулю, считаем ниже индекс n четным положительным числом.
Сделаем подстановку
х = аи sin в cos А, г = аи \]\ — е2 cos2 в,
z cos Ол/\ - £2
у = а,и sin 0 sin А, — = —== =, (4)
r л/1 - е2 cos2 в
z = cu cos в.
Поверхности u = const являются эллипсоидами S(u)
Х2+У2 , Z2 2 ,,
-2— + — = U , (5)
а2 с2
подобными внешней поверхности S = S(1). Поэтому q зависит лишь от u.
Ниже понадобятся элемент объема сплошного эллипсоида dr и элемент поверхности граничного эллипсоида dS [11]:
е
е
а
c
dr = а2 cu2 sin в du dO d\, dS = а2 sin 6л/ 1-е2 sin2 в dO dX.
(6) 517
Подставляя (4) в (3) и вычисляя интеграл по Л (а для массы и по в), получим
1
0 0 ^
1
M = 4па2 с У д(п)п2 du. (8)
о
Впрочем, для определения массы выгоднее пользоваться ее выражением через среднюю плотность Q (отношение массы к объему):
4па2с , ч
М = —(9)
В интеграле по в заменим переменную:
cos в\/1 — е2 п Р о о Л 1
Р = =, cos 0 = ; , 1-е2 cos2 в =
= Гл О о л ' /7^ о \ / -I TvTlToV ?
Vi -e2cos2e ' -е2)(1 + /32е2) ' 1+/32£2'
sin 9 d,9 =----= .
V(i-e2)(i + /32e2)3
В понятных обозначениях имеем
п 1
У 1 f Pn(№
J Vl^ J (l + /?2£2)(n+3)/2 •
0 -1
Вычисленный Лежандром [9] интеграл (10) приведен в [4, 5]:
(10)
1
Pn(0d£ 2( —l)"/2 ¡3n 2(-1)™/2 VT^e2
у (i+e2e2)(n+3)/2 (n+i)(i+e2)(n+1)/2 n+i 1
(11)
Замечание. Формула (11) доказана для вещественных чисел в2 из интервала — 1 < в2 < 1. По принципу аналитического продолжения она верна при вещественных в2 на луче —1 < в2 < го. В частности, она справедлива при 0 ^ е ^ 1.
Окончательно, будем иметь
1 \/ 1 4 -1 /„ = (_ 1)п/27п) ( [ „/„Л„."+2 1 ( 2
У g(u)un+2 duj í У e(u)u2 duj (12)
00
n + 1 \
0
1
Jn= f _ / g(u)vn+2 du. (13)
(n + 1)g J
0
In
c"
Оценки и асимптотика коэффициентов Стокса. 1. Постоянная плотность. При q = const интеграл (13) элементарен, и мы получаем результат Лежандра:
3е"
Jn = (п + 1)(п + 3) ' (14)
2. Плотность, ограниченная и отделенная от нуля. Пусть qi < q < q2,
где постоянные qs удовлетворяют неравенствам
0 < Q1 < Q < Q2 < то.
По теореме о среднем имеем
3enQ
Jn = (п+1)(п + 3)ё ' (15)
где зависящая от n величина Q заключена между qi и Q2. В частности,
A
J„ = —en, а = 2, 0 < Л* < < А* < то (16)
na
при некоторых постоянных A* , A*.
Замечание. Если плотность постоянна, следует считать 0 < q = qi = Q2 = q = q.
Получим асимптотику An при n ^ то при дополнительном условии непрерывности q слева в точке и = 1. Разобьем промежуток интегрирования в (13) на два точкой 1
ип = 1 -
л/^ТЗ '
По теореме о среднем можем записать 1
(п + 3) У ¿(и)ип+2 Ли = ¿1<+3 + ¿2 (1 - <+3) = [¿1 - ¿2] <+3 + ёа. (17)
о
Здесь зависящие от п величины д8 ограничены, причем ¿2 заключена между нижней и верхней гранью функции д(и) на отрезке ип ^ и ^ 1. Очевидно, ¿>2 ^ ¿(1) при п —> то,
/п+3
lim иП+3 = lim
1
лЛГ+з
0. (18)
В результате, получим
- = А=Щ. (19)
Пусть, например, д = ¿о(1 — 7и), 0 < 7 < 1, т. е. плотность непрерывна и возрастает с глубиной. Тогда будем иметь
1
01 = £>о(1 - 1), 92 = 90, 0 = 3 J д{и)и2 ¿и = - ^7^ £>0 ,
о
= 3£"[(1-7)п+(4-37)Ь
(п + 1)(п + 3)(п + 4)ё ' 1 }
что совпадает с (15) при
в =
(1"7)+ 7
п + 4
Непосредственно проверяются соотношения
№. (21)
01 < в < в < £2 , Иш в = 01 = 0(1), (22)
что подтверждает равенство (19).
3. Плотность, ограниченная сверху, но не снизу. Предположим, что плотность с ростом и от 0 до 1 убывает до нуля. Пусть, например, 0 = 0о(1 — и)^, V > 0. Тогда можем записать
01 = о, 02 = 00 , в = 30оВ(3, V + 1), 3еП 0
•4 = / , - В(п + 3, г/ + 1). (п + 1)в
Выразим бета-функции через гамма-функции:
_ е"Г(г/ + 4)Г(п + 3)
2(п+1)Г(п + г/ + 4) ' 1 ;
Пусть V — натуральное число. Тогда (23) равносильно
,/„ = -_-_^ + , ,_(24)
Отсюда вытекает асимптотика (19) при
* = * + 2, (25)
Пусть V — произвольное положительное число. Асимптотика гамма-функции при п ^ го подтверждает справедливость представления (19), (25) при замене (V + 3)! на Г^ + 4).
Этот результат нетрудно обобщить на случай плотности вида 0 = 0од(и)(1 — и)^, где д ограничена, непрерывна слева в точке и =1, причем $(1) = 1. Поступая так же, как при выводе формулы (17), получим
1 «л 1
У 0(и)ии+2 ¿и = 0О(д1 — в2) У (1 — иГии+2 ¿и + 0Од2 У (1 — и)"ии+2 ¿и. (26) о оо
Здесь вя ограничены, причем $2 ^ д(1) = 1 при п ^ го. Второй интеграл справа асимптотически равен
В(п + 3, V + 1) = —ттт-т;— ~ —гт~ ■
у ' Г(п + V + 4) (п)^+1
Первый интеграл справа согласно (18) убывает быстрее, как е ^. Мы приходим к асимптотике (19) при
<т = z/ + 2,
в
Для целого v ^ 0 и v раз дифференцируемой в точке и = 1 функции д справедливо
A=3gW(l) в
4■ Бесконечная в центре плотность. Плотность небесных тел в центре может быть весьма большой. Например, в центре Солнца плотность на 9 порядков выше, чем на границе фотосферы [12]. Поэтому допустима модель убывающей с ростом и плотности, бесконечной в точке и = 0, но интегрируемой по объему. Иными словами, интегрируемой должна быть функция и2 д(и). По принципу поверхностного слоя [7, 13] это не влияет на результаты пп. 2 и 3. Впрочем, это легко доказать и формально, выделяя в T малое ядро.
5■ Бесконечная на периферии плотность■ Это допущение относится только к модельным телам. Пусть тело T является тонким слоем между двумя элипсои-дами S(1 — ¿и) и S(1), а его плотность постоянна. Точнее, она не зависит от координат и, 0, Л, но зависит от параметра ¿и > 0. Согласно (6) масса элементарного параллелепипеда и его поверхностная плотность определяются равенствами
2 dM с дйи аМ = а сд sin в du ав а\, gsur t = —— = — -. (27)
dS Vi ~ £2 sin2 в С другой стороны, полная масса T согласно (8) может быть записана в виде
4п
М = — а2сд[ 1 - (1 - du)3} = Ana2cgdu + 0(du2).
Переходя к пределу д ^ то, ¿и ^ 0 при постоянной массе M, получим двумерное тело, эллипсоид S, с конечной поверхностной плотностью
^Vl-lmV (28)
Из (27), (28) вытекает сдйи = J, так что (7) переходит в
7Г , _ ч
27Га25 f о о „.„/п . / COS0A/1 — £2 *
/„ = —— / (1 - еА соъА в)п'А ъшв Рп\ , (Л9.
м У 1 ; I VI - е2 СОЗ2 в /
о
Последний интеграл дается формулами (10), (11). Окончательно, получаем
= (29)
п +1
6. Более сложная асимптотика. Асимптотика может и отличаться от вида (19). Пусть, например, д = — ¿о 1п(1 — и). Для вычисления входящих в (12) интегралов достаточно воспользоваться легко проверяемой формулой
— (к + 1) У ик 1п(1 — и) Ли = (1 — ик+1) 1п(1 — и) + + 1),
где
п 1
8 = 1 в
Отсюда получаем
1
-(к + 1) У ик 1п(1 - м) = + 1).
о
В частности, будем иметь
/11
д(и)и2 <1и = —д0,
1
о
1
[ д(и)ип+2 аи = Ь(п + 3).
} п + 3
Окончательно,
18е"Ь(п + 3) 18е" 1пп
п_ 11(п+1)(п + 3)' 11п2 ' ^
Заключение. Мы определили коэффициенты Стокса сплошного эллипсоида, поверхности равной плотности которого подобны граничному эллипсоиду. В случае ограниченной сверху и отделенной от нуля плотности коэффициенты убывают по закону (19), причем а = 2. Если плотность спадает до нуля на поверхности, то а может быть произвольным числом, большем двух, в зависимости от скорости убывания плотности. Если плотность бесконечна в центре масс, это не влияет на оценки (16), (19). Напротив, бесконечная на периферии плотность может повлечь а < 2 или более сложный закон, чем (19). Например, для распределения масс по поверхности граничного эллипсоида имеем а = 1. Это является следствием известного факта, что дифференциальные свойства потенциала простого слоя хуже свойств потенциала трехмерного тела. Соответственно, ряд Лапласа двумерных тел сходится медленнее, чем для трехмерных.
Полученные простые оценки и асимптотика позволяют найти область сходимости ряда Лапласа (2). Важно, что при любой функции £>(м) эта область одна и та же: г > ае, а при ограниченности £?(м) вблизи поверхности 5 даже г ^ ае. Впрочем, уже для градиента потенциала ряд расходится (или сходится условно) в некоторых точках сферы г = ае. Напомним, что ряд Лапласа сходится именно к потенциалу (т.е. в (2) имеет место равенство) в пересечении указанной области и внешнего к эллипсоиду пространства. Таким образом, ряд (2) представляет потенциал во всем внешнем пространстве при е ^ 1 /л/2 = 0.707107, или, что то же, с/а ^ 1 /л/2- В противном случае ряд расходится в части внешнего пространства в окрестности полюсов.
Литература
1. Пуанкаре А. Фигуры равновесия жидкой массы. М.; Ижевск: РХД, 2000. 208 с.
2. Аппель П. Фигуры равновесия вращающейся жидкости. Л.; М.: ОНТИ, 1936. 376с.
3. Сретенский Л. Н. Теория ньютоновского потенциала. М.; Л.: ГИТТЛ, 1946. 320 с.
4. Гобсон Е. В. Теория сферических и эллипсоидальных функций. М.: ИЛ, 1952. 476 с.
5. Дубошин Г.Н. Теория притяжения. М.: Физматлит, 1961. 288 с.
6. Чандрасекар С. Эллипсоидальные фигуры равновесия. М.: Мир, 1973. 289 с.
7. Антонов В. А., Тимошкова Е.И., Холшевников К. В. Введение в теорию ньютоновского потенциала. М.: Наука, 1988. 270 с.
8. Кондратьев Б. П. Теория потенциала и фигуры равновесия. М.; Ижевск: Изд. ИКИ, 2003.
624 с.
9. Legendre A. M. Memoires présentes par 1es savants etrangers. 1785. Vo1.X.
10. Kholchevnikov C. Le deve1oppement du potentie1 dans 1е cas d'une densite ana1ytique // Ce1estia1 Mechanics. 1971. Vo1.3, Ко 2. P. 232-240.
11. Фавар Ж. Курс локальной дифференциальной геометрии. М.: ИЛ, 1960. 560с.
12. Guenther D.B., Demarque P., Kim Y. C., Pinsonneault M. H. Standard so1ar mode1 // The Astrophysica1 Journa1. 1992. Vo1. 387. P. 372-393.
13. Kholshevnikov K. V., Shaidulin V. Sh. Existence of a c1ass of irregu1ar bodies with a higher convergence rate of Lap1ace series for the gravitationa1 potentia1 // Ce1est. Mech. Dyn. Astr. 2015. Vo1.122, issue 4. P. 391-403.
Статья поступила в редакцию 13 декабря 2016 г.; рекомендована в печать 30 марта 2017 г. Сведения об авторах
Холшевников Константин Владиславович —доктор физико-математических наук, профессор; [email protected]
Миланов Данила Владимирович — аспирант; [email protected]
Шайдулин Вахит Шамильевич — кандидат физико-математических наук, доцент; [email protected]
STOKES CONSTANTS OF AN OBLATE ELLIPSOID OF REVOLUTION WITH EQUIDENSITES HOMOTHETIC TO ITS SURFACE
Konstantin V. Kholshevnikov1'2, Danila V. Milanov1, Vakhit Sh. Shaidulin1
1 St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation; [email protected], [email protected], [email protected]
2 Institute of App1ied Astronomy RAS, nab. Kutuzova, 10, St. Petersburg, 191187, Russian Federation; [email protected]
Theory of the figures of equi1ibrium was deve1oped active1y during XIX century when causes making the form of observab1e massive ce1estia1 bodies (the Sun, p1anets, moons) a1most e11ipsoida1 were discovered. The existence of exact1y e11ipsoida1 figures was estab1ished. The gravitationa1 potentia1 of such figures can be presented as Lap1ace series. Its coefficients (Stoke's constants In) are defined by a certain integra1 operator. The genera1 term of the series was found in case of a homogeneous e11ipsoid, and first few terms were found for severa1 other mass distributions. Here we have found the genera1 term of the series for an arbitrary mass distributions under condition that equidensits (surfaces of equa1 density) are homothetic to the outer surface of the e11ipsoid of revo1ution. Simp1e estimates and asymptotics of In are a1so found. Refs 13.
Keywords: gravitationa1 potentia1, Lap1ace series, e11ipsoid. References
1. Poincare H., Figures d'équilibre d'une masse fluide (Gauthier-Vi11ars, Paris, 1900) [in French].
2. Appe11 P., Figures d'équilibre d'une masse liquide homogène en rotation (Gauthier-Vi11ars, Paris, 1932) [in French].
3. Sretensky L. N., Theory of Newtonian Potential (GITTL, Moscow, Leningrad, 1946, 320 p.) [in Russian].
4. Hobson E. W., The Theory of Spherical and Ellipsoidal Harmonics (Cambridge Univ. Press, Cambridge, 1931, 476p.).
5. Duboshin G.N., Theory of Attraction (Fizmatgis, Moscow, 1961) [in Russian].
6. Chandrasekhar S., Ellipsoidal Figures of Equilibrium (Dover Pub1s, 1987, 264 p.).
7. Antonov V.A., Timoshkova E. I., Kho1shevnikov K. V., Introduction to the Theory of Newtonian Potential (Nauka, Moscow, 1988, 270 p.) [in Russian].
8. Kondratiev В. P., Theory of Potential, and Figures of Equilibrium (Inst. Kosm. Res. Pub1., Moscow, Ijevsk, 2003, 624p.) [in Russian].
9. Legendre A.M., Mémoires présentées par les savants étrangers (X, 1785) [in French].
10. Kholchevnikov C., "Le developpement du potentiel dans le cas d'une densite analytique", Celestial Mechanics 3(2), 232-240 (1971).
11. Favard J., Cours de géométrie différentielle locale (Gauthier-Villars, Paris, 1957) [in French].
12. Guenther D.B., Demarque P., Kim Y. C., Pinsonneault M. H., "Standard solar model", The Astrophysical Journal 387, 372-393 (1992).
13. Kholshevnikov K. V., Shaidulin V. Sh., "Existence of a class of irregular bodies with a higher convergence rate of Laplace series for the gravitational potential", Celest. Mech. Dyn. Astr. 122, 391-403 (2015).
Для цитирования: Холшевников К. В., Миланов Д. В., Шайдулин В. Ш. Коэффициенты Стокса сжатого эллипсоида вращения, эквиденситы которого подобны его поверхности // Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 3. С. 516-524. DOI: 10.21638/11701/spbu01.2017.313
For citation: Kholshevnikov K. V., Milanov D. V., Shaidulin V. Sh. Stokes constants of an oblate ellipsoid of revolution with equidensites homothetic to its surface. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2017, vol. 4(62), issue 3, pp. 516-524. DOI: 10.21638/11701/spbu01.2017.313
