東京大学 工学系研究科 機械工学専攻 機械工学第1部 2012年8月実施 問題 2 設問 II
Author
祭音Myyura (Based on SN's answer refined with GPT 6 Astra)
Description
重力と流体の圧縮性の影響は無視してよい。
非圧縮性ニュートン流体に関する以下の設問に答えよ。ただし、流体の粘性係数を μ \mu μ とする。
(1) 単一のまっすぐの円管内の流れを考える。流れを層流とし、十分発達した定常状態にあるとする。このとき、図 2-2 に示すように、半径 r r r 、長さ Δ x \Delta x Δ x の円柱型の検査体積を考える。この検査体積に作用する圧力 p a , p b p_a,p_b p a , p b とせん断応力 τ \tau τ の関係を、検査体積に運動量法則を適用することにより求めよ。
(2) 設問 (1) において、半径 r r r の位置の検査面に働くせん断応力 τ \tau τ は、円管の流速分布 u ( r ) u(r) u ( r ) を用いて、次式で与えられる。
τ = μ d u d r \tau=\mu\frac{du}{dr} τ = μ d r d u
円管の半径を R R R として、u ( r ) u(r) u ( r ) を p a , p b , r , R , μ , Δ x p_a,p_b,r,R,\mu,\Delta x p a , p b , r , R , μ , Δ x により表せ。
(3) 設問 (2) で求めた結果を利用して、円管を流れる流量 Q Q Q を求めよ。
次に、設問 I の分岐管に関して、分岐部前後の管摩擦損失の影響が無視できない場合を考える。ただし、ここでは、分岐後の管 2、管 3 の長さは等しいとする。図 2-3 に示すように、分岐部の x x x 座標を x B x_B x B 、管端の x x x 座標を x T x_T x T とする。図に示すように位置 x T x_T x T において、x x x 軸から距離 h h h だけ離れた位置で、液体は流出しているとする。すなわち、管 2、管 3 の長さは、h 2 + ( x T − x B ) 2 \sqrt{h^2+(x_T-x_B)^2} h 2 + ( x T − x B ) 2 で与えられるとする。管 1~3 の管内の流れについて、十分発達した定常状態における層流の速度分布を仮定し、分岐部での損失は無視して、以下の設問に答えよ。
設問 I と同様に、液体は円管から大気中へ流出する。管 1 の流入部は x = 0 x=0 x = 0 、ゲージ圧は p 1 p_1 p 1 、管 1~3 の断面積は A 1 , A 2 , A 3 A_1,A_2,A_3 A 1 , A 2 , A 3 である。
(4) 管 1 の流入部の圧力 p 1 p_1 p 1 を、管 1 への流入流量 Q 1 Q_1 Q 1 を用いて表せ。
(5) 図 2-3 において、管 2 と管 3 の断面積が等しい場合を考える(A 2 = A 3 ≡ S A_2=A_3\equiv S A 2 = A 3 ≡ S )。また、A 1 > 2 S A_1>2S A 1 > 2 S とする。このとき、管 1 の流入部、管 2、管 3 の流出部の位置を固定したまま、分岐位置を変化させることを考える。図 2-3 に示した部分における圧力損失を最小にする分岐点の x x x 座標 x B x_B x B を求めよ。
题目描述
可忽略重力及流体可压缩性的影响。回答关于不可压缩牛顿流体的下列问题,流体的动力黏度为 μ \mu μ 。
(1) 考虑单根直圆管内充分发展的稳态层流。如图 2-2,取半径为 r r r 、长度为 Δ x \Delta x Δ x 的圆柱形控制体。对该控制体应用动量定律,求作用于其上的压力 p a , p b p_a,p_b p a , p b 与剪应力 τ \tau τ 之间的关系。
(2) 在 (1) 中,半径 r r r 处控制面上的剪应力与管内速度分布 u ( r ) u(r) u ( r ) 满足 τ = μ d u / d r \tau=\mu\,du/dr τ = μ d u / d r 。设管半径为 R R R ,用 p a , p b , r , R , μ , Δ x p_a,p_b,r,R,\mu,\Delta x p a , p b , r , R , μ , Δ x 表示 u ( r ) u(r) u ( r ) 。
(3) 利用 (2) 的结果求圆管中的流量 Q Q Q 。
接着考虑设问 I 的分岔管,在分岔前后各管的沿程摩擦损失不可忽略的情况。此处假设支管 2、3 等长。如图 2-3,分岔点的 x x x 坐标为 x B x_B x B ,管出口的 x x x 坐标为 x T x_T x T ;两出口均距 x x x 轴 h h h ,即支管 2、3 的长度均为 h 2 + ( x T − x B ) 2 \sqrt{h^2+(x_T-x_B)^2} h 2 + ( x T − x B ) 2 。假设各管内均具有充分发展的稳态层流速度分布,忽略分岔处局部损失。
与设问 I 相同,液体由圆管流入大气;管 1 入口位于 x = 0 x=0 x = 0 ,表压为 p 1 p_1 p 1 ;管 1~3 截面积分别为 A 1 , A 2 , A 3 A_1,A_2,A_3 A 1 , A 2 , A 3 。
(4) 用流入管 1 的流量 Q 1 Q_1 Q 1 表示管 1 入口压力 p 1 p_1 p 1 。
(5) 在图 2-3 中,设管 2、3 的截面积相等,即 A 2 = A 3 ≡ S A_2=A_3\equiv S A 2 = A 3 ≡ S ,且 A 1 > 2 S A_1>2S A 1 > 2 S 。保持管 1 入口和管 2、3 出口位置固定,改变分岔位置,求使图示部分压力损失最小的分岔点坐标 x B x_B x B 。
Kai
(1)
円柱側面に働くせん断力を 2 π r Δ x τ 2\pi r\Delta x\,\tau 2 π r Δ x τ とすると、
π r 2 ( p a − p b ) + 2 π r Δ x τ = 0. \pi r^2(p_a-p_b)+2\pi r\Delta x\,\tau=0. π r 2 ( p a − p b ) + 2 π r Δ x τ = 0.
よって
τ = p b − p a 2 Δ x r . \boxed{\tau=\frac{p_b-p_a}{2\Delta x}r}. τ = 2Δ x p b − p a r .
(2)
μ d u / d r = τ \mu\,du/dr=\tau μ d u / d r = τ を積分し、u ( R ) = 0 u(R)=0 u ( R ) = 0 を用いると、
u ( r ) = p a − p b 4 μ Δ x ( R 2 − r 2 ) . \boxed{u(r)=\frac{p_a-p_b}{4\mu\Delta x}(R^2-r^2)}. u ( r ) = 4 μ Δ x p a − p b ( R 2 − r 2 ) .
(3)
Q = 2 π ∫ 0 R u ( r ) r d r = π R 4 ( p a − p b ) 8 μ Δ x . \boxed{Q=2\pi\int_0^R u(r)r\,dr=\frac{\pi R^4(p_a-p_b)}{8\mu\Delta x}}. Q = 2 π ∫ 0 R u ( r ) r d r = 8 μ Δ x π R 4 ( p a − p b ) .
断面積 A = π R 2 A=\pi R^2 A = π R 2 の管では、長さ ℓ \ell ℓ に対する圧力損失は 8 π μ ℓ Q / A 2 8\pi\mu\ell Q/A^2 8 π μ ℓ Q / A 2 である。
(4)
分岐後の各管の長さは ℓ = h 2 + ( x T − x B ) 2 \ell=\sqrt{h^2+(x_T-x_B)^2} ℓ = h 2 + ( x T − x B ) 2 。分岐点圧力を p B p_B p B とすれば、
Q 2 = A 2 2 p B 8 π μ ℓ , Q 3 = A 3 2 p B 8 π μ ℓ , Q 1 = Q 2 + Q 3 . \begin{gathered}
Q_2=\frac{A_2^2p_B}{8\pi\mu\ell},\\
Q_3=\frac{A_3^2p_B}{8\pi\mu\ell},\\
Q_1=Q_2+Q_3.
\end{gathered} Q 2 = 8 π μ ℓ A 2 2 p B , Q 3 = 8 π μ ℓ A 3 2 p B , Q 1 = Q 2 + Q 3 .
これと管 1 の圧力損失を合わせて、
p 1 = 8 π μ Q 1 [ x B A 1 2 + h 2 + ( x T − x B ) 2 A 2 2 + A 3 2 ] . \boxed{p_1=8\pi\mu Q_1\left[\frac{x_B}{A_1^2}+\frac{\sqrt{h^2+(x_T-x_B)^2}}{A_2^2+A_3^2}\right]}. p 1 = 8 π μ Q 1 [ A 1 2 x B + A 2 2 + A 3 2 h 2 + ( x T − x B ) 2 ] .
(5)
A 2 = A 3 = S A_2=A_3=S A 2 = A 3 = S を代入して微分すると、
d p 1 d x B = 8 π μ Q 1 [ 1 A 1 2 − x T − x B 2 S 2 h 2 + ( x T − x B ) 2 ] . \frac{dp_1}{dx_B}=8\pi\mu Q_1\left[\frac1{A_1^2}-\frac{x_T-x_B}{2S^2\sqrt{h^2+(x_T-x_B)^2}}\right]. d x B d p 1 = 8 π μ Q 1 [ A 1 2 1 − 2 S 2 h 2 + ( x T − x B ) 2 x T − x B ] .
第 2 微分は正なので、内部の極小点は
x B = x T − 2 S 2 h A 1 4 − 4 S 4 . \boxed{x_B=x_T-\frac{2S^2h}{\sqrt{A_1^4-4S^4}}}. x B = x T − A 1 4 − 4 S 4 2 S 2 h .
分岐点を 0 ≤ x B ≤ x T 0\le x_B\le x_T 0 ≤ x B ≤ x T に制限し、右辺が負になる場合は境界 x B = 0 x_B=0 x B = 0 が最適となる。
Reference