微小領域を用いた流線方向のBernoulliの定理の導出
https://kakeru.app/e1e62ede4310622083fd44cf28179984 https://i.kakeru.app/e1e62ede4310622083fd44cf28179984.svg
流管の方向を$ \pmb{e}_s、密度を$ \rho、断面積を$ S(s)、流速を$ u(s,t)とする
質量収支とエネルギー収支の式を立てる
質量収支
流管の側面に垂直な成分は0なので、側面からの流量は0でいい
密度・微小領域は共に時間変化しないので、微小領域内での質量の湧き出しも考えなくていい
$ \therefore \rho S(s)u\mathrm{d}t=\rho S(s+\mathrm{d}s)u(s+\mathrm{d}s)\mathrm{d}t
$ \iff \rho\mathrm{d}(Su)=0\quad ①
エネルギー収支
運動エネルギー変化$ \mathrm{d}K
外力によるエネルギー変化$ \mathrm{d}U
体積力である重力と、面積力である圧力によるエネルギー変化がある 流管側面にかかる圧力は$ \pmb{e}_sと直交するので仕事しない
$ \therefore \frac12(\rho Su\mathrm{d}t)u^2+(\rho Su\mathrm{d}t)gz+PSu\mathrm{d}t=\frac12(\rho S(s+\mathrm{d}s)u(s+\mathrm{d}s)\mathrm{d}t)u(s+\mathrm{d}s)^2+(\rho S(s+\mathrm{d}s)u(s+\mathrm{d}s)\mathrm{d}t)gz(s+\mathrm{d}s)+P(s+\mathrm{d}s)S(s+\mathrm{d}s)u(s+\mathrm{d}s)\mathrm{d}t
$ \iff \frac12\rho\mathrm{d}(Su^3)\mathrm{d}t+\rho g\mathrm{d}(Suz)\mathrm{d}t+\mathrm{d}(PSu)\mathrm{d}t=0\quad ②
なんか変なことしてないかこれ?
どういう立式だよ
①と②を組む
$ ①\land②
$ \iff\begin{dcases}\rho\mathrm{d}(Su)=0\\Su\frac12\rho\mathrm{d}(u^2)\mathrm{d}t+u^2\frac12\rho\mathrm{d}(Su)\mathrm{d}t+\rho gSu\mathrm{d}z\mathrm{d}t+\rho gz\mathrm{d}(Su)\mathrm{d}t+Su\mathrm{d}P\mathrm{d}t+P\mathrm{d}(Su)\mathrm{d}t=0\end{dcases}
$ ①で$ \mathrm{d}(Su)を消す
$ \implies Su\frac12\rho\mathrm{d}(u^2)\mathrm{d}t+\rho gSu\mathrm{d}z\mathrm{d}t+Su\mathrm{d}P\mathrm{d}t=0
$ \iff \frac12\rho\mathrm{d}(u^2)+\rho g\mathrm{d}z+\mathrm{d}P=0
$ Su\neq0として考える
$ \iff \mathrm{d}\left(\frac12\rho u^2+\rho gz+P\right)=0
$ \underline{\iff \frac12\rho u^2+\rho gz+P=\rm Const.\quad}_\blacksquare
ここでは密度を一定として扱ったが、密度一定でなくても成立する
エネルギー収支:$ \frac12\mathrm{d}(\rho Su^3)\mathrm{d}t+g\mathrm{d}(\rho Suz)\mathrm{d}t+\mathrm{d}(PSu)\mathrm{d}t=0
$ \iff\rho Su\frac12\mathrm{d}(u^2)\mathrm{d}t+u^2\frac12\mathrm{d}(\rho Su)\mathrm{d}t+\rho Sug\mathrm{d}z\mathrm{d}t+gz\mathrm{d}(\rho Su)\mathrm{d}t+\rho Su\mathrm{d}\left(\frac{P}{\rho}\right)\mathrm{d}t+\frac{P}{\rho}\mathrm{d}(\rho Su)\mathrm{d}t=0
$ \frac{P}{\rho}が出てしまう
連続の式:$ \mathrm{d}(\rho Su)=0
エネルギー収支に連続の式を代入すると
$ \iff\rho Su\frac12\mathrm{d}(u^2)\mathrm{d}t+\rho Sug\mathrm{d}(z)\mathrm{d}t+\rho Su\mathrm{d}\left(\frac{P}{\rho}\right)\mathrm{d}t=0
$ \iff\frac12\mathrm{d}(u^2)+g\mathrm{d}z+\mathrm{d}\left(\frac{P}{\rho}\right)=0
$ \iff\mathrm{d}\left(\frac12u^2+gz+\frac{P}{\rho}\right)=0