114 年 國立臺灣大學應用力學研究所乙組《流體力學(A)》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊

第 1 題

  1. Water enters a rotating lawn sprinkler through its base at the steady rate of 1000 ml/s as sketched in the following figure. The exit area of each of the two nozzles is 30 mm². Determine the average speed of the water (m/s) leaving the nozzle, relative to the nozzle, if

(a) (15%) the rotary sprinkler head is stationary,
(b) (5%) the sprinkler head rotates at 600 rpm, and
(c) (5%) the sprinkler head accelerates from 0 to 600 rpm.

🖼️【此處有附圖,請對照原卷】The figure shows a rotating lawn sprinkler. Water enters from the base (Section 1) with a flow rate Q=1000Q = 1000 ml/s. The sprinkler head has two nozzles, each with an exit area A2=30A_2 = 30 mm². Section (3) is the control volume. W2W_2 is the velocity of the water leaving the nozzle relative to the nozzle. The sprinkler head rotates.

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查不可壓縮流體的連續方程,以及「相對於運動噴嘴的速度」定義。

對固定在灑水器上的控制體積,噴嘴出口的質量流率應以相對速度 W2W_2 計算:

Qout=2A2W2Q_{\text{out}}=2A_2W_2

由於水不可壓縮,且假設灑水器內沒有流體累積、兩個相同噴嘴平均分流,因此入口流量等於兩個出口的總流量。

解題方法

入口流量為

Q=1000 mL/s=1000×10−6=1.0×10−3 m3/sQ=1000\ \text{mL/s} =1000\times10^{-6} =1.0\times10^{-3}\ \text{m}^3/\text{s}

每個噴嘴的出口面積為

A2=30 mm2=30×10−6=3.0×10−5 m2A_2=30\ \text{mm}^2 =30\times10^{-6} =3.0\times10^{-5}\ \text{m}^2

兩個噴嘴的總出口面積為

2A2=6.0×10−5 m22A_2=6.0\times10^{-5}\ \text{m}^2

因此噴嘴相對速度為

W2=Q2A2=1.0×10−36.0×10−5=16.67 m/sW_2=\frac{Q}{2A_2} =\frac{1.0\times10^{-3}} {6.0\times10^{-5}} =16.67\ \text{m/s}

旋轉只會改變水相對地面的絕對速度:

V2=W2+U2\boldsymbol{V}_2=\boldsymbol{W}_2+\boldsymbol{U}_2

其中 U2\boldsymbol{U}_2 是噴嘴本身的切向速度;但題目要求的是相對於噴嘴的 W2W_2,故旋轉角速度不影響上述連續方程的結果。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 2 題

  1. The two-dimensional flow of a non-viscous, incompressible fluid in the vicinity of the 90° corner as shown in the following figure is described by the stream function ψ=2r2sin⁡2θ\psi = 2r^2 \sin 2\theta, where ψ\psi has units of m²/s when rr is in meters. Assume the fluid density is ρ=103\rho = 10^3 kg/m³ and the x-y plane is horizontal, that is, there is no difference in elevation between points (1) and (2).

(a) (15%) Determine the velocity potential function ϕ\phi.
(b) (10%) If the pressure at point (1) on the wall is 30 kPa, what is the pressure at point (2)?

🖼️【此處有附圖,請對照原卷】The figure shows a 90° corner (in the first quadrant, with the corner at the origin). Streamlines are shown in polar coordinates (r,θ)(r, \theta). Point (1) is on the wall at (r=1,θ=0)(r=1, \theta=0). Point (2) is at (r=0.5,θ=π/4)(r=0.5, \theta=\pi/4). The distance scale is indicated: 1 m on the x-axis, 0.5 m on the y-axis. The streamlines are curved lines originating from the corner.

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

二維不可壓縮流的極座標流函數與速度分量關係為

ur=1r∂ψ∂θ,uθ=−∂ψ∂r.u_r=\frac{1}{r}\frac{\partial \psi}{\partial \theta}, \qquad u_\theta=-\frac{\partial \psi}{\partial r}.

速度勢函數 ϕ\phi 則滿足

ur=∂ϕ∂r,uθ=1r∂ϕ∂θ.u_r=\frac{\partial \phi}{\partial r}, \qquad u_\theta=\frac{1}{r}\frac{\partial \phi}{\partial \theta}.

壓力部分使用定常、無黏、不可壓縮流的伯努力方程:

p+12ρV2+ρgz=常數.p+\frac{1}{2}\rho V^2+\rho gz=\text{常數}.

因兩點位於同一水平面,故 z1=z2z_1=z_2。


(a) 速度勢函數

已知

ψ=2r2sin⁡2θ.\psi=2r^2\sin 2\theta.

先求速度分量:

∂ψ∂θ=4r2cos⁡2θ,\frac{\partial \psi}{\partial \theta} =4r^2\cos 2\theta,

因此

ur=1r∂ψ∂θ=4rcos⁡2θ.u_r=\frac{1}{r}\frac{\partial \psi}{\partial \theta} =4r\cos 2\theta.

又

∂ψ∂r=4rsin⁡2θ,\frac{\partial \psi}{\partial r} =4r\sin 2\theta,

所以

uθ=−∂ψ∂r=−4rsin⁡2θ.u_\theta=-\frac{\partial \psi}{\partial r} =-4r\sin 2\theta.

由 ur=∂ϕ/∂ru_r=\partial\phi/\partial r,

∂ϕ∂r=4rcos⁡2θ.\frac{\partial \phi}{\partial r}=4r\cos 2\theta.

對 rr 積分得

ϕ=2r2cos⁡2θ+f(θ).\phi=2r^2\cos 2\theta+f(\theta).

再利用

uθ=1r∂ϕ∂θ,u_\theta=\frac{1}{r}\frac{\partial\phi}{\partial\theta},

可得

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 3 題

  1. Consider an incompressible Newtonian fluid with dynamic viscosity μ\mu and density ρ\rho, which is flowing steadily through a horizontal long circular straight pipe with diameter DD at an average speed VV. We can express the pressure drop per unit axial length, Δp/Δx\Delta p/\Delta x, as Δp/Δx=f(D,ρ,μ,V)\Delta p/\Delta x = f (D, \rho, \mu, V).

(a) (10%) Show that Δp/Δx\Delta p/\Delta x can be expressed in dimensionless form as ΔpρV2/D=f(Re)\frac{\Delta p}{\rho V^2 / D} = f(\text{Re}), where Re=ρVDμ\text{Re} = \frac{\rho V D}{\mu} is the Reynolds number, through the use of dimensional analysis and Buckingham pi theorem.
(b) (15%) Considered the flow as fully-developed and axisymmetric. The velocity field can be expressed as u=u(r)eru = u(r) e_r in a cylindrical coordinates system, (r,θ,x)(r, \theta, x), fixed with respect to the pipe, with the coordinate center at the central point of the inlet plane, and x-axis along the pipe axis. According to the momentum balance, we can write:
μ(d2udr2+1rdudr)=dpdx\mu \left( \frac{d^2 u}{dr^2} + \frac{1}{r} \frac{du}{dr} \right) = \frac{dp}{dx}
Solve the above equation for u(r)u(r) subject to u(D/2)=0u(D/2)=0, and derive that ΔpρV2/D=32Re\frac{\Delta p}{\rho V^2 / D} = \frac{32}{\text{Re}}.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • Buckingham Π\Pi 定理與量綱分析。
  • Reynolds number 的定義:
Re=ρVDμ\mathrm{Re}=\frac{\rho V D}{\mu}
  • 圓管內穩態、不可壓縮、牛頓流體的 fully developed laminar flow。
  • 軸向動量方程與 no-slip 邊界條件。
  • Hagen–Poiseuille 速度分布及其壓降關係。

題目中速度場寫成 u=u(r)eru=u(r)\boldsymbol e_r,但圓管內的軸向流速應沿管軸方向,因此以下以

v=u(r)ex\boldsymbol v=u(r)\boldsymbol e_x

進行推導。


(a)利用 Buckingham Π\Pi 定理進行量綱分析

已知壓力梯度可表示為

ΔpΔx=f(D,ρ,μ,V)\frac{\Delta p}{\Delta x}=f(D,\rho,\mu,V)

涉及的物理量共有五個:

ΔpΔx,D,ρ,μ,V\frac{\Delta p}{\Delta x},\quad D,\quad \rho,\quad \mu,\quad V

其基本量綱為質量 MM、長度 LL、時間 TT。

各物理量的量綱如下:

[ΔpΔx]=ML−1T−2L=ML−2T−2\left[\frac{\Delta p}{\Delta x}\right] = \frac{ML^{-1}T^{-2}}{L} = ML^{-2}T^{-2} [D]=L,[ρ]=ML−3[D]=L,\qquad [\rho]=ML^{-3} [μ]=ML−1T−1,[V]=LT−1[\mu]=ML^{-1}T^{-1},\qquad [V]=LT^{-1}

由 Buckingham Π\Pi 定理,若有 n=5n=5 個物理量、k=3k=3 個基本量綱,則獨立無因次群的數目為

n−k=5−3=2n-k=5-3=2

選取 ρ\rho、VV、DD 作為重複變數。

第一個無因次群

令

Π1=ΔpΔxρaVbDc\Pi_1= \frac{\Delta p}{\Delta x}\rho^a V^bD^c

代入量綱:

[Π1]=(ML−2T−2)(ML−3)a(LT−1)bLc[\Pi_1] = (ML^{-2}T^{-2}) (ML^{-3})^a (LT^{-1})^b L^c

比較 MM、LL、TT 的次方:

M1+aL−2−3a+b+cT−2−bM^{1+a}L^{-2-3a+b+c}T^{-2-b}

要求其為無因次,因此

1+a=0,−2−b=0,−2−3a+b+c=01+a=0,\qquad -2-b=0,\qquad -2-3a+b+c=0

解得

a=−1,b=−2,c=1a=-1,\qquad b=-2,\qquad c=1

所以

Π1=(ΔpΔx)DρV2=Δp/ΔxρV2/D\Pi_1= \frac{\left(\dfrac{\Delta p}{\Delta x}\right)D}{\rho V^2} = \frac{\Delta p/\Delta x}{\rho V^2/D}

第二個無因次群

令

Π2=μρaVbDc\Pi_2=\mu\rho^aV^bD^c

同樣比較量綱,得到

a=−1,b=−1,c=−1a=-1,\qquad b=-1,\qquad c=-1

因此

Π2=μρVD=1Re\Pi_2=\frac{\mu}{\rho V D} =\frac{1}{\mathrm{Re}}

其中

Re=ρVDμ\mathrm{Re}=\frac{\rho V D}{\mu}

由 Buckingham Π\Pi 定理,兩個無因次群之間存在函數關係:

Π1=F(Π2)\Pi_1=F(\Pi_2)

因此可寫成

Δp/ΔxρV2/D=F(1Re)\frac{\Delta p/\Delta x}{\rho V^2/D} = F\left(\frac{1}{\mathrm{Re}}\right)

將函數重新定義後,即得

Δp/ΔxρV2/D=f(Re)\boxed{ \frac{\Delta p/\Delta x}{\rho V^2/D} = f(\mathrm{Re}) }

(b)求速度分布與壓降關係

令圓管半徑為

R=D2R=\frac{D}{2}

完全發展且軸對稱時,軸向速度只與半徑 rr 有關:

u=u(r)u=u(r)

題目給定軸向動量方程:

μ(d2udr2+1rdudr)=dpdx\mu\left( \frac{d^2u}{dr^2} + \frac{1}{r}\frac{du}{dr} \right) = \frac{dp}{dx}

因為流體沿 +x+x 方向流動,壓力沿流向下降,所以

dpdx<0\frac{dp}{dx}<0

第一次積分

將方程式除以 μ\mu:

d2udr2+1rdudr=1μdpdx\frac{d^2u}{dr^2} + \frac{1}{r}\frac{du}{dr} = \frac{1}{\mu}\frac{dp}{dx}

注意左側可寫成

1rddr(rdudr)\frac{1}{r}\frac{d}{dr} \left( r\frac{du}{dr} \right)

因此

ddr(rdudr)=rμdpdx\frac{d}{dr} \left( r\frac{du}{dr} \right) = \frac{r}{\mu}\frac{dp}{dx}

積分得

rdudr=12μdpdxr2+C1r\frac{du}{dr} = \frac{1}{2\mu}\frac{dp}{dx}r^2+C_1

在管軸 r=0r=0 處,速度梯度必須有限,因此 C1=0C_1=0。於是

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 4 題

  1. The following figure shows the steady incompressible boundary layer flow over a flat plate.
    (a) (15%) Derive the momentum integral equation using the conservation of linear momentum for the flow through the control volume indicated. The result should be
    ddx∫0δρu(U−u)dy+τw=ρU2dδdx\frac{d}{dx} \int_0^\delta \rho u (U-u) dy + \tau_w = \rho U^2 \frac{d\delta}{dx}
    where δ\delta is the momentum thickness and τw\tau_w is the shear stress at wall, respectively.
    (b) (10%) Determine the boundary layer thickness and τw\tau_w, assuming a linear velocity profile across the boundary layer, i.e., u=Uyδu = U \frac{y}{\delta} for 0≤y≤δ0 \le y \le \delta and u=Uu = U for y>δy > \delta, with δ=δ(x)\delta = \delta(x) the boundary layer thickness.

🖼️【此處有附圖,請對照原卷】The figure shows a flat plate aligned with the flow. A control volume is indicated, extending from xx to x+Δxx+\Delta x and from y=0y=0 (plate surface) to y=δ(x)y=\delta(x) (boundary layer edge). The outer flow velocity is UU. Inside the boundary layer, the velocity u(x,y)u(x,y) increases from 0 at y=0y=0 to UU at y=δ(x)y=\delta(x). Streamlines are shown within the boundary layer. The shear stress at the wall is τw\tau_w. The momentum thickness is labeled as θ\theta.

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

記號釐清與核心觀念

依原卷掃描圖,第 4(a) 題的公式是

ρU2dθdx=τw,\rho U^2\frac{d\theta}{dx}=\tau_w,

其中

θ=∫0∞uU(1−uU)dy\theta=\int_0^\infty \frac{u}{U}\left(1-\frac{u}{U}\right)dy

為動量厚度;而第 4(b) 題的 δ(x)\delta(x) 才是邊界層厚度。題幹轉錄式將 θ\theta 與 δ\delta 混用;若令 δ\delta 為動量厚度,會直接導出 τw=0\tau_w=0,與壁面無滑移條件矛盾。

本題使用穩態、不可壓縮、零壓力梯度邊界層的線動量守恆。

(a) 動量積分方程式

取圖示控制體,並以單位展向長度計算。定義截面上的質量流率與 xx 方向動量流率為

m˙=∫ρu dy,Px=∫ρu2 dy.\dot m=\int \rho u\,dy, \qquad P_x=\int \rho u^2\,dy.

控制體上方位於邊界層外,流速為 UU,且無質量穿越控制面。因此穩態質量守恆給出

m˙(x+dx)−m˙(x)=0.\dot m(x+dx)-\dot m(x)=0.

由於外部流動沒有壓力梯度,控制體的 xx 方向外力只有壁面剪力。壁面對流體的作用力方向向左,因此

Px(x+dx)−Px(x)=−τw dx.P_x(x+dx)-P_x(x)=-\tau_w\,dx.

將質量守恆式乘以 UU,再由動量方程扣除,可消去外部均勻流的動量通量:

d(Px−Um˙)=−τw dx.d(P_x-U\dot m)=-\tau_w\,dx.

而

Px−Um˙=∫ρu(u−U) dy=−∫0δρu(U−u) dy.P_x-U\dot m =\int \rho u(u-U)\,dy =-\int_0^\delta \rho u(U-u)\,dy.

因此

ddx∫0δρu(U−u) dy=τw.\frac{d}{dx}\int_0^\delta \rho u(U-u)\,dy=\tau_w.

又因為

∫0δρu(U−u) dy=ρU2∫0δuU(1−uU)dy=ρU2θ,\int_0^\delta \rho u(U-u)\,dy = \rho U^2 \int_0^\delta \frac{u}{U} \left(1-\frac{u}{U}\right)dy = \rho U^2\theta,

得到原卷要求的動量積分方程式

ρU2dθdx=τw.\boxed{\rho U^2\frac{d\theta}{dx}=\tau_w}.

(b) 線性速度分布

邊界層內

u=Uyδ,0≤y≤δ.u=U\frac{y}{\delta}, \qquad 0\le y\le\delta.

動量厚度

令 η=y/δ\eta=y/\delta,則

θ=∫0δyδ(1−yδ)dy.\theta = \int_0^\delta \frac{y}{\delta} \left(1-\frac{y}{\delta}\right)dy.

計算得

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

其他考古題

114 年臺灣大學的其他科目

臺灣大學《流體力學》其他年度