---
title: "两道激波与另一种点爆炸问题"
description: "从膨胀抛射物推导包层相似律，在接触面连接两层气体，再比较独立求出的Sedov–Taylor点爆炸解。"
canonical_url: "https://sandboxphysics.com/zh/stories/supernova-remnants/two-shocks-and-a-point-blast/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/supernova-remnants/two-shocks-and-a-point-blast/"
last_modified: "2026-11-04"
---

[超新星遗迹、激波与高能粒子](https://sandboxphysics.com/zh/stories/supernova-remnants/)

超新星遗迹 · L2

# 两道激波与另一种点爆炸问题

Sandbox Physics · 2026-11-04

两道激波对应两股不同的来流：内部是高速膨胀的抛射物，外部是环境气体。把两层激波后的气体接起来，得到抛射物包层阶段的解。Sedov–Taylor点爆炸则从零抛射物质量和固定能量出发，是另一个有明确条件的问题。

## 开始前需要什么

先读强激波一章。这里会解释自相似、无量纲半径，以及怎样把球壳的质量和能量加起来。

## 膨胀抛射物带着速度标签

自由膨胀时，未激波物质近似保持各自的速度，所以 $v=r/t$。速度大的物质在更外面。设外包层密度随这个速度的 $n$ 次反比下降；同一段速度范围占据的体积随 $t^3$ 增大，密度就随 $t^{-3}$ 稀释。换成半径坐标，便得到下面的密度式。

其中 $g^n$ 是固定包层归一化的一个有量纲常数。外部环境密度 $\rho_0$ 均匀，气体起初静止。仍使用上游冷、单原子理想气体、非辐射的球对称近似。

怎样让两道激波一直按相同比例膨胀？假设三个面的半径都随 $t^m$ 变化，来流相对速度便都随 $t^{m-1}$ 变化。激波后压力大致等于上游密度乘相对速度平方，要让两侧压力维持同一种相似关系，两股来流的密度比就得固定。把半径代入包层密度，时间指数为 $n-3-mn$；令它为零，得到 $m=(n-3)/n$。

$$
\rho_{\rm ej}=g^n t^{\,n-3}r^{-n},\qquad R_{\rm c}=\left(\frac{Ag^n}{\rho_0}\right)^{1/n}t^{(n-3)/n}
$$

接触面半径 $R_{\rm c}$ 有长度单位。无量纲系数 $A$ 由两层实际激波后气体的匹配确定。

对应阅读：Chevalier (1982), Self-similar solutions for the interaction of stellar ejecta with an external medium, ApJ258,790–797

## 两层气体在接触面相遇

在正向与反向激波处，分别使用上一章的跳跃条件，再各自向接触面积分。相遇处必须有相同的接触面位置、压力，以及推动接触面的速度。两侧密度可以不同，物质也不会自由穿过这道分界。

独立球对称Euler积分给出 $n=7$ 时的 $R_{\rm f}/R_{\rm c}=1.18109$ 和 $R_{\rm r}/R_{\rm c}=0.934985$。改取 $n=12$，会得到不同的壳层相对厚度。离线计算还对照原论文表1的舍入范围，并另外积分被扫过物质的质量与能量。

理想接触面附近，某些无量纲量会出现奇异行为。数值积分在声明的正距离处停止，再比较不同容差和截断距离；密度与温度没有被直接插值穿过接触面。

交互画面采用幂律指数为七时的半径、激波紧邻两侧的跳跃值，以及已知未激波区域的密度。独立读者程序会积分激波后的壳层来做科学检查；画面没有从几个边界值臆造完整的内部密度图。

[![自有Euler积分得到的幂律指数为七、十二的两组激波半径，接触面分别归一化为一。](https://sandboxphysics.com/stories/supernova-remnants/materials/double-shock-radii.svg)](https://sandboxphysics.com/stories/supernova-remnants/materials/double-shock-radii.svg)

自有独立Euler结果。R指反向激波，C指接触面，F指正向激波。每组接触面都取一，颜色区分几何。点击SVG可放大读图。

对应阅读：Chevalier (1982), Self-similar solutions for the interaction of stellar ejecta with an external medium, ApJ258,790–797 · Sandbox Physics independent supernova-remnant Python science reader

## 幂律包层不能一直延伸到中心

包层律只描述外部抛射物。若把它延伸到零速度，中心质量就会发散。要构造有限质量的例子，可以在速度 $v_t$ 处连续接上平密度内核。利用 $r=vt$ 和 $dr=t\,dv$，对不同速度的球壳积分，时间因子会恰好消去。

设内核密度为 $Ct^{-3}$，包层为 $Ct^{-3}(v/v_t)^{-n}$。内核与包层的质量项分别含 $1/3$ 与 $1/(n-3)$，动能项则含 $1/5$ 与 $1/(n-5)$。要求 $n>5$，质量和动能才能都有限。由两种积分的比值，得到下面的转折速度。

包层相似解的适用条件是 $R_{\rm r}/t>v_t$，也就是反向激波还没有进入内核。固定年龄的实验会检查短粒子窗口的起点和终点。环境密度继续升高，最终可能破坏这个条件；那时应拒绝输入，而不是沿用包层公式继续画。

$$
v_t^2=\frac{10(n-5)E}{3(n-3)M},\qquad C=\frac{M}{4\pi v_t^3\left[1/3+1/(n-3)\right]},\qquad g^n=Cv_t^n
$$

$E$ 是给定的初始总动能，$M$ 是有限抛射物质量。这个归一化给出外包层，尚未求出以后反向激波进入内核的演化。

对应阅读：Chevalier (1982), Self-similar solutions for the interaction of stellar ejecta with an external medium, ApJ258,790–797 · Sandbox Physics independent supernova-remnant Python science reader

## 保持年龄，先做一个比例预测

对 $n=7$，能量、质量与年龄不变时，半径随 $\rho_0^{-1/7}$ 变化。把环境密度增大128倍；128恰好是二的七次方，因此三个半径都会减半。同一年龄下 $V=mR/t$，相应速度也减半。

若改成把年龄加倍，则是另一场实验。包层半径增加到 $2^{4/7}\simeq1.486$ 倍，速度变为 $2^{-3/7}\simeq0.743$ 倍。交互实验的密度预设始终保持约148.578年的同一个物理年龄，没有暗中换年龄来维持某个反向激波位置。

128倍是用于练习比例关系的例子，不保证两端都落在交互范围或包层阶段。把公式用到具体遗迹之前，还要检查反向激波与内核的关系。

对应阅读：Sandbox Physics independent supernova-remnant Python science reader

## 点爆炸从另一种守恒条件开始

独立的Sedov–Taylor问题，是把固定能量 $E$ 放进均匀密度 $\rho_0$ 的介质，忽略抛射物质量。起点没有内核，也没有反向激波。只比较单位，就能知道半径必须随 $(Et^2/\rho_0)^{1/5}$ 变化，但还差一个无量纲系数。

要确定这个系数，需要在强正向激波后积分球对称Euler方程，再把各球壳的动能和内能加起来，令总和等于给定能量。这一步提供了单靠量纲分析得不到的信息。

独立程序使用三个逐渐缩小的中心截断，单独报告未分辨中心的热能估计，得到 $\beta=1.1516664$、$\beta^5=2.0259724$。扫入的总质量和局部能流第一积分提供了另外的检查。公开读者包会重新计算这些数，不把现成结果当作重算。

包层解和点爆炸解的出发条件不同，可以比较指数，却不能把两条曲线前后接起来就称为完整演化。反向激波进入内核以后怎样过渡，必须另解动力学问题。

$$
R_{\rm ST}=\beta\left(\frac{Et^2}{\rho_0}\right)^{1/5},\qquad V_{\rm ST}=\frac{2R_{\rm ST}}{5t}
$$

这个零抛射物、能量不变的基准问题，与双激波包层解分开使用。

对应阅读：Leahy & Williams (2017), A Python Calculator for Supernova Remnant Evolution, arXiv1701.05942v2 · Sandbox Physics independent supernova-remnant Python science reader

## 练习与解释性答案

1. 在 $n=7$ 的包层解里，$V_{\rm r}/(R_{\rm r}/t)$ 是多少？为什么有反向激波？

   展开解释

   比值为 $4/7$。自由膨胀抛射物以 $R_{\rm r}/t$ 赶来，比前沿本身更快，实际流入的相对速度是 $3R_{\rm r}/(7t)$，这就是反向激波处理的来流。

2. 只用量纲分析，能确定点爆炸的全部系数吗？

   展开解释

   不能。它固定能量、密度与时间的幂次，却留下一个无量纲常数。还需要压力和速度分布，再结合总能量积分来确定。

3. 包层半径曲线与点爆炸曲线接得很光滑，为什么仍不足以证明过渡正确？

   展开解释

   一个半径不能指定整片气体的密度、速度、压力和能量分布。进入内核时的完整状态还必须满足动力学方程，视觉上的光滑连接并不够。

## 动手检验预测

### [检验同一年龄的缩放](https://sandboxphysics.com/zh/labs/supernova-remnant-shocks/)

1. 在允许范围内选两个密度，用七次方根的反比预测半径比，再比较接受后的读数。
2. 确认显示年龄相同，三个边界半径也一起变化。
3. 到方法篇重跑独立点爆炸，在纸上比较两种指数；把点爆炸结果和交互实验快照分开。

[上一篇：强激波怎样把运动变成热](https://sandboxphysics.com/zh/stories/supernova-remnants/strong-shocks-turn-motion-into-heat/) [下一篇：高能种子怎样在平均循环中继续加速](https://sandboxphysics.com/zh/stories/supernova-remnants/energetic-seeds-and-finite-acceleration/)

## 来源与计算范围

1. [Chevalier (1982), Self-similar solutions for the interaction of stellar ejecta with an external medium, ApJ258,790–797](https://ui.adsabs.harvard.edu/abs/1982ApJ...258..790C/abstract)

   原始全文与印刷公式、表1。分开求冷上游球对称Euler壳层，只用于大于五的包层阶段；没有接入内核演化或任意插值激波后密度。

2. [Leahy & Williams (2017), A Python Calculator for Supernova Remnant Evolution, arXiv1701.05942v2](https://arxiv.org/abs/1701.05942v2)

   核读第3.2.1、3.3.3节及印刷第16页系数。自有Euler点爆炸积分独立重算；只声称读过这些部分，没有声称读完此文或取得TM99全文。

3. [Sandbox Physics independent supernova-remnant Python science reader](https://sandboxphysics.com/stories/supernova-remnants/materials/supernova-remnants-reader-science-v1.zip)

   自有离线Euler n7/n12、独立零抛射物点爆炸、平均群体、角度与首次返回及Jensen检查和合成图。默认步骤与系数真正重算；不含原论文、原数据或App引擎。