ナビエ・ストークス方程式：「F = ma」の流体バージョン

「水は難しい」——なぜ流体は厄介なのか

高校物理で習う「F = ma」は、非常に強力な式です。ボールを投げれば落下軌道が計算できる。車がブレーキを踏めば停止距離が求まる。この式が成立するのは、「ボール」や「車」がひとつの物体として扱えるからです。

では、空気や水はどうでしょうか。

空気の塊に力を加えると、変形し、流れ、隣の空気の塊を押し、その押された塊がさらに隣を押す——。物体としての「境界」が存在しない。これが流体力学を難しくしている本質的な理由です。

ナビエ・ストークス方程式（以下 N-S 方程式）は、この「境界のない連続した物質」に対して F = ma を適用するために、19世紀の物理学者たちが編み出した答えです。

---

まず「F = ma」を思い出す

ニュートンの運動方程式は：

$$F = ma$$

力 $F$ を質量 $m$ で割れば加速度 $a$ が得られる、と読んでもいい：

$$a = \frac{F}{m}$$

N-S 方程式も、本質はまったく同じです。左辺に「加速度」、右辺に「力を質量で割ったもの」が並びます。ただし相手が流体なので、「加速度」の書き方が少し複雑になります。そこから説明しましょう。

---

1. 加速度の部分——「流体の加速度」はなぜ2項あるのか

(a) 時間依存の加速部分：その場所で風が強くなっていく

$$\frac{\partial \mathbf{W}}{\partial t}$$

$\partial / \partial t$ は「時間微分」、つまり「時間が経つにつれてどれだけ変化するか」です。

具体的なイメージ： あなたが公園のベンチに座ったまま、風速計を持っているとします。最初は無風だったのに、5秒後には風速 3 m/s になっていた。この「同じ場所にいながら風速が変わっていく変化率」が $\partial \mathbf{W}/\partial t$ です。

場所は固定。時間だけが動いている。これが「時間依存の加速」です。

(b) 場所依存の加速部分：流れに乗って速い場所へ運ばれていく

$$(\mathbf{W} \cdot \nabla)\mathbf{W}$$

こちらは少し直感に反します。「時刻は変わっていないのに、加速度が生じる」という項です。

具体的なイメージ： 川の流れを想像してください。上流は流れが穏やかで、岩のある場所を過ぎると急に流れが速くなる。水の粒子が「穏やかな場所」から「速い場所」へと流れ込んでいくとき、その粒子は加速しています。でもこれは時計の針とは関係ない——単に「速い場所へ移動した」から速くなったのです。

つまり、流体の加速度には2種類の原因があります。

• 「その場所の風が時間とともに強くなる」（時間依存）
• 「流れに乗って、もともと速い場所へ移動してしまう」（場所依存）

N-S 方程式の左辺は、この2つを足したものです：

$$\underbrace{\frac{\partial \mathbf{W}}{\partial t}}{\text{時間依存}} + \underbrace{(\mathbf{W} \cdot \nabla)\mathbf{W}}{\text{場所依存}}$$

これを合わせて「実質加速度（物質微分）」と呼びます。流体粒子が実際に経験する加速度の全体です。

---

2. ナブラ（∇）——「全方向への微分」を一文字で表す記号

∇ とは何者か

高校数学で習う微分 $dy/dx$ は、「$x$ 方向への変化率」を表します。しかし空間は3次元です。$x$ 方向、$y$ 方向、$z$ 方向——すべての方向への変化率を一度に扱いたい。

そのために生まれた「演算子（計算の命令記号）」が ∇（ナブラ） です。

$$\nabla \equiv \left( \frac{\partial}{\partial x},\ \frac{\partial}{\partial y},\ \frac{\partial}{\partial z} \right)$$

「$x$ 方向に微分せよ、$y$ 方向に微分せよ、$z$ 方向に微分せよ」という3つの命令を、矢印のような1つの記号にまとめたものです。∇ 単体は「命令書」であり、何かに作用して初めて意味を持ちます。

∇ の使い方：3つの顔

顔 1：スカラーに作用する「勾配（gradient）」 → $\nabla f$

気圧のような「1つの数値で表される量（スカラー）」に ∇ を作用させると、「どの方向に向かって気圧が最も急激に高くなるか」を示すベクトルが得られます。

山の地形図を思い浮かべてください。等高線が密集しているほど勾配が急です。$\nabla f$ は「この場所で、最も急な登り方向を指す矢印」です。気圧の $\nabla p$ なら、「最も気圧が高くなる方向」を指します。

顔 2：ベクトルに作用する「発散（divergence）」 → $\nabla \cdot \mathbf{W}$

風速ベクトル $\mathbf{W}$ に「$\nabla\cdot$（ナブラ・ドット）」を作用させると、スカラーが出てきます。これは「その点から流体が湧き出しているか、吸い込まれているか」を表します。

水道の蛇口を想像してください。蛇口の真上では水が四方八方に広がっていく——これが $\nabla \cdot \mathbf{W} > 0$（発散）の状態。逆に排水口の真上では水が吸い込まれていく——これが $\nabla \cdot \mathbf{W} < 0$（収束）の状態です。

顔 3：ベクトルに作用する「回転（curl）」 → $\nabla \times \mathbf{W}$

「$\nabla\times$（ナブラ・クロス）」は流体の「渦巻き具合」を表します。台風の中心では $|\nabla \times \mathbf{W}|$ が大きく、一様な風が吹いているだけの平野ではゼロに近い。

---

3. N-S 方程式の全体像——改めて F = ma として読む

以上を踏まえて、N-S 方程式を改めて眺めてみましょう：

$$\underbrace{\rho_a \left( \frac{\partial \mathbf{W}}{\partial t} + (\mathbf{W} \cdot \nabla)\mathbf{W} \right)}_{\text{「ma」の部分：流体の実質加速度 × 密度}} = \underbrace{-\nabla p}_{\text{気圧の押す力}} + \underbrace{\mu \nabla^2 \mathbf{W}}_{\text{粘性（ねばり）による力}} + \underbrace{\mathbf{f}_{ext}}_{\text{外力}}$$

右辺の各項も直感的に読めるはずです。

$-\nabla p$ は気圧の勾配、つまり「気圧の高い方から低い方へ向かう押す力」です。高気圧から低気圧へ風が吹くのはこの項のせいです。

$\mu \nabla^2 \mathbf{W}$ は粘性項です。$\nabla^2$（ラプラシアン）は「∇ を2回作用させたもの」で、「周囲との速度差を均そうとする力」を表します。蜂蜜のような粘い流体ではこの項が支配的になり、水や空気では比較的小さくなります。

$\mathbf{f}_{ext}$ はすべての外力をまとめた項です。鉄道近接飛行の文脈では、列車通過による誘起風がここに入ります。

---

まとめ：N-S 方程式が「難しい」本当の理由

N-S 方程式の構造を理解すると、「難しさ」の正体が見えてきます。

右辺——力の部分——は比較的素直です。問題は左辺の $(\mathbf{W} \cdot \nabla)\mathbf{W}$ という項です。ここには $\mathbf{W}$ が2回登場しています。求めたいものである $\mathbf{W}$ が、式の中で自分自身と掛け算されている。これを数学では「非線形項」と呼びます。

F = ma でボールの軌道を解くとき、「ボールの速度がボールの速度に影響を与える」なんてことは起きません。でも流体では起きる。速く流れているほど、周囲をより激しく引きずる。その「引きずり」がさらに流れを変える。この自己参照的な構造が、N-S 方程式を一般には解析的に解けなくしている根本的な理由です。

現代のスーパーコンピューターが膨大な計算資源を使って数値的にしか解けないのも、この非線形項のせいです。 そして「N-S 方程式の完全な解が常に存在するか」という問いは、100万ドルの懸賞がかかったミレニアム懸賞問題のひとつとして、2025年現在もなお未解決のままです。

F = ma という単純な原理が、流体という連続体に適用された瞬間に人類最難問のひとつへと変貌する——これがナビエ・ストークス方程式の持つ、知的な凄みです。

---