流体の密度伝搬は音速で伝わる。圧縮性の考慮は音速と流速の大小関係が重要

どうも(^^)

前回の続きになりますが、前回の記事では流体解析には流れの性質をよく理解しておく必要があることを述べました。

特に、流体解析において「非圧縮性流体」「圧縮性流体」に区別は特に重要であり、両者の区別はマッハ数がとても重要なパラメータになることを示しました。

マッハ数\(M=\frac{u}{a}\)
※\(a\)：音速、\(u\)：代表の流速

マッハ数はこんな感じで書けるのですが、分母に音速\(a\)というのがあります。

この時点では、音速\(a\)が何者かわかりませんよね（笑）

でもこの音速\(a\)が「非圧縮性流体」「圧縮性流体」の判別にとても重要な物理量なのですが、なぜ重要なのかを本記事で示したいと思います。

本記事の目標

「非圧縮性流体」「圧縮性流体」の判別に音速が重要であることを示す。

理解したい内容をもう一度整理しておくと、

• 非圧縮性流体・圧縮性の判別をしたい

↓

• 密度変化が無視できるかどうかを見積もれば良い

↓

• マッハ数が見積もりに使われるようだ

↓

• 流速と音速の大小が重要である：本記事の内容

このフローを逆に辿って理解していけば、非圧縮性流体・圧縮性流体の判別を理解できます。（少しづつ投稿していきます。）

密度波

非圧縮性流体と考えれるかどうかは、密度の変化を考慮するかどうかを考えれば良いわけですよね。

ということは、流体における密度伝搬がどういった速度で伝わっているかを考えれば良いということになります。

イメージ

y
↑→x

イメージはこんな感じでしょうか。
密度が疎密波で伝搬する（密度波）といった感じです。
そこを流体が流速\(u\)で流れているというイメージです。

流体力学の基礎方程式は、基本的には流速\(u\)を求めるためのものです。
ですので、この密度波が考慮した方が良いのかどうかというのは、「密度伝搬の速度」と「流速」の大小を比較することに他なりません。

この密度伝搬の速度のことを音速と言います。

密度伝搬の速度と流速が同程度なら密度変化を考慮しなくはいけないということになりますね(^^)
この密度伝搬の速度を見積れば良いということになりますので、早速やってみましょう(^^)/

密度伝搬の速度：微小じょう乱

非圧縮性流体の基礎方程式：1次元

簡単のために一次元で考えましょう(^^)

粘性も熱拡散も無視します。(\(\nu=0,\alpha=0\))

※粘性は流速の垂直方向の速度勾配(\(\nu=\mu/\rho=\frac{1}{\rho}\frac{\partial u(x)}{\partial y}\))があるときに重要なので一様な流れを考えている場合は考慮しなくても良いでしょう。
温度に関しても温度を一様としているので熱拡散はないものとしています。
しかし、微小じょう乱が入り一様で無くなる場合は本来考慮しなくてはいけません。
粘性\(\nu\frac{\partial^2 u}{\partial x^2}\)も熱拡散も\(\alpha\frac{\partial^2 T}{\partial x^2}\)もいくら係数が小さくても微小じょう乱によって変動する場合は二階微分が大きくなるのでさすがに無視できなくなります（実現象と異なってきます）。

いきなり複雑に扱ってもよくわからなくなるかという思いで、ここはいったん粘性も熱拡散も無視しましょう。

そうしたときの、流体の基礎方程式を列挙します。

質量保存則

\begin{align*}\frac{\partial \rho}{\partial t}+\frac{\partial}{\partial x}\big(\rho u\big)=0\cdot\cdot\cdot (1)\end{align*}

 

ナビエストークス方程式

\begin{align*}\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}=-\frac{1}{\rho}\frac{\partial P}{\partial x}+\nu\frac{\partial^2 u}{\partial x^2}\cdot\cdot\cdot (2)\end{align*}

↓簡単に考えるため、粘性を無視します\(\nu=0\)

\begin{align*}\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}=-\frac{1}{\rho}\frac{\partial P}{\partial x}\cdot\cdot\cdot (2)\end{align*}

 

エネルギー保存則

\begin{align*}\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=\alpha\frac{\partial^2 T}{\partial x^2}\cdot\cdot\cdot (3)\end{align*}

↓簡単に考えるため、拡散項を無視します\(\alpha\)

\begin{align*}\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0\cdot\cdot\cdot (3)\end{align*}

 

状態方程式

系の状態と言うのは、圧力\(P\)、密度\(\rho\)、温度\(T\)で決まるのですが、熱平衡状態にある場合は、

\begin{align*}f(P,\rho,T)=0\cdot\cdot\cdot (4)\end{align*}

のような関係式（これを状態方程式と言う）があり、独立な変数は2つだけです。（未知数３つと(4)式なので自由度は２つ）

さらに、(4)式から温度の輸送は無視できるので温度は(4)式から除くことができ、

\begin{align*}f(P,\rho)=0\cdot\cdot\cdot (5)\end{align*}

と書けます。

つまり、\(P\)が\(\rho\)に依存しているという形ですね。

ですので、(5)式を

\begin{align*}P=P(\rho)\cdot\cdot\cdot (6)\end{align*}

と置けます。

微小じょう乱を導入

下記のように一様に流れている流体（速度\(u_{0}\)、密度\(\rho_{0}\)、圧力\(P_{0}\)）に、微小じょう乱\({u}’、{\rho}’、{P}’\)を考えます。

そうすると、速度、密度、圧力はそれぞれ

• 速度：\(u=u_{0}+{u}’\)・・・(7)
• 密度：\(\rho=\rho_{0}+{\rho}’\)・・・(8)
• 圧力：\(P=P_{0}+{P}’\)・・・(9)

となります。

ここで考えたいのは、微小じょう乱\(\rho\)がどのような速度で伝播するかです。

さて、(7)(8)(9)式を(1)(2)に代入してみましょう。

結果

\begin{align*}\frac{\partial {\rho}’}{\partial t}+\rho_{0}\frac{\partial {u}’}{\partial x}=0\cdot\cdot\cdot (10)\end{align*}

\begin{align*}\frac{\partial {u}’}{\partial t}+\frac{1}{\rho_{0}}\big(\frac{dP}{d\rho}\big)_{0}\frac{\partial {\rho}’}{\partial x}=0\cdot\cdot\cdot (11)\end{align*}

いきなり結果を示しましたが、気になる方は下記を参考にやってみてください。

(1)式に(7)(8)式を代入

これは簡単です(^^)

\begin{align*}\frac{\partial (\rho_{0}+{\rho}’)}{\partial t}+\frac{\partial (\rho_{0}+{\rho}’)(u_{0}+{u}’)}{\partial x}=0\end{align*}

\begin{align*}\frac{\partial (\rho_{0}+{\rho}’)}{\partial t}+\frac{\partial (\rho_{0}u_{0}+\rho_{0}{u}+{\rho}’u_{0}+{\rho}'{u}’)(u_{0}+{u}’)}{\partial x}=0\end{align*}

微小じょう乱の2次の項は無視すると、

\begin{align*}\frac{\partial {\rho}’}{\partial t}+\rho_{0}\frac{\partial {u}’}{\partial x}=0\cdot\cdot\cdot (10)\end{align*}

(2)式に(7)(8)(9)式を代入

\begin{align*}\frac{\partial (u_{0}+{u}’)}{\partial t}+ (u_{0}+{u}’)\frac{\partial (u_{0}+{u}’)}{\partial x}=-\frac{1}{\rho_{0}+{\rho}’}\frac{\partial (P+{P}’)}{\partial x}\end{align*}

右辺の分母については、\(\rho_{0}>>{\rho}’\)なので\(\rho_{0}+{\rho}=\rho_{0}\)。

さらに微小じょう乱の2次の項は無視すると、

\begin{align*}\frac{\partial {u}’}{\partial t}=-\frac{1}{\rho_{0}}\frac{\partial (P+{P}’)}{\partial x}\end{align*}

ここで、右辺の\(P+{P}’\)は(6)式より「密度変動したことによって圧力変動が起こった」と考えるはずなので、

\(P+{P}’=P(\rho_{0}+{\rho}’)\)として、これを\(\rho=\rho_{0}\)のまわりでテーラー展開しておきましょう。

\begin{align*}P(\rho_{0}+{\rho}’)=P_{0}(\rho_{0})+\big(\frac{\partial P}{\partial \rho}\big)_{0}{\rho}’\end{align*}

そして、xで偏微分を行って、

\begin{align*}\frac{\partial P(\rho_{0}+{\rho}’)}{\partial x}=\frac{\partial P_{0}(\rho_{0})}{\partial x}+\big(\frac{\partial P}{\partial \rho}\big)_{0}\frac{\partial {\rho}’}{\partial x}\end{align*}

よって、

\begin{align*}\frac{\partial {u}’}{\partial t}+\frac{1}{\rho_{0}}\big(\frac{dP}{d\rho}\big)_{0}\frac{\partial {\rho}’}{\partial x}=0\cdot\cdot\cdot (11)\end{align*}

さてもう一息です。

\(\frac{\partial}{\partial t}\)×(10)ー\(\rho_{0}\)\(\frac{\partial}{\partial x}\)×(11)として、\({u}’\)を消しましょう。

そうすると、

\begin{align*}\frac{\partial^2 {\rho}’}{\partial t^2}-\big(\frac{\partial P}{\partial \rho}\big)_{0}\frac{\partial^2 {\rho}’}{\partial x^2}=0\end{align*}

となりました!(^^)!

ここで、\(a_{0}=\sqrt{\big(\frac{\partial P}{\partial \rho}\big)_{0}}\)と置いたらどうでしょう。

\begin{align*}\frac{\partial^2 {\rho}’}{\partial t^2}-a_{0}^2\frac{\partial^2 {\rho}’}{\partial x^2}=0\cdot\cdot\cdot (12)\end{align*}

これは波動方程式ですね!(^^)!

つまり、微小な密度変動は\(a_{0}=\sqrt{\big(\frac{\partial P}{\partial \rho}\big)_{0}}\)で伝播することになります。

この密度伝播を音波と言います。

\begin{align*}a_{0}=\sqrt{\big(\frac{\partial P}{\partial \rho}\big)_{0}}\end{align*}を音速と言います。

ですので、密度変動は音速で伝播するわけですので、密度変動を考慮するかどうかというのはこの音速の値が相対的に無視できるかどうかということにつながりますね。

音速が流速に比べてものすごく大きい場合は、密度変動が一瞬で端まで伝播するかですから、それは密度変化がないというのと等価であると考えてよいですし、逆に音速が流速と同等である場合は、密度変動は決して無視できないということになります。

理想気体の状態方程式から空気や水の音速

結果だけ示しておきましょう(^^)/

水の音速が大きいですね(笑)

物理量の伝播は疎密波で伝わります。つまり、縦波で伝播していくので伝播するために媒介してくれるものがある方がより速く伝播することになりますよね。

ということを考えると、水は空気に比べて密度が大きいので音速が速いという理解になります。

ですので、多くの場合「水は非圧縮性流体（密度変動が無視できる流体）」として扱っても良いことになります。

でも音速と同等なくらいの流速であれば、水であっても圧縮性を考慮しなくてはいけないです（「水＝非圧縮性流体」が必ず成り立つわけではないことに注意です）。

さらに、音速も温度依存しますので、いつも上記の速度で伝播するということではないということに注意しておきましょう(^^)/

まとめ

微小な密度変動は、\begin{align*}a_{0}=\sqrt{\big(\frac{\partial P}{\partial \rho}\big)_{0}}\end{align*}で伝播する。

音速

\begin{align*}a_{0}=\sqrt{\big(\frac{\partial P}{\partial \rho}\big)_{0}}\end{align*}