липання [2, 9]. Пропоновані можливі механізми злипання частинок пилу: ван дер Ваальсових сили; різні типи В«радіаційногоВ» спікання [2], стор 413; ефект холодної зварювання [9], стр. 139 та інші. Чи відбудеться злипання або дроблення частинок при зіткненні залежить від їх відносної швидкості. За даними робіт [3, 7] в протопланетному диску Сонця частинки досягають розподілу за розмірами, в якому є і дрібні частинки розміром близько 1 мкм, що підтримують високу непрозорість речовини диска, і великі близько 1 см. Середній розмір пилових частинок становить кілька десятків мікрометрів.
Так як в даних дослідженнях вивчалася еволюція протопланетного диска на її початковому етапі, то було прийнято, що пилові частинки однорідно розподілені за обсягом протопланетного диска, а їх розмір не сильно відрізнявся від розміру пилових частинок в протосонячній туманності.
Для опису поведінки протопланетного диска необхідно знати рівняння стану його середовища.
2.1 Рівняння стану середовища протопланетного диска
Отже, припустимо, що початковий склад середовища протопланетного диска має розподіл щільності пиловий компоненти, близьке до однорідного, і вміст пилу за масою не перевищує декількох відсотків. За цих умов можна показати, що усереднені параметри даної середовища з достатньою точністю описуються рівнянням стану ідеального газу [15]. p align="justify"> Так, за даними роботи [3], якщо на пилові частинки припадає по масі близько 1,5% речовини сонячного складу, то для такого середовища молекулярний вага дорівнює 2,53, а показник адіабати - 1, 43. Опис протопланетної хмари у наближенні ідеального газу дає досить точну картину поведінки його деяких середніх характеристик, а саме тих, які локально визначаються газової компонентою, навіть у тому випадку, коли пилова компонента конденсується і перетворюється у тверду речовину. br/>
3. Чисельна двовимірна модель протопланетного газопилового диска
Описати з достатньою точністю еволюцію протопланетного диска Сонця можна тільки за допомогою адекватних багатовимірних чисельних моделей. Побудована чисельна двовимірна газодинамічна модель, яка відповідає осесиметричними руху газового середовища в гравітаційному полі. Ця модель може використовуватися на початковій стадії еволюції протопланетного диска, коли пилова компонента розподілена відносно рівномірно за обсягом диска. br/>
3.1 Рівняння газової динаміки у формі законів збереження для криволінійної системи координат
Використовувана чисельна методика є варіантом методики для чисельного рішення двовимірних газодинамічних течій в областях складної форми з рухомими межами, розробленої Годуновим С.К. і Забродін А.В. з співавторами [19]. Переваги цього методу полягають у тому, що він дозволяє виділити межі областей газу...
