112 年 國立中央大學土木工程學系碩士班戊：水資源組《流體力學》

此題考驗對水躍 (hydraulic jump) 的觀念以及其在矩形渠道中的應用。水躍是發生在開放水域中，由超臨界流 (supercritical flow) 轉變為次臨界流 (subcritical flow) 的一種流動現象，其過程中動量增加而能量減少。

(1) 證明下游處會產生水躍：
水躍產生的條件是上游的流況為超臨界流，而下游的流況為次臨界流，且下游的水位高到足以支撐水躍的發生。由題目可知，閘門上游的水深為 $y_1 = 4 \text{ m}$，閘門通過後的水深為 $y_2 = 0.5 \text{ m}$。
首先，我們計算上游的弗汝德數 (Froude number, $Fr$) 來判斷流態。對於矩形渠道，弗汝德數定義為 $Fr = \frac{V}{\sqrt{gy}}$，其中 $V$ 為平均流速，$g$ 為重力加速度，$y$ 為水深。
由於題目未給定流量或渠道寬度，我們假設渠道寬度為 $B$。設上游流量為 $Q$。則上游平均流速 $V_1 = \frac{Q}{B y_1}$。
上游的弗汝德數 $Fr_1 = \frac{V_1}{\sqrt{gy_1}} = \frac{Q/(B y_1)}{\sqrt{gy_1}} = \frac{Q}{B y_1 \sqrt{gy_1}}$。
我們知道，對於矩形渠道，當 $Fr < 1$ 時為次臨界流，當 $Fr > 1$ 時為超臨界流。
題目提到閘門上游水深為 $4 \text{ m}$，通過閘門後水深變為 $0.5 \text{ m}$。這表示閘門對水流產生了阻礙作用，使得水深減小，流速增加，這通常發生在超臨界流的條件下，例如在寬頂堰 (broad-crested weir) 或閘門 (sluice gate) 下方。
若閘門後的水深 $y_2 = 0.5 \text{ m}$ 是由超臨界流產生的，我們需要證明下游水深 $y_2$ 確實對應於超臨界流。
一個關鍵的觀念是，在特定流量下，次臨界流和超臨界流分別對應不同的水深。水躍發生時，代表流態從超臨界轉變為次臨界。
根據能量守恆和動量守恆，當流體從水深 $y_1$ 變到水深 $y_2$ 時，若 $y_2 < y_1$ 且流態從超臨界轉到次臨界，則會發生水躍。
在沒有能量損失的情況下，我們可以利用單位寬度的能量方程和動量方程來分析。
單位寬度的流量 $q = Q/B$。
單位寬度的能量 $E = y + \frac{V^2}{2g} = y + \frac{q^2}{2gy^2}$。
單位寬度的動量 $M = \frac{q^2}{gy} + \frac{y^2}{2}$。
若閘門造成水深從 $y_1 = 4 \text{ m}$ 減小到 $y_2 = 0.5 \text{ m}$，且假設在此過程中能量損失不大，則 $E_1 \approx E_2$。
$y_1 + \frac{V_1^2}{2g} \approx y_2 + \frac{V_2^2}{2g}$。
由於 $y_2 < y_1$，為了保持能量大致相等，必須 $V_2 > V_1$。
若上游是超臨界流 ($Fr_1 > 1$)，下游是次臨界流 ($Fr_2 < 1$)，且 $y_2$ 是由 $y_1$ 經過某種結構（如閘門）產生的，那麼水躍就可能發生。
水躍發生的充要條件是：從超臨界流到次臨界流的轉變，並且下游的水位足以維持此種轉變。
我們知道，對於給定的流量 $q$，存在一個臨界水深 $y_c = (\frac{q^2}{g})^{1/3}$。當 $y < y_c$ 時為超臨界流，當 $y > y_c$ 時為次臨界流。
水躍發生時，水深會從 $y_1$ 跳到 $y_2$，其中 $y_1$ 是超臨界水深， $y_2$ 是次臨界水深。
對於矩形渠道，水躍前後的水深 $y_1$ 和 $y_2$ 滿足動量方程：
$\frac{q^2}{gy_1} + \frac{y_1^2}{2} = \frac{q^2}{gy_2} + \frac{y_2^2}{2}$。
這可以改寫成 $y_2$ 與 $y_1$ 的關係。
另一種判斷水躍是否發生的方法是看 $y_2$ 是否大於 $y_1$ 的共軛水深 (conjugate depth)。
對於矩形渠道，由動量方程可知，若水深由 $y_1$ 變為 $y_2$，則：
$\frac{q^2}{g} (\frac{1}{y_1} - \frac{1}{y_2}) = \frac{y_2^2 - y_1^2}{2}$
$\frac{q^2}{g} \frac{y_2 - y_1}{y_1 y_2} = \frac{(y_2 - y_1)(y_2 + y_1)}{2}$
若 $y_1 \neq y_2$，則
$\frac{q^2}{g y_1 y_2} = \frac{y_1 + y_2}{2}$
$2 q^2 = g y_1 y_2 (y_1 + y_2)$
$q^2 = \frac{g y_1 y_2 (y_1 + y_2)}{2}$
我們也可以將動量方程改寫成 $y_2$ 與 $y_1$ 的關係式：
$y_2^3 - \frac{2q^2}{g} y_2 - \frac{2q^2}{g} y_1 - y_1^2 y_2 = 0$ (這不是標準形式)
標準的水躍共軛水深公式為：
$y_2 = \frac{y_1}{2} \left( \sqrt{1 + \frac{8q^2}{gy_1^3}} - 1 \right)$
或
$y_1 = \frac{y_2}{2} \left( \sqrt{1 + \frac{8q^2}{gy_2^3}} - 1 \right)$
這個公式是基於假設 $y_1$ 是超臨界水深，$y_2$ 是次臨界水深，且動量守恆。
題目給定 $y_1 = 4 \text{ m}$ 和 $y_2 = 0.5 \text{ m}$。
根據題目描述，閘門通過後水深變為 $0.5 \text{ m}$。這表示閘門對水流產生了很大的阻礙，使得水深急劇減小，流速急劇增加。這種情況下，上游的流況通常是超臨界流。
若我們假設 $y_1 = 4 \text{ m}$ 是超臨界水深，那麼我們可以計算出它的共軛水深 $y_2'$。如果實際的下游水深 $y_2 = 0.5 \text{ m}$ 小於 $y_2'$，則水躍會發生。
然而，題目似乎是反過來給定，閘門上游是 $4 \text{ m}$，通過閘門後是 $0.5 \text{ m}$。這表示閘門的作用是降低水深，增加流速。
若將 $y_1 = 4 \text{ m}$ 作為超臨界水深，我們需要計算 $q$。
題目中提到「請證明下游處會產生水躍」。這意味著我們需要說明為什麼會發生這種從較深水流到較淺水流的轉變，並且下游水位會導致水躍。
一個可能的情境是，閘門前的水流是次臨界流 ($y_1=4 \text{ m}$)，通過閘門後，由於閘門的結構，水深急劇減小到 $y_2=0.5 \text{ m}$，使得流速急劇增加，變為超臨界流。然後，在下游某處，由於水流的慣性或下游水位條件，流態又從超臨界變回次臨界，發生水躍。
但題目描述「閘門上游的水深為 $4 \text{ m}$，而當閘門些微抬起時。水通過閘門後產生一較淺水深為 $50 \text{ cm}$」，這暗示著 $4 \text{ m}$ 是閘門前的「正常」水深，而 $0.5 \text{ m}$ 是閘門後的「正常」水深。
通常，閘門（如 weirs 或 sluice gates）的作用是產生超臨界流。所以，我們假設 $4 \text{ m}$ 是次臨界水深，而 $0.5 \text{ m}$ 是超臨界水深。
如果 $y_1 = 4 \text{ m}$ 是次臨界水深 ($Fr_1 < 1$)，而 $y_2 = 0.5 \text{ m}$ 是超臨界水深 ($Fr_2 > 1$)，那麼這是不符合水躍從超臨界到次臨界的定義。
水躍的定義是從超臨界流轉變為次臨界流。因此，我們應該將 $y_1 = 4 \text{ m}$ 視為閘門前的「正常」水深，而 $y_2 = 0.5 \text{ m}$ 是閘門處產生的水深。
題目說「當閘門些微抬起時」，這表示閘門的開度是可調的。
標準的水躍現象是指，當超臨界流遇到阻礙或下游水位升高時，會發生水躍，將流態轉變為次臨界流。
在此題中，閘門的作用是減小水深並增加流速。所以，我們假設閘門前的流況是次臨界流，水深 $y_1 = 4 \text{ m}$。通過閘門後，流速急劇增加，水深減小到 $y_2 = 0.5 \text{ m}$，這變成了超臨界流。
如果下游的水位（由 $y_2$ 和某個下游水位共同決定）能夠支撐超臨界流發展，並在某處發生回歸為次臨界流，則會形成水躍。
題目問「證明下游處會產生水躍」。這表示我們需要說明為什麼會發生這種情況。
通常，水躍是發生在超臨界流進入次臨界流的時候。
設閘門前的渠道寬度為 $B$。設流量為 $Q$。
上游水深 $y_1 = 4 \text{ m}$。
通過閘門後的水深 $y_2 = 0.5 \text{ m}$。
若 $y_1 = 4 \text{ m}$ 是次臨界水深，則 $Fr_1 < 1$。
若 $y_2 = 0.5 \text{ m}$ 是超臨界水深，則 $Fr_2 > 1$。
這種從次臨界流到超臨界流的轉變，通常是由於流道變窄或有結構導致加速。閘門的作用就是如此。
水躍是從超臨界流轉變為次臨界流。所以，應該是某個超臨界水深 $y_a$ 經過水躍變為次臨界水深 $y_b$。
題目給定的 $4 \text{ m}$ 和 $0.5 \text{ m}$ 是閘門前後的水深。
最合理的解釋是：閘門前的水流是次臨界流 ($y_1 = 4 \text{ m}$)，通過閘門後，水流被壓縮，水深減小到 $y_2 = 0.5 \text{ m}$，流速急劇增加，變成了超臨界流。然後，在下游某處，由於能量損失和下游水位的作用，流態從超臨界轉變回次臨界，這就是水躍。
題目問「證明下游處會產生水躍」。這句話的重點在於「下游處」。
水躍的發生需要滿足動量守恆和能量損失。
對於給定的流量 $q$，存在一個共軛水深關係：
$y_2 = \frac{y_1}{2} \left( \sqrt{1 + \frac{8q^2}{gy_1^3}} - 1 \right)$ （若 $y_1$ 是超臨界水深）
$y_1 = \frac{y_2}{2} \left( \sqrt{1 + \frac{8q^2}{gy_2^3}} - 1 \right)$ （若 $y_2$ 是超臨界水深）
或者，若 $y_1$ 和 $y_2$ 是共軛水深，則 $y_1+y_2 = \frac{2q^2}{gy_1y_2}$。
我們需要確定 $q$ 的值。題目沒有直接給出。
但是，題目說「證明下游處會產生水躍」。這暗示著，閘門後的 $0.5 \text{ m}$ 水深，雖然是超臨界流，但下游的水位（或是某種阻礙）導致了水躍的發生。
假設閘門後的 $0.5 \text{ m}$ 水深是超臨界流，我們稱之為 $y_a = 0.5 \text{ m}$。
為了證明水躍發生，我們需要證明存在一個次臨界水深 $y_b$ 與 $y_a$ 是共軛水深，並且下游的水位至少能達到 $y_b$。
共軛水深公式：$y_b = \frac{y_a}{2} (\sqrt{1 + \frac{8q^2}{gy_a^3}} - 1)$。
我們不知道 $q$。但是，我們知道閘門前的 $y_1 = 4 \text{ m}$。
如果 $y_1 = 4 \text{ m}$ 是次臨界流，那麼 $Fr_1 < 1$。
通過閘門後，水深變為 $y_2 = 0.5 \text{ m}$。這通常意味著流速增加，變為超臨界流，即 $Fr_2 > 1$。
所以，閘門的作用是將次臨界流轉變為超臨界流。
然後，這個超臨界流 ($y_2 = 0.5 \text{ m}$) 在下游會發生水躍，變為次臨界流。
要證明水躍發生，我們需要證明：

1. 閘門後的流況是超臨界流。
2. 下游的水位足夠高，以至於超臨界流無法維持，必須轉變為次臨界流。
   我們仍然不知道 $q$。但題目給了 $y_1=4 \text{ m}$ 和 $y_2=0.5 \text{ m}$。
   如果我們假設這兩個水深是與某個流量 $q$ 相關的，例如，通過閘門前後，流量是守恆的。
   若 $y_1 = 4 \text{ m}$ 是次臨界水深，則 $Fr_1 < 1$。
   若 $y_2 = 0.5 \text{ m}$ 是超臨界水深，則 $Fr_2 > 1$。
   水躍是從超臨界流到次臨界流的轉變。
   因此，我們需要證明，在 $y_2 = 0.5 \text{ m}$ 的超臨界流之後，會發生一個水躍，形成一個較大的次臨界水深。
   如何「證明」？
   一個方法是利用能量和動量守恆。
   在水躍之前，流體是超臨界流，水深為 $y_1'$。在水躍之後，流體是次臨界流，水深為 $y_2'$。
   題目給定的 $y_1 = 4 \text{ m}$ 和 $y_2 = 0.5 \text{ m}$ 是閘門前後的水深。
   最關鍵的點是：水躍是從超臨界流轉變到次臨界流。

核心觀念

本題依序使用：

1. 連續方程式求 $30\text{ cm}$ 與 $15\text{ cm}$ 管內平均流速。
2. 擴充伯努力方程式，計算沿程損失與局部損失。
3. 線性動量方程式，求水流對管線及支撐 $C$ 的水平作用力。

取 $B$ 為向大氣排出的出口，因此 $p_B=0$（表壓）；依圖示，$A$ 比 $B$ 高 $0.5\text{ m}$。

---

解題方法與基本流速

大管徑與小管徑分別為

$$D_1=0.30\text{ m},\qquad D_2=0.15\text{ m}$$

截面積：

$$A_1=\frac{\pi D_1^2}{4}=0.07069\text{ m}^2$$ $$A_2=\frac{\pi D_2^2}{4}=0.01767\text{ m}^2$$

由連續方程式 $Q=AV$：

$$V_1=\frac{0.0707}{0.07069}\approx 1.00\text{ m/s}$$ $$V_2=\frac{0.0707}{0.01767}\approx 4.00\text{ m/s}$$

因此

$$\frac{V_1^2}{2g}=\frac{1^2}{2(9.81)} =0.05099\text{ m}$$ $$\frac{V_2^2}{2g}=\frac{4^2}{2(9.81)} =0.8155\text{ m}$$

圖中 $30\text{ cm}$ 直管總長度為

$$L=1+0.5+2=3.5\text{ m}$$

且共有兩個 $90^\circ$ 彎管。

---

(1) $A$ 至 $B$ 的總水頭損失

沿程損失：

$$h_f=f\frac{L}{D_1}\frac{V_1^2}{2g}$$ $$h_f =0.018\left(\frac{3.5}{0.30}\right)(0.05099) =0.0107\text{ m}$$

兩個 $90^\circ$ 彎管損失：

$$h_{\text{bend}} =2(0.9)\frac{V_1^2}{2g} =2(0.9)(0.05099) =0.0918\text{ m}$$

閥門損失：

$$h_{\text{valve}} =0.2\frac{V_1^2}{2g} =0.0102\text{ m}$$

束縮損失係數以束縮後的小管徑流速為基準：

$$h_{\text{contraction}} =0.1\frac{V_2^2}{2g} =0.1(0.8155) =0.0816\text{ m}$$

所以總水頭損失為

$$h_{L,AB} =h_f+h_{\text{bend}}+h_{\text{valve}}+h_{\text{contraction}}$$ $$h_{L,AB} =0.0107+0.0918+0.0102+0.0816 \approx 0.194\text{ m}$$

---

(2) 管線入口 $A$ 處所需壓力

核心觀念

本題屬於流體靜力學中經典的**浮體平衡與浮性穩定度（Stability of Floating Bodies）**分析，主要涵蓋兩大核心定理與觀念：

1. 阿基米德原理（浮力平衡）：
   浮體在靜止液體中呈垂直力平衡時，浮體總重力 $W$ 等於其所排開流體的總重量（即浮力 $F_B$）：
   $W = F_B = \rho_w g V_{\text{sub}}$
   其中 $\rho_w$ 為液體密度，$V_{\text{sub}}$ 為浸水體積（排開水體積）。

2. 定傾高（Metacentric Height, $\overline{GM}$）穩定性判別法：
   浮體受小角度傾側微擾時，其穩定性取決於定傾中心 $M$ 與重心 $G$ 的相對位置：

   • $\overline{GM} = \overline{KB} + \overline{BM} - \overline{KG}$
   • $\overline{KB}$：浮心 $B$（排開水體之幾何形心）距底面的垂直距離。
   • $\overline{KG}$：整體木塊重心 $G$ 距底面的垂直距離。
   • $\overline{BM} = \dfrac{I}{V_{\text{sub}}}$：定傾半徑（Metacentric Radius），$I$ 為浮體在水面處截面（水線面）繞傾側軸之面積二次矩（慣性矩）。
   • 穩定平衡條件：$\overline{GM} > 0$（定傾中心 $M$ 必須在重心 $G$ 之上方）。

---

解題方法

以木塊底部為座標原點基準線（$y = 0$），垂直向上為正。

(1) 求木塊之吃水深度 $d$

1. 計算木塊總質量與總重力：
   木塊底面積 $A = B \times B = B^2$，總高度為 $H$。

   • 下半部（高度 $0 \le y \le \dfrac{H}{2}$）：密度為 $\rho_s$，體積 $V_1 = B^2 \left(\dfrac{H}{2}\right)$。
     質量 $M_1 = \rho_s \cdot \dfrac{B^2 H}{2}$。
   • 上半部（高度 $\dfrac{H}{2} \le y \le H$）：密度為 $\dfrac{\rho_s}{2}$，體積 $V_2 = B^2 \left(\dfrac{H}{2}\right)$。
     質量 $M_2 = \left(\dfrac{\rho_s}{2}\right) \cdot \dfrac{B^2 H}{2} = \dfrac{\rho_s B^2 H}{4}$。
   • 木塊總重力 $W$：
     $W = (M_1 + M_2)g = \left(\frac{1}{2} + \frac{1}{4}\right)\rho_s B^2 H g = \frac{3}{4} \rho_s g B^2 H$

2. 建立浮力平衡方程式：
   設吃水深度為 $d$。不論吃水深度是否超過中線 $\dfrac{H}{2}$，排開的流體均為均勻的水（密度 $\rho_w$），排開水體積為均勻柱體 $V_{\text{sub}} = B^2 d$。
   浮力為：
   $F_B = \rho_w g V_{\text{sub}} = \rho_w g B^2 d$
   由垂直靜力平衡 $W = F_B$：
   $\frac{3}{4} \rho_s g B^2 H = \rho_w g B^2 d$
   解得吃水深度 $d$：
   $d = \frac{3}{4}\frac{\rho_s}{\rho_w} H$

---

(2) 求木塊穩定不翻覆之幾何條件 $(B/H)$

1. 求重心位置 $\overline{KG}$：
   以木塊底部 $y = 0$ 為基準，利用力矩原理計算複合材料的質心位置：

   • 下半物質心位置：$y_1 = \dfrac{H}{4}$
   • 上半物質心位置：$y_2 = \dfrac{H}{2} + \dfrac{H}{4} = \dfrac{3H}{4}$
     $\overline{KG} = \frac{M_1 y_1 + M_2 y_2}{M_1 + M_2} = \frac{\left(\dfrac{\rho_s B^2 H}{2}\right)\left(\dfrac{H}{4}\right) + \left(\dfrac{\rho_s B^2 H}{4}\right)\left(\dfrac{3H}{4}\right)}{\dfrac{3}{4}\rho_s B^2 H}$
     分子分母同除以 $\rho_s B^2 H$：
     $\overline{KG} = \frac{\dfrac{1}{8}H + \dfrac{3}{16}H}{\dfrac{3}{4}} = \frac{\dfrac{5}{16}H}{\dfrac{3}{4}} = \frac{5}{12} H$

2. 求浮心位置 $\overline{KB}$：
   浮心為浸水體積（長、寬皆為 $B$、深度為 $d$ 的均勻水體柱）之形心，其距離底面的高度為：

此題考驗旋轉座標系下的流體力學，特別是圓柱座標系下的流動分析。

題目給定：
圓桶直徑 $D$，半徑 $R = D/2$。
圓桶深 $H$。
初始水深 $h_0 = H/2$。
圓筒以角速度 $\omega$ 轉動。
流速 $V(r) = \omega r$ (沿圓周方向)。
流線為同心圓。

(1) 解釋此流況是否為均勻流 (uniform flow) 或旋轉流 (rotational flow)?
均勻流 (Uniform flow) 指的是流場中速度向量在空間上不隨位置變化。即 $\frac{\partial \vec{v}}{\partial x} = \frac{\partial \vec{v}}{\partial y} = \frac{\partial \vec{v}}{\partial z} = 0$ (在笛卡爾座標系下)。
旋轉流 (Rotational flow) 指的是流體質點具有角速度，即渦量 (vorticity) 不為零。在圓柱座標系 $(r, \theta, z)$ 下，渦量 $\vec{\omega}_{rot} = \nabla \times \vec{v}$。對於軸對稱流動，主要的渦量分量是 $\omega_{z,rot} = \frac{1}{r} \frac{\partial}{\partial r} (r V_\theta) - \frac{1}{r} \frac{\partial V_r}{\partial \theta}$。
在本題中，速度向量為 $\vec{v} = V_\theta(r) \hat{e}_\theta = \omega r \hat{e}_\theta$。
速度 $V_\theta$ 只隨 $r$ 變化，不隨 $\theta$ 或 $z$ 變化。
然而，速度向量本身在空間上是變化的。例如，在不同的 $r$ 值處，速度大小不同。在不同的 $\theta$ 值處，速度方向不同。所以，這不是均勻流。

判斷是否為旋轉流：
渦量 $\vec{\omega}_{rot} = \nabla \times \vec{v}$。
在圓柱座標系下，速度 $\vec{v} = V_r \hat{e}_r + V_\theta \hat{e}_\theta + V_z \hat{e}_z$。
本題中，$V_r = 0$, $V_\theta = \omega r$, $V_z = 0$。
渦量的 $\theta$ 分量為 $\omega_{\theta,rot} = \frac{1}{r} \frac{\partial V_z}{\partial \theta} - \frac{\partial V_\theta}{\partial z} = 0 - 0 = 0$。
渦量的 $z$ 分量為 $\omega_{z,rot} = \frac{1}{r} \frac{\partial}{\partial r} (r V_\theta) - \frac{1}{r} \frac{\partial V_r}{\partial \theta} = \frac{1}{r} \frac{\partial}{\partial r} (r (\omega r)) - 0 = \frac{1}{r} \frac{\partial}{\partial r} (\omega r^2) = \frac{1}{r} (2 \omega r) = 2 \omega$。
由於 $\omega_{z,rot} = 2 \omega \neq 0$ (假設 $\omega \neq 0$)，所以此流況是旋轉流。
這個流動是一個自由渦 (free vortex) 的模型，其中速度隨半徑增大而增大，並且具有非零的渦量。

(2) 推求水深 $h(r)$ 沿徑向 $(r$ 方向) 之變化。
當圓筒旋轉時，由於離心力的作用，水面會形成一個拋物面。
我們可以利用歐拉方程或考慮質點所受的力來推導。
考慮一個質量為 $dm$ 的水微元，在半徑為 $r$ 處。
它受到重力 $dm \cdot g$ (向下)，以及壓力梯度力 $dp \cdot dA$ (徑向)。
在旋轉參考系中，還需要考慮科氏力 (Coriolis force) 和離心慣性力 (centrifugal inertial force)。
然而，題目說「水流速度 $V(r) = \omega r$」，這已經給定了速度場。
我們需要考慮的是，在這種速度場下，水面會形成什麼形狀。
水面是一個自由表面，其上的壓力是恆定的 (通常為大氣壓 $P_{atm}$)。
我們可以使用伯努利方程來分析。
在旋轉參考系中，考慮一個在 $(r, z)$ 平面上的水微元。
作用在水微元上的力是：壓力梯度力、重力、離心力。
考慮徑向方向的力平衡。
對於一個在半徑 $r$ 處，高度 $z$ 的水微元，其受力分析：
徑向壓力梯度力 $dp/dr$。
離心力 $\rho \omega^2 r$。
在穩態情況下，合力為零（或者說，壓力梯度平衡了離心力）。
$\frac{dp}{dr} = \rho \omega^2 r$。
積分得到壓力 $p(r)$：
$p(r) = \int \rho \omega^2 r dr = \frac{1}{2} \rho \omega^2 r^2 + C_1$。
其中 $C_1$ 是積分常數。

現在考慮垂直方向 (z 方向)。
在垂直方向上，壓力梯度平衡重力。
$\frac{dp}{dz} = -\rho g$。
積分得到壓力 $p(z)$：
$p(z) = -\rho g z + C_2$。
其中 $C_2$ 是積分常數。

水面是一個等壓面。
所以，壓力 $p$ 在水面上是恆定的。
我們可以將 $p(r, z)$ 寫成 $p(r, z) = \frac{1}{2} \rho \omega^2 r^2 - \rho g z + C'$ (將兩個積分常數合併)。
水面是一個等壓面，例如 $p = P_{atm}$。
$\frac{1}{2} \rho \omega^2 r^2 - \rho g z + C' = P_{atm}$。
$z = \frac{\frac{1}{2} \rho \omega^2 r^2 + C' - P_{atm}}{\rho g}$。
$z = \frac{\omega^2}{2g} r^2 + \text{constant}$。
這表示水面是拋物面。
設水面方程為 $z = f(r)$。
我們需要確定這個常數。
初始狀態：靜止時，水深 $h_0 = H/2$。水面是平的， $z = H/2$ (相對於桶底)。
當圓筒旋轉時，水的總體積保持不變 (如果沒有溢出)。
設水面的高度為 $z(r)$。
桶底為 $z=0$。桶壁為 $r=R=D/2$。
桶的總高度為 $H$。
水面方程為 $z(r) = z_0 + \frac{\omega^2 r^2}{2g}$，其中 $z_0$ 是圓桶中心處的水面高度。
水的總體積 $V = \int_0^R 2\pi r z(r) dr$ (假設桶壁是垂直的)。
$V = \int_0^R 2\pi r (z_0 + \frac{\omega^2 r^2}{2g}) dr$
$V = 2\pi [z_0 \frac{r^2}{2} + \frac{\omega^2}{2g} \frac{r^4}{4}]_0^R$
$V = 2\pi [z_0 \frac{R^2}{2} + \frac{\omega^2 R^4}{8g}]$。
初始體積 $V_0 = \pi R^2 h_0 = \pi R^2 (H/2)$。