Navier-Stokes方程式の厳密解：Poiseuille流、Couette流、Stokes拡散

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

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

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

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

ナビエ–ストークス方程式の厳密解とは、速度場 $u$ と圧力 $p$ の組であって、

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

を数値近似なしに閉じた形で満たすものです。

鍵となるメカニズムは、非線形移流項 $(u \cdot \nabla)u$ の消失または簡略化です。平行流 — デカルト座標で $u = (u(y,t),\, 0,\, 0)$、円筒座標で $u = (0,\, 0,\, u(r,t))$ の形を持つ流れ — では、速度場は自動的に発散なしであり、移流項は恒等的に消えます。速度が流れ方向に成分を持たず、流れ方向に沿って変化しないためです。運動量方程式は線形PDEに帰着されます — 定常流の場合はODEに。

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

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

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

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

ナビエ–ストークス方程式の実際の閉形式解を探しているなら、まずはこれらの古典的な解から始めるのがよいでしょう。以下の各公式は、特定の幾何学と境界条件の仮定のもとで、非圧縮方程式の簡約版を解いたものです。

重要な注意点は、これらが一般解ではないということです。三次元のナビエ–ストークス問題を難しくしている非線形な挙動の大部分を幾何学的設定が取り除いてくれるからこそ、厳密に解けるのです。

この表は、滑りなし境界条件のもとでの非圧縮ニュートン流体のナビエ–ストークス方程式について、標準的な閉形式解をまとめたものです。ここで $\mu$ は粘性係数、$\nu=\mu/\rho$ は動粘性係数、$R$ は管の半径、$h$ は平板間隔、$dp/dx$ は一定の流れ方向圧力勾配です。

これらの公式は、バチェラーのAn Introduction to Fluid Dynamics、ランダウ＝リフシッツのFluid Mechanics、シュリヒティングのBoundary-Layer Theoryといった大学院レベルの流体力学の教科書では標準的です。以下の導出では、公式を引用結果として頼るのではなく、帰着を直接示します。

ポアズイユ流れ（ハーゲン・ポアズイユ流れ）とは、まっすぐな円管内を流れるニュートン流体の、定常かつ十分に発達した流れです。速度分布は放物線で、壁面ではゼロ、中心軸上で最大になります。

$Q=\pi R^4\Delta p/(8\mu L)$ は流量と圧力降下を結び、流量が半径に非常に敏感な理由を示します。順を追った導出、仮定、計算例、適用限界については、ハーゲン・ポアズイユ方程式の詳しい導出をご覧ください。

軸対称で十分に発達した仮定形 $\mathbf{u}=u_z(r)\mathbf{e}_z$ では移流項が消え、軸方向方程式を積分すると $u_z(r)=G(R^2-r^2)/(4\mu)$ となります。断面積分から $Q=\pi R^4\Delta p/(8\mu L)$ が得られます。

この厳密解には、非圧縮ニュートン流体、まっすぐで剛体の円管、定常かつ十分に発達した流れ、滑りなしなどの条件が必要です。ポアズイユ流れの完全な導出では、境界条件、$R^4$ 則、壁面せん断応力、流れの状態と形状による適用限界を詳しく扱います。

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

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

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

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

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

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

単純クエット流。上板が速度 $U$ で $x$ 方向に動き、底板は静止、$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}.$$

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

相似変数 $\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（直接数値シミュレーション）コードの標準的な検証ケースとして機能します。
• ジェフリー–ハーメル流。定常・純径方向流 $u_r(r,\theta)$（半角 $\alpha$ のくさび内）。流れ関数 $\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$ 上の解が、滑らかで発散なし、急減少する初期データに対して、すべての時刻で $C^\infty(\mathbb{R}^3 \times [0,\infty))$ に留まるかどうかです。厳密解はこれに対処しません。なぜなら、それらは完全な非線形ダイナミクスが決して関与しない解空間の部分空間に属するからです。

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

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

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

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

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

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

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