ナビエ–ストークス問題はなぜ難しいか
立ちはだかる核心的な数学的障害
非線形性の罠
物理学のほとんどの方程式は線形です — 出力は入力に比例します。原因を倍にすれば、効果も倍になります。線形方程式はよく理解されており、(比較的)解きやすいです。
ナビエ–ストークス方程式は非線形です。流体の速度が自身の変化率に影響する — 流体が自分自身を押し動かします。群衆の中のすべての人の動きが他の全員の動きに依存するとき、群衆がどこに向かうかを予測するようなものです。
この自己相互作用項 $(u \cdot \nabla)u$ が方程式を非常に困難にしています。小さな擾乱が大きなものに増幅されるフィードバックループを作り出し、流体の乱流がこれほど複雑になる理由です(詳細は部分問題をご覧ください)。
対流的非線形性 $(u \cdot \nabla)u$ が根本的な障害です。渦度の定式化 $\omega = \nabla \times u$ では、方程式は
$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$
となります。渦伸長項 $(\omega \cdot \nabla)u$ は符号を持たず — 渦度を際限なく増幅できます。2次元では、$\omega$ が流れに垂直なスカラーであるためこの項は消失し、2次元の大域正則性が知られている理由です(ラディジェンスカヤ、1969年)。3次元では、渦伸長が有限時間爆発の主要な候補メカニズムです。
決定的に重要なのは、非線形性が $u$ について2次であることです:$H^1$ エネルギー評価は $\|\nabla u\|_{L^2}$ を与えますが、$(u \cdot \nabla)u$ を $L^2$ で制御するには $u \in L^\infty$ または少なくとも $u \in L^3$ が必要です — エネルギークラスが提供しない情報です。
超臨界性:スケーリングギャップ
重要な洞察:ナビエ–ストークス方程式はスケーリング対称性を持っています。解を拡大する(適切な量ですべてをより小さく、より速くする)と、別の有効な解が得られます。
これが問題です。制御できる唯一の量 — 流体の全エネルギー — が「間違ったスケール」にあるからです。全体像については教えてくれますが、爆発が形成されるであろう非常に小さなスケールで何が起きているかについては教えてくれません。
都市の総電力使用量を監視して1つの火花を検出するようなものです。測定は実在し有用ですが、心配している事象を捉えるには十分に細かくありません。それが研究者が橋渡しする必要のあるギャップです。
自然なスケーリング $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$ の下で、臨界ソボレフ空間は $\dot{H}^{1/2}(\mathbb{R}^3)$(同値的に $L^3$)です。量は以下のように分類されます:
- 劣臨界:再スケーリング下でノルムが縮小 — 大スケールの振る舞いを捉えるが、小スケールの集中を見逃す — 例:$\|u\|_{L^2}$
- 臨界:スケール不変 — 例:$\|u\|_{L^3}$、$\|u\|_{\dot{H}^{1/2}}$
- 超臨界:スケーリング下で縮小(小スケールを優先)
エネルギー不等式は $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ の制御を与えます。両成分とも劣臨界です:
$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$
これは、エネルギー評価が小スケールの制御を提供しないことを意味します — 非線形性は原理的に細かいスケールで散逸を圧倒できます。劣臨界エネルギークラスから臨界ノルムへの橋渡しが中心的な困難です。
乱流とエネルギーカスケード
川を見たことがある人なら誰でも、流体の運動がカオス的 — 乱流的 — になることを知っています。大きな渦が小さな渦に砕け、それらがさらに小さな渦に砕け、粘性がようやく滑らかにする微視的スケールまで続きます。
このエネルギーカスケード(1941年にコルモゴロフが記述)は、ナビエ–ストークス方程式によって美しく捉えられます。しかし、それは危険も示唆しています:エネルギーが粘性が散逸できるよりも速く、ますます小さな領域に集中したらどうなるか? それが爆発です。
これが実際に起こりうるのか — あるいは粘性が最終的に常に勝つのか — がまさに未解決の問題です。レイノルズ数からこの小スケールの描像へのより物理的な橋渡しについては、レイノルズ数、乱流、そして小スケールが重要な理由をご覧ください。
コルモゴロフのK41理論は、慣性領域 $k_f \ll k \ll k_\eta$ でエネルギースペクトル $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ を予測します。ここで $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ はコルモゴロフ散逸波数です。この領域ではエネルギーフラックスはスケール間で一定です。
正則性問題は、このカスケードが退化しうるかを問います:散逸スケール $k_\eta^{-1}$ が有限時間でゼロに縮小しうるか? これには $\|\nabla u\|_{L^2} \to \infty$(散逸率の爆発)が必要ですが、全エネルギーは有限のままです。
散逸異常予想(オンサーガー、1949年)は、消失粘性極限 $\nu \to 0$ でエネルギー散逸が持続することを示唆します — オイラーの弱解がエネルギーを散逸できます。これは $1/3$ 未満のヘルダー指数に対して確認されています(イゼット、2018年;バックマスターら、2018年)が、ナビエ–ストークスの正則性との関連は不明のままです。このセクションのレジームレベルの直感については、レイノルズ数、乱流、そして小スケールが重要な理由をご覧ください。
圧力の問題
ナビエ–ストークス方程式において、圧力は奇妙な役割を果たします。独立変数ではなく — 拘束条件(流体は非圧縮であり、押しつぶせない)を通じて速度によって完全に決定されます。
これにより圧力は非局所的になります:ある点での速度の変化がどこでも瞬時に圧力に影響します。流体のあらゆる部分が見えないバネで他のすべての部分と結ばれているようなものです。
この非局所性により方程式の解析がはるかに困難になります。1つの点で何が起きているかを、流体全体を同時に考慮せずに研究することはできません。
非圧縮拘束 $\nabla \cdot u = 0$ はポアソン方程式を通じて圧力を決定します
$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j)$$
したがって $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$ であり、リース変換 — 特異積分作用素 — を含みます。圧力は速度の非局所的な汎関数であり、この非局所性が点ごとの、または空間的に局所的な評価に対する主要な障害です。
特に、標準的な最大値原理の議論は失敗します:粘性項 $\nu \Delta u$ は散逸的であるにもかかわらず、圧力勾配 $-\nabla p$ は遠方の領域からエネルギーを集中させることができます。カファレリ–コーン–ニーレンバーグ理論は放物型シリンダー上の局所エネルギー不等式でこれに対処しますが、これらから点ごとの正則性を抽出することが困難なステップのままです。
なぜ3次元が特別なのか
2次元では、ナビエ–ストークス問題は解決されています — 滑らかな解はいつでも存在します。(ラディジェンスカヤが1969年に証明しました。)
では3次元で何がうまくいかないのか? 鍵となる違いは渦伸長です。2次元では、渦は回転や合体はできますが、伸びることはできません。3次元では、流体が渦管をますます細く引き延ばし、すべてのエネルギーを無限に細いフィラメントに集中させる可能性があります。
この伸長過程が有限時間で無限大に暴走しうるのか — あるいは粘性が常に介入してそれを止めるのか — が百万ドルの問題です。
2次元と3次元の二分法は鋭いです:
- 2次元:渦度 $\omega$ は $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$ を満たすスカラーです。最大値原理は $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$ を与え、BKMは大域正則性を含意します。渦伸長項 $(\omega \cdot \nabla)u$ は恒等的にゼロです。
- 3次元:渦度 $\omega \in \mathbb{R}^3$ は $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$ を満たします。伸長項 $(\omega \cdot \nabla)u$ は $|\omega|$ を超線形的に増幅できます(ビオ–サバールを通じて形式的に $\sim |\omega|^2$)、そして最大値原理は利用できません。
エンストロフィー $\|\omega\|_{L^2}^2$ は
$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$
を満たしますが、$\|\nabla u\|_{L^\infty}$ を制御するには $\omega \in L^\infty$ が必要であり、既存のいかなる手法も破ることのできない循環依存関係を生み出しています。
さらに探索する
この記事は未解決問題の一部です。
これらの障害により、数学者たちは問題を部分問題に分解し、各問題に対して専門的なアプローチを開発してきました。
粘性項がなぜ役に立つが十分ではないかという文脈については、オイラー方程式 vs. ナビエ–ストークス方程式をご覧ください。
この記事は未解決問題の一部です。
正則性問題の取り組み可能な構成要素への分解については、部分問題をご覧ください。これらの障害に対処するために開発された解析ツールについては、アプローチをご覧ください。粘性項がなぜ役に立つがギャップを閉じないかという文脈については、オイラー方程式 vs. ナビエ–ストークス方程式をご覧ください。