三维引力碰撞系统的数值模拟

本文研究了一个三维引力碰撞质点系统的模型,即“椭球体模型”,其中所有质点是在一个旋转椭球中心体的引力作用下运动,且它们相互之间的碰撞是非弹性的。我们用BCY语言茬TQ-6机上对这个模型的碰撞演化进行数值模拟。计算结果证实了潘加莱的结论(1911)。计算表明该系统在很快变平以后,趋于准平衡态且形成一个有限厚度的盘。盘缓慢地伸展,中心凝聚逐渐形成。这些结论与宇宙中实际的盘状结构相符。