Navier-Stokes方程式の解き方：「解く」とは本当は何を意味するのか

ナビエ–ストークス方程式は解けるのかと問われたとき、正直な答えはこうです：「解く」の意味によります。エンジニアは航空機を設計し天気を予測するために、毎日数値的に解いています。数学者は1840年代から、単純な幾何学に対する紙と鉛筆による厳密解を知っています。しかし、あの問題 — 解が常に存在し、良い振る舞いをすることの証明 — を解いた者は誰もいません。

三つの異なるレベルがあります：

• 厳密解。特別な幾何学では、対称性が方程式を大幅に単純化するため、答えを閉じた形で書き下せます。
• 数値解。現実世界の幾何学に対しては、コンピュータがセルごと・時間ステップごとに流れを近似します。これが計算流体力学（CFD）です。
• 存在の問題。3次元のあらゆる滑らかな初期流れは、永遠に滑らかなままでいるのか？それがミレニアム賞問題であり、未解決です。

このページの残りでは、各レベルと — より重要なことに — それらの間のギャップを順に見ていきます。興味深い数学が住んでいるのは、そのギャップの中です。

非圧縮系 $$\partial_t u + (u\cdot\nabla)u = -\frac{1}{\rho}\nabla p + \nu\,\Delta u + f, \qquad \nabla\cdot u = 0.$$ を固定します。これを「解く」ことには、互いに同値でない三つの意味がありえます：

1. 閉じた形：系を厳密に満たす場 $u(x,t)$ と $p(x,t)$ を明示的に提示すること。実際には、対称性による簡約でPDE系を線形ODEや熱型方程式に変換することが必要です。
2. 離散近似：近似解 $u_h$ をメッシュ幅 $h$ の格子上で構成し、スキームの安定性と整合性の制約のもとで、$u_h$ が $h\to 0$ の極限で解に収束するようにすること。
3. 適切設定性：クレイ数学研究所の定式化は、滑らかで発散なし・急減少の任意の初期データ $u_0$ が $\mathbb{R}^3$ 上に与えられたとき（またはトーラス上の滑らかなデータに対して）、滑らかな解 $(u,p)$ がすべての $t>0$ でエネルギー有界のまま存在するかを問います。

各レベルは前のレベルよりも少ない構造を前提とし、より大きな一般性を要求します。最初の二つは確立された実践です。三つ目が未解決問題です。

幾何学が非常に単純なため、方程式の最も困難な部分が消えてしまうことがあります。真っ直ぐな管の中や二枚の平行平板の間では、流体は一方向に動き、別の方向に沿って変化します。この配置は非線形項を完全に消し去り、残るのは手で積分できる方程式です。

すべての流体力学の講義で扱われる古典的な族：

• ポアズイユ流：管や流路内の圧力駆動流で、放物線速度分布を持つ
• クエット流：動く板と固定された板の間でせん断される流体
• ストークスの第一問題：突然動き出した壁によって運動を始める流体
• ラム–オセーン渦：粘性の作用でゆっくりと広がっていく単一の渦
• よどみ点流：壁に正面から衝突する流れ

これらは博物館の展示品ではありません。あらゆるCFDコードがテストされるベンチマークであり、粘性が流れをどう形作るかについての本物の直観を運んでいます。厳密解のページでは、これらの公式を一つずつ導いています。ナビエ・ストークス方程式から速度分布と流量までの完全な簡約は、ポアズイユ流れの導出で確認できます。

平行流 $u=(u(y,t),0,0)$ に対しては、発散なし条件が自動的に成り立ち、移流項は恒等的に消えます：$(u\cdot\nabla)u = 0$。運動量方程式は $$\frac{\partial u}{\partial t} = -\frac{1}{\rho}\frac{dp}{dx} + \nu\,\frac{\partial^2 u}{\partial y^2},$$ という線形方程式に縮退します。定常状態では直接積分で解ける2階のODEになり、$y=0$ と $y=h$ の板の間では平面ポアズイユ分布 $$u(y)=\frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h-y).$$ が得られます。

このパターンは一般化できます：すべての古典的厳密解が存在するのは、何らかの対称性（平行な流線、径方向対称性、自己相似性）が、速度場とその自身による輸送との間の非線形フィードバックを除去または飼い慣らすからです。対称性を破れば閉じた形は失われます。だからこそ厳密解は、価値あるものではあっても、一般の3次元問題についてはほとんど何も語らないのです。ポアズイユ流れの厳密な導出では、この消去と半径方向問題を閉じる境界条件を詳しく示しています。

真っ直ぐな管より複雑なもの（翼、血管、暴風雨のシステム）には、公式はありません。代わりにCFDは、空間を何百万もの小さなセルに刻み、各セルの速度と圧力が隣接セルにどう応答するかを計算しながら、流れを時間的に前進させます。解いているのは元の連続な方程式ではありません。セルが縮小するにつれて答えが真の流れに近づくように構築された、離散的な代替物を解いているのです。

主な離散化の族：

• 有限差分法：格子上で微分を差分に置き換える
• 有限体積法：セル面を通るフラックスを追跡し、保存量が保存され続けるようにする（商用CFDの多くはこの方式）
• 有限要素法：単純な基底関数から解を組み立てる
• スペクトル法：流れをフーリエモードで展開する。周期領域では極めて高精度
• 格子ボルツマン法：平均が流体方程式を回復するような、簡略化された粒子集団をシミュレートする

乱流はさらにモデリングの選択の層を加えます — すべてを解像する方式（DNS）からすべてを平均する方式（RANS）まで。どこまで細かく解像する必要があるかはレイノルズ数に依存し、どの方程式を離散化するかは流れが非圧縮か圧縮性かに依存します。

中心的な構造的困難は、圧力と速度の結合です。非圧縮ナビエ–ストークス方程式には $p$ の発展方程式がありません。圧力は $\nabla\cdot u = 0$ を強制するラグランジュ乗数として作用します。標準的な取り扱いには、暫定速度を進めてから圧力ポアソン解法で発散なし場へ射影する射影法（Chorin, 1968）、SIMPLE/SIMPLEC/PISO系の圧力補正スキーム、完全連成型ソルバーがあります。

移流の離散化はトレードオフを強います：風上スキームは安定ですが数値拡散を導入し、中心スキームはより正確ですが偽の振動を起こしやすい。時間積分は、陽的スキーム（1ステップあたり安価だがCFL制限あり）、陰的スキーム（安定だが各ステップで非線形ソルブが必要）、およびIMEXの組み合わせに分かれます。

乱流に対する階層はDNS、LES、RANSです。直接数値シミュレーションは動的に活性なすべてのスケールを解像しますが、コストがレイノルズ数とともに急激に増大するため、中程度の $Re$ に限られます。ラージ・エディ・シミュレーションは大きなスケールを解像し、サブグリッド応力をモデル化します。レイノルズ平均定式化は変動場全体をモデル化し、得られるモーメント方程式を経験的モデルで閉じます。

ここにひねりがあります。エンジニアは毎日流体の流れを計算し、物理学者は方程式を完全に信頼しています。それなのに数学者は、方程式が常に良い振る舞いをする答えを与えることを証明できないのです。3次元の任意の滑らかな流れのパターンから出発して、発展させてみてください。永遠に滑らかなままでしょうか、それともエネルギーを激しく集中させ、有限時間で速度が無限大になりうるのでしょうか？誰も知りません。その問い — 存在と滑らかさ — は七つのミレニアム賞問題の一つで、100万ドルの価値があります。

すでに確定していること：

• 2次元では答えはイエスです：滑らかな解はすべての時刻で存在します。これは20世紀半ばに解決されました。
• 3次元の短時間では、滑らかな解が存在します。未知なのは、「短時間」を常に「永遠」に延長できるかどうかです。
• 3次元のすべての時刻では、より弱い種類の解（ルレイの弱解）が常に存在します。その弱解が実際に滑らかかどうかが、まさに未解決の問いです。

2次元が抵抗しないのに3次元が抵抗する理由は、渦の伸長とエネルギーが小スケールへ移る仕方に帰着します。問題が難しい理由のページで、さらに深く掘り下げています。

ルレイ（1934）は、発散なしの初期データ $u_0\in L^2(\mathbb{R}^3)$ に対する弱解の時間大域的存在を証明しました：分布の意味で方程式とエネルギー不等式を満たすが、各点での正則性は保証されない場です。これらの解の一意性は未知です。

クレイ問題はこう問います：$f\equiv 0$ かつ滑らかで発散なし・急減少の $u_0$ が $\mathbb{R}^3$ 上に与えられたとき、組 $(u,p)$ であって、$u\in C^\infty(\mathbb{R}^3\times[0,\infty))$ かつエネルギーが大域的に有界であるものは存在するか？大域正則性の証明か、許容されるデータに対する有限時間破綻の実証か、いずれかが問題を解決します。

次元による分裂は明確です。2次元の大域正則性は、ルレイ、ラジゼンスカヤ、リオンス、プロディの一連の仕事によって確立されました。鍵は、2次元の渦度が伸長なしに輸送され、強い事前評価による制御が得られることです。3次元では、滑らかな解の時間局所的存在は古典的であり、あるスケール臨界ノルムが有限に留まる限り滑らかさが持続することも知られていますが、一般のデータに対する大域正則性はどちらの方向にも未解決のままです：証明もなく、反例もありません。

コンピュータが好きな流れを何でもシミュレートできるなら、なぜ数学的な問いが重要なのでしょうか？シミュレーションと証明は異なる問いに答えるからです。シミュレーションは格子上の数値を与えます。証明は方程式そのものについての保証を与えます：解が存在し、一意であり、初期データの変化に連続的に応答するという保証です。その保証がなければ、流体力学の基礎には論理的なギャップが残ります。「方程式はうまく機能する」と「方程式がうまく機能することを証明できる」は同じ主張ではありません。

このギャップには実務的な牙があります。シミュレーションは時に破綻しますが、破綻したとき、それが数値的アーティファクト（メッシュが粗すぎる、時間ステップが大きすぎる）なのか、方程式そのものの真の特異点に触れたのかを、常に見分けられるとは限りません。正則性理論があれば、どの失敗が本物かを決着させられるでしょう。このギャップを埋めることを目指す戦略は、アプローチと部分問題のページにあります。

数値スキームの収束定理は条件付きです：典型的には、極限の解がミレニアム問題がまさに確立を求めている正則性を持つことを仮定しています。それがなければ、メッシュ細分化が生み出す列の極限は、せいぜい弱解であることしか保証されません。そして弱解のクラスには驚きがあります。凸積分法（オンサーガー予想の解決でIsettおよびBuckmaster、De Lellis、Székelyhidi、Vicolによりオイラー方程式向けに発展し、BuckmasterとVicolによりルレイ–ホップより弱いクラスのナビエ–ストークスに拡張された手法）は、分布解がひどく非一意でありうることを示しています。したがって「計算された解」という概念は、見かけほど無害ではありません。

解像度のギャップもあります。DNSは現在、大型計算機上で $Re\sim 10^4\text{-}10^5$ に達しますが、大気の流れは $Re\sim 10^9$ で動いています。工学の実践は、厳密な適切設定性理論による裏付けなしに、スケール分離の数桁にわたって方程式への信頼を外挿しています。これらはいずれも、実務におけるCFDの信頼性を損なうものではありません。証明が何を付け加えるのかを、正確に位置づけているのです。

現在の状況：

• 厳密解：対称的な流れに対するカタログは確立済みで、周辺部でゆっくりと成長し続けている
• 数値解：工学用途では成熟し目覚ましく高精度だが、乱流モデリングの負担が増えるにつれて信頼性は低下する
• 存在の問題：未解決。大域的な滑らかさの証明はなく、爆発の反例もなく、懸賞金は未請求

「解く」ことの階層は、理解の階層を映しています。ほぼあらゆる流れを計算でき、予算の許す限り細かく近似できますが、3次元方程式についての最も基本的な適切設定性の主張をまだ証明できていません。部分的な結果は届き続けており、進捗のページがそれらを追跡しています。

3次元で最も強い無条件の結果は、依然としてカファレリ–コーン–ニーレンバーグの部分正則性です：適切な弱解に対して、特異集合の1次元放物型ハウスドルフ測度はゼロです。条件付き正則性判定条件（渦度を用いるBeale–Kato–Majda — 元はオイラー方程式向けに定式化されナビエ–ストークスに適応されたもの — と、Prodi–Serrin–Ladyzhenskayaのスケール臨界可積分性条件）は、単一の臨界量の制御を完全な滑らかさに変換します。

活発な方向には、凸積分とルレイ–ホップ級より下での非一意性の構造、定量的正則性と爆発レートの結果、自己相似爆発候補シナリオの計算機支援解析、確率的正則化の問題が含まれます。スケーリング・ギャップ（既知の事前評価が方程式の自然なスケーリングに対して優臨界であること）が依然として中心的な障害であり、漸進的な改良が問題を閉じられていない理由です。