ナビエ–ストークス方程式 — 流体力学の支配方程式

流体力学の中心に鎮座する方程式がナビエ–ストークス方程式だ。オイラーの運動方程式が粘性を無視した理想流体の式だったのに対し、ナビエ–ストークス方程式は粘性の効果を含めた「本物の」流体の運動方程式である。天気予報、航空機の設計、血流のシミュレーション、海流の予測——流体が関わるほぼすべての科学技術がこの方程式に依存している。

方程式の全貌

非圧縮性のニュートン流体に対するナビエ–ストークス方程式は次のように書ける。

$ρ (\frac{\partial v}{\partial t} + (v \cdot \nabla) v) = - \nabla p + μ \nabla^{2} v + f$

ここで $ρ$ は密度、 $v$ は速度場、 $p$ は圧力、 $μ$ は粘性係数（動粘度ではなく粘度）、 $f$ は体積力だ。オイラーの運動方程式と見比べると、右辺に $μ \nabla^{2} v$ という項が加わっただけだと気づく。

オイラーの運動方程式

$ρ \frac{D v}{D t} = - \nabla p + f$ 　——粘性なし

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

$ρ \frac{D v}{D t} = - \nabla p + μ \nabla^{2} v + f$ 　——粘性あり

たった一つの項の追加だが、この項が方程式の振る舞いを劇的に変える。

各項の物理的意味

方程式の左辺と右辺に現れる各項を整理しよう。

慣性項

ρ (v \cdot \nabla) v

流体粒子が空間的に速度の異なる場所へ移動することで生じる加速度。この項の非線形性が、流体力学の解析を困難にしている最大の原因だ。

圧力勾配

- \nabla p

圧力が高い側から低い側へ流体を押す力。流れを駆動する主要な要因の一つ。

粘性項

μ \nabla^{2} v

速度場のラプラシアン（2 階微分）に粘性係数を掛けたもの。流体内部の摩擦を表し、速度分布を滑らかにする方向に働く。はちみつのようにドロドロした流体ほど $μ$ が大きい。

粘性項 $μ \nabla^{2} v$ は拡散方程式に現れるラプラシアンと同じ構造を持つ。熱が温度差を均す方向に伝わるように、粘性は速度差を均す方向に作用する。この類似性は偶然ではなく、運動量の拡散という統一的な視点から理解できる。

連続の方程式との組み合わせ

ナビエ–ストークス方程式だけでは未知数が足りない。非圧縮性流体の場合、速度場の 3 成分 $(v_{x}, v_{y}, v_{z})$ と圧力 $p$ の計 4 つの未知関数に対して、運動方程式は 3 本（各成分に 1 本ずつ）しかない。4 本目の方程式として連続の方程式が必要になる。

$\nabla \cdot v = 0$

この 4 本の方程式と適切な境界条件・初期条件を与えれば、原理的には流れの完全な時間発展が決まる。

レイノルズ数 — 慣性と粘性の綱引き

ナビエ–ストークス方程式の振る舞いを支配する最も重要な無次元量がレイノルズ数 $R e$ だ。代表的な流速 $U$ 、代表的な長さ $L$ 、動粘度 $ν = μ / ρ$ を用いて、

$R e = \frac{UL}{ν}$

と定義される。レイノルズ数は慣性力と粘性力の比を表しており、その大きさによって流れの性質が劇的に変わる。

Re の範囲	流れの特徴
$R e ≪ 1$	粘性支配（ストークス流）
$R e \sim 1$	慣性と粘性が拮抗
$R e ≫ 1$	慣性支配（乱流が発生しうる）

レイノルズ数が小さい極限では非線形の慣性項が無視でき、方程式は線形化される。これがストークス方程式であり、微生物の遊泳やコロイド粒子の沈降など、低レイノルズ数の世界を記述する。一方、レイノルズ数が大きくなると非線形項が支配的になり、流れは不安定になって乱流へ遷移する。

未解決のミレニアム問題

ナビエ–ストークス方程式は 19 世紀前半にクロード＝ルイ・ナビエとジョージ・ガブリエル・ストークスによって定式化された。200 年近い歴史を持つにもかかわらず、この方程式には数学的に未解決の根本問題が残っている。

3 次元のナビエ–ストークス方程式に対して、滑らかな初期条件から出発した解が有限時間内に「爆発」（速度や渦度が無限大に発散）しないかどうかは、証明されていない。この問題はクレイ数学研究所が選定した 7 つのミレニアム懸賞問題の一つであり、解決には 100 万ドルの賞金がかかっている。

現時点では、2 次元の場合については解の滑らかさと一意性が証明されているが、3 次元では大域的な正則性が未解決のままだ。

任意の時刻まで解が滑らかに存在し続けることの証明。有限時間で特異点が生じるかどうかすらわかっていない。

数値解法の重要性

解析的に解ける場合は極めて限られるため、実用上はコンピュータによる数値シミュレーション（CFD: Computational Fluid Dynamics）が不可欠だ。有限差分法、有限体積法、有限要素法、スペクトル法などの手法が開発され、現代の工学設計はほぼすべて CFD に依存している。

ただし、乱流の完全な直接数値シミュレーション（DNS）はレイノルズ数の 3 乗に比例する計算コストを必要とし、実用的な高レイノルズ数の流れに対しては依然として近似モデル（RANS やラージ・エディ・シミュレーション）に頼らざるを得ない。ナビエ–ストークス方程式は書き下すこと自体は簡単だが、その解を求めることは物理学と数学の最前線に位置する挑戦であり続けている。