ナビエ–ストークス方程式とは何か？

流体運動の偏微分方程式への明快な入門 — 簡単な直感から数学的形式まで

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

ナビエ–ストークス方程式は、水、空気、血液などの粘性流体の運動を記述する偏微分方程式の系です。

管の中の水、飛行機の翼の周りの空気、動脈の中の血液、その他無数の流れを記述するために使われています。

大まかに言えば、この方程式は次のことを述べています：流体は、圧力、粘性、そして作用する力によって運動を変化させる。圧力は流体を押し動かし、粘性は運動の急激な差を滑らかにし、重力のような外力が流れを駆動します。

これらの方程式は単なる物理のスローガンではありません。流体力学、工学、数値シミュレーションの大部分で使われている実用的な言語です。

ナビエ–ストークス方程式は、粘性ニュートン流体の運動量収支方程式です。非圧縮の設定では、速度場 $u(x,t)$ と圧力 $p(x,t)$ を非線形PDE系を通じて結びつけます。

厳密に言えば、ナビエ–ストークスは単一の方程式ではなく、方程式の系です。粘性応力が対称速度勾配に比例するという構成法則とともに運動量保存をモデル化しています。非圧縮性は、体積が局所的に保存されるという拘束条件を追加します。

クレイ・ミレニアム問題および本サイトの大部分において、関連する設定は3次元非圧縮系です。

ナビエ–ストークス方程式の簡単な説明

最も簡単な一般的形式では、方程式は次のようになります：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

$$\nabla \cdot u = 0$$

ここで：

  • $u$ は流体の速度
  • $p$ は圧力
  • $\nu$ は粘性
  • $f$ は重力のような外力

左辺は速度が時間とともにどう変化するか、そして流体がどのように自身の運動を輸送するかを記述しています。右辺には流れを押し動かし滑らかにする力が含まれています。

$\mathbb{R}^3$ 上の非圧縮ナビエ–ストークス系は

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

$$\nabla \cdot u = 0, \qquad u(x,0)=u_0(x).$$

各項の標準的な解釈は以下の通りです：

  • $\partial_t u$：局所的な時間変化
  • $(u \cdot \nabla)u$：非線形移流、すなわち流体が自身の速度を輸送する
  • $-\nabla p$：圧力による力
  • $\nu \Delta u$：粘性拡散
  • $f$：外力
  • $\nabla \cdot u = 0$：非圧縮拘束条件

これはミレニアム問題、困難さ、証明戦略に関する本サイトのページ全体で使われている形式です。

ナビエ–ストークス方程式の由来

方程式は単純なアイデアから生まれます：ニュートンの第二法則 — 力は質量かける加速度に等しい — を流体の微小な塊に適用します。その塊には質量があり、力が作用します。圧力はあらゆる方向から押し、内部摩擦が引っ張り、重力が引き下げます。ニュートンの法則によれば、正味の力が塊の加速度を決定します。

剛体とは異なり、流体の塊は移動しながら変形できます。流れに乗って伸びたり歪んだりします。そのため「加速度」は単一の物体を追跡するほど単純ではありません。塊が流れに乗る様子を追う必要があります。

この収支 — すべての力は質量かけるこの複合的な加速度に等しい — を流体のあらゆる点で同時に書き下すことが、ナビエ–ストークス方程式の導出の出発点です。完全なステップバイステップの導出については、ナビエ–ストークス方程式の導出をご覧ください。

導出は、連続体に対するニュートンの第二法則である運動量保存から始まります。物質体積上での線形運動量の収支を書き、恒等式を局所化してPDEを得ます。

ニュートン流体の非圧縮流体では、コーシー応力テンソルは

$$T = -pI + 2\mu D(u), \qquad D(u)=\frac{1}{2}(\nabla u + \nabla u^T),$$

の形を持ちます。ここで $\mu$ は動粘性係数です。この構成法則を運動量方程式に代入し、密度で割ると、動粘性係数 $\nu = \mu/\rho$ を持つよく知られた非圧縮系が得られます。

非圧縮性は一定密度に対応し、発散なし条件 $\nabla \cdot u = 0$ を与えます。運動量収支とニュートン構成法則からの完全な導出については、導出をご覧ください。

ナビエ–ストークス方程式が難しい理由

困難な部分は非線形項 $(u \cdot \nabla)u$ です。流体は外力に応答するだけでなく、自分自身を押し動かします。このフィードバックが乱流やカオス的に見える運動を可能にします。

2次元空間では、方程式ははるかに良く振る舞います。3次元では、すべての滑らかな初期流れが永遠に滑らかであり続けるかどうかは、まだわかっていません。

これが、これらの方程式が工学を超えて有名な理由です：ナビエ–ストークス・ミレニアム問題に直結するのです。

主要な解析的困難は、自然なエネルギー評価が3次元方程式のスケーリングに対して弱すぎることです。大まかに言えば、標準的な $L^2$ 制御は非常に小さなスケールでの集中を排除するのに十分ではありません。

これが、大域弱解について知られていることと、大域的な滑らかさを証明するために必要なこととの間のギャップの源です。非線形移流項は、伝播したい正則性に対してエネルギー臨界です。

より詳しい議論については、なぜ難しいのかとアプローチをご覧ください。

何に使われているか

ナビエ–ストークス方程式は、科学と工学で毎日使われています。代表的な応用には：

  • 翼や車両周りの気流
  • 気象・気候モデル
  • 海洋循環
  • 産業用流体輸送
  • 血流やその他の生体輸送問題

実用上、人々は追加のモデリング仮定を含むこれらの方程式の数値近似を解いています。この実用的成功が、残された数学的問題をさらに際立たせています。

応用研究は通常、特定のレジーム下でのナビエ–ストークスまたは関連モデルの数値近似を使用します：非圧縮流、圧縮流、乱流クロージャ、境界層近似、縮約モデルなど。

直接数値シミュレーション、ラージエディシミュレーション、レイノルズ平均クロージャはすべて同じ連続体PDEフレームワークに遡りますが、3次元の基礎的な正則性問題を解消するものではありません。

実用的有効性と不完全な理論のこの乖離が、この主題を非常に魅力的にしている一因です。

次に読むべきもの

問題が解決されているかどうかが主な疑問であれば、ナビエ–ストークス問題は解決されたのか？から始めてください。

広範な数学的な意義を知りたい場合は、ミレニアム問題に進んでください。

ナビエ–ストークスが非粘性のオイラー方程式とどう異なるか、粘性がなぜ重要かを理解したい場合は、オイラー方程式 vs. ナビエ–ストークス方程式をお読みください。

乱流と小スケールの物理的直感を得たい場合は、レイノルズ数、乱流、そして小スケールが重要な理由をお読みください。

主な障害を知りたい場合は、なぜ難しいのかに進んでください。

本サイトでの自然な次のステップ：