• 学問 -top-

2.3.3.2 非圧縮性流体（密度$\rho$ [kg/m$^3$ ]は一定）

密度変化を考慮しないため、質量保存式（連続の式(2.38) $^{\text{p.\pageref{eq-mass}}}$ ）より $\nabla \cdot \bm{v} = \dfrac{\partial u}{\partial x} + \dfrac{\partial v}{\partial y} + \dfrac{\partial w}{\partial z} = 0$ となり、式(2.65) $^{\text{p.\pageref{eq-sigma-com}}}$ は次のように表される。

$$\left. \begin{array}{ccc} \sigma_x & = & - P + 2 \mu \dfrac{\part... ... \sigma_z & = & - P + 2 \mu \dfrac{\partial w}{\partial z} \end{array} \right\}$$ (2.76)

$\sigma$ の式(2.77) $^{\text{p.\pageref{eq-sigma}}}$ と$\tau$ の式(2.66) $^{\text{p.\pageref{eq-tau}}}$ から、それぞれの面に掛かる力を求める。力は外への方向が正となっている。

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

$$\left( \begin{array}{c} \sigma_{x, {x -}} dydz \\ \tau_{xy, {x -}... ...c{\partial w}{\partial x} \right\vert _ {{x -}} \bigg) \end{array} \right) dydz$$ (2.77)

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

$$\left( \begin{array}{c} \sigma_{x, {x +}} dydz \\ \tau_{xy, {x +}... ...rtial x} \right\vert _ {{x +}} \bigg) \vspace{.5em} \\ \end{array} \right) dydz$$

$$= \left( \begin{array}{c} - \bigg( P_{x -} + \left. \dfrac{\parti... ...x^2 } \right\vert _ {{x -}} dx \bigg) \vspace{.5em} \\ \end{array} \right) dydz$$ (2.78)

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

$$\left( \begin{array}{c} \tau_{yx, {y -}} dzdx \vspace{.5em} \\ \s... ...c{\partial w}{\partial y} \right\vert _ {{y -}} \bigg) \end{array} \right) dzdx$$ (2.79)

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

$$\left( \begin{array}{c} \tau_{yx, {y +}} dzdx \vspace{.5em} \\ \s... ...rtial y} \right\vert _ {{y +}} \bigg) \vspace{.5em} \\ \end{array} \right) dzdx$$

$$= \left( \begin{array}{c} \mu \bigg( \left. \dfrac{\partial u}{\p... ...y^2 } \right\vert _ {{y -}} dy \bigg) \vspace{.5em} \\ \end{array} \right) dzdx$$ (2.80)

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

$$\left( \begin{array}{c} \tau_{zx, {z -}} dxdy \vspace{.5em} \\ \t... ...dfrac{\partial w}{\partial z} \right\vert _ {{z -}} \\ \end{array} \right) dxdy$$ (2.81)

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

$$\left( \begin{array}{c} \tau_{zx, {z +}} dxdy \vspace{.5em} \\ \t... ... w}{\partial z} \right\vert _ {{z +}} \vspace{.5em} \\ \end{array} \right) dxdy$$

$$= \left( \begin{array}{c} \mu \bigg( \left. \dfrac{\partial u}{\p... ...z^2 } \right\vert _ {{z -}} dz \bigg) \vspace{.5em} \\ \end{array} \right) dxdy$$ (2.82)

上六式から、まずそれぞれの軸に沿った作用を求める。$xyz$ の各軸にそってコントロールボリューム内を正の方向へ加速させる力を正、負の方向へ加速させる力を負となるように符号を加え、向かい合う面を足し合わせる^2.14。

$x$ 軸に沿った作用

$-$ 式(2.78)$+$ 式(2.79)

$$\left( \begin{array}{c} - \sigma_{x, {x -}} + \sigma_{x, {x +}} \... ...xy, {x +}} \\ - \tau_{xz, {x -}} + \tau_{xz, {x +}} \\ \end{array} \right) dydz$$

$$= \left( \begin{array}{c} - \left\{ - P_{x -}+ 2 \mu \left. \dfra... ...ght\vert _ {{x -}} dx \bigg) \right\} \vspace{.5em} \\ \end{array} \right) dydz$$

$$= \left( \begin{array}{c} - \left. \dfrac{\partial P}{\partial x}... ...artial x^2 } \right\vert _ {{x -}} dx \vspace{.5em} \\ \end{array} \right) dydz$$

$$= \left( \begin{array}{c} - \left. \dfrac{\partial P}{\partial x}... ... x^2 } \right\vert _ {{x -}} \bigg) \vspace{.5em} \\ \end{array} \right) dxdydz$$ (2.83)

$y$ 軸に沿った作用

$-$ 式(2.80)$+$ 式(2.81)

$$\left( \begin{array}{c} - \tau_{yx, {y -}} + \tau_{yx, {y +}} \\ ... ... y^2 } \right\vert _ {{y -}} \bigg) \vspace{.5em} \\ \end{array} \right) dxdydz$$ (2.84)

$z$ 軸に沿った作用

$-$ 式(2.82)$+$ 式(2.83)

$$\left( \begin{array}{c} - \tau_{zx, {z -}} + \tau_{zx, {z +}} \\ ... ...partial z^2 } \right\vert _ {{z -}} \vspace{.5em} \\ \end{array} \right) dxdydz$$ (2.85)

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

$$\left( \begin{array}{c} - \dfrac{\partial P}{\partial x} + 2 \mu ... ...dfrac{\partial ^2 w}{\partial z^2 } \vspace{.5em} \\ \end{array} \right) dxdydz$$

$$= \left( \begin{array}{c} - \dfrac{\partial P}{\partial x} + \mu ... ...partial ^2 w}{\partial z^2 } \bigg) \vspace{.5em} \\ \end{array} \right) dxdydz$$

$$= \left( \begin{array}{c} - \dfrac{\partial P}{\partial x} + \mu ... ...\mu \bm{\nabla} \cdot \bm{\nabla} w \vspace{.5em} \\ \end{array} \right) dxdydz$$

$$= \left\{ - \bm{\nabla} P + \mu \bm{\nabla} (\bm{\nabla} \cdot \bm{v}) + \mu \bm{\nabla} \cdot \bm{\nabla} \bm{v} \right\} dxdydz$$

非圧縮性流体（密度$\rho$ [kg/m$^3$ ]は一定）では、式(2.38) $^{\text{p.\pageref{eq-mass}}}$ より $\bm{\nabla} \cdot \bm{v} = 0$ であるので、次式となる。

$$\left( - \bm{\nabla} P + \mu \bm{\nabla} \cdot \bm{\nabla} \bm{v} \right) dxdydz$$ (2.86)

• 学問 -top-

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