名古屋大学 工学研究科 機械航空系 2017年8月実施 専門部門 問題2
Author
祭音Myyura (co-authored with GPT 6 Astra)
Description
問題2.
(1) 図1に示すように,上部が大気に開放された断面積 S A S_A S A の容器 A と断面積 S B S_B S B の容器 B が断面積 σ ( σ ≪ S A , S B ) \sigma\ (\sigma\ll S_A,S_B) σ ( σ ≪ S A , S B ) の小さな孔のあいた仕切板で隔てられている。容器 A,B には孔中心から測ってそれぞれ高さ z A , z B z_A,z_B z A , z B まで液体が入っている。液体は非圧縮かつ非粘性で密度を ρ \rho ρ とする。重力加速度を鉛直下向きに g g g とする。液面は孔から十分に離れており,液面は水平を保ったまま上下するとみなせる。また,孔断面では流れは一様とみなせるとする。水平方向右向きに正方向をとる。このとき以下の問いに答えよ。
1) 容器 A の液面高さが Δ z A \Delta z_A Δ z A 変化するとき,容器 B の液面高さの変化 Δ z B \Delta z_B Δ z B を求めよ。
2) z A > z B z_A>z_B z A > z B のとき,孔を通過する液体の速度を求めよ。
3) 時刻 t ( > 0 ) t\ (>0) t ( > 0 ) での容器 A,B の液面高さの差 z A ( t ) − z B ( t ) z_A(t)-z_B(t) z A ( t ) − z B ( t ) を求めよ。ただし,t = 0 t=0 t = 0 での容器 A,B の液面高さをそれぞれ z A ( 0 ) = h A , z B ( 0 ) = h B ( h A > h B ) z_A(0)=h_A,\ z_B(0)=h_B\ (h_A>h_B) z A ( 0 ) = h A , z B ( 0 ) = h B ( h A > h B ) とする。
(2) 図2に示すように,水平方向から角度 θ \theta θ 傾いた斜面を考える。斜面に沿って下向きに x x x 軸,斜面に対して垂直上向きに y y y 軸をとり,重力加速度 g g g が鉛直下向きに働くとする。液体が一定の厚さ h h h および体積流量 Q Q Q で層流状態を保ち,x x x 軸に平行な速度成分のみをもって流れ落ちている。液体は密度 ρ \rho ρ ,粘度 μ \mu μ の非圧縮粘性流体(ニュートン流体)とする。図に示す2次元平面での運動のみを扱うこととし,y = h y=h y = h の面は自由表面であるとする。このとき以下の問いに答えよ。
1) 液体に作用する重力と液体中の摩擦応力 τ \tau τ のつり合いを考え,d τ / d y d\tau/dy d τ / d y を求めよ。
2) 図2のように自由表面での液体の流速を v ( h ) = V v(h)=V v ( h ) = V とする。流速分布 v ( y ) v(y) v ( y ) が,
d τ d y = μ d 2 v d y 2 \frac{d\tau}{dy}=\mu\frac{d^2v}{dy^2} d y d τ = μ d y 2 d 2 v
の微分方程式に従うとき,
v ( y ) = V y h ( 2 − y h ) v(y)=V\frac yh\left(2-\frac yh\right) v ( y ) = V h y ( 2 − h y )
と表されることを示せ。
3) 速度分布を y y y 方向に関して積分することで紙面垂直方向単位長さあたりの体積流量 Q Q Q を求めよ。また液体層の厚さ h h h は体積流量 Q Q Q の何乗に比例するか答えよ。
题目描述
问题2.
(1) 如图1,上方与大气相通的容器 A、B 的截面积分别为 S A , S B S_A,S_B S A , S B ,两容器之间由开有小孔的隔板隔开,孔面积为 σ \sigma σ ,且 σ ≪ S A , S B \sigma\ll S_A,S_B σ ≪ S A , S B 。从孔中心计,容器 A、B 中液体高度分别为 z A , z B z_A,z_B z A , z B 。液体不可压缩、无黏,密度为 ρ \rho ρ ;重力加速度 g g g 竖直向下。液面距孔足够远,可视为始终保持水平地上下移动;孔截面内的流动可视为均匀。取水平向右为正方向。回答下列问题。
1) 当容器 A 的液面高度变化 Δ z A \Delta z_A Δ z A 时,求容器 B 液面高度的变化 Δ z B \Delta z_B Δ z B 。
2) 当 z A > z B z_A>z_B z A > z B 时,求流过孔的液体速度。
3) 求 t > 0 t>0 t > 0 时两容器液面高度差 z A ( t ) − z B ( t ) z_A(t)-z_B(t) z A ( t ) − z B ( t ) 。初始 t = 0 t=0 t = 0 时,液面高度分别为 z A ( 0 ) = h A , z B ( 0 ) = h B z_A(0)=h_A,\ z_B(0)=h_B z A ( 0 ) = h A , z B ( 0 ) = h B ,且 h A > h B h_A>h_B h A > h B 。
(2) 如图2,考虑与水平方向夹角为 θ \theta θ 的斜面。沿斜面向下取 x x x 轴,垂直斜面向上取 y y y 轴,重力加速度 g g g 竖直向下。液体保持恒定厚度 h h h 、体积流量 Q Q Q 及层流状态,仅具有平行于 x x x 轴的速度分量并沿斜面流下。液体为密度 ρ \rho ρ 、黏度 μ \mu μ 的不可压缩黏性牛顿流体。只讨论图示二维平面内的运动,y = h y=h y = h 为自由表面。回答下列问题。
1) 考虑液体所受重力与液体内剪应力 τ \tau τ 的平衡,求 d τ / d y d\tau/dy d τ / d y 。
2) 如图2,设自由表面液体流速为 v ( h ) = V v(h)=V v ( h ) = V 。当速度分布 v ( y ) v(y) v ( y ) 满足微分方程
d τ d y = μ d 2 v d y 2 , \frac{d\tau}{dy}=\mu\frac{d^2v}{dy^2}, d y d τ = μ d y 2 d 2 v ,
证明其可表示为
v ( y ) = V y h ( 2 − y h ) . v(y)=V\frac yh\left(2-\frac yh\right). v ( y ) = V h y ( 2 − h y ) .
3) 沿 y y y 方向积分速度分布,求垂直纸面方向单位长度上的体积流量 Q Q Q ,并回答液层厚度 h h h 与体积流量 Q Q Q 的多少次方成正比。
Kai
(1)
1), 2)
質量保存から
S A Δ z A + S B Δ z B = 0 . \boxed{S_A\Delta z_A+S_B\Delta z_B=0}. S A Δ z A + S B Δ z B = 0 .
自由表面の速度を小孔速度に比べて無視し、圧力差 ρ g ( z A − z B ) \rho g(z_A-z_B) ρ g ( z A − z B ) を動圧に等置すると
v = 2 g ( z A − z B ) . \boxed{v=\sqrt{2g(z_A-z_B)}}. v = 2 g ( z A − z B ) .
S A z ˙ A = − σ v S_A\dot z_A=-\sigma v S A z ˙ A = − σ v 、S B z ˙ B = σ v S_B\dot z_B=\sigma v S B z ˙ B = σ v より
h ˙ = − σ ( 1 S A + 1 S B ) 2 g h . \dot h=-\sigma\left(\frac1{S_A}+\frac1{S_B}\right)\sqrt{2gh}. h ˙ = − σ ( S A 1 + S B 1 ) 2 g h .
h > 0 h>0 h > 0 の間、変数分離して
h ( t ) = [ h A − h B − σ ( S A + S B ) S A S B g 2 t ] 2 , 0 ≤ t ≤ t ∗ . \boxed{h(t)=\left[\sqrt{h_A-h_B}-\frac{\sigma(S_A+S_B)}{S_AS_B}\sqrt{\frac g2}\,t\right]^2},
\quad0\le t\le t_*. h ( t ) = [ h A − h B − S A S B σ ( S A + S B ) 2 g t ] 2 , 0 ≤ t ≤ t ∗ .
ここで
t ∗ = S A S B σ ( S A + S B ) 2 ( h A − h B ) g . t_* =\frac{S_AS_B}{\sigma(S_A+S_B)}\sqrt{\frac{2(h_A-h_B)}g}. t ∗ = σ ( S A + S B ) S A S B g 2 ( h A − h B ) .
この準定常モデルでは t ≥ t ∗ t\ge t_* t ≥ t ∗ で h = 0 h=0 h = 0 である。二乗式を t ∗ t_* t ∗ より先まで延長してはならない。
(2)
微小要素の x x x 方向の力のつり合いより
d τ d y + ρ g sin θ = 0 . \boxed{\frac{d\tau}{dy}+\rho g\sin\theta=0}. d y d τ + ρ g sin θ = 0 .
τ = μ v ′ \tau=\mu v' τ = μ v ′ 、τ ( h ) = 0 \tau(h)=0 τ ( h ) = 0 、v ( 0 ) = 0 v(0)=0 v ( 0 ) = 0 を用いて二回積分する。
τ = ρ g sin θ ( h − y ) , v = ρ g sin θ μ ( h y − y 2 2 ) . \tau=\rho g\sin\theta(h-y),\qquad
v=\frac{\rho g\sin\theta}{\mu}\left(hy-\frac{y^2}{2}\right). τ = ρ g sin θ ( h − y ) , v = μ ρ g sin θ ( h y − 2 y 2 ) .
従って
V = ρ g sin θ 2 μ h 2 , v = V y h ( 2 − y h ) . V=\frac{\rho g\sin\theta}{2\mu}h^2,\qquad
\boxed{v=V\frac yh\left(2-\frac yh\right)}. V = 2 μ ρ g sin θ h 2 , v = V h y ( 2 − h y ) .
Q = ∫ 0 h v d y = 2 3 V h = ρ g sin θ 3 μ h 3 , \boxed{Q=\int_0^h v\,dy=\frac23Vh=\frac{\rho g\sin\theta}{3\mu}h^3}, Q = ∫ 0 h v d y = 3 2 Vh = 3 μ ρ g sin θ h 3 ,
h = ( 3 μ Q ρ g sin θ ) 1 / 3 ∝ Q 1 / 3 . \boxed{h=\left(\frac{3\mu Q}{\rho g\sin\theta}\right)^{1/3}\propto Q^{1/3}}. h = ( ρ g sin θ 3 μ Q ) 1/3 ∝ Q 1/3 .
Reference