多体问题的计算机模拟.pdfVIP

多体问题的计算机模拟 蒋绍周 云天梁 2001012474 2001012503 （清华大学物理系物11 班） 摘要：借助数学软件 MATLAB ，用数值解微分方程，模拟三星系统在各种初条件下的 演化，并进行分析。 关键词：三体问题，数学软件MATLAB 宇宙中单恒星的存在相对来说是少数的，大量存在的是双星系统，还有一定 的三星系统。双星的运动问题比较简单，早就有完美的解析分析。而由于三体问 题没有完美的解析解，受初始条件的影响很大，所以有必要用计算机进行数值模 拟，给出其演化规律。下面就给出几种典型的演化过程。 在模拟时，不考虑恒星自身的演化，只考虑其它星的引力影响。 当速度与三星组成的平面共面时，下图是一种演化的可能情况 图1 三星中有两颗星在经过一段时间后发生碰撞合成一颗星，上图显示了碰撞的 瞬间。碰撞之前，各恒星做类似螺线的运动，大质量的恒星的回旋半径较大，周 期较长；小质量恒星的回旋半径较小，周期较短。分别对于上图中绿色和蓝色的 轨迹。发生碰撞后，剩下的两恒星的运动如下图。 图2 实际上这就成了两体问题，由于初速度的选取使质心速度不为零，所以双星 还有平动。 如果初始条件选得好，使得三星的运动过程中没有相互碰撞的现象，这得到 的是另一景象，如下图 图3 开始的时候红色和绿色轨迹的星距离较近，互相缠绕着转动，蓝色轨迹的星 在距离较远在外围转动。经过一段时间后如下图， 图4 三颗星的距离又靠近，然后变为另两颗星距离较近相互转动，如下图 图5 这时，是绿色和蓝色轨迹的星距离较近，缠绕着运动，红色的星距离较远。 过了一段时间以后，变为红色与蓝色的轨迹的星距离较近，如下图 图6 如此重复，在每一阶段，都有两颗星的距离较近，另外一颗星距离较远。当 那颗距离较远的星转回去时，对另两颗星产生影响，最后缠住一颗星，使另一颗 星远离。结果是每颗星都有远离的机会，每次回来时，换了另一颗星远离，如此 循环。 以上模拟的是速度与三星的位置共面时的情况，如果速度与三星位置不共 面，下图是一种可能的演化方式 图7 三颗星开始的运动由于受另外两颗星的作用都很大，规律不是很明显。但是 经过一段时间以后，有两颗星碰到了一起( 图中为红蓝两色轨迹的星碰到了一 起) ， 图8 还有一种情况是由于初速度太大，各星相互散开了，这种三星系统不稳定， 不再作讨论了。还有其它各种演化方式，由初条件不同而不同，如两颗星相互围 着转，再于第三颗星围着转等等，就不一一列举了，具体的演示可参看我们的程 序。 在模拟的过程中还发现，大多数情况下三星系统都有两颗星碰到一起，最后 形成双星系统，形成稳定的三星系统的情况较少，这可能是宇宙中双星系统较多 的原因之一。