Устойчивость протопланетного облака системы "Земля - Луна" часть 2

Вторая модель

(Черновые тезисы)

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

• Работа состояла из двух частей: Компьютерного 2D моделирования эксперимента и поиска аналитического решения.

• Эксперимент и аналитика дали хорошее совпадение результатов. Рисунки представлены ниже.

рис. 1. Прохождение "лучей" через неиспаряющееся облако.

рис. 2. Прохождение "лучей" в испаряющемся облаке.

рис. 3. Сумма двух вкладов.

рис. 4. Эта зависимость от расстояния при заданной концентрации.

Об аналитическом решении

Чтобы найти закон распространения радиационного излучения в облаке, воспользуемся простой моделью. Некоторый объём пространства представим в виде цилиндра с основанием [math]H[/math] и длиной [math]L[/math]. В этом цилиндре находятся частицы-сферы радиуса r.

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

Наша задача ответить на вопрос: сколько в среднем лучей достигнет противоположной стенки цилиндра?

Эта задача равносильна следующей: На некоторой поверхности [math]H[/math] нанесена в случайном месте точка. На поверхность набрасываются случайным образом [math]N[/math] окружностей площадью [math]S[/math] [math](S=\pi r^2)[/math]. Какова вероятность того, что точка останется непокрытой? Будем решать её.

Ответ легко получить из закона распределения Пуассона.

Для начала-формальное определение:

Пусть производится [math]n[/math] независимых испытаний. Если число испытаний [math]N[/math] достаточно велико, а вероятность появления события [math]А[/math] в каждом испытании мало ([math]p\in(0,1)[/math]), то вероятность появления события [math]А[/math] [math]k[/math] раз находится следующим образом:

[math]P_n(k)=\frac{\lambda^k e^{-\lambda}}{k!}[/math]

Сделаем важное допущение – произведение [math]np[/math] сохраняет постоянное значение: [math]np=\lambda[/math] Практически это допущение означает, что среднее число появления события в различных сериях испытаний (при разном [math]N[/math]) остается неизменным.

Очевидно, что вероятность того, что точка покроется одной окружностью есть отношение площади окружности к площади поверхности:[math]p=\frac{S}{H}[/math]. Тогда для [math]N[/math] окружностей эта вероятность: [math]\lambda=\frac{n\cdot S}{H}[/math]

Вероятность, что точка будет покрыта ровна [math]к[/math] окружностями из [math]n[/math] даёт формула Пуассона.

Нас интересует случай, когда [math]k=0[/math]. Случай ,когда точка не покрылась ни одной окружностью, т.е. случай, когда луч прошёл через цилиндр со сферами ни встретив ни одну на своём пути.

Подставляя [math]k=0[/math], в формулу Пуассона, находим [math]P_n(0)=e^{-\lambda}[/math].

Поэтому из [math]N[/math] лучей в среднем пройдёт

[math]N_{ex}=N\cdot e^{-\lambda}=N\cdot e^{-\frac{n\cdot S}{H}}[/math] .......................................(1)

лучей.

Зависимость числа пройдённых лучей от концентрации в некотором объёме проиллюстрирована на рис. 1.

Теперь рассмотрим другую ситуацию.

Тот же объём и те же шарики. Но теперь радиационное излучение не падает на стенки цилиндра, а каждый шарик излучает в произвольном направлении луч. Спрашивается, сколько лучей пройдёт через основание цилиндра?

Сперва откажемся от требования произвольности направления излучения. Пусть луч выходит из центра шарика в направлении основания [math]H[/math]. Тогда для луча с шарика, находящегося на расстоянии [math]x[/math] от поверхности [math]H[/math], вероятность прохожднения через оставшейся объём одного луча есть

[math]N_j^{*}=e^{-\frac{n\cdot S\cdot x}{H\cdot L}}[/math].................................................................(2)

где [math]L[/math] это длинна цилиндра. Соответственно, если излучается [math]I[/math] лучей, то последнее выражение надо умножить на [math]I[/math].

Чтобы найти количество лучей вышедших из бесконечно тонкого объёма, параллельного [math]H[/math], и находящегося на расстоянии [math]x[/math], очевидно необходимо (2) умножить на [math]N/L[/math].

Для всех шариков в цилиндре:

[math]N_j=\frac{N\cdot I}{L}\int_{0}^{L} e^{-\frac{n\cdot S\cdot x}{H\cdot L}} dx =\frac{N\cdot I\cdot H}{n\cdot S}(1-e^{-\frac{n\cdot S}{H}} ) [/math].................................................(3)

Для того, чтобы учесть случайность направления луча, следует выражение (3) умножить на вероятность того, что луч попадёт в угловой сегмент с основанием [math]H[/math] и вершиной в центре частички.

Иллюстрацию этого закона можно видеть на рис. 3.

К сожалению, это выражение столь сложно, что интеграл "руками" не берётся. Поэтому проводилось только численное интегрирование.

[math]N_j=\int_0^l I\cdot n\cdot \left(2\cdot arctg(\frac{H}{x})\cdot H- x\cdot ln(1+\frac{H^2}{x^2})\right)\cdot exp(-n \cdot S \cdot x)[/math]......................................................(4)

, где [math]n[/math]-концентрация пылинок

[math]I[/math]- интенсивность испарения.

Оценки для устойчивости

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

Сила гравитационного притяжения:

[math]F_g=G\frac{m^2}{r^2}=\frac{16 \pi^2 G\rho_p^2 r_p^6}{9 r^2}[/math].......................................(5)

Сила радиационного отталкивания в пустом пространстве есть

[math]F_v=\frac{N m_0 v r_p^2}{4 r^2} [/math]

При выводе этой формулы, столкновения между молекулами не учитывалось.

Если ввести

[math]Nm_0=4\pi r_p^2 \nu[/math]- количество массы испаряющееся с пылинки за единицу времени (величина независящая от размеров частицы), то предыдущая формула перепишется, как

[math]F_v=\frac{\pi \nu v r_p^4 }{4 r^2} [/math]..................................................................(6)

[math]k=\frac{f_v}{f_g}= \frac {9\nu v}{16\pi G r_p^2 \rho^2} [/math]

При [math]k=1[/math], силы газодинамического отталкивания, полностью уравновешивают силу гравитационного притяжения (рис.5).

Для пространства заполненного другими частицами, создающими экранирующий эффект

[math]F_v=\frac{\pi \nu v r_p^4 exp(-n \cdot \pi r_p^2 \cdot r)}{r^2}[/math]..................................................................(7)

,Где [math]n[/math]-концентрация пылинок в рассматриваемом объёме.

При равенстве этих сил облако будет находиться в равновесии.

Между тем, понятно, что при любой, отличной от нуля концентрации и при любом, отличном от нуля размере пылинок будет присутствовать эффект экранирования. Взаимодействие далёких частиц будут всё более отходить от закона обратных квадратов.

Будем считать, что облако может находиться в равновесии, если отношение сил радиационного отталкивания и гравитационного притяжения меньше одного процента.

[math]k=\frac{F_v}{F_g}\le 0.99[/math] ....................................................(8)

[math]k=\frac{F_v}{F_g}=\frac{9\pi\nu v \cdot exp(-n \cdot \pi r_p^2 \cdot r)}{16 \pi^2 G r_p^2\rho^2 }[/math]

Среднею скорость примем равной среднему значению модуля скорости идеального газа:

[math]v=\sqrt{\frac{8RT}{\pi}}[/math]

Радиус облака равен среднему расстоянию между Землёй и Луной ( [math]L=3,84\cdot 10^{10}[/math] см ).

Пусть на расстоянии [math]2L[/math] выполняется условие (8).

Если принять [math]k=\frac{F_v}{F_g}=\frac{9\pi\nu v }{16 \pi^2 G r_p^2\rho^2 } =1 [/math], то значит , то [math] exp(-n \cdot \pi r_p^2 \cdot 2L)=0.99[/math]

Следовательно:

[math] n \cdot \pi r_p^2 \cdot 2L=-ln(0.99)=0.01 =\gt n \cdot \pi r_p^2=2.6\cdot 10^{-13} [/math]

Тогда построим зависимость концентрации от площади, при условии сохранения постоянным их произведения.

рис.5

Радиус порядка [math]10^{-6}[/math] см, всего в 100 раз больше характерных размеров молекул водяного пара.

На рис. 6 приведена зависимость времени жизни частицы от радиуса. Интенсивность испарения расчитыалась из предпосылки исходного равновесия облака.

рис.6

На рис. 7 Масса протопланетного облака от размеров частиц. То что она растёт линейно, не удивительно. Ибо масса частиц хоть и увеличивается как куб радиуса, зато их концентрация падает, как квадрат радиуса.

Значение Итоговой Массы в [math]10^{14}/math\gt много меньше массы Земли (\lt math\gt 10^{24}[/math])

С другой стороны понятно, что требовать, чтобы каждая частица взаимодействовала с каждой, глупо, так как на больших расстояниях случайные флуктуации космологических факторов (давление света, плотность газа, случайные грав. сгустки.) будут оказывать на частицы большее воздействие. Поэтому необходимо ввести радиус обрезания грав. сил.

Чтоб решить эту проблему пойдём от обратного. Предположим помимо устойчивости облака то, что его масса порядка 1/10 массы системы Земля-Луна [math]6.05 \cdot 10^{27}[/math] грамм, то есть [math]6\cdot 10^{26}[/math] грамм.

Отсюда, произведение концентрации на массу частицы есть отношение Массы облака к его объёму.

[math]n\cdot m=M/V=\frac{6\cdot 10^{26}}{4/3 \pi R^3}[/math]

Подстановка значений даёт:

[math]n\cdot m=1.43\cdot 10^{-4}\frac{g}{sm^3}[/math]...............................................(9)

Задача состоит в том, чтобы найти такой радиус обрезания [math]l[/math] рассматриваемых сил, чтобы выполнялось (9). Коэффециент [math]k=0.99[/math] оставим без изменения.

[math]n\cdot m=\frac{-ln(0.99)}{\pi r_p^2 \cdot l}\cdot \frac{4\pi r_p^3}{3}=\frac{-4 ln(0.99)\cdot r_p}{3\cdot l}[/math]

Таким образом:

[math]l=\frac{-ln(0.99)\cdot 4 \cdot r_p}{3\cdot 1.43\cdot 10^{-4}} \rightarrow l=93.24 r [/math]

Что получилось:

рис.7

рис.8

Ссылки по теме

• Снимки зарождающихся планетных систем в туманности Ориона, полученные космическим телескопом Хаббл
• А. Мурачев. Некоторые замечания по модели образования системы Земля-Луна в результате ротационного коллапса газопылевого облака (черновые наброски). Скачать презентацию, pptx: 2741 kb