• 学問 -top-

2.3.5.1 圧縮性流体（密度$\rho$ [kg/m$^3$ ]は変化する）

式(2.35)へ、式(2.36)、式(2.47)、式(2.71)、式(2.83)を入れると、以下の式が得られる。

$$\bigg( \rho \frac{\partial \bm{v}}{\partial t} + \bm{v} \frac{\partial \rho}{\partial t} \bigg) dxdydx = - \left( \begin{array}{c} \bm{\nabla} \cdot (\rho u \bm{v}) ... ...} \mu \bm{\nabla}(\bm{\nabla} \cdot \bm{v}) \right) dxdydz + \rho \bm{g} dxdydz$$

$$\rho \frac{\partial \bm{v}}{\partial t} + \bm{v} \frac{\partial \rho}{\partial t} = - \left( \begin{array}{c} \bm{\nabla} \cdot (\rho u \bm{v}) ... ...2 \bm{v} + \dfrac{1}{3} \mu \bm{\nabla}(\bm{\nabla} \cdot \bm{v}) + \rho \bm{g}$$

$$\rho \frac{\partial \bm{v}}{\partial t} + \bm{v} \frac{\partial \rho}{\partial t} = - \left( \begin{array}{c} \rho u \bm{\nabla} \cdot \bm{v} + \rh... ...2 \bm{v} + \dfrac{1}{3} \mu \bm{\nabla}(\bm{\nabla} \cdot \bm{v}) + \rho \bm{g}$$ (2.83)

質量保存式(2.33)の両辺に速度ベクトル$\bm{v}$ [m/s]を掛けると次式となる。

$$\bm{v} \frac{\partial \rho}{\partial t} = - \bm{v} \{\bm{\nabla} \cdot ( \rho \bm{v})\}$$

$$= - \left( \begin{array}{c} u v w \end{array} \right) \{\bm{\nabla} \cdot ( \rho \bm{v})\}$$

$$= - \left( \begin{array}{c} \rho u \bm{\nabla} \cdot \bm{v} + u \... ... \bm{\nabla} \cdot \bm{v} + w \bm{v} \cdot \bm{\nabla} \rho \end{array} \right)$$

上式を先ほど求めた運動量保存の式(2.84)から引き密度$\rho$ [kg/m$^3$ ]で割ると次式が得られる。ここで、$\nu$ [m$^2$ /s]は動粘性係数であり粘性係数$\mu$ [Pa$\cdot$ s]により $\nu = \mu / \rho$ と表される。

$$\frac{\partial \bm{v}}{\partial t} = - \left( \begin{array}{c} \bm{v} \cdot \bm{\nabla} u \bm{v} ... ...bla}^2 \bm{v} + \dfrac{1}{3} \nu \bm{\nabla}(\bm{\nabla} \cdot \bm{v}) + \bm{g}$$ (2.84)

$$\left( \begin{array}{c} \dfrac{\partial u}{\partial t} \vspace{.5... ...\partial t} \vspace{.5em} \dfrac{\partial w}{\partial t} \end{array} \right) = \left( \begin{array}{c} - \bm{v} \cdot \bm{\nabla} u - \dfrac{1... ...\partial z} \big( \bm{\nabla} \cdot \bm{v} \big) \bigg\} \end{array} \right)$$

• 学問 -top-

この図を含む文章の著作権は著者にあり、クリエイティブ・コモンズ 表示 - 非営利 - 改変禁止 3.0 非移植 ライセンスの下に公開する。