ナビエ–ストークス方程式

1. なぜ重要なのか:ミレニアム問題

船の周囲を流れる水と翼の周囲を流れる空気は、同じ局所的なつり合いに従う。質量は保存され、力は運動量を変化させる。ナビエ–ストークス方程式は、この考えを連続体のモデルにしたものである。方程式を書き下すことは、そのすべての解がどう振る舞うかを証明することより、はるかに容易である。

研究状況の確認日:2026年10月13日。 OpenAIは9月8日、滑らかな外力を伴う三次元非圧縮ナビエ–ストークス方程式について、有限時間で発散する解の構成を発表し、論文とLeanのコードを公開した。クレイ数学研究所のミレニアム問題の選択肢 CとD を対象とする、大きな進展を主張する発表である。ただし、外力なしのナビエ–ストークス方程式で発散が起こることを示したものではない。別に発表された外力なしオイラー方程式の結果は粘性がゼロの場合を扱っており、混同してはならない。OpenAIの発表(9月10日更新)。

ここでは、結果の発表、形式化された命題の検査、人が読んで理解できる証明、数学界での受容を区別することが大切になる。第14節では、10月に入ってからの検討も含め、評価の状況を日付とともに説明する。その前に、エネルギーが有界でも局所的な運動が激しくなりうる理由を、物理から見ていこう。


2. 連続体の場と保存則

質量密度を 、速度を 、圧力を 、コーシー応力を とする。単位質量当たりの体積力(例えば重力)を とする。

場は、各位置と時刻に値を割り当てるものであり、多数の分子について平均した量を表す。このモデルでは、対象とする長さのスケールが、その平均化や局所的な応力の法則に意味があるほど大きいと仮定する。数学的な特異点は連続体方程式に関する概念であり、現実の分子が無限の速さに達するという予言ではない。

2.1 質量保存

第1項は局所的にどれだけ質量が蓄積するかを、発散の項は外向きの正味の質量流束を表す。密度が一定なら、質量保存は速度の発散がゼロという条件になる。流体の小さな塊は、体積を保ちながら伸びたり形を変えたりできる。

2.2 運動量保存(コーシーの式)

応力を、等方的な圧力と粘性による部分に分解する。

物質微分 は、移動する流体の塊を追いかける微分である。速度場が時間に依存しなくても、塊は加速しうる。例えば狭い通路に入る水は、速度の異なる領域へ移動する。この効果が、非線形な移流項 に対応する。


3. 構成則:ニュートン流体

等方的なニュートン流体では、

ここで は粘性係数、 は第2粘性係数であり、体積粘性係数は となる。

静止した下側の板と動く上側の板の間では、定常層流の速度は高さとともに直線的に増える。粘性によるせん断と滑りなし条件を示している。

粘性係数が空間的に一定なら、これを運動量のつり合いに代入して、圧縮性ナビエ–ストークス方程式を得る。

が一定のとき、動粘性係数を と定義する。

粘性は、一様な並進ではなく変形に応答する。その次元は 、 である。板の間の流れの図では、間隔を 、上側の板の速さを とすると、、せん断応力は となる。この簡単な厳密解は滑らかである。特異点が存在するという主張は、特定のデータに対する存在命題であって、すべての流れについての主張ではない。


4. 非圧縮流れの形

密度 が一定で、

が成り立つとき、次式を得る。

この式は、加速度が「圧力による加速度、粘性拡散、外から加えた加速度」の和に等しい、と読める。

項意味
位置を固定して測った時間変化
流体の塊が速度場の中を移動するために生じる変化
圧力による力が運動量を再配分する効果
粘性が空間的な変動をならす効果
外から加えられる単位質量当たりの力

非圧縮性は、速度や加速度がゼロであることを要求しない。圧力は速度の発散をゼロに保つ拘束として働く。この一定密度のモデルでは、圧力は独立した状態方程式で決まる量ではない。

運動量方程式の発散を取ると、適切な境界条件のもとで圧力のポアソン方程式が得られる。


5. 渦度とヘリシティ

渦度を と定義する。 非圧縮流れ()では、

成分で書くと、

の項は渦の伸長を表す。これは三次元でのみ現れる。 ヘリシティ は、渦線のトポロジーを調べる指標となる。

渦管が軸方向に引き伸ばされると、細くなり、回転が強まることがある。渦度方程式では、伸長は回転を増幅しうる一方、粘性はそれを拡散させる。しかし、それぞれの効果を単独で考えるだけでは、結び付いた運動全体が滑らかに保たれるかどうかは決まらない。

真に平面的な流れ では、渦度は平面に垂直な方向を向き、場はその方向には変化しない。したがって である。これにより二次元では主要な増幅機構の一つがなくなるが、だからといって三次元の渦が必ず発散するわけではない。

平面内の反時計回りの循環に対応する渦度は、紙面の手前を向く。丸の中の点でその方向を示している。


6. エネルギー収支:非圧縮流れ

滑らかな解について、運動量方程式と の内積を取り、静止した滑りなし境界をもつ領域 上で積分する。周期領域、または十分な減衰がある全空間でも、同じ恒等式が成り立つ。

の場合、運動エネルギーは単調に減少し、粘性による散逸率は で与えられる。

ここでのエネルギーは一定密度で規格化している。物理的な運動エネルギーは、速度の二乗の積分に を掛けたものである。移流と圧力の寄与は境界を通る流束だけになり、上記の仮定では消える。部分積分によって、非負の粘性項が得られる。壁が動く場合や境界が開いている場合には、仕事や流束の項を追加する必要がある。

6.1 有限エネルギーだけでは滑らかさを保証できない理由

積分値の上界は運動の総量を制御するが、局所的な最大値までは制御しない。スケーリングの例として、滑らかで台がコンパクトな発散ゼロの場 を一つ固定し、次のように定める。

すると、 は で大きくなる一方、 は小さくなる。これらは縮尺を変えた場の例であり、方程式の解でも、発散の証明でもない。ますます小さい体積の中で速さがますます大きくなっても、エネルギーの上限とは矛盾しないことを示している。導関数や集中を制御することが、正則性を考える際の追加の難題となる。


7. 無次元化と支配パラメータ

代表的な長さ と速さ を選び、次のように定義する。 、、、。

非圧縮流れの方程式は、次の形になる。

ここでレイノルズ数は、

また、 である。拡散の時間スケールは 、移流の時間スケールは であり、その比が になる。レイノルズ数が大きいとは、選んだスケールにおいて、粘性が変動をならす速さに比べて移流が速いことを意味する。普遍的な乱流の閾値でも、特異点の判定条件でもない。形状、擾乱、どのスケールを測るかも重要である。例えば、、、 なら、 となる。

必要に応じて、ほかの無次元数も用いる。

  • マッハ数 (圧縮性)。
  • フルード数 (重力)。
  • プラントル数 (運動量拡散と熱拡散の比)。

8. 圧縮性ナビエ–ストークス–フーリエ系の概観

内部エネルギー または温度 と、フーリエの熱流束 を含める。ここで は熱伝導率であり、上で用いた熱拡散率 とは異なる。

状態方程式(例えば理想気体の式) を加えることで、方程式系が閉じる。

ここでは圧縮による仕事と、運動エネルギーから内部エネルギーへの変換が明示的に重要となる。非圧縮流れに関するミレニアム問題の結果によって、この別の系における衝撃波形成、熱輸送、正則性の理論が解決するわけではない。


9. 数学的理論:存在と一意性

局所的な存在とは初期時刻の後のある時間区間で解が存在することであり、大域的な存在とは将来のすべての時刻で解が存在することである。古典解は、各点で方程式を満たすのに十分な導関数をもつ。ここで滑らかとは、任意の階まで微分可能であることをいう。弱解は、微分を滑らかな試験関数に移した積分形式を満たす。これらは異なる意味での解の存在である。

9.1 弱解(ルレイ–ホップ)

標準的な全空間または周期領域の設定では、発散ゼロで有限エネルギーの初期データと許容される外力に対し、ルレイ–ホップ理論はエネルギー不等式を満たす大域的弱解を与える。滑らかな解は、滑らかである間、対応するエネルギークラスの中で一意である(弱解と強解の一意性)。したがって、滑らかさが失われたとしても、すべての広義の解が存在しなくなることを意味しない。フェファーマンによる公式の問題解説。

一意性は、解のクラス、初期データ、外力に依存する。Albritton、Brué、Colombo は特定の外力に対する非一意性を示している。Hou、Wang、Yang は、特異な初期データに対し、外力なしでの非一意性を計算機援用で示したと報告している。こうした結果を、滑らかな初期データから外力なしで発散が起こるという結果に読み替えてはならない。

9.2 三次元の正則性条件

  • 渦度による解の延長条件(ビール–加藤–マイダ型): 十分な正則性を仮定した通常の全空間または周期領域の設定で、外力が時刻 を含めて十分に正則であるとき、古典解が で破綻するならば、次式が必要となる。

Luoによるナビエ–ストークス方程式の解の延長条件の解説には、外力なしの条件が述べられている。上記で外力の正則性を要求するのは、特異な外部入力によって引き起こされる破綻を除くためである。

  • プロディ–セリン条件: 適切な外力のもとで、エネルギークラスの弱解が を満たし、

であれば、その解は 上で滑らかである。

  • カファレッリ–コーン–ニーレンバーグの定理: 定理で課される外力の仮定のもとで、適切な弱解の特異点集合は、一次元放物型ハウスドルフ測度がゼロになる。これは部分正則性であり、特異点集合が空であることの証明ではない。公式の問題解説にある背景と参考文献。

これらは条件付きの命題である。どのような制御があれば破綻を防げるかを指定するものであり、その制御を三次元のすべての流れに対して与えるわけではない。上の指数は、自然なスケーリング を反映している。三次元では、エネルギーノルムだけでは臨界的なスケーリングにならない。

9.3 二次元での大域的正則性

標準的な全空間または周期領域の設定で、条件に適合する滑らかなデータと十分に正則な外力があれば、二次元非圧縮流れはすべての有限時刻で滑らかに保たれる。外力がなければ、エンストロフィー は増加しない。外力があれば生成項が加わる。渦の伸長がないため、その制御は三次元よりもはるかに強い。

9.4 圧縮性流れ

一定の条件のもとで、 則に従う気体には有限エネルギーの大域的弱解が存在する。しかし、正則性と一意性の多くは未解決である。

9.5 クレイの四つの選択肢は何を問うのか

公式の定式化では、正の粘性のもとで、次の四つの命題のいずれか一つを証明することが求められる。

選択肢領域外力と問い
A外力ゼロ。許容されるすべての滑らかな初期データについて、有限エネルギーの滑らかな大域解が存在するか
B周期領域 外力ゼロ。許容されるすべての滑らかな周期的初期データについて、滑らかな大域解が存在するか
C許容されるある滑らかなデータと、滑らかで急減衰する外力について、一様に有界なエネルギーをもつ滑らかな大域解が存在しないか
D周期領域 許容されるある滑らかな周期的データと外力について、滑らかな大域解が存在しないか

全空間では、初期データと外力は指定された減衰条件を満たす必要がある。周期的な外力にも時間方向の減衰条件が課される。CとDはゼロでない外力を許すため、どちらかを証明しても、外力なしの命題AやBを論理的に否定することにはならない。「クレイの選択肢の一つを解決する」ことと「正則性のあらゆる問題を解決する」ことは、異なる主張である。


10. 乱流とスペクトル

高 の一様・等方性乱流に関するコルモゴロフの1941年理論(K41)では、慣性領域のエネルギースペクトルは、

ここで は平均散逸率、 はコルモゴロフ定数である。

コルモゴロフの微小スケール:

K41の見積もり を使うと、 となる。このスケールを三次元で解像するには、時間積分の計算量を数える前の段階で、格子点数が 程度必要になると見積もられる。この考察が、異なる解像度とモデル化を選ぶ直接数値シミュレーション(DNS)、ラージ・エディ・シミュレーション(LES)、**レイノルズ平均ナビエ–ストークス(RANS)**の手法を動機付ける。LESとRANSには、解像されない相関を扱う閉鎖仮定が必要である。

のスペクトルは、慣性領域に関する統計的なモデルであり、すべての流れやスケールを記述するものではない。間欠性や境界は、この描像を複雑にする。特別に構成された特異解が得られても、普遍的な乱流モデルが手に入るわけではなく、数値誤差の評価が不要になるわけでもない。


11. 境界条件と領域の形状

  • 滑りなしの壁: 。
  • 自由滑り: 、。
  • 理論解析やDNSで用いる周期領域(トーラス)。
  • 適切な応力・圧力条件を指定した流入・流出境界。

境界条件の選択は、問題の適切性、エネルギー収支、数値安定性に影響する。

周期的な箱では向かい合う面を同一視し、物理的な壁は存在しない。全空間では、無限遠での減衰条件が壁の条件に代わる。これらの幾何での結果を、航空機の表面や管内の流れにそのまま適用することはできない。境界条件も数学的な問題の一部だからである。


12. 数値計算の要点

  • 射影法: を時間発展させ、 のポアソン方程式を解き、 を満たす空間へ射影する。
  • エネルギー安定な計算法: 対流項を歪対称形にすると、離散系のエネルギーを保存できる。
  • 移流・拡散を陽的に扱う計算法のCFL条件: 。
  • 発散の制御: 発散ゼロのスペクトル基底、またはinf–sup条件を満たす速度・圧力空間の組を用いた混合有限要素法。厳密に発散ゼロとするには適切な空間が必要であり、inf–sup安定性だけでは保証されない。

シミュレーションでピークの成長は観察できるが、有限の格子だけで無限大への極限を確立することはできない。本当の集中と数値的不安定性を見分けるには、格子を細かくしたときの検討、離散化誤差の制御、解析的な議論が必要である。計算機援用証明は厳密な評価を与え、形式証明支援系は符号化された定義と仮定からの推論を検査する。どちらも単なる高解像度の動画ではない。


13. 関連する理論と拡張

  • オイラー方程式: とすると粘性拡散がなくなる。極限 を取ること自体も、特に壁の近くでは繊細な問題となる。OpenAIが別に発表した外力なしオイラー方程式の結果は、この非粘性の方程式に関するものであり、外力なしで正の粘性をもつ問題の結果ではない。発表と別論文のオイラー方程式の結果。
  • 相対論的流体力学: 素朴なナビエ–ストークス理論は因果律に反する。2次のイスラエル–スチュワート理論は、輸送係数と流体の状態に適切な条件を課すことで因果性を回復できる(体積粘性を含む厳密な例)。
  • 量子流体: 二流体模型と量子渦を用いた、超流動におけるナビエ–ストークス方程式の類似系。

14. 研究の状況

2026年10月13日時点: 滑らかな外力を伴う三次元の構成が発表され、検討が進んでいる。外力なし三次元ナビエ–ストークス方程式の正則性は未解決である。二次元の正則性や、滑らかな流れの解が日常の応用で役立つことが覆されたわけではない。

14.1 正確な主張と構成の仕組み

OpenAIの論文の定理1.1は、任意の に対し、空間と時間についてコンパクトな台をもつ滑らかな外力と、静止状態から出発する 上の滑らかな非圧縮流れが存在し、次を満たすと述べている。

時刻1より前では、速度と圧力は固定されたコンパクトな空間的台をもち、外力はその時刻を越えて滑らかに延長される。論文は選択肢Cと、系10.6の周期的な構成によるDを主張している。論文の圧力の定義では、一定密度の因子を圧力に吸収している。

提案された仕組みでは、内向きにらせん運動をし、軸方向に伸びる渦を集中させる。振動的な補正とその運動量輸送が、そのままでは特異になる残りの力を打ち消す。発散する速度を好きに選び、方程式から外力を定義するだけでは不十分である。外力とそのすべての導関数が、特異時刻を越えて滑らかでなければならない。これは主張された構成の概念的な説明であり、独立した証明ではない。

「任意の粘性に対して外力が存在する」という量化では、外力が粘性に依存してよい。一つの固定された外力があらゆる粘性に打ち勝つこと、任意に小さな外力で十分なこと、あるいは静止した流体がエネルギーを与えられずに発散することを意味しない。

14.2 Lean、数学的な検討、懸賞金

公開されたLeanリポジトリには、全空間と周期領域の証明対象、およびComparatorを用いて独立に検査する手順が含まれる。形式検査器は、正確に符号化された命題と仮定に対する推論を検証する。その定義が意図した問題を表しているか、また形式的な議論が文章で書かれた証明をどのように説明しているかは、別に検討すべき問いである。本稿は、Leanプロジェクトを再ビルドしたとも、独立に検証したとも主張しない。

10月6日、Bastounis、Circelli、Hansen は、自然言語による説明とLeanの対応関係を批判的に検討した論文を公開した。例として、補題8.6で必要とする導関数の階数の違い、補題10.5における圧力評価の議論の違いを挙げている。この指摘は、形式化によって文章の証明のすべての段階が自動的に認証されたと扱うことに疑問を投げかける。ただし、それだけで外力付き発散の最終定理が反証されたり、最終的な形式命題が誤っていると示されたりするわけではない。

9月11日のクレイ数学研究所の声明は、画期的な進展と見られる発表を歓迎する一方、評価を意図的に急がず進めることを強調した。懸賞金の規則では、条件を満たす出版、その出版から少なくとも2年の経過、数学界での広い受容が、審査に先立って必要となる。発表日だけで、懸賞金に向けた所定の期間が始まったとは判断できない。ここでは受賞が決まったとは述べていない。OpenAIは懸賞金を請求する意向はないとしている。

14.3 後続研究が示すこと、まだ示さないこと

9月の後続研究についても、適用範囲に注意が必要である。PetrilloとGlimm は、外力付きの結果と、残された外力なしの問題を明確に区別している。LeiとRenによる解説はプロファイルの構成を扱い、その先の振動による相殺は後続の議論に委ねているため、独立した完全な検証ではない。Constantin、Ignatova、Vicol は、幾何、増大率、外力に特定の仮定を課したもとで正則性を示し、特異点形成を阻む結果を得ている。これはこの機構から外力を取り除く試みに関係するが、解析的な外力をもつすべての流れについての大域的正則性の定理ではない。

今後の実りある問いは明確である。主張された仮定と結論は検証されているか、仕組みを見通しよく説明できるか、外力をなくしたときに何が残るか。乱流の予測、あらゆる解のクラスでの弱解の一意性、圧縮性流れは、それぞれ別の問題である。

15. 理解を確かめよう

  1. 一様な並進: 速度と圧力が一定で外力がゼロのとき、非圧縮方程式のどの項が消えるか。空間微分も時間微分もすべて消えるので、どの正の粘性とも両立する運動である。
  2. エネルギーと最大速度: 第6.1節で と変数変換しよう。なぜエネルギーは に比例するのか。振幅の二乗が 、三次元の体積が を与えるからである。
  3. スケールを変える: と を固定して を半分にすると、レイノルズ数はどうなるか。半分になり、小さいスケールで粘性の相対的な重要性が増す。
  4. 定理を読む: 選択肢Aが未決着でもCが成り立ちうるのはなぜか。Cは選ばれたゼロでない滑らかな外力を許すのに対し、Aは外力ゼロの場合だけを扱うからである。
  5. 証明検査器: Leanが証明を受け入れた後、何を確認する必要があるか。符号化された定義、仮定、最終的な結論が意図した数学的な問いと一致していること、さらに別に書かれた説明との対応も検討する必要がある。