重力井戸模拟器 返回
交互式模拟器

重力井戸模拟器

点击放置重力源,在时空橡皮膜类比中体验网格扭曲。通过测试粒子描绘的轨道,直观掌握广义相对论的概念。

点击放置重力源 右键删除 牛顿重力 + Verlet积分
引力势阱
参数
新建质量
重力常数 G
网格密度
操作
预设
显示设置
网格显示
粒子轨迹
力矢量
统计
计算结果
0
重力源
0
粒子数
0
总动能
60
FPS
50
新质量
1.0
G
理论·主要公式

$$F = -\frac{G m_1 m_2}{r^2}$$

万有引力定律(牛):$G = 6.674 \times 10^{-11}$ N·m²/kg²,$r$ 为质点间距离(m)。

$$U = -\frac{G m_1 m_2}{r}$$

重力势能(焦):当 $r \to \infty$ 时为零,接近时变为负值。

$$v_{esc} = \sqrt{\frac{2GM}{r}}$$

逃逸速度(m/s):从半径 $r$ 处逃离引力范围所需的最小速度。

重力井戸模拟器简介

🙋
网格变形表示什么呢?
🎓
这是爱因斯坦"橡皮膜类比"的可视化。简而言之,质量置于时空时,时空本身会弯曲,这是广义相对论的观点。在2D薄膜上呈现。上面的"重力源质量"滑块改变时,可以看到网格沉陷越来越深。这就是"重力井戸"。
🙋
原来如此。那么在屏幕上单击运动的粒子,为什么会旋转或被吸收呢?
🎓
粒子试图在弯曲的网格(时空)上直线行进。但因为舞台本身弯曲,轨道看起来是弯曲的。这就是重力"下落"和"公转"的本质。例如地球绕太阳公转也是这个原理。改变"粒子初速度"滑块,轨道会变成椭圆或抛物线。
🙋
"软化"参数有什么作用呢?
🎓
很好的观察!这在实务(模拟)中非常重要。重力在距离接近零时理论上发散到无穷,计算会崩溃。为了防止这种情况,分母加上小值来稳定计算。试试把右侧的"软化"滑块调得太小,粒子会撞上重力源,计算爆炸。在实际星体碰撞模拟等工程中,这样的调整是必需的。

常见问题

多个重力源时,粒子受到各重力源引力矢量和的作用。轨道不再是简单椭圆,可能产生复杂的混沌运动。例如两个重力源间的"8字"轨道,或被一个吸引后又被抛出的摆动现象都可以观察到。
网格各顶点的高度(变形深度)代表该位置重力势的总和。质量越大或越靠近重力源,沉陷越深。这是广义相对论中时空弯曲的橡皮膜类比表达。
当前版本中,可通过拖动粒子释放时的方向和距离赋予初速。拖动距离越长初速越大。也可在设置面板中直接输入初速度矢量(x、y分量)的数值。
重力源或粒子过多会增加计算负荷。首先右击删除不需要的重力源。也可在设置菜单降低"网格分辨率"(如高→中)减少绘制负荷。仍然卡顿的话,刷新浏览器标签回到初始状态。

现实应用

天体物理学·太空探索:向其他行星发送探测器时的轨道设计(摆动旁路)需要这样的重力场运动计算。模拟器中改变初速创建抛物线轨道的操作是基础概念。

太空望远镜定位:哈勃太空望远镜和詹姆斯·韦布太空望远镜等被放置在地球重力井戸内的特定点(拉格朗日点),实现稳定观测。重力场理解有助于轨道确定。

高精度定位系统(GPS):GPS卫星在地球重力井戸外高速运行。广义相对论效应(时空弯曲导致的时间流速差异)必须补正,否则定位误差数分钟内会达数公里。

数值相对论模拟:黑洞双星合并等强重力场极端现象的计算"数值相对论"领域,将时空弯曲表示在格子上,在其上计算物理量。本模拟器是这类基础概念的初步体验。

常见误解和注意点

首先要明确这个模拟器是2D类比。实际时空弯曲发生在3维空间加时间的4维空间中,这个"橡皮膜"仅是理解辅助。特别是粒子"下落"不是膜的"下"方向,而是弯曲几何本身中的运动这是本质。

其次是"软化"参数的设置。设得过小(如0.1或0.01)时,粒子容易撞上重力源导致计算发散。实务中根据模拟对象的尺度调整。例如星系碰撞模拟中星间距很大,软化参数相对设大来稳定计算。反之精密计算小尺度系时则设小,但相应需要减小时间步长。本模拟器中提高"粒子初速度"时轨道出现震荡,是因为时间步长固定,实务中这是需要注意的。

最后是"Verlet积分"的特性。该方法能量保守性好,但速度不直接计算(从位置差分求得),后来添加速度依赖的力(如空气阻力)时需特别处理。本模拟器仅纯重力计算,可简单使用,应用时要牢记这一点。

使用指南

  1. 用滑块"质量"(sl-mass)设置重力源质量。可在太阳质量相当的10⁸kg~10³⁰kg范围调整。
  2. "重力常数"(sl-g)以牛顿常数G=6.67×10⁻¹¹为基准,可10~1000倍缩放以强调效果。
  3. 在模拟器上任意位置单击放置重力源,周围网格(tog-grid)立即变形,测试粒子(tog-trail)描绘曲线轨道。
  4. "网格分辨率"(sl-grid)改变网格分割数,调整时空扭曲的细节程度。

具体计算例

质量M=1.989×10³⁰kg(太阳质量)、软化参数h=0.15au、初速v=42.1km/s(地球公转速度)时,用Verlet积分、计算时间步Δt=0.001年运行约365步,测试粒子描绘椭圆轨道(长半轴a≈1au)。从重力加速度a=GM/r²,1au处得a≈5.9mm/s²,网格变形以半径r²反比强度可视化。

实务注意点

  1. Verlet积分能量保守性优秀,即使数万步长期轨道计算也最小化轨道衰减。但软化系数h(内部参数)过小会导致尺度不稳定,应设为网格间隔的5~10%。
  2. 多个重力源时,若源间距离小于软化长度,会抑制特异性,物理上接近时防止发散。
  3. 时空网格变形曲率依赖局部质量密度,高密度区有助于理解黑体辐射计算的施瓦茨施尔德半径rs=2GM/c²的概念。与标准求解器(odeint等)结果对比时,相对误差应以1%为目标。