A model of a multiphase medium based on the closure of moment chains for the Vlasov-Poisson equations

技术摘要：基于 Vlasov-Poisson 方程矩链闭合的多相介质模型

问题陈述
从描述无碰撞等离子体的动力学 Vlasov-Poisson 模型向流体动力学模型的过渡，是电子气建模中的一个基本挑战。虽然通过速度平均矩（密度、速度、压力）推导出的流体动力学方程可以得出守恒律，但所得系统本质上是不封闭的。具体而言，低阶矩的演化依赖于高阶矩（例如代表热通量的三阶矩）。传统的闭合方法试图将这些高阶矩表示为低阶矩的函数，但往往无法保留动力学模型的本质属性，如兰道阻尼（Landau damping）、湍流或扰动的有限传播速度。此外，将矩链扩展到更高阶往往会导致所得系统失去双曲性，从而使其对于标准数值方法而言是不适定的。

研究方法
本文针对一维无碰撞 Vlasous-Poisson 方程提出了一种替代性的保守矩闭合方案。该方法依赖于以下数学和物理构造：

1. 矩链构建： 作者定义了源自分布函数 $F$ 的 Hamburger 矩 $M_i(t, x) = \int v^i F(t, x, v) dv$。动力学方程产生了一个矩方程链，其中 $M_i$ 的时间导数取决于 $M_{i+1}$ 的空间导数以及电场 $E$。
2. 基于 Hölder 不等式的闭合原理： 本文方法的核心是在奇数步 $j = 2k - 1$ 处闭合链。闭合条件是通过 Hölder 不等式的等号成立情况导出的：$(M_{2k-1})^2 = M_{2k}M_{2k-2}$。这意味着将项 $\lambda_j$（代表偏离等号的情况）设为零。
3. 递归定义： 一旦链在第 $j$ 步闭合，所有后续的高阶矩（$M_{j+1}, M_{j+2}, \dots$）都作为低阶矩的逻辑后果被递归地定义，而不是被丢弃。
4. 多相解释： 闭合系统被解释为一种多组分（多相）流体。每个“相” $i$ 对应于一对关于质量平衡和线性动量的方程，其特征为密度 $n_i$ 和速度 $U_i$。至关重要的是，这些相之间的耦合是通过电场实现的。闭合原理规定，最后一相（以及所有后续相）是无压力的（$P_k = 0$）。

主要贡献与结果

• 双曲性保持： 作者证明了所得方程组在每一个闭合步骤中都是非严格双曲的。系统的系数矩阵具有实谱，其谱由一个零值（重数为 $j-1$）和一个非零值（重数为 2）组成。然而，该系统缺乏完备的特征向量集，这意味着它是非严格双曲的。尽管如此，实谱确保了扰动的有限传播速度，这是其他闭合尝试中经常丢失的属性。
• 保守形式： 闭合系统被成功地改写为保守形式（记作系统 MF）。对于第 $k$ 步闭合，系统由 $k$ 对流体方程及一个电场方程组成。
• 无限守恒律： 本文证明了闭合系统允许存在无限级数守恒律。这些守恒律直接从矩方程和电场演化中导出，确保了特定矩组合与场能的守恒。
• 热通量表达式： 该方法提供了热通量 $q$（第三阶中心矩）关于闭合系统矩的显式表达式。这使得可以通过求解导出的流体动力学系统来计算热通量，而无需进行人为假设。
• 特定情形：
  • 第一次闭合 ($k=1$)： 退化为冷等离子体方程，其分布函数为 Dirac $\delta$ 函数。
  • 第二次闭合 ($k=2$)： 产生一个涉及两个流体组分的系统，其中第二组分为无压力的。
  • 分布假设： 文中指出，如果假设分布函数为高斯分布（热核），则闭合将导致零热通量的等熵气体动力学。

意义与主张
作者将这项工作定位为从动力学建模向流体动力学建模过渡的一种替代方法。他们明确声明，并不声称该方法完全解决了在流体动力学框架下重现所有动力学特性（如兰道阻尼或湍流）的问题。

然而，本文声称该方法具有以下有用的特征：

1. 它允许在任何奇数步进行闭合，同时保持方程的双曲性质。
2. 闭合过程并未丢弃无限的高阶矩方程序列；相反，它将它们确立为低阶系统的逻辑后果。
3. 它提供了一个严谨的数学框架（基于 Hölder 不等式的等号成立），并可以赋予物理阐释，即通过电场相互作用的多相介质。

这项工作贡献了一个数学一致的闭合方案，该方案保持了方程的双曲结构（尽管是非严格的），为在传统闭合方法难以维持适定性的无碰撞等离子体建模提供了一条潜在路径。