恒星系统与三体运动稳定性的模拟探讨-科展群杰厅.PDF

恆星系統與三體運動穩定性的模擬探討 高中組　第一名 縣　市：高雄市 校　名：高雄中學 　 作　　者：陳威尹 謝俊駒 　 　　　　　 　 　 指導教師：蔡國泰 　 　 　 　 　 　 　 　 　 　 我是澎湖人，自幼在馬公長大，國三時有幸考取雄中，於是寄宿高雄。 從小性好探尋數理之奇與自然之妙，常在閒暇時思考數理問題自娛。因興趣使然，於中學時作了不少小型研究。竟日苦 思、徹夜長考之後，解決難題的喜悅和一絲絲的自豪實乃人生至樂。 願此至樂得與同好共享。 關鍵詞：恆星系統、三體運動、模擬、數值分析、小行星、混沌、雙星、拉格朗日點 　 一、研究動機： 我們在討論天體運動時，當只有一個天體，就行等速度運動；兩個天體時，軌道呈圓錐曲線；但是三個天體的運動分析 並不能算出解析解，而只能以數值分析來算出近似值（早期對於太陽系──簡化的三體運動問題──使用攝動理論（先 考慮恆星重力，再加入行星的影響）做短期模擬（因為長期的攝動共振可以使運動軌跡完全改觀））。 一般而言，想要進行預測星體軌道，人們只要解開相關的微分方程式就可以了。雖然如此，經由牛頓的計算，證實了想 解出只有三個物體的微分方程式都是極度困難的。它的運動定律指出三體運動的模式沒有解析解。要預測三體運動，唯 一的方法就是利用攝動的理論計算近似值。但是現在藉由電腦的的輔助，我們可以用數值分析的方法去做三體運動的模 擬。只要計算夠精確、模擬時間不要太長，都可以得到令人滿意的結果。既然有了良好的工具，就可以探究一些三體運 動的問題了。 說 明 三體運動的最大特徵是「蝴蝶效應」和「重力彈弓」。蝴蝶效應的本意是指「對初始條件的敏感」──微小的初始條件 變化會造成極大影響的效應，也就是有混沌現象（一種介於規律與混亂的現象），像雙擺、彈珠臺等等物理系統；而重 力彈弓是一種三體運動中可能發生的的定性現象──意思是：重力系統中有一物體會受到其他物體的重力影響使它的速 度超過脫離速度而彈射出去，永不回頭。像太陽系中被甩出去的小行星。 拉格朗日點是指當太陽質量遠大於行星質量遠大於小行星質量時，小行星可以穩定存在的點，總共有５個解。本文中只 討論穩定的L4與L5。它是在以太陽、行星連線為一邊的正三角形的頂點上。 二、研究目的： 1.觀察雙星系統中行星的行為。 2.觀察三恆星系統的運動。 3.行星系統的小行星是否能穩定存在。 4.拉格朗日點附近的小行星運動。 　 三、研究設備和器材： ●個人電腦 ●Microsoft Visual Basic 6.0 四、研究過程： 為了實際研究三體運動的細節，我用 Visual Basic 6.0 寫了幾個模擬三體運動的程式。它可以將輸入的多個星體之質量、位 置、和初速度進行運算，得到往後的發展情形，並畫出運行圖。 程式內容省略。 ＊說明1 這三個程式都是用數值積分的方法，但是演算法不同。 程式１是先算出速度再用 的公式計算位置。 程式２是三十八屆物理組第二名科展作品中的方法。他是先用 的公式算出位置，再算出速度。 程式３利用n階階差法來作計算。因為５階計算最精確，而且超過5階就沒有明顯的減少誤差，所以以後提到程式３ 就是指５階階差法。 ＊說明２ 本程式是將複雜的運動切割為很多段的運動，每一個小段有時間間隔Tstep，在本文中以T值表示。因為個人電腦的運算 速度有限，所以其時間間隔T不能取太小，故計算會有一些誤差。為了使精確度和運算時間達到平衡，用以下實驗衡量 誤差和比較三種演算法的優劣。 為了方便設計數值，以下所有實驗的重力常數都為１。 四-(一)、衡量誤差 四-(一)-1、行星圓周運動的軌跡誤差 問　　題：實驗T值的改變對計算誤差的影響。 初始條件：模擬行星系統的變化，其中有一行星對恆星作圓周運動，半徑為１。恆星不動，行星沒有質量。 行星的 實驗方法：以無質量的星體軌道和正圓形比較，計算誤差。 　 　 演算法１： T值設定對誤差最大值做圖（對數座標軸） 結　　論：在圓周運動中誤差最大值和T值成正比。 演算法２：→與演算法１的實驗比起來誤差明顯大很多，而且是發