2026年10月11日日曜日

ナビエ−ストークス方程式(2)

ナビエ−ストークス方程式(1)からの続き


図:ナビエ−ストークス方程式(ncu_inoue_labから引用)

Q(koshix):
外力の有無,粘性の有無,圧縮性の有無,その他の条件分類の中で,流体力学の方程式の難易度はどうなっていて,ナビエ−ストークス方程式はどこに位置していて,オイラー方程式その他はどこに対応するのかが,混乱しているのでできるだけ,整理してわかりやすく説明してください。

A(Claude):

1.整理の枠組み:条件は「2値」ではなく「連続的な階段」

混乱の最大の原因は、外力・粘性・圧縮性を「アリ/ナシ」の2値だと思ってしまうことです。実際には外力の正則性クラスも散逸の強さも連続パラメータで、難易度はその上で滑らかに変化します。整理の軸は7つあります。

|軸――――――|パラメータ――――――|効果――――――――――――――――|
|① 空間次元  |2D / 3D       |最大の分岐。渦伸長項の有無     |
|② 散逸の強さ|$(-\Delta)^\alpha$ の $\alpha$(0=Euler, 1=NS)|大きいほど易しい|
|③ 外力   |$f$ の正則性クラス(超臨界→…→$C^\infty$急減少→$f\equiv0$)|滑らかにするほど難しい|
|④ 圧縮性  |非圧縮 / 圧縮性    |困難の種類が変わる         |
|⑤ 境界   |アリ / ナシ      |境界アリの方が易しい        |
|⑥ 解のクラス|古典解・強解 / Leray–Hopf弱解 / 一般の弱解|弱くするほど「易しい(が非一意)」|
|⑦ 対称性   |軸対称(swirlなし) など |対称性が強いほど易しい       |

そして、これら全部を貫くマスターキーが「スケーリング臨界性」です。


2.支配原理:超臨界性がすべてを決める

分数階散逸 $(-\Delta)^\alpha$ を持つ一般化NS方程式は、
$$u_\lambda(x,t) = \lambda^{2\alpha-1}u(\lambda x,\ \lambda^{2\alpha}t)$$
で不変です。唯一のア・プリオリ量であるエネルギー $\|u\|_{L^2}^2$ はこの変換で $\lambda^{2(2\alpha-1)-n}$ 倍になる。$n=3$ で指数がゼロになるのは
$$2\alpha - 1 = \tfrac{3}{2} \quad\Longrightarrow\quad \boxed{\alpha = \tfrac{5}{4}}$$
これが臨界点です。これより上なら勝ち、下なら(現在の道具では)負け。

|$\alpha$―――――――|名称―――――――――|状態―――――|
|$\alpha \ge 5/4$    |超散逸(hyperdissipative)|大域正則性 証明済(Ladyzhenskaya, Lions)|
|$5/4$ をわずかに下回る|対数的超臨界     |証明済(Tao 2009, Barbato–Morandin–Romito ほか)|
|$\alpha = 1$     |標準のNS方程式    |未解決(ミレニアム問題)|
|$0 < \alpha < 1$   |劣散逸(hypodissipative) |粗い外力付きなら爆発 証明済|
|$\alpha = 0$     |Euler方程式      |未解決(ただし爆発が強く予想)|

ここが決定的に重要な感覚です。標準のNS方程式は、臨界点からわずか $1/4$ だけ超臨界側に落ちている。ほんの少し粘性を強く(非局所的に)してやれば、90年来の難問は既に解けている。その $1/4$ の隙間に、ミレニアム賞100万ドルが埋まっているわけです。


3.軸① 次元:2Dはなぜ易しいのか

3次元の渦度方程式には渦伸長項 $\omega\cdot\nabla u$ があります。
$$\partial_t\omega + (u\cdot\nabla)\omega = \underbrace{(\omega\cdot\nabla)u}_{\text{3Dのみ}} + \nu\Delta\omega$$
2次元では渦度がスカラーになりこの項が消え、$\omega$ は単に流れに沿って運ばれるだけになる。したがって $\|\omega\|_{L^\infty}$ が保存され、爆発しようがない。これが2D Eulerの大域可解性(Yudovich 1963)の本質であり、2D NSも自動的に大域正則です。

渦伸長こそが、非圧縮流体の全困難の単一の源だと言ってよい。逆に言えば、渦伸長を持つ模型方程式はどれも同じ難しさを共有します。


4.軸② 粘性:Eulerの方が難しいのか易しいのか(重要な逆説)

ここが最も混乱しやすい点です。答えは「目標による」。

・爆発を証明したい → Eulerの方が易しい。粘性は特異点を潰そうと働くので、爆発を作るには粘性という敵と戦わねばならない。実際、爆発の証明はほぼすべてEuler側(ないし劣散逸側)で先に達成されています。
・大域正則性を証明したい → NSの方が望みがある。粘性のおかげでエネルギー散逸、Leray弱解の存在、Caffarelli–Kohn–Nirenbergの部分正則性といった構造が使える。Eulerにはこれらが一切ない。

だから「NSはEulerより難しい」も「EulerはNSより難しい」も、どちらも文脈次第で正しい。前回触れた「無限小の摩擦を入れただけで挙動が本質的に変わる」というのは、この非対称性のことです。

なおEuler方程式はクレイの賞金対象外ですが、Fefferman自身が公式記述の中で、これらの問題はEuler($\nu=0$)についても未解決かつ非常に重要だと明記しています。


5.軸③ 外力:ここが「階段」であることが最大の要点

外力を一切許さないのが本丸(A/B)ですが、外力を許す場合もその正則性で難易度が全く違う。これが今回のニュースの核心なので、階段として並べます。

|段|外力 $f$ のクラス―――――――|難易度――――――|状況――――――|
|0|超臨界($u$ から逆算した $f$)|ほぼ自明|爆発する $u$ を勝手に決めて $f:=\partial_t u+(u\cdot\nabla)u-\nu\Delta u+\nabla p$ と置くだけ。$f$ も一緒に爆発するので無意味|
|1|$L^1_t C^{1,\epsilon}_x \cap L^\infty_t L^2_x$ 程度 | 難しい|Córdoba–Martínez-Zoroa(3D Euler, 2023)、Córdoba–Martínez-Zoroa–Zheng(劣散逸NS, ARMA 2026)。well-posedクラス内での爆発=方程式の力学が生む爆発であり、外力が押し潰しているのではない|
|2|空間的に滑らか(時間方向は粗い)|さらに難しい | Córdoba–Martínez-Zoroa の IPM(2024)|
|3|時空一様に $C^\infty$ かつ全導関数が急減少(Clayの条件(5))|非常に難しい|Alpöge–Buckmaster–Coiculescu(IPM)→ 2026年8月のBoussinesq・3D Euler → OpenAIのNS主張|
|4|$f \equiv 0$           |本丸|3D Euler:OpenAIが主張(未検証)/3D NS(A/B):完全に未解決|

段0と段3の落差を掴むのが肝心です。「外力を許せば簡単」というのは段0の話であって、Clayが要求する段3は、$u$ が無限大に発散する最中も $f$ が最後まで完璧に滑らかで急減少していなければならない。各項が同時に巨大化しながら、外力として残る部分だけが精妙に相殺して滑らかに留まる——この設計が2023年から2026年までの3年をかけて詰められた部分です。


6.軸④ 圧縮性:難しさの「種類」が変わる(同じ物差しに乗らない)

ここも混同されがちですが、圧縮性は難易度の上下というより問題の性格を変える軸です。

|――――――――――|特異点の種類―――――――――|状態――――――――――|
|圧縮性Euler(3D) |衝撃波(微分の爆発、速度自体は有界)|有限時間で必ず形成されることが証明済(Sideris 1985、Christodoulouの精密理論)  |
|圧縮性Euler(1D系)|衝撃波             |エントロピー解の枠組みで存在・一意性理論が整備(Glimm, Bressan)|
|圧縮性NS     |真空・密度の退化        |大域弱解の存在は証明済(P.-L. Lions, Feireisl)/正則性は未解決/真空を含む古典解は爆発しうる(Xin 1998)|
|非圧縮Euler / NS  |渦の自己増幅(速度そのものが無限大)|未解決|

つまり圧縮性流体では「特異点が起きること」は驚きでも何でもなく、衝撃波として日常的に起きる。しかもそれは弱解(エントロピー解)の理論で数学的に乗り越え済みです。ミレニアム問題が問うている非圧縮の特異点は、衝撃波とは全く別種の、渦度が自己増幅して速度が発散するタイプであり、しかもそれを乗り越える弱解理論が存在しない(後述の非一意性のため)。

もう一つの構造的な違い:非圧縮では圧力が $-\Delta p = \nabla\cdot((u\cdot\nabla)u)$ で非局所的に決まり、情報が無限速度で伝わる。圧縮性では有限伝播速度で双曲型。解析の道具立てが根本的に異なります。


7.軸⑤ 境界:あると易しい

境界の役割はFeffermanが明確に述べています。境界のあるモデルで特異点が見つかっても、それは流体が境界との相互作用のおかげで特異点を作れることを示すにすぎない。境界なしなら、無茶をしているのは流体自身だ、と。

・Hou–Luo 2013:円筒の上下半分を逆回転させるとEulerが爆発する(数値+誤差評価)。「最初の本格的な特異点の主張」。
・Chen–Hou 2022:それを計算機援用で厳密証明(2D Boussinesq+境界あり3D軸対称Euler、滑らかな初期データ)。
・Elgindi 2021(Annals):境界なし $\mathbb{R}^3$ で爆発を証明。ただし初期データは $C^{1,\alpha}$ にとどまり、$C^\infty$ ではない。

つまり3D Eulerの爆発は「境界を許せば滑らかなデータで」「滑らかさを緩めれば境界なしで」既に証明済みで、両方同時(境界なし+滑らかなデータ+外力なし)だけが残っていた。そこがOpenAIが主張している地点です。


8.軸⑥ 解のクラス:弱くすると「易しくなるが無意味になる」

これも独立した軸です。

・Scheffer(1993)・Shnirelman(1997):2D Eulerに、時空コンパクト台を持つ弱解が存在する。静止した流体が外部刺激なしに勝手に動き出し、また止まる。Euler弱解の一意性は劇的に偽。
・Buckmaster–Vicol(2019):convex integrationにより3D NSの弱解も非一意(有限エネルギーだがLeray–Hopfクラス外)。
・Isett(2018)ほか:Onsager予想の解決。$C^{1/3}$ 以下ではエネルギー散逸する弱解が作れる。
・Vishik(2018):粗い外力を許すと2D Eulerですら非一意になる。

要するに、弱解の層では「何でも作れる」。だから「NSの弱解を作ってから正則性を示す」という古典的プログラムの後半が機能しない。ミレニアム問題が古典解($C^\infty$)を要求しているのはこのためです。


9.軸⑦ 対称性と、周辺モデル方程式の難易度階段

爆発現象の研究は、易しいモデルから順に登ってきました。

|モデル――――――――|性格――――――――|状態――――――――|
|Constantin–Lax–Majda  |1D、非局所項のみ |明示的に爆発(解ける)|
|De Gregorio / gCLM   |CLM+移流    |ほぼ解決|
|3D軸対称NS/Euler(swirlなし)|$\omega_\theta/r$ が輸送される|大域正則性 証明済|
|IPM(多孔質媒体、Darcy則)|2D、密度輸送  |外力付き爆発 証明済|
|2D Boussinesq     |3D軸対称Eulerの境界近傍と同型|境界ありで爆発 証明済(Chen–Hou)|
|SQG         |$\alpha=1/2$ 臨界散逸なら大域正則(Kiselev–Nazarov–Volberg, Caffarelli–Vasseur)。非粘性は未解決       |一部解決|
|3D Euler       |本命の一つ      |前述のとおり条件付きで証明済|
|3D NS        |最終目標       |A/B 未解決|
|平均化NS(Tao 2014)|エネルギー恒等式は保つ改変NS|爆発 証明済 ← 「エネルギー法では絶対に決着しない」ことの証明|

3D軸対称でswirlがない場合が解けているのに、swirlを入れた瞬間に未解決になる——この一点にも渦伸長の効き方が凝縮しています。


10.全体マップ:何がどこにあるか

最後に、すべてを一枚に並べます。上ほど易しい。

|設定―――――――――――――――|状態―――――――――――――――|
|2D NS/2D Euler(滑らかデータ、外力なし) |✅ 大域正則性 証明済|
|3D 超散逸NS $\alpha\ge5/4$(対数的にわずか下も可) |✅ 大域正則性 証明済|
|3D 軸対称 NS/Euler、swirlなし     |✅ 大域正則性 証明済|
|3D NS、初期値が小さい/短時間    |✅ 証明済|
|圧縮性Euler(3D)          |✅ 衝撃波形成 証明済(別種の特異点)|
|3D NS Leray弱解の存在        |✅ 証明済(1934)/一意性は ❓|
|3D NS 部分正則性(特異集合の1次元放物測度=0)|✅ CKN 1982。これ以上進むのが極めて困難|
|平均化NS/dyadicモデルの爆発     |✅ 証明済(エネルギー法の限界の証明)|
|3D Euler、境界あり、滑らかデータ、外力なし|✅ Chen–Hou 2022(計算機援用)|
|3D Euler、境界なし、$C^{1,\alpha}$ データ、外力なし|✅ Elgindi 2021|
|3D Euler/劣散逸NS、滑らかデータ、粗い外力|✅ Córdoba–Martínez-Zoroa 系列(2023–2026)|
|IPM・Boussinesq・3D Euler、時空一様に滑らかな外力|🆕 Buckmaster–Alpöge 2026(Lean検証済、査読前)|
|3D NS、滑らかな外力(Clay C/D)     |⚠️ OpenAI 2026 が主張(論文・Lean公開、独立検証待ち、Clay未受理) |
|3D Euler、境界なし、滑らかデータ、外力なし|⚠️ OpenAI が主張(未検証)|
|3D NS、外力なし(Clay A/B)       |❌ 完全に未解決。本丸|
|3D NS Leray–Hopf弱解の一意性      |❌ 未解決|
|圧縮性NSの正則性            |❌ 未解決|


11.まとめると

(1) 難易度の主軸は「超臨界性の度合い」であり、次元・散逸指数・外力の正則性はすべてその一つの物差しの上に並ぶ。
(2) 3D NSは臨界点から $1/4$ だけ外れている。この微小なギャップが90年の壁。
(3) 「外力アリ/ナシ」は2値ではなく5段の階段で、Clayが要求する段3(時空一様に滑らか+急減少)と段0(逆算した外力)は天と地ほど違う。2026年9月のニュースは、段3に到達したという主張。
(4) 粘性ナシ(Euler)は「爆発を示すには易しく、正則性を示すには難しい」という逆説的な位置にある。
(5) 圧縮性は難易度の軸ではなく、問題の種類を変える軸。圧縮性の特異点(衝撃波)は証明済みかつ理論的に処理済みで、ミレニアム問題とは別世界の話。
(6) 残る最難関は依然として 3D NS、外力なし、境界なし、滑らかデータ、古典解(A/B)。そしてそこへ至る最大の技術的障壁は、外力を消したときに粘性散逸が成長機構を食い潰してしまう点にある。

0 件のコメント:

コメントを投稿