点击放置重力源,在时空橡皮膜类比中体验网格扭曲。通过测试粒子描绘的轨道,直观掌握广义相对论的概念。
$$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$ 处逃离引力范围所需的最小速度。
天体物理学·太空探索:向其他行星发送探测器时的轨道设计(摆动旁路)需要这样的重力场运动计算。模拟器中改变初速创建抛物线轨道的操作是基础概念。
太空望远镜定位:哈勃太空望远镜和詹姆斯·韦布太空望远镜等被放置在地球重力井戸内的特定点(拉格朗日点),实现稳定观测。重力场理解有助于轨道确定。
高精度定位系统(GPS):GPS卫星在地球重力井戸外高速运行。广义相对论效应(时空弯曲导致的时间流速差异)必须补正,否则定位误差数分钟内会达数公里。
数值相对论模拟:黑洞双星合并等强重力场极端现象的计算"数值相对论"领域,将时空弯曲表示在格子上,在其上计算物理量。本模拟器是这类基础概念的初步体验。
首先要明确这个模拟器是2D类比。实际时空弯曲发生在3维空间加时间的4维空间中,这个"橡皮膜"仅是理解辅助。特别是粒子"下落"不是膜的"下"方向,而是弯曲几何本身中的运动这是本质。
其次是"软化"参数的设置。设得过小(如0.1或0.01)时,粒子容易撞上重力源导致计算发散。实务中根据模拟对象的尺度调整。例如星系碰撞模拟中星间距很大,软化参数相对设大来稳定计算。反之精密计算小尺度系时则设小,但相应需要减小时间步长。本模拟器中提高"粒子初速度"时轨道出现震荡,是因为时间步长固定,实务中这是需要注意的。
最后是"Verlet积分"的特性。该方法能量保守性好,但速度不直接计算(从位置差分求得),后来添加速度依赖的力(如空气阻力)时需特别处理。本模拟器仅纯重力计算,可简单使用,应用时要牢记这一点。
质量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²反比强度可视化。