Equilibria of the Pressureless Euler System in Dimension 1

技术摘要：一维无压 Euler 系统的平衡态

问题陈述
本文研究了一维无压 Euler 系统（PE）的长期行为，该系统最初由 Zeldovich 提出，用于模拟大规模结构的形成。该系统由连续性方程和动量方程组成：
$$\partial_t \rho + \partial_y(\rho v) = 0, \quad \partial_t(\rho v) + \partial_y(\rho v^2) = 0$$
其中 $\rho$ 是质量分布（具有有限二阶矩的概率测度），$v$ 是速度场。作者特别关注粘性粒子解（SPS），这是一类在碰撞时质量粒子会粘在一起的分布解。虽然这类解的存在性已通过投影公式（Natile-Savaré）得到确立，但这些解在何种条件下收敛于平衡态、平衡态的性质以及坍缩所需的时间，仍是悬而未决的问题。与吸引型 Euler-Poisson 系统不同（后者通常收敛于单个质点质量），无压 Euler 系统允许存在丰富的平衡分布。

方法论
分析在很大程度上依赖于 [10] 中推导出的投影公式，该公式将问题从物理上的“位置空间”转化为“质量空间”。通过利用最优传输映射 $N_t: (0,1) \to \mathbb{R}$（其中 $N_t \# \mathcal{L} = \rho_t$），系统的演化可以通过将线性轨迹 $N_0 + tV_0$ 投影到 $L^2(0,1)$ 中的非减函数凸锥 $K$ 上来表征：
$$N_t = \text{proj}_K(N_0 + tV_0)$$
其中 $V_0 = v_0 \circ N_0$。作者利用凸分析工具，特别是凸包（$\text{conv}$）的性质和右导数，来分析 $N_t$ 的渐近行为。他们定义了一个原函数 $F(x) = \int_0^x V_0(w) dw$，并利用系统的伽利略不变性将总动量归一化为零。

核心贡献与结果

1. 平衡态的充要条件：
   本文建立了一个仅基于初始数据 $(N_0, V_0)$ 的精确条件，记作 $(\star)$，用于判定 SPS 是否收敛于平衡态。该条件要求：

   • 总动量为零：$\int_0^1 V_0(w) dw = 0$（等价于 $F(1)=0$）。
   • 部分动量为非负：对于所有 $x \in (0,1)$，有 $\int_0^x V_0(w) dw \geq 0$（等效于 $F(x) \geq 0$）。

   作者证明，当且仅当条件 $(\star)$ 成立时，SPS 在 1-Wasserstein 度量下收敛于平衡态。如果对于任何 $x$ 都有 $F(x) < 0$，则解会发散（质量扩散至无穷远）。

2. 平衡态的特征描述：
   当满足条件 $(\star)$ 时，本文对极限分布 $\rho_\infty$ 提供了完整的描述。平衡态由函数 $F$ 的零点决定：

   • 静止区域： 在 $F$ 恒为零的任何区间上，质量分布在所有时间内保持静止（$N_t(x) = N_0(x)$）。
   • 坍缩区域： 在任何满足 $F(a)=F(b)=0$ 且在 $(a,b)$ 内部 $F > 0$ 的区间 $(a, b)$ 上，该区间内的质量会坍缩为位于该初始“块”质心处的单个 Dirac 测度（点质量）。
     这意味着平衡态是静止质量与坍缩点质量的叠加。

3. 坍缩时间：
   作者推导出了坍缩时间（$T_c$）的显式公式，其定义为使得对于所有 $t \geq T$，都有 $\rho_t = \rho_\infty$ 的最小时间 $T$。如果解收敛，则 $T_c$ 由 $F$ 为正的所有区间 $(a,b)$ 的上确界给出：
   $$T_c = \sup_{[a,b] \in I} \sup_{x \in (a,b)} \frac{\int_x^b N_0(w) dw + \frac{x-b}{b-a}\int_a^b N_0(w) dw}{\int_0^x V_0(w) dw}$$
   如果该值为无穷大，则解渐近地趋向于平衡态，但永远不会在有限时间内达到平衡。

4. 收敛度量：
   研究结果主要基于 1-Wasserstein 距离（$W_1$）进行构建。作者指出，如果解的支撑集保持有界（这发生在初始支撑集有界且满足 $(\star)$ 时），则在 $W_p$ 下的收敛性随之成立。作者讨论了向 $W_2$ 收敛的扩展，指出虽然在质量区间边界附近的特定单调性条件下 $W_2$ 收敛成立，但由于凸包操作的全局性质，仅依赖于初始数据的 $W_2$ 收敛的一般性条件仍是一个开放性的难题。

意义与主张
本文声称，通过仅使用初始条件，完整地刻画了一维无压 Euler 系统粘性粒子解的平衡行为。它解决了关于平衡态存在性、平衡态的具体形式以及坍缩时间的三个核心问题。作者强调，其方法避免了以往文献（如 [8]）中可能需要的关于随时间保持紧支撑的假设。这项工作通过区分于排斥型 Euler-Poisson 系统（后者由于病态行为导致此类预测方法失效），突显了无压 Euler 系统的动量驱动特性使得这种精确判定成为可能。研究结果依赖于投影公式和凸分析，在初始速度分布的几何结构与系统渐近状态之间建立了严谨的联系。