ナビエ–ストークス方程式の導出

すべてのナビエ–ストークス方程式の導出は、同じアイデアから始まります：ニュートンの第二法則 — 力は質量かける加速度に等しい — を流体の微小な塊に適用するのです。

水、空気、その他の流体の小さな塊を取り出してください。それには質量があります。力が作用しています：圧力があらゆる方向から押し、内部摩擦が引っ張り、重力が引き下げます。ニュートンの法則によれば、正味の力がこの塊の加速度を決定します。

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

日常的な言葉で言えば：川の上のカヌーに座っていると想像してください。あなたの加速度は、あなたの場所での流れの時間変化だけでなく、流れがあなたをより速いあるいは遅い場所へ運んでいるという事実にも依存します。両方の効果が速度の変化に寄与します。

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

ナビエ–ストークス方程式の導出は、流体とともに移動する物質体積 $\Omega(t)$ にニュートンの第二法則の連続体力学的形式であるコーシーの運動量方程式から始まります。

密度 $\rho$、速度場 $u(x,t)$ を持つ流体に対し、線形運動量の保存は

$$\frac{d}{dt}\int_{\Omega(t)} \rho\, u\, dV = \int_{\partial\Omega(t)} T\, n\, dS + \int_{\Omega(t)} \rho\, f\, dV,$$

を述べます。ここで $T$ はコーシー応力テンソル、$n$ は外向き単位法線、$f$ は単位質量あたりの体積力（典型的には重力）です。

レイノルズ輸送定理で局所化すると、微分形

$$\rho\frac{Du}{Dt} = \nabla \cdot T + \rho f,$$

が得られます。ここで物質微分

$$\frac{Du}{Dt} = \partial_t u + (u \cdot \nabla)u$$

は、局所的な時間変化率と対流加速度の両方を捉えています。これは連続体のあらゆる点における $ma = F$ の厳密なアナログです。

ナビエ–ストークス方程式の導出には3種類の力が現れます：

1. 圧力。流体はあらゆる方向から塊を押します。片側の圧力がもう片側より高ければ、塊は低圧側に押されます。この不均衡は圧力勾配 $-\nabla p$ で捉えられます：高圧から低圧を指すベクトルであり、流体にどちらに加速するかを教えます。

2. 粘性力。異なる速度で動く隣接する流体層は互いに引き合います — 手を擦り合わせるようなものですが、流体全体で連続的に起きます。速い層は遅い隣を引っ張り、遅い層は速い層を引き留めます。この内部摩擦は速度差を滑らかにします。水のような単純な流体では、この摩擦の強さは一つの数値：粘性 $\mu$ で決まります。

3. 外力。外部から流体に作用するもの — 最も一般的には重力。これらは体積力です：表面だけでなく体積内のすべての流体に作用します。

ナビエ–ストークス方程式は、「質量×加速度 = 圧力 + 粘性力 + 外力」をあらゆる点で書いたものです。

コーシー応力テンソル $T$ はすべての内部接触力を符号化します。任意の流体に対して、等方的な圧力部分と偏差（粘性）部分に分解されます：

$$T = -pI + \tau,$$

ここで $p$ は力学的圧力（圧縮流に対しては $-\tfrac{1}{3}\mathrm{tr}\,T$ として定義）、$\tau$ は粘性応力テンソルです。

圧力勾配。$-pI$ に発散定理を適用すると、単位体積あたりの力 $-\nabla p$ が得られます。これは応力の発散の等方的部分です。

粘性応力。テンソル $\tau$ は、応力と変形速度テンソルを関連づける特定の構成法則に依存します。その発散 $\nabla \cdot \tau$ が単位体積あたりの粘性力を与えます。この段階では $\tau$ の形はまだ特定されていません — それは次のセクションのニュートン流体の仮定から来ます。

体積力。重力のような外力は単位体積あたり $\rho f$ として入ります。応力の分解を運動量方程式に代入すると

$$\rho\frac{Du}{Dt} = -\nabla p + \nabla \cdot \tau + \rho f.$$

これはニュートン流体であるかどうかを問わず、すべての単純流体に対する一般的な運動量方程式です。系を閉じるには、$\tau$ を特定する構成法則が必要です。

すべての流体が応力のもとで同じように振る舞うわけではありません。蜂蜜は水とは異なる仕方で運動に抵抗します。ケチャップは振ると流れやすくなります。水溶きコーンスターチは叩くと硬くなります。

ナビエ–ストークス方程式は特定の仮定を置きます：流体はニュートン流体であるということです。これは、内部摩擦（粘性応力）が流体の変形速度に直接比例することを意味します。変形速度を倍にすれば、応力も倍になります。構成法則は線形です：粘性応力は瞬間的な変形速度に比例します。

水、空気、多くの一般的な流体はニュートン流体として非常によくモデル化されます。しかしこれは真にニュートンの法則からの帰結ではなく仮定です。導出にはこれが必要です。これがなければ、異なるクラスの方程式（非ニュートン流体モデル）が得られます。

比例定数は動粘性率 $\mu$：流体がせん断に抵抗する程度を測る材料定数です。水は低い粘性を持ち、蜂蜜は高い粘性を持ちます。

第二の粘性パラメータ $\lambda$ もあり、第二粘性係数と呼ばれることもあります。これは流体が圧縮または膨張するときに重要になります。さらなる一般的な簡略化 — ストークスの仮説 — は $\lambda = -\frac{2}{3}\mu$ と設定します。これは定理ではなく追加の仮定です。

ニュートン流体の構成法則は、粘性応力 $\tau$ がひずみ速度テンソル $D(u) = \tfrac{1}{2}(\nabla u + \nabla u^T)$ の線形等方関数であることを仮定します。等方テンソル関数の表現定理により、最も一般的なそのような法則は

$$\tau = 2\mu\, D(u) + \lambda (\nabla \cdot u)\, I,$$

等価的に

$$\tau = \mu(\nabla u + \nabla u^T) + \lambda(\nabla \cdot u)\, I,$$

と書けます。ここで $\mu > 0$ は動粘性率（せん断粘性）、$\lambda$ は第二粘性係数です。

これがニュートン流体の定義的仮定です。第一原理から導出されるものではなく、多くの一般的な流体に対して経験的に検証された構成法則的仮説です。

ストークスの仮説はさらに $\lambda = -\frac{2}{3}\mu$ を仮定し、体積粘性 $\kappa = \lambda + \frac{2}{3}\mu$ をゼロにします。この簡略化は広く用いられていますが、独立した仮定であり — 熱力学やニュートン流体の仮説自体の帰結ではありません。古典的な運動論は理想化された仮定のもとで単原子理想気体に対してゼロの体積粘性を予測しますが、多くの実際の気体や液体に対してストークスの仮説は近似にすぎません。

クラウジウス–デュエムの不等式からの熱力学的拘束は $\mu \geq 0$ と $3\lambda + 2\mu \geq 0$（等価的に $\kappa \geq 0$）のみを要求します。

ここでパーツを組み合わせます。以下があります：

• 左辺に質量かける加速度（速度の物質微分）
• 右辺に圧力、粘性摩擦、重力
• 摩擦と変形速度を結びつけるニュートン流体の規則

代入して整理すると、圧縮性ナビエ–ストークス運動量方程式が得られます：

$$\rho\Big(\frac{\partial u}{\partial t} + (u \cdot \nabla) u\Big) = -\nabla p + \mu\, \Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u) + \rho\, f$$

各項には物理的意味があります：

• $\rho\,\partial_t u$ — 固定点での速度の時間変化
• $\rho(u \cdot \nabla) u$ — 流体が自身の速度を場所から場所へ運ぶこと（移流）
• $-\nabla p$ — 高圧から低圧への押し
• $\mu\,\Delta u$ — 粘性が速度差を滑らかにする
• $(\mu + \lambda)\nabla(\nabla \cdot u)$ — 流体が圧縮または膨張するときにのみ重要な追加の粘性項
• $\rho f$ — 重力のような外力

これが完全な運動量方程式です。質量保存と状態方程式（圧力を密度に結びつける）とともに、圧縮性ナビエ–ストークス系を形成します。

ニュートン構成法則 $\tau = 2\mu D(u) + \lambda(\nabla \cdot u)I$ を運動量方程式 $\rho \frac{Du}{Dt} = -\nabla p + \nabla \cdot \tau + \rho f$ に代入し、$\tau$ の発散を計算すると：

$$\nabla \cdot \tau = \nabla \cdot \big[\mu(\nabla u + \nabla u^T)\big] + \nabla\big[\lambda(\nabla \cdot u)\big].$$

$\mu$ と $\lambda$ が定数であれば（多くの導出における標準的な仮定）、これは

$$\nabla \cdot \tau = \mu\,\Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u),$$

と簡略化されます。ここでベクトル恒等式 $\nabla \cdot (\nabla u^T) = \nabla(\nabla \cdot u)$ を用いました。圧縮性ナビエ–ストークス運動量方程式は

$$\rho\big(\partial_t u + (u \cdot \nabla)u\big) = -\nabla p + \mu\,\Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u) + \rho f.$$

これは質量保存（連続の方程式）

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

と状態方程式 $p = p(\rho, \theta)$（またはエネルギー方程式）と結合して系を閉じる必要があります。ストークスの仮説 $\lambda = -\frac{2}{3}\mu$ では、係数 $\mu + \lambda = \frac{1}{3}\mu$ となります。非圧縮系との詳細な比較は、非圧縮 vs. 圧縮ナビエ–ストークスをご覧ください。

多くの流れ — 管内の水、遅い空気の流れ、海洋循環 — は密度が実質的に一定である流体を含みます。この簡略化を行うと、一般的な方程式はより明快で扱いやすい系に変換されます。

一定密度は $\rho$ がどこでもいつでも同じであることを意味します。すると質量保存から、速度場は $\nabla \cdot u = 0$ を満たさなければなりません：流体はどこでも圧縮も膨張もしません。これが非圧縮性拘束条件です。

$\nabla \cdot u = 0$ により、追加の粘性項 $(\mu + \lambda)\nabla(\nabla \cdot u)$ は完全に消えます。第二粘性係数 $\lambda$ は脱落します — 流体が圧縮できないとき、その役割はありません。運動量方程式は次のようになります：

$$\rho\Big(\frac{\partial u}{\partial t} + (u \cdot \nabla)u\Big) = -\nabla p + \mu\,\Delta u + \rho f$$

両辺を $\rho$ で割り、$\nu = \mu / \rho$（動粘性係数）と書き、$1/\rho$ の因子を吸収するように圧力を再定義すると、標準形が得られます：

$$\frac{\partial u}{\partial t} + (u \cdot \nabla)u = -\nabla p + \nu\,\Delta u + f$$

$$\nabla \cdot u = 0$$

これがクレイ・ミレニアム問題で研究されている系です。本サイト全体、および多くの数学的なナビエ–ストークスの取り扱いで見かけるバージョンです。この形式と圧縮系の詳細な比較については、非圧縮 vs. 圧縮ナビエ–ストークスをご覧ください。

$\nu = 0$（粘性が全くない場合）とするとオイラー方程式が得られることに注意してください — 関連はありますが重要な意味で異なる系です。

一般的な非圧縮の特殊化は、流れ全体で一定密度 $\rho$ を仮定します。連続の方程式 $\partial_t \rho + \nabla \cdot(\rho u) = 0$ から $\nabla \cdot u = 0$ が強制されます。より一般的には、非圧縮性は $\nabla \cdot u = 0$（連続の式を通じて等価的に $D\rho/Dt = 0$）で符号化され、原理的には可変密度も許容しますが、一定密度の場合がクレイ問題の標準的な設定です。

非圧縮性のもとで、項 $(\mu + \lambda)\nabla(\nabla \cdot u)$ は恒等的に消えます。第二粘性係数 $\lambda$ は無関係になります — ストークスの仮説もその他の $\lambda$ に関する仮定も非圧縮系には影響しません。運動量方程式は

$$\rho\big(\partial_t u + (u \cdot \nabla)u\big) = -\nabla p + \mu\,\Delta u + \rho f.$$

$\rho$ で割り、動粘性係数 $\nu = \mu/\rho$ を定義すると（$p$ は $1/\rho$ の因子を吸収するように再定義）、$\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).$$

これがクレイ・ミレニアム問題（Fefferman, 2000）で研究されている系です。公式の問題は $\mathbb{R}^3$（または $\mathbb{T}^3$）上の滑らかで発散なしの初期データに対する大域正則性に関するもので、クレイの標準的定式化（Fefferman, 2000）に従います。

圧力 $p$ は独立した動的変数ではありません：運動量方程式の発散をとり $\nabla \cdot u = 0$ を用いることで決定され（定数を除いて）、ポアソン方程式 $\Delta p = -\nabla \cdot [(u \cdot \nabla)u] + \nabla \cdot f$（または $f$ が発散なしまたは不在のとき $\Delta p = -\partial_i \partial_j(u_i u_j)$）を得ます。この楕円型結合は非圧縮系の特徴的な特徴です。

$\nu = 0$ とすると非圧縮オイラー方程式が回復されます。これら二つの系の関係は、数学的流体力学の多くの未解決問題の中心にあります。

上の導出はナビエ–ストークス方程式を与えます。方程式が何であるか、そしてなぜその形を取るかを教えてくれます。すべての項は物理的原理または明示的な仮定に遡ります。

しかし方程式を導出することと、その解を理解することは同じではありません。導出は以下に答えません：

• 解はすべての時刻で常に存在するか？
• 滑らかに始まれば、滑らかであり続けるか？
• 速度は有限時間で無限大に爆発しうるか？

これらは方程式の物理的起源ではなく、数学的振る舞いに関する問いです。そして3次元では、未解決のままです。このギャップがクレイ・ミレニアム問題につながります：滑らかな3次元非圧縮ナビエ–ストークスのデータが常に大域滑らかな解を生成するか、それとも有限時間で特異点が形成されうるか。

クレイ数学研究所はその解決に百万ドルを提示しています。175年以上の歳月を経て、問題は立ち続けています。

導出はナビエ–ストークス方程式を、運動量保存とニュートン構成法則に基づく根拠のあるPDE系として確立します。しかし数学的適切性の中心的な問いには触れません。

具体的に、導出は以下について沈黙しています：

• 大域存在：滑らかな初期データに対して $t > 0$ のすべてで滑らかな解が持続するかどうか。
• 正則性：解が $C^\infty(\mathbb{R}^3 \times [0,\infty))$ に留まるか、特異点を発展させうるか。
• 一意性：さまざまな関数空間の設定での解が一意であるかどうか。

ルレイの理論（1934）は $L^2$ における弱解の大域存在を保証しますが、ルレイ–ホップ解の一意性と正則性は3次元で未証明のままです。利用可能なエネルギー評価と非線形性のスケーリングとのギャップが、困難の解析的核心です。

クレイ・ミレニアム問題は、$\mathbb{R}^3$ における大域正則性の証明または反証を正確に求めています。導出はPDE系を動機づけますが、3次元における存在、正則性、一意性を解決するものではありません。

方程式の由来がわかったところで、次のステップはこちらです：

• ナビエ–ストークス方程式とは何か？ — 方程式そのものと項ごとの説明
• 非圧縮 vs. 圧縮ナビエ–ストークス — 密度が変化するとき何が変わるか、非圧縮の場合がなぜ特別か
• オイラー方程式 vs. ナビエ–ストークス方程式 — 粘性を完全に取り除くとどうなるか
• ミレニアム問題 — これらの方程式を有名にした未解決問題

さらなる研究のための関連ページ：

• ナビエ–ストークス方程式とは何か？ — 各項の詳細な議論を含む標準形
• 非圧縮 vs. 圧縮ナビエ–ストークス — 完全な圧縮系、順圧クロージャ、低マッハ数極限
• オイラー方程式 vs. ナビエ–ストークス方程式 — 非粘性極限、$\nu$ の役割、粘性消失問題
• ナビエ–ストークスの存在と滑らかさ — クレイの正確な定式化、既知の結果、未解決問題