流体力学入門：水や空気の振る舞いを支配する法則を理解する

はじめに

私たちの身の回りには、水や空気といった「流体」が満ちています。蛇口から流れる水、窓から入ってくる風、飛行機が空を飛ぶ仕組み、血液が体内を循環するメカニズム──これらすべてを理解するための学問が流体力学です。

流体力学は、機械工学、航空宇宙工学、土木工学、化学工学、生体工学など、工学のあらゆる分野で必要とされる基礎学問です。また、気象学、海洋学、天体物理学といった自然科学の分野でも不可欠な知識となっています。

この記事では、流体力学の基礎から応用まで、理系の学生が押さえておくべき内容を体系的に解説していきます。数式の導出過程も丁寧に追いながら、物理的な意味を理解することを目指します。

---

第1章：流体とは何か

1.1 流体の定義

流体（fluid） とは、せん断応力を受けると連続的に変形し続ける物質のことです。これは固体との本質的な違いです。

固体にせん断力を加えると、ある程度変形した後に静止します（弾性変形の場合）。しかし流体は、どんなに小さなせん断応力であっても、それが作用し続ける限り変形を続けます。

より厳密に言えば、流体は以下のように定義されます。

流体の定義
静止状態においてせん断応力に抵抗できない物質。言い換えれば、静止流体内部ではせん断応力は常にゼロである。

この定義から直ちに導かれる重要な帰結があります。静止流体内の任意の点において、応力は等方的（すべての方向に同じ大きさ）でなければなりません。これが圧力という概念につながります。

1.2 連続体仮説

流体は実際には分子の集合体です。水1cm³には約3.3×10²²個もの分子が含まれています。しかし、これらの分子一つ一つの運動を追跡することは現実的ではありません。

そこで流体力学では連続体仮説（continuum hypothesis） を採用します。

連続体仮説
流体を、その物理量（密度、速度、圧力など）が空間的に連続的に分布する連続体として扱う。

この仮説が成り立つためには、考える現象の特性長さスケールLが、分子の平均自由行程λに比べて十分大きい必要があります。この比をクヌーセン数（Knudsen number） と呼びます。

$$\text{Kn} = \frac{\lambda}{L}$$

通常の流体力学が適用できるのは、Kn ≪ 1 の場合です。

• Kn < 0.01：連続体として扱える（通常の流体力学）
• 0.01 < Kn < 0.1：希薄気体効果が現れ始める
• Kn > 0.1：分子運動論的な取り扱いが必要

標準大気圧・常温の空気では、λ ≈ 0.07 μm です。したがって、ミリメートルスケール以上の現象であれば、連続体仮説は十分に成り立ちます。

1.3 流体の分類

圧縮性による分類

非圧縮性流体（incompressible fluid）：密度が一定とみなせる流体。液体は一般にこの近似が良く成り立ちます。また、気体でも流速がマッハ数 Ma < 0.3 程度であれば、非圧縮性近似が使えます。

圧縮性流体（compressible fluid）：密度変化を考慮する必要がある流体。高速気流、音波の伝播、衝撃波などを扱う場合に重要です。

粘性による分類

粘性流体（viscous fluid）：内部摩擦（粘性）を持つ実在の流体。

完全流体/理想流体（perfect/ideal fluid）：粘性がゼロと仮定した仮想的な流体。数学的取り扱いが容易なため、流体力学の基礎理論では重要です。

応力と変形速度の関係による分類

ニュートン流体（Newtonian fluid）：せん断応力がせん断変形速度に比例する流体。水、空気、多くの油などが該当します。

$$\tau = \mu \frac{du}{dy}$$

ここで、τはせん断応力、μは粘度（粘性係数）、du/dyは速度勾配です。

非ニュートン流体（non-Newtonian fluid）：上記の線形関係が成り立たない流体。血液、ペンキ、マヨネーズ、高分子溶液などが該当します。

---

第2章：流体静力学

2.1 静止流体中の圧力

静止している流体内部では、せん断応力は存在せず、垂直応力（圧力）のみが作用します。まず、静止流体中の圧力の基本的性質を確認しましょう。

パスカルの原理

パスカルの原理
静止流体中の任意の点において、圧力はすべての方向に等しく作用する（圧力は等方的なスカラー量である）。

これを証明するために、静止流体中の微小な直角三角形要素を考えます。

        /|
       / |
      /  | dy
     /   |
    /θ   |
   /_____|
     dx

斜面の法線方向を n、その面の圧力を p_n、水平面・垂直面の圧力をそれぞれ p_x、p_y とします。斜面の面積を dA とすると、

• 斜面：$p_n \cdot dA$（面に垂直な方向）
• 水平面：$p_x \cdot dA \sin\theta$
• 垂直面：$p_y \cdot dA \cos\theta$

体積力（重力）は体積に比例するため、微小要素では面に働く力に比べて高次の微小量となり無視できます。

x方向の力の釣り合い：

$$p_n dA \sin\theta = p_x dA \sin\theta$$

y方向の力の釣り合い：

$$p_n dA \cos\theta = p_y dA \cos\theta$$

したがって、$p_n = p_x = p_y$ となり、圧力は方向によらないことが示されました。

2.2 圧力の基本方程式

静止流体中で、重力の作用下における圧力分布を求めましょう。鉛直上向きにz軸をとります。

微小直方体要素（体積 dx·dy·dz）に作用する力を考えます。

z方向の力の釣り合い：

$$p \cdot dx \cdot dy - \left(p + \frac{\partial p}{\partial z}dz\right) dx \cdot dy - \rho g \cdot dx \cdot dy \cdot dz = 0$$

整理すると、

$$\frac{\partial p}{\partial z} = -\rho g$$

同様に、x方向、y方向についても、

$$\frac{\partial p}{\partial x} = 0, \quad \frac{\partial p}{\partial y} = 0$$

これをベクトル形式でまとめると、

静止流体の基本方程式

$$\nabla p = -\rho \mathbf{g} = -\rho g \mathbf{k}$$

または、重力方向（鉛直下向き）に深さ h をとると、

$$\frac{dp}{dh} = \rho g$$

非圧縮性流体の場合

密度ρが一定の場合、上式を積分すると、

$$p = p_0 + \rho g h$$

ここで、p₀は基準点（通常は自由表面）の圧力、hは基準点からの深さです。

これは有名な水圧の公式です。深さ10mごとに約1気圧（≈ 0.1 MPa）ずつ圧力が増加します。

圧縮性流体（等温大気）の場合

気体では密度が圧力に依存します。理想気体の状態方程式 $p = \rho R T$ を用いると、等温（T = 一定）の場合、

$$\frac{dp}{dz} = -\rho g = -\frac{pg}{RT}$$

これを解くと、

$$p = p_0 \exp\left(-\frac{gz}{RT}\right) = p_0 \exp\left(-\frac{z}{H}\right)$$

ここで、$H = RT/g$ はスケールハイトと呼ばれ、圧力が1/eになる高度を表します。地球大気（T ≈ 288 K）では H ≈ 8.4 km です。

2.3 マノメーター（圧力計）

圧力の測定には、静止流体中の圧力分布を利用したマノメーターが用いられます。

U字管マノメーター

U字管に測定流体と基準液体（通常は水または水銀）を入れ、液面の高さの差から圧力を求めます。

    ┌───┐        ┌───┐
    │   │        │   │
    │ A │        │ B │
    │   │        │   │
    │   └────────┘   │
    │       ▲        │
    │       h        │
    │       ▼        │
    └────────────────┘
         水銀

点Aと点Bの圧力差は、

$$p_A - p_B = \rho_{Hg} g h - \rho_{fluid} g h$$

測定流体が気体の場合、$\rho_{fluid} \ll \rho_{Hg}$ なので、

$$p_A - p_B \approx \rho_{Hg} g h$$

2.4 浮力とアルキメデスの原理

流体中に沈められた物体には、その物体が排除した流体の重さに等しい上向きの力（浮力）が作用します。

アルキメデスの原理

$$F_b = \rho_{fluid} V_{submerged} g$$

F_b：浮力、V_submerged：物体の流体中に沈んでいる部分の体積

この原理の証明は、物体表面に作用する圧力を積分することで得られます。物体表面の微小面積 dA に作用する圧力 p による力を、物体全体で積分すると、

$$\mathbf{F}_b = -\oint_S p \mathbf{n} dA$$

ガウスの発散定理を用いると、

$$\mathbf{F}_b = -\int_V \nabla p dV = \int_V \rho_{fluid} \mathbf{g} dV = \rho_{fluid} V \mathbf{g}$$

浮体の安定性

浮遊物体が平衡位置から傾いたとき、元の位置に戻ろうとするか（安定）、さらに傾こうとするか（不安定）は、メタセンタの位置で決まります。

• 浮心（B）：物体が排除した流体の体積の重心
• 重心（G）：物体自体の重心
• メタセンタ（M）：物体が微小角度傾いたときの、新しい浮力の作用線と元の対称軸との交点

メタセンタ高さ GM は、

$$GM = \frac{I}{V} - BG$$

ここで、I は水面での物体断面の断面二次モーメント、V は排除流体の体積、BG は浮心と重心の距離です。

• GM > 0：安定（傾いても元に戻る）
• GM < 0：不安定（傾きが増大する）

船舶設計では、この安定性の確保が極めて重要です。

---

第3章：流体運動学

3.1 流れの記述方法

流体の運動を記述する方法には、大きく分けて2つのアプローチがあります。

ラグランジュ記述（Lagrangian description）

個々の流体粒子を追跡し、その位置・速度・加速度を時間の関数として記述する方法です。

$$\mathbf{x} = \mathbf{x}(\mathbf{x}_0, t)$$

ここで、$\mathbf{x}_0$は時刻 t = 0 における粒子の初期位置（ラグランジュ座標）、$\mathbf{x}$は時刻 t における位置です。

速度は、

$$\mathbf{v} = \frac{d\mathbf{x}}{dt}\bigg|_{\mathbf{x}_0}$$

この方法は粒子追跡に適しており、質点力学の延長として理解しやすいですが、流体全体の運動を把握するには不向きな場合があります。

オイラー記述（Eulerian description）

空間に固定した点で、その点を通過する流体の速度・圧力などの物理量を観測する方法です。

$$\mathbf{v} = \mathbf{v}(\mathbf{x}, t)$$

流体力学では、主にこのオイラー記述が用いられます。実験でも、特定の位置に計測器を設置して測定することが多いため、実用的です。

物質微分（material derivative）

オイラー記述において、流体粒子に乗って移動しながら観測する変化率を表すのが物質微分（または実質微分、ラグランジュ微分）です。

任意のスカラー量 φ について、

物質微分

$$\frac{D\phi}{Dt} = \frac{\partial \phi}{\partial t} + (\mathbf{v} \cdot \nabla)\phi$$

• 第1項：局所的時間変化率（その場所での時間変化）
• 第2項：対流変化率（流れによる移動に伴う変化）

速度についても同様に、

$$\frac{D\mathbf{v}}{Dt} = \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v}$$

これが流体粒子の加速度を表します。

具体例

温度 T(x, t) を考えます。川を流れる水の温度を測定する場合：

1. 岸に立って温度計を川に入れておくと、その場所での $\partial T/\partial t$ が測定されます
2. ボートに乗って流れと一緒に移動しながら測定すると、$DT/Dt$ が測定されます

仮に川の上流は冷たく下流は暖かいとします。流れが定常（$\partial T/\partial t = 0$）でも、ボートに乗っていると温度上昇を感じます。これが対流項の寄与です。

3.2 流れの可視化

流線（streamline）

流線は、各点で速度ベクトルに接する曲線です。瞬間の流れの様子を表します。

流線の微分方程式は、

$$\frac{dx}{v_x} = \frac{dy}{v_y} = \frac{dz}{v_z}$$

または、

$$\mathbf{v} \times d\mathbf{s} = 0$$

流跡線（pathline）

流跡線は、一つの流体粒子が時間とともに描く軌跡です。

$$\frac{d\mathbf{x}}{dt} = \mathbf{v}(\mathbf{x}, t)$$

を積分して得られます。

流脈線（streakline）

流脈線は、ある瞬間に、過去のある時刻からその瞬間までに特定の点を通過したすべての粒子を結んだ線です。煙突からの煙や、水面に落とした染料が描く線がこれに相当します。

重要な関係
定常流（$\partial \mathbf{v}/\partial t = 0$）では、流線、流跡線、流脈線はすべて一致します。非定常流では、一般にこれらは異なります。

3.3 変形と回転

流体要素の運動は、以下の4つの基本的な運動の重ね合わせとして理解できます。

1. 並進：全体が平行移動
2. 回転：全体が剛体的に回転
3. 線形変形（伸縮）：体積が変化
4. せん断変形：体積を保ったまま形状が変化

速度勾配テンソル

速度場 $\mathbf{v}(\mathbf{x})$ の空間微分を表すテンソルを考えます。

$$L_{ij} = \frac{\partial v_i}{\partial x_j}$$

このテンソルは対称部分と反対称部分に分解できます。

$$L_{ij} = D_{ij} + \Omega_{ij}$$

変形速度テンソル（対称部分）：

$$D_{ij} = \frac{1}{2}\left(\frac{\partial v_i}{\partial x_j} + \frac{\partial v_j}{\partial x_i}\right)$$

回転速度テンソル（反対称部分）：

$$\Omega_{ij} = \frac{1}{2}\left(\frac{\partial v_i}{\partial x_j} - \frac{\partial v_j}{\partial x_i}\right)$$

渦度（vorticity）

回転速度テンソルの情報は、渦度ベクトルに集約されます。

渦度の定義

$$\boldsymbol{\omega} = \nabla \times \mathbf{v} = \text{curl } \mathbf{v}$$

渦度は流体要素の局所的な角速度の2倍に等しく、

$$\boldsymbol{\omega} = 2\boldsymbol{\Omega}$$

ここで $\boldsymbol{\Omega}$ は流体要素の角速度ベクトルです。

渦なし流れ（irrotational flow）：$\boldsymbol{\omega} = 0$ となる流れ。この場合、速度はスカラーポテンシャルの勾配として表せます。

$$\mathbf{v} = \nabla \phi$$

このφを速度ポテンシャルと呼びます。

---

第4章：流体力学の基礎方程式

4.1 質量保存則（連続の式）

積分形

検査体積 V（空間に固定）を考えます。質量保存則は、「検査体積内の質量の時間変化率 = 検査表面を通って流入する質量流量」と表されます。

$$\frac{d}{dt}\int_V \rho dV = -\oint_S \rho \mathbf{v} \cdot \mathbf{n} dA$$

ここで S は検査体積の表面、$\mathbf{n}$ は外向き単位法線ベクトルです。

微分形

ガウスの発散定理を用いて、上式を微分形に変換します。

$$\int_V \frac{\partial \rho}{\partial t} dV = -\int_V \nabla \cdot (\rho \mathbf{v}) dV$$

これが任意の検査体積について成り立つので、

連続の式（質量保存則）

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$

または、

$$\frac{D\rho}{Dt} + \rho \nabla \cdot \mathbf{v} = 0$$

非圧縮性流体の場合

密度が一定（$D\rho/Dt = 0$）の場合、連続の式は単純化されます。

$$\nabla \cdot \mathbf{v} = 0$$

これは、非圧縮性流体では速度場がソレノイダル（発散なし）であることを意味します。

直交座標では、

$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$$

4.2 運動量保存則（ナビエ・ストークス方程式の導出）

ニュートンの第2法則の適用

流体粒子に対してニュートンの第2法則を適用します。

$$\rho \frac{D\mathbf{v}}{Dt} = \text{（単位体積あたりの力）}$$

流体に作用する力は、体積力と表面力に分類されます。

体積力（body force）：重力、電磁力など、流体の体積に比例して作用する力

$$\mathbf{f}_b = \rho \mathbf{g}$$

表面力（surface force）：圧力、粘性応力など、流体要素の表面に作用する力

応力テンソル

流体中の任意の点における応力状態は、応力テンソル $\boldsymbol{\sigma}$ で記述されます。

$$\sigma_{ij} = \begin{pmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{pmatrix}$$

対角成分は垂直応力、非対角成分はせん断応力です。角運動量保存則から、応力テンソルは対称（$\sigma_{ij} = \sigma_{ji}$）です。

流体では、応力テンソルを圧力部分と粘性応力部分に分解するのが便利です。

$$\sigma_{ij} = -p\delta_{ij} + \tau_{ij}$$

ここで、$\delta_{ij}$はクロネッカーのデルタ、$\tau_{ij}$は偏差応力テンソル（粘性応力テンソル）です。

ニュートン流体の構成方程式

ニュートン流体では、粘性応力は変形速度に比例します。

$$\tau_{ij} = 2\mu D_{ij} + \lambda (\nabla \cdot \mathbf{v}) \delta_{ij}$$

ここで、μは動粘度（せん断粘度）、λは第二粘度（体積粘度）です。

非圧縮性流体（$\nabla \cdot \mathbf{v} = 0$）では、第二項は消えます。

ナビエ・ストークス方程式

以上をまとめると、運動量保存則は次のようになります。

ナビエ・ストークス方程式（非圧縮性ニュートン流体）

$$\rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \rho \mathbf{g}$$

成分形（直交座標、重力は-z方向）：

$$\rho \left( \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z} \right) = -\frac{\partial p}{\partial x} + \mu \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)$$

（v, w成分も同様）

各項の物理的意味：

1. 慣性項（左辺）：流体粒子の加速度に対応

   • $\partial \mathbf{v}/\partial t$：非定常項（時間変化）
   • $(\mathbf{v} \cdot \nabla)\mathbf{v}$：対流項（移流による加速）

2. 圧力項：$-\nabla p$

   • 圧力勾配による力

3. 粘性項：$\mu \nabla^2 \mathbf{v}$

   • 粘性による運動量の拡散

4. 体積力項：$\rho \mathbf{g}$

   • 重力などの外力

オイラー方程式

粘性を無視した理想流体の運動方程式をオイラー方程式と呼びます。

$$\rho \frac{D\mathbf{v}}{Dt} = -\nabla p + \rho \mathbf{g}$$

高レイノルズ数の流れで、壁面から離れた領域では、この近似が有効です。

4.3 エネルギー保存則

完全な流体力学系を記述するには、質量・運動量に加えてエネルギーの保存則も必要です。特に圧縮性流体では必須です。

$$\rho \frac{D}{Dt}\left( e + \frac{v^2}{2} \right) = -\nabla \cdot (p\mathbf{v}) + \nabla \cdot (\boldsymbol{\tau} \cdot \mathbf{v}) + \nabla \cdot (k \nabla T) + \rho \mathbf{g} \cdot \mathbf{v} + Q$$

ここで、e は単位質量あたりの内部エネルギー、k は熱伝導率、Q は体積熱源です。

4.4 境界条件

偏微分方程式を解くには、適切な境界条件が必要です。

壁面での境界条件

滑りなし条件（no-slip condition）：粘性流体では、固体壁面で流体の速度が壁の速度と一致します。

$$\mathbf{v}_{fluid} = \mathbf{v}_{wall}$$

静止壁では $\mathbf{v} = 0$ となります。これは実験的に確立された条件です。

滑り条件（slip condition）：理想流体では、壁面に垂直な速度成分のみがゼロ。

$$\mathbf{v} \cdot \mathbf{n} = 0$$

自由表面での境界条件

• 運動学的条件：流体粒子が表面上にとどまる
• 力学的条件：表面での応力の連続性

---

第5章：無次元数と相似則

5.1 無次元化の重要性

流体力学では、問題を無次元化することで、多くの重要な知見が得られます。

• 現象を支配するパラメータを明確化
• 異なるスケールの現象の比較（相似則）
• 数値計算の効率化
• 実験結果の一般化

5.2 ナビエ・ストークス方程式の無次元化

代表長さ L、代表速度 U、代表時間 L/U、代表圧力 ρU² を用いて無次元変数を定義します。

$$\mathbf{x}^* = \frac{\mathbf{x}}{L}, \quad \mathbf{v}^* = \frac{\mathbf{v}}{U}, \quad t^* = \frac{tU}{L}, \quad p^* = \frac{p}{\rho U^2}$$

非圧縮性ナビエ・ストークス方程式を無次元化すると、

$$\frac{\partial \mathbf{v}^*}{\partial t^*} + (\mathbf{v}^* \cdot \nabla^*)\mathbf{v}^* = -\nabla^* p^* + \frac{1}{Re} \nabla^{*2} \mathbf{v}^* + \frac{1}{Fr^2} \mathbf{e}_g$$

ここで現れた無次元数がレイノルズ数とフルード数です。

5.3 主要な無次元数

レイノルズ数（Reynolds number）

レイノルズ数

$$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$$

• 慣性力と粘性力の比
• ν = μ/ρ：動粘度

レイノルズ数は流体力学で最も重要な無次元数です。

• Re ≪ 1：粘性力支配（ストークス流れ、クリープ流れ）
• Re ≫ 1：慣性力支配（しかし境界層内は粘性重要）
• Re が臨界値を超えると層流から乱流へ遷移

流れの種類	臨界レイノルズ数
円管内流れ	Re_cr ≈ 2300
平板上の境界層	Re_cr ≈ 5×10⁵
円柱まわりの流れ	Re_cr ≈ 3×10⁵

フルード数（Froude number）

$$Fr = \frac{U}{\sqrt{gL}}$$

慣性力と重力の比を表します。自由表面を持つ流れ（開水路流れ、船の造波抵抗）で重要です。

• Fr < 1：常流（亜臨界流）
• Fr = 1：臨界流
• Fr > 1：射流（超臨界流）

マッハ数（Mach number）

$$Ma = \frac{U}{c}$$

流速と音速の比。圧縮性の指標です。

• Ma < 0.3：非圧縮性近似が有効
• 0.3 < Ma < 0.8：亜音速流
• 0.8 < Ma < 1.2：遷音速流
• 1.2 < Ma < 5：超音速流
• Ma > 5：極超音速流

その他の重要な無次元数

ストローハル数（Strouhal number）：

$$St = \frac{fL}{U}$$

周期的な流れ現象の特性。カルマン渦列では St ≈ 0.2。

オイラー数（Euler number）：

$$Eu = \frac{\Delta p}{\rho U^2}$$

圧力降下の無次元表示。

ウェーバー数（Weber number）：

$$We = \frac{\rho U^2 L}{\sigma}$$

慣性力と表面張力の比。液滴、気泡の現象で重要。

グラスホフ数（Grashof number）：

$$Gr = \frac{g\beta \Delta T L^3}{\nu^2}$$

自然対流における浮力と粘性力の比。

プラントル数（Prandtl number）：

$$Pr = \frac{\nu}{\alpha} = \frac{c_p \mu}{k}$$

運動量拡散と熱拡散の比。流体の物性のみに依存。

5.4 相似則と模型実験

二つの流れが力学的に相似であるためには、

1. 幾何学的相似：形状が相似
2. 運動学的相似：速度場が相似
3. 力学的相似：支配的な無次元数が等しい

例えば、高レイノルズ数の流れを風洞で再現する場合、実機と模型でレイノルズ数を一致させる必要があります。

$$Re_{model} = Re_{prototype}$$

$$\frac{U_m L_m}{\nu_m} = \frac{U_p L_p}{\nu_p}$$

模型が1/10スケールで同じ流体を使う場合、流速を10倍にする必要があります。実際には風洞の能力に限界があるため、圧縮空気や低温空気を使って動粘度を変えることもあります。

---

第6章：理想流体の流れ

6.1 ベルヌーイの定理

定常・非粘性流れに対するベルヌーイの定理

オイラー方程式から出発して、ベルヌーイの定理を導出します。

定常（$\partial/\partial t = 0$）のオイラー方程式：

$$(\mathbf{v} \cdot \nabla)\mathbf{v} = -\frac{1}{\rho}\nabla p + \mathbf{g}$$

ベクトル恒等式 $(\mathbf{v} \cdot \nabla)\mathbf{v} = \nabla(v^2/2) - \mathbf{v} \times (\nabla \times \mathbf{v})$ を用いると、

$$\nabla\left(\frac{v^2}{2}\right) - \mathbf{v} \times \boldsymbol{\omega} = -\frac{1}{\rho}\nabla p - \nabla(gz)$$

ここで、重力 $\mathbf{g} = -g\nabla z$ としました。

この式を流線に沿って積分すると、$\mathbf{v} \times \boldsymbol{\omega}$ は流線に垂直なので積分に寄与せず、

ベルヌーイの定理（定常・非粘性・非圧縮性流れ）

$$\frac{v^2}{2} + \frac{p}{\rho} + gz = \text{const. along a streamline}$$

または、

$$\frac{1}{2}\rho v^2 + p + \rho gz = \text{const.}$$

各項の物理的意味（単位体積あたり）：

• $\frac{1}{2}\rho v^2$：動圧（運動エネルギー）
• $p$：静圧（圧力エネルギー）
• $\rho gz$：位置圧（位置エネルギー）

これらの和を全圧または淀み圧と呼びます。

渦なし流れでのベルヌーイの定理

流れ全体で渦度がゼロ（$\boldsymbol{\omega} = 0$）の場合、ベルヌーイの定理は流線に限らず、流れ場全体で成り立ちます。

$$\frac{v^2}{2} + \frac{p}{\rho} + gz = \text{const. throughout the flow}$$

ベルヌーイの定理の応用例

ピトー管：流速の測定に使用。淀み点（v = 0）と流れの中の点を比較。

$$p_{stagnation} = p_{static} + \frac{1}{2}\rho v^2$$

$$v = \sqrt{\frac{2(p_{stag} - p_{static})}{\rho}}$$

ベンチュリ管：流量計。断面積の変化に伴う圧力変化を測定。

$$v_1 A_1 = v_2 A_2$$

$$p_1 + \frac{1}{2}\rho v_1^2 = p_2 + \frac{1}{2}\rho v_2^2$$

これらから流量 Q = v₁A₁ を求められます。

トリチェリの定理：容器からの流出速度。

$$v = \sqrt{2gh}$$

ここで h は水面から流出口までの高さ。

6.2 ポテンシャル流れ

速度ポテンシャル

渦なし流れでは、速度場は勾配場として表せます。

$$\mathbf{v} = \nabla \phi$$

φを速度ポテンシャルと呼びます。

連続の式（$\nabla \cdot \mathbf{v} = 0$）を代入すると、

$$\nabla^2 \phi = 0$$

これはラプラス方程式です。調和関数の理論を流体力学に適用できます。

流れ関数

2次元非圧縮性流れでは、連続の式を自動的に満たす流れ関数 ψ を導入できます。

$$u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}$$

流れ関数の等値線（ψ = const.）は流線に一致します。

渦なし流れでは、ψ もラプラス方程式を満たします。

$$\nabla^2 \psi = 0$$

複素ポテンシャル

2次元ポテンシャル流れでは、複素関数論を活用できます。

複素速度ポテンシャル：

$$F(z) = \phi + i\psi$$

ここで $z = x + iy$ は複素座標です。F(z) が解析関数であれば、φとψは自動的にラプラス方程式を満たします（コーシー・リーマンの関係）。

複素速度：

$$W(z) = \frac{dF}{dz} = u - iv$$

基本的なポテンシャル流れ

一様流：

$$F(z) = Uz, \quad \phi = Ux, \quad \psi = Uy$$

湧き出し（source）：原点から流体が湧き出す。強さ Q（単位時間あたりの体積流量/単位奥行き）。

$$F(z) = \frac{Q}{2\pi}\ln z$$

吸い込み（sink）：Q が負の湧き出し。

渦糸（vortex）：原点まわりの循環 Γ を持つ流れ。

$$F(z) = -\frac{i\Gamma}{2\pi}\ln z$$

二重湧き出し（doublet）：湧き出しと吸い込みを近づけた極限。

$$F(z) = \frac{\kappa}{2\pi z}$$

円柱まわりの流れ

一様流と二重湧き出しを重ね合わせると、円柱まわりの流れが得られます。

$$F(z) = U\left(z + \frac{a^2}{z}\right)$$

ここで a は円柱の半径。

円柱表面上（$|z| = a$）で $\psi = 0$（流線）となることが確認できます。

循環を加えると、

$$F(z) = U\left(z + \frac{a^2}{z}\right) - \frac{i\Gamma}{2\pi}\ln z$$

この場合、円柱には揚力（クッタ・ジュコフスキーの定理）が発生します。

$$L = \rho U \Gamma$$

---

第7章：粘性流体の流れ

7.1 層流と乱流

流体の流れは、大きく層流（laminar flow） と乱流（turbulent flow） に分類されます。

層流：流体粒子が整然と層状に流れる。流線は交わらず、速度の時間変動がない（または小さい）。低レイノルズ数で発生。

乱流：流体粒子が不規則に混合しながら流れる。速度・圧力が激しく時間変動。渦構造が様々なスケールで存在。高レイノルズ数で発生。

レイノルズの実験（1883年）

レイノルズは円管内流れに着色した水を注入し、流れの様子を観察しました。

• 低流速：着色水は直線状のまま流れる（層流）
• 高流速：着色水は不規則に乱れ、管全体に広がる（乱流）
• 遷移は Re ≈ 2300 で起こる

7.2 管内層流（ハーゲン・ポアズイユ流れ）

円管内の定常な完全発達層流を考えます。

支配方程式

軸対称性から、速度は軸方向成分 $v_z(r)$ のみ。ナビエ・ストークス方程式のz成分：

$$0 = -\frac{dp}{dz} + \mu \frac{1}{r}\frac{d}{dr}\left(r\frac{dv_z}{dr}\right)$$

圧力勾配は軸方向に一定（$dp/dz = -\Delta p / L$）なので、

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{dv_z}{dr}\right) = -\frac{\Delta p}{\mu L}$$

速度分布

境界条件：

• $r = 0$ で速度が有限（軸上で特異点なし）
• $r = R$ で $v_z = 0$（滑りなし条件）

積分すると、

ハーゲン・ポアズイユ流れの速度分布

$$v_z(r) = \frac{\Delta p}{4\mu L}(R^2 - r^2) = v_{max}\left(1 - \frac{r^2}{R^2}\right)$$

• 放物線分布
• 最大速度：$v_{max} = \frac{\Delta p R^2}{4\mu L}$（中心軸上）
• 平均速度：$\bar{v} = \frac{v_{max}}{2}$

体積流量と圧力損失

体積流量：

$$Q = \int_0^R v_z(r) \cdot 2\pi r dr = \frac{\pi R^4 \Delta p}{8\mu L}$$

これをハーゲン・ポアズイユの法則と呼びます。流量は半径の4乗に比例するため、管径のわずかな変化が流量に大きく影響します。

圧力損失：

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

管摩擦係数

管摩擦係数（ダルシー・ワイスバッハの摩擦係数）f は、

$$\Delta p = f \frac{L}{D} \frac{\rho \bar{v}^2}{2}$$

層流の場合、

$$f = \frac{64}{Re}$$

7.3 境界層理論

境界層の概念

レイノルズ数が大きい流れでも、固体壁近傍では速度がゼロ（滑りなし条件）から主流速度まで変化する必要があります。この薄い領域を境界層（boundary layer） と呼びます。

プラントル（1904年）の境界層理論は、流体力学における最も重要な概念の一つです。

境界層の特徴：

• 厚さ δ は主流方向の距離 x に比べて非常に薄い（δ ≪ x）
• 境界層内では粘性効果が重要
• 境界層外では非粘性流れ（ポテンシャル流れ）として扱える

境界層方程式

2次元定常境界層方程式（プラントルの境界層方程式）：

$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = U\frac{dU}{dx} + \nu\frac{\partial^2 u}{\partial y^2}$$

$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$

ここで、U(x) は境界層外縁での主流速度です。

境界条件：

• $y = 0$：$u = v = 0$（壁面、滑りなし）
• $y \to \infty$：$u \to U(x)$（主流との接続）

平板上の層流境界層（ブラジウス解）

一様流中に置かれた平板上の境界層を考えます（$U = U_\infty = $ 一定）。

ブラジウスは相似解を求めました。相似変数：

$$\eta = y\sqrt{\frac{U_\infty}{\nu x}}$$

境界層厚さ（99%厚さ）：

$$\delta_{99} \approx 5.0 \sqrt{\frac{\nu x}{U_\infty}} = \frac{5.0 x}{\sqrt{Re_x}}$$

ここで $Re_x = U_\infty x / \nu$ は局所レイノルズ数。

壁面せん断応力：

$$\tau_w = 0.332 \rho U_\infty^2 \sqrt{\frac{\nu}{U_\infty x}}$$

局所摩擦係数：

$$C_f = \frac{\tau_w}{\frac{1}{2}\rho U_\infty^2} = \frac{0.664}{\sqrt{Re_x}}$$

排除厚さと運動量厚さ

境界層の影響を特徴づける積分量として、以下が重要です。

排除厚さ（displacement thickness）：

$$\delta^* = \int_0^\infty \left(1 - \frac{u}{U}\right) dy$$

主流が境界層の存在によって外側に押しやられる量を表します。

運動量厚さ（momentum thickness）：

$$\theta = \int_0^\infty \frac{u}{U}\left(1 - \frac{u}{U}\right) dy$$

壁面摩擦による運動量損失を表します。

境界層の剥離

逆圧力勾配（$dp/dx > 0$、すなわち $dU/dx < 0$）が存在すると、境界層内の流体は減速されます。十分に強い逆圧力勾配では、壁面近傍の流体が逆流し始め、境界層剥離（separation） が起こります。

剥離点では、壁面せん断応力がゼロになります。

$$\tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_{y=0} = 0$$

剥離は、物体の抵抗増大、揚力低下（失速）、振動・騒音の原因となります。

---

第8章：乱流

8.1 乱流の特徴

乱流は以下の特徴を持ちます。

1. 不規則性（irregularity）：速度・圧力が時間的・空間的にランダムに変動
2. 拡散性（diffusivity）：運動量・熱・物質の混合が層流より格段に大きい
3. 大きなレイノルズ数：慣性力が粘性力を大きく上回る
4. 3次元性：本質的に3次元的な渦構造
5. 散逸性（dissipation）：運動エネルギーが最終的に熱に変換される

8.2 レイノルズ分解

乱流の解析では、瞬時値を時間平均と変動成分に分解するレイノルズ分解が基本となります。

$$u = \bar{u} + u', \quad v = \bar{v} + v', \quad p = \bar{p} + p'$$

ここで、オーバーバーは時間平均（またはアンサンブル平均）、プライムは変動成分を表します。

定義から、変動成分の平均はゼロです。

$$\overline{u'} = 0$$

しかし、変動成分の積の平均は一般にゼロではありません。

$$\overline{u'v'} \neq 0$$

8.3 レイノルズ平均ナビエ・ストークス方程式（RANS）

ナビエ・ストークス方程式にレイノルズ分解を適用し、時間平均を取ると、

$$\bar{u}_j \frac{\partial \bar{u}_i}{\partial x_j} = -\frac{1}{\rho}\frac{\partial \bar{p}}{\partial x_i} + \nu \frac{\partial^2 \bar{u}_i}{\partial x_j \partial x_j} - \frac{\partial \overline{u'_i u'_j}}{\partial x_j}$$

最後の項 $-\rho\overline{u'_i u'_j}$ はレイノルズ応力（Reynolds stress） と呼ばれ、乱流による運動量輸送を表します。

レイノルズ応力テンソルは対称テンソルで、6つの独立成分を持ちます。しかし、方程式は4つ（連続の式1つ + 運動量3つ）しかないため、未知数の方が多い状態です。これが乱流の閉じ込め問題（closure problem） です。

8.4 乱流モデル

閉じ込め問題を解決するため、レイノルズ応力を平均流の量で表現する乱流モデルが開発されてきました。

渦粘性モデル

最も基本的なアプローチは、レイノルズ応力を平均ひずみ速度テンソルに比例させる渦粘性仮説（ブシネスク仮説）です。

$$-\overline{u'_i u'_j} = \nu_t \left( \frac{\partial \bar{u}_i}{\partial x_j} + \frac{\partial \bar{u}_j}{\partial x_i} \right) - \frac{2}{3}k\delta_{ij}$$

ここで、$\nu_t$ は渦粘性係数（渦動粘度）、$k = \frac{1}{2}\overline{u'_i u'_i}$ は乱流運動エネルギーです。

代表的な乱流モデル

混合距離モデル（Prandtl）：壁面近傍の流れに対する最も単純なモデル。

$$\nu_t = l_m^2 \left| \frac{\partial \bar{u}}{\partial y} \right|$$

$l_m = \kappa y$（κ ≈ 0.41 はカルマン定数）

k-εモデル：産業界で最も広く使われる2方程式モデル。乱流運動エネルギー k とその散逸率 ε の輸送方程式を解く。

$$\nu_t = C_\mu \frac{k^2}{\varepsilon}$$

k-ωモデル：境界層内の流れに適したモデル。ε の代わりに比散逸率 ω = ε/k を使用。

SSTモデル：k-ω と k-ε のハイブリッド。壁近傍では k-ω、主流では k-ε の特性を活かす。

レイノルズ応力モデル（RSM）：レイノルズ応力各成分の輸送方程式を解く。異方性を捉えられるが、計算コストが高い。

8.5 乱流境界層

乱流境界層は、壁面に近い順に以下の領域に分けられます。

粘性底層（viscous sublayer）：y⁺ < 5
壁のごく近傍で、乱流変動が壁面に抑えられ、粘性が支配。

$$u^+ = y^+$$

バッファ層：5 < y⁺ < 30
粘性と乱流の両方が重要な遷移領域。

対数則層（log-law region）：30 < y⁺ < 0.2δ⁺
慣性力が支配。普遍的な対数則が成り立つ。

$$u^+ = \frac{1}{\kappa} \ln y^+ + B$$

κ ≈ 0.41、B ≈ 5.0

外層（wake region）：y⁺ > 0.2δ⁺
境界層外縁に近い領域。普遍的な法則からずれる。

ここで、壁面単位（wall unit）として、

$$u^+ = \frac{u}{u_\tau}, \quad y^+ = \frac{y u_\tau}{\nu}$$

$u_\tau = \sqrt{\tau_w/\rho}$ は摩擦速度。

8.6 管内乱流

乱流円管流れの摩擦係数は、経験式で与えられます。

ブラジウスの式（滑らかな管、Re < 10⁵）：

$$f = \frac{0.316}{Re^{0.25}}$$

コールブルック・ホワイトの式（任意の粗さ）：

$$\frac{1}{\sqrt{f}} = -2\log_{10}\left( \frac{\varepsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}} \right)$$

ε は壁面粗さ。この陰的な式はムーディ線図としてまとめられています。

---

第9章：圧縮性流体力学

9.1 音速と圧縮性

音波は圧力・密度の微小な擾乱として伝播します。等エントロピー過程（断熱・可逆）を仮定すると、音速 c は、

$$c = \sqrt{\left(\frac{\partial p}{\partial \rho}\right)_s}$$

理想気体では、

$$c = \sqrt{\gamma R T}$$

ここで、γ = cp/cv は比熱比（空気では γ ≈ 1.4）。

標準大気（T = 288 K）での空気の音速は約 340 m/s です。

9.2 圧縮性流れの基礎方程式

圧縮性流れでは、密度変化を考慮する必要があります。

質量保存：

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$

運動量保存：

$$\rho \frac{D\mathbf{v}}{Dt} = -\nabla p + \nabla \cdot \boldsymbol{\tau} + \rho \mathbf{g}$$

エネルギー保存：

$$\rho \frac{D}{Dt}\left( e + \frac{v^2}{2} \right) = -\nabla \cdot (p\mathbf{v}) + \nabla \cdot (\boldsymbol{\tau} \cdot \mathbf{v}) + \nabla \cdot (k \nabla T)$$

状態方程式：

$$p = \rho R T$$

9.3 等エントロピー流れの関係式

断熱・可逆の流れ（等エントロピー流れ）では、重要な関係式が得られます。

$$\frac{T_0}{T} = 1 + \frac{\gamma - 1}{2}Ma^2$$

$$\frac{p_0}{p} = \left( 1 + \frac{\gamma - 1}{2}Ma^2 \right)^{\frac{\gamma}{\gamma-1}}$$

$$\frac{\rho_0}{\rho} = \left( 1 + \frac{\gamma - 1}{2}Ma^2 \right)^{\frac{1}{\gamma-1}}$$

ここで、添字0は淀み点（Ma = 0）の値を表します。

空気（γ = 1.4）の場合の数値例：

Ma	T/T₀	p/p₀	ρ/ρ₀
0	1.000	1.000	1.000
0.5	0.952	0.843	0.885
1.0	0.833	0.528	0.634
2.0	0.556	0.128	0.230
3.0	0.357	0.027	0.076

9.4 ノズル流れ

面積-速度関係

定常・1次元・等エントロピー流れにおける断面積 A と流速 v の関係は、

$$\frac{dA}{A} = (Ma^2 - 1)\frac{dv}{v}$$

この式から重要な結論が導かれます。

• 亜音速流れ（Ma < 1）：

  • 断面積が減少 → 速度が増加（先細ノズル）
  • 断面積が増加 → 速度が減少（ディフューザ）

• 超音速流れ（Ma > 1）：

  • 断面積が減少 → 速度が減少
  • 断面積が増加 → 速度が増加（末広ノズル）

• Ma = 1：dA = 0 でなければならない（のど部）

したがって、亜音速から超音速へ加速するには、先細末広ノズル（収束発散ノズル、ラバールノズル） が必要です。

チョーク（閉塞）

ノズルのど部で Ma = 1 になると、流量はそれ以上増加しません。これをチョークと呼びます。チョーク状態での質量流量は、

$$\dot{m}_{max} = A^* \frac{p_0}{\sqrt{T_0}} \sqrt{\frac{\gamma}{R}} \left( \frac{2}{\gamma+1} \right)^{\frac{\gamma+1}{2(\gamma-1)}}$$

$A^*$ はのど部面積。

9.5 衝撃波

超音速流れが障害物に遭遇したり、圧力が急激に上昇する状況では、衝撃波（shock wave） が発生します。衝撃波を横切ると、流体の状態が不連続的に変化します。

垂直衝撃波の関係式

垂直衝撃波の前後の状態（添字1が上流、2が下流）は、ランキン・ユゴニオの関係で結ばれます。

質量保存：$\rho_1 v_1 = \rho_2 v_2$

運動量保存：$p_1 + \rho_1 v_1^2 = p_2 + \rho_2 v_2^2$

エネルギー保存：$h_1 + \frac{v_1^2}{2} = h_2 + \frac{v_2^2}{2}$

これらから、マッハ数で整理すると、

$$Ma_2^2 = \frac{Ma_1^2 + \frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1}Ma_1^2 - 1}$$

$$\frac{p_2}{p_1} = \frac{2\gamma Ma_1^2 - (\gamma-1)}{\gamma+1}$$

$$\frac{\rho_2}{\rho_1} = \frac{(\gamma+1)Ma_1^2}{(\gamma-1)Ma_1^2 + 2}$$

重要な特徴：

• 衝撃波の前後では必ず Ma₁ > 1 → Ma₂ < 1（超音速→亜音速）
• 圧力、密度、温度は上昇
• 全圧は減少（不可逆過程、エントロピー増大）

---

第10章：応用と発展的話題

10.1 抗力と揚力

抗力の成分

物体に作用する流れ方向の力（抗力）は、以下の成分からなります。

摩擦抗力（friction drag）：壁面せん断応力による寄与

圧力抗力（pressure drag, form drag）：前後の圧力差による寄与

$$D = \underbrace{\int_S \tau_w \cos\theta \, dA}_{\text{摩擦抗力}} + \underbrace{\int_S p \sin\theta \, dA}_{\text{圧力抗力}}$$

抗力係数：

$$C_D = \frac{D}{\frac{1}{2}\rho U^2 A}$$

A は代表面積（翼なら翼平面積、鈍頭物体なら前面投影面積）。

揚力

流れに垂直な方向の力が揚力です。

$$L = \frac{1}{2}\rho U^2 A \cdot C_L$$

揚力の発生メカニズムは、翼上下面での圧力差です。これは循環Γと関連付けられます（クッタ・ジュコフスキーの定理）。

$$L = \rho U \Gamma$$

抗力係数のレイノルズ数依存性

球の抗力係数 C_D は、レイノルズ数に強く依存します。

• Re ≪ 1：ストークス域、$C_D = 24/Re$
• 1 < Re < 1000：遷移域
• 10³ < Re < 2×10⁵：ニュートン域、$C_D ≈ 0.44$
• Re > 2×10⁵：抗力危機、境界層が乱流に遷移し $C_D$ が急減少

ゴルフボールのディンプルは、境界層の乱流遷移を促進して剥離を遅らせ、圧力抗力を減少させる効果があります。

10.2 ターボ機械

ポンプ・ファン・圧縮機

流体に仕事を加える機械です。

遠心ポンプ：インペラ（羽根車）の回転により、流体に運動エネルギーと圧力を与える。

オイラーのターボ機械方程式：

$$W = \omega(r_2 v_{\theta 2} - r_1 v_{\theta 1})$$

W は単位質量あたりの仕事、ω は角速度、添字1,2は入口・出口。

比速度：ポンプの特性を表す無次元数。

$$N_s = \frac{N\sqrt{Q}}{H^{3/4}}$$

N：回転数、Q：流量、H：揚程

水車・タービン

流体から仕事を取り出す機械です。

ペルトン水車：衝撃式。高落差・小流量向け。

フランシス水車：反動式。中落差向け。最も広く使われる。

カプラン水車：反動式。低落差・大流量向け。

揚力型風車

風のエネルギーを回転エネルギーに変換。ベッツ限界により、風のエネルギーの最大59.3%しか取り出せない。

$$P_{max} = \frac{16}{27} \cdot \frac{1}{2}\rho A U^3$$

10.3 数値流体力学（CFD）

計算機の発展により、ナビエ・ストークス方程式を数値的に解く数値流体力学（CFD: Computational Fluid Dynamics） が実用的な設計ツールとなっています。

基本的な離散化手法

有限差分法（FDM）：微分を差分で近似。構造格子に適用しやすい。

$$\frac{\partial u}{\partial x} \approx \frac{u_{i+1} - u_{i-1}}{2\Delta x}$$

有限体積法（FVM）：保存則の積分形を離散化。保存性が良く、複雑形状に適用可能。最も広く使われる。

有限要素法（FEM）：変分原理に基づく。非構造格子に適し、構造解析との連成に向く。

代表的なCFD手法

直接数値シミュレーション（DNS）：乱流の全スケールを解像。最も正確だが、計算コストが膨大（Re³に比例）。研究用途に限定。

ラージエディシミュレーション（LES）：大きな渦は直接計算、小さな渦はモデル化。DNS と RANS の中間的手法。

レイノルズ平均（RANS）：時間平均方程式を解く。乱流モデルが必要。産業界で最も広く使用。

10.4 特殊な流れ

多相流

気液二相流、固液二相流など、複数の相が混在する流れ。

• ボイド率（気相の体積分率）
• 流動様式（気泡流、スラグ流、環状流など）
• 相間の運動量・熱・物質移動

非ニュートン流体

血液、高分子溶液、食品、塗料など。

• 擬塑性流体（shear-thinning）：せん断速度が増すと見かけ粘度が減少
• ダイラタント流体（shear-thickening）：せん断速度が増すと見かけ粘度が増加
• ビンガム流体：降伏応力を持つ

希薄気体力学

高高度大気、真空装置内、MEMS など、クヌーセン数が大きい流れ。

• 滑り流れ領域：壁面で滑り境界条件
• 遷移領域：ボルツマン方程式または DSMC 法
• 自由分子流：分子間衝突が無視できる

---

おわりに

流体力学は、古典力学の一分野でありながら、現代においても活発に研究が続けられている生きた学問です。乱流の完全な理解はまだ達成されておらず、ナビエ・ストークス方程式の解の存在と滑らかさは、ミレニアム懸賞問題の一つとして未解決のままです。

本記事では、流体力学の基礎概念から応用まで幅広く紹介しました。実際に流体力学を使いこなすには、本記事の内容をさらに深く学び、多くの例題・演習問題に取り組む必要があります。

参考文献・推薦図書

入門レベル

• 今井功『流体力学（前編）』（物理学選書）裳華房 - 定評のある古典的教科書
• 巽友正『流体力学』培風館 - 丁寧な解説で定評
• 日野幹雄『流体力学』朝倉書店 - 工学寄りの実用的な教科書

中級〜上級レベル

• Kundu, Cohen, Dowling『Fluid Mechanics』Academic Press - 世界的に使われている標準的教科書
• Batchelor『An Introduction to Fluid Dynamics』Cambridge - 理論的に厳密な古典
• White『Viscous Fluid Flow』McGraw-Hill - 粘性流体の詳細な解説
• Pope『Turbulent Flows』Cambridge - 乱流理論の決定版

数値流体力学

• Ferziger, Perić『Computational Methods for Fluid Dynamics』Springer
• Anderson『Computational Fluid Dynamics: The Basics with Applications』McGraw-Hill