跳到主要内容

名古屋大学 工学研究科 機械航空系 2018年度 専門部門 問題2

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

(1)​

大気に開いた断面積 SA,SBS_A,S_B の容器 A、B を、底部の面積 σ\sigma の小孔で接続する。σ≪SA,SB\sigma\ll S_A,S_B とし、密度一定の非粘性液体が A から B へ流れる。自由表面は水平で、小孔の中心を基準とする高さを zA,zBz_A,z_B とする。初期値は zA(0)=hA>hB=zB(0)z_A(0)=h_A>h_B=z_B(0)、重力加速度は gg。

  1. 微小な水位変化 ΔzA,ΔzB\Delta z_A,\Delta z_B の関係を求めよ。
  2. 水位差を用いて小孔流速 vv を求めよ。
  3. 水位差 h=zA−zBh=z_A-z_B の時間変化を求めよ。

(2)​

水平面から角度 θ\theta 傾いた平板上を、密度 ρ\rho、粘度 μ\mu の非圧縮ニュートン流体が厚さ hh の液膜として定常層流で流下する。斜面下向きを xx、壁面から垂直外向きを yy とし、速度は xx 方向の v(y)v(y) のみである。壁面では滑りなし、自由表面のせん断応力はゼロとする。

  1. せん断応力 τ(y)\tau(y) の微分方程式を導け。
  2. 自由表面速度を V=v(h)V=v(h) として v=V(y/h)(2−y/h)v=V(y/h)(2-y/h) を示せ。
  3. 単位幅当たり流量 QQ を求め、hh と QQ の関係を示せ。

連通容器と斜面上の液膜

题目描述​

(1)两个向大气敞开的容器 A、B 截面积分别为 SA,SBS_A,S_B,由底部面积为 σ≪SA,SB\sigma\ll S_A,S_B 的小孔连通。无黏、恒密度液体从 A 流向 B,液面保持水平;以孔中心为基准,高度为 zA,zBz_A,z_B,初值为 hA>hBh_A>h_B,重力加速度为 gg。求两侧微小液位变化的关系、小孔速度,以及液位差 h=zA−zBh=z_A-z_B 随时间的变化。(2)密度 ρ\rho、黏度 μ\mu 的不可压牛顿流体以厚度 hh 在倾角 θ\theta 的平板上稳态层流流动;xx 沿下坡方向,yy 为向外法向,速度为 v(y)v(y),壁面无滑移,液面剪应力为零。推导剪应力方程,以液面速度 VV 表示速度分布,并求单位宽流量 QQ 及膜厚与流量的关系。

Kai​

(1)​

1), 2)​

質量保存から

SAΔzA+SBΔzB=0.\boxed{S_A\Delta z_A+S_B\Delta z_B=0}.

自由表面の速度を小孔速度に比べて無視し、圧力差 ρg(zA−zB)\rho g(z_A-z_B) を動圧に等置すると

v=2g(zA−zB).\boxed{v=\sqrt{2g(z_A-z_B)}}.

3)​

SAz˙A=−σvS_A\dot z_A=-\sigma v、SBz˙B=σvS_B\dot z_B=\sigma v より

h˙=−σ(1SA+1SB)2gh.\dot h=-\sigma\left(\frac1{S_A}+\frac1{S_B}\right)\sqrt{2gh}.

h>0h>0 の間、変数分離して

h(t)=[hA−hB−σ(SA+SB)SASBg2 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_*.

ここで

t∗=SASBσ(SA+SB)2(hA−hB)g.t_* =\frac{S_AS_B}{\sigma(S_A+S_B)}\sqrt{\frac{2(h_A-h_B)}g}.

この準定常モデルでは t≥t∗t\ge t_* で h=0h=0 である。二乗式を t∗t_* より先まで延長してはならない。

(2)​

1)​

微小要素の xx 方向の力のつり合いより

dτdy+ρgsin⁡θ=0.\boxed{\frac{d\tau}{dy}+\rho g\sin\theta=0}.

2)​

τ=μv′\tau=\mu v'、τ(h)=0\tau(h)=0、v(0)=0v(0)=0 を用いて二回積分する。

τ=ρgsin⁡θ(h−y),v=ρgsin⁡θμ(hy−y22).\tau=\rho g\sin\theta(h-y),\qquad v=\frac{\rho g\sin\theta}{\mu}\left(hy-\frac{y^2}{2}\right).

従って

V=ρgsin⁡θ2μh2,v=Vyh(2−yh).V=\frac{\rho g\sin\theta}{2\mu}h^2,\qquad \boxed{v=V\frac yh\left(2-\frac yh\right)}.

3)​

Q=∫0hv dy=23Vh=ρgsin⁡θ3μh3,\boxed{Q=\int_0^h v\,dy=\frac23Vh=\frac{\rho g\sin\theta}{3\mu}h^3},
h=(3μQρgsin⁡θ)1/3∝Q1/3.\boxed{h=\left(\frac{3\mu Q}{\rho g\sin\theta}\right)^{1/3}\propto Q^{1/3}}.

Reference​