非圧縮と圧縮のNavier-Stokes方程式：密度、圧力、正則性

「非圧縮 vs. 圧縮」とは、実は密度に関する問いです。流体が流れるとき、その密度は場所ごとに同じままでしょうか、それとも変化するでしょうか？

注射器に水を入れてプランジャーを押してみてください。水は動きますが、圧縮はされません — 同じ空間により多くの水を押し込むことはできません。水はほぼ非圧縮です。今度は注射器に空気を入れて先端を塞いでください。プランジャーを押すと空気が圧縮されるのを感じることができます。同じ質量の空気がより小さな体積を占めます。これが圧縮流です。

非圧縮流では、密度 $\rho$ は流体全体で一定です。流体の微小な塊はすべて、移動しても体積を保ちます。圧縮流では、密度は変数になります — 場所ごと、時刻ごとに変化しうるのです。ジェットエンジン周りの空気、爆発中のガス、大規模な大気：これらは密度変化が物理を駆動する圧縮流を含みます。

この区別が重要なのは、ナビエ–ストークス方程式の見た目、予測するもの、解く難しさを変えるからです。

非圧縮の仮定は、密度場 $\rho$ を流れ領域全体で正の定数に設定します。物理的には、すべての物質体積要素が流れ写像のもとでその体積を保存することを意味します。圧縮の設定では、$\rho(x,t)$ を完全な未知関数に格上げし、独自の発展方程式で支配されます。

これらは同じ系の二つの表記法ではありません。未知数の数、拘束方程式の構造、そして圧力の性格が異なります。非圧縮系は $d+1$ 個のスカラー未知数（$u$ と $p$、$\mathbb{R}^d$ において）を持ちます。圧縮系は通常 $d+2$ 個の主要な場（速度、密度、内部エネルギーまたは温度）を持ち、圧力は状態方程式を通じて決定されます。

この区別は圧力のPDEの型も変えます：非圧縮の場合は楕円型（圧力は瞬時に調整される）ですが、圧縮の場合は系全体が混合双曲型–放物型であり、圧力の擾乱は有限速度で伝播します。これは技術的な細部ではなく、流体中を情報がどのように伝わるかを制御しています。

非圧縮ナビエ–ストークス方程式は、密度が一定の流体を記述します。クレイ・ミレニアム問題に登場するバージョンであり、本サイトが焦点を当てているバージョンです。

系は二つの部分からなります。運動量方程式：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

と、非圧縮性拘束条件：

$$\nabla \cdot u = 0$$

拘束条件 $\nabla \cdot u = 0$ は、速度場が発散なしであることを意味します：流体はどこでも溜まりも薄くもなりません。微小な領域に流入するものは、同じ速度で流出しなければなりません。この単一の条件が密度方程式全体を置き換えます — 密度が変化しないので、それを追跡する方程式は不要です。

この系における圧力は特殊な役割を果たします。圧力は（理想気体の法則のような）熱力学的法則によって決まるのではありません。むしろ、流れを発散なしに保つためにどこでも瞬時に調整されます。数学的には、$p$ は拘束条件から導出されるポアソン方程式を解きます。これは圧力変化が無限の速度で伝播することを意味します — 非圧縮流には「音速」がありません。

非圧縮ナビエ–ストークス系は二つの未知場を持ちます：速度場 $u$ と圧力場 $p$。この相対的な単純さは欺瞞的です — 非線形項 $(u \cdot \nabla)u$ が、3次元でこの系を途方もなく困難にしています。

$\mathbb{R}^3$ 上の、動粘性係数 $\nu > 0$ を持つ非圧縮ナビエ–ストークス系：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

発散なし条件 $\nabla \cdot u = 0$ は各点での拘束条件であり、発展方程式ではありません。局所的な体積保存を符号化しています：流れ写像 $\Phi_t$ は $\det(D\Phi_t) = 1$ をすべての $t$ で満たします。

運動量方程式に発散演算子を適用し、非圧縮性を用いると、圧力のポアソン方程式が得られます：

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

これは各固定時刻での $p$ に対する楕円型方程式です。圧力は独立した熱力学的変数ではなく、速度場によって大域的かつ瞬時的に決定されます。情報は圧力を通じて無限の速度で伝播します — 圧縮系との根本的な構造的差異です。

未知数は $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ と $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$ です。クレイの定式化（Fefferman, 2000）では、$u_0 \in C^\infty(\mathbb{R}^3)$ は発散なしであり、問いは $u$ がエネルギー有界で $C^\infty(\mathbb{R}^3 \times [0,\infty))$ に留まるかどうかです。

圧縮性ナビエ–ストークス方程式は、密度が変化する流れを支配します。より多くの未知数とより多くの方程式を持つ、より大きくより複雑な系です。

運動量方程式はまだありますが、今度は密度 $\rho$ が明示的に現れます：

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

拘束条件 $\nabla \cdot u = 0$ は消えます。その代わりに、密度がどう発展するかを追跡する連続の方程式が登場します：

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

これは質量が保存されることを述べています：流れが流体の塊を圧縮または膨張させることにより密度が変化します。

系にはさらにエネルギー方程式と状態方程式 — 圧力を密度と温度に結びつける $p = \rho R T$（理想気体の法則）のような熱力学的関係 — が必要です。圧力はもはや拘束条件の受動的な執行者ではありません。音響的な成分は音速で伝わり、さらに完全な粘性系には拡散効果も含まれます。

圧縮系は、高速空気力学、天体物理学的気体力学、燃焼、および密度変化が重要なあらゆる流れに不可欠です。しかしそれは、非圧縮方程式とは真に異なる数学的対象です — より多くの未知数、より多くの方程式、異なるPDE構造。

圧縮性ナビエ–ストークス系は、速度 $u(x,t)$、密度 $\rho(x,t)$、圧力 $p(x,t)$、比内部エネルギー $e(x,t)$（または温度 $\theta$）を結合します。保存形では：

連続の方程式： $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

運動量： $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

エネルギー： $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

ここで $E = e + \tfrac{1}{2}|u|^2$ は全比エネルギー、$\tau$ は粘性応力テンソル（ニュートン流体では $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$、体積粘性 $\lambda$）、$\kappa$ は熱伝導率です。

閉じるには状態方程式が必要であり、例えば理想気体では $p = (\gamma - 1)\rho e$（$\gamma$ は断熱指数）です。

圧縮系の非粘性・双曲型部分では、音響擾乱は有限速度で伝播し、特性的な音速は $c = \sqrt{\partial p / \partial \rho |_s}$ です。完全な粘性ナビエ–ストークス–フーリエ系は双曲型と放物型の混合であり、粘性拡散と熱拡散は数学的には無限伝播を導入します。圧縮方程式は衝撃波、膨張波、接触不連続を許容しますが、これらは非圧縮流にはアナログがありません。

実用的には、圧縮流と非圧縮流の境界線はマッハ数です：

$$\text{Ma} = \frac{|u|}{c}$$

ここで $|u|$ は流速、$c$ は音速です。マッハ数は、圧力擾乱が伝わる速度に対して流れがどれだけ速く動いているかを測定します。

経験則：$\text{Ma} < 0.3$ のとき、密度変化は約5%未満であり、非圧縮方程式は優れた近似です。部屋の中の空気、管の中の水、建物周りの風 — これらはすべて低マッハ流です。空気は技術的には圧縮性ですが、非圧縮であるかのように振る舞います。

$\text{Ma} \approx 0.3$ を超えると、圧縮効果が大きくなります。$\text{Ma} \approx 1$ 付近の遷音速領域では、局所的な超音速ポケットが現れ、衝撃波が形成され始めます。戦闘機、ロケットのノズル、再突入する宇宙船がこの領域に存在します。

これはバイナリスイッチではなくスペクトルです。低マッハ領域は、非圧縮ナビエ–ストークス方程式が物理モデルとして適用される場所です。日常の流体流のほとんど — そしてクレイ・ミレニアム問題 — はこの領域の真っ只中にあります。

マッハ数 $\text{Ma} = |u|/c$ は圧縮性の重要度をパラメータ化します。ここで $c = \sqrt{\partial p / \partial \rho |_s}$ は等エントロピー音速です。形式的には、非圧縮方程式は圧縮系の低マッハ極限として出現します。

$\text{Ma}^2$ のべきでの漸近展開（Klainerman & Majda, 1981, 1982; Schochet, 1986参照）は、適切な初期データで $\text{Ma} \to 0$ のとき、圧縮解が非圧縮解に収束することを示しています。標準的な低マッハ無次元化で準備された初期データでは、圧力は概ね空間的に一様な熱力学的背景と、極限で非圧縮拘束条件を強制するより小さな動的補正として書かれることが多いです。

これは特異極限です：音速 $c \to \infty$ であり、圧縮系の双曲型の性格が非圧縮流の楕円型圧力方程式に退化します。音響モードは無限に速くなり、渦度ダイナミクスから切り離されます。

$\text{Ma} < 0.3$ という領域は工学的な経験則です。相対的な密度変動 $\delta\rho / \rho \sim \text{Ma}^2 / 2$ がこの範囲で $\sim$5% 以下に留まるという経験的観察を反映しています。数学的な正当化は低マッハ極限の収束定理であり、これは準備された初期データと適合する境界条件を要求します。

クレイ・ミレニアム問題は正確な問いを発しています：$\mathbb{R}^3$ 上の滑らかで発散なしの初期速度が与えられたとき、非圧縮ナビエ–ストークス系は常にすべての時刻で存在する滑らかな解を生成するか？

なぜ特に非圧縮なのか？ 三つの理由があります。

第一に、それだけでも十分に困難です。非圧縮3次元方程式は、1934年のルレイの基礎的研究以来、大域正則性の証明に抵抗してきました。可変密度、熱力学、衝撃波の複雑さを加えれば、問題ははるかに困難になり、より扱いやすくはなりません。

第二に、困難は純粋な流体力学です。非圧縮系は、中核的な数学的挑戦 — 非線形移流 $(u \cdot \nabla)u$ と粘性散逸 $\nu \Delta u$ の相互作用 — を、熱力学的・音響的合併症なしに分離しています。非圧縮の場合は、圧縮性、熱力学、衝撃波を追加する前に正則性の問題を分離するのに十分困難です。

第三に、物理がきれいです。非圧縮方程式は、最も一般的な日常の流れをモデル化しています。滑らかなデータから特異点を生成しうるかどうかを理解することは、古典流体力学の数学的整合性に関する根本的な問いです。

圧縮系にも独自の深い未解決問題 — 大データでの大域解の存在、衝撃波の形成と相互作用 — がありますが、これらは異なる構造を持つ異なる問題です。クレイ賞は、3次元ナビエ–ストークスに対するクレイの定式化した特定の正則性問題であるため、非圧縮の場合を対象としています。

クレイの公式定式化（Fefferman, 2000）は $\mathbb{R}^3$ 上の非圧縮系を指定しています：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

$u_0 \in C^\infty(\mathbb{R}^3)$ は発散なしで、問いは $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ であって $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ がすべての $t \geq 0$ で有界であるような解が存在するかどうかです。

非圧縮系の選択は数学的に動機づけられています。主要な未解決問題 — ルレイ–ホップ弱解（グローバルに存在するが一意性も滑らかさも保証されない）と古典的な滑らかな解（局所的に存在するが爆発しうる）のギャップ — は、非圧縮3次元方程式に固有のものです。2次元では、滑らかな非圧縮ナビエ–ストークス解の大域正則性は知られています。3次元の場合は未解決のままです。

圧縮系は質的に異なる困難を導入します：衝撃波形成（滑らかなデータからのオイラー方程式でも起きる）、真空状態（$\rho \to 0$）、渦度と音響モードの結合。これらは重要な未解決問題ですが、非圧縮正則性の問題とは構造的に異なります。

非圧縮問題は、エネルギー超臨界的な非線形性と粘性散逸の間の競合を分離します。自然なエネルギー評価は $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ を与え、これは3次元でスケーリング臨界空間 $L^\infty_t \dot{H}^{1/2}_x$ に半微分分だけ足りません。このギャップを閉じること — あるいは閉じられないことを証明すること — がミレニアム問題の核心です。

完全な非圧縮系と各項の意味については、ナビエ–ストークス方程式とは何か？をお読みください。

方程式がニュートンの第二法則と応力テンソルからどう構築されるかについては、ナビエ–ストークス方程式の導出をご覧ください。

粘性を完全に取り除いてオイラー方程式を得るとどうなるかについては、オイラー方程式 vs. ナビエ–ストークス方程式をご覧ください。

ミレニアム問題の正確な記述と「存在と滑らかさ」が実際に何を意味するかについては、ナビエ–ストークスの存在と滑らかさの問題に進んでください。

ここからの推奨パス：

• ナビエ–ストークス方程式とは何か？ — 完全な非圧縮系と項ごとの解析
• ナビエ–ストークス方程式の導出 — 連続体力学の公理からPDE系へ
• オイラー方程式 vs. ナビエ–ストークス方程式 — 非粘性極限 $\nu \to 0$ と正則性における粘性の役割
• ナビエ–ストークスの存在と滑らかさの問題 — クレイの正確な定式化、ルレイ–ホップ理論、弱解と強解のギャップ