Научная статья на тему 'Распределение плотности в газовой оболочке планеты'

Распределение плотности в газовой оболочке планеты Текст научной статьи по специальности «Техника и технологии»

CC BY
65
4
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
Газовая оболочка планеты / самогравитация / барометрическая формула / уравнение Пуассона / плотность газа в газовых скоплениях / радиус Хилла / The planet gaseous environment / self-gravity / barometric formula / Poisson equation / gas density in gaseous accumulations / Hill radius

Аннотация научной статьи по технике и технологии, автор научной работы — Е. Н. Бодунов, А. О. Семенов

Цель: Теоретическое исследование распределения плотности газа вокруг планеты с учетом самогравитации. Методы: Предложенные в настоящей работе уравнения для плотности газовой оболочки планеты с соответствующими граничными условиями решаются аналитически и численно методом Рунге — Кутта. Результаты: Впервые с помощью метода подобия проведен анализ численного решения уравнений для всего пространства, в котором преобладает гравитационное влияние планеты — в диапазоне расстояний от поверхности планеты до радиуса Хилла. Вблизи планеты решение совпадает с классической барометрической формулой, на промежуточных расстояниях — с барометрической формулой, учитывающей зависимость ускорения свободного падения от расстояния до планеты, на больших расстояниях — с зависимостью плотности для сингулярной изотермической газовой сферы, обусловленной самогравитацией. Практическая значимость: На основе полученного решения создана единая картина распределения плотности газовой оболочки планеты. Представленные в работе результаты могут быть полезны как для преподавателей физики вузов, так и для научных сотрудников, занимающихся астрофизикой.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по технике и технологии , автор научной работы — Е. Н. Бодунов, А. О. Семенов

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Density Distribution in the Gaseous Envelope of the Planet

Purpose: Theoretical study of gas density distribution around the planet given self-gravity. Methods: Proposed in this paper equations for the planet gaseous envelope density with corresponding boundary conditions are solved analytically and numerically by Runge-Kutta method. Results: For the first time, the analysis of numerical solution of equations for all the space where the planet gravitational influence prevails in the distance range from the planet surface till Hill radius has been pursued using similarity method. Near the planet, the solution coincides with classical barometric formula, at intermediary distances, with barometric formula which takes into consideration the dependence of free fall acceleration from the distance till the planet, at large distances, with the dependence of density for singular isothermal gas sphere due to self-gravity. Practical significance: On the basis of the solution obtained, the unified picture of the planet gaseous envelope density distribution was analyzed. The results presented in the paper can be useful both for university physics professors and for the researchers involved in astrophysics.

Текст научной работы на тему «Распределение плотности в газовой оболочке планеты»

УДК 533 + 524.7

Распределение плотности в газовой оболочке планеты Е. Н. Бодунов, А. О. Семенов

Петербургский государственный университет путей сообщения Императора Александра I, Российская Федерация, 190031, Санкт-Петербург, Московский пр., 9

Для цитирования: Бодунов Е. Н., Семенов А. О. Распределение плотности в газовой оболочке планеты // Известия Петербургского университета путей сообщения. — СПб.: ПГУПС, 2023. — Т. 20. — Вып. 1. — С. 101-109. DOI: 10.20295/1815-588Х-2023-1-101-109

Аннотация

Цель: Теоретическое исследование распределения плотности газа вокруг планеты с учетом самогравитации. Методы: Предложенные в настоящей работе уравнения для плотности газовой оболочки планеты с соответствующими граничными условиями решаются аналитически и численно методом Рунге — Кутта. Результаты: Впервые с помощью метода подобия проведен анализ численного решения уравнений для всего пространства, в котором преобладает гравитационное влияние планеты — в диапазоне расстояний от поверхности планеты до радиуса Хилла. Вблизи планеты решение совпадает с классической барометрической формулой, на промежуточных расстояниях — с барометрической формулой, учитывающей зависимость ускорения свободного падения от расстояния до планеты, на больших расстояниях — с зависимостью плотности для сингулярной изотермической газовой сферы, обусловленной самогравитацией. Практическая значимость: На основе полученного решения создана единая картина распределения плотности газовой оболочки планеты. Представленные в работе результаты могут быть полезны как для преподавателей физики вузов, так и для научных сотрудников, занимающихся астрофизикой.

Ключевые слова: Газовая оболочка планеты, самогравитация, барометрическая формула, уравнение Пуассона, плотность газа в газовых скоплениях, радиус Хилла.

Введение

В учебниках по общей физике зависимость давления газар(г) от высоты г = г - Яс над поверхностью планеты (Яс — радиус планеты, г — расстояние от центра планеты) обычно описывается барометрической формулой, имеющей вид:

рМ=е-z / н

(1)

Ро

где р0 — давление на поверхности планеты (при г = 0);

Н0 — шкала высот (высота, на которой давление газа в атмосфере уменьшается в е раз) или же высота однородной атмосферы (толщина, которую имела бы атмосфера, если всюду в ней давление было бы такое же, как на поверхности):

н 0 = kT

kT0 R

mg0 GMcm

(2)

В формулах (1) и (2) k — постоянная Боль-цмана, g0 — ускорение свободного падения, g0 = ОЫс/Яс2 (О — гравитационная постоянная, Мс — масса планеты), Т0 — температура. Формула (1) получена при следующих предположениях: атмосферный газ считается идеальным, состоит из молекул массы т, находится при постоянной температуре Т = Т0 в гравитационном поле планеты, ускорение свободного падения g не зависит от г и равно ускорению свободного падения g0 на поверхности планеты (при г = 0). При выводе уравнения (1) обычно постулируется, что давление на высоте г обусловлено

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

Несмотря на свою простоту и сделанные при ее выводе предположения, формула (1) достаточно хорошо передает зависимость давления от высоты в нижних слоях земной атмосферы. Отметим, что так как идеальный газ подчинятся уравнению р = ркТ / т (р — плотность газа, р = пт, п — концентрация газа), то барометрическая формула (1) описывает также распределение концентрации и плотности газа в изотермической атмосфере, так как п^)/п0 = р^)/р0 = р^)/р0 (п0 и р0 — концентрация и плотность газа на поверхности планеты).

Если учесть зависимость ускорения свободного падения от высоты z ^ = ОЫс/(Яс + z)2 = ОЫс/г2), то барометрическая формула приобретает вид [1, 2]:

Р (r) Р п (г) —— = — = = exp

Р Р п

f GM m r

J0 p0 f

l0

kZ

1 1

r

\\

J )

= exp

R f1 - R

r

v Ho V

(3)

Уравнение (3) является следствием распределения Больцмана, согласно которому в состоянии термодинамического равновесия отношение концентраций газа на разных высотах (Яс и г) в потенциальном поле определяется разностью потенциальных энергий молекул газа ЩДс) -Ш(г), деленной на

п

(r )

п

= exp

exp

W (Rc ) W (r)

kTn

kTn

0

GM m 1 GM m 1

- + -

(4)

kT0 R„

kT

Несмотря на давнюю историю, барометрическая формула продолжает интересовать препо-

давателей физики и научных работников с точки зрения ее использования в учебной работе, обобщения на случай неидеального газа, наличия центробежных сил и так далее [1-7], применения в различных приложениях, начиная с метеорологических [8-11] и заканчивая физикой плазмы [12], ядерной физикой [13, 14] и астрофизикой [15-18].

Недостатком обобщенной барометрической формулы (3), согласно, например, классическому учебнику Л. Д. Ландау и Е. М. Лифшица [2], является стремление концентрации газа п при г ^ ю к постоянной величине, п(г ^ ю) = п0ехр(-^с/#0), что некорректно, так как конечное количество газа не может быть распределено по бесконечному объему с отличной от нуля концентрацией. Этот недостаток барометрической формулы (3) можно исправить, если предположить уменьшение температуры газа до нуля с увеличением высоты быстрее, чем (1 + z/RcУ2 [1, 7]. Однако известно, что температура газа на больших высотах в атмосфере отлична от нуля.

Дополнительным фактором, не учтенным в (3) и приводящим к удержанию газа вблизи планеты, является самогравитация [15-18]: газ на расстоянии г удерживается не только притяжением со стороны планеты, но и со стороны слоев газа, расположенных ниже — на расстояниях, меньших г. Как известно [15-18], самогравитация используется в астрофизике для объяснения пространственного строения звезд, галактик и газо-вопылевых структур.

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

Система уравнений, численное решение и обсуждение результатов

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

сферически симметричной самогравитирующей оболочкой из идеального изотермического газа, находящегося в гидростатическом равновесии. Условие гидростатического равновесия такой оболочки можно записать в виде системы из двух уравнений [19, 20]:

dp dr dM dr

GM --— P>

r

= 4 nr 2p,

(5)

где М — масса, заключенная в сфере радиуса г.

Дифференцируя первое уравнение (5) по г и подставляя в него второе уравнение, с учетом уравнения состояния идеального газа ^ = р^ / m) получаем уравнение Пуассона:

d_ r2 dr

dlnp dr

m

= -4nG — p. kT

(6)

Это уравнение можно представить также в интегральном виде:

Р(r)_

GmM

kT

1 R

1

-ir J R

Gm

kT

X"T(iRR 4np(y)y2dy)dx

= exp

(7)

В приближении малой массы атмосферы по сравнению с массой планеты (т. е. когда вклад интегрального слагаемого в экспоненте мал) уравнение (7) можно решать итерационным методом. Из уравнения (7), эквивалентного уравнению (6), отчетливо видно, что на малых расстояниях от планеты и при малых плотностях газа (р ^ 0), когда самогравитацией можно пренебречь (второе слагаемое под знаком экспоненты равно нулю), распределение плотности газа имеет вид барометрической формулы (3). На больших расстояниях, где самогравитация преобладает над притяжением планеты (первым слага-

емым под знаком экспоненты в (7) можно пренебречь), уравнения (6) и (7) имеют аналитическое решение:

p(r):

kT0 1 _ 2 Pc HR

lnGm r2 3 p0

(8)

Введем безразмерную координату х и два безразмерных параметра:

x = r / Rc, a = Rc / Ho, b = Po/Pc

(9)

где рс — плотность ядра (планеты), Мс = (4п/3) рсЯс3. Параметр а — безразмерный радиус планеты, выраженный в единицах Н0. Параметр Ь есть отношение плотности газа р0 у поверхности ядра (планеты) к плотности ядра рс (планеты). Используя параметры (9), уравнение (8) преобразуем к виду:

p(r) _ 2 1 1

Po

3 ab x2

(10)

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

Проинтегрировав второе уравнение системы (5) от 0 до г > Я , получим:

M = Mc + ir 4nr 2pdr.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(11)

Для уменьшения количества параметров задачи (с целью изучения свойств решения системы уравнений (5)) используем метод подобия, успешно применяющийся в теории планетных атмосфер (см., например, [21, 22]).

Введем безразмерные функции:

u = р/Ро, f = M /Mc. (12)

С помощью безразмерных величин (12) и параметров (9) система уравнений (5) записывается в виде:

p

0

Параметры Земли, Юпитера и протоюпитера (в скобках указаны единицы измерения)

Планета Земля Юпитер Протоюпитер1*

Плотность планеты рс (кг/м3) 5500 1326 5500

Плотность газовой оболочки у поверхности планеты р0 (кг/м3) 1,22 0,16 ~2002*

Средняя молекулярная масса т (а. е. м) 28,5 2,3 2,3

Температура у поверхности Т0 (К) 288 129 5000

Масса планеты Мг (МЕ) 3* 1 318 1-254*

Радиус планеты Яг = (3М74гсрс)1/3 (ЯЕ) 5* 1 11 1-3

Ускорение свободного падения g0 = (м/с2) 9,8 25,9 9,8-28,7

Шкала высот Н0 = kT0/mg0 (км) 8,6 18 630-1840

« = КН 744 3886 3,46-29,6

Ь = Р0/Рс 0,0002 3 ■ 10-5 0,036

Примечания:

'"Под протоюпитером понимается твердое ядро будущей планеты на месте современного Юпитера в протопла-нетном облаке на раннем этапе образования планет Солнечной системы; 2* порядок величины взят из [20]; 3* М = 5,97217 • 10-24 кг — масса Земли;

Е ' '

4* массы твердых ядер внутри современных планет—гигантов Солнечной системы оцениваются как 15-25 МЕ [23]; 5* ЯЕ = 6371 км — радиус Земли.

af

x

du dx

f = 3bx2u. dx

(13)

Система уравнений (13) решалась численно методом Рунге — Кутта со следующими граничными условиями на поверхности планеты:

u(1) = 1 и Д1) = 1. (14)

Расчеты производились для набора параметров a и b (9), характерных для Земли, Юпитера и протоюпитера (таблица).

Результаты численных расчетов приведены на рисунке (синие, зеленые и красные кривые). Пунктирной серой кривой добавлено аналитическое решение уравнений (13) без учета самогравитации (при f = 1 ):

Ро

: u = exp

a

( x -1)

(15)

Это уравнение есть уравнение (3), записанное через безразмерные параметры. На больших высотах (х >> 1) функция (15) стремится к постоянному (асимптотическому) значению:

Р / Ро = exP (-a).

(16)

С учетом самогравитации численные решения отклоняются от асимптотики (16) (степень отклонения зависит от параметра Ь, см. рисунок) и приближаются к другому предельному случаю — к выражению (10) с зависимостью р/р0 ~ 1/х2, соответствующей газовой оболочке с массой, значительно большей массы планеты. Приравнивая асимптотики (16) и (10), можно оценить безразмерную высоту х , начиная с которой становится существенно влияние самогравитации на плотность газовой оболочки:

' 3ab

(17)

10° 101 102 103 104 ю5

г/В с

Зависимость плотности газовой оболочки планеты от безразмерного расстояния х = г / Яс при различных значениях параметров а и Ь. Красные кривые — а = 30, зеленые — а = 10, синие — а = 3,5, сплошные кривые — Ь = 1, штриховые — Ь = 0,036, штрихпунктирные — Ь = 0,002. Серая пунктирная кривая — расчет по барометрической формуле (15), серая штриховая кривая — расчет по формуле (10), вертикальная черная штриховая линия — безразмерный радиус Хилла для системы Юпитер и Солнце (18)

Отметим, что область гравитационного влияния планеты ограничена радиусом Хилла [24] гн — радиусом сферической области вокруг планеты, внутри которой гравитационное влияние планеты преобладает над гравитационным влиянием звезды:

D 3

M '3M.

(18)

где Б — расстояние от планеты до звезды; М* — масса звезды.

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

Предполагая, что масса М в основном определяется массой планеты М формулу (18) можно переписать в виде:

H _

R

(19)

где Я* — радиус звезды;

рс и р* — плотности планеты и звезды соответственно.

Считая, что плотности планеты рс и звезды р* примерно равны, и учитывая, что рс и р* входят в уравнение (19) под знаком кубического корня, можно считать, что отношение гн/Яс

определяется только отношением D/R*, т. е. значение rH/Rc является фиксированным для заданного положения планеты и не зависящим от параметров a и b.

На рисунке вертикальной черной штриховой линией показан безразмерный радиус Хилла (19) для системы Юпитер и Солнце (rH/Rc й 1,06 • 103). Из рисунка следует, что для параметров a и b, характерных для протоюпитера, внутри радиуса Хилла наблюдаются все три предельных случая зависимости плотности газовой оболочки, соответствующих барометрической формуле (3), выходу на плато (асиптотика (16)) и обратной степенной зависимости (10).

Для Земли (rH/Rc ~ 0,21-103) и Юпитера, для которых параметр a очень большой (744 и 3886 соответственно), а плотность р0 мала, плотность газовой оболочки р в области плато (в соответствии с формулой (15)) чрезвычайно мала, и расстояние, начиная с которого следует учитывать самогравитацию, становится значительно больше радиуса Хилла. Поэтому для этих планет на всех разумных расстояниях (меньших радиуса Хилла) справедливы формулы (3) и (15).

Без сомнения, ситуация, подобная протоюпи-теру, возможна и для других планетных систем.

Заключение

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

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

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

Библиографический список

1. Berberan-Santos M. N. On the Barometric Formula / M. N. Berberan-Santos, E. N. Bodunov, L. Pogliani // Amer. J. Phys. — 1997. — Vol. 65. — P. 404. — DOI: 10.1119/ 1.18555.

2. Landau L. D. Statistical Physics / L. D. Landau,

E. M. Lifshiz. — Pergamon Press, 1969. — Vol. 5. — P. 108.

3. Berberan-Santos M. N. Liquid — Vapor Equilibrium in a Gravitational Field / M. N. Berberan-Santos, E. N. Bodunov, L. Pogliani // Amer. J. Phys. — 2002. — Vol. 70. — P. 438. DOI: 10.1119/1.1424264.

4. Berberan-Santos M. N. On the Barometric Formula inside the Earth / M. N. Berberan-Santos, E. N. Bodunov, L. Pogliani // J. Math. Chem. — 2010. — Vol. 47. — P. 990. — DOI: 10.1007/s10910-009-9620-7.

5. Rodrigues D. S. Analyzing Atmospheric Pressure Variations in Time and Height: a Didactic Proposal Employing a Smartphone Barometer / D. S. Rodrigues,

F. J. Arnold // Revista Brasileira de Ensino de Física. — 2022. — Vol. 44. — P. e20210422. — DOI: 10.1590/1806-9126-RBEF-2021-0422.

6. Slisko J. The Physical Cause of Atmospheric Pressure: Weight of Air or Molecular Motion and Impacts? / J. Slisko, T. M. Topalovic, M. Bozic // The Physics Teacher. — 2021. — Vol. 59. — P. 470. — DOI: 10.1119/10.0006132.

7. Bodunov E. N. Barometric Formula for Non-Isothermal Atmosphere / E. N. Bodunov, G. G. Khokhlov // J. Phys.: Conf. Ser. — 2021. — Vol. 2131. — P. 022053. — DOI: 10.1088/1742-6596/2131/2/022053.

8. Wang J. Improving the Vertical Modeling of Tropospheric Delay / J. Wang, K. Balidakis, F. Zus et al. // Geophys. Research Lett. — 2022. — Vol. 49(5). — P. e2021GL096732. — DOI: 10.1029/2021GL096732.

9. Michalak P. Impact of Air Density Variation on a Simulated Earth-to-Air Heat Exchanger's Performance / P. Michalak // Energies. — 2022. — Vol. 15(9). — P. 3215. — DOI: 10.3390/en15093215.

10. Stenner C. Development and Persistence of Hazardous Atmospheres in a Glaciovolcanic Cave System — Mount Rainier, Washington, USA. / C. Stenner, A. Pflitsch, L. J. Florea et al. // J. Cave & Karst Studies. — 2022. — Vol. 84(2). — P. 66. — DOI: 10.4311/2021ex0102.

11. Chkeir S. Nowcasting extreme rain and extreme wind speed with machine learning techniques applied to different input datasets / S. Chkeir, A. Anesiadou, A. Mascitelli et al. // Atmospheric Research. — 2023. — Vol. 282. — P. 106548. — DOI: 10.1016/j.atmosres.2022.106548.

12. Dubinov A. E. Mathematical Tricks for Pseudopotentials in the Theories of Nonlinear Waves in Plasmas / A. E. Dubinov // Physics of Plasmas. — 2022. — Vol. 29. — P. 020901. — DOI: 10.1063/5.0078573.

13. Dubinov A. E. Barometric Formula for Ultrarelativistic Degenerate Fermi-Gases / A. E. Dubinov // Astrophysics. — 2020. — Vol. 63(4). — P. 580. — DOI: 10.1007/s10511-020-09660-1.

14. Moon J. Design of Air-cooled Waste Heat Removal System with String Type Direct Contact Heat Exchanger and Investigation of Oil Film Instability / J. Moon, Y. H. Jeong, Y. Addad // Nuclear Engineering and Technology. — 2020. — Vol. 52(4). — P. 734. — DOI: 10.1016/j.net.2019.10.010.

15. Shu F. N. The Physics of Astrophysics / F. N. Chu. — California: Univer. Sci. Books, 1992. — Vol. II — P. 246.

16. Ткаченко Р. В. Гравитационная неустойчивость газопылевых околоядерных дисков близких галактик / Р. В. Ткаченко, В. И. Корчагин, Б. Б. Жмайлов. // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. — 2022. — Т. 9(67). — № 3. — С. 561. — DOI: 10.21638/spbu01.2022.316.

17. Patra N. N. Theoretical modelling of two-component molecular discs in spiral galaxies / N. N. Patra // Astronomy & Astrophysics. — 2020. — Vol. 638. — P. A66. — DOI: 10.1051/0004-6361/201936483.

18. Bodmer D. Asymptotic tracking position control with active oscillation damping of a multibody Mars vehicle using two artificial augmentation approaches / D. Bodmer, M. Krenmayr, F. Holzapfel // CEAS Space J. — 2022. — Vol. 14. — P. 125. — DOI: 10.1007/s12567-021-00364-6.

19. Тейлер Р. Строение и эволюция звезд / Р. Тейлер. — М.: Мир, 1973. — 276 с.

20. Mizuno H. Formation of the Giant Planets / H. Mizuno // Progress of Theoretical Physics. — 1980. — Vol. 64(2). — P. 544. — DOI: 10.1143/PTP.64.544.

21. Semenov A. O. Upper Thermal Boundary Layer of Planetary Atmosphere: An Attempt of Developing a General Model / A. O. Semenov, G. M. Shved // Icarus. — 2008. — Vol. 194(1). — P. 290. — DOI: 10.1016/j.ica-rus.2007.08.040.

22. Shved G. M. The Standard Problem of Nonlocal Thermodynamic Equilibrium Radiative Transfer in the Rovibrational Band of the Planetary Atmosphere / G. M. Shved, A. O. Semenov // Solar System Research. — 2001.—Vol. 35.—P. 212. — DOI: 10.1023/A:1010478906172.

23. Bodenheimer P. Calculations of the Accretion and Evolution of Giant Planets: The Effects of Solid Cores / P. Bodenheimer, J. B. Pollack // Icarus. — 1986.—Vol. 67(3). — P. 391. — DOI: 10.1016/0019-1035(86)90122-3.

24. Pecnik B. Giant Planet Formation. A First Classification of Isothermal Protoplanetary Equilibria / B. Pecnik, G. Wuchterl // Astronomy & Astrophysics. — 2005. — Vol. 440. — P. 1183. — DOI: 10.1051/0004-6361:200500005.

Дата поступления: 23.01.2023 Решение о публикации: 14.02.2023

Контактная информация:

БОДУНОВ Евгений Николаевич — д-р физ.-мат. наук, проф., зав. кафедрой «Физика»; evgeny.bodunov@inbox.ru

СЕМЕНОВ Алексей Олегович — канд. физ.-мат. наук, доц.; aleksem@mail.ru

Density Distribution in the Gaseous Envelope of the Planet E. N. Bodunov, A. O. Semenov

Emperor Alexander I St. Petersburg State Transport University, 9, Moskovsky pr., Saint Petersburg, 190031, Russian Federation

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

For citation: Bodunov E. N., Semenov A. O. Density Distribution in the Gaseous Envelope of the Planet //

Proceedings of Petersburg Transport University, 2023, vol. 20, iss. 1, pp. 101-109. (In Russian). DOI: 10.20295/1815-588X-2023-1-101-109

Summary

Purpose: Theoretical study of gas density distribution around the planet given self-gravity. Methods: Proposed in this paper equations for the planet gaseous envelope density with corresponding boundary conditions are solved analytically and numerically by Runge-Kutta method. Results: For the first time, the analysis of numerical solution of equations for all the space where the planet gravitational influence prevails - in the distance range from the planet surface till Hill radius - has been pursued using similarity method. Near the planet, the solution coincides with classical barometric formula, at intermediary distances, - with barometric formula which takes into consideration the dependence of free fall acceleration from the distance till the planet, at large distances, - with the dependence of density for singular isothermal gas sphere due to self-gravity. Practical significance: On the basis of the solution obtained, the unified picture of the planet gaseous envelope density distribution was analyzed. The results presented in the paper can be useful both for university physics professors and for the researchers involved in astrophysics.

Keywords: The planet gaseous environment, self-gravity, barometric formula, Poisson equation, gas density in gaseous accumulations, Hill radius.

References

1. Berberan-Santos M. N., Bodunov E. N., Pogliani L. On the Barometric Formula. Amer. J. Phys, i997, vol. б5, p. 404. DOI: in.iii9/i.i8555.

2. Landau L. D., Lifshiz E. M. Statistical Physics. Pergamon Press, ^9, vol. 5, p. 108.

3. Berberan-Santos M. N., Bodunov E. N., Pogliani L. Liquid — Vapor Equilibrium in a Gravitational Field, 2nn2, vol. 70, p. 438. DOI: Ш.Ш9/и4242б4.

4. Berberan-Santos M. N., Bodunov E. N., Pogliani L. On the Barometric Formula inside the Earth. J. Math. Chem, 2010, vol. 47, p. 990. DOI: Ш.ЮОШДОЮ^^П^.

5. Rodrigues D. S., Arnold F. J. Analyzing Atmospheric Pressure Variations in Time and Height: a Didactic Proposal Employing a Smartphone Barometer. Revista Brasileira de Ensino de Física, 2022, vol. 44, p. e2n2in422. DOI: 10.1590/1806-9126-RBEF-2021-0422.

6. Slisko J., Topalovic T. M., Bozic M. The Physical Cause of Atmospheric Pressure: Weight of Air or Molecular

Motion and Impacts? The Physics Teacher, 2021, vol. 59, p. 470. DOI: 10.1119/10.0006132.

7. Bodunov E. N., Khokhlov G. G. Barometric Formula for Non-Isothermal Atmosphere. J. Phys.: Conf. Ser., 2021, vol. 2131, p. 022053. DOI: 10.1088/17426596/2131/2/022053.

8. Wang J., Balidakis K., Zus F. et al. Improving the Vertical Modeling of Tropospheric Delay. Geophys. Research Lett, 2022, vol. 49(5), p. e2021GL096732. DOI: 10.1029/2021GL096732.

9. Michalak P. Impact of Air Density Variation on a Simulated Earth-to-Air Heat Exchanger's Performance. Energies, 2022, vol. 15(9), p. 3215. DOI: 10.3390/ en15093215.

10. Stenner C., Pflitsch A., Florea L. J. et al. Development and Persistence of Hazardous Atmospheres in a Glaciovolcanic Cave System—Mount Rainier, Washington, USA. J. Cave & Karst Studies, 2022, vol. 84(2), p. 66. DOI: 10.4311/2021ex0102.

11. Chkeir S., Anesiadou A., Mascitelli A. et al. Nowcasting extreme rain and extreme wind speed with machine learning techniques applied to different input datasets. Atmospheric Research, 2023, vol. 282, p. 106548. DOI: 10.1016/j.atmosres.2022.106548.

12. Dubinov A. E. Mathematical Tricks for Pseudopotentials in the Theories of Nonlinear Waves in Plasmas. Physics of Plasmas, 2022, vol. 29, p. 020901. DOI: 10.1063/5.0078573.

13. Dubinov A. E. Barometric Formula for Ultrarelativistic Degenerate Fermi-Gases. Astrophysics, 2020, vol. 63(4), p. 580. DOI: 10.1007/s10511-020-09660-1.

14. Moon J., Jeong Y. H., Addad Y. Design of Air-cooled Waste Heat Removal System with String Type Direct Contact Heat Exchanger and Investigation of Oil Film Instability. Nuclear Engineering and Technology, 2020, vol. 52(4), p. 734. DOI: 10.1016/j.net.2019.10.010.

15. Shu F. N. The Physics of Astrophysics. California: Univer. Sci. Books, 1992, vol. II, p. 246.

16. Tkachenko R. V, Korchagin V I., Zhmaylov B. B. Gravitatsionnaya neustoychivost' gazopylevykh okolo-yadernykh diskov blizkikh galaktik [Gravitational instability of gas and dust circumnuclear disks of nearby galaxies]. Vestnik Sankt-Peterburgskogo universiteta. Matematika. Mekhanika. Astronomiya [Bulletin of St. Petersburg University. Mathematics. Mechanics. Astronomy]. 2022, vol. 9(67), Iss. 3, pp. 561. DOI: 10.21638/spbu01.2022.316. (In Russian)

17. Patra N. N. Theoretical modelling of two-component molecular discs in spiral galaxies. Astronomy & Astrophysics, 2020, vol. 638, p. A66. DOI: 10.1051/0004-6361/201936483.

18. Bodmer D., Krenmayr M., Holzapfel F. Asymptotic tracking position control with active oscillation damping of a multibody Mars vehicle using two artificial augmentation

approaches. CEAS Space J., 2022, vol. 14, p. 125. DOI: 10.1007/s12567-021-00364-6.

19. Teyler R. Stroenie i evolyutsiya zvezd [Structure and evolution of stars]. Moscow: Mir Publ., 1973, 276 p. (In Russian)

20. Mizuno H. Formation of the Giant Planets. Progress of Theoretical Physics, 1980, vol. 64(2), p. 544. DOI: 10.1143/PTP.64.544.

21. Semenov A. O., Shved G. M. Upper Thermal Boundary Layer of Planetary Atmosphere: An Attempt of Developing a General Model. Icarus, 2008, vol. 194(1), p. 290. DOI: 10.1016/j.icarus.2007.08.040.

22. Shved G. M., Semenov A. O. The Standard Problem of Nonlocal Thermodynamic Equilibrium Radiative Transfer in the Rovibrational Band of the Planetary Atmosphere. Solar System Research, 2001, vol. 35, p. 212. DOI: 10.1023/A:1010478906172.

23. Bodenheimer P., Pollack J. B. Calculations of the Accretion and Evolution of Giant Planets: The Effects of Solid Cores. Icarus. 1986, vol. 67(3), p. 391. DOI: 10.1016/0019-1035(86)90122-3.

24. Pecnik B., Wuchterl G. Giant Planet Formation. A First Classification of Isothermal Protoplanetary Equilibria. Astronomy & Astrophysics. 2005, vol. 440, p. 1183. DOI: 10.1051/0004-6361:200500005.

Received: January 23, 2023 Accepted: February 14, 2023

Author's information:

Evgeniy N. BODUNOV — Dr. Sci. in Physics and Mathematics, Professor, Head of the Department "Physics"; evgeny.bodunov@inbox.ru Aleksey O. SEMENOV — PhD in Physics and Mathematics, Associate Professor; aleksem@mail.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.