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

Текущая версия на 14:52, 28 мая 2013

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

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

Группа: 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 тысяч итераций

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

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

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

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

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

• Справочная информация
• Курсовые проекты 2012
• Курсовые проекты 2013