Редактирование: КП: Моделирование образования Земли и Луны в газопылевом облаке с учётом солнечного притяжения(моделирование эллипсоида)

[[А.М. Кривцов]] > [[Теоретическая механика: физико-механический факультет|Теоретическая механика]] > [[Курсовые проекты ТМ 2013|Курсовые проекты 2013]] > '''Моделирование образования Земли и Луны в газопылевом облаке с учётом солнечного притяжения (моделирование эллипсоида)''' <HR>

[[Файл:GazObl.jpg|справа|350px]]

'''''Курсовой проект по [[Теоретическая механика: физико-механический факультет|Теоретической механике]]'''''

'''Исполнитель:''' [[Лапин Руслан]]

'''Группа:''' [[Группа 07|07]] (20510)

'''Семестр:''' весна 2013

== Аннотация проекта ==
В основу данного курсового проекта положен [[Проект "Земля - Луна"|проект "Земля-Луна"]], выполненный [[А.М.Кривцов]]ым и [[Ле-Захаров А.А.|А.А. Ле-Захаровым]]. Исследуется возможное появление Солнечной системы из газопылевого облака путем моделирования облака (эллипсоида) частиц в силовом поле массивного центрального тела (Солнца), которое расположено в центре системы.

== Постановка задачи ==
Смоделировать систему облака частиц в центральном поле массивного тела (Солнца). Оценить правдоподобность гипотезы [http://elementy.ru/trefil/21175| о появлении Солнечной системы из газопылевого облака].

== Общие сведения по теме ==
Внесение массивного центрального тела влечет за собой необходимость пересчета начальных данных, в частности угловой скорости вращения облака <math>\omega</math>. Это имеет важное значение в дальнейшем исследовании, так как интуитивно понятно, что, если задать её слишком маленькой, получится так, что все частицы "упадут" на Солнце. Если выбрать её слишком большой, то получится так, что все частицы разлетятся. Значит, существует интервал значений угловой скорости, при которых облако будет стабильно вращаться и появится возможность образования кластеров. Моделирование при различных угловых скоростях подтверждает эти предположения.

== Решение ==
[[Файл:NachUslLapiRuslan.jpg|thumb|350px|Начальное состояние, 10 тысяч частиц]]

Начальными условиями в данной задаче являются следующие характеристики системы:

<math>N</math> — число частиц

<math>d_o</math> — среднее расстояние между частицами

<math>m</math> — масса частицы

<math>\frac{mN}{M}</math> —  отношение массы системы, к массе массивного тела

<math>\frac{V_{xy}}{V}</math> — отношение рандомной составляющей скорости к линейной скорости вращения


Без учета солнца на частицу действуют следующие силы:

<math>\overline{F} = \overline{F}_{gr} + \overline{F}_{dis} + \overline{F}_{react}</math>

где <math>F_{gr} = - G\frac{m^2}{r^2}</math> гравитационная составляющая.

<math>F_{dis} = -G\frac{m^{2}a^{11}\beta r'}{r^{14}}</math> диссипативная составляющая.

<math>F_{react} = G\frac{m^{2}a^{11}}{r^{13}}</math> реактивная составляющая.

а <math>G</math> -гравитационная постоянная.

<math>m</math> - масса частицы

<math>a</math> - равновесное расстояние

<math>r</math> - расстояние между частицами

<math>\beta</math> - коэффициент диссипации

При внесении массивного тела массы <math>M</math>, возникает центральное поле и на частицы начинает действовать сила 

<math>F = G\frac{mM}{R^2}</math>, где <math>m</math> - масса частицы, <math>R</math> - расстояние до частицы.

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

<math> m\omega^2 R = G\frac{mM}{R^2}</math>

То есть угловая скорость от появления центрального тела, должна быть 

<math>\omega = \sqrt{G\frac{M}{R^3}}</math>

А тогда учитывая еще и твердотельное вращение получаем искомую угловую скорость
 
<math>\omega = \sqrt{G\frac{M}{R^3}} + \sqrt{\frac{3 \pi GNm}{4R^3}}</math>

== Обсуждение результатов и выводы ==

[[Файл:ResultLapinRuslan.bmp|thumb|350px|10 тысяч частиц, результат после 15 тысяч итераций ]]
                                        
В результате моделирования получилось стабильное облако частиц, которое вращается вокруг солнца. Начали образовываться небольшие кластеры, что свидетельствует о правильности модели и метода моделирования. Далее возможно исследование системы в зависимости от начальных параметров. Например, увеличение числа частиц, изменение скоростей частиц, и наблюдение за тем, как будут меняться результаты. 

Параллельно моделированию был проведен [[КП: Динамика молекулы (расчет)|аналитический анализ данной задачи]]. В дальнейшем планируется подведение результатов моделирования к аналитическим.

== Ссылки по теме == 
[[Проект "Земля - Луна"]]
== См. также ==

* [[Справка:Содержание|Справочная информация]]
* [[Курсовые проекты ТМ 2012|Курсовые проекты 2012]]
* [[Курсовые проекты ТМ 2013|Курсовые проекты 2013]]


[[Category: Студенческие проекты]]
@@ Строка 71: / Строка 71: @@
 <math>\omega = \sqrt{G\frac{M}{R^3}} + \sqrt{\frac{3 \pi GNm}{4R^3}}</math>
-где
-<math>\omega_0 = \sqrt{\frac{3 \pi GNm}{4R^3}}</math>
-составляющая возникающая из-за присутствия твердотельного вращения (частицы движутся, как твердое тело).
 == Обсуждение результатов и выводы ==