ナビエ–ストークス方程式の弱解・強解・滑らかな解
ミレニアム賞が求めるのは滑らかな解です。存在が証明できるのは弱解だけです。そのギャップこそが問題のすべてです。
弱解とは何か?
状況を整理しましょう。ミレニアム賞は、ナビエ–ストークス方程式が常に滑らかな解を持つことを証明すれば100万ドルを贈呈します。滑らかとは無限回微分可能、折れ目なし、爆発なし、予想外の振る舞いなし、という意味です。しかし、ほぼ1世紀に及ぶ挑戦を経て、誰かが証明できた最良の存在結果は、ナビエ–ストークス方程式の弱解の存在だけです。
では弱解とは何でしょうか?近似ではありません。「だいたい正しい」ものでもありません。弱解は方程式の厳密な解ですが、ルールが緩和されています。速度場が直接微分できるほど滑らかであることを要求する代わりに、方程式にテスト関数を掛けて積分し、部分積分によってすべての微分をテスト関数に移します。その結果がすべてのテスト関数に対して成り立てば、弱解が得られたことになります。
こう考えてみてください。古典解は試験のすべての問題を途中式をきちんと示しながら解く学生です。弱解は途中式を見せることはできませんが、あなたが出しうるどんな問題に対しても最終的な答えが正しいことが証明できる学生です。微分する過程は見えませんが、積分した結果は必ず辻褄が合います。
なぜわざわざそんなことを? 古典解には方程式が荒すぎる場合があるからです。流体が速度の急激な変化を起こし、通常の意味での微分が存在しなくなる領域を発生させることがあります。弱解は古典解が行き詰まるところでも先に進めます。それは安全ネットなのです。
落とし穴は:弱解は一意でないかもしれないことです。同じ初期条件から複数の弱解が存在しうるのです。物理学は流体が一つの特定の振る舞いをすべきだと言っているのに、複数が出てくるのは問題です。そして弱解は滑らかでないかもしれません。滑らかさこそミレニアム賞が要求するものであり、誰も証明できていないものです。
ナビエ–ストークス方程式の弱解は、各点でのPDEを分布的定式化に置き換えたものです。$\mathbb{R}^3 \times (0,T)$ 上の非圧縮系を考えます:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$$
ベクトル場 $u \in L^2_{\mathrm{loc}}(\mathbb{R}^3 \times [0,T))$ が弱解であるとは、すべての滑らかでコンパクト台を持つ発散なしテスト関数 $\varphi$ に対して以下が成り立つことです:
$$\int_0^T \!\int_{\mathbb{R}^3} \bigl[-u \cdot \partial_t \varphi - (u \otimes u) : \nabla \varphi + \nu \, u \cdot \Delta \varphi \bigr] \, dx \, dt = \int_{\mathbb{R}^3} u_0 \cdot \varphi(x,0) \, dx.$$
$\varphi$ が発散なしであるため、圧力はこの定式化から消えます。部分積分によって、すべての微分は $u$ から $\varphi$ に移されています。要点は:$u$ は古典的に微分可能である必要がなく、これらの積分が収束するのに十分な可積分性があればよいということです。
弱定式化は近似ではありません。PDEを各点で満たす古典解も、すべてのテスト関数に対して弱定式化を満たします(逆方向に部分積分すればよい)。逆は成り立ちません:弱解がPDEを各点で満たすほど滑らかである必要はありません。
これが重要なのは以下の理由です:(1) $L^2$ 初期データに対して弱解は大域的に存在する(Leray 1934)が、古典解は3Dでは時間的に局所的にしか存在が知られていない;(2) 弱解は一般に一意であることが知られていない;(3) 弱解が常に滑らかであるかどうかという問題は、適切な初期データに対するクレイ・ミレニアム問題と同値である。
Lerayと最初の存在証明(1934年)
1934年、Jean Lerayはこの分野を今なお規定する偉業を成し遂げました。Acta Mathematicaに掲載された73ページの一本の論文「Sur le mouvement d'un liquide visqueux emplissant l'espace」において、3次元の発散なし初期データに対するナビエ–ストークス方程式の弱解がすべての時刻にわたって存在することを証明したのです。これは3D方程式に対する最初の大域的存在結果であり、90年以上が経過した今日でも最も強い無条件の存在定理です。
彼の手法は美しく具体的でした。非線形性をモリファイ(滑らかに近似)して方程式を正則化し、正則化された方程式を解き(最悪の非線形相互作用が抑えられているため容易)、次に平滑化パラメータをゼロに近づけて極限を取り、コンパクト性の議論によって極限が存在し弱定式化を満たすことを保証します。
しかしLerayが証明しなかったことがあります。一意性です。彼の手法は少なくとも一つの弱解を生成しますが、同じ初期データから始まる他の弱解が存在するかもしれず、それを排除できませんでした。滑らかさも証明していません。彼の解は有限エネルギーを持ち、エネルギー不等式を満たします(エネルギーは粘性を通じて散逸しうるが自然に増加はしない)。それだけです。それ以上のことはありません。
Leray自身は特異点が形成されると考えていました。彼はその姿を描きました:流体が特定のスケーリング構造を持ちながら一点に向かって崩壊し、ますます速く、エネルギーを無限小の領域に集中させる自己相似的爆発です。1996年、Nečas、Růžička、Šverákは、この正確な自己相似的シナリオが起こりえないことを証明しました。潜在的な爆発の形に関するLerayの推測は間違っていたのです。爆発がどんな形であれ実際に起こるのか?誰にもわかりません。
1951年、Eberhard HopfはLerayの構成を有界領域に拡張し、得られたクラスはLeray-Hopf弱解として知られるようになりました:エネルギー不等式を満たす弱解です。これが標準的な概念です。研究者が追加の修飾なしに「弱解」と言うとき、ほぼ常にこれを意味しています。
もう一点。Leray-Hopf弱解の中にも、適切な弱解(suitable weak solutions)と呼ばれる厳密により小さなクラスがあり、大域的なエネルギー不等式だけでなく局所的なエネルギー不等式を満たします。Caffarelli、Kohn、Nirenbergは1982年の有名な部分的正則性の結果を、まさにこの小さなクラスに対して証明しました。二つを混同しないでください:CKNが適用されるのは適切な弱解であり、すべてのLeray-Hopf解ではありません。
Lerayの1934年の論文は以下を確立しました:任意の発散なし $u_0 \in L^2(\mathbb{R}^3)$ に対して、$\mathbb{R}^3 \times (0,\infty)$ 上のナビエ–ストークス方程式の弱解 $u$ が少なくとも一つ存在し、以下を満たす:
- $u \in L^\infty(0,\infty; L^2(\mathbb{R}^3)) \cap L^2(0,\infty; \dot{H}^1(\mathbb{R}^3))$
- エネルギー不等式:ほとんどすべての $t > 0$ に対して $\|u(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u(s)\|_{L^2}^2 \, ds \leq \|u_0\|_{L^2}^2$
構成はモリファイによって進められます。非線形性 $(u \cdot \nabla)u$ を $(u_\varepsilon \cdot \nabla)u$($u_\varepsilon = J_\varepsilon * u$ は空間的モリファイ)に置き換えます。正則化された系は大域的に滑らかな解を持ちます(モリファイが最悪の非線形相互作用を除去するため)。Lerayは正則化された解に対する一様なエネルギー評価を得て、弱収束する部分列を抽出しました。その極限は弱定式化とエネルギー不等式を満たします。
Lerayは一意性を確立しませんでした。コンパクト性の議論は少なくとも一つの集積点の存在を与えますが、異なる部分列は異なる極限に収束する可能性があります。3DにおけるLeray-Hopf弱解の一意性は今日まで未解決です。
Hopf (1951)はGalerkin近似(有限次元部分空間への射影)を用いて、Dirichlet境界条件を持つ有界領域 $\Omega \subset \mathbb{R}^3$ に構成を適応させました。得られたクラス——エネルギー不等式を満たす弱解——は両方の名前を持ちます:Leray-Hopf弱解。
Caffarelli-Kohn-Nirenbergの定理(1982)は厳密により小さなクラスに関するものです:適切な弱解は、以下の形の局所エネルギー不等式を追加的に満たします
$$\partial_t \left(\frac{|u|^2}{2}\right) + \nabla \cdot \left(\left(\frac{|u|^2}{2} + p\right)u\right) + \nu |\nabla u|^2 \leq \nu \Delta \left(\frac{|u|^2}{2}\right)$$
分布の意味で成り立ちます。CKNは、任意の適切な弱解に対して、時空における特異点集合の1次元放物型Hausdorff測度がゼロであることを証明しました。これは特異点が存在するとしても極めてまばらであることを意味します(時空の曲線を埋め尽くすことはできない)。しかしこの定理は特異点が実際に発生するかどうかについては何も述べておらず、適切な弱解にのみ適用され、すべてのLeray-Hopf解には適用されません。
Leray自身は $u(x,t) = (T-t)^{-1/2} U(x / (T-t)^{1/2})$ の形の自己相似的爆発の存在を予想しました。Nečas、Růžička、Šverák (1996)は $L^3(\mathbb{R}^3)$ に属する解に対してそのような自己相似的爆発が存在しないことを証明し、Tsai (2009)はこれを漸近的に自己相似的な爆発の排除にまで拡張しました。潜在的な特異点の形は、もし存在するとしても、未知のままです。
強解と正則性
弱解は大域的に存在します。しかし一意でないかもしれず、滑らかでないかもしれません。もっと良い結果を得られるでしょうか?
はい、ただし局所的にだけです。ナビエ–ストークス方程式の強解は、PDEが分布的な意味だけでなく古典的(各点的)な意味で成り立つのに十分な正則性を持つ解です。3Dの滑らかな初期データに対して、強解は短い時間だけ存在します。どれだけ短いかはデータに依存します:初期速度が大きく荒いほど、保証される存在区間は短くなります。そしてこれらの強解が最終的に爆発しないことは誰にも証明できていません。
美しい部分的結果があります。1962年、James Serrinは、Leray-Hopf弱解がたまたま $L^p_t L^q_x$($2/p + 3/q \leq 1$、$q > 3$)に属していれば、自動的に滑らかであることを証明しました。そして弱解-強解の一意性により、その時間区間においてその初期条件を持つ唯一のLeray-Hopf弱解となります。これは十分条件であって必要条件ではありません。Leray-Hopf弱解がSerrin条件を満たせば、それを「昇格」させることができます:最初からずっと正則だったのであり、同じデータを持つ他のLeray-Hopf弱解はそれと異なることはできません。条件を満たさなければ、何も結論できません。
Serrin条件は条件付き正則性の結果です。「もし解が荒すぎなければ、完全に良い振る舞いをする」と教えてくれます。困難の全体は「もし」の部分を証明することにあります。
2次元では、エネルギー評価が十分に強く、すべてのLeray-Hopf弱解が自動的にSerrin型の条件を満たします。だからこそ2Dは解決済みなのです。3Dでは、エネルギー評価が必要なものより半微分だけ不足しており、そのギャップを埋めることがすべてです。
他の条件付き結果も存在します。オイラー方程式に対しては、Beale-Kato-Majda基準により、渦度が時間に関して可積分である限り滑らかな解は継続します。ナビエ–ストークスの類似物が存在し、これを精密化しています。Escauriaza-Seregin-Šverák基準(2003)はこれを端点まで押し進めます:$u \in L^\infty_t L^3_x$ ならば爆発は起きません。これらのそれぞれは、大域的正則性の問題を一つのアプリオリ評価に変換します。これらの評価のいずれか一つを無条件に証明すればミレニアム問題は解決されます。誰も成功していません。人々が試みてきたさまざまな証明戦略の概観は、専用のページがあります。
ナビエ–ストークス方程式の強解とは、PDEが(ほとんどいたるところで)各点的に成り立ち、非線形項 $(u \cdot \nabla)u$ が単なる分布ではなく関数として定義できるだけの正則性を持つ解です。典型的には $u \in L^\infty(0,T; H^1(\mathbb{R}^3)) \cap L^2(0,T; H^2(\mathbb{R}^3))$ を意味します。関連するが異なる枠組みとしてFujita-Kato (1964)があり、臨界空間における局所的な温和解を構成します。
初期データ $u_0 \in H^s(\mathbb{R}^3)$($s \geq 1/2$)に対する強解の局所存在は、温和解の形でのFujita-Katoの不動点論法から従います:
$$u(t) = e^{\nu t \Delta} u_0 - \int_0^t e^{\nu(t-s)\Delta} \mathbb{P} \nabla \cdot (u \otimes u)(s) \, ds,$$
ここで $\mathbb{P}$ は発散なし場へのLeray射影です。この積分方程式は適切な関数空間における縮小写像原理によって一意な局所解を持ちます。臨界空間($L^3$、$\dot{H}^{1/2}$、$BMO^{-1}$)における小さなデータに対しては、解は大域的です。
問題は、大きなデータの強解がすべての時刻にわたって持続するかどうかです。鍵となる条件付き結果はSerrin (1962)のものです:Leray-Hopf弱解が $u \in L^p(0,T; L^q(\mathbb{R}^3))$ を満たし
$$\frac{2}{p} + \frac{3}{q} \leq 1, \qquad q > 3,$$
ならば、$u$ は $(0,T] \times \mathbb{R}^3$ 上で滑らかであり、与えられた初期データを持つ一意のLeray-Hopf弱解です。これらはProdi-Serrin条件と呼ばれます(Prodi 1959が関連する結果を確立しました)。
端点 $q = 3$($p = \infty$)はEscauriaza、Seregin、Šverák (2003)によって解決されました:$u \in L^\infty(0,T; L^3(\mathbb{R}^3))$ ならば、$u$ は時刻 $T$ で爆発しません。これはLeray-Hopfの枠組みにおける最も強い既知の継続基準です。
2次元では、エネルギー有界性 $u \in L^\infty_t L^2_x \cap L^2_t H^1_x$ とLadyzhenskaya不等式 $\|f\|_{L^4}^2 \leq C\|f\|_{L^2}\|\nabla f\|_{L^2}$(2D特有)を組み合わせることで $u \in L^4_t L^4_x$ が得られ、これはSerrin条件 $2/4 + 2/4 = 1$(2D版で $2/p + 2/q \leq 1$)を満たします。3Dでは、エネルギー有界性はSobolev埋め込みにより $u \in L^{10/3}_t L^{10/3}_x$ を与え、$2/(10/3) + 3/(10/3) = 3/2 > 1$ となります。Serrin条件は超臨界性ギャップに正確に対応するマージンで成り立ちません。この構造的障害の詳細については、なぜナビエ–ストークスは難しいのかを参照してください。
滑らかな解とミレニアム問題
滑らかな解は最高基準です。空間と時間の両方で無限回微分可能、特異点なし、角なし、不連続性なし。何度微分しても常にきちんとした答えが得られます。
クレイ・ミレニアム賞問題は、2000年にCharles Feffermanによって定式化され、正確な問いを発しています:$\mathbb{R}^3$ 上の滑らかで発散なしの、適切な減衰を持つ初期データが与えられたとき、ナビエ–ストークス方程式は有界なエネルギーを持ちすべての時刻に存在する滑らかな解を常に持つか?あるいは解が爆発する初期データを見つけることができるか?
どちらの答えも100万ドルの価値があります。大域的な滑らかさを証明するか、爆発を構成するか。どちらも歴史的な成果です。
現状はこうです。3Dで滑らかな解が常に大域的に存在することは誰も証明していません。爆発を構成した人もいません。その中間で行き詰まっており、Lerayの1934年の論文以来ずっと行き詰まっています。90年以上にわたって答えがわからないのです。
局所的な描像は問題ありません。滑らかな初期データに対して、方程式はある時間区間で滑らかな解を生成します。流体は動き始め、数学は機能し、すべてがきれいです。問題はその後に何が起こるかです。解は永遠に滑らかなままか、それとも最終的に速度(またはその微分)が無限大になる特異点を発生させるか?
もし滑らかなままであれば、その滑らかな解は自動的に弱解でもあります(すべての古典解は弱定式化を満たすため)。弱解-強解の一意性原理により、Leray-Hopf弱解の中でも自動的に一意です:強解が存在すれば、他の弱解はそれと競合できません。つまり大域的な滑らかさを証明すれば、階層全体が崩壊します。弱、強、滑らか:すべて同じ対象です。だからこそこの問題はこれほど魅力的であり、これほど難しいのです。存在が証明できるもの(弱解)と求めるもの(滑らかな解)のギャップこそが、問題のすべてです。
クレイ・ミレニアム問題(Fefferman 2000)は次のように問います:任意の発散なし $u_0 \in C^\infty(\mathbb{R}^3)$ で、すべての $\alpha, K$ に対して $|\partial^\alpha u_0(x)| \leq C_{\alpha K}(1 + |x|)^{-K}$ を満たすものに対して、ナビエ–ストークス方程式を満たす $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ と $p \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ で、すべての $t \geq 0$ に対して $\int_{\mathbb{R}^3} |u(x,t)|^2 \, dx < C$ となるものは存在するか?
もう一つの選択肢(同様に賞に値する):上記のクラスに属する $u_0$ で、そのような滑らかな解が存在しないものを見つけよ。
知られていること:
- 局所存在:$u_0 \in H^s(\mathbb{R}^3)$($s \geq 1/2$)に対して、$\|u_0\|_{H^s}$ に依存するある $T^* > 0$ について $[0,T^*)$ 上に一意な局所温和解/強解が存在し、すべての正の時刻 $t > 0$ で滑らかです。$T^* < \infty$ ならば $t \to T^*$ のとき $\|u(t)\|_{H^s} \to \infty$ です。
- 小さなデータの大域的存在:$\|u_0\|_{\dot{H}^{1/2}}$(あるいは $\|u_0\|_{L^3}$、$\|u_0\|_{BMO^{-1}}$)が普遍定数より小さければ、解は大域的で滑らかです(Fujita-Kato 1964、Koch-Tataru 2001)。
- 弱解-強解の一意性:$[0,T]$ 上に強解が存在すれば、同じ初期データを持つすべてのLeray-Hopf弱解は $[0,T]$ 上でそれと一致します。これはSerrin (1962)によって確立され、後続の研究で精密化されました。正則性を証明することがLeray-Hopfクラス内での一意性も決定することを意味します。
階層は上方に崩壊します:滑らか $\Rightarrow$ 強 $\Rightarrow$ 弱、そして弱解-強解の一意性は、滑らかな解が存在するならそれが一意のLeray-Hopf弱解であることを意味します。したがってミレニアム問題は次と同値です:滑らかな初期データを持つすべてのLeray-Hopf弱解はそれ自身が滑らかか?「弱解が存在する」(Leray 1934)と「滑らかな解が存在する」(未解決)のギャップこそが、まさに賞の問いです。
なぜこの区別が重要なのか
弱解が存在し流体を記述するのなら、なぜ滑らかさを気にするのでしょうか?
理由は三つあります。
第一に、一意性。物理学は一つの答えを要求します。流体の初期状態を与えられたら、次に何が起こるかを正確に、厳密に、曖昧さなく言えるべきであり、可能性のメニューを差し出して肩をすくめるべきではありません。弱解は一意性を保証しません。同じ出発点から複数の弱解が現れる可能性があり、物理的な流体がどれに従うかはわかりません。一意性が崩れれば、方程式はそのレベルの正則性での決定論的予測力を失います。
第二に、数値計算の信頼性。地球上のすべてのCFDシミュレーションは、これらの方程式を近似的に解いています。一般的な滑らかさと一意性の理論がなければ、格子の細分化がすべての3D設定で単一の物理的に妥当な解に収束することを保証する定理はありません。あるマシンでは一つの弱解に収束し、別のマシンでは別の弱解に収束する可能性があります。
第三に、極端な物理。もし特異点が形成されうるなら、それは現実の流体が最も小さなスケールで何をするかについて何かを教えてくれます:方程式が無限の速度や無限の圧力勾配を予測し、連続体モデル自体が破綻する合図です。
これは技術的な細事ではありません。すべてを貫く断層線です。一方の側に:存在(Leray、1934年、完了)。もう一方の側に:正則性(未解決、100万ドル、そして生きているすべての解析学者の集合的な苛立ち)。なぜ渡るのがこれほど難しいのか? 3Dのエネルギー評価が必要なものからちょうど半微分だけ不足しており、90年の努力と数百の論文、キャリア全体をかけても、そのギャップは埋まっていないからです。
現在追求されているすべての証明戦略は、この分水嶺を橋渡しする試みです。弱解が滑らかであることを証明する。あるいはそうでないことを証明する。二語か三語か。どちらの答えも教科書を書き換えます。
3D非圧縮ナビエ–ストークス方程式の解の階層は以下の通りです:
$$\text{滑らか} \subset \text{強} \subset \text{Leray-Hopf弱} \subset \text{分布的弱}$$
各レベルで知られていること:
| クラス | 大域的存在 | 一意性 | 正則性 |
|---|---|---|---|
| 分布的弱 | あり(エネルギー制御なし) | なし | 非常に粗い可能性がある |
| Leray-Hopf弱 | あり(Leray 1934、Hopf 1951) | 未解決 | $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ |
| 適切な弱 | あり | 未解決 | CKN:特異点集合の $\mathcal{H}^1 = 0$ |
| 強/温和 | 未解決(小さなデータでは可) | あり(存在区間内で) | 存在すれば $C^\infty$ |
| 滑らか($C^\infty$) | 未解決(= ミレニアム問題) | あり(弱解-強解の一意性) | 定義上 $C^\infty$ |
ミレニアム問題は、大域的Leray-Hopf弱解の存在と大域的滑らかな正則性のギャップに位置しています。核心的困難:エネルギー不等式は $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ を与え、これは臨界スケーリング $\dot{H}^{1/2}$ より半微分低いです。この超臨界性ギャップを橋渡しすることは、大域的正則性を証明することと同値です。
弱解-強解の一意性は、滑らかな解が大域的に存在すれば階層が崩壊することを意味します。より正確には:与えられた $u_0 \in C^\infty$ に対して $[0,T]$ 上に滑らかな解 $u$ が存在すれば、同じデータを持つすべてのLeray-Hopf弱解は $[0,T]$ 上で $u$ に等しいです。したがって大域的滑らかさ $\Rightarrow$ Leray-Hopfクラス内での大域的一意性。
逆に、Leray-Hopf弱解の非一意性は、滑らかな解が大域的に持続しないことを意味します(滑らかな解は一意性を強制するため)。凸積分に関する最近の研究(オイラー方程式に対するDe Lellis-Székelyhidiに基づき、Buckmaster-Vicol 2019がLeray-Hopf正則性以下のナビエ–ストークスの非一意弱解を構成するように拡張)は、分布的弱解が高度に非一意でありうることを示しています。この非一意性がLeray-Hopfクラスに及ぶかどうかは、ミレニアム問題に直接的な含意を持つ重要な未解決問題です。
このギャップに挑む証明戦略の現状、およびそれがかくも困難である構造的理由については、なぜナビエ–ストークスは難しいのかを参照してください。