非圧縮性流れ

非圧縮性流れ（ひあっしゅくせいながれ）とは流体力学において、流体粒子の内部で密度が一定の流体である。縮まない流体とも呼ばれる^[1]^[2]。連続体力学における非圧縮性の概念を流体に適用したものである。

言い換えると、非圧縮性とは流体の速度の発散が 0 になることである（この表現が等価である理由は後述）。

非圧縮性流れは、流体自体が非圧縮性であることを意味するものではない。圧縮性流体でも（適切な条件の下で）良い近似で非圧縮性流れとしてモデル化できる。非圧縮性流れは流体と同じ速度で移動する流体粒子の中で密度が一定であることを意味する。

非圧縮性流れに対して、密度が変化する流れを圧縮性流れという。厳密な意味での非圧縮性流れは自然界には存在しないが、一般的に流れのマッハ数（局所音速と流速との比）が小さい流れに対しては圧縮性の影響は無視できる。マッハ数が0.3を超えるか、または流体が非常に大きな圧力変化を受ける場合に、圧縮性の影響は考慮される。

速度の発散が 0 になることの導出

非圧縮性流れのための基本的な要件は、密度ρが、速度v で移動する流体の微小体積dV 内で一定であることである。数学的にはこの非圧縮性流れの条件は、密度の物質微分（後述）が 0 にならなければならないことを意味する。この条件を導入する前に、質量保存則から必要な関係式を導く。ある領域（コントロールボリューム）V の中にある流体の質量m は、密度ρの体積積分によって計算される。

m=\iiint _{V}\rho \,dV

質量保存則より、コントロールボリュームV 内の質量の時間変化はその境界面S を通る質量流束J に等しくなければならない。この質量流速は数学的には面積分で表される。

{\frac {\partial m}{\partial t}}=-\iint _{S}\!\!\!\!\!\!\!\!\!\!\!\subset \!\supset {\boldsymbol {J}}\cdot d{\boldsymbol {S}}

ここで負号は、面積ベクトルdS を外向きに定義していることより、コントロールボリュームから出る流れによってボリューム内の質量は時間的に減少することを意味する。この式の右辺に発散定理を用い、質量流束と密度の時間変化の間の関係が導かれる。

{\begin{aligned}&\iiint _{V}{\frac {\partial \rho }{\partial t}}\,dV=-\iiint _{V}\operatorname {div} {\boldsymbol {J}}\,dV\\&\therefore \quad {\frac {\partial \rho }{\partial t}}=-\operatorname {div} {\boldsymbol {J}}\end{aligned}}

左辺の密度の時間微分は非圧縮性流れを保証するためには 0 になる必要はない。この場合の密度の時間微分とは、固定位置のコントロールボリューム内でのこの変化率のことを言っている。密度の時間微分が 0 にならないことを許しても、非圧縮性流体は制限されない。なぜなら固定位置のコントロールボリュームを観察していれば、流れによってそこの密度は変化することが許されているからである。このアプローチは一般的に言え、密度の時間微分が 0 になってもよいことは、圧縮性流体が依然として非圧縮性流れになりうることを示している。

我々が興味を持っていることは、流体速度v とともに移動するコントロールボリュームの密度の変化である。流束J は、次の式で流体速度v と関連している：

{\boldsymbol {J}}=\rho {\boldsymbol {v}}

したがって質量保存則は

{\frac {\partial \rho }{\partial t}}+\operatorname {div} (\rho {\boldsymbol {v}})\equiv {\frac {\partial \rho }{\partial t}}+\operatorname {grad} \rho \cdot {\boldsymbol {v}}+\rho (\operatorname {div} {\boldsymbol {v}})=0

と表される。この式は連続の式と呼ばれている。密度の全微分を考える（そして連鎖律を適用する）ことにより、

{\frac {d\rho }{dt}}={\frac {\partial \rho }{\partial t}}+{\frac {\partial \rho }{\partial x}}{\frac {dx}{dt}}+{\frac {\partial \rho }{\partial y}}{\frac {dy}{dt}}+{\frac {\partial \rho }{\partial z}}{\frac {dz}{dt}}

であるから、流体と同じ速度で動いている（すなわちv = (dx /dt , dy /dt , dz /dt ) となる）コントロールボリュームをとれば、上式は物質微分：

{\frac {D\rho }{Dt}}:={\frac {\partial \rho }{\partial t}}+\operatorname {grad} \rho \cdot {\boldsymbol {v}}

を用いて簡潔に表される。連続の式に代入すれば、

{\frac {D\rho }{Dt}}=-\rho (\operatorname {div} {\boldsymbol {v}})

を得る。密度の時間微分は、流体が圧縮または膨張（一定体積dV に含まれる質量が変化）することを示しているが、それは禁止されている。密度の物質微分が 0 になることを要求することは、（密度が 0 でなければ）流体速度の発散が 0 でなければならないことと等価である。

\operatorname {div} {\boldsymbol {v}}=0

結局、質量保存則と、流体とともに移動する体積内の密度が一定であるという制約条件から、それに等価な非圧縮性流れの必要条件は流体速度の発散が 0 になることであることが示された。

圧縮率との関係

いくつかの分野では、流れの非圧縮性の尺度には、圧力変化による密度の変化がとられる。これは圧縮率βという用語で表現される。

\beta :={\frac {1}{\rho }}{\frac {d\rho }{dp}}

圧縮率が十分小さければ、流れは非圧縮性であるとされる。圧縮率が 0 であることを非圧縮性流体の定義とすることもある^[1]。

ソレノイド場との関係

非圧縮性流れはソレノイド（solenoidal）な速度場によって記述される。しかし、ソレノイド場は発散が 0 である以外に、0 でない回転を有する意味合いを持っている。

あるいは、非圧縮性流れの回転が 0 、つまり非回転なら、速度場は実際にはラプラス場（英語版）である。