Navier–Stokes 方程式 III：外力つき方程式の有限時間爆発という主張（OpenAI, 2026）を読む

0. この記事の要点

静かな水に、外から力を加えて動かし始めるとします。力の加え方をうまく選ぶと、水全体の運動エネルギーは有限のままなのに、ある決まった時刻に、ある一点で流れの速さがいくらでも大きくなってしまう。2026 年に公表された原稿は、そうした力と流れが存在すると主張しました。この記事は、その主張が正確には何を言っているのか、何を言っていないのか、そして確かめられたことと確かめられていないことを分けて読みます。

• 2026 年 9 月 8 日、OpenAI は原稿 “Finite time blowup for Navier–Stokes”（166 ページ）を公表しました。主張は「任意の粘性 $\nu\gt0$ に対し、初期値 $0$ と、時空でコンパクト台（ある有界な範囲の外では $0$）の滑らかな外力（流体を外から押す力）で、運動エネルギーを有界に保ったまま時刻 $1$ で速度が非有界になる解がある」です。この記事の執筆時点（2026-10-08）で、独立した査読を経た結果ではありません。以下では「主張」として扱い、正否の判定はしません。
• 主張が扱うのは Clay の問題文の命題 (C)(D)（外力を許す「解の破綻」）です。外力なし（$f=0$）の命題 (A)(B) の成否、すなわち「滑らかな初期値だけから特異点が生じるか」は、この主張では扱われていません。 主張が正しくても (A)(B) は未解決のままです。
• 構成の骨格は次のとおりです。任意の発散なし（非圧縮）の速度場は、運動量の残差（方程式の左辺と右辺のずれ）を外力と定めれば方程式の解になる。問題は、速度が爆発するのに残差が全導関数とともに特異時刻を越えて滑らかに延びるようにできるか、であり、原稿は (1) 半径方向が軸方向より速く縮む軸対称の渦、(2) その外縁の環状領域に残る特異な残差、(3) 背景のせん断で成長する 2 族の振動パルスの運動量流束（運動量が単位時間にどれだけ運ばれるか）による残差の相殺、(4) 逐次補正と切り落とし、の順でこれを行うと述べています。
• この記事で証明するのは、主張の正否に依らない検算（残差と外力、粘性の変換、周期化）、主張を仮定に明記した導出（(D) と、条件つきの (C)）、原稿のスケールが既知の定理（爆発の速さの下界、LPS の条件）と矛盾しないことの確認です。
• Lean による形式化が公開されています。定理文が Fefferman の条件と一致するかは、この記事で 1 行ずつ照合しました。形式証明そのもの（ビルドと検査）と、原稿の記述と形式証明の対応は確認していません。
• 一次資料で確認したことと、報道だけに基づくことを分けて書きます。確認日はすべて 2026-10-08 です。

この記事で使う事前知識

• 外力: 流体を外から押し続ける力 $f$。方程式の右辺に入る。外力なしの場合が Clay の問題の命題 (A)(B)、外力を許す場合が (C)(D)。→ 命題 (A)〜(D)
• Leray–Hopf の弱解: 方程式を積分の形で満たし、エネルギーの不等式を満たす、滑らかとは限らない解。→ Leray–Hopf の弱解
• コンパクト: 台がコンパクトとは、有界な範囲の外で $0$ になること。距離空間では、どの点列も収束する部分列をもつ性質。→ コンパクト
• 発散: ベクトル場の湧き出しの強さ。発散が $0$ の場は非圧縮の流れを表す。→ 発散
• 熱方程式: 熱が広がる様子を表す線形の方程式で、渦の粘性による拡散にも同じ形が現れる。→ 熱方程式
• 正則性の判定条件: 解が滑らかなままでいられるための条件（LPS の条件など）。→ Navier–Stokes 方程式 II

1. 動機と、確認した資料

第 9 章 で見たとおり、Clay の問題は 4 つの命題のどれか 1 つの証明を求めています。(C)(D) は外力を許すので、外力の設計に自由度があります。第 10 章 の最後で見た Córdoba と Martínez-Zoroa の系列は、この自由度を使って特異性を作る方針で、物差しは「外力をどこまで滑らかにできるか、散逸をどこまで強くできるか」でした。2026 年 9 月の主張は、外力が $C^\infty$・散逸が通常の $\nu\Delta$ という、その物差しの端に当たります。

この記事の目的は、主張の内容を、読める範囲で正確に述べることです。この記事が依拠する資料と、読んだ範囲を先に示します。

資料	所在	読んだ範囲
原稿 “Finite time blowup for Navier–Stokes”（著者表記 OpenAI、166 ページ、PDF の作成日時 2026-09-08 19:06 UTC）	cdn.openai.com/pdf/…/navier-stokes.pdf	要旨、§1〜§3（1〜24 ページ）、§10（116〜126 ページ）を通読。§4〜§9 は Theorem 4.6、Proposition 5.5・7.5・9.6・9.9、Lemma 7.4 の文と各節の冒頭。§4〜§9 と付録 A〜C の証明（約 130 ページ）は検証していない
Lean リポジトリ openai/NavierStokesAndEuler	github.com/openai/NavierStokesAndEuler	README.md、formalization.yaml、lean-toolchain、ComparatorChallenges/README.md・NavierStokes.json・NavierStokes.lean（定理文。全文）、NavierStokes/ComparatorSolution.lean、NavierStokes/PeriodicPaperTheorem.lean の冒頭。証明本体は読んでおらず、ビルドも検査も実行していない
OpenAI の発表ページ	openai.com/index/navier-stokes-solution/	読めなかった（HTTP 403）。このページの内容は、報道を通じてしか確認できていない
原稿 “Finite time blowup for the Euler equation”（同時に公表）	cdn.openai.com	読んでいない（この記事では扱わない）
Clay の問題文（Fefferman）	claymath.org の PDF	全文（第 9 章で条件を書き出した）
Clay の問題ページ・規則のページ	claymath.org	問題ページの表示と、規則の要点
T. Buckmaster の声明（PDF の作成日時 2026-09-08 02:38 UTC）	cims.nyu.edu/~tristanb/statement.pdf	全文（4 ページ）
Lean リポジトリ tristanbuckmaster/fluid_lean	github.com/tristanbuckmaster/fluid_lean	euler-blowup/README.md・boussinesq-blowup/README.md の定理文、コミットの日時
L. Alpöge・T. Buckmaster・M. P. Coiculescu の原稿（多孔質媒質方程式）	arXiv:2609.16470	要旨のページ。Boussinesq と Euler の原稿（著者のページに置かれている PDF）は読んでいない

本稿の執筆時点で確認できたのは上の範囲までです。 原稿の §4〜§9 と付録について以下に書くことは、原稿自身の §2・§3 の説明と各命題の文に基づく要約で、その証明が正しいことの確認ではありません。

2. 準備：記号

原稿の記号に合わせます。特異時刻は $t=1$、特異点は空間の原点で、$\tau=1-t$ は特異時刻までの残り時間です。

定義 2.1（運動量の残差）↗

滑らかな速度場 $u$ と圧力 $p$、粘性 $\nu\gt0$ に対し

$$R_\nu(u,p)=\partial_tu+(u\cdot\nabla)u-\nu\Delta u+\nabla p$$

を運動量の残差という。$\nu=1$ のとき $R(u,p)$ と書く（原稿の式 (3.1)）。

定義 2.2（円柱座標・相似座標・指数）↗

• 円柱座標 $x=(r\cos\theta,r\sin\theta,z)$ と、単位ベクトル $e_r,e_\theta,e_z$。速度の成分を $u_r,u_\theta,u_z$ と書く。$\theta$ によらない場を軸対称という。$u_\theta$ を旋回（swirl）、$u_\theta,u_z$ を接線成分という。
• 指数 $0\lt h\lt\frac1{100}$ を固定し、$A=\frac12+h$、$D=\frac12-h$ とおく。
• 相似座標（原稿の式 (3.2)）: $\tau\gt0$ と $z$ に対し、$q\gt0$ と $-1\lt\eta\lt1$ を

$$\tau=q(1-\eta^2),\qquad z=q^D\eta$$

で定め、$X=\dfrac{r^2}{2q}$ とおく。$q$ を集中のスケールという。

• 速さ $U$、長さ $\ell$、粘性 $\nu$ に対し Reynolds 数を $U\ell/\nu$ とする。粘性による拡散の時間 $\ell^2/\nu$ と、流れが距離 $\ell$ を進む時間 $\ell/U$ の比である。

相似座標が定まることと、$q$ の大きさを確かめておきます。爆発点の近くでは、半径方向と軸方向で長さの縮む速さが違うため、時間と高さ $(\tau,z)$ のままでは見通しが悪くなります。そこで両方をまとめて 1 つのスケール $q$ で測り直す座標を導入するのですが、そのためには $(\tau,z)$ から $q$ が一意に決まることと、$q$ が $\tau$ と $|z|$ で上下から抑えられることが要ります。次の補題はその確認で、証明は中間値の定理を使うだけです。

補題 2.3（集中のスケール q の一意性と大きさ）✓

$0\lt h\lt\frac12$ とする。$\tau\gt0$、$z\in\mathbb{R}$ に対し、定義 2.2 の $(q,\eta)$ はただ 1 組存在し、

$$\tau+|z|^{1/D}\le2q,\qquad q\le2\tau+2^{1/(2D)}|z|^{1/D}$$

が成り立つ。とくに $z=0$ では $q=\tau$、$\tau\downarrow0$ では $q\to|z|^{1/D}$ である。

証明（補題 2.3）

$\eta=zq^{-D}$ を代入すると、条件は $g(q):=q-z^2q^{2h}=\tau$（$1-2D=2h$）かつ $|\eta|\lt1$、すなわち $q\gt|z|^{1/D}$ です。$q_0=|z|^{1/D}$ では $g(q_0)=q_0-q_0^{2D+2h}=0$ です。$q\gt q_0$ で $g'(q)=1-2hz^2q^{2h-1}=1-2h\eta^2\ge1-2h\gt0$ なので $g$ は狭義単調増加で、$g(q)\ge(1-2h)(q-q_0)\to\infty$ です。中間値の定理により $g(q)=\tau$ となる $q\gt q_0$ がただ 1 つあります。

$\tau=q(1-\eta^2)\le q$、$|z|^{1/D}=q|\eta|^{1/D}\le q$ から第 1 の不等式が出ます。第 2 の不等式は場合を分けます。$1-\eta^2\ge\frac12$ なら $q\le2\tau$ です。そうでなければ $|\eta|\gt2^{-1/2}$ で、$|z|=q^D|\eta|\gt q^D2^{-1/2}$ から $q\lt(2^{1/2}|z|)^{1/D}=2^{1/(2D)}|z|^{1/D}$ です。

∎

$z=0$ の平面では $q=\tau$、$X=r^2/(2\tau)$ で、半径方向の長さは $\tau^{1/2}$ です。$|\eta|$ が $1$ より小さい一定の範囲では $|z|\le q^D\approx\tau^{1/2-h}$ なので、軸方向の長さは $\tau^{1/2-h}$ で、半径方向よりゆっくり縮みます。$|\eta|\to1$ は「特異時刻に、原点以外の高さ $z\ne0$ に近づく」ことに当たります。

3. 主張の正確な文

注意 3.1（主張（OpenAI, 2026。独立の査読は未了））

出典: OpenAI, “Finite time blowup for Navier–Stokes”（原稿。上の表の PDF。確認日 2026-10-08）の Theorem 1.1、Proposition 10.1、Corollary 10.6。

主張 1（原稿 Theorem 1.1）。 任意の $\nu\gt0$ に対し、外力 $f\in C_c^\infty(\mathbb{R}^3\times(0,\infty);\mathbb{R}^3)$、コンパクト集合 $K\subset\mathbb{R}^3$、$\mathbb{R}^3\times[0,1)$ 上の滑らかな速度場 $u$ と圧力 $p$ が存在して、次を満たす。

• (a) $\partial_tu+(u\cdot\nabla)u-\nu\Delta u+\nabla p=f$、$\operatorname{div}u=0$、$u(\cdot,0)=0$。
• (b) すべての $0\le t\lt1$ で $\operatorname{supp}u(\cdot,t)\cup\operatorname{supp}p(\cdot,t)\subset K$。さらに、十分小さいすべての $t\ge0$ で $u(\cdot,t)=p(\cdot,t)=0$（Proposition 10.1 の式 (10.2)）。
• (c) $\sup_{0\le t\lt1}\|u(t)\|_{L^2(\mathbb{R}^3)}\lt\infty$ かつ $\limsup_{t\uparrow1}\|u(t)\|_{L^\infty(\mathbb{R}^3)}=\infty$。
• (d) 同じ外力と初期値に対して、運動エネルギーが一様に有界（$\sup_{t\ge0}\frac12\int_{\mathbb{R}^3}|v(x,t)|^2dx\lt\infty$）な $\mathbb{R}^3\times[0,\infty)$ 上の滑らかな解 $(v,P)$ は存在しない。

主張 2（原稿 Corollary 10.6）。 任意の $\nu\gt0$ に対し、$\mathbb{T}^3\times[0,\infty)$ 上の滑らかで時間についてコンパクト台の外力があって、粘性 $\nu$・初期速度 $0$ の周期的な Navier–Stokes 方程式の解 $(U,P_{\mathrm{per}})$ は $[0,1)$ で滑らかで $\limsup_{t\uparrow1}\|U(t)\|_{L^\infty(\mathbb{T}^3)}=\infty$ を満たし、同じ初期値と外力に対する大域的な滑らかな周期解は存在しない。

原稿は、主張 1 が Fefferman の問題文の命題 (C) を、主張 2 が命題 (D) を成り立たせると述べている。

これは 2026-10-08 の時点で、査読誌に掲載された結果でも、独立した数学者による検証が公表された結果でもない。 以下、この Remark の内容を「主張」と呼び、定理としては扱わない。

主張の各要素を、関数空間とノルムの観点で確認します。

• 粘性: すべての $\nu\gt0$。構成は $\nu=1$ で行い、変換で一般の $\nu$ に移す（命題 4.2）。
• 初期値: $u_0=0$（静止状態）。
• 外力: 時空の $C^\infty$ 級で、台は $K\times[0,2]$ に含まれる（原稿 Lemma 10.3）。$t=0$ の近くでは $0$。特異時刻 $t=1$ を越えて滑らかである点が本質で、$t\lt1$ だけで滑らかな外力なら 命題 4.1 により自明に作れます。外力は発散なしとは限りません（原稿 §10.1）。
• 解: $t\lt1$ で $C^\infty$ 級、空間の台は固定されたコンパクト集合 $K$ の中。
• 爆発のノルム: $L^\infty$（速度の大きさそのもの）。爆発は 1 点（原点、時刻 $1$）で起こり、原稿 Theorem 3.1 (iv) によれば、$z=0$ の平面上の半径 $\sqrt{2X_{\mathrm{in}}\tau}$ の円周に沿って $u_\theta=\tau^{-1/2-h}(e_0+O(\tau^{2h}))$、$e_0\gt0$ です。
• 有界に保たれる量: 運動エネルギー（$L^2$ ノルム）と、散逸の時間積分 $\int_0^1\|\nabla u\|^2dt$（原稿 Lemma 10.4）。

注意 3.2（この主張が扱っていないこと）

Clay の 4 つの命題との対応は次のとおりです（命題の定義は第 9 章の 定義 2.3［Navier–Stokes 方程式 I］）。

命題	外力	主張との関係
(A) 全空間で滑らかな大域解が存在	$f\equiv0$	扱っていない。主張が正しくても真偽は決まらない
(B) 周期的な場合に滑らかな大域解が存在	$f\equiv0$	扱っていない
(C) 全空間で解が破綻するデータが存在	$f\ne0$ を許す	主張 1 が成り立たせると原稿が述べる
(D) 周期的な場合に解が破綻するデータが存在	$f\ne0$ を許す	主張 2 が成り立たせると原稿が述べる

主張の外力は恒等的に $0$ ではありません（$0$ なら 命題 5.1 (i) のエネルギーの評価から $u=0$ になります）。爆発は外力が流体に運動量を注入し続けることで起こり、外力なしで静止状態から何かが起こるわけではありません。第 9 章の注意（注意 2.4［Navier–Stokes 方程式 I］）のとおり、(C) は (A) の否定を導かないので、「外力なしの方程式で、滑らかな初期値から有限時間で特異点が生じるか」という、物理的に最も関心のある問いは、主張が正しいとしても開いたままです。問題文は 4 つの命題のどれか 1 つの証明を求めているので、(C) の証明は問題文の意味では解決に当たりますが、それと (A)(B) の解明とは別のことです。

4. 主張の正否に依らない 3 つの検算

この節の命題は、主張の構造を支える初等的な事実で、主張が正しいかどうかに関係なく成り立ちます。主張の「外力」が何を意味するのか（残差を外力と呼び替えるだけで何が変わらないのか）、粘性の大きさは本質的でないこと、全空間の解を周期的な問題へ移せることを、計算で確かめるのがねらいです。

命題 4.1（残差を外力と定めれば方程式は成り立つ）✓

(i) $I$ を区間、$u$ を $\mathbb{R}^3\times I$ 上の滑らかな発散なしベクトル場、$p$ を滑らかな関数とする。$f:=R_\nu(u,p)$ とおくと、$(u,p)$ は外力 $f$ の Navier–Stokes 方程式の滑らかな解である。

(ii) さらに $w$ を滑らかな発散なしベクトル場、$\pi$ を滑らかな関数とすると、$\nu=1$ で

$$R(u+w,p+\pi)=R(u,p)+L_u(w,\pi)+\operatorname{div}(w\otimes w),\qquad L_u(w,\pi)=\partial_tw+(u\cdot\nabla)w+(w\cdot\nabla)u-\Delta w+\nabla\pi,$$

ここで $(\operatorname{div}(w\otimes w))_j=\partial_i(w_iw_j)$ である（原稿 §3.3 の増分の恒等式）。

証明（命題 4.1）

(i) $f$ の定義がそのまま方程式 $\partial_tu+(u\cdot\nabla)u-\nu\Delta u+\nabla p=f$ で、$\operatorname{div}u=0$ は仮定です。

(ii) $R$ のうち $u$ と $p$ について 1 次の項は和に分かれます。2 次の項は $((u+w)\cdot\nabla)(u+w)=(u\cdot\nabla)u+(u\cdot\nabla)w+(w\cdot\nabla)u+(w\cdot\nabla)w$ で、$\operatorname{div}w=0$ から $(w\cdot\nabla)w_j=w_i\partial_iw_j=\partial_i(w_iw_j)$ です。

∎

(i) により、問題は次のように言い換えられます。

$\mathbb{R}^3\times[0,1)$ 上の滑らかな発散なしの $u$ と $p$ で、$u$ は $t\uparrow1$ で非有界になり、しかも残差 $R_\nu(u,p)$ とそのすべての時空の導関数が、$t=1$ を越えて滑らかに延びるものを作れ。

残差の個々の項（$\partial_tu$、$(u\cdot\nabla)u$、$\Delta u$、$\nabla p$）は発散してかまいません。それらの和だけが滑らかであればよい、というのが構成の自由度です。逆に、この条件がどれほど強いかは 例 11.2 で確かめます。(ii) は、速度に補正 $w$ を加えたときの残差の変化を、$w$ について線形な部分（背景の中での $w$ の発展）と 2 次の部分（$w$ 自身の運動量流束の発散）に分ける式です。原稿は、2 次の部分の平均で背景の残差を打ち消します。

命題 4.2（粘性の変換）✓

$(u,p)$ を $\mathbb{R}^3\times[0,T)$ 上の滑らかな発散なしの場と関数とし、$f=R_1(u,p)$ とおく。$\nu\gt0$ に対し

$$u_\nu(x,t)=\sqrt\nu\,u\Bigl(\frac{x}{\sqrt\nu},t\Bigr),\qquad p_\nu(x,t)=\nu\,p\Bigl(\frac{x}{\sqrt\nu},t\Bigr),\qquad f_\nu(x,t)=\sqrt\nu\,f\Bigl(\frac{x}{\sqrt\nu},t\Bigr)$$

とおくと、$R_\nu(u_\nu,p_\nu)=f_\nu$、$\operatorname{div}u_\nu=0$ である。時間は変わらず、

$$\|u_\nu(t)\|_{L^\infty}=\sqrt\nu\,\|u(t)\|_{L^\infty},\qquad\|u_\nu(t)\|_{L^2}^2=\nu^{5/2}\|u(t)\|_{L^2}^2,\qquad\operatorname{supp}u_\nu(\cdot,t)=\sqrt\nu\,\operatorname{supp}u(\cdot,t)$$

である。

証明（命題 4.2）

$y=x/\sqrt\nu$ とおくと、$x$ についての微分 1 回ごとに $\nu^{-1/2}$ が掛かります。各項は

$$\partial_tu_\nu=\sqrt\nu\,\partial_tu,\quad(u_\nu\cdot\nabla_x)u_\nu=\sqrt\nu\cdot\sqrt\nu\cdot\nu^{-1/2}(u\cdot\nabla_y)u=\sqrt\nu\,(u\cdot\nabla_y)u,$$$$\nu\Delta_xu_\nu=\nu\cdot\sqrt\nu\cdot\nu^{-1}\Delta_yu=\sqrt\nu\,\Delta_yu,\quad\nabla_xp_\nu=\nu\cdot\nu^{-1/2}\nabla_yp=\sqrt\nu\,\nabla_yp$$

（右辺はすべて $(y,t)$ での値）です。よって $R_\nu(u_\nu,p_\nu)(x,t)=\sqrt\nu\,R_1(u,p)(y,t)=f_\nu(x,t)$ です。$\operatorname{div}_xu_\nu=(\operatorname{div}_yu)(y,t)=0$。ノルムは $dx=\nu^{3/2}dy$ から $\int|u_\nu|^2dx=\nu\cdot\nu^{3/2}\int|u|^2dy$ です。

∎

この変換は、第 9 章のスケール変換（粘性を変えない $u_\lambda=\lambda u(\lambda x,\lambda^2t)$）とは別のもので、長さの単位だけを取り替えて粘性を $1$ に正規化します。原稿 §3 の冒頭と式 (10.22)(10.23) がこれに当たり、上の計算はその検算です。$f\in C_c^\infty$、台のコンパクト性、$t=1$ での爆発は、どれもこの変換で保たれます。したがって $\nu=1$ で主張 1 の (a)〜(c) が成り立てば、すべての $\nu\gt0$ で成り立ちます。

命題 4.3（コンパクト台の解の周期化）✓

$\nu\gt0$ とする。$K\subset\mathbb{R}^3$ をコンパクト集合、$f$ を $\mathbb{R}^3\times[0,\infty)$ 上の滑らかな関数で $\operatorname{supp}f\subset K\times[0,T_f]$ となるもの、$(u,p)$ を $\mathbb{R}^3\times[0,1)$ 上の滑らかな解（粘性 $\nu$、外力 $f$）で、$\operatorname{supp}u(\cdot,t)\cup\operatorname{supp}p(\cdot,t)\subset K$、かつある $t_1\gt0$ について $t\le t_1$ で $u=p=f=0$ とする。$Q_0=(-\frac12,\frac12)^3$ とし、$\lambda\gt1$ を $\lambda^{-1}K\subset Q_0$ となるように取って $t_0=1-\lambda^{-2}$ とおく。$t\ge t_0$ で

$$\tilde u(x,t)=\lambda u(\lambda x,\lambda^2(t-t_0)),\quad\tilde p(x,t)=\lambda^2p(\lambda x,\lambda^2(t-t_0)),\quad\tilde f(x,t)=\lambda^3f(\lambda x,\lambda^2(t-t_0))$$

とし、$t\lt t_0$ では $0$ とする。このとき

$$U(x,t)=\sum_{k\in\mathbb{Z}^3}\tilde u(x+k,t),\qquad P_{\mathrm{per}}(x,t)=\sum_{k\in\mathbb{Z}^3}\tilde p(x+k,t),\qquad F_{\mathrm{per}}(x,t)=\sum_{k\in\mathbb{Z}^3}\tilde f(x+k,t)$$

は、$U,P_{\mathrm{per}}$ が $\mathbb{R}^3\times[0,1)$ 上、$F_{\mathrm{per}}$ が $\mathbb{R}^3\times[0,\infty)$ 上の滑らかな周期関数で、$(U,P_{\mathrm{per}})$ は同じ粘性 $\nu$・外力 $F_{\mathrm{per}}$・初期値 $0$ の周期的な Navier–Stokes 方程式の解である。$F_{\mathrm{per}}$ の時間の台は $[t_0,t_0+\lambda^{-2}T_f]$ に含まれ、$0\le s\lt1$ について $\|U(t_0+\lambda^{-2}s)\|_{L^\infty}=\lambda\|u(s)\|_{L^\infty}$ である。

証明（命題 4.3）

$t\ge t_0$ では、時間の平行移動と第 9 章のスケール変換（命題 3.5［Navier–Stokes 方程式 I］ (i)）により、$(\tilde u,\tilde p)$ は粘性 $\nu$・外力 $\tilde f$ の解です。$u,p,f$ は $t\le t_1$ で $0$ なので、$\tilde u,\tilde p,\tilde f$ は $t_0\le t\le t_0+\lambda^{-2}t_1$ で $0$ であり、$t\lt t_0$ への $0$ による延長は滑らかで、方程式は $t\lt t_0$ でも（すべての項が $0$ として）成り立ちます。元の時刻 $s=\lambda^2(t-t_0)$ が $1$ に達するのは $t=t_0+\lambda^{-2}=1$ のときです。

$\tilde u(\cdot,t),\tilde p(\cdot,t),\tilde f(\cdot,t)$ の台は $\lambda^{-1}K$ に含まれます。$\lambda^{-1}K$ は開集合 $Q_0$ のコンパクト部分集合なので、境界から正の距離 $\rho$ だけ離れており、平行移動 $\lambda^{-1}K-k$（$k\in\mathbb{Z}^3$）は互いに $2\rho$ 以上離れています。したがって和は局所有限で、各点の近傍で高々 1 つの項だけが $0$ でなく、微分と和が交換できます。とくに非線形項について

$$(U\cdot\nabla)U(x,t)=\sum_k(\tilde u\cdot\nabla)\tilde u(x+k,t)$$

が成り立ち（異なる平行移動の積は恒等的に $0$）、方程式の各項が項別に和になるので、$R_\nu(U,P_{\mathrm{per}})=F_{\mathrm{per}}$、$\operatorname{div}U=0$ です。周期性は和の定義から、初期値が $0$ であることは $t_0\gt0$ から従います。最後の等式は $\sup_x|U(x,t)|=\sup_x|\tilde u(x,t)|=\lambda\sup_y|u(y,\lambda^2(t-t_0))|$ です。

∎

非線形な方程式で解の和がまた解になるのは、台が互いに離れているからです。これは原稿 Corollary 10.6 の証明の前半の検算です。空間を縮める（$\lambda\gt1$）と時間は $\lambda^{-2}$ 倍に縮むので、特異時刻を $1$ に保つために時間を $t_0$ だけずらしています。

5. 主張から従うこと

命題 5.1（主張を仮定したときの (D) と (C)）✓

(i) （主張に依らない）$f$ を $\mathbb{R}^3\times[0,\infty)$ 上の滑らかな関数で $\sup_t\|f(t)\|_{L^2}\lt\infty$ となるもの、$(u,p)$ を $\mathbb{R}^3\times[0,1)$ 上の滑らかな解（粘性 $\nu$、外力 $f$、$u(\cdot,0)=0$）で、$u,p$ の空間の台が固定されたコンパクト集合に含まれるものとする。このとき $\|u(t)\|_{L^2}\le\int_0^t\|f(s)\|_{L^2}ds$（$0\le t\lt1$）である。

(ii) （データの条件）$u_0=0$ は条件 (4) を満たす。$f\in C_c^\infty(\mathbb{R}^3\times(0,\infty))$ は条件 (5) を満たす。命題 4.3 の $F_{\mathrm{per}}$ は条件 (8)(9) を満たす。

(iii) 注意 3.1 の主張 1 の (a)(b)(c) が $\nu\gt0$ について正しいと仮定する。 このとき、その $\nu$ について命題 (D) が成り立つ。すなわち、命題 4.3 の $F_{\mathrm{per}}$ と初期値 $0$ に対して、条件 (10)(11) と圧力の周期性を満たす $\mathbb{R}^3\times[0,\infty)$ 上の滑らかな解は存在しない。

(iv) 同じ仮定のもとで、同じ $f$ と初期値 $0$ に対して、$\mathbb{R}^3\times[0,\infty)$ 上の滑らかな解 $(v,P)$ で、エネルギーが一様に有界、かつすべての $T\lt1$ について

$$\int_0^T\|v(t)\|_{L^3}^3dt\lt\infty,\qquad\int_0^T\|P(t)-p(t)\|_{L^{3/2}}^{3/2}dt\lt\infty$$

を満たすものは存在しない。

(v) 主張 1 の (d) も正しいと仮定すれば、命題 (C) が成り立つ。

証明（命題 5.1）

(i) $T\lt1$ とします。$[0,T]$ 上で $u,p$ は滑らかでコンパクト台をもつので、第 9 章の減衰する解の条件を満たし、エネルギー等式（命題 3.3［Navier–Stokes 方程式 I］）の微分形 $\frac12\frac{d}{dt}\|u\|^2+\nu\|\nabla u\|^2=\langle f,u\rangle\le\|f\|\,\|u\|$ が成り立ちます。$\delta\gt0$ に対し $\frac{d}{dt}(\|u\|^2+\delta^2)^{1/2}=\frac{\frac{d}{dt}\|u\|^2}{2(\|u\|^2+\delta^2)^{1/2}}\le\|f\|$ なので、$u(0)=0$ から積分して $(\|u(t)\|^2+\delta^2)^{1/2}\le\delta+\int_0^t\|f\|$、$\delta\to0$ として主張を得ます。

(ii) $0$ はすべての評価を満たします。$\operatorname{supp}f\subset B_\rho\times[0,T_f]$ とし、$M_{\alpha m}=\sup|\partial_x^\alpha\partial_t^mf|$ とおくと、台の上では $1+|x|+t\le1+\rho+T_f$ なので、任意の $K\gt0$ について

$$|\partial_x^\alpha\partial_t^mf(x,t)|\le M_{\alpha m}\,(1+\rho+T_f)^K\,(1+|x|+t)^{-K}$$

で、これが (5) です。$F_{\mathrm{per}}$ は滑らかで空間について周期 $1$、時間の台がコンパクトなので、同じ計算で (9) を満たします。(8) は周期性です。

(iii) 主張 1 の (b) の後半と $f\in C_c^\infty(\mathbb{R}^3\times(0,\infty))$ から、ある $t_1\gt0$ について $t\le t_1$ で $u=p=f=0$ です。$f$ の台はコンパクトなので、$K$ を $f$ の空間の台も含むコンパクト集合に取り直しておきます。命題 4.3 により周期的な解 $(U,P_{\mathrm{per}})$ が $[0,1)$ 上にあり、主張 1 の (c) から $\limsup_{t\uparrow1}\|U(t)\|_{L^\infty}=\lambda\limsup_{s\uparrow1}\|u(s)\|_{L^\infty}=\infty$ です。同じデータに対して (10)(11) と圧力の周期性を満たす大域的な滑らかな解 $(V,P')$ があったとします。任意の $T\lt1$ について、$(U,P_{\mathrm{per}})$ と $(V,P')$ は $\mathbb{T}^3\times[0,T]$ 上の滑らかな周期解で初期値と外力が等しいので、第 9 章の 定理 8.2［Navier–Stokes 方程式 I］ (ii) により $[0,T]$ で $U=V$ です。よって $[0,1)$ で $U=V$ です。$V$ は $\mathbb{R}^3\times[0,\infty)$ で連続かつ周期的なので、コンパクト集合 $[0,1]^3\times[0,1]$ 上で有界で、$\sup_{t\lt1}\|U(t)\|_{L^\infty}\lt\infty$ となり、矛盾です。(ii) と合わせて (D) です。

(iv) そのような $(v,P)$ があったとします。$T\lt1$ とします。$u$ は $\mathbb{R}^3\times[0,T]$ で滑らかで台が $K\times[0,T]$ に含まれるので、$\nabla u$ は有界、$\sup_t\|u(t)\|\lt\infty$、$\int_0^T\|u\|_3^3\lt\infty$ です。$v$ についての仮定と合わせて、定理 8.2［Navier–Stokes 方程式 I］ (iii) の仮定がすべて満たされ、$[0,T]$ で $v=u$ です。よって $[0,1)$ で $v=u$ で、$\|u(t)\|_{L^\infty}=\sup_K|u(\cdot,t)|=\sup_K|v(\cdot,t)|\le\max_{K\times[0,1]}|v|\lt\infty$ となり、主張 1 の (c) に反します。

(v) (ii) により $u_0=0$ と $f$ は (4)(5) を満たす滑らかなデータで、(d) は「(1)(2)(3)(6)(7) を満たす解が無い」ことそのものです。

∎

注意 5.2（(iv) と (v) の差：原稿 Lemma 10.5）

(iv) の 2 つの可積分性の仮定は、Fefferman の条件 (6)(7) には含まれません。原稿はこの差を Lemma 10.5 で埋めています。そこでの議論は、原稿に書かれているとおりに要約すると、次の 3 段です。第 1 に、差 $w=v-u$ の圧力の勾配が、$g_{ij}=v_iv_j-u_iu_j$ の Riesz 変換 $\pi_*=\sum R_iR_jg_{ij}$ の勾配に一致することを、$\nabla\pi-\nabla\pi_*$ が調和で、時間について平均すると $H^{-3}$ に属する（したがって $0$ である）ことから示す。第 2 に、半径 $R$ の切り落としを掛けた圧力の流束を、重みつきの $L^6$ ノルムと Riesz 変換の $L^p$ 有界性、交換子の核の評価で抑え、重みつきの散逸に吸収する。第 3 に、差のエネルギーに Gronwall の不等式を使い、$R\to\infty$ とする。第 9 章の 注意 8.3［Navier–Stokes 方程式 I］ で述べたとおり、Riesz 変換の $L^p$ 有界性は Rikai に未収録なので、この記事では Lemma 10.5 を再証明せず、(d) を主張の一部として仮定に置きました。(i) により、主張 1 の (c) のうちエネルギーの有界性は、(a)(b) と外力の滑らかさから自動的に従います（原稿 Lemma 10.4 と同じ評価です）。

6. 証明の構成：原稿に書かれている範囲

この節は、原稿の §2「Physical description of the blowup」と §3「Proof outline」、および各節の冒頭と主要な命題の文に書かれていることの要約です。式は原稿のものを記号を合わせて書き、原稿に無い説明は足しません。各段の評価の中身（§4〜§9 と付録）は、この記事では検証していません。

6.1. 全体の流れ

原稿は $\nu=1$ で、次の順に場を作ります（原稿の図 5・図 6 に当たる流れ）。

段	何をするか	原稿の箇所
1	軸対称な主部 $(u^{(0)},p^{(0)})$ のプロファイル $E,U,\Pi$ を作る。残差は環状領域に台をもつ応力 $T$ の発散の形になり、外側は熱方程式の厳密解になる	§4（Theorem 4.6）、付録 A・B・C
2	$q^{2nh}$（$n=1,2,\ldots$）のべきで補正し、滑らかな軸対称の背景 $(u_B,p_B)$ を作る。残差は「環状領域の応力の発散」と「平坦な余り」の和になる	§5（Proposition 5.5、式 (5.41)）
3	振動を記述する補助のトーラス変数を導入し、異なるパルスの台を分離する	§6（Lemma 6.1）
4	2 族の振動パルスを作る。背景のせん断で成長し粘性で減衰する。振幅の 2 乗を、平均した運動量流束が応力 $T$ に一致するように選ぶ	§7（Lemma 7.1・7.4・7.7、Proposition 7.5）
5	角度平均した残差を打ち消す、コンパクト台の平均流の補正と、5 つの半径方向のモーメントの方程式	§8
6	補正のサイクルを繰り返して残差の減衰を改善し、縮む切り落としを掛けて総和して、局所的な場 $(u_{\mathrm{loc}},p_{\mathrm{loc}})$ を得る	§9（Proposition 9.6・9.9。まとめが Theorem 3.1）
7	空間と時間で切り落とし、残差を $t\ge1$ へ滑らかに延ばして外力とし、エネルギーの評価・比較・粘性の変換・周期化を行う	§10（Proposition 10.1、Lemma 10.2〜10.5、Corollary 10.6）

時刻を固定した半径–軸方向の断面の模式図（原稿の図 3(a) に当たる配置）。横軸は r/√τ、縦軸は z/τ^D で、どちらの実際の長さも τ → 0 で縮む。内側の芯と外側の熱の解の間の環状領域に、応力とパルスが置かれる。

6.2. 主部：2 つの長さで縮む渦

主部は軸対称で、相似座標 $(X,\eta)$ の固定された関数 $E,U,\Pi$（プロファイル）を使って

$$u^{(0)}_\theta=q^{-A}E(X,\eta),\qquad u^{(0)}_z=q^{-A}U(X,\eta),\qquad r\,u^{(0)}_r=V_0(X,\eta),\qquad p^{(0)}=q^{-2A}\Pi(X,\eta)$$

と書かれます（原稿の式 (4.3)）。$V_0$ は非圧縮条件 $\frac1r\partial_r(ru_r)+\partial_zu_z=0$ と軸での正則性から決まり、$\Pi$ は半径方向の釣り合い $\partial_rp^{(0)}=(u^{(0)}_\theta)^2/r$（遠心力と圧力勾配）と無限遠での正規化から決まります。自由に選べるのは $E$ と $U$ です。

原稿 §2.1 の物理的な説明では、芯の中心部で流体は軸に向かって渦を巻きながら流れ込み、$z=0$ に近い層の上下で軸方向に流れ出します。内向きの流れは角運動量を小さい半径に運んで旋回を強め、軸方向の流出が非圧縮性と両立させます。これは第 9 章の Burgers 渦（例 9.3［Navier–Stokes 方程式 I］）と同じ型の流れです。

スケールは、芯 $C_\tau$（相似座標の固定された範囲 $0\le X\le X_c$、$|\eta|\le\eta_c$）の上で

$$\ell_r\asymp\tau^{1/2},\qquad\ell_z\asymp\tau^{1/2-h},\qquad|u^{(0)}_\theta|,|u^{(0)}_z|\asymp\tau^{-1/2-h},\qquad|u^{(0)}_r|=O(\tau^{-1/2})$$

です（$\asymp$ は $\tau$ によらない正の定数倍で上下から抑えられること）。半径方向の方が速く縮み、芯は細長い柱になって原点に集中します。$h=0$ なら第 10 章の自己相似解（定義 6.1［Navier–Stokes 方程式 II］）と同じ等方的なスケールですが、$h\gt0$ なのでそれとは異なります。原稿は、半径方向の拡散と半径・軸方向の輸送が同じ速さ $\tau^{-1}$ で釣り合い、軸方向の拡散はそれより $\ell_r^2/\ell_z^2\asymp\tau^{2h}$ 倍だけ弱いことを述べ、この $q^{2h}$ を背景の補正の展開パラメータに使います。また原稿は、軸方向のプロファイルを $z=0$ について完全には対称にせず、わずかに上向きに偏らせると述べています（対称だと $z=0$ でパルスの増幅に要るせん断が消えるため）。

芯の外縁に必要なせん断を与えることと、次に述べる応力の条件のために、プロファイルの構成（§4 と付録 A〜C。本文と付録を合わせて約 60 ページ）が要ります。原稿 §3.2 によれば手順は、(1) 純粋に方位角方向の外側のプロファイルを作り、その裾を熱方程式の厳密解（“heat exterior”）に取り替える、(2) 軸の近くで解析的な内側の解を作る、(3) 両者をつなぎ、圧力・半径方向の速度・応力を決める 5 つの半径方向の累積積分を一致させる、(4) 環状領域で応力の向きを調整する、の 4 段です。

外側の熱の解は、原稿 Theorem 3.1 (iii) の式で

$$u=K(r,\tau)\,e_\theta,\qquad K(r,\tau)=r^{-1-2h}H_{\mathrm{ext}}\Bigl(\frac{\tau}{r^2}\Bigr),\qquad-\partial_\tau K=\partial_r^2K+\frac1r\partial_rK-\frac{1}{r^2}K$$

です（$\partial_t=-\partial_\tau$ なので、これは時間 $t$ について前向きの熱方程式）。圧力を $\partial_rp=K^2/r$ で定めれば、この場の残差は恒等的に $0$ です（例 11.1 で同じ型の厳密解を確かめます）。$r\gt0$ を固定すると、$\tau\to0$ で $K$ とその導関数は滑らかな極限をもち、特異時刻の速度は $r^{-1-2h}$ に比例します。

6.3. 残差を応力の発散として書く

主部だけでは方程式を満たしません。芯の内側ではプロファイルの方程式が主要な釣り合いを与え、外側は厳密解ですが、両者をつなぐ環状領域 $X_a\lt X\lt X_b$ に残差が残ります。原稿 §3.2 は、主部の接線方向の残差 $R^{(0)}_\theta,R^{(0)}_z$（半径方向の粘性を残し、軸方向の粘性を除いたもの）を、応力の成分 $T=(T_{r\theta},T_{rz})$ を使って

$$R^{(0)}_\theta=-\Bigl(\partial_r+\frac2r\Bigr)T_{r\theta},\qquad R^{(0)}_z=-\Bigl(\partial_r+\frac1r\Bigr)T_{rz}$$

と表します。$(\partial_r+\frac2r)T=r^{-2}\partial_r(r^2T)$、$(\partial_r+\frac1r)T=r^{-1}\partial_r(rT)$ なので、軸で正則な解は

$$T_{r\theta}(r)=-\frac{1}{r^2}\int_0^rs^2R^{(0)}_\theta(s)\,ds,\qquad T_{rz}(r)=-\frac1r\int_0^rsR^{(0)}_z(s)\,ds$$

です。$T$ が環状領域の外で $0$ になるように、残差を内側と外側で $0$ にし、モーメントの条件 $\int_0^\infty r^2R^{(0)}_\theta dr=0$、$\int_0^\infty rR^{(0)}_zdr=0$ を課します。Theorem 4.6 (ii) は、応力のプロファイルが $0\le X\le X_a$ と $X\ge X_b$ で $0$、その間のすべての点で $0$ でないことを述べています。

背景だけでは、この残差（＝背景を維持するのに要る外力）は $t\uparrow1$ で非有界です。応力の大きさは $q^{-A-1/2}=q^{-1-h}$、その半径方向の発散は $q^{-3/2-h}$ の大きさです。

6.4. 振動パルスが運動量を運ぶ

この節のねらいは、前節で残った残差を、外力ではなく流体自身の細かい振動で打ち消す仕組みを見ることです。細かい振動は平均すると $0$ でも、2 つの成分の積は平均が $0$ にならず、運動量を一方向に運ぶ働きをします。この働きが、背景の側に残った分をちょうど埋めます。

残差を打ち消すのは外力ではなく、流体自身の運動です。命題 4.1 (ii) の式で、背景 $u_B$ に発散なしの振動 $w$ を加えると、残差に $\operatorname{div}(w\otimes w)$ が加わります。$w$ の各成分の角度平均が $0$ でも、積 $w_rw_\theta$、$w_rw_z$ の平均は $0$ とは限りません（$\cos^2$ の平均は $\frac12$）。原稿は、角度と補助変数について平均した流束が

$$\begin{pmatrix}\langle w_rw_\theta\rangle\\\langle w_rw_z\rangle\end{pmatrix}=T+(\text{高次の項})$$

となるように $w$ を設計します。こうすると、平均した流束の半径方向の発散が、背景の応力の発散をちょうど打ち消します。

原稿 §2.2 と §3.3 に書かれているパルスの性質は次のとおりです。

• 形: 半径・高さ・時間について局在し、方位角方向には 1 周する（輪の形の）振動。各パルスの種は指数関数的に小さい外力が与え、その後の成長は背景のせん断が与える。
• 成長と減衰: パルスは背景のせん断からエネルギーを取り出して成長する。原稿は単純な例として、回転する流れの中で方位角速度の過剰が外向きの加速を生み、角速度が半径とともに十分速く減るならその外向きの運動が過剰をさらに増す、という機構を挙げている（演習 12.1 で初等的な模型を計算します）。同じせん断は振動の模様を変形して半径方向の波長を縮め、増幅を弱めるとともに粘性による減衰を強めるので、パルスは「成長してから減衰する」。時間の切り落としは指数関数的に小さい裾にだけ掛かる（Lemma 7.1・7.4）。
• 2 つの族: 半径方向の角運動量の流束と軸方向運動量の流束の比が異なる 2 族を使う。各点で 2 族の流束の方向 $v_1,v_2\in\mathbb{R}^2$ は 1 次独立で、目標の応力を正の係数で $T=c_1v_1+c_2v_2$（$c_1,c_2\gt0$）と表す。係数は振幅の 2 乗なので正でなければならず、これが主部のプロファイルへの制約（原稿の “admissible stress cone condition”。Lemma 4.5、Theorem 4.6 (iii)、Proposition 7.5）になる。
• スケール: 振幅 $A_{\mathrm{wave}}\asymp q^{-1/2-h/2}$、波長 $\ell_{\mathrm{wave}}\asymp q^{1/2+h/2}$（対数因子と境界での重みを除く）。背景の速度 $q^{-1/2-h}$ と芯の半径 $q^{1/2}$ に対して、どちらも $q^{h/2}$ 倍だけ小さい。振幅の 2 乗 $q^{-1-h}$ が応力の大きさに一致する。
• 非圧縮性: 波はベクトルポテンシャルの回転として作るので、厳密に発散なし（Lemma 7.7）。
• 配置: 特異時刻に近づくにつれ、順に細かい空間・時間のスケールにパルスを並べる（2 進的な帯 $q\asymp2^{-\ell}$ ごと）。台が重なる異なるパルスは補助変数の中で台を分け、積が各点で $0$ になるようにする（Lemma 6.1）。

残差を高周波の振動の 2 次の相互作用で打ち消す、という点は、第 10 章で触れた凸積分と共通です（原稿も、振動で指定された応力を実現する先行例として Daneri–Székelyhidi を挙げています）。違いは、ここでは振動が線形化方程式の解として力学的に成長・減衰するように作られ、得られる解が滑らかで、粘性項を含むすべての評価が粘性つきの方程式について行われる、と原稿が述べている点です。

6.5. 補正のサイクルと総和

最初のパルスで打ち消せるのは、残差の主要部分だけです。残りには、角度方向の高調波、補助変数についての平均、半径方向の積分の不整合があります。原稿 §3.4 と Proposition 9.6 は、これを 4 つの操作からなるサイクルで順に処理します。

1. 台をもつ $0$ でない角度方向の Fourier モードについて、非同次の波の振幅の方程式を解く。
2. 波の振幅を変えて、主パルスとの交差項で平均の応力の補正を与える。
3. 角度平均した残差のうち補助変数の平均が $0$ の部分を、速い補助時間についての微分の逆で打ち消す。
4. 5 つの半径方向のモーメントの方程式を解く（2 つは角運動量と軸方向の流量の積分を $0$ に保ち、3 つは圧力と接線方向の運動量の不整合を打ち消す）。

各操作のたびに圧力を半径方向の方程式から作り直し、残差全体を計算し直します。1 サイクルごとに、残差の減衰の指数 $\sigma_j$ が

$$\sigma_0=\frac15,\qquad\sigma_{j+1}=\sigma_j+\frac1{10},\qquad\sigma_j=\frac15+\frac{j}{10}\to\infty$$

と改善し、$m$ 階の導関数については $q^{h\sigma_j-K_m}$（$K_m$ は $j$ によらない）の形の評価になる、と原稿は述べています。

無限個の補正の和を取るために、各補正のベクトルポテンシャルと圧力に、$q=0$ の近くで $1$ で台が段とともに縮む切り落としを掛けてから回転を取ります（非圧縮性が保たれる）。和は $q\gt0$ で局所有限で、$t\lt1$ で滑らかな場 $(u_{\mathrm{loc}},p_{\mathrm{loc}})$ を定めます。結論が原稿の Theorem 3.1 で、要点は 4 つです。(i) $u=\operatorname{curl}\mathsf{A}+\mathsf{B}e_\theta$ の形で厳密に発散なし。(ii) $q$ が正の範囲では、すべての導関数が $\tau=0$ まで一様に有界。(iii) 残差が平坦: すべての導関数とすべての $N$ について $|\partial_x^\alpha\partial_t^bR(u,p)|\le Cq^N$（$X$ が有界な範囲で $q\downarrow0$）、$X\ge X_{\mathrm{ext}}$ では残差は恒等的に $0$ で場は熱の解。(iv) 内側の点列に沿った旋回の増大 $u_\theta=\tau^{-A}(e_0+O(\tau^{2h}))$。

残差は $0$ になるのではなく、特異点で無限次に消える、という形で制御されます。補題 2.3 により $q\le2\tau+2^{1/(2D)}|z|^{1/D}$ なので、$q^N$ の評価は「残差とその導関数が、時空の点 $(0,1)$ に近づくとき任意の多項式より速く $0$ に近づく」ことを意味し、これが残差を $t=1$ を越えて滑らかに延ばせることの核心です。

6.6. 切り落としと外力の延長

最後の段（§10）は、この記事で通読した部分です。

• 切り落とし（Proposition 10.1）: ポテンシャル $\mathsf{A}$、係数 $\mathsf{B}$、圧力に、空間の切り落とし $\chi_x$（軸対称、コンパクト台 $K$）と時間の切り落とし $\chi_t$（$t\le1-\tau_0$ で $0$）の積を掛け、そのあとで回転を取る。$u=\operatorname{curl}(c\mathsf{A})+c\mathsf{B}e_\theta$ は厳密に発散なしで、$(0,1)$ の近くでは局所的な場に一致し、初期時刻の近くでは $0$ になる。
• 外力の極限（Lemma 10.2）: $t\lt1$ で $f=R(u,p)$ と定める。切り落としで生じる項は $(0,1)$ から離れた所に台をもつ。原点の近くでは平坦性、原点から離れた $q\gt0$ の所では Theorem 3.1 (ii)、$r\gt0$ で $q\to0$ の所では熱の解の明示式から、$f$ のすべての導関数が $t\uparrow1$ で $\mathbb{R}^3$ 上一様に収束し、極限 $F_j=\lim\partial_t^jf$ は滑らかで台が $K$ に入る。
• $t\ge1$ への延長（Lemma 10.3）: $t=1$ での「Taylor 係数」$F_j$ を実現する関数を、$f(x,1+\sigma)=\sum_j\chi_0(b_j\sigma)\frac{\sigma^j}{j!}F_j(x)$ で作る。切り落としの幅 $1/b_j$ を速く縮めて級数を収束させる（任意の Taylor 係数をもつ滑らかな関数を作る Borel の補題と同じ方法）。これで $f\in C_c^\infty(\mathbb{R}^3\times(0,\infty))$、台は $K\times[0,2]$ に入る。$t\gt1$ での $f$ の値は延長の仕方によるが、主張に使うのは $t\lt1$ での方程式と、$f$ が全体で滑らかであることだけである。
• エネルギー（Lemma 10.4）: 命題 5.1 (i) と同じ評価。
• 比較（Lemma 10.5）: 注意 5.2 のとおり。
• 粘性と周期化: 命題 4.2 と 命題 4.3。

7. スケールの検算

原稿 §2 のスケールが、第 10 章で証明した定理と矛盾しないことを確かめます。これは主張の内部の整合性の確認で、主張の証明ではありません。

命題 7.1（原稿のスケールと既知の定理の整合）✓

$0\lt h\lt\frac1{100}$、$c_1,C_1\gt0$、$0\lt\tau_0\lt1$ とする。$\mathbb{R}^3\times[0,1)$ 上の滑らかなベクトル場 $u$ と、$0\lt\tau\lt\tau_0$ に対する可測集合 $A_\tau\subset C_\tau\subset\mathbb{R}^3$ が、$t=1-\tau$ で

• (S1) $|A_\tau|\ge c_1\tau^{3/2-h}$ かつ $A_\tau$ 上で $|u(x,t)|\ge c_1\tau^{-1/2-h}$、
• (S2) $|C_\tau|\le C_1\tau^{3/2-h}$ かつ $C_\tau$ 上で $|u(x,t)|\le C_1\tau^{-1/2-h}$、$|\nabla u(x,t)|\le C_1\tau^{-1-h}$

を満たすとする。このとき次が成り立つ。

(i) $\int_{C_\tau}|u|^2dx\le C_1^3\tau^{1/2-3h}\to0$、$\int_{C_\tau}|\nabla u|^2dx\le C_1^3\tau^{-1/2-3h}$ で、後者は $\tau$ について $(0,\tau_0)$ 上可積分である。

(ii) $\|u(t)\|_{L^\infty}\ge c_1\tau^{-1/2-h}$ で、$\int_0^1\|u(t)\|_{L^\infty}^2dt=\infty$。

(iii) LPS の指数 $(p,q)$（$\frac2p+\frac3q=1$、$3\lt q\lt\infty$）について

$$\|u(t)\|_{L^q}\ge c_1^{1+1/q}\,\tau^{-\frac{q-3}{2q}-h(1+\frac1q)},\qquad\int_0^1\|u(t)\|_{L^q}^pdt=\infty,\qquad(1-t)^{\frac{q-3}{2q}}\|u(t)\|_{L^q}\to\infty\ (t\uparrow1).$$

(iv) $\|u(t)\|_{L^3}\ge c_1^{4/3}\tau^{-4h/3}\to\infty$。

証明（命題 7.1）

(i) $\int_{C_\tau}|u|^2\le|C_\tau|\sup_{C_\tau}|u|^2\le C_1\tau^{3/2-h}\cdot C_1^2\tau^{-1-2h}=C_1^3\tau^{1/2-3h}$。同様に $\int_{C_\tau}|\nabla u|^2\le C_1\tau^{3/2-h}\cdot C_1^2\tau^{-2-2h}=C_1^3\tau^{-1/2-3h}$ です。$h\lt\frac16$ なので $\frac12-3h\gt0$、$\frac12+3h\lt1$ で、前者は $0$ に収束し、後者は $\int_0^{\tau_0}\tau^{-1/2-3h}d\tau=\tau_0^{1/2-3h}/(\frac12-3h)\lt\infty$ です。

(ii) (S1) で $|A_\tau|\gt0$ なので上限は $c_1\tau^{-1/2-h}$ 以上です。$\int_{1-\tau_0}^1\|u\|_\infty^2dt\ge c_1^2\int_0^{\tau_0}\tau^{-1-2h}d\tau=\infty$。

(iii) $\int|u|^q\ge|A_\tau|\,(c_1\tau^{-1/2-h})^q\ge c_1^{q+1}\tau^{3/2-h-q/2-qh}$ の $q$ 乗根を取ると、指数は $\frac{3}{2q}-\frac hq-\frac12-h=-\frac{q-3}{2q}-h(1+\frac1q)$ です。$p=\frac{2q}{q-3}$ 乗すると指数は $-1-hp(1+\frac1q)\le-1$ で、$\tau=0$ の近くで可積分ではありません。$(1-t)^{\frac{q-3}{2q}}$ を掛けると $\tau^{-h(1+1/q)}\to\infty$ です。

(iv) (iii) の第 1 式は $q=3$ でも同じ計算で成り立ち、指数は $-h(1+\frac13)=-\frac43h$ です。

∎

仮定 (S1)(S2) は、原稿 §2.1・§3.1 が主部 $u^{(0)}$ と芯 $C_\tau$ について述べるスケール（芯の体積 $\tau^{3/2-h}$、接線成分の大きさ $\tau^{-1/2-h}$、半径方向の微分 1 回で $\tau^{-1/2}$）を書き直したものです。$A_\tau$ としては、$E$ が正の下界をもつ相似座標の範囲を取ります。実際の解は主部に補正を加えたものなので、(S1)(S2) が解そのものについて成り立つかどうかは原稿の評価に依り、この記事では確かめていません（$L^\infty$ の下界 (ii) だけは Theorem 3.1 (iv) に明示されています）。

結論を既知の定理と照らします。

• (ii)(iii) と LPS の条件: 第 10 章の 定理 3.1［Navier–Stokes 方程式 II］ は外力つきでも成り立ち、爆発するなら臨界な混合ノルムはすべて発散しなければなりません。(ii)(iii) はそのとおりになっています。
• (iii) と爆発の速さ: 系 3.3［Navier–Stokes 方程式 II］ (ii) の下界 $(1-t)^{-(q-3)/(2q)}$ に対し、主張のスケールは $\tau^{-h(1+1/q)}$ 倍だけ速く発散します。下界と矛盾せず、しかも「ちょうど自己相似」の速さではありません。$\|\nabla u\|^2$ についても、(i) の上界 $\tau^{-1/2-3h}$ は 系 3.3［Navier–Stokes 方程式 II］ (i) の下界 $\tau^{-1/2}$ より大きく、両立します。
• (iv) と端点の定理: 定理 3.5［Navier–Stokes 方程式 II］ は外力なしの方程式の定理なので、この主張に直接は適用できません。それでも (iv) は $\|u(t)\|_{L^3}\to\infty$ を示しており、「爆発するなら $L^3$ ノルムが非有界」という型の結論と同じ向きです。$h=0$ なら $L^3$ ノルムは有界にとどまる計算になるので、$h\gt0$ はここでも効いています。
• (i) とエネルギー: 速度が発散しても芯のエネルギーは $0$ に近づきます。これは第 10 章の 命題 4.1［Navier–Stokes 方程式 II］ で見た優臨界性の具体例で、エネルギーの有界性は集中を妨げません。散逸の時間積分が有限であるための条件が $h\lt\frac16$ です。
• Reynolds 数: 芯の半径を長さに取ると、$\mathrm{Re}_\theta=|u_\theta|\ell_r/\nu\asymp\tau^{-1/2-h}\tau^{1/2}/\nu=\tau^{-h}/\nu\to\infty$、$\mathrm{Re}_r=|u_r|\ell_r/\nu=O(\tau^{-1/2}\tau^{1/2})=O(1)$ です。旋回は粘性に対して際限なく強くなり、半径方向の流入は粘性と同じ大きさで競い続けます。
• 部分正則性: 特異点は時空の 1 点 $(0,1)$ で、1 点の $\mathcal{P}^1$ 測度は $0$ なので、定理 5.6［Navier–Stokes 方程式 II］ の型の評価とも矛盾しません（原稿も、この評価は孤立した特異点を許すと述べています）。

8. 既知の結果との関係

第 10 章の結果と並べます。

結果	方程式	外力	解のクラス	示されたこと
Tao（定理 7.2［Navier–Stokes 方程式 II］）	非線形項を平均化した方程式	なし	Schwartz 級の初期値	有限時間爆発
Buckmaster–Vicol	Navier–Stokes	なし	エネルギー不等式を満たすとは限らない弱解	非一意性
Albritton–Brué–Colombo（定理 7.7［Navier–Stokes 方程式 II］）	Navier–Stokes	$L^1_tL^2_x$、初期時刻で特異	Leray–Hopf の弱解	非一意性
Elgindi	Euler	なし	$C^{1,\alpha}$ の初期値	有限時間特異性
Córdoba–Martínez-Zoroa ほか（定理 7.10［Navier–Stokes 方程式 II］）	Euler、弱散逸 Navier–Stokes	$C^{1,\frac12-\epsilon}$ など、有限の正則性	特異時刻まで古典解（Euler の結果では原点以外で滑らか）	有限時間特異性
2026 年の主張（注意 3.1。査読は未了）	Navier–Stokes	$C_c^\infty$（特異時刻を越えて滑らか）	特異時刻まで滑らか、コンパクト台	有限時間爆発

• Tao の平均化方程式との違い。 Tao の結果は方程式を変えており、元の方程式については「エネルギー恒等式と調和解析の評価だけでは (A) を証明できない」という障害を示すものでした（系 7.4［Navier–Stokes 方程式 II］）。主張は方程式を変えず、代わりに外力を入れます。爆発の仕組みも、スケールからスケールへエネルギーを順送りする回路ではなく、1 つの縮む渦と、それを支える流束を運ぶパルスの組です。
• 非一意性の結果との違い。 凸積分の解は滑らかでなく、Albritton–Brué–Colombo の外力は初期時刻で特異です。主張の外力は全時刻で滑らかで、解は特異時刻まで滑らかです。そのため第 9 章の弱強一意性が使え、「この解以外に滑らかな解は無い」と結論できます（命題 5.1）。非一意性ではなく、解の破綻の主張です。
• Córdoba–Martínez-Zoroa の系列との関係。 外力を許して特異性を作るという方針は、この系列が先に確立したものです。原稿は文献 [6][7][8] としてこれらを挙げ、「大きいスケールのひずみが小さいスケールの渦度を増幅する」構成と、外力の正則性を制御する戦略に言及したうえで、自分の構成では増幅された擾乱の役割が異なる（縮む背景の渦に欠けている力を、パルスの平均の運動量流束が供給する）と述べています。
• 判定条件・部分正則性との関係。 前節のとおり、矛盾しません。

9. 検証の状況

9.1. Lean による形式化：一次資料で確認したこと

リポジトリ openai/NavierStokesAndEuler について、読んで確かめた事実は次のとおりです（確認日 2026-10-08）。

• ツールは Lean v4.34.0-rc2（lean-toolchain）と Mathlib。ライセンスは Apache-2.0。GitHub の記録ではコミットは 2 個で、日時は 2026-09-08 と 2026-09-10。
• 検査の対象は ComparatorChallenges/NavierStokes.json に定められており、定理名は NavierStokes.Comparator.navier_stokes_breakdown_R3 と NavierStokes.Comparator.navier_stokes_breakdown_periodic、許される公理は propext・Quot.sound・Classical.choice（Lean の標準の 3 公理）です。
• 定理文（挑戦文）のファイル ComparatorChallenges/NavierStokes.lean は、Formal Conjectures プロジェクト（google-deepmind/formal-conjectures）の Navier–Stokes の問題文の形式化を、依存を Mathlib だけにする形で改変したものだと、ファイルの冒頭と ComparatorChallenges/README.md に書かれています。挑戦文の証明は sorry で、これは「この文を証明せよ」という形式です。
• formalization.yaml は、主結果の sorry の個数を 0、使用公理を上の 3 つ、形式化の方法を「エージェント（モデル GPT-6 Astra、枠組み Codex）」、レビューの状態を “self-assessed”（自己評価） と記載しています。これはリポジトリ自身の申告です。
• NavierStokes/ComparatorSolution.lean は、挑戦文と同じ文の定理を、プロジェクト内の別の定理（ComparatorBridge.…）から 1 行で導いています。その先の証明本体は読んでいません。

9.2. 定理文と Fefferman の条件の照合

挑戦文のファイルを全文読み、第 9 章の 定義 2.2［Navier–Stokes 方程式 I］ と 1 行ずつ対応させました。Lean では速度を v x t、圧力を p x t、外力を f x t（位置が先、時刻が後）と書き、空間は EuclideanSpace ℝ (Fin 3) です。

Fefferman	Lean の定義（挑戦文のファイル）	照合の結果
(1) 運動方程式	navier_stokes: すべての x と t ≥ 0 で、derivWithin (v x ·) (Set.Ici 0) t + fderiv ℝ (v · t) x (v x t) = nu • Δ (v · t) x - gradient (p · t) x + f x t	項ごとに一致。時間微分は $[0,\infty)$ 内での微分（$t=0$ では片側）。fderiv … x (v x t) は $x$ での微分を方向 $v$ に適用したもので $(v\cdot\nabla)v$
(2) 非圧縮	div_free: すべての x、t ≥ 0 で divergence (v · t) x = 0。divergence は fderiv のトレース	一致
(3) 初期条件	initial_condition: v x 0 = u₀ x	一致
$u_0$ が滑らかで発散なし	InitialVelocityCondition: ContDiff ℝ ∞ u₀ と divergence u₀ x = 0	一致
(4) $u_0$ の減衰	decay: すべての m : ℕ、K : ℝ に対し C があって ‖iteratedFDeriv ℝ m u₀ x‖ ≤ C / (1 + ‖x‖) ^ K	同値。原文は多重指数 $\alpha$ ごとの偏導関数、Lean は $m$ 階導関数（多重線形写像）のノルム。有限次元なので互いに定数倍で抑え合う。K は負も動くが、負の K の条件は正の K の条件から従うので差は無い
(5) $f$ の減衰	ForceConditionDecay.decay: ‖iteratedFDerivWithin ℝ m (↿f) (Set.univ ×ˢ Set.Ici 0) (x, t)‖ ≤ C / (1 + ‖x‖ + t) ^ K（t ≥ 0）	同値。時空の $m$ 階導関数全体のノルムで、$\partial_x^\alpha\partial_t^m$ をすべて含む
(6)(11) 解の滑らかさ	velocity_smooth・pressure_smooth: ContDiffOn ℝ ∞ (↿v) (Set.univ ×ˢ Set.Ici 0)	一致。閉じた半空間 $\mathbb{R}^3\times[0,\infty)$ 上の $C^\infty$
(7) エネルギー有界	integrable: 各 t ≥ 0 で MemLp (‖v · t‖) 2。globally_bounded_energy: ∃ E, ∀ t ≥ 0, (∫ x, ‖v x t‖ ^ 2) < E	一致。可積分性を別に課しているのは、Lean の積分が可積分でない関数に $0$ を返すため
(8) 周期性	IsOnePeriodic u₀、および各 t ≥ 0 で IsOnePeriodic (f · t)（各座標方向の単位ベクトルを足しても値が同じ）	一致
(9) $f$ の時間減衰	ForceConditionPeriodic.decay: 上と同じ形で右辺が C / (1 + t) ^ K	一致
(10) と正誤表	isOnePeriodic_velocity、isOnePeriodic_pressure	一致。正誤表の圧力の周期性を解の条件に入れている
(C)	navier_stokes_breakdown_R3 (nu : ℝ) (hnu : nu > 0): ∃ u₀ f, InitialVelocityConditionDecay u₀ ∧ ForceConditionDecay f ∧ ¬ (∃ v p, NavierStokesExistenceAndSmoothnessRn nu u₀ f v p)	一致
(D)	navier_stokes_breakdown_periodic: 同じ形で …Periodic の条件	一致

照合できた範囲。 上の表のとおり、挑戦文の 2 つの定理は、命題 (C)(D) を素直に形式化したものと読めます。量化子の順序（すべての $\nu\gt0$ に対し、データが存在して、解が存在しない）も原文と同じです。(D) で圧力の周期性を解の条件に入れることは、存在しないと主張する解の範囲を狭める（主張を弱める）向きの変更ですが、問題文の正誤表に従ったものです。

照合できなかった範囲。

• Mathlib の Δ（Laplacian）、gradient、iteratedFDerivWithin、ContDiffOn の定義が、この版で通常の意味（成分ごとのラプラシアン、勾配、片側の高階導関数、閉集合上の $C^\infty$）と一致することは、標準の定義であるという前提で読んでおり、Mathlib のソースで確かめていません。
• 解答側の証明（ComparatorBridge より先）は読んでいません。Lean のビルドも、Comparator による検査も実行していません。したがって「挑戦文が 3 公理だけから導かれている」ことは、リポジトリの記載として知っているだけです。
• 形式証明が、166 ページの原稿に書かれた構成（§6 の各段）と同じ議論を形式化したものかどうかは確かめていません。formalization.yaml は Theorem 1.1 と Corollary 10.6 を挑戦文の 2 定理に対応させていますが、挑戦文は「解が存在しないデータがある」という存在の主張で、爆発の時刻・台・スケールといった原稿の定理の詳細は挑戦文には現れません。

9.3. Lean の検査が保証すること、しないこと

Comparator による検査が通ることが意味するのは、形式化された 2 つの命題が、Lean の論理と 3 つの公理のもとで証明されていることです。証明の途中の議論がどれほど長くても、人が読めなくても、この保証は変わりません。これは通常の査読では得られない種類の確かさです。

一方で、次のことは Lean の検査の外にあります。

• 形式化された命題が、意図した数学的な命題（Fefferman の (C)(D)）と同じであること。これは定義を人が読んで確かめるしかなく、前項の表はその試みです。
• 原稿の記述と形式証明が対応していること、原稿の証明が人に理解できる形で正しいこと。
• 検査に使う Lean の核、Mathlib、検査の道具自体の正しさ（独立な核による再検査 enable_nanoda が設定されているのは、この点を補強するためです）。

9.4. 査読・掲載・Clay 数学研究所

• 査読・掲載: 原稿にも Lean リポジトリにも、査読誌への投稿・掲載の記載はありません。独立した数学者による検証結果の公表は、確認した範囲では見つかっていません。原稿の版の更新・訂正・正誤表も、確認した範囲ではありません。
• Clay 数学研究所の問題ページ: 2026-10-08 に確認した時点で、この主張への言及は無く、問題の状態の表示は “Active” です。
• ミレニアム賞の規則: Clay 数学研究所の規則のページ（2018 年 9 月 26 日に理事会が採択した改訂版と記されている）によれば、研究所が解決を検討する条件は、(1) 規則が定める媒体での公表、(2) 公表から少なくとも 2 年の経過、(3) 世界の数学界で広く受け入れられていること、です。研究所は直接の投稿を受け付けないとも書かれています（規則のページは要点を確認したもので、全文の精読ではありません）。
• 報道だけが出典の事柄: OpenAI は賞金を請求する意向が無いと述べている、Lean による形式化と検証に GPT-6 Astra で 17 時間を要した、と報じられています（The Next Web、2026-09-09）。OpenAI の発表ページを読めていないので、一次資料では確認していません。モデルの稼働時間やエージェントの数など、その他の数字はこの記事には書きません。

10. 公表の経緯：一次資料で確認できた事実

同じ時期に、外力つきの特異性について別のグループの結果が公表されています。誰に先取権があるか、誰の説明が正しいかは、この記事では判断しません。出典ごとに、確認できた事実を日付順に並べます。

先行研究（arXiv の記録）。

• 2023-09-15: Córdoba・Martínez-Zoroa、外力つき 3 次元 Euler 方程式の特異性（arXiv:2309.08495）。
• 2024-07-09: Córdoba・Martínez-Zoroa・Zheng、弱散逸 Navier–Stokes 方程式（arXiv:2407.06776）。
• 2024 年: Córdoba・Martínez-Zoroa、滑らかな源つきの 2 次元非圧縮多孔質媒質方程式（arXiv:2410.22920。OpenAI の原稿の文献 [7] による）。

Buckmaster の声明（PDF の作成日時は 2026-09-08 02:38 UTC。本文に日付の記載は無い。本人による記述であり、以下は「声明にそう書かれている」という事実）。

• 声明の当日（本文の表現は「本日」）、Alpöge と共同で 3 つの結果を公開した: 滑らかな外力つきの有限時間爆発を、非圧縮多孔質媒質方程式、Boussinesq 方程式、3 次元非圧縮 Euler 方程式について。
• 弱散逸 Navier–Stokes 方程式についても爆発を得たと考えているが、Lean による検証が終わっていないので公開していない。
• この研究計画の基本的な着想は Córdoba と Martínez-Zoroa のもので、両氏の粗い外力での結果を、言語モデルの助けを借りて滑らかな外力と Euler 方程式まで進めた。
• Boussinesq と Euler の結果を得たのは 8 月 15 日、Lean で検証したのは 8 月 22 日。
• 9 月 3 日に OpenAI の数学者に連絡し、9 月 6 日の通話で、OpenAI の内部のモデルが外力つき Navier–Stokes 方程式の有限時間爆発の証明を出したこと、その対象は $\mathbb{R}^3$ と $\mathbb{T}^3$ の滑らかな外力つきの方程式（Fefferman の (C)(D)）であること、証明は約 100 ページであることを伝えられた。声明の筆者は OpenAI の証明を見ていない。

Buckmaster・Alpöge 側の公開物（確認できたもの）。

• Lean リポジトリ tristanbuckmaster/fluid_lean の最初のコミットは 2026-09-08（GitHub の記録）。euler-blowup と boussinesq-blowup の README に定理文があり、Euler については「$[0,T)\times\mathbb{R}^3$ 上で滑らか、初期値はコンパクト台、外力は $[0,T)\times\mathbb{R}^3$ で滑らかですべての導関数が有界、固定された球の外で $0$、渦度の上限が $t\uparrow T$ で発散し $\int_0^T\|\omega\|_\infty dt=\infty$」という内容です。
• 2026-09-15: Alpöge・Buckmaster・Coiculescu、多孔質媒質方程式について外力を時空で滑らかにした結果（arXiv:2609.16470）。
• Boussinesq と Euler の原稿は著者のページに PDF があるとされますが、この記事では読んでいません。

OpenAI 側の公開物（確認できたもの）。

• 原稿 PDF の作成日時は 2026-09-08、Lean リポジトリの最初のコミットも 2026-09-08（GitHub の記録）。
• 原稿は先行研究として Córdoba と Martínez-Zoroa の 3 編を挙げて構成の違いを述べている。原稿の本文と文献表に、Buckmaster・Alpöge の 2026 年の結果への言及は無い（文献表にあるのは Buckmaster–Vicol の 2019 年の論文）。
• 同時に公表された Euler 方程式の原稿は、Lean リポジトリの README によれば外力なしの方程式についてのものです。

OpenAI 側の説明（発表ページを読めなかったため、報道による）。

• 報道によれば、OpenAI の発表ページは Alpöge と Buckmaster の同時期の外力つき Euler 方程式の仕事に言及し、OpenAI が扱った Euler 方程式は外力なしの版であると述べている（The Next Web、2026-09-09）。
• 報道によれば、OpenAI の研究者は、両氏の仕事に依拠したことを否定し、相手のプロンプトや証明は使っていないと述べた（Scientific American、2026-09-08）。

両者が公開した結果の数学的な違いを表にします。

	Buckmaster・Alpöge（および Coiculescu）	OpenAI
方程式	多孔質媒質方程式、2 次元 Boussinesq、3 次元 Euler（いずれも非粘性）	3 次元 Navier–Stokes（すべての $\nu\gt0$）、3 次元 Euler
外力	あり（滑らか）	Navier–Stokes はあり（$C_c^\infty$）、Euler はなし
Clay の問題との関係	粘性が無いので問題文の命題そのものではない	Navier–Stokes の結果が (C)(D) の主張
公開された形式化	Euler と Boussinesq の Lean（多孔質媒質は原稿に形式化済みと記載）	Navier–Stokes と Euler の Lean

11. 例

例 11.1（Lamb–Oseen 渦：熱方程式を満たす旋回）

純粋に方位角方向で高さによらない流れ $u=V(r,t)e_\theta$ を考えます。$\operatorname{div}u=\frac1r\partial_\theta V=0$ で、円柱座標の Navier–Stokes 方程式は

$$\text{半径方向: }-\frac{V^2}{r}=-\partial_rp,\qquad\text{方位角方向: }\partial_tV=\nu\Bigl(\partial_r^2V+\frac1r\partial_rV-\frac{V}{r^2}\Bigr)$$

になります（移流項 $(u\cdot\nabla)u=-\frac{V^2}{r}e_r$ は向心加速度だけ）。圧力を $\partial_rp=V^2/r$ で定めれば、$V$ が右の線形方程式を満たす限り、外力なしの厳密解です。軸方向の渦度 $\omega=\frac1r\partial_r(rV)$ で書くと、右の方程式は 2 次元の熱方程式 $\partial_t\omega=\nu\Delta\omega$ の回転対称な場合です。点渦 $\omega_0=\Gamma\delta$ から出発すると、熱核（定理 2.3［熱方程式］ の基本解）により

$$\omega=\frac{\Gamma}{4\pi\nu t}e^{-r^2/(4\nu t)},\qquad V=\frac{\Gamma}{2\pi r}\Bigl(1-e^{-r^2/(4\nu t)}\Bigr)$$

で、これが Lamb–Oseen 渦です。原稿の「熱の外側」（§6.2 の $K$）は同じ方程式の別の解で、$r^{-1}$ の代わりに $r^{-1-2h}$ の減衰をもち、時間は特異時刻に向かって進みます。非線形項が圧力に吸収されて線形方程式に帰着するので、この部分には外力が要りません。

例 11.2（ちょうど自己相似な爆発は、滑らかな外力では作れない（1 次元 Burgers 方程式））

命題 4.1 によれば、どんな $u$ も残差を外力にすれば解です。それでも問題が難しいことを、1 次元の $\partial_tu+u\partial_xu-\partial_x^2u=f$ で確かめます。滑らかで有界な $W\not\equiv0$ を使って $u(x,t)=\tau^{-1/2}W(x/\sqrt\tau)$（$\tau=1-t$）とおくと、第 10 章の 命題 6.2［Navier–Stokes 方程式 II］ と同じ計算で、残差は

$$f(x,t)=\tau^{-3/2}G\Bigl(\frac{x}{\sqrt\tau}\Bigr),\qquad G(y)=\frac12W+\frac12yW'+WW'-W''$$

です。$f$ が点 $(x,t)=(0,1)$ の近傍まで連続に延びるとします。$y$ を固定して曲線 $x=y\sqrt\tau$ に沿って $\tau\to0$ とすると、点は $(0,1)$ に近づき、$f=\tau^{-3/2}G(y)$ が有界であるためには $G(y)=0$ でなければなりません。$y$ は任意なので $G\equiv0$、すなわち $f\equiv0$ です。つまり、ちょうど自己相似な形の爆発を連続な外力で支えることはできず、外力なしのプロファイル方程式を解くしかありません。Burgers 方程式ではそれは不可能です（有界な解についての最大値原理により $\sup|u(t)|\le\sup|u(0)|$ となるはずですが、$\sup|u(t)|=\tau^{-1/2}\sup|W|$ は増加します）。

同じ議論は 3 次元の Navier–Stokes 方程式でも成り立ち、ちょうど自己相似な $(u,p)$ の残差が $(0,T)$ まで連続に延びるなら、プロファイルは Leray の方程式を満たし、プロファイルが $L^3(\mathbb{R}^3)$ に属する場合は 定理 6.3［Navier–Stokes 方程式 II］ によって排除されます。これはこの記事での観察で、原稿に書かれていることではありませんが、原稿の構成がちょうど自己相似ではなく（2 つの指数 $\frac12$ と $\frac12-h$、相似座標に収まらないパルス）、残差を「$0$」ではなく「無限次に平坦」にする形になっていることと符合します。

12. 演習

演習 12.1標準

非粘性で、定常な旋回流 $u=V(r)e_\theta$ に、軸対称で $z$ だけに依存する（半径方向の波数が $0$ の）小さい擾乱 $(u_r',u_\theta',u_z')(z,t)$ を加える。$r$ を固定して係数を定数とみなす近似で、線形化した方程式が

$$\partial_tu_r'-\frac{2V}{r}u_\theta'=0,\qquad\partial_tu_\theta'+\Bigl(V'+\frac Vr\Bigr)u_r'=0$$

となることを認めて、$\partial_t^2u_r'=-\Phi\,u_r'$、$\Phi=\frac{1}{r^3}\frac{d}{dr}(rV)^2$ を導け。$V=cr^{-1-2h}$（$c\ne0$、$h\gt0$）ではどうなるか。

解答

第 1 式を $t$ で微分して第 2 式を代入すると $\partial_t^2u_r'=\frac{2V}{r}\partial_tu_\theta'=-\frac{2V}{r}\bigl(V'+\frac Vr\bigr)u_r'$ です。$\frac{1}{r^3}\frac{d}{dr}(r^2V^2)=\frac{2rV^2+2r^2VV'}{r^3}=\frac{2V}{r}\bigl(\frac Vr+V'\bigr)$ なので $\Phi$ に一致します。$\Phi\gt0$ なら振動、$\Phi\lt0$ なら $e^{\sqrt{-\Phi}\,t}$ で成長します（Rayleigh の判定条件: 角運動量の 2 乗 $(rV)^2$ が外向きに減るとき不安定）。$V=cr^{-1-2h}$ では $(rV)^2=c^2r^{-4h}$ が減少するので $\Phi=-4hc^2r^{-4-4h}\lt0$ で、成長率は $\sqrt{-\Phi}=2\sqrt h\,|c|\,r^{-2-2h}$ です。原稿が「角速度が半径とともに十分速く減ると、方位角速度の過剰が外向きの運動を生み、それが過剰を増す」と述べる機構の、最も単純な模型です。$r^{-1}$（Lamb–Oseen 渦の遠方）では $\Phi=0$ でちょうど中立で、$h\gt0$ が不安定の側に倒します。原稿の実際のパルスは軸対称ではなく、軸方向のせん断、波数の変形、粘性による減衰を含むので（Lemma 7.1・7.4）、これはその解析の代わりにはなりません。

演習 12.2標準

特異時刻 $t=1$ で、外側の場が $|u|=c\,r^{-1-2h}$（$0\lt h\lt\frac12$、$D=\frac12-h$）で、それが成り立つ範囲が $|z|\le1$、$|z|^{1/(2D)}\le r\le1$ であるとする（補題 2.3 により $\tau=0$ では $q=|z|^{1/D}$ で、$X=r^2/(2q)$ が一定以上という条件が $r\ge\text{定数}\cdot|z|^{1/(2D)}$ になることに対応する）。この範囲のエネルギー $\int|u|^2dx$ が有限であるための $h$ の条件を求めよ。

解答

円柱座標で $dx=2\pi r\,dr\,dz$ なので

$$\int|u|^2dx=2\pi c^2\int_{-1}^1\int_{|z|^{1/(2D)}}^1r^{-2-4h}\,r\,dr\,dz=2\pi c^2\int_{-1}^1\frac{|z|^{-4h/(2D)}-1}{4h}\,dz$$

です。$\int_{-1}^1|z|^{-2h/D}dz$ が有限である条件は $\frac{2h}{D}\lt1$、すなわち $2h\lt\frac12-h$、$h\lt\frac16$ です。命題 7.1 (i) で散逸の時間積分が有限である条件として出た $h\lt\frac16$ と同じ値です。特異時刻の速度場は原点で非有界ですが、エネルギーは有限です。

演習 12.3易

注意 3.1 の主張 1 で、外力 $f$ が $\mathbb{R}^3\times[0,1)$ 上で恒等的に $0$ ではありえないことを示せ。また、$f$ が $t\lt1$ でだけ滑らかでよいなら、主張 1 の (a)(b)(c) を満たす例が初等的に作れることを説明せよ。

解答

前半: 命題 5.1 (i) により $\|u(t)\|_{L^2}\le\int_0^t\|f(s)\|ds$ です。$f\equiv0$ なら $u\equiv0$ で、$\limsup\|u(t)\|_\infty=\infty$ に反します。後半: $\phi\in C_c^\infty(\mathbb{R}^3;\mathbb{R}^3)$ を $\operatorname{curl}\phi\ne0$ となるベクトル場、$g$ を $[0,1]$ 上の滑らかな関数で $t\le\frac14$ で $0$、$t\ge\frac12$ で $1$ となるものとし、

$$u(x,t)=g(t)\,\tau^{-1/2}\,(\operatorname{curl}\phi)\bigl(x/\tau^{1/2}\bigr),\qquad p=0$$

とおきます。$u$ は $t\lt1$ で滑らか、発散なしで、台は $\operatorname{supp}\phi$ を含む固定された球に入り、$t\le\frac14$ で $0$ です。$t\ge\frac12$ で $\|u(t)\|_\infty=\tau^{-1/2}\sup|\operatorname{curl}\phi|\to\infty$、$\|u(t)\|_2^2=\tau^{-1}\tau^{3/2}\|\operatorname{curl}\phi\|^2=\tau^{1/2}\|\operatorname{curl}\phi\|^2$ は有界です。$f=R_\nu(u,0)$ とすれば $t\lt1$ で方程式が成り立ち（命題 4.1）、(a)(b)(c) の形の条件はすべて満たされます。しかしこの $f$ は $t\ge\frac12$ で $\tau^{-3/2}G(x/\sqrt\tau)$ の形で、例 11.2 と同じ理由により、$G\equiv0$ でない限り $t=1$ を越えて連続に延びません（$G\equiv0$ は、コンパクト台の $\operatorname{curl}\phi$ が Leray のプロファイル方程式を満たすことを意味し、第 10 章の 定理 6.3［Navier–Stokes 方程式 II］ により起こりません）。主張の内容のすべては、「$f\in C_c^\infty(\mathbb{R}^3\times(0,\infty))$」の一語にあります。

参考文献

確認日はすべて 2026-10-08 です。

1. OpenAI, “Finite time blowup for Navier–Stokes”（原稿、166 ページ、PDF の作成日 2026-09-08）。https://cdn.openai.com/pdf/32d9f210-8b73-45e0-91bc-82a30aef8a9a/navier-stokes.pdf — 査読を経ていない原稿。
2. OpenAI, Lean リポジトリ NavierStokesAndEuler。https://github.com/openai/NavierStokesAndEuler
3. C. L. Fefferman, “Existence and smoothness of the Navier–Stokes equation”, Clay Mathematics Institute の問題文（正誤表つき）。https://www.claymath.org/wp-content/uploads/2022/06/navierstokes.pdf
4. Clay Mathematics Institute, “Navier–Stokes Equation”（問題ページ）https://www.claymath.org/millennium/navier-stokes-equation/、および規則のページ https://www.claymath.org/millennium-problems/rules/。
5. D. Córdoba, L. Martínez-Zoroa, “Blow-up for the incompressible 3D-Euler equations with uniform $C^{1,\frac12-\epsilon}\cap L^2$ force”, arXiv:2309.08495。
6. D. Córdoba, L. Martínez-Zoroa, F. Zheng, “Finite time blow-up for the hypodissipative Navier Stokes equations with a force in $L^1_tC^{1,\epsilon}_x\cap L^\infty_tL^2_x$”, arXiv:2407.06776。
7. T. Buckmaster の声明（PDF の作成日 2026-09-08）。https://cims.nyu.edu/~tristanb/statement.pdf — 当事者による記述。
8. Lean リポジトリ fluid_lean。https://github.com/tristanbuckmaster/fluid_lean
9. L. Alpöge, T. Buckmaster, M. P. Coiculescu, “Extending the Córdoba-Martínez-Zoroa IPM Blow-Up to Uniformly Space-Time Smooth Forcing”, arXiv:2609.16470。
10. 報道（二次情報）: The Next Web（2026-09-09）https://thenextweb.com/news/openai-navier-stokes-proof-published-millennium-prize、Scientific American（2026-09-08）https://www.scientificamerican.com/article/openai-claims-blockbuster-math-breakthrough-amid-swirl-of-controversy/。