Редактирование: Устойчивость протопланетного облака системы "Земля - Луна"
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 77: | Строка 77: | ||
==Случай дискообразного протопланетного облака== | ==Случай дискообразного протопланетного облака== | ||
+ | Выделим в облаке окружность с центром в точке (0,0) некоторого радиуса <math>r</math>. Концентрация молекул газа, а следовательно и оказываемое ими давление в силу симметрии будет одинаковым. Найдём эту концентрацию: Сначала учтём вклад "внутренней" области выделенной окружности. Представим ещё одну сферу радиуса <math>a</math>, где <math>a<r</math> и найдём влияние этой сферы на бесконечно маленький объём <math>G</math> на сфере радиуса <math>r</math>. Для этого, очевидно надо воспользоваться формулой [1], где <math>r</math> пробегает по всем расстояниям (с каждой точки поверхности сферы до G). | ||
− | + | <math>L</math>-расстояние от точки <math>G</math> до точек правой полуокружности. | |
− | + | <math>\alpha</math>-угол между радиус-вектором <math>r</math> и радиусом <math>a</math> | |
− | + | Смотри рисунок '''internal.''' | |
+ | [[Файл: internal.jpg|thumb|right|540px|internal]] | ||
+ | |||
+ | Легко сообразить для левой полуокружности: | ||
+ | |||
+ | <math>L_1=\sqrt{r^2+a^2-2ra\cos (\alpha+3\frac\pi2)}=\sqrt{r^2+a^2+2ra\cos \alpha}</math> | ||
+ | |||
+ | и | ||
+ | |||
+ | <math>L_2=\sqrt{r^2+a^2-2ra\cos (\frac \pi2 +\alpha) }=\sqrt{r^2+a^2-2ra\cos (\alpha) }</math><math>Вставьте сюда формулу</math> | ||
+ | |||
+ | для правой. | ||
+ | |||
+ | Интегрирование по <math>L</math> сводится к интегрированию по <math>\alpha</math> от 0 до <math>\frac \pi 2</math> для "верхнего левого" и "нижнего правого" сектора. А потом результат удваивается. | ||
+ | |||
+ | |||
+ | <math>n(r)_{internal}=2\int_0^r da \left( | ||
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2-2ra\cos \alpha}}+ | ||
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2+2ra\cos \alpha}}\right)</math> | ||
+ | |||
+ | |||
+ | Теперь можно оценить вклад внешних слоёв на наш выделенный объём G . | ||
+ | Смотри рисунок '''external'''. Можно заметить, что тут всё то же самое, что и в предыдущем примере, лишь точка <math>G</math> <<переехала>> с одного конца отрезка на другой, а <math>a</math> и <math>r</math> поменялись местами. Поэтому повторяя предыдущие рассуждения, можно записать | ||
+ | [[Файл: external.jpg|thumb|right|540px|external]] | ||
+ | |||
+ | <math>n(r)_{internal}=2\int_r^R da \left( | ||
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2-2ra\cos \alpha}}+ | ||
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2+2ra\cos \alpha}}\right)</math> | ||
+ | |||
+ | |||
+ | Тогда влияние вклад в концентрацию от испарений всего облака будет | ||
+ | |||
+ | |||
+ | |||
+ | <math>n(r)=n(r)_{external}+n(r)_{internal}=</math> | ||
− | |||
− | |||
− | |||
− | + | <math>2\int_0^R da \left( | |
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2-2ra\cos \alpha}}+ | ||
+ | \int_0^{\frac \pi 2 } d\alpha | ||
+ | \frac{\dot N\cdot n(a)^*}{4\pi D\sqrt{r^2+a^2+2ra\cos \alpha}}\right)</math> | ||
− | + | [[Файл:Points.png|600px|thumb|right|concentration]] | |
− | + | Значение концентрации пылинок | |
− | + | <math>w(a)=\sqrt{1-\frac{a^2}{R^2}}</math> | |
− | + | Задача свелась к вычислению данных интегралов. | |
− | |||
− | |||
− | <math> | + | Но этот интеграл не берётся в элементарных функциях. Поэтому он был вычислен численно. На рисунке concentration, изображен результат численного интегрирования. Здесь <math>R=1</math>, и константа перед интегралом тоже взята за 1 (Точный вид этой константы зависит от константы диффузии <math>D</math>). |
+ | Что удивительно, при <math>r=R</math>, значение интеграла не равно нулю. Это можно объяснить тем, что у реального протопланетного диска нет чёткой границы. | ||
==Уравнение равновесия.== | ==Уравнение равновесия.== | ||
Строка 119: | Строка 158: | ||
==Испарение пылинок в вакуум== | ==Испарение пылинок в вакуум== | ||
− | Интенсивность испарения | + | Интенсивность испарения [<math>g/cm^2\cdot sek</math>] определяется формулой Ленгмюра. |
<math> | <math> | ||
Строка 132: | Строка 171: | ||
Эта формула выведена для абсолютного вакуума, поэтому реальная скорость испарения в космическом пространстве будет меньше расчётной. | Эта формула выведена для абсолютного вакуума, поэтому реальная скорость испарения в космическом пространстве будет меньше расчётной. | ||
− | |||
− | |||
− | |||
<math> | <math> | ||
− | \dot m =4\pi r^ | + | \dot m =4\pi r^3 \nu => \left(\frac{4}{3}\pi r^3 \rho \right)'=4\nu\pi r^3 => \dot r = \frac{\nu}{\rho} |
</math> | </math> | ||
Строка 167: | Строка 203: | ||
Теперь остался главный вопрос о концентрации и скорости пылинок. | Теперь остался главный вопрос о концентрации и скорости пылинок. | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
== См. также == | == См. также == | ||
− | * [[Устойчивость протопланетного облака системы "Земля - Луна" часть | + | * [[Устойчивость протопланетного облака системы "Земля - Луна" часть 2]]. (2012 г.) |
− | |||
[[Category: Проект "Земля - Луна"]] | [[Category: Проект "Земля - Луна"]] | ||
[[Category: Студенческие проекты]] | [[Category: Студенческие проекты]] |