设万维读者为首页 万维读者网 -- 全球华人的精神家园 广告服务 联系我们 关于万维
 
首  页 新  闻 视  频 博  客 论  坛 分类广告 购  物
搜索>> 发表日志 控制面板 个人相册 给我留言
帮助 退出
小鹰的博客  
小鹰的博客  
https://blog.creaders.net/u/24357/ > 复制 > 收藏本页
网络日志正文
小鹰:三体问题 (附多幅照片及视频) 2024-05-12 16:15:19

由于相关的小说和影视作品上市,最近“三体”这个词儿在中国火了。的确,物理学的三体问题,因其解的混沌性质,容易令人遐想,或墬入迷思。

2006年,出于好奇,我曾研读了大卫·艾奇逊的《从微积分到混沌》[1]一书,其中在二维平面空间中,讨论了三体的引力运动问题。

现将笔记整理如下。

我们知道,如果只有二体存在,只要给出初始条件,根据牛顿定律,即可以精确解出二体的运动轨迹和速度。

然而,在三体情形下,如图一所示,每个粒子受到其它二粒子的引力作用,根据牛顿第二定律,固然也可以写出各粒子在X和Y方向上的六个二阶运动方程,但这些方程式彼此耦合,没有解析解。

图一

3BD_1.jpg

通常的求解办法是:将其降阶为十二个一阶微分方程,经去量纲化(non-dimensionalization)後,用可变步长的数字叠代法(the Runge-Kutta method)计算,近似求出每个粒子瞬时的坐标和速度。

图二a显示,如果没有第三粒子存在,当粒子m1和m2的质量相同,等于0.5,初始位置分别在X轴的两侧等距之(0,-0.3)和(0,0.3)点处,初速度大小相等,一个向南,一个向北,即为(0,-0.3)和(0,0.3)时,两个粒子将各自绕系统重心做稳定的椭圆运动。宇宙中观察到的联星或复星(binary star)现象,即是如此。

但当有第三粒子介入时,每个粒子的轨迹会变得极为复杂。

例如,具有同样质量的m3,从(-0.1,0.75)的位置,以(0,-0.3)的初速度向南进入,如图二b、c、d所示,粒子1与粒子3很快接近,一起向西北方向行,粒子2则被排斥到东南方。不久,粒子2又转回来,与粒子3缠绕,把粒子1踢了出去。

图二

3BD_2.jpg

这种相互纠缠,不断变换“伙伴”的关系,在下面的图三中可以看得更清楚,此时系统其它条件与图二相同,只是m3的初速改为(0,-0.2)。

图三

3BD_3.jpg

该图显示,粒子3与彼此纠缠的粒子1和2最後可能再次相遇,那之後的结果又会是怎样?

笔者按照Acheson书中提供的程序,进一步演算,并将全过程拍成影片,後经友人建议和帮助,又改用录屏技术,获得更清晰的视频。有兴趣的读者可以点击这里http://www.azcolabs.com/3BD.mp4观看MP4模式的视频链接:三体动态图。如有读者无法观看MP4模式,可试点击AVI模式的原摄影片视频链接 http://www.azcolabs.com/3BD.avi。该短片长约三分钟,但可按播放软件中的快进键或选择播放速度,加快播放。

影片结尾时显示,粒子3(黑)与粒子1/2(红与绿)相会之後,再次交换位置,最终,粒子2/3(绿与黑)结伴,类似一对双星(double star),一路向西,粒子1(红)则形同一颗脱离联星的“流浪行星”(a rogue planet),在太空深渊中独自投东,踽踽而行。二者背道而驰,渐行渐远,达到了三体衰变(decay)或解构崩溃的结局。

有些读者可能会问,在我们生活的多体大千世界,怎麽没有见过如此诡异的行径?太阳系的九大行星不是秩序井然地运行多年,安然无恙麽?而且其轨迹可以预测,人们准确知道月蚀、日蚀会在什麽时候发生,也可推算出哈雷彗星回归的日子。

我的理解是:我们这里说的三体混沌现象会发生,关键在于三体间的相互作用量级相当,彼此影响,形成强耦合状态,以致其运动完全无法预测。

但是,如果一个体系,只有一个物体占绝对优势,如太阳系中太阳的质量占99.86%,其余众多小物体,如诸行星及其卫星,则各自围绕它做稳定的二体运动;至于行星彼此间的吸引,由于质量小且相距较远,属于高阶小项微扰,可能会引起行星的些许进动,或者可以忽略不计。

此时,若有更小的物体,以一定的初始条件,闯入二体之间,结果会怎么样?

以下四张图片,对此做了说明。

初始条件与图二相同,第三体质量m 0。

第一种情形,如图四显示,第三体从远处(0,10)的位置,以初速(0,-0.3)出发,受二体的吸引加速向下。如果加速时间够长,使第三体获得足够的动能,或当切入时另二体恰在远端附近,那它就有可能直接穿越二体体系。

图四

3bd_a.jpg

第二种情形,如图五显示,第三体从近处(0,1)的位置,同样以初速(0,-0.3)出发,因受二体吸引加速时间较短,第三体速度不够高,且刚好与第二体接近,结果被第二体俘获,成为绕其旋转的卫星。然而,经过相当的时间之後,第三体又解脱出来,向西南方向飞去。

图五

3bd_b.jpg

第三种情形,如图六显示,第三体从稍远的(0,3)位置,同样以初速(0,-0.3)出发,因受二体吸引加速时间略长,第三体速度居中,它既不被俘获,又无力直接穿越,结果在二体的势井中,混沌漫游一番之後,终于得以向东北约60度方向逃逸。

图六

3bd_c.jpg

第四种情形,如图七显示,第三体从偏离Y轴的(-0.1,0.75)位置,同样以初速(0,-0.3)出发,结果与上相似。只是这一次,它没能幸运逃逸,而是直接撞上了第二体。图中显示二粒子相撞的位置在屏幕的(318,120)点处。

图七

3bd_d.jpg

从月球上的无数陨石坑和地球上夜间常见到的流星,可知这种相撞情形并不罕见。

太阳系存在已有几十亿年,该撞的或能逃的天体大概多已被淘汰了,现在剩下的,可以说,是一个基本稳定的结构。但多体耦合导致的混沌法则依然存在,从大尺度来看,一切仍旧可能改变,包括月球每年正悄悄地远离地球约一英寸多,无数小行星和碎石块仍在太阳系里杂乱飞舞,随时威胁着地球上的生命存在。

而放眼银河系,那里大约有4万亿个不隶属于任何恒星的“流浪行星”(the rogue planets)在星际空间中游荡,其大小如地球到木星不等,有些还带有自己的卫星环绕。

如前文所述,这些流浪行星很可能是从联星或复星天体系统(binary stellar systems)中因混沌被排斥出来的。

以上只是些数字模拟,试图形象地给出有关三体混沌运动中的一些基本概念,以及可能发生的情形。

浩瀚的宇宙是和谐秩序与紊乱混沌的统一体,然而,无论如何,质量相近的三体天体极不稳定,应不适合生命起源和进化的漫长过程,遑论产生造就高级智慧生物。 

[1],David Acheson,From Calculus to Chaos,Oxford,1997。

To read English version please click: http://www.azcolabs.com/xy_3BD_Eng.html


写于2024年4月28日。


浏览(3216) (1) 评论(0)
发表评论
我的名片
小鹰
注册日期: 2021-09-07
访问总量: 186,994 次
点击查看我的个人资料
Calendar
最新发布
· 小鹰:三体问题 (附多幅照片及
· 小鹰:从北大的一幅挽联谈起 (
· 小鹰:熵与状态
· 小鹰:援“穆”求愚──五答穆白硕(
· 小鹰:辩论刘少奇
· 小鹰:四答穆白硕(雨村)
· 小鹰:三答穆白硕(雨村)
分类目录
【文革与溯源】
· 小鹰:援“穆”求愚──五答穆白硕(
· 小鹰:辩论刘少奇
· 小鹰:四答穆白硕(雨村)
· 小鹰:三答穆白硕(雨村)
· 小鹰:看看马克思还讲了些什么?
· 小鹰:再答穆白硕(雨村)
· 小鹰:答穆白硕(雨村)
· 小鹰:四、结 语
· 小鹰:三、“文革”反动了什么?
· 小鹰:二、“文革”是什么?
【随笔与杂评】
· 小鹰:三体问题 (附多幅照片及
· 小鹰:从北大的一幅挽联谈起 (
【鲁迅与文艺】
· 小鹰:叫不醒的“书呆子”里,还有
· 小鹰转发:《圣诞节》注解
· 小鹰(整理):对《转了“基因”的
· 小鹰:转了“基因”的“钱式鲁迅”
· 小鹰:对谈钱理群先生的“鲁迅观”
· 小鹰:再谈解读《影的告别》及其
· 小鹰:如何解读《影的告别》?──
· 小鹰:沉默的鲁迅,还是呐喊的鲁
· 小鹰:评戴克刚的《那个…鲁迅》
· 小鹰:从钱理群的 “XXX鲁迅”谈起
【追忆与思考】
· 小鹰:从北大的一幅挽联谈起 (
· 小鹰:纪念丹尼尔•艾尔斯
· 小鹰:有两拨人
· 小鹰:从悼念唐伟所想到的
· 小鹰:闭卷判分 名人归零
· 小鹰:争论的价值──犹太文化的启
· 小鹰:文革中我参加过的“相对论
存档目录
2024-05-12 - 2024-05-12
2023-11-06 - 2023-11-06
2023-10-12 - 2023-10-20
2023-09-04 - 2023-09-27
2023-08-05 - 2023-08-22
2023-07-05 - 2023-07-29
2023-06-17 - 2023-06-29
2023-02-07 - 2023-02-07
2022-03-15 - 2022-03-15
2022-02-12 - 2022-02-26
2022-01-02 - 2022-01-20
2021-12-01 - 2021-12-23
2021-11-05 - 2021-11-19
2021-10-02 - 2021-10-28
2021-09-09 - 2021-09-27
 
关于本站 | 广告服务 | 联系我们 | 招聘信息 | 网站导航 | 隐私保护
Copyright (C) 1998-2024. Creaders.NET. All Rights Reserved.