太陽系n體模擬的流形改正方法及其應(yīng)用.pdf

1、可靠的數(shù)值方法是非線性動力學(xué)研究的基礎(chǔ)。由于傳統(tǒng)數(shù)值算法引入人工耗散，而辛算法應(yīng)用有限，因此尋找新的可靠的數(shù)值算法成為當(dāng)前非線性動力學(xué)研究的重要課題。 本文基于Nacozy提出的流形改正思想以及積分不變關(guān)系，設(shè)計了三種流形改正方案用于保持太陽系n體問題中每個天體的運動緩變量，提高與之相關(guān)軌道根數(shù)的精度。 首先，分別構(gòu)造了速度因子方法和坐標(biāo)因子方法，用于保持系統(tǒng)中每一體的開普勒能量。也就是在哈密頓正則方程的數(shù)值解中，對速度

2、分量或者坐標(biāo)分量乘以一改正因子，其值在嚴格的開普勒能量關(guān)系中求出，下一步的計算都以改正后的值為初始值。通過在理想二體和日木土三體問題中數(shù)值驗證，無論是速度因子方法還是坐標(biāo)因子方法，都取得了與前人一致的結(jié)果，極大的提高了軌道半長徑的精度。一般情況下，速度因子方法的求解要比坐標(biāo)因子方法簡單，故推薦在實際應(yīng)用中采用速度因子方法。 其次，考慮同時對開普勒能量和拉普拉斯積分的改正，用于提高與之對應(yīng)的軌道半長徑和偏心率的精度。將開普勒能量和

3、拉普拉斯積分對速度分量的偏導(dǎo)數(shù)組成一2×3矩陣，根據(jù)Nacozy流形改正原理，求出速度改正因子。數(shù)值實驗表明，新方案取得了預(yù)期的效果，并與Fukushima的雙因子法相比，其更適合于高偏心率情況。 再次，對太陽系n體問題中與每一體6個軌道根數(shù)直接或間接相關(guān)的開普勒能量、3個角動量分量、拉普拉斯積分z分量同時改正。與第二個工作相比，新方法考慮的積分量更多，而且對坐標(biāo)和速度分量都要求導(dǎo)。通過數(shù)值實驗驗證，此方案與Fukushima的

4、線性轉(zhuǎn)換方法在對軌道根數(shù)精度的改進上效果一致，所有的軌道根數(shù)的精度均保持在雙精度的量級。 最后，將速度因子改正方法應(yīng)用到共形耦合標(biāo)量場閉宇宙模型當(dāng)中。一方面保持了能量面，另一方面，也為該模型非線性動力學(xué)性質(zhì)演化提供了高精度的數(shù)值解，避免了數(shù)值誤差帶來的影響，使得結(jié)果更為真實可靠。通過混沌識別方法Poincare截面分析，表明宇宙學(xué)常數(shù)對宇宙非線性的演化有一定的影響。當(dāng)宇宙學(xué)常數(shù)和自相互作用系數(shù)都小于1的時候，它們對系統(tǒng)非線性沒有