Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

技术摘要：三维可压缩 Navier–Stokes–Korteweg 方程中内爆奇点的形成

问题陈述
本文研究了具有密度依赖粘性和毛细系数的三维可压缩 Navier–Stokes–Korteweg (NSK) 系统在有限时间内的奇点形成问题。该系统定义为：
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \text{div}(2\mu(\rho)D(u)) + \nabla(\lambda(\rho)\text{div} u) + \text{div} K, \end{cases}$$
其中系数满足特定的幂律关系：
$$\mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha - 1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}.$$
之前的研究（Gu–Huang–Meng–Zhou [19] 以及 Huang–Lei–Zhou [21]）已经证明，当指数 $\alpha$ 处于特定适宜范围内（例如在 3D 中 $\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}, 1)$）时，对于远离真空的任意大初始数据，存在全局强解。本文探讨的核心问题是：当 $\alpha$ 落在该范围之外（特别是对于较小的正指数 $\alpha < 1/2$）时，解是全局存在的还是会在有限时间内发生爆破。在此机制下，有效体积粘度变为负值（$\lambda(\rho) + \frac{2}{3}\mu(\rho) < 0$），这一特征与非平衡声学介质和活性物质相关，对标准的全局存在性理论提出了挑战。

方法论
作者采用了一种基于构造平滑的、源自可压缩 Euler 方程的自相似内爆剖面的摄动方法。该方法通过以下几个关键步骤进行：

1. 有效速度重构： 为了处理三阶 Korteweg 张量，系统通过有效速度 $v = u + d\rho^{\alpha-2}\nabla\rho$（其中 $d = \nu - \sqrt{\nu^2-\varepsilon^2}$）进行重构。这种变换将系统转化为关于 $(\rho, v)$ 的抛物型-抛物型结构，揭示了质量方程中隐藏的退化扩散项 $d\Delta(\rho^\alpha)$。
2. 自相似缩放： 问题被转化到自相似坐标 $(\tau, y)$ 中，其中 $\tau = -\log(T-t)/\Lambda$ 且 $y = x/(T-t)^{1/\Lambda}$。在这些坐标下，只要缩放参数 $\Lambda$ 满足特定条件（$f_1 > 0$），粘性和毛细效应表现为指数衰减的摄动项（$e^{-f_1\tau}$）。
3. 摄动的稳定性分析： 证明的核心在于分析系统在平滑的、球对称的自相似 Euler 剖面 $(\bar{S}, \bar{V})$ 周围的非线性稳定性。通过推导摄动方程并构建 Bootstrap 论证来控制解。
4. 能量估计与谱分解：
   • 加权能量估计： 建立了高阶加权能量估计以控制摄动的衰减。作者利用特定的权重函数 $\phi$ 和参数层级来管理传输与耗散之间的相互作用。
   • Bootstrap 论证： 设计了一个严谨的 Bootstrap 方案来封闭低阶导数（特别是密度下界）和高阶导数的估计。一个关键的创新是同时利用传输结构和耗散结构来建立密度的下界，从而克服了密度方程除了传输结构外还具有耗散结构的困难。
   • 不稳定模控制： 将线性算子分解为稳定子空间和不稳定子空间。通过选择特定的初始数据系数来控制对应于有限维空间的成分。这种选择确保了解保持在一个导致爆破的有限余维流形上。

主要贡献与结果
本文证明了当 $\alpha$ 足够小（$\alpha < 1/2$）时，存在密度均匀远离真空的平滑初始数据，其演化出 3D NSK 系统的有限时间内爆奇点。

• 主定理 (Theorem 1.1)： 对于满足以下条件的参数：
  • (P1) $1 < \gamma < 1 + 2/\sqrt{3}$, $0 < \alpha < 1/2$, 且 $f_1 > 0$；或
  • (P2) $\gamma \in (1, 1 + 2/\sqrt{3}) \setminus \mathcal{J}$（其中 $\mathcal{J}$ 是一个可数例外集）且 $0 < \alpha < \alpha^*(\gamma)$，
    存在 $C^\infty$ 初始数据 $(\rho_0, u_0)$，使得对应的解 $(\rho, u)$ 在 $[0, T)$ 上存在，并在 $t=T$ 时产生内爆奇点。
• 奇点特征： 在爆破时刻 $T$：
  • 密度在原点处变为无穷大：$\lim_{t \to T^-} \rho(t, 0) = +\infty$。
  • 有效速度 $v$ 在原点的每个邻域内都是无界的。
  • 在经过重缩放后，解收敛于预设的自相似内爆剖面。
• 推论 1.1： 原系统产生了有限时间奇点，具体表现为密度爆破（$\lim_{t \to T^-} \rho(t, 0) = +\infty$）。然而，由于梯度项 $\nabla\rho$ 在原点附近也可能变得奇异，且无法排除在关系式 $v = u + d\rho^{\alpha-2}\nabla\rho$ 中出现抵消的可能性，因此不能直接从有效速度的爆破推导出物理速度 $u$ 的爆破。

意义
本文声称补充了现有的 NSK 系统全局存在性理论。虽然之前的研究展示了在特定 $\alpha$ 范围内全局适定性，但这项工作识别出了“小 $\alpha$”机制下一种截然不同的有限时间爆破机制。该机制的特征是负有效体积粘度结构，这一现象与非平衡物理相关。研究结果表明，即使在远离真空的平滑初始数据下，密度依赖的粘性、毛细作用与负体积粘度的相互作用，也会导致密度和有效速度场本身发生爆破（而非仅仅是其导数发生冲击波形成）的内爆奇点。该构造依赖于自相似 Euler 剖面的稳定性，将已知于 Euler 方程的内爆机制扩展到了具有粘性和毛细作用的 NSK 设定中。