跳到主要内容
Sandbox Physics
从阳光推力读到真实彗尾

L2 · 彗尘动力学

继承的速度,怎样把光压变成尘尾

知道力比,还不能确定尘粒走哪条轨道;还要知道它在哪里、带着什么速度释放。保留这组初始条件,就能看出向内的力为什么也允许粒子逃逸,以及尘尾为什么会记住释放时间。

把释放时刻写清楚

设观察时刻为 t,尘粒年龄为 τ\tau,先求彗核在更早的 t−τt-\tau 时刻的位置与速度。尘粒从那里出发,速度等于彗核速度加上脱离气体拖曳后的释放速度;随后在减弱引力下走完整个年龄。彗核则在正常太阳引力下走到同一个观察时刻。

只有两条轨道到达同一时刻,才能相减求尘粒相对彗核的位置。把粒子放在今天的彗核,再向前走一个释放年龄,会画错几何。把纯太阳模型的粒子直接减去含完整摄动力的 JPL 彗核,也可能在力比为零时造出假尾巴。真实图像教程会专门处理这种不一致。

rd(t−τ)=rc(t−τ),vd(t−τ)=vc(t−τ)+vej,Δr(t)=rd(t)−rc(t)\boldsymbol r_d(t-\tau)=\boldsymbol r_c(t-\tau),\quad \boldsymbol v_d(t-\tau)=\boldsymbol v_c(t-\tau)+\boldsymbol v_{\rm ej},\quad \Delta\boldsymbol r(t)=\boldsymbol r_d(t)-\boldsymbol r_c(t)
d、c 分别指尘埃与彗核。所有三维向量使用同一太阳原点、坐标框架与时间单位。把力比和释放速度都设为零,必须恢复彗核轨道;这是很有力的检验。

对应阅读:Finson & Probstein · JPL Horizons API and observer/vector conventions

从能量求逃逸门槛

先令释放速度为零。母椭圆半长轴为 aca_c 时,速度平方为 μ(2/r−1/ac)\mu(2/r-1/a_c),其中 μ=GM⊙\mu=GM_\odot。尘粒继承这个速度,势能却要用减小后的 μ(1−β)\mu(1-\beta)。把动能减去引力势的绝对值,再整理。

单位质量能量为负,对应受束缚轨道;等于零是抛物极限,为正便是双曲逃逸。门槛取决于释放点在母椭圆的哪里。同样大小的粒子,在远处释放可能受束缚,在近日点释放却可能逃逸。叠加释放速度后,动能会变,门槛还会移动。

εd=vc22−μ(1−β)r=μ(βr−12ac),βescape=rrel2ac\varepsilon_d=\frac{v_c^2}{2}-\frac{\mu(1-\beta)}{r}=\mu\left(\frac{\beta}{r}-\frac{1}{2a_c}\right),\qquad \beta_{\rm escape}=\frac{r_{\rm rel}}{2a_c}
此式要求母轨道为受束缚开普勒轨道、释放速度为零。它不同于合力向外所需的力比大于一。教学椭圆偏心率为零点六五,近日点门槛为零点一七五,远日点为零点八二五。

对应阅读:Burns · Finson & Probstein

完整算一次圆轨道释放

从一天文单位圆轨道释放力比为零点一的尘粒,不另加释放速度。母天体速度平方等于太阳引力参数除以一天文单位;尘粒保留这个速度,也保留释放瞬间的单位质量角动量。先从能量求新半长轴,再从角动量求偏心率。

新半长轴为一点一二五天文单位,偏心率为九分之一。释放点恰是新轨道的近日点,仍在一天文单位;远日点则到一点二五。力比为二分之一时,原有圆轨道速度恰好等于减弱引力下的逃逸速度;再大便逃逸,但力比尚未到一以前,加速度仍然向内。

ad=1−β1−2βrrel,ed=β1−β(0≤β<1/2)a_d=\frac{1-\beta}{1-2\beta}r_{\rm rel},\qquad e_d=\frac{\beta}{1-\beta}\quad(0\le\beta<1/2)
这组表达式只适用于圆形母轨道与零释放速度。它来自能量和角动量关系。偏心彗星应先计算实际释放状态,再从能量与角动量求轨道。

对应阅读:Burns · Finson & Probstein

从一颗粒子,扩展到两类轨道族

等时线固定释放时刻、改变力比;等粒径线固定力比、改变释放时刻。在零释放速度、材料性质固定的近似下,两组曲线展示了简单尘尾图像的主要自由度。投影到天空后,本来分开的空间曲线可能交叉,也可能几乎重合。

交互场景另加了气体加速结束后的释放速度,按粒径和释放日心距离的平方根反比缩放,再用球面四个方向作稀疏三维采样。每个点的视觉权重相同,因此画面点密度不是表面亮度预测。若要预测亮度,还需尘埃生产历史、散射相函数和粒径分布。

vej=v0s0s(rrelAU)−1/2v_{\rm ej}=v_0\sqrt{\frac{s_0}{s}}\left(\frac{r_{\rm rel}}{\mathrm{AU}}\right)^{-1/2}
这是教学释放定律,不是 ISON 气体拟合。青色参考粒子不加释放速度,因此拖动粒径能单独观察辐射压响应;下方彩色族线也不加释放速度。

对应阅读:Finson & Probstein · Li et al.

留下可复算的结果,再检查画法

轨道位置为六十五度时,教学彗核距太阳约一点二七天文单位。释放历史覆盖一百二十日、中心半径三微米、速度尺度每秒一百二十米时,采样粒子最远距彗核约零点七一一天文单位。下表记录刚进入页面、时钟还未推进时的代表结果,数值由前述离散粒子族计算后舍入。

若要复算这组离散例题,母轨道半长轴取二点八天文单位,偏心率零点六五,倾角三十五度、升交点经度三十五度、近日点幅角二十五度。年龄下标从一至二十八,粒径下标从零至六;四个释放方向按下式在教学惯性轴中取点。释放速度先换成每天多少天文单位,再与轨道速度相加。

圆锥轨道传播器用通用变量穿过受束缚和逃逸两类状态。独立龙格—库塔积分比较圆轨道、近抛物线与双曲线,另检验能量、角动量和反向传播。取景另算几何:太阳与所有粒子都进入视野,两条图像轴保持相同距离比例。

τj=120j28 d,sk=3 22k/6−1 μm\tau_j=\frac{120j}{28}\,\mathrm d,\qquad s_k=3\,2^{2k/6-1}\,\mu\mathrm m
此默认例题取二十八档年龄、七档半径。换参数时,半径限制在一至一百微米之间。
zℓ=1−2(ℓ+1/2)4,ϕℓ=ℓπ(3−5),u^ℓ=(1−zℓ2cos⁡ϕℓ,1−zℓ2sin⁡ϕℓ,zℓ)z_\ell=1-\frac{2(\ell+1/2)}4,\quad\phi_\ell=\ell\pi(3-\sqrt5),\quad\hat{\boldsymbol u}_\ell=(\sqrt{1-z_\ell^2}\cos\phi_\ell,\sqrt{1-z_\ell^2}\sin\phi_\ell,z_\ell)
方向下标从零至三,单位向量乘上释放速度定律,便组成默认七百八十四个粒子。速度为零时只需一个方向,重复方向会给出同一条轨道。
情形中心力比释放历史
默认0.19128 档年龄、7 档半径、4 个方向
细尘0.574同样年龄;小粒子分离更强
粗尘0.00574180 日;速度尺度每秒 5 米
关闭光力与释放速度0全部粒子与彗核重合

对应阅读:Finson & Probstein · Burns

亲手观察

比较同一年龄与同一粒径

  1. 暂停“一次喷发,多种粒径”,观察窄曲线;换成“同一粒径,不同年龄”,比较弯曲。
  2. 从零开始增大释放速度尺度,转向接近轨道平面的方向,看厚度怎样出现。
  3. 把彗核移到近日点附近,再移到远处。同一个力比,不保证相同的轨道能量。

练习与答案

力比零点二、零释放速度的粒子,在教学轨道近日点释放会逃逸吗?远日点呢?

门槛分别为零点一七五与零点八二五。零点二超过近日点门槛,能量为正;低于远日点门槛,仍受束缚。两处的力比都小于一,所以中心力都向内。

零光力、零释放速度的粒子,相对 JPL 彗核却画出一条尾巴,先查哪里?

先查粒子和参考母轨道是否用了相同的力、原点与时刻。纯太阳粒子与完整摄动 JPL 彗核未必重合。先让粒子与母体在同一太阳模型、零力比下重合,再把相对偏移加到共同参考上。另一种做法是给两者一致加入缺失的力。

一条投影亮脊贴合某条等时线,就能唯一确定喷发日期吗?

不能只凭亮脊判断。其他粒径、年龄、释放速度和视线方向可能给出相近投影。需比较多期图像、多个方向特征和亮度模型,才能进一步把曲线解释为尘埃生产事件。

来源与计算范围

  1. Finson & Probstein, A Theory of Dust Comets I (1968)

    ApJ 154,327—352,DOI 10.1086/149761。第 I—III 节,特别是第 329、334—335 页与图 4:脱离气体拖曳后的尘粒、减弱的太阳引力、等时线与等粒径线。教学场景另加明确写出的释放速度假设。

  2. Burns, Lamy & Soter, Radiation Forces on Small Particles (1979)

    Icarus 40,1—48,DOI 10.1016/0019-1035(79)90050-2。第 III—IV 节推导一阶辐射力;第 VI 节式 (19) 与图 7 讨论压力效率;第 VII—VIII 节、式 (47)—(50) 讨论减弱引力的轨道与轨道平均 PR 演化。核读作者托管的 48 页全文。

  3. Li et al., HST observations of Comet ISON, arXiv:1311.0826v1

    第 2 节给出 2013 年 4 月 HST 观测序列与原处理;向日内喷流与延伸尘尾促使教程在粗略方向诊断中排除中心两角秒。教程使用档案 2026 年处理的单帧,与论文合成图不同。

  4. JPL Horizons API and observer/vector conventions

    API 响应版本 1.2,API 文档版本 1.3。固定彗星记录 90004196、解 JPL#57(2025-Nov-21_13:42:48)、DE441。对应 TIME_TYPE、TLIST、ICRF/FRAME 向量约定和观测量 1、19、20、21、24、27、30;保存的请求与响应保留时标和几何选择。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策