[论文解读] Self-similarity of the Third Type in the Strong Explosion Problem

Andrei Gruzinov|arXiv (Cornell University)|Mar 11, 2003

Marine and environmental studies参考文献 1被引用 4

一句话总结

本文提出了一种新型自相似解——III型，用于球对称、幂律密度介质（指数 $\alpha$ 在 3 到 3.26 之间）中的强爆炸问题。与 I 型（衰减激波）和 II 型（需特征值解的加速波）不同，该解具有恒定速度激波（$\beta = 1$），由非自相似活塞驱动，尽管在后激波区域收敛性有限，数值上仍被确认为长时间渐近解。

ABSTRACT

Propagation of a blast wave due to strong explosion in the center of a power-law-density ($ρ\propto r^{-α}$) spherically symmetric atmosphere is studied. For adiabatic index of 5/3, the solution was known to be self-similar, (of type I) for $α<3$, self-similar (of type II) for $α>3.26$, and unknown in between. We find a self-similar solution for $3

研究动机与目标

• 解决在 $3<\alpha<3.26$ 范围内，强爆炸问题在幂律密度介质中自相似解的缺失问题，该区域此前未被解释。
• 提出一类新的自相似解——III 型，与 I 型（能量主导，$\beta<1$）和 II 型（特征值主导，$\beta>1$）截然不同。
• 证明当 $\alpha \in (3, 3.26)$ 时，激波以恒定速度（$\beta=1$）传播，这是由非自相似气体产生的活塞效应所致。
• 提供初步的数值证据，表明该 III 型解是强爆炸的真正长时间渐近解，尽管存在数值不稳定性。

提出的方法

• 使用无量纲变量 $\xi = r / R(t)$，其中 $R(t) = t^\beta$，从流体动力学方程推导出密度 $\rho(\xi)$、速度 $v(\xi)$ 和压力 $p(\xi)$ 的常微分方程。
• 引入变换变量 $V(\xi) = v(\xi)/\xi$ 和 $C(\xi)$，使得 $p(\xi) = \rho(\xi) (C(\xi) \xi)^2$，将系统简化为 $V(\xi)$ 和 $C(\xi)$ 的两个常微分方程。
• 通过从 $\xi=1$（对应激波前缘）反向数值积分求解常微分方程组，初始条件为 $V(1)=C(1)=1$，并检查在 $H=0$ 处奇异性是否被消除。
• 对于 III 型，解在 $V=4/3$，$C=0$（即 $H=0$）处终止，且该点处 $F=0$ 的条件可确保正则性，无需求解特征值问题。
• 使用一阶球对称欧拉格式模拟幂律大气中的强爆炸，通过最多 20 次自适应网格加倍以达到大半径。
• 将数值得到的密度、速度和压力分布与预测的自相似解进行比较，重点关注 $\alpha \approx 3.15$ 以检验 III 型区域。

实验结果

研究问题

• RQ1 在 $3 < \alpha < 3.26$ 的幂律密度介质中，强爆炸问题是否存在自相似解？该区域此前既不适用 I 型也不适用 II 型解。
• RQ2 在这一中间区域，何种物理解释机制可实现恒定速度（$\beta=1$）激波？其与能量主导或特征值主导自相似性的区别何在？
• RQ3 所提出的 III 型解是否如活塞驱动机制所暗示的那样，是强爆炸的真实长时间渐近解？
• RQ4 为何在此区域激波后部的数值收敛性会失效？这是否影响对解有效性的信心？

主要发现

• 在 $3 < \alpha < 3.26$ 范围内发现了一类新型自相似解——III 型，其特征为恒定激波速度（$\beta = 1$），既非能量主导也非特征值决定。
• 解在 $V=4/3$，$C=0$（即 $H=0$）处终止，且仅当该点处 $F=0$ 时奇异性可被消除，而 $\beta=1$ 时该条件自然满足，无需求解特征值问题。
• 对于 $\alpha \approx 3.15$，使用高达 $10^5$ 个网格点和 18 次网格加倍的数值模拟，在密度、速度和压力分布上与预测的 III 型分布具有良好一致性。
• 数值方法成功再现了已知的 I 型和 II 型解作为基准，证实了其在测试案例中的可靠性。
• 尽管分布匹配良好，但在高分辨率下激波后部区域出现数值不稳定性，表明可能存在弱不稳定性或收敛性问题，略微降低了对渐近性声明的信心。
• 该解在物理解释上是合理的：激波后方的非自相似气体作为活塞，维持了恒定速度的自相似前缘。