• 学問 -top-

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

対流により各面で時間当たりに出入する量は式(2.193) $^{\text{\pageref{eq-component-adv}}}$ と各面の質量流量の式(2.17)-(2.22) $^{\text{p.\pageref{eq-mass_f_x-com}}}$ より次のように求まる。途中で相対する面の関係を式(2.11) $^{\text{p.\pageref{eq-diffplus}}}$ を用いて求める。

$x$ 軸に垂直 $yz$ 面左

$$\omega_{i,{x -}} \dot{m}_{x -}= \rho_{x -}\omega_{i,{x -}} u_{x -}dydz$$ (2.193)

$x$ 軸に垂直 $yz$ 面右

$$\omega_{i,{x +}} \dot{m}_{x +} = \rho_{x +}\omega_{i,{x +}} u_{x +}dydz$$

$$= \bigg( \rho_{x -} + \left. \frac{\partial \rho}{\partial x} \righ... ...-} + \left. \frac{\partial u}{\partial x} \right\vert _ {{x -}} d x \bigg) dydz$$

$$= \bigg( \rho_{x -}\omega_{i,{x -}} u_{x -}+ \rho_{x -}\omega_{i,{x... ...,{x -}} u_{x -}\left. \frac{\partial \rho}{\partial x} \right\vert _ {{x -}} dx$$

$$\underbrace{+ \rho_{x -}\left. \fra... ...視する（\ref{sec-DifferentialTerm}節 {\rm p}. \pageref{sec-DifferentialTerm}）}$$

$$\underbrace{ + \left. \frac{\partia... ...{sec-DifferentialTerm}節 {\rm p}. \pageref{sec-DifferentialTerm}）} \bigg) dydz$$

$$= \bigg( \rho_{x -}\omega_{i,{x -}} u_{x -}+ \rho_{x -}\omega_{i,{x... ... -}\left. \frac{\partial \rho}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$

$$= \bigg( \rho_{x -}\omega_{i,{x -}} u_{x -}+ \left. \frac{\partial (\rho \omega_{i} u)}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$ (2.194)

$y$ 軸に垂直 $zx$ 面下

$$\omega_{i,{y -}} \dot{m}_{y -}= \rho_{y -}\omega_{i,{y -}} v_{y -}dzdx$$ (2.195)

$y$ 軸に垂直 $zx$ 面上

$$\omega_{i,{y +}} \dot{m}_{y +} = \rho_{y +}\omega_{i,{y +}} v_{y +}dzdx$$

$$= \bigg( \rho_{y -}\omega_{i,{y -}} v_{y -}+ \left. \frac{\partial (\rho \omega_{i} v)}{\partial y} \right\vert _ {{y -}} dy \bigg) dzdx$$ (2.196)

$z$ 軸に垂直 $xy$ 面後

$$\omega_{i,{z -}} \dot{m}_{z -}= \rho_{z -}\omega_{i,{z -}} w_{z -}dxdy$$ (2.197)

$z$ 軸に垂直 $xy$ 面前

$$\omega_{i,{z +}} \dot{m}_{z +}= \rho_{z +}\omega_{i,{z +}} w_{z +}dxdy$$

$$= \bigg( \rho_{z -}\omega_{i,{z -}} w_{z -}+ \left. \frac{\partial (\rho \omega_{i} w)}{\partial z} \right\vert _ {{z -}} dz \bigg) dxdy$$ (2.198)

上六式から、まずそれぞれの軸に沿った出入を求める。$xyz$ の各軸にそってコントロールボリュームから出て行く方向が負、入る方向が正になるように符号を加え、向かい合う面を足し合わせる^2.26（2.1.6節 $^{\text{p.\pageref{sec-Difference}}}$ ）。

$x$ 軸に沿った出入

式(2.194)$-$ 式(2.195)

$$\rho_{x -} \omega_{i,{x -}} u_{x -}dydz - \bigg( \rho_{x -}\omega_{i,{x -}} ... ...{\partial (\rho \omega_{i} u)}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$

$$= - \left. \frac{\partial (\rho \omega_{i} u)}{\partial x} \right\vert _ {{x -}} dxdydz$$ (2.199)

$y$ 軸に沿った出入

式(2.196)$-$ 式(2.197)

$$\rho_{y -} \omega_{i,{y -}} v_{y -}dzdx - \bigg( \rho_{y -}\omega_{i,{y -}} ... ...{\partial (\rho \omega_{i} v)}{\partial y} \right\vert _ {{y -}} dy \bigg) dzdx$$

$$= - \left. \frac{\partial (\rho \omega_{i} v)}{\partial y} \right\vert _ {{y -}} dxdydz$$ (2.200)

$z$ 軸に沿った出入

式(2.198)$-$ 式(2.199)

$$\rho_{z -} \omega_{i,{z -}} w_{z -}dxdy - \bigg( \rho_{z -}\omega_{i,{z -}} ... ...{\partial (\rho \omega_{i} w)}{\partial z} \right\vert _ {{z -}} dz \bigg) dxdy$$

$$= - \left. \frac{\partial (\rho \omega_{i} w)}{\partial z} \right\vert _ {{z -}} dxdydz$$ (2.201)

xyz軸での出入の総和、式(2.200)$+$ 式(2.201)$+$ 式(2.202)をとる。ここで、全ての項がコントロールボリュームの体積($dxdydz$ )で括られているため各境界面での区別はせず(2.1.7節 $^{\text{p.\pageref{sec-DifferentialTerm}}}$ )、下付きを外す。コントロールボリューム全体での対流による成分の質量の出入は次式で表される^2.27。

$$- \dfrac{\partial (\rho \omega_i u)}{\partial x} dxdydz - \dfrac{... ... v)}{\partial y} dxdydz - \dfrac{\partial (\rho \omega_i w)}{\partial z} dxdydz$$

$$= - \bigg( \rho \omega_{i} \dfrac{\partial u}{\partial x} + \rho ... ...a_i}{\partial z} + \omega_{i} w \dfrac{\partial \rho}{\partial z} \bigg) dxdydz$$

$$= - \bigg\{ \rho \bigg( u \dfrac{\partial \omega_i}{\partial x} +... ...l \rho}{\partial y} + w \dfrac{\partial \rho}{\partial z} \bigg) \bigg\} dxdydz$$

$$= - \big( \rho \bm{v} \cdot \bm{\nabla} \omega_i + \rho \omega_i \bm{\nabla} \cdot \bm{v} + \omega_i \bm{v} \cdot \bm{\nabla}\rho \big) dxdydz$$

$$= - \bm{\nabla} \cdot (\rho \omega_i \bm{v}) dxdydz$$ (2.202)

• 学問 -top-

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