Редактирование: КП: Динамика молекулы (моделирование)
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 1: | Строка 1: | ||
[[А.М. Кривцов]] > [[Теоретическая механика: физико-механический факультет|Теоретическая механика]] > [[Курсовые проекты ТМ 2013|Курсовые проекты 2013]] > '''Моделирование образования Земли и Луны в газопылевом облаке с учётом солнечного притяжения (моделирование)''' <HR> | [[А.М. Кривцов]] > [[Теоретическая механика: физико-механический факультет|Теоретическая механика]] > [[Курсовые проекты ТМ 2013|Курсовые проекты 2013]] > '''Моделирование образования Земли и Луны в газопылевом облаке с учётом солнечного притяжения (моделирование)''' <HR> | ||
− | + | ||
'''''Курсовой проект по [[Теоретическая механика: физико-механический факультет|Теоретической механике]]''''' | '''''Курсовой проект по [[Теоретическая механика: физико-механический факультет|Теоретической механике]]''''' | ||
Строка 12: | Строка 12: | ||
== Аннотация проекта == | == Аннотация проекта == | ||
− | Данный проект является продолжением работы над проектом "Земля-Луна" | + | Данный проект является продолжением работы над проектом "Земля-Луна".Исследуется поведение газопылевого облака в гравитационном поле, создаваемом Солнцем. Уже получены результаты для облака, имеющего форму эллипсоида. Результатом является образование двойной системы путём одновременного возникновения двух небесных тел, одно из которых в несколько раз больше другого по размеру. Для этого случая были проведены исследования, которые выявили наиболее оптимальные исходные данные, при которых образование двойной системы происходит с большой вероятностью. |
== Постановка задачи == | == Постановка задачи == | ||
Строка 20: | Строка 20: | ||
== Общие сведения по теме == | == Общие сведения по теме == | ||
− | Для | + | Для этого случая необходимо учитывать многое из того, что не учитывается в случае, когда облако задается в виде эллипсоида. Скорость твердотельного вращения должна считаться исходя из факта влияния на облако как самих частиц,его составляющих, так и центрального тела. Таким образом, задав скорость вращения облака, можно добиться того, что частицы пыли не будут падать в центр или разлетаться. |
− | |||
− | |||
'''Облако в виде тора в поле центрального тела.''' | '''Облако в виде тора в поле центрального тела.''' | ||
Строка 29: | Строка 27: | ||
Также необходимо учитывать тот факт, что из-за более сложной геометрии тора по сравнению с геометрией эллипсоида, понять где и как будут образовываться необходимые кластеры гораздо сложнее. В случае с эллипсоидом, можно было предполагать образование системы в центральной части, тогда как для тора оно может произойти где угодно. Причем в данном случае может получиться так, что образуется несколько кластеров, подобных друг другу. Таким образом, для вероятности образования двойной системы, необходимо увеличивать толщину тора. | Также необходимо учитывать тот факт, что из-за более сложной геометрии тора по сравнению с геометрией эллипсоида, понять где и как будут образовываться необходимые кластеры гораздо сложнее. В случае с эллипсоидом, можно было предполагать образование системы в центральной части, тогда как для тора оно может произойти где угодно. Причем в данном случае может получиться так, что образуется несколько кластеров, подобных друг другу. Таким образом, для вероятности образования двойной системы, необходимо увеличивать толщину тора. | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
== Решение == | == Решение == | ||
Строка 56: | Строка 40: | ||
Если учитывать, что угловая скорость твердотельного вращения без центрального тела имеет вид | Если учитывать, что угловая скорость твердотельного вращения без центрального тела имеет вид | ||
− | <math> \omega_{s0} =\sqrt{\frac{3 \pi GNm_{0}}{4R_{0}^3}}</math>, где <math>N</math> | + | <math> \omega_{s0} =\sqrt{\frac{3 \pi GNm_{0}}{4R_{0}^3}}</math>, где <math>N</math> - число частиц |
то получим, что искомая уловая скорость твердотельного вращения будет иметь вид | то получим, что искомая уловая скорость твердотельного вращения будет иметь вид | ||
Строка 63: | Строка 47: | ||
Задав отношение <math>\frac{\omega_{0}}{\omega_{s}}</math>, получаем значение угловой скорости облака <math>\omega_{0}</math> | Задав отношение <math>\frac{\omega_{0}}{\omega_{s}}</math>, получаем значение угловой скорости облака <math>\omega_{0}</math> | ||
− | Это отношение угловых скоростей позволяет | + | Это отношение угловых скоростей позволяет сохронить подобие при изменении числа частиц. |
− | Зная <math>\omega_{0}</math> можно найти скорость частицы <math>V_{0} = [\omega,r]</math>,где <math>r</math> | + | Зная <math>\omega_{0}</math> можно найти скорость частицы <math>V_{0} = [\omega,r]</math>,где <math>r</math> - расстояние до центра. |
+ | Результат работы программы при: | ||
+ | <math>\frac{\omega_{0}}{\omega_{s}}</math> = 1.05 | ||
− | |||
− | |||
<math>N</math> = 10000 | <math>N</math> = 10000 | ||
[[Файл: Prog2.png]] | [[Файл: Prog2.png]] | ||
− | Видим 12 кластеров, движущиеся по орбитам. Этот результат получен после 24804 итераций. То есть | + | Видим 12 кластеров, движущиеся по орбитам. Этот результат получен после 24804 итераций. То есть пи данном отношении угловых скоростей облако устойчиво. |
== Обсуждение результатов и выводы == | == Обсуждение результатов и выводы == | ||
− | |||
− | |||
== Ссылки по теме == | == Ссылки по теме == | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
== См. также == | == См. также == |