公共包层释放的能量去了哪里

先问热能来不来得及做功

伴星在巨星包层中运动，压缩气体并产生激波。留下来的热可以提高压力，让气体膨胀；若辐射在此之前就把能量送到别处，同一团气体受到的压力推动会变弱。于是，轨道释放能量相近的两次计算，仍可能有很不同的气体分布和抛出质量。

这也是标量能量终态留下的问题。一个效率参数汇总了多少轨道能用于某项成本，却没有追踪能量在哪一层、什么时候转成压力功，或以什么方式离开。阅读辐射流体模拟时，先找出这条能量路线，再看最终数字。

本次阅读固定为Lau等人的 $mathrm{arXiv}:2503.20506mathrm{v2}$ 。它比较12倍太阳质量红超巨星与3倍太阳质量点质量伴星的三维计算：一组绝热，一组含辐射扩散，另一组在辐射计算中另加核心附近的热输入。主要分辨率为两百万个SPH粒子。

先读第2节，记住这些并不是交互场景那颗学校巨星的“更高分辨率版”。论文的供体不自转，初始包层用平熵构造，并未驱动原先就存在的平衡对流；内侧结构由软化点质量代替。我们另有的15倍太阳质量学校剖面，质量和结构都不同，不能因为也叫红超巨星就拿来对齐成同一条演化序列。

额外加热组用于比较中心附近的热结构：它帮助维持密度和熵，却并未被作者选为最真实的默认组。初始对流输运被抑制，而热又从中心加入，这个热结构本身需要解释。中心加热、对流与表面冷却，需要共同定义一个可解释的热初态。

图6画包层中被计为非束缚的质量比例，图7画两致密中心的间距。先比较相同曲线颜色，再沿时间分清快速卷入和之后较慢的旋入。若只挑最后一帧，很容易把一个仍在演化的慢阶段当成永久终态。

图6实线用机械比能为正作为非束缚判据；虚线把气体和辐射内能也算进去。二者不能混着比较：

e_{\rm k}+e_{\rm p}>0,\qquad\hbox{或}\qquad e_{\rm k}+e_{\rm p}+e_{\rm gas}+e_{\rm rad}>0

同一团气体可以满足第二种判据，却还没把内能转成足够的向外机械运动。继续冷却或做功，会改变后续命运。

论文给出的后期例子中，绝热组非束缚比例约20.8%，辐射组约11.0%，加热的辐射组约8.9%。这些是特定时段、特定判据下的比较。“约减半”可以帮助抓住量级，但还应回到相同时间的曲线读数，不能把它当作每个时刻都严格成立的比例。

与这种显著的质量差异相比，快速卷入结束后最初一段时间的间距演化较接近。气体留下多少热、分布在哪里，与双星当时缩到哪里，不是同一个输出。这个结果把问题推向较慢、非绝热的后续阶段，还需要继续追踪包层和轨道。

图8把密度、局部扩散时间、机械比能、不透明度和氢电离状态放在同样的空间切片里。先选同一列演化阶段，再横向比较位置：哪里气体致密，哪里辐射能量梯度陡，哪里可能在气体运动以前就重新分配辐射能？

t_{\rm diff}={\kappa\rho\over c\lambda_{\rm FLD}}\left({E_{\rm rad}\over|\nabla E_{\rm rad}|}\right)^2

这是论文使用的局部辐射能输运时间。

lambda_{ m FLD}

是通量限制因子，与上一章描述包层重力成本的lambda不是同一个量。

在光学厚的地方，辐射反复相互作用，扩散近似有清楚的物理图像。在光学薄的地方，通量限制避免让能量流超过自由传播的尺度，但它并不恢复完整的方向性辐射场。气体压力占主导时，辐射能的搬运时间也不自动等于整团气体的热调整时间。

氢部分电离和不透明度的图很吸引人，仍要回到第2.3节检查状态方程：这些计算没有加入原子、分子复合所释放的能量。它们可以讨论复合区附近的光学厚度，却不能因此改名为“复合能驱动的抛射模拟”。

图注还有一处符号不一致：固定v2的图8把“bound material”旁写成 $e_{ m kp}>0$ ，与第3.2节正文正机械能表示非束缚的定义不一致。这里按正文定义解释能量，不从这个标签矛盾推断其实际代码已经反向筛选。把它记为固定版本的图注疑点，再沿正文的能量定义继续读图。

现在翻到附录A和图9。作者报告，到约8年时，两组辐射计算的总能量漂移分别达到10.7%和6.9%；后一项已把明确注入的加热另作处理。主要漂移发生在卷入以后，与隐式辐射求解器不严格守能有关。

同一对计算的角动量保存却好得多，分别约0.020%和0.022%。轨道角动量与总能量各自保留自己的误差记录，辐射求解的误差还要沿相应的能量交换项追踪。也不能把总能量漂移直接当作每个输出都具有相同百分比误差，再给非束缚质量做一项统一修正；哪些气体被误差影响，需要空间与时间信息。

交互能量方程的残差可以接近机器精度，因为它只检验一个代数等式。论文则在解气体、引力和辐射的耦合问题。把二者的小数位排在一起比较“谁更准”，会抹掉它们真正计算的问题。

附录B把主要计算同少十倍粒子的试验比较。先看图10的间距：快速卷入后最初约一至两年，几个计算保持相近的变化趋势；但辐射组约8.5年出现的晚期明显下降，在较低分辨率里没有出现。作者保留了随机湍动或数值伪影的可能，不能把这个下降单独当成确定的长期轨道命运。

图11再比较非束缚质量。晚期羽流阶段，两种分辨率的差异扩大。与此同时，两种分辨率的绝热组都在卷入后持续抛出质量，加入辐射输运以后，抛出量都明显减弱。这份证据支持有范围的物理趋势；晚期细节仍需要更细分辨率与更长演化的比较。

读到这里，可以把自己的结论写得更有内容：在这组初态、状态方程和输运近似下，辐射改变气体保留热并膨胀的能力；较早的卷入后间距对照与晚期细节具有不同的数值证据。若要讨论其他质量、旋转、真实对流初态和复合能，仍要做新的比较。