Global existence and optimal decay for a three-dimensional penalized Navier--Stokes system with biharmonic damping

技術的要約：3次元における双調減衰を伴うペナルティ付きナビエ・ストークス系の全域的存在性と最適減衰

問題設定
本論文は、非圧縮ナビエ・ストークス方程式の超粘性およびペナルティ近似として機能する、3次元放物型系を調査している。このモデルは、3次元における任意の初期データに対する古典的なナビエ・ストークス方程式（1.1）のグローバルな正則性と一意性の未解決問題に対処するものである。提案された系（1.3）は、以下の3つのメカニズムを組み合わせている。

1. 古典的な粘性拡散 ($-\nu\Delta u$)。
2. 双調（超粘性）減衰 ($\beta\Delta^2 u$)：高次の散逸を提供する。
3. 発散ペナルティ ($-\frac{1}{\varepsilon}\nabla \text{div} u$)：$\varepsilon \to 0$ のときに近非圧縮性を強制することにより、厳密な非圧縮制約 ($\text{div} u = 0$) を緩和する。

決定的なことに、非線形対流項はテマム型の歪対称補正を用いて修正されている：$N(u) = (u \cdot \nabla)u + \frac{1}{2}u \text{div} u$。この補正は、ペナルティ導入によって生じる弱い圧縮性を補償し、圧縮形式では通常失われる輸送項の $L^2$ キャンセレーション特性を回復させる。

手法
解析は、エネルギー法、半群理論、およびフーリエ分裂技術の組み合わせに基づいており、特に $\varepsilon > 0$ に関して**一様（uniform）**な推定値を導出することに焦点を当てている。

• 弱解（定理 1.3）： 任意の $L^2$ 初期データに対する弱解の全域的存在性は、フリードリヒス型のガラーキン近似を通じて確立される。著者らは、ルレー射影（Leray projector）ではなく、カットオフ演算子 $J_n$ を用いて方程式をフーリエ球へ投影する。この区別は極めて重要である。ルレー射影を使用するとペナルティ項が消滅し、極限問題が変化してしまうからである。一様なア・プリオリ推定はエネルギー恒等式から導かれ、コンパクト性はオービンの・リオン・シモン（Aubin–Lions–Simon）補題によって得られる。
• 強解（定理 1.4）： $H^2(\mathbb{R}^3)$ における小さな初期データに対して、線形半群 $S_\varepsilon(t) = e^{tL_\varepsilon}$ に関連する緩和形式を用いた不動点引数を用いて、大域的な強解が証明される。線形作用素 $L_\varepsilon = \nu\Delta - \beta\Delta^2 + \frac{1}{\varepsilon}\nabla\text{div}$ は、平滑化特性を持つ半群を生成する。証明の鍵は、非線形項が粘性および双調の散逸によって吸収される、一様な $H^2$ エネルギー推定にある。
• 減衰推定（定理 1.6）： フーリエ分裂法（Schonbek）を用いて、最適な長時間減衰率が導出される。解析では低周波と高周波を分離する。低周波の挙動は、ドゥアメル（Duhamel）公式と非線形の特定の分解を通じて制御され、高周波は微分エネルギー不等式とガリオルディ・ニレンバーグ（Gagliardo–Nirenberg）補間を通じて管理される。

主要な貢献と結果

1. 弱解の全域的存在： 本論文は、任意の初期データ $u_0 \in L^2(\mathbb{R}^3)$ に対して弱解の全域的存在を証明している。解は、$L^\infty(0, T; L^2)$, $L^2(0, T; H^2)$, および $L^4(0, T; L^4)$ を含む特定の正則性クラスを満たす。
2. 一様な小ささを伴う大域的強解： $H^2(\mathbb{R}^3)$ における十分に小さな初期データに対して、系は一意な大域的強解を持つ。重要な貢献は、小ささの閾値 $\delta_0$ および関連するエネルギー境界が $\varepsilon$ に依存しないことである。この一様性は、結果が固定された単一の事例だけでなく、ペナルティ付き系の全家族に適用されることを保証する。
3. 最適減衰率： 初期データが $L^1(\mathbb{R}^3) \cap H^2(\mathbb{R}^3)$ である場合、著者らは以下の最適な減衰推定を導出している：
   $$\|u(t)\|_{L^2} \leq C(1+t)^{-3/4}$$
   $$\|\nabla^k u(t)\|_{L^2} \leq C(1+t)^{-3/4 - k/2}, \quad k=1, 2$$
   これらの速度は、3次元における古典的な熱方程式の漸近的減衰と一致する。解析は、双調項が高周波の散逸を支配する一方で、長時間の減衰は低周波の熱流挙動によって支配されることを確認している。
4. $\varepsilon$ に関する一様性： 解析の決定的な特徴は、すべての定数と推定値が $\varepsilon \to 0$ としても一様であることである。ペナルティ項 $\frac{1}{\varepsilon}\|\text{div} u\|^2$ は、保持されるか破棄されるかのいずれかであり、$1/\varepsilon$ を含む量によって上限が抑えられることはない、非負の散逸項として扱われる。

意義と主張
本論文は、非圧縮流のペナルティ近似に対して「安定かつ厳密な解析的基礎」を提供すると主張している。双調減衰とテマム補正を組み込むことで、本モデルは、無条件の $L^2$ エネルギー推定や高次の平滑化に関して、古典的なナビエ・ストークス方程式よりも優れたウェル・ポーズド（適切に設定された）結果を実現している。

著者らは、一様な推定値から、$\varepsilon \to 0$ となるにつれて、強解が非圧縮超粘性ナビエ・ストークス系を解く極限へと収束することを示唆している。しかし、極限系の特定および $\varepsilon$ に関する収束速度の証明については、今後の課題として明示的に残されている。さらに、非圧縮超粘性系（Lionsによる）については大きなデータに対する大域的正則性が既知であるが、著者らは、歪対称な圧縮構造が高次のエネルギー推定を修正するため、ペナルティ付き系において大きなデータに対する結果を直接拡張することはできないとして、本件に関する主張を限定している。

本研究は、吸収項を伴う減衰ナビエ・ストークス方程式（例）の研究とは異なり、より高次の減衰メカニズムと、圧力項を排除するペナルティ法に焦点を当てている。これにより、テスト関数に対する発散フリーの制約を取り除くことで、弱形式を簡素化している。