2026年10月3日更新：ナビエ・ストークス方程式の「速度が発散してもエネルギーは有限」という主張は、どう理解すればよいのでしょうか。OpenAIが9月8日に公表した解決の主張を題材に、式の読み方、有限時間の発散、速度と体積の計算を順に解説します。発表の経緯と検証状況は解決発表の確認状況を整理した関連記事にまとめています。本記事では数学的な条件と量の違いを確かめます。
結論を先に述べると、外力があるという理由だけでCMIの問題から外れるとは言えません。一方、この発表を「どんな流れも計算できるようになった」「現実の水が無限の速さになると実証された」と読むのも誤りです。以下では、式の記号から発散の意味、問題文との対応へ順番に進みます。
確認した一次資料はOpenAIの公式発表、公開論文、公式のLeanリポジトリ、クレイ数学研究所（CMI）の問題文と賞の規則です。本記事は論文全166ページの証明を独立に検証したものではなく、Leanのビルドや証明検査も実行していません。
ナビエ・ストークス方程式とは何を表す式か
流体の小さな部分に運動方程式を当てはめる
ボールなら、位置と速度を追えば運動を考えられます。水や空気では、場所ごとに動く向きや速さが異なります。そこで、空間の各点に速度を割り当てた「速度場」を使います。川の中央は速く、岸の近くは遅い、という違いも一つの速度場で表せます。
ナビエ・ストークス方程式は、こうした流体の小さな部分について、加速度と力の釣り合いを記述する式です。ここで扱うのは、密度が一定で粘性のある非圧縮性流体です。空気を含め、あらゆる状況をこの条件で表せるわけではありません。圧縮や温度変化が重要なら、それに対応する別の条件や式も必要になります。
記号と各項を日本語で読む
密度で割った圧力をpと書くと、対象の式は次の形になります。uとfは三つの方向の成分を持つベクトルです。
\[ \begin{aligned}\frac{\partial\mathbf u}{\partial t}+(\mathbf u\cdot\nabla)\mathbf u &=-\nabla p+\nu\Delta\mathbf u+\mathbf f,\\ \nabla\cdot\mathbf u &=0.\end{aligned} \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| x、t | 位置、時刻 | m、s |
| u、|u| | 速度ベクトル、その大きさ（速さ） | m/s |
| P、ρ、p=P/ρ | 圧力、一定密度、密度で割った圧力 | Pa、kg/m³、m²/s² |
| ν | 動粘性係数 | m²/s、正の定数 |
| f | 単位質量あたりの外力 | m/s² |
| ∇、Δ | 空間微分、各方向の2階微分の和 | それぞれ長さの逆数、逆数の2乗を伴う |
| ∂u/∂t、(u・∇)u | 同じ場所での速度変化、移流による変化 | ともにm/s² |
| −∇p、νΔu | 圧力差、粘性による加速度 | ともにm/s² |
| ∇・u | 速度の発散（流出入の差） | 1/s、非圧縮なら0 |
加速度と非圧縮条件の読み方
1行目は「流体と一緒に移動して見た加速度＝圧力差、粘性、外力による加速度」です。2行目は、その小さな流体部分の体積が局所的に増えたり減ったりしない条件です。
移流項を、流れの速い場所へ運ばれる小さな浮きで考えましょう。川の流れが時間的に変化していなくても、浮きは移動先で速くなります。「同じ場所での変化」だけでは、この加速を数えられません。そのため左辺には二つの項があります。これは移流を理解する例であり、川そのものを今回の論文の解とみなしているわけではありません。
粘性があれば、なぜ必ず滑らかになるとは言えないのか
粘性には隣り合う流れの速度差をならす働きがあります。しかし移流項では、速度uが速度の変化にも掛かります。同じ空間分布の速度場を定数倍して2倍にすると、その空間微分も2倍になり、移流項は4倍になります。したがって、式全体を比例計算だけで扱えません。このような非線形性が、複数の流れを足して解を作ることを難しくします。
「ならす働きがある」ことと、「どんな細かい構造ができても、全時間にわたって変化を抑えられる」ことには隔たりがあります。3次元で必要な条件の下、解の滑らかさをどう保証し、あるいは破綻をどう示すかが問題です。νを0にしたオイラー方程式では粘性項がなくなるので、その証明をそのままナビエ・ストークス方程式の証明として扱うこともできません。
「有限時間で発散する」を図と数値で理解する
まず「非常に大きい」と「発散する」を区別しましょう。速さが100万でも、それだけなら有限です。有限時間の発散とは、ある有限の時刻へ近づく途中で、どんな有限の上限を決めても、それを超える値が現れることです。単に長い時間をかけて増え続ける場合とも違います。
1/(1−t)という説明用の関数
無次元の説明用関数を考えます。次のqは、時刻を表す無次元の量tが1に近づくと大きくなります。比較用のrは同じ範囲で有限です。ここでのtは実験の秒数ではなく、qも論文の流速ではありません。
\[ q(t)=\frac{1}{1-t},\qquad r(t)=1+t,\qquad 0\le t\lt1. \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| t | 説明用の時刻 | 無次元、0以上1未満 |
| q、r | 比較する関数の値 | 無次元。qは発散、rは2に近づく |
図の青線は上式のrです。こちらはt=1で2になり、有限時間で発散しません。橙線は上式のqで、t=1には定義されません。図は途中までしか描けないので、「線の上端が無限大」という意味ではありません。描画範囲の外でも値が増え続けることを、式の分母から確認します。
| t | 1−t | q=1/(1−t) |
|---|---|---|
| 0 | 1 | 1 |
| 0.9 | 1−0.9=0.1 | 1÷0.1=10 |
| 0.99 | 1−0.99=0.01 | 1÷0.01=100 |
| 0.999 | 1−0.999=0.001 | 1÷0.001=1,000 |
数表から任意の上限を超える議論へ
この表はすべて無次元の計算です。何行増やしても、有限個の数値だけで無限の振る舞いを証明したことにはなりません。任意の上限をBとすると、次の関係が成り立ちます。分母を正のまま1/Bより小さくできるので、その逆数はBを超えます。どれほど大きな上限にも同じ議論が使える点が、大量計算との違いです。\[ B>1,\quad 1-\frac1B\lt t\lt1\quad\Longrightarrow\quad q(t)>B. \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| B | 任意に決める上限 | 無次元、1より大きい |
| t、q(t) | 直前の関数の時刻と値 | 無次元、tは1より小さい |
発散する量を確かめる
流体の論文では、速度そのもの、速度の微分、渦度など、どの量に注目するかで主張が変わります。「特異点ができる」とだけ書かれていたら、その対象まで読む必要があります。また、通常の意味で微分できる滑らかな解と、積分の形で式を満たす弱い解も区別します。滑らかな解の破綻から、あらゆる意味の解がなくなるとは結論できません。
OpenAIの公開論文は何を主張しているか
公開論文「Finite time blowup for Navier–Stokes」の定理1.1は、正の粘性を固定した3次元の非圧縮性方程式について、静止した初期状態と滑らかな外力を使う構成を述べています。主張の要点は、運動エネルギーに対応する量に共通の有限な上限がある一方、tが1へ近づくと速度の最大の大きさが有界でなくなる、というものです。
「静止から始める」は、開始時の速度がすべての場所で0という意味です。ここで発散時刻を1とするのは数学上の尺度の選び方であり、現実の流体が1秒後に発散するという予測ではありません。外力まで0という意味ではありません。初期条件は開始時の状態、外力はその後も方程式に入る条件です。この二つを混同すると、何が流れを駆動しているかを読み違えます。
外力を好きに選べば簡単なのでは？
式を外力について並べ替えると、次のようになります。\[ \mathbf f=\frac{\partial\mathbf u}{\partial t}+(\mathbf u\cdot\nabla)\mathbf u+\nabla p-\nu\Delta\mathbf u. \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| u、t | 速度、時刻 | m/s、s |
| p、ν | 密度で割った圧力、動粘性係数 | m²/s²、m²/s |
| ∇、Δ | 1階の空間微分、2階微分の和 | 先の式と同じ演算 |
| f | 残りの項から求めた外力 | m/s²、全項で単位が一致 |
形式的には速度と圧力を選び、残りを外力と名付けられます。しかし速度を発散させた結果、外力まで特異になれば、滑らかな外力という要求を満たしません。「方程式の形が合う」ことと「外力の条件を満たす」ことは別です。
論文の第2節は、各項の特異な部分を相殺し、外力を滑らかに保つことが構成の課題だと説明しています。したがって読者が注目すべきなのは、迫力のある渦の図だけでなく、残る外力とその微分が要求された条件を本当に満たすかです。ここは本記事の例題から保証できる部分ではありません。
最大速度が増えてもエネルギーが有限な理由
ピークと全体の合計は別の量
質量と速さから物体の運動エネルギーを求める式を、流体の小さな部分に適用します。密度が一定なら、小さな体積の質量は密度と体積の積です。すべての部分を積分して足し合わせると、流体全体の運動エネルギーになります。
\[ \begin{aligned}K&=\frac12mv^2,\quad dm=\rho\,dV,\\ E(t)&=\frac{\rho}{2}\int_{\mathbb R^3}|\mathbf u(x,t)|^2\,dV.\end{aligned} \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| K、m、v | 物体の運動エネルギー、質量、速さ | J、kg、m/s |
| ρ、dm、dV | 一定密度、微小質量、体積要素 | kg/m³、kg、m³ |
| E(t)、u(x,t) | 流体全体のエネルギー、速度 | J、m/s |
| x、t、ℝ³ | 位置、時刻、3次元全空間 | m、s、積分範囲。定密度を仮定 |
積分記号∫は「空間全体で足し合わせる」操作です。最大速度は最も高いピークを測りますが、エネルギーは速さの2乗を体積付きで合計します。したがって、ピークが高くなると同時に、その領域が十分に狭くなる場合を考える必要があります。
速さ10倍・体積1,000分の1の計算
一定の速さを持つ領域で、速さを基準のA倍、体積を基準のAの3乗分の1にする説明用モデルを考えます。基準エネルギーと比べると、密度と基準値が約分され、倍率だけで比較できます。
\[ U=AU_0,\quad V=\frac{V_0}{A^3},\quad E_0=\frac12\rho U_0^2V_0, \]
\[ \begin{aligned}\frac{E}{E_0}&=\frac{\rho(AU_0)^2(V_0/A^3)}{\rho U_0^2V_0}\\ &=A^2A^{-3}=\frac1A.\end{aligned} \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| A | 速さの倍率 | 無次元、1以上 |
| U₀、U | 基準と比較対象の速さ | m/s、領域内で一定 |
| V₀、V | 基準と比較対象の体積 | m³、高速な空間領域を比較 |
| ρ | 共通の一定密度 | kg/m³ |
| E₀、E | 基準と比較対象のエネルギー | J、この領域だけの値 |
速さの2乗でA²倍になっても、体積がA³分の1になるため、差し引きで1/A倍です。単位を省略した数だけを掛けているのではなく、同じ基準に対する比を計算しています。密度は一定という前提です。
| 速さの倍率A | 体積の比V/V₀ | エネルギーの比E/E₀ |
|---|---|---|
| 1 | 1 | 1 |
| 10 | 0.001 | 100×0.001=0.1 |
| 100 | 0.000001 | 10,000×0.000001=0.01 |
身近なたとえなら、狭い場所の山の高さと、山全体の体積を比べる違いに似ています。ただしエネルギーで足すのは速さ自体ではなく、その2乗です。ピークが高いだけで、全体の合計が必ず大きくなるとは言えない点を押さえましょう。
単位を入れてエネルギーを検算する
単位も確かめましょう。密度を1,000 kg/m³、基準の速さを1 m/s、体積を1 L（0.001 m³）とします。速さを10倍、体積を1,000分の1にすると、比較する領域は10 m/s、1 mL（0.000001 m³）です。エネルギーは次のように0.5 Jから0.05 Jになります。
\[ \begin{aligned}E_0&=\frac12\times1000\times1^2\times0.001\ \mathrm J=0.5\ \mathrm J,\\ E&=\frac12\times1000\times10^2\times0.000001\ \mathrm J=0.05\ \mathrm J.\end{aligned} \]
| 記号 | 意味 | 単位・条件 |
|---|---|---|
| E₀、E | 基準と比較対象のエネルギー | J=kg·m²/s²、比は0.05÷0.5=0.1 |
| ρ | 代入した密度 | 1,000 kg/m³ |
| U₀、U | 代入した速さ | 1 m/s、10 m/s |
| V₀、V | 代入した体積 | 0.001 m³、0.000001 m³ |
この計算が証明していないこと
図と表は、量の違いを理解するための独自の説明例です。実際の論文の解や数値シミュレーションではなく、Aと体積の関係も説明のために選んでいます。この領域を同じ流体粒子の集まりと考え、非圧縮なのに押し縮めたと解釈しないでください。ここで比べているのは、速さが大きい空間領域の大きさです。
反対に、速さが10倍でも体積が10分の1にしか減らなければ、エネルギーは100×0.1で10倍に増えます。「領域が狭くなる」という言葉だけでは足りず、その割合が重要です。さらに非圧縮条件、運動方程式、滑らかな外力を同時に満たすかは、この倍率計算からは決まりません。
CMIのA〜Dと今回の主張を照合する
CMIの正式な問題文には、空間と証明の方向を分けたA〜Dの選択肢があります。まず外力と空間の条件を照合しましょう。次の表は概略で、滑らかさ・減衰・エネルギーなどの厳密な条件は原文で確認する必要があります。
| 選択肢 | 空間と外力 | 示す内容 |
|---|---|---|
| A | 3次元全空間／外力なし | 条件を満たす初期値すべてについて、要求された滑らかな解の存在。 |
| B | 3次元の周期的設定／外力なし | 周期条件の下で同様の存在。 |
| C | 3次元全空間／条件付きで外力あり | 要求された滑らかな解が存在しない初期値・外力の例。 |
| D | 3次元の周期的設定／条件付きで外力あり | 周期条件の下で破綻する例。 |
「すべて」と「ある例」、全空間と周期空間
「すべての条件で滑らか」と「ある条件で破綻」は、示す方向が違います。ただしAとCでは外力の条件も異なるので、CをそのままAの否定と呼んではいけません。OpenAIはC・Dに対応すると主張しており、それを外力なしの全問題を分類したという意味へ広げることはできません。
全空間は、壁のない3次元空間を考える設定です。周期的な設定は、一つの箱の向かい合う面をつなぎ、同じ流れが繰り返されるように考えます。実際の水槽の壁に流体が付着する条件とは違います。どちらも数学的に指定された舞台であり、家庭の配管や翼の周りをそのまま解いたことにはなりません。
AIとLeanの役割をどう評価するか
証明を探す作業と検査する作業
OpenAIは、多数のAIエージェントによる探索とLeanによる形式化を報告しています。探索は証明の候補を作る作業、形式化は定義や推論を形式言語で記述する作業です。その記述を検査器が受理するかという確認も、さらに区別します。自然言語の説明がもっともらしいことだけでは、数学的な正しさは決まりません。
公式リポジトリには、Lean 4による形式化、ビルド方法、独立した証明検査への案内があります。公開されていることは検証可能性の入口です。本記事では実行結果を確認していないため、「第三者として機械検証に成功した」とは述べません。
形式化を読む際は、検査が通るかに加え、最終定理にどんな仮定が含まれるか、形式化された「滑らか」「外力」「エネルギー」が原問題に対応するかを確かめます。狭い条件の定理を正しく証明しても、それが目的の問題を覆っているとは限りません。この確認はAI製か人間製かに関係なく必要です。
計算規模を成果そのものと取り違えない
エージェント数や所要時間は研究過程の情報ですが、大きな数字だけで証明の正否は判断できません。また企業の報告値と、独立した再現結果は区別します。AI一般の仕組みを知りたい方は、関連記事の機械学習では何を学習するのかも参照してください。予測や文章生成の仕組みと、個別の数学的主張の正しさは別の確認対象です。
未確認の点と続報を読む手順
2026年10月3日に確認したCMIの問題ページは「Active」と表示されています。この表示だけで証明が誤りだとも、単に更新が遅れているだけだとも断定できません。本記事で確認した資料から、数学界全体の受容や賞の認定を確認したとは言えない、というのが現在の範囲です。
発表、数学的評価、賞の手続きを分ける
CMIの賞の規則では、所定の媒体での出版、出版から少なくとも2年の経過、世界の数学界での一般的な受容が、提案された解決を検討する条件です。「2年待てば自動的に受賞」ではありません。また、OpenAI自身は賞金を請求するつもりはないと発表しています。請求の意向は証明の正否を判定する根拠にはなりません。
9月28日にはLeiとRenによる解説論文Part Iも公開されました。扱う構成部分と未完の部分を区別する必要があり、この公開だけで証明全体の第三者検証や賞の認定が完了したとは扱いません。続報は、まず公開論文の版や訂正を確認し、次に独立した検証がどの定理・補題を扱ったかを読み、最後にCMIの公式な扱いを確認する順で追えます。「検証済み」という一語ではなく、誰が、何を、どの方法で確かめたのかを見ると、確認範囲がはっきりします。
実用上、何がすぐ変わるのか
特定の条件で解が破綻するという存在の結果は、任意の条件の流れを高速・高精度に求める計算法とは違います。したがって、この発表だけで天気予報が急に正確になったり、設計シミュレーションが不要になったりするとは言えません。実用上の改善を評価するなら、計算法の誤差や計算時間、対象条件について別の証拠が必要です。
また連続体モデルの速度発散は、現実の水で無限速度を測定したという意味ではありません。数学では指定した方程式と仮定を厳密に調べ、現実への適用ではその仮定が有効な範囲も調べます。この二段階を保つことで、成果の価値を過大にも過小にも扱わずに読めます。
まとめ：式・ピーク・エネルギーを分けて読む
今回の発表を理解する鍵は、ナビエ・ストークス方程式の「外力」「正の粘性」「3次元」という条件と、発散する量を落とさないことです。最大速度とエネルギーは別の尺度なので、エネルギーが有界でも速度のピークが有界とは限りません。特に、速さ10倍と体積1,000分の1からエネルギー0.1倍が得られる計算が、その違いを示します。図と計算例はこの区別を理解する助けになりますが、論文の証明の代わりにはなりません。
OpenAIの解決の主張は、原問題との対応、公開された証明と形式化の検証、数学界での評価を追って判断する対象です。「AIだから正しい」「AIだから信用できない」のどちらにも寄せず、確認できた条件と未確認の範囲を明示して読むことが、このニュースから理解を深める方法です。
コメント