An example of convection. Though the flow may be steady (time-independent), the fluid decelerates as it moves down the diverging duct (assuming incompressible or subsonic compressible flow), hence there is an acceleration happening over position.
← ホームへ戻る

ナビエ・ストークス方程式流体力学の根幹をなす連続体力学の解析

🗓 2026年8月11日

私たちの身の回りに存在する空気や水などの流体の挙動を数学的に記述しようとする際、避けて通れないのがナビエ・ストークス方程式です。これは連続体力学(物質を離散的な分子の集まりではなく、連続した広がりを持つものとして扱う学問)に基づいた方程式であり、流体の運動量保存則を表現しています。本記事では、この複雑な方程式の基礎から、圧縮性と非圧縮性の違い、そして具体的な座標系での展開までを詳しく解説します。

Key Facts

  • 物理的本質:流体におけるニュートンの第2法則(運動量保存則)を記述したものである。
  • 主要変数:密度 ($\rho$)、流速 ($\mathbf{u}$)、圧力 ($p$)、粘性係数 ($\mu$) などで構成される。
  • 非圧縮性流体:密度が一定と見なされ、方程式が簡略化される(液体などに適用)。
  • 圧縮性流体:密度が変化することを考慮し、熱力学的な状態方程式が必要となる(高速の気体などに適用)。
  • 非線形性:対流項(流速が自身の速度に影響を与える項)により、乱流などの複雑な現象が発生する。

流体運動の基礎概念と保存則

ナビエ・ストークス方程式を理解するためには、まず連続体力学の視点が必要です。流体は、質量、運動量、エネルギーの3つの保存則に従います。特に運動量保存則は、流体に働く力(圧力勾配や粘性力、外力)が流体の加速を引き起こすことを示しています。

ここで重要な概念が対流加速(Convective acceleration)です。これは、流体粒子が空間的に移動することで、速度が変化する現象を指します。例えば、管の形状が変化して流路が広がると、流速が低下し、位置による加速(減速)が生じます。

An example of convection. Though the flow may be steady (time-independent), the fluid decelerates as it moves down the diverging duct (assuming incompressible or subsonic compressible flow), hence there is an acceleration happening over position.

Cauchy応力テンソルと構成方程式

流体内部で発生する応力は、Cauchy応力テンソルを用いて記述されます。ニュートン流体の場合、応力は速度勾配(ひずみ速度)に比例するという線形な関係を持ちます。この関係式を運動量保存則に組み込むことで、ナビエ・ストークス方程式が導出されます。

圧縮性と非圧縮性のモデル化

解析対象となる流体の性質によって、方程式の形態は大きく2つに分かれます。

圧縮性流体(Compressible Flow)

気体のように密度が圧力や温度で変化する場合、圧縮性流体として扱います。この場合、方程式には密度の時間変化や空間変化が含まれ、より複雑な項(体積粘性係数 $\xi$ など)が登場します。熱力学的な影響を強く受けるため、エネルギー方程式との連立解析が不可欠です。

非圧縮性流体(Incompressible Flow)

水などの液体や、低速の気流では密度を一定と見なす「非圧縮性」の仮定が用いられます。これにより、連続の式は単純な速度場の発散ゼロ ($\nabla \cdot \mathbf{u} = 0$) となり、計算負荷が大幅に軽減されます。このとき、圧力場はポアソン方程式に従うことになります。

多様な流動パターンと解析手法

ナビエ・ストークス方程式は、境界条件や流れの対称性に応じて様々な解を持ちます。代表的な例として、流線が平行に流れる平行流や、一点から放射状に広がる放射流が挙げられます。

Visualization of (a) parallel flow and (b) radial flow

また、より高度な数学的アプローチとして、ホップファイブレーション(Hopf fibration)に沿った流線を持つ3次元の定常渦解などの特殊解も研究されています。これらは流体のトポロジー的な性質を理解する上で重要な役割を果たします。

Wire model of flow lines along a Hopf fibration

座標系による表現の最適化

解析する問題の形状に合わせて、方程式は異なる座標系で記述されます。直交座標系(デカルト座標系)が基本ですが、円筒形や球形の物体を扱う場合は、それぞれ円筒座標系球座標系に変換することで、対称性を利用して変数を削減し、計算を効率化できます。

流体力学の主要概念まとめ

流体力学における主要な法則と概念の整理
カテゴリー 主要な法則・概念 物理的意味
基本保存則 質量・運動量・エネルギー保存 物理量の総量は不変であり、変化は流入出と生成消滅で決まる
流体特性 粘性 (Viscosity) 流体の「ねばりけ」。内部摩擦として運動を抑制する力
代表的な法則 ベルヌーイの定理 / パスカルの原理 エネルギー保存や圧力伝播に関する簡略化された法則
拡散現象 フィックの法則 (Fick's laws) 濃度勾配による物質の拡散速度を記述
特殊流体 スマート流体 (電磁流体など) 外部電場や磁場によって粘性が変化する流体

Frequently Asked Questions

ナビエ・ストークス方程式が「難しい」と言われる理由は何ですか?

最大の原因は、方程式に含まれる非線形項(対流項)にあります。これにより、わずかな条件の変化が劇的な結果の違いを生む「乱流」という現象が発生し、数学的な一般解を得ることが極めて困難であるためです。

ニュートン流体と非ニュートン流体の違いは何ですか?

ニュートン流体は、せん断応力がひずみ速度に正比例する流体(水や空気など)です。一方、非ニュートン流体は、加える力によって粘性が変化する流体(血液やケチャップなど)を指します。

非圧縮性の仮定はどのような時に利用できますか?

一般的に、流体の流速が音速よりも十分に遅い(マッハ数が小さい)場合や、液体のように外部圧力による体積変化が無視できる場合に適用されます。

連続体力学とは具体的にどのような考え方ですか?

物質を個々の分子や原子の集まりとしてではなく、空間的に連続して分布している「連続体」として扱う手法です。これにより、微分方程式を用いて密度や速度などの物理量を滑らかに記述することが可能になります。

粘性係数 $\mu$ は流体の挙動にどう影響しますか?

粘性係数が大きいほど、流体内部の摩擦が強くなり、流れは安定して層流(ラミナーフロー)になりやすくなります。逆に粘性が低いと、慣性力が支配的になり、複雑な渦を巻く乱流へと移行しやすくなります。

References

  1. Refer to the mathematical operator represented by the nabla ( {\displaystyle \nabla } ) symbol.

📸 フォトギャラリー

An example of convection. Though the flow may be steady (time-independent), the fluid decelerates as it moves down the diverging duct (assuming incompressible or subsonic compressible flow), hence there is an acceleration happening over position.
Visualization of (a) parallel flow and (b) radial flow
Wire model of flow lines along a Hopf fibration