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

技術要約：3次元圧縮性Navier–Stokes–Korteweg方程式における内破特異点の形成

問題設定
本論文は、密度依存型の粘性と毛管現象係数を持つ3次元圧縮性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$ が特定の適切な範囲（例：3次元において $\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$）となり、これは非平衡音響媒体やアクティブマターに関連する特徴であり、標準的な大域存在理論を困難にするものである。

手法
著者らは、圧縮性オイラー方程式から導出された滑らかな自己相似内破プロファイルに基づく摂動的アプローチを採用している。その手法は以下の主要なステップを通じて進行する：

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. 摂動の安定性解析： 証明の核心は、滑らかな球対称自己相似オイラープロファイル $(\bar{S}, \bar{V})$ の周りでの系の非線形安定性の解析である。摂動方程式が導かれ、解を制御するためのブートストラップ引数が構築される。
4. エネルギー評価とスペクトル分解：
   • 重み付きエネルギー評価： 摂動の減衰を制御するために、高次の重み付きエネルギー評価が確立される。著者らは、輸送と散逸の相互作用を管理するために、特定の重み関数 $\phi$ とパラメータの階層を利用する。
   • ブートストラップ引数： 低次の導関数（特に密度の下限）と高次の導関数の両方の評価を閉じるために、厳密なブートストラップ・スキームが設計されている。重要な革新は、密度方程式が輸送構造に加えて散逸構造を持つという困難を克服するために、輸送構造と散逸構造の両方を用いて密度の下限を確立した点にある。
   • 不安定モードの制御： 線形作用素は安定部分空間と不安定部分空間に分解される。不安定成分は有限次元の空間に対応しており、特定の初期データ係数を選択することによって制御される。この選択により、解が爆発へと至る有限余次元多様体に留まることが保証される。

主な貢献と結果
本論文は、$\alpha$ が十分に小さい（$\alpha < 1/2$）とき、密度が真空から一様に離れた滑らかな初期データが存在し、それが有限時間の内破特異点へと進化することを証明している。

• 主定理 (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)$
    対応する解 $(\rho, u)$ が $[0, T)$ 上に存在し、時刻 $t=T$ で内破特異点を形成するような $C^\infty$ 初期データ $(\rho_0, u_0)$ が存在する。
• 特異性の特性： 爆発時刻 $T$ において：
  • 密度は原点で無限大となる：$\lim_{t \to T^-} \rho(t, 0) = +\infty$。
  • 有効速度 $v$ は、原点のあらゆる近傍において非有界となる。
  • 解は、再スケーリング後、規定された自己相似内破プロファイルに収束する。
• 系図 1.1 (Corollary 1.1): 元の系は有限時間の特異性を発達させる。具体的には、密度爆発（$\lim_{t \to T^-} \rho(t, 0) = +\infty$）によって特徴付けられる。ただし、有効速度の爆発から物理的速度 $u$ の爆発を直接結論付けることはできない。なぜなら、原点付近で勾配項 $\nabla\rho$ も特異になる可能性があり、関係式 $v = u + d\rho^{\alpha-2}\nabla\rho$ における相殺の可能性を排除できないためである。

意義
本論文は、NSK系に関する既存の大域存在理論を補完することを目的としている。これまでの結果は、特定の $\alpha$ の範囲において大域的なウェルポーズドネスを示してきたが、本研究は「小さな $\alpha$」の領域における明確な有限時間爆発メカニズムを特定している。この領域は、負の有効体積粘性構造によって特徴付けられ、非平衡物理学において関連のある現象である。これらの結果は、滑らかな初期データであっても、密度依存の粘性、毛管現象、および負の体積粘性の相互作用によって、単なる導関数の爆発（衝撃波形成）ではなく、密度や有効速度場自体が爆発する内破特異点が生じ得ることを示している。この構成は、自己相似オイラープロファイルの安定性に依拠しており、オイラー方程式で知られている内破メカニズムを、粘性と毛管現象を持つNSKの設定へと拡張したものである。