よみもの:音速

	よみものページ
	この「よみもの」は、Imuzakが作成したものです。 他ユーザーの編集も歓迎いたします。

音を物理的に定義し、ついでに音速を計算してみましょう。

使うもの

参考になりそうなwebページにリンクしてあります。

• 高校レベルの数学の知識(公式をたくさん暗記する必要はない)
• 偏微分の知識(高校数学でいう「微分」は常微分と呼ばれる)　偏微分
• テイラー展開の知識　テイラー展開
• 古典力学の知識　力学
• 熱力学の基礎知識　熱力学

以上が用意できた方はこのページに沿ってやってみましょう!

「何も知らない」、「何から勉強したらいいか全くわからない」という方は、「関数って一体何者?」という問いかけから微分を勉強してみるといいと思います。

音の定義と音速の導出

図1のように、断面積が$S$の円柱状の領域について、この領域がある物質(弾性体)で満たされていて、この領域の側面で物質の出入りがないような状況を考える。
この領域内の物質に密度の偏りがないとき、図1の斜線部のような部分にある物質の質量を$dm$とおくと、この部分は体積が$Sdx$であるから、このときのこの部分の質量密度は
${\rho }_0=\frac{dm}{Sdx}$…(1)

である。さて、図1の状況から$x$軸方向に物質がわずかに移動して物質に密度の偏りが出来たとき、図1の斜線部分が図2の斜線部分のように変化したとき、変化後の斜線部分の幅$d{x}^{\prime }$を求める。斜線部分の両端の変位を図1の状況での位置$x$と時刻$t$の関数として$u=u(x,t)$と表わすと、
$d{x}^{\prime }=dx+u(x+dx,t)-u(x,t)$

となる。ここから

$d{x}^{\prime }=dx(1+\frac{u(x+dx,t)-u(x,t)}{dx})$

と式変形できる。ここで、$dx$が微小量だとすると、偏微分の定義から

$d{x}^{\prime }=dx(1+\frac{\partial u}{\partial x})$…(2)

となる。また、このとき斜線部分の質量密度は$\rho =\frac{dm}{Sd{x}^{\prime }}$…(3)　である。(3)に(2)を代入すると、

$\rho =\frac{dm}{Sdx(1+\frac{\partial u}{\partial x})}$…(4) となる。また、$\frac{\partial u}{\partial x}$が位置$x$と時刻$t$の関数なので、(4)が成り立つなら質量密度$\rho$も

$\rho =\rho (x,t)$…(5)　と表わせる

(4)に(1)を代入して整理すると、

$\rho ={\rho }_0{(1+\frac{\partial u}{\partial x})}^{-1}$…(6)　となる。

ここで、${(1+\frac{\partial u}{\partial x})}^{-1}$を$\frac{\partial u}{\partial x}=0$のまわりでテイラー展開すると、

${(1+\frac{\partial u}{\partial x})}^{-1}=1-\frac{\partial u}{\partial x}+{\left(\frac{\partial u}{\partial x}\right)}^2-2{\left(\frac{\partial u}{\partial x}\right)}^3+\cdots$

となる。ここで、$\left|\frac{\partial u}{\partial x}\right|$が1より充分小さいと仮定すると、

${(1+\frac{\partial u}{\partial x})}^{-1}=1-\frac{\partial u}{\partial x}$…(7)　と近似できる。(6)に(7)を代入すると、

$\rho ={\rho }_0(1-\frac{\partial u}{\partial x})$　となる。両辺を$x$で偏微分すると、

$\frac{\partial \rho }{\partial x}=-{\rho }_0\frac{{\partial }^{2}u}{\partial x^2}$…(8)　と表わせることがわかる。

ここで、この状況を力学的な視点で考察しよう。今扱っている物質の性質として、圧力$P$が質量密度$\rho$の関数として、$P=P(\rho )$と表わせると仮定すると、ここに(5)を代入すると、

$P=P(\rho (x,t))$と表わせることもわかる。

さて、質量が$dm$であるような微小部分についての運動方程式は、この微小部分の加速度が$\frac{{\partial }^{2}u}{\partial t^2}$と表わせることと、図3のような力が掛かっていることから、

$dm\frac{{\partial }^{2}u}{\partial t^2}=P(\rho (x,t))S-P(\rho (x+dx,t))S$…(9)

とわかる。ここに(1)を代入すると、

${\rho }_{0}Sdx\frac{{\partial }^{2}u}{\partial t^2}=P(\rho (x,t))S-P(\rho (x+dx,t))S$

となる。よって

${\rho }_0\frac{{\partial }^{2}u}{\partial t^2}=-\frac{P(\rho (x+dx,t))-P(\rho (x,t))}{dx}$

となる。さらに、

${\rho }_0\frac{{\partial }^{2}u}{\partial t^2}=-\frac{P(\rho (x+dx,t))-P(\rho (x,t))}{\rho (x+dx,t)-\rho (x,t)}\frac{\rho (x+dx,t)-\rho (x,t)}{dx}$

と式変形できる。ここで微分の定義より、

${\rho }_0\frac{{\partial }^{2}u}{\partial t^2}=-\frac{dP}{d\rho }\frac{\partial \rho }{\partial x}$…(10)

と表わせることがわかる。(10)に(8)を代入して整理すると、

$\frac{{\partial }^{2}u}{\partial t^2}=\frac{dP}{d\rho }\frac{{\partial }^{2}u}{\partial x^2}$…(11)

となる。(11)のような形の方程式は一次元の波動方程式と呼ばれていて、この一般解は

$u=f_1(x-\sqrt{\frac{dP}{d\rho }}t)+f_2(x+\sqrt{\frac{dP}{d\rho }}t)$ 　(ただし、$f_1,f_2$の関数形は任意)

である。これは、$u$は関数形が変化せずに$x$軸正の向きと負の向きに速さ$\sqrt{\frac{dP}{d\rho }}$で伝わる関数の和で表せることを意味している。

これは波の定義を満たしている。この波を音と定義する。

空気中の音速

ここまで考えてきた「物質」が理想気体として扱えるような気体だとしたら、音の速さはどのように計算できるだろうか。熱力学を用いて考えてみよう。
その気体中の音の速さが気体中の熱伝導の速さよりも充分大きいと仮定すると、音によって気体の微小部分は断熱変化すると考えることができる。
その気体の定圧比熱を$C_V$、定積比熱を$C_P$として、比熱比$\gamma =\frac{C_V}{C_P}$を用いて微小部分における圧力$P$と微小部分の体積$V$の関係式は

$P{V}^{\gamma }=(const)$…(12)

と表わせる。これは熱力学から導かれたポアソンの公式である。また、この気体の平均分子量を$M$とし、この微小部分に含まれる気体の物質量を$n$とおくと、この気体の質量密度$\rho$は

$\rho =\frac{Mn}{V}$…(13)　と表わせる。よって、

$V=\frac{Mn}{\rho }$…(14)　といえる。(14)を(12)に代入して整理すると、

$P=\frac{(const)}{(Mn{)}^{\gamma }}{\rho }^{\gamma }$…(15)　となる。両辺を$\rho$で微分すると、

$\frac{dP}{d\rho }=\frac{(const)}{(Mn{)}^{\gamma }}\gamma {\rho }^{\gamma -1}$　となる。式変形すると、

$\frac{dP}{d\rho }=\frac{\gamma }{\rho }\frac{(const)}{(Mn{)}^{\gamma }}{\rho }^{\gamma }$…(16)　となる。(16)に(15)を代入すると、

$\frac{dP}{d\rho }=\frac{\gamma }{\rho }P$…(17) と表わせる。また、気体定数を$R$とし、この気体の温度を$T$とすると、この気体の状態方程式は

$PV=nRT$　である。式変形すると

$P=\frac{nRT}{V}$…(18)　となる。(17)に(13)と(18)を代入すると、

$\frac{dP}{d\rho }=\frac{\gamma }{\left(\frac{Mn}{V}\right)}\left(\frac{nRT}{V}\right)$ となる。整理すると、

$\frac{dP}{d\rho }=\frac{\gamma RT}{M}$ となる。

前節で$\sqrt{\frac{dP}{d\rho }}$が音の速さであるとわかっているので、

$\sqrt{\frac{dP}{d\rho }}=\sqrt{\frac{\gamma RT}{M}}$ が空気中の音の速さである。

この式に、比熱比$\gamma =1.403$(一般的な組成の空気の比熱比)、気体定数$R=8.314\mathrm{J}/\mathrm{(}\mathrm{K}\mathrm{\cdot }\mathrm{m}\mathrm{o}\mathrm{l}\mathrm{)}$、温度$T=293\mathrm{K}$　(20℃)　、分子量$M=28.8\mathrm{g}/\mathrm{m}\mathrm{o}\mathrm{l}$を代入すると、音の速さは344m/sとわかる。これは同じ条件での音速の実測値と非常に近い値である。
日常的に使いやすい音速の近似式としては、気温を$t$とおくと、常温($T_0=293K$、つまり20℃)に近い値において、

$\sqrt{\frac{\gamma R(T_0+\mathrm{\Delta }T)}{M}}=\sqrt{\frac{\gamma R{T}_0}{M}}+\frac{1}{2}\sqrt{\frac{\gamma R}{M{T}_0}}(t-20^{\circ }\text{C})$


$=(344+0.636\times (t-20))\mathrm{m}/\mathrm{s}$と表すことができる。