Редактирование: Устойчивость протопланетного облака системы "Земля - Луна" часть 2
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 1: | Строка 1: | ||
− | |||
− | |||
==Вторая модель== | ==Вторая модель== | ||
(Черновые тезисы) | (Черновые тезисы) | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | * | + | * Вопрос на который должна дать данная часть работы, можно сформулировать так: "Как изменится концентрация испарившихся молекул после прохождения через испаряющейся участок пылевого облака" |
* Работа состояла из двух частей: Компьютерного 2D моделирования эксперимента и поиска аналитического решения. | * Работа состояла из двух частей: Компьютерного 2D моделирования эксперимента и поиска аналитического решения. | ||
Строка 17: | Строка 9: | ||
* Эксперимент и аналитика дали хорошее совпадение результатов. Рисунки представлены ниже. | * Эксперимент и аналитика дали хорошее совпадение результатов. Рисунки представлены ниже. | ||
− | [[Файл: Nonrad.png|thumb|left| | + | [[Файл: Nonrad.png|thumb|left|500px|рис. 1. Прохождение "лучей" через неиспаряющееся облако.]] |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
+ | [[Файл: Radiation.png|thumb|right|500px|рис. 2. Прохождение "лучей" в испаряющемся облаке.]] | ||
+ | [[Файл: Summ.png|thumb|left|500px|рис. 3. Сумма двух вкладов.]] | ||
+ | [[Файл: lengh.png|thumb|right|500px|рис. 4. Эта зависимость от расстояния при заданной концентрации.]] | ||
== Об аналитическом решении== | == Об аналитическом решении== | ||
Строка 95: | Строка 32: | ||
Для начала-формальное определение: | Для начала-формальное определение: | ||
− | ''Пусть производится <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> | + | Сделаем важное допущение – произведение <math>np</math> сохраняет постоянное значение: <math>np=\lambda</math> Практически это допущение означает, что среднее число появления события в различных сериях испытаний (при разном <math>N</math>) остается неизменным.'' |
− | Очевидно, что вероятность того, что точка покроется одной окружностью есть отношение площади окружности к площади поверхности:<math>p=\frac{S}{H}</math>. Тогда для <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>. Случай ,когда точка не покрылась ни одной окружностью, т.е. случай, когда луч прошёл через цилиндр со сферами ни встретив ни одну на своём пути. | ||
Строка 113: | Строка 48: | ||
Поэтому из <math>N</math> лучей в среднем пройдёт | Поэтому из <math>N</math> лучей в среднем пройдёт | ||
− | <math> | + | <math>N_{ex}=N\cdot e^{-\lambda}=N\cdot e^{-\frac{n\cdot S}{H}}=N\cdot exp(-n\cdot S \cdot l) </math>.......................................(1) |
лучей. | лучей. | ||
− | |||
− | |||
Зависимость числа пройдённых лучей от концентрации в некотором объёме проиллюстрирована на рис. 1. | Зависимость числа пройдённых лучей от концентрации в некотором объёме проиллюстрирована на рис. 1. | ||
Строка 127: | Строка 60: | ||
Сперва откажемся от требования произвольности направления излучения. Пусть луч выходит из центра шарика в направлении основания <math>H</math>. Тогда для луча с шарика, находящегося на расстоянии <math>x</math> от поверхности <math>H</math>, вероятность прохожднения через оставшейся объём одного луча есть | Сперва откажемся от требования произвольности направления излучения. Пусть луч выходит из центра шарика в направлении основания <math>H</math>. Тогда для луча с шарика, находящегося на расстоянии <math>x</math> от поверхности <math>H</math>, вероятность прохожднения через оставшейся объём одного луча есть | ||
− | <math> | + | <math>N_j^{*}=e^{-\frac{n\cdot S\cdot x}{H\cdot L}}</math>.................................................................(2) |
− | </math> | ||
− | + | где <math>L</math> это длинна цилиндра. Соответственно, если излучается <math>I</math> лучей, то последнее выражение надо умножить на <math>I</math>. | |
− | Чтобы найти количество лучей вышедших из бесконечно тонкого объёма, параллельного | + | Чтобы найти количество лучей вышедших из бесконечно тонкого объёма, параллельного <math>H</math>, и находящегося на расстоянии <math>x</math>, очевидно необходимо (2) умножить на <math>N/L</math>. |
− | + | Для всех шариков в цилиндре: | |
− | <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>H</math> и вершиной в центре частички. |
− | + | Иллюстрацию этого закона можно видеть на рис. 3. | |
− | + | К сожалению, это выражение столь сложно, что интеграл "руками" не берётся. Поэтому проводилось только численное интегрирование. | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | <math> | + | <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>- интенсивность испарения. | |
− | |||
== Оценки для устойчивости== | == Оценки для устойчивости== | ||
Строка 186: | Строка 105: | ||
<math>F_v=\frac{\pi \nu v r_p^4 }{4 r^2} </math>..................................................................(6) | <math>F_v=\frac{\pi \nu v r_p^4 }{4 r^2} </math>..................................................................(6) | ||
− | <math>k=\frac{ | + | <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>k=1</math>, силы газодинамического отталкивания, полностью уравновешивают силу гравитационного притяжения (рис.5). | ||
Строка 311: | Строка 230: | ||
[[Файл:res2.png|800px|thumb|right|рис.7]] | [[Файл:res2.png|800px|thumb|right|рис.7]] | ||
− | [[Файл: | + | [[Файл:Newtime.png|800px|thumb|right|рис.8]] |
Строка 366: | Строка 285: | ||
<math>\nu=1.38\cdot 10^{-11} \frac{kg}{cm^2\cdot sek}</math> | <math>\nu=1.38\cdot 10^{-11} \frac{kg}{cm^2\cdot sek}</math> | ||
− | ==Другой подход к | + | ==Другой подход к от<math>ы</math>сканию условия равновесия протопланетного облака== |
Идея: [[А.М.Кривцов]] | Идея: [[А.М.Кривцов]] | ||
Строка 372: | Строка 291: | ||
Обозначения: | Обозначения: | ||
− | <math>V</math> | + | <math>V</math>-объём протооблока. |
− | <math>V_{\sigma}</math> | + | <math>V_{\sigma}</math>-объём частицы. |
− | <math>V_0</math> | + | <math>V_0</math>-суммарный объём всех частиц, т.е. объём конденсированной фазы всего вещества протооблака. |
− | <math>R</math> | + | <math>R</math>-радиус протооблака. |
− | <math>R_0</math> | + | <math>R_0</math>-радиус конденсированной фазы вещ-ва протооблака. |
− | <math>\rho</math> | + | <math>\rho</math>-радиус частицы. Частицы сферические. |
− | <math>\nu</math> | + | <math>\nu</math>-концентрация частиц. |
− | <math>r</math> | + | <math>r</math>-Отношение объёма облака к площади поверхности. Характеризует его форму. |
− | <math>\sigma</math> | + | <math>\sigma</math>- площадь поверхности частицы. |
− | <math>N</math> | + | <math>N</math>-количество частиц в протооблаке. <math>N=\nu\cdot V</math>. |
Рассмотрим разряженное протопланетное облако. Эффект экранирования возникает, когда между двумя рассматриваемыми частицами существует третья частица. Рассмотрим цилиндр- с площадью оснований, равных площади окружности, на которых построена сферческая частица | Рассмотрим разряженное протопланетное облако. Эффект экранирования возникает, когда между двумя рассматриваемыми частицами существует третья частица. Рассмотрим цилиндр- с площадью оснований, равных площади окружности, на которых построена сферческая частица | ||
Строка 414: | Строка 333: | ||
<math>k=\frac{3}{4}\frac{R_0^2}{R^2}\frac{R_0}{\rho}</math> | <math>k=\frac{3}{4}\frac{R_0^2}{R^2}\frac{R_0}{\rho}</math> | ||
− | И | + | И от сюда: |
<math>\frac{R}{R_0}=\sqrt{\frac{3}{4}\frac{R_0}{k \rho}}</math> | <math>\frac{R}{R_0}=\sqrt{\frac{3}{4}\frac{R_0}{k \rho}}</math> | ||
Строка 426: | Строка 345: | ||
Для того чтоб свести одну задачу к другой, нужно в формуле | Для того чтоб свести одну задачу к другой, нужно в формуле | ||
− | <math>n \cdot S \cdot L= 6 \left( \frac{r_p}{L}\right)^2 \left( \frac{R}{r_p}\right)^3</math>, под <math>L | + | <math>n \cdot S \cdot L= 6 \left( \frac{r_p}{L}\right)^2 \left( \frac{R}{r_p}\right)^3</math>, под <math>L</math>, подразумевать радиус протооблака. |
Тогда вместо 6, будет стоять 3. | Тогда вместо 6, будет стоять 3. | ||
− | <math>3 \left( \frac{r_p}{L | + | <math>3 \left( \frac{r_p}{L}\right)^2 \left( \frac{R}{r_p}\right)^3=-ln(k)=k^{*}</math> |
Проделаем элементарные алгебраические преобразования и получим ответ. | Проделаем элементарные алгебраические преобразования и получим ответ. | ||
Строка 442: | Строка 361: | ||
Разобранный пример позволяет связать искусственно подбираемый в '''Задаче 1''' малый параметр с параметрами протопланетного облака. | Разобранный пример позволяет связать искусственно подбираемый в '''Задаче 1''' малый параметр с параметрами протопланетного облака. | ||
− | == | + | ==Задача 4== |
− | + | '''Получить конкретные параметры протопланетного облака, которое будет жить достаточно долго.''' | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | : | + | Как оказалось, учёт взаимодействия каждой частицы с каждой, не привёл к удовлетворительным результатам, поэтому будем использовать радиус обрезания сил взаимодействия. Для этого введём самый большой радиус обрезания из возможных: пусть разность между силами радиационного отталкивания и гравитационного притяжения по модулю равна силе светового давления на расстоянии 1 а.е. |
− | :<math> | + | Как [http://ru.wikipedia.org/wiki/%C4%E0%E2%EB%E5%ED%E8%E5_%FD%EB%E5%EA%F2%F0%EE%EC%E0%E3%ED%E8%F2%ED%EE%E3%EE_%E8%E7%EB%F3%F7%E5%ED%E8%FF известно] давление солнечного света на перпендикулярную свету зеркальную поверхность, находящуюся в космосе в районе Земли равно примерно <math>F_s=4.5 \cdot 10^{-6} \frac{N}{m^2}</math> |
− | + | <math>F_g-F_v=F_s</math>.........................(9) | |
− | + | Предположим, что параметры облака, таковы, что в отсутствии экранирования <math>F_g=F_v</math>. Необходимые соотношения обсуждались выше. | |
− | + | При "включении" экранирования, условие (9) запишется, как | |
− | |||
− | |||
− | --- | + | <math>F_g-F_v=F_g-F_g\cdot(1-k)=F_g\cdot k=F_s</math> |
− | |||
− | <math> | + | <math>F_g=\frac{16 \pi^2 G\rho_p^2 r_p^6}{9 r^2}</math> |
− | + | Поэтому | |
− | <math> | + | <math>r=\sqrt{\frac{16 \pi^2 G\rho_p^2 r_p^6 k}{9 F_s}}</math> |
− | + | ==Ссылки по теме== | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
* [http://www.spacetelescope.org/news/heic0917/ Снимки зарождающихся планетных систем в туманности Ориона, полученные космическим телескопом Хаббл] | * [http://www.spacetelescope.org/news/heic0917/ Снимки зарождающихся планетных систем в туманности Ориона, полученные космическим телескопом Хаббл] | ||
* [[А. Мурачев]]. '''Некоторые замечания по модели образования системы Земля-Луна в результате ротационного коллапса газопылевого облака (черновые наброски)'''. Скачать презентацию, pptx: [http://195.209.230.53:8088/Presentations/Student_seminar/2011_10_21_Murachev.pptx 2741 kb] | * [[А. Мурачев]]. '''Некоторые замечания по модели образования системы Земля-Луна в результате ротационного коллапса газопылевого облака (черновые наброски)'''. Скачать презентацию, pptx: [http://195.209.230.53:8088/Presentations/Student_seminar/2011_10_21_Murachev.pptx 2741 kb] | ||
− | + | ||
− | |||
− | |||
[[Category: Проект "Земля - Луна"]] | [[Category: Проект "Земля - Луна"]] | ||
[[Category: Студенческие проекты]] | [[Category: Студенческие проекты]] |