跳到主要内容

東京大学 工学系研究科 機械工学専攻 機械工学第1部 2012年8月実施 問題 2 設問 II

Author​

祭音Myyura (Based on SN's answer refined with GPT 6 Astra)

Description​

重力と流体の圧縮性の影響は無視してよい。

非圧縮性ニュートン流体に関する以下の設問に答えよ。ただし、流体の粘性係数を μ\mu とする。

(1) 単一のまっすぐの円管内の流れを考える。流れを層流とし、十分発達した定常状態にあるとする。このとき、図 2-2 に示すように、半径 rr、長さ Δx\Delta x の円柱型の検査体積を考える。この検査体積に作用する圧力 pa,pbp_a,p_b とせん断応力 τ\tau の関係を、検査体積に運動量法則を適用することにより求めよ。

(2) 設問 (1) において、半径 rr の位置の検査面に働くせん断応力 τ\tau は、円管の流速分布 u(r)u(r) を用いて、次式で与えられる。

τ=μdudr\tau=\mu\frac{du}{dr}

円管の半径を RR として、u(r)u(r) を pa,pb,r,R,μ,Δxp_a,p_b,r,R,\mu,\Delta x により表せ。

(3) 設問 (2) で求めた結果を利用して、円管を流れる流量 QQ を求めよ。

図2-2:円柱型検査体積の原図

次に、設問 I の分岐管に関して、分岐部前後の管摩擦損失の影響が無視できない場合を考える。ただし、ここでは、分岐後の管 2、管 3 の長さは等しいとする。図 2-3 に示すように、分岐部の xx 座標を xBx_B、管端の xx 座標を xTx_T とする。図に示すように位置 xTx_T において、xx 軸から距離 hh だけ離れた位置で、液体は流出しているとする。すなわち、管 2、管 3 の長さは、h2+(xT−xB)2\sqrt{h^2+(x_T-x_B)^2} で与えられるとする。管 1~3 の管内の流れについて、十分発達した定常状態における層流の速度分布を仮定し、分岐部での損失は無視して、以下の設問に答えよ。

設問 I と同様に、液体は円管から大気中へ流出する。管 1 の流入部は x=0x=0、ゲージ圧は p1p_1、管 1~3 の断面積は A1,A2,A3A_1,A_2,A_3 である。

(4) 管 1 の流入部の圧力 p1p_1 を、管 1 への流入流量 Q1Q_1 を用いて表せ。

(5) 図 2-3 において、管 2 と管 3 の断面積が等しい場合を考える(A2=A3≡SA_2=A_3\equiv S)。また、A1>2SA_1>2S とする。このとき、管 1 の流入部、管 2、管 3 の流出部の位置を固定したまま、分岐位置を変化させることを考える。図 2-3 に示した部分における圧力損失を最小にする分岐点の xx 座標 xBx_B を求めよ。

図2-3:分岐位置と固定した流出部の原図

题目描述​

可忽略重力及流体可压缩性的影响。回答关于不可压缩牛顿流体的下列问题,流体的动力黏度为 μ\mu。

(1) 考虑单根直圆管内充分发展的稳态层流。如图 2-2,取半径为 rr、长度为 Δx\Delta x 的圆柱形控制体。对该控制体应用动量定律,求作用于其上的压力 pa,pbp_a,p_b 与剪应力 τ\tau 之间的关系。

(2) 在 (1) 中,半径 rr 处控制面上的剪应力与管内速度分布 u(r)u(r) 满足 τ=μ du/dr\tau=\mu\,du/dr。设管半径为 RR,用 pa,pb,r,R,μ,Δxp_a,p_b,r,R,\mu,\Delta x 表示 u(r)u(r)。

(3) 利用 (2) 的结果求圆管中的流量 QQ。

接着考虑设问 I 的分岔管,在分岔前后各管的沿程摩擦损失不可忽略的情况。此处假设支管 2、3 等长。如图 2-3,分岔点的 xx 坐标为 xBx_B,管出口的 xx 坐标为 xTx_T;两出口均距 xx 轴 hh,即支管 2、3 的长度均为 h2+(xT−xB)2\sqrt{h^2+(x_T-x_B)^2}。假设各管内均具有充分发展的稳态层流速度分布,忽略分岔处局部损失。

与设问 I 相同,液体由圆管流入大气;管 1 入口位于 x=0x=0,表压为 p1p_1;管 1~3 截面积分别为 A1,A2,A3A_1,A_2,A_3。

(4) 用流入管 1 的流量 Q1Q_1 表示管 1 入口压力 p1p_1。

(5) 在图 2-3 中,设管 2、3 的截面积相等,即 A2=A3≡SA_2=A_3\equiv S,且 A1>2SA_1>2S。保持管 1 入口和管 2、3 出口位置固定,改变分岔位置,求使图示部分压力损失最小的分岔点坐标 xBx_B。

Kai​

(1)​

円柱側面に働くせん断力を 2πrΔx τ2\pi r\Delta x\,\tau とすると、

πr2(pa−pb)+2πrΔx τ=0.\pi r^2(p_a-p_b)+2\pi r\Delta x\,\tau=0.

よって

τ=pb−pa2Δxr.\boxed{\tau=\frac{p_b-p_a}{2\Delta x}r}.

(2)​

μ du/dr=τ\mu\,du/dr=\tau を積分し、u(R)=0u(R)=0 を用いると、

u(r)=pa−pb4μΔx(R2−r2).\boxed{u(r)=\frac{p_a-p_b}{4\mu\Delta x}(R^2-r^2)}.

(3)​

Q=2π∫0Ru(r)r dr=πR4(pa−pb)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}}.

断面積 A=πR2A=\pi R^2 の管では、長さ ℓ\ell に対する圧力損失は 8πμℓQ/A28\pi\mu\ell Q/A^2 である。

(4)​

分岐後の各管の長さは ℓ=h2+(xT−xB)2\ell=\sqrt{h^2+(x_T-x_B)^2}。分岐点圧力を pBp_B とすれば、

Q2=A22pB8πμℓ,Q3=A32pB8πμℓ,Q1=Q2+Q3.\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}

これと管 1 の圧力損失を合わせて、

p1=8πμQ1[xBA12+h2+(xT−xB)2A22+A32].\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]}.

(5)​

A2=A3=SA_2=A_3=S を代入して微分すると、

dp1dxB=8πμQ1[1A12−xT−xB2S2h2+(xT−xB)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].

第 2 微分は正なので、内部の極小点は

xB=xT−2S2hA14−4S4.\boxed{x_B=x_T-\frac{2S^2h}{\sqrt{A_1^4-4S^4}}}.

分岐点を 0≤xB≤xT0\le x_B\le x_T に制限し、右辺が負になる場合は境界 xB=0x_B=0 が最適となる。

Reference​