Кіріспе

Күн жүйесінің сандық моделі – уақытқа байланысты планеталардың шамамен орналасуын анықтайтын математикалық теңдеулер жиынтығы. Мұндай модельді жасау әрекеттері аспан механикасының дамуына негіз болды. Осы модельдеу нәтижелерін өткен өлшемдермен салыстырып, дұрыстығын тексеруге және болашақ орналасуларын болжауға болады. Оның басты қолданысы – альманахтарды дайындау.

Бұрынғы еңбектер

Симуляциялар декартық немесе сфералық координаттарда жүргізілуі мүмкін. Алғашқысы оңай, бірақ есептеулерге өте көп күш жұмсалады және тек электрондық компьютерде ғана тиімді. Сондықтан бұрын тек соңғысы қолданылған. Шындығында, соңғысы да есептеулерге аса көп күш жұмсауды қажет етпеді, бірақ бастапқыда қарапайым шамалаудан бастап, қажетті дәлдікке жету үшін қажет болғанша бұзылыстарды қосуға болатын. Бұл Күн жүйесінің математикалық симуляциясы, ол N дене проблемасының бір түрі болып табылады. N символы денелердің санын көрсетеді, ол Күнді, 8 планетаны, ондаған серіктерді және сансыз кішкентай планеталарды, кометаларды және т.б. есепке алғанда өте үлкен болуы мүмкін. Алайда, Күннің кез келген басқа денеге әсері өте зор, ал қалған денелердің бір-біріне әсері өте шағын болғандықтан, мәселені аналитикалық шешімі бар 2 дене проблемасына дейін тоғытуға болады. Әр планета үшін нәтиже – орбита, уақыт функциясы ретінде оның орналасуының қарапайым сипаттамасы. Бұл шешім табылғаннан кейін, серіктер мен планеталардың бір-біріне тигізетін әсері шағын түзетулер ретінде қосылады. Бұл түзетулер толыққанды планеталық орбитамен салыстырғанда өте кішкентай. Кейбір түзетулер әлі де бірнеше градусқа жете алады, ал өлшеулер 1 секундтан артық дәлдікпен жүргізілуі мүмкін. Бұл әдіс қазіргі уақытта симуляциялар үшін қолданылмаса да, шамамен эфемерида табу үшін әлі де пайдалы, өйткені салыстырмалы түрде қарапайым негізгі шешімді алып, ең ірі бұзылыстардың бірнешеуін қосып, планетаның қалаған орналасуын көп еңбекті жұмсамай табуға болады. Бірақ бұл әдістің кемшілігі – бұзылыстар теориясы өте күрделі математика саласы.

Қазіргі заманғы әдіс

Қазіргі заманғы әдіс 3 өлшемді кеңістікте сандық интеграциядан тұрады. Біріншіден, қатысатын әрбір дененің орны (x, y, z) және жылдамдығы (vx, vy, vz) үшін жоғары дәлдік мәнімен басталады. Әрбір дененің массасы белгілі болғанда, Ньютонның гравитация заңы бойынша оның үдеуі (ax, ay, az) есептелуі мүмкін. Әрбір дене басқа денелерді тартады, ал жалпы үдеу – осы тартылыстардың қосындысы. Келесі қадам – Δt кішкентай уақыт аралығын таңдап, Ньютонның екінші қозғалыс заңын қолдану. Үдеуді Δt-ға көбейту жылдамдыққа түзету енгізеді. Жылдамдықты Δt-ға көбейту орнына түзету енгізеді. Бұл процедура басқа барлық денелер үшін қайталанады. Нәтижесінде барлық денелердің орны мен жылдамдығы үшін жаңа мән шығады. Содан кейін осы жаңа мәндерді пайдаланып келесі Δt уақыт аралығы үшін есептеудің барлық кезеңдері қайта басталады. Бұл процедураны жеткілікті түрде қайталаса, уақыт өте келе барлық денелердің орналасуын сипаттауға болады. Бұл әдістің артықшылығы – компьютер үшін бұл өте оңай жұмыс, және ол барлық денелер үшін бір уақытта жоғары дәлдік нәтижелер береді, бұл бұзылыстарды анықтаудың күрделі және қиын процедураларын қажет етпейді. Кемшілігі – бастапқыда өте дәл мәндермен жұмыс істеу қажет, әйтпесе нәтижелер уақыт өте келе шындықтан алшақтайды; x, y, z координаттарын аламыз, олар көбінесе пайдалану алдында практикалық эклиптикалық немесе экваторлық координаттарға түрлендіріледі; және бұл «бәрі немесе ештеңе» тәсілі. Егер сіз бір планетаның белгілі бір уақыттағы орнын білгіңіз келсе, онда басқа планеталардың және аралық уақыт кезеңдерінің есептеулері де қажет.

Асқынулар

Егер ең жоғары дәлдік қажет болса, есептеулер әлдеқайда күрделі болады. Кометалардың жағдайында, радиациялық қысым және газдың қарсыласуы сияқты гравитациялық емес күштерді ескеру қажет. Меркурий және басқа планеталардың жағдайында, ұзақ мерзімді есептеулер үшін релятивистік эффектілерді назардан тыс қалдыруға болмайды. Сондай-ақ, жалпы энергия енді тұрақты емес (өйткені төрт векторлық энергиясы мен сызықтық импульс өзара байланысты). Жарық жылдамдығының шекті болуы, классикалық және релятивистік эффектілердің екеуінің де жарық уақытының әсерін ескеруді маңызды етеді. Планеталарды енді бөлшектер ретінде қарастыруға болмайды, олардың пішіні мен тығыздығы да ескерілуі керек. Мысалы, Жердің экваторлық сығылуы прецессияны тудырады, ол барлық планеталардың ұзақ мерзімді қозғалысына әсер ететін осьтік бұрылыстың өзгеруіне алып келеді. Күн жүйесінің тұрақсыздығына байланысты, ондаған миллион жылдан астам мерзімді модельдер жасау мүмкін емес.