自由膨胀状态的演化方程

在等离子体物理中，**自由膨胀（free expansion）**指的是等离子体在没有外部约束或驱动的情况下向真空或低密度介质自行展开的过程。该过程的动力学由一组耦合的演化方程描述，通常采用流体模型或动力学模型（如 Vlasov 方程）进行分析。下面系统阐述自由膨胀状态的基本假设、核心方程以及典型的自相似解。
在多数实验和数值研究中，针对自由膨胀过程常作以下简化：

• 无碰撞近似：等离子体尺度远大于碰撞长度，采用碰撞稀疏的 Vlasov 方程或无粘性流体方程。
• 单温度近似：电子与离子温度相同（或分别处理），忽略温度不平衡导致的热流。
• 一维平面对称：膨胀方向设为 (x) 轴，所有物理量仅随 (x) 与时间 (t) 变化。
• 电中性近似：在尺度远大于 Debye 长度时，等离子体整体保持准电中性 ((n_e \approx Z n_i))。

在上述假设下，可分别采用 流体模型（连续方程、动量方程、能量方程）或 动力学模型（Vlasov‑Poisson 方程）进行描述。下面先给出流体模型的完整形式。

2. 流体演化方程

2.1 连续方程

[
\frac{\partial n}{\partial t} + \frac{\partial}{\partial x}(n u) = 0
]

其中 (n(x,t)) 为粒子密度，(u(x,t)) 为流体速度。

2.2 动量方程

[
\frac{\partial}{\partial t}(n u) + \frac{\partial}{\partial x}!\left(n u^{2} + p\right) = 0
]

(p) 为热压强，若采用理想气体状态方程 (p = n k_{\mathrm{B}} T)。

2.3 能量方程（等熵近似）

在自由膨胀过程中，系统对外不做功且无热传导，可采用 等熵（adiabatic）关系：

[
\frac{\partial}{\partial t}!\left(\frac{p}{\gamma-1}\right) + \frac{\partial}{\partial x}!\left(\frac{p}{\gamma-1}u + p u\right) = 0
]

其中 (\gamma) 为绝热指数（电子-离子单温度等离子体常取 (\gamma = 5/3)）。

2.4 归一化

为便于数值求解，常将变量归一化：

[
\tilde{x} = \frac{x}{\lambda_D},\quad \tilde{t} = \omega_{pi} t,\quad \tilde{u} = \frac{u}{c_s},\quad \tilde{n} = \frac{n}{n_0},\quad \tilde{p} = \frac{p}{n_0 k_{\mathrm{B}} T_0}
]

其中 (\lambda_D) 为 Debye 长度，(\omega_{pi}) 为离子等离子体频率，(c_s = \sqrt{\gamma k_{\mathrm{B}} T_0/m_i}) 为声速。

3. 自相似解的构造

自由膨胀常呈现 自相似（self‑similar）结构，即物理量仅依赖于无量纲变量 (\xi = x / (c_s t))。将连续、动量、能量方程代入自相似假设，可得到一组常微分方程：

[
\begin{aligned}
&(1-\xi) \frac{d\tilde{n}}{d\xi} + \tilde{n}\frac{d\tilde{u}}{d\xi}=0,\
&(1-\xi) \frac{d\tilde{u}}{d\xi} + \frac{1}{\tilde{n}}\frac{d\tilde{p}}{d\xi}=0,\
&(1-\xi) \frac{d\tilde{p}}{d\xi} + \gamma \tilde{p}\frac{d\tilde{u}}{d\xi}=0.
\end{aligned}
]

利用状态方程 (\tilde{p}= \tilde{n}^{\gamma}) 消去 (\tilde{p})，可得到闭合的二阶方程，进而求得解析解：

[
\tilde{u}(\xi)=\frac{2}{\gamma+1}\bigl(\xi-1\bigr),\qquad
\tilde{n}(\xi)=\left[1-\frac{\gamma-1}{\gamma+1}\xi\right]^{\frac{2}{\gamma-1}},\qquad
\tilde{p}(\xi)=\tilde{n}^{\gamma}.
]

该解在 (\xi \le \xi_{\max}= \frac{2}{\gamma-1}) 区域有效，超出该范围密度趋于零，表现为等离子体前缘的 稀疏前沿。

示例：电子-离子单温度等离子体（(\gamma=5/3)）

代入 (\gamma=5/3)：

[
\tilde{u}(\xi)=\frac{3}{4}(\xi-1),\quad
\tilde{n}(\xi)=\left(1-\frac{1}{4}\xi\right)^{3},\quad
\tilde{p}(\xi)=\tilde{n}^{5/3}.
]

当 (\xi=0)（即膨胀起点）时，(\tilde{u}= -\frac{3}{4})，说明等离子体内部仍保持向左的微弱流动；而在 (\xi=\xi_{\max}=4) 处，(\tilde{n}=0)，对应膨胀前沿的真空界面。

4. 动力学视角：Vlasov‑Poisson 方程

在更高精度的描述中，需要考虑粒子分布函数 (f_s(x,v,t))（(s=e,i)），其演化满足 Vlasov 方程：

[
\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E,\frac{\partial f_s}{\partial v}=0,
]

其中电场 (E(x,t)) 由 Poisson 方程给出：

[
\frac{\partial E}{\partial x}= \frac{e}{\varepsilon_0}\bigl(Z n_i - n_e\bigr),\qquad
n_s(x,t)=\int f_s,dv.
]

在自由膨胀的极限（(E) 主要出现在稀疏前沿），可采用 薄层近似（thin‑sheath approximation），将电场视为局部常数，从而得到类似流体模型的自相似结构，但能够捕捉 离子加速 与 电子逃逸 的细节。

5. 数值实现要点

• 网格划分：采用自适应网格或拉格朗日跟踪，以保证稀疏前沿的分辨率。
• 守恒性：使用有限体积或有限差分守恒格式，确保质量、动量、能量的数值守恒。
• 边界条件：左端（等离子体内部）设为固定密度/温度，右端（真空）采用开放边界或吸收层，防止非物理反射。

# 简单的 1D Euler 代码框架（Python + NumPy）
import numpy as np

Nx = 400
x = np.linspace(0, 10, Nx)
dx = x[1] - x[0]
dt = 0.001
t_end = 5.0

# 初始条件：左侧均匀等离子体，右侧真空
n = np.where(x < 2.0, 1.0, 1e-6)
u = np.zeros_like(x)
p = np.where(x < 2.0, 1.0, 1e-6)

gamma = 5/3

def flux(n, u, p):
    return np.array([n*u,
                     n*u**2 + p,
                     (p/(gamma-1) + p)*u])

while t < t_end:
    # 计算通量
    F = flux(n, u, p)
    # Lax‑Friedrichs 更新
    n[1:-1] = 0.5*(n[2:] + n[:-2]) - dt/(2*dx)*(F[0,2:] - F[0,:-2])
    u[1:-1] = 0.5*(u[2:] + u[:-2]) - dt/(2*dx)*(F[1,2:] - F[1,:-2])/n[1:-1]
    p[1:-1] = 0.5*(p[2:] + p[:-2]) - dt/(2*dx)*(F[2,2:] - F[2,:-2])
    t += dt

上述代码演示了最简化的欧拉方程求解框架，实际应用时需加入 限制器（Limiter）以抑制数值振荡，并结合 自相似变量 进行后处理，以验证数值解与解析自相似解的一致性。

6. 小结

• 自由膨胀状态的演化方程在流体模型下由连续、动量、能量三组守恒方程组成，配合等熵关系即可得到完整闭合系统。
• 通过自相似变量 (\xi = x/(c_s t)) 可导出解析的自相似解，直观展示密度、速度、压强随膨胀前沿的空间分布。
• 对于更高精度的需求，Vlasov‑Poisson 方程提供了动力学层面的描述，能够捕捉电子逃逸和离子加速等微观效应。
• 数值实现时应关注守恒性、网格自适应以及合适的边界条件，以确保模拟结果可靠并与理论解相吻合。

掌握上述演化方程及其求解技巧，是进一步研究等离子体膨胀、激光等离子体相互作用以及等离子体推进器等前沿课题的基础。