半径有限薄圆盘引力场的数值计算.DOC

半径有限薄圆盘引力场的的（福州原创物理研究所Cognitech Computer Inc, USA） 内容摘要 本文计算机数值计算方法证明半径有限物质均匀分布引力场的旋转曲线是自然向上弯曲的，并在边界处达到极大。在物质的情况下，半径有限引力场的旋转曲线平坦的或向上弯曲的。处以外的外盘物质会对半径为处以内的内盘物质产生向外的引力。在漩涡星系的演化过程中，由于某些偶然的原因，内盘某个小区域的物质密度比外盘物质密度小时，该区域的物质受到的总引力的方向是指向外盘的。结果会使物质向外盘方向运动，内盘中心区域会形成物质空洞。同时证明这种盘的结构是不稳定，内盘的物质会不断地向外盘运动，使中心的空洞不断扩大，最终导致结构消散。由此就能较合理地解释著名的环状星系形成机制疑难问题，同时还能解释宇宙中环状星系为什么非常稀少的原因。 关键词 牛顿引力，漩涡星系，环状星系，旋转曲线，数值计算方法，天体物理学 一 前 言 由于缺乏的了解。天体力学中常用质量球对称分布代替盘状分布，近似地估计漩涡星系的的球面以外的物质对该球面以内物质的引力不产生影响。但物质盘状分布时，半径为盘面以外的物质对该半径以内物质的引力会产生影响，因此用球状代替盘状分布会产生较大的误差。 本文计算机数值计算方法证明半径有限物质均匀分布引力场的旋转曲线是自然向上弯曲的，并在边界处达到极大。在物质的情况下，半径有限引力场的旋转曲线平坦的或向上弯曲的。盘面以外的物质对该半径以内物质会产生引力。在盘面中心区域的密度小于外部盘的密度时，物质受到的引力是指向盘外，而不是指向盘心的。结果使物质向外盘运动，并在圆盘中心区域形成物质空洞。同时证明这种盘的结构是不稳定，会使内盘物质不断地向外盘运动，使中心的空洞不断扩大，最终导致结构消散。由此就能较合理地解释著名的环状星系形成机制疑难问题，同时还能解释宇宙中环状星系为什么非常稀少的原因。 二 半径有限密度均匀薄圆盘的引力和旋转曲线 2. 1 半径有限密度均匀圆盘的引力场： 径有限密度均匀圆盘的引力场设盘的密度为，半径为，在任意空间点，引力势为【1】： （） 式中，，，，和分别是第三类完全椭圆积分： （2） （） 对上式进行微分，得到在方向单位质量受到的引力： （） （） （） 是Heuman’s Lambda函数。在平面上，，： （） （） 在边界点上，，得： （） （1） 按（2）式有： （） 在边界点上无穷大，这个结果至今没有被注意。本文的目的就是通过计算机数值计算，找到一般性的方法来处理圆盘的引力场问题。 2. 2 半径有限密度均匀薄圆盘引力场的数值计算 假设厚度不计物质均匀分布的薄圆盘半径为，面密度为常数。如图1所示，圆盘内任意点上面积元的质量对盘内点上单位质量的引力为： （1） 从对称性考虑，沿轴方向的引力分量相互抵消，沿轴方向的引力分量为： （1） 整个薄圆盘质量对点单位质量的引力就为： （14）当时有，容易证明此时。因此在圆盘外区间，上式中被积函数是正常的。在的盘内点上被积函数无穷大，（1）式是一致的。 图 1密度均匀半径有限薄圆盘引力设星系盘半径，厚度，质量，体密度。若将它看成薄圆盘，面密度。利用计算机数值方法（1）式计算，考虑速度与力的关系，得到图的旋转曲线。结果表明，薄圆盘的引力和旋转速度都随圆盘半径的增大而增大，且在盘边缘处达到极大。原则上这个峰值应当是无穷大，但由于数值计算采用的步长不可能无穷小，得到的峰不无穷大。 图2 是半径有限密度均匀分布薄圆盘的旋转曲线 是半径有限密度递减薄圆盘的旋转曲线是质量球对称均匀分布的旋转曲线图消除盘边界点上旋转曲线尖锐峰参数选择为了消除边界上无穷大，令盘质量递减在盘边界上密度为零，即： （1）取足够大时是一条于1的直线，图 是时的曲线时足够了。和，半径的圆盘。用数值方法计算得到曲线在边界点上无穷大消失，曲线变得光滑，在边界内仍然向上弯曲 或 （1） 同样的质量球，，密度为，半径则为。按（1）式计算旋转曲线，得到图中的曲线。可知质量和体密度相同时，曲线与有很大的差别。的峰值点比的峰值点明显右移，盘外区域的旋转速度也比球状分布的旋转速度大。密度指数衰减薄圆盘引力场旋转曲线 3. 1半径无限薄盘面密度分布指数衰减时，引力场有解析解。因此天体力学中较精确的方法是，将漩涡星系看成边界无限的薄圆盘，求出盘的引力势，然后通过微分计算求引力。令：