• 学問 -top-

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

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

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

$$\dot{m}_{x -}c_v T_{x -}= \rho_{x -}c_v T_{x -}u_{x -}dydz$$ (2.139)

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

$$\dot{m}_{x +} c_v T_{x +}= \rho_{x +}c_v T_{x +}u_{x +}dydz$$

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

$$= c_v \bigg( \rho_{x -}T_{x -}u_{x -}... ...視する（\ref{sec-DifferentialTerm}節 {\rm p}. \pageref{sec-DifferentialTerm}）}$$

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

$$= c_v \bigg( \rho_{x -}T_{x -}u_{x -}+ \rho_{x -}T_{x -}\left. \fra... ... -}\left. \frac{\partial \rho}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$

$$= c_v \bigg( \rho_{x -}T_{x -}u_{x -}+ \left. \frac{\partial (\rho T u)}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$ (2.140)

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

$$\dot{m}_{y -}c_v T_{y -}= \rho_{y -}c_v T_{y -}v_{y -}dzdx$$ (2.141)

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

$$\dot{m}_{y +} c_v T_{y +}= \rho_{y +}c_v T_{y +}v_{y +}dzdx$$

$$= c_v \bigg( \rho_{y -}T_{y -}v_{y -}+ \left. \frac{\partial (\rho T v)}{\partial y} \right\vert _ {{y -}} dy \bigg) dzdx$$ (2.142)

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

$$\dot{m}_{z -}c_v T_{z -}= \rho_{z -}c_v T_{z -}w_{z -}dxdy$$ (2.143)

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

$$\dot{m}_{z +} c_v T_{z +}= \rho_{z +}c_v T_{z +}w_{z +}dxdy$$

$$= c_v \bigg( \rho_{z -}T_{z -}w_{z -}+ \left. \frac{\partial (\rho T w)}{\partial z} \right\vert _ {{z -}} dz \bigg) dxdy$$ (2.144)

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

$x$ 軸に沿った出入

式(2.140)$-$ 式(2.141)

$$\rho_{x -} c_v T_{x -}u_{x -}dydz - c_v \bigg( \rho_{x -}T_{x -}u_{x -}+ \left. \frac{\partial (\rho T u)}{\partial x} \right\vert _ {{x -}} dx \bigg) dydz$$

$$= - c_v \left. \frac{\partial (\rho T u)}{\partial x} \right\vert _ {{x -}} dxdydz$$ (2.145)

$y$ 軸に沿った出入

式(2.142)$-$ 式(2.143)

$$\rho_{y -} c_v T_{y -}v_{y -}dzdx - c_v \bigg( \rho_{y -}T_{y -}v_{y -}+ \left. \frac{\partial (\rho T v)}{\partial y} \right\vert _ {{y -}} dy \bigg) dzdx$$

$$= - c_v \left. \frac{\partial (\rho T v)}{\partial y} \right\vert _ {{y -}} dxdydz$$ (2.146)

$z$ 軸に沿った出入

式(2.144)$-$ 式(2.145)

$$\rho_{z -} c_v T_{z -}w_{z -}dxdy - c_v \bigg( \rho_{z -}T_{z -}w_{z -}+ \left. \frac{\partial (\rho T w)}{\partial z} \right\vert _ {{z -}} dz \bigg) dxdy$$

$$= - c_v \left. \frac{\partial (\rho T w)}{\partial z} \right\vert _ {{z -}} dxdydz$$ (2.147)

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

$$- c_v \bigg( \dfrac{\partial (\rho Tu)}{\partial x} + \dfrac{\partial (\rho Tv)}{\partial y} + \dfrac{\partial (\rho Tw)}{\partial z} \bigg) dxdydz$$

$$= - c_v \bigg( \rho T \dfrac{\partial u}{\partial x} + \rho u \df... ...c{\partial T}{\partial z} + T w \dfrac{\partial \rho}{\partial z} \bigg) dxdydz$$

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

$$= - c_v ( \rho T \bm{\nabla} \cdot \bm{v} + \rho \bm{v} \cdot \bm{\nabla}T + T \bm{v} \cdot \bm{\nabla}\rho ) dxdydz$$

$$= - c_v \bm{\nabla} \cdot (\rho T \bm{v}) dxdydz$$ (2.148)

• 学問 -top-

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