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

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Общие сведения по теме)
(Обсуждение результатов и выводы)
 
(не показано 8 промежуточных версий этого же участника)
Строка 18: Строка 18:
  
 
== Общие сведения по теме ==
 
== Общие сведения по теме ==
Внесение массивного центрального тела влечет за собой необходимость пересчета начальных данных, в частности угловой скорости вращения облака <math>\omega</math>. Это имеет важную роль в дальнейшем исследовании, так как интуитивно понятно, если задать слишком маленькой, получится так, что все частицы "упадут" на Солнце. Выбрав её слишком большой, то получится так, что все частицы разлетятся. Значит, существует интервал значений угловой скорости, при которых облако будет стабильно вращаться и появится возможность образования кластеров.
+
Внесение массивного центрального тела влечет за собой необходимость пересчета начальных данных, в частности угловой скорости вращения облака <math>\omega</math>. Это имеет важное значение в дальнейшем исследовании, так как интуитивно понятно, что, если задать её слишком маленькой, получится так, что все частицы "упадут" на Солнце. Если выбрать её слишком большой, то получится так, что все частицы разлетятся. Значит, существует интервал значений угловой скорости, при которых облако будет стабильно вращаться и появится возможность образования кластеров. Моделирование при различных угловых скоростях подтверждает эти предположения.
  
 
== Решение ==
 
== Решение ==
Строка 25: Строка 25:
 
Начальными условиями в данной задаче являются следующие характеристики системы:
 
Начальными условиями в данной задаче являются следующие характеристики системы:
  
<math>N</math> - число частиц
+
<math>N</math> число частиц
  
<math>d_o</math> - среднее расстояние между частицами
+
<math>d_o</math> среднее расстояние между частицами
  
<math>m</math> масса частицы
+
<math>m</math> масса частицы
  
<math>\frac{mN}{M}</math> отношение массы системы, к массе массивного тела
+
<math>\frac{mN}{M}</math> —  отношение массы системы, к массе массивного тела
  
<math>\frac{V_{xy}}{V}</math> отношение рандомной составляющей скорости к линейной скорости вращения
+
<math>\frac{V_{xy}}{V}</math> отношение рандомной составляющей скорости к линейной скорости вращения
  
  
Строка 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>G</math> -гравитационная постоянная.
  
 
<math>m</math> - масса частицы
 
<math>m</math> - масса частицы
Строка 82: Строка 82:
 
[[Файл:ResultLapinRuslan.bmp|thumb|350px|10 тысяч частиц, результат после 15 тысяч итераций ]]
 
[[Файл:ResultLapinRuslan.bmp|thumb|350px|10 тысяч частиц, результат после 15 тысяч итераций ]]
 
                                          
 
                                          
В результате моделирования получилось стабильное облако частиц, которое вращается вокруг солнца. Начали образовываться небольшие кластеры, что свидетельствует о правильности модели и метода моделирования. Далее возможно исследование системы в зависимости от начальных параметров. Например увеличение числа частиц, изменение скоростей частиц, и наблюдение как будут меняться результаты.
+
В результате моделирования получилось стабильное облако частиц, которое вращается вокруг солнца. Начали образовываться небольшие кластеры, что свидетельствует о правильности модели и метода моделирования. Далее возможно исследование системы в зависимости от начальных параметров. Например, увеличение числа частиц, изменение скоростей частиц, и наблюдение за тем, как будут меняться результаты.
 +
 
 +
Параллельно моделированию был проведен [[КП: Динамика молекулы (расчет)|аналитический анализ данной задачи]]. В дальнейшем планируется подведение результатов моделирования к аналитическим.
  
 
== Ссылки по теме ==  
 
== Ссылки по теме ==  

Текущая версия на 14:52, 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]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]

где

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

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

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

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

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

Параллельно моделированию был проведен аналитический анализ данной задачи. В дальнейшем планируется подведение результатов моделирования к аналитическим.

Ссылки по теме[править]

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

См. также[править]