--- title: "从两条轨道推导一小块星环的运动" description: "推导 Kepler 剪切与 Hill 运动，解出自由周转振荡，再看卫星冲量和非弹性碰撞怎样改变局部环片。" canonical_url: "https://sandboxphysics.com/zh/stories/planetary-rings/shearing-orbits/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/planetary-rings/shearing-orbits/" last_modified: "2026-10-28" --- [土星环的波：从轨道到星光](https://sandboxphysics.com/zh/stories/planetary-rings/) AC10 · L2 # 从两条轨道推导一小块星环的运动 Sandbox Physics · 2026-10-28 先从两条圆轨道开始，它们的相对速度给出剪切。再允许粒子偏离圆轨道，把转动系中的科里奥利加速度留下来，就能得到一个可解的局部模型。卫星与碰撞随后再加进去。 ## 开始前需要什么 需要牛顿引力、圆周运动、导数与简谐振动。坐标沿径向向外为正，方位沿顺行方向为正；参考系随卫星的参考圆轨道转动。 ## 先减去参考轨道，再做展开 在参考半径处，圆轨道平衡给出下式中的角频率。把邻近半径写成参考半径加上一个很小的向外位移，再对角频率做一阶展开。在转动系中，将角速度差换成局部方位速度时，使用参考半径这个周长因子；丢掉的修正相对于环片尺度与轨道半径之比是二阶小量。 负号有明确方向：向外位移为正的粒子，在顺行坐标中相对落后。拖动卫星相对采样环条的位置，可以让坐标原点跟随另一条卫星轨道。模型的径向坐标始终从这条卫星圆轨道量起，画面上的偏移不表示把一个天体强行固定在违背剪切的位置。 $$ \Omega(r)=\sqrt{\frac{GM_S}{r^3}},\qquad \frac{d\Omega}{dr}=-\frac{3\Omega}{2r} $$ 先对半径的幂求导，再取局部近似。 $$ \dot y\simeq r_0[\Omega(r_0+x)-\Omega(r_0)]=-\frac32\Omega_0x $$ 一条圆轨道在局部坐标中变成直线剪切轨迹。 对应阅读：Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2 · NASA NSSDCA Saturn Fact Sheet ## 把转动系加速度一起留下 偏离参考位置的粒子还可以有径向和方位速度，科里奥利加速度把两者联系起来。将中心引力与离心加速度围绕参考圆轨道展开，留下一个径向潮汐项。合在一起就是下面的 Hill 方程；右边的力表示卫星等局部作用造成的额外加速度，行星引力已经包含在左边的结构里。 浏览器只保留平面运动，不计算竖直方向、扁率修正、粒子互引力、自转和摩擦。这让机制容易读懂，却不能定量模拟自引力尾迹，也不能标定土星环粘性。场景长度单位是一百米，时间单位是参考角频率的倒数。 $$ \ddot x=2\Omega_0\dot y+3\Omega_0^2x+F_x,\qquad \ddot y=-2\Omega_0\dot x+F_y $$ 径向潮汐与科里奥利项必须一起保留。把速度相关力当作普通位置力处理，会损坏简单蛙跳积分的表现。 $$ r_0=1.30\times10^8\,\mathrm m,\quad GM_S=3.7931208\times10^{16}\,\mathrm{m^3\,s^{-2}},\quad\Omega_0\simeq1.31\times10^{-4}\,\mathrm{s^{-1}} $$ 这组给定参考值对应每个模拟时间单位约 2.11 小时，屏幕播放经过加速。 对应阅读：Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2 · NASA NSSDCA Saturn Fact Sheet ## 先解自由运动，再加入卫星 令额外作用为零。方位方程说明下式中的组合保持常数。用它表示方位速度，再代入径向方程，剩下的就是以导向半径为中心的简谐振动，而不是一个会单独把粒子甩出去的径向加速度。 解里同时有振荡与导向中心的匀速漂移。圆轨道初始条件给出零径向速度和对应的剪切方位速度，振荡幅度恰好为零。卫星的一次冲量通常会改变那个常数组合与径向速度，后续轨迹便一边振荡，一边漂移。 积分器精确推进这段无外力解，在前后各放半次软化卫星冲量。独立的小步长 Runge–Kutta 轨迹核对受迫运动；无卫星剪切和自由周转能量则提供解析检查。这些结果检验的是局部近似，没有检验尚未加入的整环物理。 $$ C=\dot y+2\Omega_0x,\qquad\ddot x+\Omega_0^2x=2\Omega_0C,\qquad x_g=\frac{2C}{\Omega_0} $$ 先对方位方程积分一次，再从径向方程消去方位速度。 $$ x(t)=x_g+(x_0-x_g)\cos(\Omega_0t)+\frac{\dot x_0}{\Omega_0}\sin(\Omega_0t) $$ 在 Kepler 极限中，周转振荡频率等于圆轨道频率。 $$ y(t)=y_0-3Ct-2(x_0-x_g)\sin(\Omega_0t)-\frac{2\dot x_0}{\Omega_0}[1-\cos(\Omega_0t)] $$ 将径向解代入方位速度后积分，别漏掉导向中心的长期漂移。 对应阅读：Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2 ## 把一条轨道和一次碰撞算到底 把粒子放在参考轨道外侧一百米，它的局部方位速度约每秒负 0.0197 米，一小时后向后漂移约七十一米。没有卫星与碰撞时，径向位置保持不变。屏幕上的短暂运动对应很长的物理时间，因此实验明确标注了时间加速。 无摩擦的等质量碰撞只改变沿接触法线的相对速度。动量守恒固定质心运动，恢复系数决定碰撞后的法向相对速度。用碰撞前后的相对动能相减，就得到损失的动能。 举个数值例子：两个各一千克的物体，法向接近速度每秒 0.02 米，恢复系数 0.35。这次碰撞损失约 87.8 微焦耳，切向相对速度保持不变。这里是在计算碰撞定律，没有给每个显示粒子指定一千克的真实质量或真实尺寸。 $$ u_n'=-\epsilon u_n,\qquad\Delta K=\frac{m}{4}(1-\epsilon^2)u_n^2 $$ 两个等质量物体的约化质量为单个物体的一半；恢复系数为一时，法向相对动能不损失。 $$ \Delta K=\frac{1\,\mathrm{kg}}4(1-0.35^2)(0.02\,\mathrm{m\,s^{-1}})^2=8.78\times10^{-5}\,\mathrm J $$ 先说明速度是法向相对速度，再代入能量损失。 对应阅读：Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2 · Salo & Mondino-Llermanos, Photometric modelling of self-gravity wakes and overstable oscillations ## 这套计算能形成哪一种图样 局部模型包含给定卫星、方位周期相接带来的重复经过，以及有限数目的碰撞示踪粒子；径向开放。它可以产生偏折、周转聚集和碰撞重排。高斯密度核随后测量同一批位置。数值碰撞半径与光学权重是独立教学选择，数画出来的圆盘不能测出真实粒径。 整环共振密度波需要集体响应，以及与轨道共振相连的图样速度；自引力尾迹需要粒子互相吸引。额外画几道装饰条纹不能代替这些作用。下一章先把局部光学柱连到探测器，再引入单独的整环色散关系。 对应阅读：Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2 · Goldreich & Tremaine, The excitation and evolution of density waves · Salo & Mondino-Llermanos, Photometric modelling of self-gravity wakes and overstable oscillations ## 练习与解释性答案 1. 将同一圆轨道粒子移到参考轨道内侧一百米，一小时后向哪边漂移？ 展开解释 径向位移变号，剪切速度也变为每秒正 0.0197 米。它在顺行方向领先约七十一米，无额外作用时仍保持原来的半径。 2. 关闭额外作用，但给粒子一次径向速度推动，它会沿径向直线走下去吗？ 展开解释 不会。推动激发周转振荡，径向位置改变时，科里奥利加速度也改变方位速度。径向坐标会围绕导向中心振荡。 3. 一对粒子碰撞后动量守恒，动能却降低，是否违反守恒律？ 展开解释 非弹性碰撞允许这种结果。质心动量守恒，法向相对动能转为未显式计算的内能。可用恢复系数的能量损失式核对；完全弹性时不应损失。 ## 动手检验预测 ### [分开剪切、强迫与耗散](https://sandboxphysics.com/zh/labs/planetary-ring-waves/) 1. 比较无卫星与有卫星预设，再预测偏离参考轨道的粒子会往哪边漂移。 2. 比较弹性与非弹性碰撞，用局部模型的范围解释出现的聚集。 [上一篇：一颗小卫星，怎样在土星环上留下条纹？](https://sandboxphysics.com/zh/stories/planetary-rings/moving-stripes/) [下一篇：从流动的环粒子，算到缺失的恒星光](https://sandboxphysics.com/zh/stories/planetary-rings/columns-and-waves/) ## 来源与计算范围 1. [Rein & Tremaine, Symplectic integrators in the shearing sheet, arXiv v2](https://arxiv.org/pdf/1103.1376v2) 2011 年 4 月 19 日第二版；第 2–3 节的 Hill 方程与精确周转分裂，第 4 节和图 1–4 的自由/受扰轨迹检查。本站实现平面精确自由运动加对称冲量，以及明确软化的给定卫星，没有复现完整多体计算。 2. [NASA NSSDCA Saturn Fact Sheet](https://nssdc.gsfc.nasa.gov/planetary/factsheet/saturnfact.html) 页面更新于 2025 年 3 月 18 日，2026 年 10 月 28 日核读。引力参数与半径提供尺度参照；实验把给定引力参数和十三万千米参考半径作为固定教学输入，没有拟合历表或更新引力场。来源页面只列舍入后的引力参数。 3. [Goldreich & Tremaine, The excitation and evolution of density waves](https://authors.library.caltech.edu/records/ef5pg-y1e07) 1978 年 6 月 15 日发表于 ApJ 222,850–858；核读作者档案全文，第二部分式 4–17，尤其式 9、11–15：薄自引力环片、紧卷色散、共振与卷绕。课程另行推导 Kepler 极限展开，保留压力、扁率和非线性失效范围。 4. [Salo & Mondino-Llermanos, Photometric modelling of self-gravity wakes and overstable oscillations](https://doi.org/10.1051/0004-6361/202452328) 2025 年 3 月 3 日发表于 A&A 695,A37。核读原文索引的 4.2 节/图 5（弱引力不对称与相关检验）、图 4（测光区间）、8.3 节/附录 B（输运贡献）及混合模型范围。通过原文索引读到所选正文，没有复现二十九页数值研究。