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

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(→‎Общие сведения по теме)
(→‎Решение)
Строка 42: Строка 42:
 
где <math>F_{gr} = - G\frac{m^2}{r^2}</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_{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>F_{react} = G\frac{m^{2}a^{11}}{r^{13}}</math> реактивная составляющая.
  
 
а <math>\gamma</math> -гравитационная постоянная.
 
а <math>\gamma</math> -гравитационная постоянная.

Версия 14:36, 28 мая 2013

А.М. Кривцов > Теоретическая механика > Курсовые проекты 2013 > Моделирование образования Земли и Луны в газопылевом облаке с учётом солнечного притяжения (моделирование эллипсоида)
GazObl.jpg

Курсовой проект по Теоретической механике

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

Группа: 07 (20510)

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

Аннотация проекта

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

Постановка задачи

Смоделировать систему облака частиц в центральном поле массивного тела (Солнца). Оценить правдоподобность гипотезы о появлении Солнечной системы из газопылевого облака.

Общие сведения по теме

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

Решение

Начальное состояние, 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]\gamma[/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]

где

[math]\omega_0 = \sqrt{\frac{3 \pi GNm}{4R^3}}[/math]

составляющая возникающая из-за присутствия твердотельного вращения (частицы движутся, как твердое тело).

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

10 тысяч частиц, результат после 15 тысяч итераций

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

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

Проект "Земля - Луна"

См. также

Источник — «http://tm.spbstu.ru/?title=КП:_Моделирование_образования_Земли_и_Луны_в_газопылевом_облаке_с_учётом_солнечного_притяжения(моделирование_эллипсоида)&oldid=19809»