連続の方程式

カテゴリ：流体力学

水というのはたとえ流れの中にあっても量は変化しません。

100mlの水は、どこであろうと100mlだということです。

この質量保存則を記述したものが連続の方程式というものです。

ここでは、この連続の方程式を導出します。

連続の方程式の導出

まず、以下のように空間上のある点 $(x,y,z)$ にある体積 $\delta x\delta y\delta z$ の立方体を考えます。

ここでの流体の速度は $\mathbf{v}(x,y,z,t)=(u,v,w)$ です。

でははじめに、x軸に対して垂直な面から出入りする流体の量を計算します。

$x-\frac{\delta x}{2}$ から微小時間 $\delta t$ の間に流入する流体の質量は

　 $\displaystyle\rho\left(x-\frac{\delta x}{2},y,z,t\right)u\left(x-\frac{\delta x}{2},y,z,t\right)\delta y\delta z\delta t$

となります。これをテイラー展開すると

　 $\displaystyle\left(\rho(x,y,z,t)-\frac{\partial\rho}{\partial x}\frac{\delta x}{2}+O(\delta x^{2})\right)\left(u(x,y,z,t)-\frac{\partial u}{\partial x}\frac{\delta x}{2}+O(\delta x^{2})\right)\delta y\delta z\delta t\\=\left(\rho(x,y,z,t)u(x,y,z,t)-\frac{\partial\rho}{\partial x}\frac{\delta x}{2}u-\rho\frac{\partial u}{\partial x}\frac{\delta x}{2}+O(\delta x^{2})\right)\delta y\delta z\delta t\\=\left(\rho(x,y,z,t)u(x,y,z,t)-\left(\frac{\partial\rho}{\partial x}u+\rho\frac{\partial u}{\partial x}\right)\frac{\delta x}{2}+O(\delta x^{2})\right)\delta y\delta z\delta t\\=\left(\rho(x,y,z,t)u(x,y,z,t)-\frac{\partial}{\partial x}\left(\rho u\right)\frac{\delta x}{2}+O(\delta x^{2})\right)\delta y\delta z\delta t$

となります。二乗の項 $O(\delta x^{2})$ を無視すれば

　 $\displaystyle\left(\rho(x,y,z,t)u(x,y,z,t)-\frac{\partial}{\partial x}\left(\rho u\right)\frac{\delta x}{2}\right)\delta y\delta z\delta t$ ・・・①

が得られます。

次に、 $x+\frac{\delta x}{2}$ から微小時間 $\delta t$ の間に流出する流体の質量を計算します。すると

　 $\displaystyle\rho\left(x+\frac{\delta x}{2},y,z,t\right)u\left(x+\frac{\delta x}{2},y,z,t\right)\delta y\delta z\delta t$

となります。よって上記と同じように計算すれば

　 $\displaystyle\left(\rho(x,y,z,t)u(x,y,z,t)+\frac{\partial}{\partial x}\left(\rho u\right)\frac{\delta x}{2}\right)\delta y\delta z\delta t$ ・・・②

が得られます。

では、体積 $\delta x\delta y\delta z$ の増加量を考えます。まず、x軸に垂直な面から流入する質量は式①と②より

　 $\displaystyle\left(\rho(x,y,z,t)u(x,y,z,t)-\frac{\partial}{\partial x}\left(\rho u\right)\frac{\delta x}{2}\right)\delta y\delta z\delta t-\left(\rho(x,y,z,t)u(x,y,z,t)+\frac{\partial}{\partial x}\left(\rho u\right)\frac{\delta x}{2}\right)\delta y\delta z\delta t\\=-\frac{\partial}{\partial x}\left(\rho u\right)\delta x\delta y\delta z\delta t$

となります。あとは同様に $y,z$ 軸でやり総和をとれば

　 $\displaystyle-\left(\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)\right)\delta x\delta y\delta z\delta t$

が得られます。体積 $\delta x\delta y\delta z$ における質量の時間変化は

　 $\displaystyle\frac{\partial\rho}{\partial t}\delta x\delta y\delta z\delta t$

となるので、

　 $\displaystyle\frac{\partial\rho}{\partial t}\delta x\delta y\delta z\delta t=-\left(\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)\right)\delta x\delta y\delta z\delta t\\\Leftrightarrow \frac{\partial\rho}{\partial t}=-\left(\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)\right)\\\Leftrightarrow\frac{\partial\rho}{\partial t}+\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)=0$

となります。この方程式が連続の方程式です。密度が一定で変化しない流体の場合(水とかはほとんど密度一定なのでこの近似でいいでしょう)は

　 $\displaystyle\frac{\partial \rho}{\partial t}=0$

より、

　 $\displaystyle\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)=0$

が連続の方程式となります。

以上まとめると以下のようになります。

連続の方程式

密度 $\rho$ 、速度 $(u,v,w)$ における連続の方程式は

　 $\displaystyle\frac{\partial\rho}{\partial t}+\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)=0$

となる。非圧縮性流体の場合は

　 $\displaystyle\frac{\partial}{\partial x}\left(\rho u\right)+\frac{\partial}{\partial y}\left(\rho v\right)+\frac{\partial}{\partial z}\left(\rho w\right)=0$

となる。

方程式だけ見ると難しそうですが、ただ流体の質量は変わらないことを強調しているだけの単純な式です。

その単純さを理解するために、式の導出方法はきちんと理解しておくといいでしょう。

著者：安井真人(やすい　まさと)
@yasui_masatoさんをフォロー