---
title: "读懂 PR 拖曳与热反冲各自的时钟和证据"
description: "对照固定版辐射力论文与雅可夫斯基综述，读方程、定义热模型，再从图中分清轨道拟合与力的解释。"
canonical_url: "https://sandboxphysics.com/zh/stories/comet-dust/reading-slow-radiation-forces/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/comet-dust/reading-slow-radiation-forces/"
last_modified: "2026-10-28"
---

[从阳光推力读到真实彗尾](https://sandboxphysics.com/zh/stories/comet-dust/)

L4 · 彗尘动力学

# 读懂 PR 拖曳与热反冲各自的时钟和证据

Sandbox Physics · 2026-10-28

把 Burns、Lamy、Soter 的 1979 年论文与 Vokrouhlický 等 arXiv:1502.01249v1（2015 年 2 月 4 日）放在旁边。先跟着前者读能量损失，再读后者的热变量和图 2—3，判断轨道漂移究竟在检验什么力学模型。这里导读固定原文，没有重跑小行星天体测量拟合，也不汇总此后新探测。

## 开始前需要什么

上一章给出 PR 能量与圆轨道时间推导。热惯量描述材料在周期受热时抵抗温度变化的性质；倾角是自转角动量与轨道角动量的夹角。

## Burns 第 VIII 节，先认出平均的是哪条轨道

读式 (47)—(50) 时，把前面的力比定义和一阶辐射力一起翻开。标出受束缚轨道、半长轴、偏心率，以及对一周公转取平均这一步。瞬时加速度与轨道平均漂移的单位、用途都不同。先恢复上一章的圆轨道极限，再使用偏心轨道结果。

可以在旁边重算十微米、一天文单位例题：力比约零点零五七四，圆轨道形式时间约六千九百七十四年。这只时钟属于定义清楚的受束缚粒子，不能直接贴在任意逃逸的 ISON 尘粒上；为方便观看而改光速，也会改变物理模型。

对应阅读：Burns · JPL Solar System Dynamics astronomical parameters

## 先有温度模型，才有热推力

对密度、比热和热导率固定的均匀固体，用热扩散方程描述温度。表面外法向记为 n，吸收的阳光分成热辐射和流入内部的热量。求解前还须选定自转状态与周期内部边界条件。下面给出最小示例模型，足以说明热推力还需要哪些输入。

C 是单位质量比热，k 是热导率，A 是邦德反照率，发射率无量纲，S 是当地太阳能流。遮阴、粗糙度、变热物性或升华出气都需要补充项。得到温度图后，再假设朗伯热发射，对全表面辐射动量积分，才得到反冲力。

$$
\rho C\frac{\partial T}{\partial t}=k\nabla^2T,\qquad k\frac{\partial T}{\partial n}=(1-A)S\max(\hat{\boldsymbol n}\cdot\hat{\boldsymbol s},0)-\epsilon\sigma T^4
$$

太阳方向从表面指向光源，边界各项单位都是每平方米多少瓦。采用此符号，表面较热时，热量导入内部。

$$
\boldsymbol F_{\rm th}=-\frac{2}{3c}\int_{\rm surface}\epsilon\sigma T^4\hat{\boldsymbol n}\,dA
$$

三分之二来自对向外半球的朗伯光子发射积分。球面温度均匀时，各处法向相消，静止系合力为零；另行推导的一阶 PR 参考系修正仍然存在。

对应阅读：Vokrouhlický et al. · Burns

## 综述第 2.1 节，先读符号再读大小

对照式 (1)—(3) 里的自转倾角 $\gamma$ 与热响应 W。沿用文中经典假设，昼夜项随自转方向可正可负，季节项则让半长轴减小。式 (3) 后的定义把自转频率接到热扩散上；单改密度，还不能指定温度滞后。

直接的思想检验是把自转轴翻转：简单昼夜项的轨道漂移应反向；若只保留季节项，漂移仍向内。再与 PR 拖曳比较，它的能量损失并不依靠下午侧较热。

$$
\Gamma=\sqrt{k\rho C},\qquad \ell_\nu=\sqrt{\frac{k}{\rho C\nu}},\qquad \Theta_\nu=\frac{\Gamma\sqrt\nu}{\epsilon\sigma T_*^3}
$$

依次为热惯量、扩散长度与比较周期加热和辐射冷却的热参数。昼夜项取自转频率，季节项取公转频率；这里沿用固定综述的定义。

对应阅读：Vokrouhlický et al. · Burns

## 图 2 与图 3，找出真正测量的量

图 2 把 Golevka 的轨道预测放在雷达距离和距离变化率坐标中。先找纯引力预测与雷达测量，再看加入热力后的预测。图中不确定度属于轨道预测，不是直接拍摄热光子的照片。

图 3 比较 1992 BF 的旧照相底片轨迹与回推天区位置。先认清赤经、赤纬方向，再比较模型轨迹。两图都利用累积的位置后果检验力模型；长期半长轴漂移由多期位置与动力学共同推断。

对应阅读：Vokrouhlický et al. · JPL Horizons API and observer/vector conventions

## 第 2.2 节与表 1，追问解释依赖什么

第 2.2 节讨论更细的热与形状计算。表 1 的轨道拟合误差是在采用的观测和动力学下回答统计问题；要作物理解释，还需合理的尺寸、密度、自转和热响应。比较天体的精度前，先看观测弧段和雷达覆盖。

读下一篇论文时，记下版本日期、实测量、力学模型、材料与自转信息，以及其他可能加速度。横向拟合项与物理成因要分开：彗核的 Marsden 出气参数和小行星热漂移都可能含横向项，描述的却不是同一种力。ISON 教程将很小的 A1/A2 出气采样贡献，与热物理分开检验。

对应阅读：Vokrouhlický et al. · JPL Horizons API and observer/vector conventions · Li et al.

## 亲手观察

### [把 Lab 作为光压与 PR 参考](https://sandboxphysics.com/zh/labs/comet-dust-tails/)

1. 选定粒径，用上一章公式算圆轨道时间，再对照下方图的真实年份。
2. 比较论文热模型前，先指出交互场景明确省略了哪些项。

## 练习与答案

模型只有力比和一个“热反冲”开关，还缺哪些物理输入？

还缺温度演化、表面边界条件、自转、形状、热物性与发射动量定律，才能预测反冲向量。力比比较辐射压力与引力，并不给出温度分布；缺少这些内容，开关便没有定义清楚的热预测。

这条课程实际复算了哪些结果？

配套教程复算固定 HST 方向诊断、JPL 投影与明确采用的太阳相对释放几何；PR 章节复算圆轨道解析参考，并用直接积分比较。综述里的小行星探测和热形状反演在此作原文阅读，没有重跑其原拟合。

[上一篇：把固定 JPL 尘埃几何投影到真实 HST 的 ISON 图像](https://sandboxphysics.com/zh/stories/comet-dust/ison-image-projection-tutorial/)

## 来源与计算范围

1. [Burns, Lamy & Soter, Radiation Forces on Small Particles (1979)](https://ntrs.nasa.gov/citations/19790070469)

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

2. [Li et al., HST observations of Comet ISON, arXiv:1311.0826v1](https://arxiv.org/abs/1311.0826v1)

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

3. [JPL Horizons API and observer/vector conventions](https://ssd.jpl.nasa.gov/horizons/manual.html)

   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；保存的请求与响应保留时标和几何选择。

4. [JPL Solar System Dynamics astronomical parameters](https://ssd.jpl.nasa.gov/astro_par.html)

   天文单位、日、光速以及相对力检验所用行星系统引力参数。太阳引力参数 132712440041.93938 立方千米每平方秒、光度三点八二八乘十的二十六次方瓦固定自 Horizons 太阳元数据。教学密度与压力效率是设定，不是 ISON 尘粒实测值。

5. [Vokrouhlický et al., The Yarkovsky and YORP Effects, arXiv:1502.01249v1](https://arxiv.org/abs/1502.01249v1)

   固定 2015 年 v1，共 24 页。第 2.1 节式 (1)—(3) 给出昼夜与季节性半长轴漂移；第 2.2 节讨论热物性和形状；第 3.1 节及图 2—3 讨论 Golevka、1992 BF 等天体测量探测。本页导读原文，不添加新热模拟，也不声称汇总最新探测。