---
title: "泡壁怎样推动等离子体"
description: "把泡壁分成爆燃、混合与爆轰三类，求出 Jouguet 速度，算出释放的能量有多少变成流体运动。"
canonical_url: "https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/how-bubble-walls-push-the-plasma/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/phase-transition-gravitational-waves/how-bubble-walls-push-the-plasma/"
last_modified: "2026-11-12"
---

[宇宙相变引力波](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/)

# 泡壁怎样推动等离子体

Sandbox Physics · 2026-11-12

泡壁一边推进，一边把真空能变成等离子体的热和整体运动。有多少变成运动，取决于相变强度$\alpha$，也取决于泡壁比声速快还是慢。取$\alpha=0.2$、$\xi_w=0.9$，单个泡泡的效率是$\kappa=0.274$，模板里保留的动能占比是$K=0.0274$。

## 开始前需要什么

先读成核那一章。需要知道气体的声速、相对速度，会做带平方根的代数。流体方程只引用、不推导，解是数值求出的。

## 相变强度与袋模型

把等离子体看成相对论性气体，压强是能量密度的三分之一，声速是$c_s=1/\sqrt3\approx0.577$。袋模型里，假真空比破缺相多带一份恒定的能量密度$\epsilon$。相变强度拿它和泡壁前方的辐射能量密度比：$\alpha=\epsilon/(a_+T_+^4)$，$a_+$是定义辐射能量的常数。

泡壁以速度$\xi_w$穿过远处静止的等离子体。除此之外没有别的尺度，所以单个膨胀泡泡周围的流动只依赖$\xi=r/t$，是自相似的。一条剖面解出来，就描述了泡泡在碰到别的泡泡之前的每一时刻。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, arXiv:1910.13125v2 / JCAP (2020), doi:10.1088/1475-7516/2020/03/024

## 三种泡壁

比声速慢的泡壁叫爆燃。它推着前方的等离子体走，一道激波以略高于声速的速度跑在前面，激波和泡壁之间的等离子体向外流。

比某个极限速度、即 Jouguet 速度$\xi_J$还快的泡壁叫爆轰。前方的等离子体完全没被扰动，泡壁撞上它、拖着它走；后面跟着一道稀疏波，让流体到$\xi=c_s$时重新慢下来。

速度在$c_s$和$\xi_J$之间时，解是混合型：前有激波，后有稀疏波。$\alpha=0.2$时 Jouguet 速度是$0.828$，泡壁速度 0.9 就是爆轰；$\alpha=0.01$时它只有$0.653$。

$$
\xi_J(\alpha)=\frac{\sqrt{2\alpha/3+\alpha^2}+\sqrt{1/3}}{1+\alpha}
$$

$\alpha\to0$时它趋于$c_s$；相变越强，这条分界越往光速靠。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

## 多少能量推动了流体

效率$\kappa$是释放的真空能里，最终变成单个泡泡流动的整体动能的那一份。它是$w\,v^2\gamma^2\xi^2$在自相似剖面上的积分，$w$是焓密度，$v$是流速，$\gamma$是洛伦兹因子。Espinosa、Konstandin、No 和 Servant 在 2010 年解出了这些剖面，并给出拟合公式，说明拟合与解的偏差小于 15%。

流体动能还要换算成占总能量密度的比例。2024 年的 LISA 模板再乘上 0.6，这个系数来自多泡泡模拟：泡泡多了以后，整体运动比孤立泡泡要小。

算个例子。$\alpha=0.2$、$\xi_w=0.9$时，拟合给出$\kappa=0.2742$，我们自己解剖面得到$0.2740$，于是$K=0.6\times0.2742\times0.2/1.2=0.0274$。慢泡壁算例$\alpha=0.1$、$\xi_w=0.4$，拟合给$0.1877$，我们的解是$0.1890$，$K=0.0102$。

$$
K\simeq0.6\,\frac{\kappa\,\alpha}{1+\alpha}
$$

0.6 是 arXiv:2403.03723v2 采用的多泡泡修正。2020 年的 LISA 综述 arXiv:1910.13125v2 写的是$K=\kappa\alpha/(1+\alpha)$，没有这个系数；同样的$\alpha$和$\xi_w$，两篇文章给出的$K$不一样。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

## 我们的剖面解，以及它在哪里停下

我们自己积分了自相似流体方程，在泡壁和激波处用衔接条件，覆盖从 0.01 到 0.5 的六个强度、从 0.2 到 0.95 的泡壁速度。离 Jouguet 速度不太近时，和已发表拟合的最大相对偏差是 7.4%，出现在$\alpha=0.01$、$\xi_w=0.2$。

有两类点没有剖面。离 Jouguet 速度 0.01 以内，我们的打靶法不收敛，这些点就空着，不用拟合值去填。泡壁慢、相变又强时，根本不存在爆燃解：$\xi_w=0.4$时有解的最大强度是$\alpha=0.590$，$\xi_w=0.2$时是$0.396$。超过这个界，实验室不画谱。

声壳厚度也来自剖面。模板用$\xi_{\rm shell}=|\xi_w-c_s|$，这对亚声速爆燃和爆轰成立。混合型的声壳更薄，模板没有为它标定，所以实验室把从$c_s-0.05$到$\xi_J$这一段标为超出适用范围。

![三种相变强度下，单个泡泡的动能效率随泡壁速度的变化](https://sandboxphysics.com/_astro/efficiency.DeId2Pne.svg)

$\alpha = 0.03$、$0.1$ 和 $0.3$ 时的效率 $\kappa$。曲线是 Espinosa 等人 2010 年发表的拟合，圆点和方块是我们在袋模型里自己解出的剖面。虚线是声速 $c_s = 1/\sqrt{3}$。靠近 Jouguet 速度、我们的打靶法不收敛的地方没有画点。 [查看完整图像](https://sandboxphysics.com/_astro/efficiency.DeId2Pne.svg)

对应阅读：Original benchmark tutorial materials (standard-library Python) · Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010)

## 练习与解释性答案

1. $\alpha=0.05$时，速度$\xi_w=0.7$的泡壁是爆轰吗？$\alpha=0.01$呢？

   展开解释

   算$\xi_J$：$\alpha=0.05$时是$0.730$，$\alpha=0.01$时是$0.653$。0.7 的泡壁在前一种情况是混合型，在后一种情况是爆轰。同样的速度换了类型，是因为相变越强，泡壁要跑得越快才能不扰动前方的等离子体。

2. 用$\kappa=0.1545$算$\alpha=0.1$、$\xi_w=0.9$时的$K$，和同一强度的慢泡壁比一比。

   展开解释

   $K=0.6\times0.1545\times0.1/1.1=0.00843$，比$\xi_w=0.4$时的$0.0102$略低。快爆轰并不自动更高效：很多能量用来加热泡壁后面的等离子体。效率在 Jouguet 速度附近最高。

3. 为什么实验室不能直接对$\alpha=0.5$、$\xi_w=0.25$用拟合公式？

   展开解释

   拟合公式照样给出一个数，可袋模型里那里不存在稳定的爆燃解；我们的求解器给出$\xi_w=0.25$时的界是$\alpha=0.429$。没有物理流动支撑的数字会误导人，实验室就给出提示、不画谱。

## 动手检验预测

### [把泡壁速度扫过三种类型](https://sandboxphysics.com/zh/labs/phase-transition-gravitational-waves/)

1. 先预测$\alpha=0.2$时效率曲线的峰在哪里，再把泡壁速度从 0.4 慢慢拖到 0.95。
2. 看画面标签从“激波跑在泡壁前面”变成“稀疏波拖在泡壁后面”，中间有一段不画谱。
3. 把泡壁速度放在 0.3，再把$\alpha$调到 0.5，读一读出现的提示。

[上一篇：泡泡怎样冒出来，又怎样填满空间](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/how-bubbles-appear-and-fill-space/) [下一篇：从声波寿命算到今天的引力波谱](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/from-sound-lifetime-to-the-spectrum-today/)

## 来源与计算范围

1. [Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020](https://arxiv.org/abs/2403.03723v2)

   读的是 arXiv v2 的 LaTeX 源文件：摘要与引言、第 2.2.1 节（声波双折幂律、振幅、平均泡间距与红移因子）、第 3 节和附录 A（双参数 LISA 噪声与响应），以及第 5 节开头关于参数重建的部分。噪声公式把光学测量拐点印成 2e-3 mHz；我们的资料用 2 mHz，并同时给出两种结果。图 2 只读了图注。

2. [Caprini et al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, arXiv:1910.13125v2 / JCAP (2020), doi:10.1088/1475-7516/2020/03/024](https://arxiv.org/abs/1910.13125v2)

   只读了 LaTeX 源的两处：声波信号一节（signal.tex 约第 180 至 260 行），以及均方根流体速度和动能占比的定义（prelim.tex 第 171 至 182 行）。只用来对比早期的激波形成时间定义。

3. [Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010)](https://arxiv.org/abs/1004.4187v1)

   读了单个自相似泡泡的流体方程、爆燃、混合与爆轰三类剖面，以及附录 B 的效率拟合和其“优于 15%”的精度说明。我们在袋模型里自己解出的剖面独立检验这些拟合。

4. [Original benchmark tutorial materials (standard-library Python)](https://sandboxphysics.com/stories/phase-transition-gravitational-waves/materials/README.md)

   我们的脚本、三个固定基准、四个失败案例和预期输出。它计算已发表的参数化模板、效率拟合和噪声模型；不求解爱因斯坦方程或流体方程，也不含浏览器引擎代码、论文原文或 LISA 数据。

## 继续阅读

宇宙相变引力波 · 阅读路径

### 从泡泡成核到 LISA 可能看到的谱

1. [文章 **宇宙早期若冒过泡泡，会留下怎样的引力波**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/what-a-bubbling-early-universe-leaves-in-gravitational-waves/)
2. [文章 **泡泡怎样冒出来，又怎样填满空间**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/how-bubbles-appear-and-fill-space/)
3. [当前页面 **泡壁怎样推动等离子体**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/how-bubble-walls-push-the-plasma/)
4. [文章 **从声波寿命算到今天的引力波谱**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/from-sound-lifetime-to-the-spectrum-today/)
5. [文章 **独立复算一份基准声波谱**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/recompute-a-benchmark-sound-wave-spectrum/)
6. [文章 **读懂 LISA 声波模板和它的边界**](https://sandboxphysics.com/zh/stories/phase-transition-gravitational-waves/read-the-lisa-sound-wave-template-and-its-limits/)