ナビエ–ストークス方程式の厳密解

ナビエ–ストークス方程式は悪名高いほど非線形です。では、どうやって正確に解くことができるのでしょうか？

答えは対称性です。流れの幾何学が十分に単純なとき — 真っ直ぐな管、二枚の平板、無限平面 — 速度は一方向のみを向き、一つまたは二つの座標に沿ってのみ変化します。多くの古典的な対称的設定では、非線形項が消えるか、十分に単純化されて方程式が線形PDEに帰着されます — 定常の場合にはしばしば紙と鉛筆で解けるODEになります。

こう考えてください：ナビエ–ストークス方程式はすべての可能な流体運動を記述しています。しかし、流体を非常に秩序だった状況 — 旋回なし、カオスなし、すべてが一方向に行進 — に強制すれば、方程式の複雑さの大部分は無関係になります。ナビエ–ストークスの難しい部分は、流体が自分自身を押し回すフィードバックループです。これらの対称的な流れでは、押し返すものがなく、非線形自己移流項が脱落するのです。

これらの厳密解は好奇心の対象ではありません。流体力学教育の基礎であり、数値コードのベンチマークであり、実際の流れがいつどのように理想化から逸脱するかを理解する出発点です。

ナビエ–ストークス方程式の厳密解とは、

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0$$

を、数値近似なしに閉じた形で満たす速度場 $u$ と圧力 $p$ です。

鍵となるメカニズムは、非線形移流項 $(u \cdot \nabla)u$ の消失または簡略化です。平行流 — デカルト座標で $u = (u(y,t),\, 0,\, 0)$、円筒座標で $u = (0,\, 0,\, u(r,t))$ の形を持つ流れ — では、速度場は自動的に発散なしであり、移流項は恒等的に消えます：

$$(u \cdot \nabla)u = u\,\partial_x u\, \hat{e}_x = 0$$

$u$ が流れ方向座標 $x$ に依存しないためです。運動量方程式は線形PDEに帰着されます — 定常流の場合はODEに。

流れが定常（$\partial_t u = 0$）であれば、残るのは

$$\nu \Delta u = \nabla p,$$

ポアソン方程式であり、その解が粘性流の古典的厳密解です：ポアズイユ流、クエット流、およびそれらの派生。

非定常平行流（$\partial_t u \neq 0$ だが移流は消える）の場合、方程式は拡散方程式 $\partial_t u = \nu\, \partial_{yy} u$ となり、ストークスの第一・第二問題を与えます。

ポアズイユ流（ハーゲン–ポアズイユ流とも呼ばれる）は、ナビエ–ストークス方程式の最も重要な厳密解であり、多くのエンジニアが最初に出会うものです。

水が長い真っ直ぐな管を定常的に流れている様子を想像してください。両端の一定の圧力差が流れを駆動します。管壁は静止しているので、壁に触れている流体はゼロ速度で固定されています（すべりなし条件）。壁から遠い管の中心の流体が最も速く動きます。

得られる速度分布は放物線です：壁でゼロ、中心で最大、その間は滑らかな曲線。管を切断して断面を見ると、速度は逆さまの椀のように見えるでしょう。

この放物線プロファイルには驚くべき帰結があります。管を通る総流量は管の半径の4乗に比例します。半径を倍にすると流量は16倍になります。これがハーゲン–ポアズイユの法則であり、わずかな動脈の狭窄が血流を劇的に減少させうる理由を説明しています。

仮定：ポアズイユ流は、流体が非圧縮でニュートン的（一定の粘性）、流れが定常で完全に発達（入口からまだ加速中ではない）、そして流れが層流 — 滑らかで秩序ある — であることを仮定します。実際には、管内流はレイノルズ数約2,300で乱流に遷移します。この数値は理論から導出されたものではなく、経験的な観察です。

半径 $R$ の円管内の定常・完全発達・軸対称流を考えます。一様な圧力勾配 $dp/dx < 0$ が軸方向 $x$ に駆動します。円筒座標 $(r, \theta, x)$ で速度場は $u = (0,\, 0,\, u(r))$ の形を取ります。

非圧縮条件 $\nabla \cdot u = 0$ は自動的に満たされます。$u$ が $x$ に依存しないので移流項は消えます。軸方向運動量方程式は

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

あるいは等価的に

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

$dp/dx$ が定数なので、これは $r$ のODEです。境界条件 $u(R) = 0$（すべりなし）と $du/dr|_{r=0} = 0$（対称性）で二回積分すると、放物線速度分布が得られます：

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

中心線での最大速度は $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$、平均速度は $\bar{u} = u_{\max}/2$ です。

断面にわたって積分すると、体積流量のハーゲン–ポアズイユの法則が得られます：

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

ここで $\Delta p > 0$ は管長 $L$ にわたる圧力降下です。

$R^4$ 依存性が特徴的です：半径の小さな変化が流量に大きな変化を生じます。血行力学では、これが動脈狭窄に対する血流の感度の基盤です。

妥当性：この解は、非圧縮・ニュートン・定常・完全発達・層流に適用されます。乱流への遷移は実験的に $\text{Re} = \bar{u} D / \nu \approx 2{,}300$（$D = 2R$ は管径）で起きます。この閾値はナビエ–ストークス方程式から予測する定理ではなく、経験的観察です。レイノルズ数と層流–乱流遷移の議論については、レイノルズ数、乱流、そして小スケールが重要な理由をご覧ください。

クエット流は、二枚の平行な板の間に閉じ込められた流体が、片方の板が動きもう片方が静止しているときの厳密解です。

テーブルの上に平らに置いたトランプの束を想像してください。一番上のカードを横にずらすと、下のカードも動きます — それぞれが上のカードよりも少しだけ少なく。速度は底のゼロから上の板の速度まで直線的に変化します。

これが単純クエット流です：直線的な速度分布。底の板に接する流体は静止し、上の板に接する流体は板とともに動き、その間のすべてが直線的に補間されます。圧力勾配は不要で — 運動は純粋に動く境界によって駆動されます。

流路方向にも圧力勾配を加えると、せん断駆動流と圧力駆動流の組み合わせになります。速度分布は直線分布の上に放物線が重畳されたものになります。この複合流は平面ポアズイユ–クエット流と呼ばれることがあります。

純粋な圧力駆動の場合 — 両板が静止し、圧力差で駆動される流れ — は平面ポアズイユ流であり、管内流の平板版です。分布は放物線で、中央が最も速く、両壁でゼロです。

間隔 $h$ で隔てられた二枚の無限平行板の間の定常流を考えます。$y$ を板に垂直な座標とし、底板を $y = 0$、上板を $y = h$ とします。流れは平行です：$u = (u(y),\, 0,\, 0)$。

単純クエット流。上板が $x$ 方向に速度 $U$ で動き、底板は静止、$dp/dx = 0$。運動量方程式は $d^2u/dy^2 = 0$ に帰着し、

$$u(y) = U\frac{y}{h}.$$

これはナビエ–ストークス方程式の最も単純な非自明な厳密解です：境界条件のみによって駆動される直線速度分布。

平面ポアズイユ流。両板が静止し、一定の圧力勾配 $dp/dx < 0$ が流れを駆動します。運動量方程式は

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

$u(0) = u(h) = 0$ の境界条件で。解は

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

中心線 $y = h/2$ について対称な放物線分布です。

一般的なクエット–ポアズイユ流。動く上板と圧力勾配の組み合わせは重ね合わせ

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

第一項はせん断駆動の寄与、第二項は圧力駆動の寄与です。$dp/dx$ の符号と大きさに応じて、$\mu U / h^2$ に対する比率によって、分布は単調、内部に極大を持つ、あるいは片方の壁面近くで逆流を示すことさえあります。

ポアズイユ流とクエット流は定常です — 時間とともに何も変化しません。ストークスの問題は最も単純な非定常厳密解であり、根本的なことを明らかにします：粘性は、熱が固体中を拡散するのと同様に、運動量を流体中に拡散させるのです。

ストークスの第一問題（レイリー問題とも呼ばれる）：平板の上に広大で静止した流体のプールがあると想像してください。時刻ゼロで、板が突然一定速度で横に滑り始めます。板のすぐ隣の流体は即座に引きずられますが、遠くの流体は気づくまでに時間がかかります。滑らかな境界層が板から外側へ成長し、時間とともに厚くなります。

板の上の任意の高さでの速度は、その高さと特性的な拡散長 $\sqrt{\nu t}$ の比に依存します。ここで $\nu$ は粘性、$t$ は経過時間です。流体の粘性が大きいほど、運動はより速く上方に広がります。

ストークスの第二問題：同じ設定ですが、板が一定速度で動く代わりに正弦波的に前後に振動します。振動は流体中に有限の距離しか浸透しません — より上方では、流体はほとんど気づきません。運動の振幅は高さとともに指数関数的に減衰し、薄い振動境界層を形成します。深いプールの水面を撹拌しても、表面近くの浅い領域しか乱されない理由がこれです。

両方のストークスの問題は、$y = 0$ の平板の上にある半無限流体（$y > 0$）を含みます。流れは平行です：$u = (u(y,t),\, 0,\, 0)$。移流項は消え、支配方程式は1次元拡散方程式になります：

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

ストークスの第一問題（レイリー問題）。初期条件：$y > 0$ で $u(y,0) = 0$。境界条件：$t > 0$ で $u(0,t) = U$。遠方場：$y \to \infty$ で $u \to 0$。

相似変数 $\eta = y / (2\sqrt{\nu t})$ がPDEをODEに帰着させます。解は

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

ここで $\operatorname{erfc}$ は相補誤差関数です。境界層の厚さは $\delta \sim \sqrt{\nu t}$ として成長し、拡散的な広がりの特徴です。

ストークスの第二問題。板が振動します：$u(0,t) = U\cos(\omega t)$。$u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$ の形の解を求めると

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

したがって

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

振幅は浸透深さ $\delta_s$ で指数関数的に減衰し、位相は直線的に遅れます — エネルギーがすべて粘性により散逸する横波。高い周波数ほど浸透が浅くなります。

ポアズイユ、クエット、ストークスの流れは最も有名ですが、唯一の厳密解ではありません。文献に頻繁に登場するものがいくつかあります：

• テイラー–グリーン渦 — 2次元における減衰する渦のパターン。真の渦構造を持つ有名な厳密解であり、計算流体力学コードのテストのベンチマークとして広く使われています。
• ジェフリー–ハーメル流 — 収束または発散するくさび形の流路内の流れ。流体が狭まるギャップに加速する、あるいは広がるギャップで減速する様子を捉えます。
• ヒーメンツよどみ点流 — 平板に正面から衝突する流体。流体は壁面でゼロに減速し、横に逸れます。風が建物の正面に当たるとき、あるいはジェットが表面に衝突するときに起きることをモデル化しています。

これらはそれぞれ、方程式を扱いやすくするために特定の幾何学的対称性を利用しています。専門的な文脈で重要ですが、入門的な流体力学の主力はポアズイユ流とクエット流です。

平行流を超えて、いくつかの他の厳密解の族が特定の対称性や自己相似構造を利用しています：

• テイラー–グリーン渦。2次元では $u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$ で $aA = bB$（非圧縮性）、指数的な時間減衰 $e^{-\nu(a^2+b^2)t}$。これは非線形項を含む完全な2次元ナビエ–ストークス方程式の厳密解です（非線形項は勾配になり圧力に吸収されます）。DNS（直接数値シミュレーション）コードの標準的な検証ケースとして機能します。
• ジェフリー–ハーメル流。半角 $\alpha$ のくさび内の定常・純径方向流 $u_r(r,\theta)$。流れ関数 $\psi = \nu f(\theta)$ は $\theta$ の3次非線形ODEを満たします。収束・発散流路の両方に解が存在し、高レイノルズ数で豊かな分岐構造を持ちます。
• ヒーメンツよどみ点流。平面壁に衝突する流れに対する2次元相似解。流れ関数 $\psi = \sqrt{a\nu}\, x\, f(\eta)$、$\eta = y\sqrt{a/\nu}$ がナビエ–ストークスをヒーメンツODEに帰着させます：$f''' + ff'' - f'^2 + 1 = 0$、$f(0)=f'(0)=0$、$f'(\infty)=1$。

これらすべての場合にナビエ–ストークス方程式を正確に解けるなら、なぜまだ百万ドルの未解決問題があるのでしょうか？

すべての厳密解が特別なトリックに頼っているからです。幾何学が非常に注意深く選ばれているため、方程式の最も困難な部分 — 非線形項 — が完全に消えるか、扱いやすいものに帰着されるのです。管内流は実質的に1次元です。クエット流は直線です。テイラー–グリーン渦さえも、その非線形性を圧力の中に隠しています。

ミレニアム賞問題は、特別な対称性のない一般的な3次元流れについて問うています。任意の初期条件、簡略化する幾何学なし、すべてのスケールで完全な非線形相互作用。その設定では、解が常に滑らかであり続けることは誰にも証明されておらず、そうではないことも証明されていません。

これらの厳密解は、強い対称性の仮定のもとで方程式が明示的な滑らかな解を許容することを示しています。未解決の問いは、方程式が常に機能するか — あるいは、3次元乱流の完全な一般性の中で、何か壊滅的なことが起きうるかどうかです。

上で議論したすべての厳密解は、非線形移流項 $(u \cdot \nabla)u$ を排除または自明化することで扱いやすさを達成しています。平行流では恒等的に消えます。テイラー–グリーン渦では勾配になり圧力に吸収されます。相似解では変数変換により低次元の問題に帰着されます。

クレイ・ミレニアム賞問題は、滑らかで急減少する任意の初期データを持つ3次元非圧縮ナビエ–ストークス方程式の初期値問題に関するものです — まさにこれらの簡略化がまったく適用されない領域です。

問いは

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

の $\mathbb{R}^3$ 上の解が、滑らかな初期データ $u_0 \in C_c^\infty(\mathbb{R}^3)$ に対してすべての時刻で $C^\infty(\mathbb{R}^3 \times [0,\infty))$ に留まるかどうかです。厳密解はこれに対処しません。なぜなら、それらは完全な非線形ダイナミクスが決して関与しない解空間の部分空間に属するからです。

要約すれば：厳密解は高度に対称的な設定で明示的な滑らかな解を提供します。未解決問題は、適切設定性が制限のない3次元流れに拡張されるかどうかです。正確な定式化については、ナビエ–ストークスの存在と滑らかさをご覧ください。

方程式そのものと各項の意味を理解するには、ナビエ–ストークス方程式とは何か？から始めてください。

方程式が第一原理からどのように構築されるかを見るには、ナビエ–ストークス方程式の導出をお読みください。

ポアズイユ流やクエット流のような層流がいつ乱流に崩壊するかを理解するには、レイノルズ数、乱流、そして小スケールが重要な理由をお読みください。

厳密解では答えられない百万ドルの問いを理解するには、ミレニアム賞問題：存在と滑らかさをお読みください。

本サイトでの推奨される次のページ：

• ナビエ–ストークス方程式とは何か？ — 支配方程式系とその各項
• ナビエ–ストークス方程式の導出 — 運動量収支とニュートン構成法則からPDEへ
• レイノルズ数、乱流、そして小スケールが重要な理由 — ここで記述した層流領域から完全に乱流への遷移
• ミレニアム賞問題：存在と滑らかさ — 厳密解が解決しない未解決問題の正確な定式化