112 年 國立中央大學地球科學學系地球物理碩士班在職生《流體力學》

此題考驗對溫度梯度（temperature gradient）概念的理解，以及如何結合相對運動與絕對變化來判斷梯度。

題目給了幾個資訊：

1. 觀察者（車）的運動：沿南北向公路往北走，車速 $v_{car} = 100 \text{ km/hr}$。
2. 車上溫度計記錄的變化率：$\frac{dT_{car}}{dt} = 5 \text{ °C/hr}$。這是車內觀察者感受到的溫度變化率，包含了環境溫度變化和車輛移動造成的溫度變化。
3. 氣象報告的環境溫度變化率：各地氣溫隨時間下降的速率為 $2 \text{ °C/hr}$。這代表環境溫度隨時間的變化，即 $\frac{\partial T}{\partial t} = -2 \text{ °C/hr}$。
4. 目標：求南北方向的溫度梯度大小與方向。溫度梯度是空間上的變化率，即 $\nabla T$。

我們知道，觀察者感受到的溫度變化率 $\frac{dT_{car}}{dt}$ 可以用鏈式法則表示為：
$\frac{dT_{car}}{dt} = \frac{\partial T}{\partial t} + \mathbf{v}_{car} \cdot \nabla T$
其中：

• $\frac{\partial T}{\partial t}$ 是環境溫度隨時間的變化率。
• $\mathbf{v}_{car}$ 是觀察者（車輛）的速度向量。
• $\nabla T$ 是溫度梯度向量。

我們已知：

• $\frac{dT_{car}}{dt} = 5 \text{ °C/hr}$
• $\frac{\partial T}{\partial t} = -2 \text{ °C/hr}$ (氣象報告說氣溫下降，所以是負值)
• 車輛沿南北向往北走，所以 $\mathbf{v}_{car}$ 的方向是往北。若我們定義北方為正向，則 $\mathbf{v}_{car} = (100, 0) \text{ km/hr}$ (假設在 x-y 平面，x 代表東西向，y 代表南北向，或者直接將南北向視為一個方向)。
• 我們要找的溫度梯度 $\nabla T$ 是南北向的。假設溫度梯度只有南北向的分量，即 $\nabla T = (0, \frac{\partial T}{\partial y})$，其中 $\frac{\partial T}{\partial y}$ 是溫度在南北方向的變化率。

將這些代入鏈式法則公式：
$5 \text{ °C/hr} = -2 \text{ °C/hr} + \mathbf{v}_{car} \cdot \nabla T$
$7 \text{ °C/hr} = \mathbf{v}_{car} \cdot \nabla T$

由於車輛往北，我們假設 $\mathbf{v}_{car}$ 的方向是往北。同時，我們也假設溫度梯度 $\nabla T$ 的主要分量也是南北向的。令 $\nabla T$ 在南北向的分量為 $G_y = \frac{\partial T}{\partial y}$。
如果車輛往北 $\mathbf{v}_{car}$，而我們關心的溫度梯度也是南北向 $\nabla T$，那麼 $\mathbf{v}_{car} \cdot \nabla T = |\mathbf{v}_{car}| |\nabla T| \cos \theta$。
題目問的是「南北的溫度梯度」。這暗示我們只需要考慮溫度在南北方向的變化。

假設 $\mathbf{v}_{car}$ 是往北，速度大小為 $v_{car} = 100 \text{ km/hr}$。
假設溫度梯度 $\nabla T$ 在南北方向的分量為 $G_{north}$。

此題考驗對流體運動基本概念的理解，包括恆穩性（steady）、維度（dimension）、不可壓縮性（incompressibility），以及全加速度（total acceleration）的計算。

首先，我們分析速度場 $\mathbf{v} = (u,v) = \langle u(x,y,t), v(x,y,t), w(x,y,t) \rangle$。
題目給的速度向量為 $\mathbf{v} = (u,v) = \langle x^2t^2 + (-2xy-1), y^2t^2 + (2xy-1) \rangle$。
這裡只給定了 $u$ 和 $v$ 分量，而 $w$ 分量未給出。通常，如果速度向量只包含 $x$ 和 $y$ 分量，且這兩個分量不依賴於 $z$，我們就稱之為二維流動。

1. 恆穩性 (Steady Flow)？
一個流動是恆穩的，如果流場的任何物理量（如速度、壓力等）在空間的固定點上隨時間不變。也就是說，速度向量 $\mathbf{v}(x,y,t)$ 的時間偏導數為零：$\frac{\partial \mathbf{v}}{\partial t} = \mathbf{0}$。
我們來看速度場的兩個分量：
$u(x,y,t) = x^2t^2 - 2xy - 1$
$v(x,y,t) = y^2t^2 + 2xy - 1$

計算時間偏導數：
$\frac{\partial u}{\partial t} = \frac{\partial}{\partial t}(x^2t^2 - 2xy - 1) = 2x^2t$
$\frac{\partial v}{\partial t} = \frac{\partial}{\partial t}(y^2t^2 + 2xy - 1) = 2y^2t$

因為 $\frac{\partial u}{\partial t}$ 和 $\frac{\partial v}{\partial t}$ 一般不為零（除非 $x=0$ 或 $y=0$ 或 $t=0$），所以速度場隨時間變化。因此，這個流動 不是恆穩的 (unsteady)。

2. 二維流動還是三維流動？
題目給的速度向量只有 $u$ 和 $v$ 兩個分量，且這兩個分量都不依賴於 $z$ 座標。通常，這意味著 $w=0$ 且 $\frac{\partial}{\partial z}$ 對所有物理量都為零。因此，我們可以將此流動視為 二維流動。
如果題目意指速度向量 $\mathbf{v} = \langle u(x,y,t), v(x,y,t), w(x,y,t) \rangle$ 且 $w$ 函數未給定，但從表達式看，它只依賴 $x,y,t$，那它就不是三維流動。若 $w$ 真的存在且依賴 $z$，則需要額外資訊。基於題目給的表達式，我們推斷為二維。

3. 不可壓縮流 (Incompressible Flow)？
一個流動是不可壓縮的，如果其質量守恆方程式（連續方程式）滿足：
$\nabla \cdot \mathbf{v} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$
我們已經假設這是二維流動，所以 $w=0$ 且 $\frac{\partial w}{\partial z}=0$。我們只需要計算 $\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}$。

計算空間偏導數：
$\frac{\partial u}{\partial x} = \frac{\partial}{\partial x}(x^2t^2 - 2xy - 1) = 2xt^2 - 2y$
$\frac{\partial v}{\partial y} = \frac{\partial}{\partial y}(y^2t^2 + 2xy - 1) = 2yt^2 + 2x$

計算散度：
$\nabla \cdot \mathbf{v} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = (2xt^2 - 2y) + (2yt^2 + 2x) = 2xt^2 - 2y + 2yt^2 + 2x$
這個散度一般不為零（除非 $x=0, y=0, t=0$ 或者特定的組合）。因此，這個流動 不是不可壓縮流。

4. 計算全加速度 (Total Acceleration)
全加速度 $\mathbf{a}$ 的表達式為：
$\mathbf{a} = \frac{D\mathbf{v}}{Dt} = \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v}$
其中 $(\mathbf{v} \cdot \nabla)\mathbf{v}$ 是對流加速度（convective acceleration）。

我們已經計算了時間導數部分：
$\frac{\partial u}{\partial t} = 2x^2t$
$\frac{\partial v}{\partial t} = 2y^2t$
所以，時間加速度向量為 $\mathbf{a}_t = \langle 2x^2t, 2y^2t \rangle$。

現在計算對流加速度部分 $(\mathbf{v} \cdot \nabla)\mathbf{v}$。
$\mathbf{v} \cdot \nabla = u \frac{\partial}{\partial x} + v \frac{\partial}{\partial y}$
$(\mathbf{v} \cdot \nabla)\mathbf{v} = \left( u \frac{\partial}{\partial x} + v \frac{\partial}{\partial y} \right) \langle u, v \rangle = \left\langle u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y}, u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} \right\rangle$

我們需要計算這些偏導數：
$\frac{\partial u}{\partial x} = 2xt^2 - 2y$ (已計算)
$\frac{\partial v}{\partial y} = 2yt^2 + 2x$ (已計算)

$\frac{\partial u}{\partial y} = \frac{\partial}{\partial y}(x^2t^2 - 2xy - 1) = -2x$
$\frac{\partial v}{\partial x} = \frac{\partial}{\partial x}(y^2t^2 + 2xy - 1) = 2y$

現在代入對流加速度的計算：
x-分量：
$u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = (x^2t^2 - 2xy - 1)(2xt^2 - 2y) + (y^2t^2 + 2xy - 1)(-2x)$
$= (2x^3t^4 - 2x y t^2 - 2x t^2) - (2x^2y t^2 - 4xy^2 - 2y) - 2xy^2t^2 - 4x^2y + 2x$
$= 2x^3t^4 - 2x y t^2 - 2x t^2 - 2x^2y t^2 + 4xy^2 + 2y - 2xy^2t^2 - 4x^2y + 2x$

y-分量：
$u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} = (x^2t^2 - 2xy - 1)(2y) + (y^2t^2 + 2xy - 1)(2yt^2 + 2x)$
$= 2x^2yt^2 - 4xy^2 - 2y + (2y^3t^4 + 2xy^2t^2) + (4xy^3t^2 + 4x^2y) - 2yt^2 - 2x$
$= 2x^2yt^2 - 4xy^2 - 2y + 2y^3t^4 + 2xy^2t^2 + 4xy^3t^2 + 4x^2y - 2yt^2 - 2x$

我們需要計算在 $(x,y) = (4, -2)$ 的全加速度。
為了簡化計算，我們將 $x=4$ 和 $y=-2$ 代入。
首先，代入 $u$ 和 $v$ 的值：
$u(4, -2, t) = (4)^2t^2 - 2(4)(-2) - 1 = 16t^2 + 16 - 1 = 16t^2 + 15$
$v(4, -2, t) = (-2)^2t^2 + 2(4)(-2) - 1 = 4t^2 - 16 - 1 = 4t^2 - 17$

代入時間加速度：
$\frac{\partial u}{\partial t} = 2x^2t = 2(4^2)t = 32t$
$\frac{\partial v}{\partial t} = 2y^2t = 2((-2)^2)t = 8t$
$\mathbf{a}_t(4, -2, t) = \langle 32t, 8t \rangle$

現在代入對流加速度的偏導數，並在 $(x,y)=(4,-2)$ 計算：
$\frac{\partial u}{\partial x} = 2xt^2 - 2y = 2(4)t^2 - 2(-2) = 8t^2 + 4$
$\frac{\partial v}{\partial y} = 2yt^2 + 2x = 2(-2)t^2 + 2(4) = -4t^2 + 8$
$\frac{\partial u}{\partial y} = -2x = -2(4) = -8$
$\frac{\partial v}{\partial x} = 2y = 2(-2) = -4$

計算對流加速度在 $(4, -2)$ 的 x-分量：
$u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = (16t^2 + 15)(8t^2 + 4) + (4t^2 - 17)(-8)$
$= (128t^4 + 64t^2 + 120t^2 + 60) - 32t^2 + 136$
$= 128t^4 + 184t^2 + 60 + 136$
$= 128t^4 + 184t^2 + 196$

計算對流加速度在 $(4, -2)$ 的 y-分量：

此題考驗對質量守恆定律（Continuity Equation）不同形式的理解及其物理意義。質量守恆定律是流體力學中最基本且重要的定律之一，它指出在一個封閉系統中，質量既不能被創生也不能被消滅。

題目給了三個方程式，分別是積分形式、微分形式（張量形式）以及針對質點的物質導數形式。

第一個方程式：積分形式 (Integral Form)
$\frac{\partial}{\partial t} \int_{V} \rho dV = - \oint_{S} \rho \mathbf{v} \cdot d\mathbf{S}$

• 左邊: $\frac{\partial}{\partial t} \int_{V} \rho dV$
  • $\int_{V} \rho dV$ 代表一個固定體積 $V$ 內所含有的總質量。
  • $\rho$ 是流體的密度（mass per unit volume）。
  • $\frac{\partial}{\partial t} \int_{V} \rho dV$ 代表這個固定體積 $V$ 內總質量隨時間的變化率。
• 右邊: $- \oint_{S} \rho \mathbf{v} \cdot d\mathbf{S}$
  • $S$ 是包圍體積 $V$ 的閉合表面。
  • $d\mathbf{S}$ 是表面 $S$ 上的一個微小面積向量，其方向垂直於表面並指向外。
  • $\rho \mathbf{v}$ 是質量通量向量（mass flux vector），表示單位時間單位面積通過的質量。
  • $\rho \mathbf{v} \cdot d\mathbf{S}$ 是通過微小面積 $d\mathbf{S}$ 的淨質量流量。如果 $\mathbf{v}$ 和 $d\mathbf{S}$ 同向（流體流出），則為正；如果反向（流體流入），則為負。
  • $\oint_{S} \rho \mathbf{v} \cdot d\mathbf{S}$ 是通過整個閉合表面 $S$ 的淨質量流量（net mass flow rate out of the control volume）。
  • 負號 $-\oint_{S} \rho \mathbf{v} \cdot d\mathbf{S}$ 代表流入體積 $V$ 的淨質量流量。
• 整體意義: 這個方程式說明，一個固定控制體積 $V$ 內總質量的變化率，等於通過該體積表面 $S$ 流入的淨質量率。簡單來說，就是「進入多少質量，質量就增加多少；出去多少質量，質量就減少多少」。這是對質量守恆定律的一種「控制體積」表述。

第二個方程式：微分形式 (Differential Form)
$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$
這個方程式是從第一個積分形式，利用高斯散度定理（Gauss's Divergence Theorem）轉換而來的，適用於流場中的每一個點。

• 左邊第一項: $\frac{\partial \rho}{\partial t}$
  • $\rho$ 是密度。
  • $\frac{\partial \rho}{\partial t}$ 代表在空間中某固定點上的密度隨時間的變化率。這是由於質量在該點累積或耗散所致。
• 左邊第二項: $\nabla \cdot (\rho \mathbf{v})$
  • $\rho \mathbf{v}$ 是質量通量向量（mass flux vector）。
  • $\nabla \cdot (\rho \mathbf{v})$ 是質量通量的散度（divergence of mass flux）。它表示在空間中某固定點上，從該點發散出去的淨質量通量（單位體積）。如果散度為正，表示有淨質量流出；如果為負，表示有淨質量流入。
• 整體意義: 這個方程式說明，在空間的任何一點，密度的時間變化率（$\frac{\partial \rho}{\partial t}$）加上該點的質量通量的散度（$\nabla \cdot (\rho \mathbf{v})$），其總和為零。
  • 若 $\frac{\partial \rho}{\partial t} > 0$ (密度增加)，則必須有 $\nabla \cdot (\rho \mathbf{v}) < 0$ (淨質量流入)，以保持總和為零。
  • 若 $\frac{\partial \rho}{\partial t} < 0$ (密度減少)，則必須有 $\nabla \cdot (\rho \mathbf{v}) > 0$ (淨質量流出)。
  • 這就是質量守恆定律的微分形式，它描述了密度在空間中如何隨著質量通量的分佈而變化。

第三個方程式：物質導數形式 (Material Derivative Form)
$\frac{d\rho}{dt} = - \rho \nabla \cdot \mathbf{v}$
這個方程式是從第二個微分形式推導出來的，它描述了一個隨流體運動的質點（material particle）的密度如何變化。

• 左邊: $\frac{d\rho}{dt}$
  • $\frac{d}{dt}$ 是物質導數（material derivative 或 total derivative），表示沿著流體質點的運動軌跡的變化率。
  • $\frac{d\rho}{dt}$ 代表一個隨流體運動的質點的密度隨時間的變化率。
• 右邊: $- \rho \nabla \cdot \mathbf{v}$
  • $\rho$ 是流體密度。
  • $\nabla \cdot \mathbf{v}$ 是速度場的散度（divergence of velocity）。它表示單位體積的流體在某點的膨脹率。

此題考驗對環流量（Circulation）及其變化率的理解，以及正壓（barotropic）與斜壓（baroclinic）流體的區別，並應用至海陸風和渦度-環流量關係的分析。

首先，題目中給出的環流量變化方程式應為動量守恆（Navier-Stokes equation）在環流量上的應用。標準的環流量定理（Kelvin's Circulation Theorem 或 Reynolds Transport Theorem for Circulation）是：
$\frac{d\Gamma}{dt} = \int_{S} \frac{\partial \mathbf{v}}{\partial t} \cdot d\mathbf{A} + \oint_{C} \mathbf{v} \cdot d\mathbf{r} \text{ (此處有誤，應為 } \oint_{C} \frac{Dp}{Dt} \cdot d\mathbf{r} \text{ 或相關項)}$
題目給的公式形式較為複雜，可能是一個特定的形式，或是將動量方程的時間項、對流項、壓力梯度項、重力項合併。
我們假設題目給的公式是有效的，它描述了環流量 $\Gamma$ 的變化率 $\frac{d\Gamma}{dt}$。
環流量 $\Gamma$ 定義為沿一個閉合曲線 $C$ 的速度積分：$\Gamma = \oint_C \mathbf{v} \cdot d\mathbf{r}$。
題目給的公式中的積分涉及到面積 $S$，這暗示了 $\Gamma$ 的變化率與通過某個曲面 $S$ 的相關物理量有關。

標準的環流量定理（Kelvin's Theorem）
$\frac{d\Gamma}{dt} = \frac{\partial \Gamma}{\partial t} + \oint_{C} \frac{Dp}{Dt} \cdot d\mathbf{r} + \oint_{C} \nu \nabla^2 \mathbf{v} \cdot d\mathbf{r} + \oint_{C} \nabla (\Phi) \cdot d\mathbf{r}$
其中 $\nu$ 是運動黏滯係數，$\Phi$ 是力勢。
對於無黏滯、無外力、壓力與密度單值函數的流體，$\frac{d\Gamma}{dt} = 0$，即環流量守恆。

題目給的公式 $\frac{\partial}{\partial t} \int_{V} \rho dV = - \oint_{S} \rho \mathbf{v} \cdot d\mathbf{S}$ 是質量守恆。
題目給的公式 $\frac{d\rho}{dt} = - \rho \nabla \cdot \mathbf{v}$ 是質量守恆。
題目問題 4 中的公式：
$\frac{d\Gamma}{dt} = \int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A} + \int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A} + \int_{S} \nabla p \cdot d\mathbf{A} + \int_{S} \rho \mathbf{g} \cdot d\mathbf{A}$
這個公式看起來是將動量守恆方程（N-S 方程）右邊的每一項都積分在一個曲面上，並與環流量的變化率關聯。
通常，環流量的變化率與力（如壓力梯度、重力）和黏滯力有關。

我們將逐一回答問題。

---

(1) 什麼是正壓流體？針對正壓流體改寫上式並說明其意義。

• 正壓流體 (Barotropic Fluid):
  正壓流體是指其密度 $\rho$ 僅是壓強 $p$ 的函數，即 $\rho = \rho(p)$。這意味著在流場中，等密度線（isopycnals）與等壓強線（isobars）是重合的。在這種流體中，壓力梯度力 $\nabla p$ 和密度梯度力 $\nabla \rho$ 是平行的（或反平行的）。

• 改寫上式:
  如果流體是正壓的，則 $\rho = \rho(p)$。
  那麼，$\nabla \rho = \frac{d\rho}{dp} \nabla p$。
  題目給的公式是：
  $\frac{d\Gamma}{dt} = \int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A} + \int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A} + \int_{S} \nabla p \cdot d\mathbf{A} + \int_{S} \rho \mathbf{g} \cdot d\mathbf{A}$
  這個公式的結構看起來不是標準的環流量定理。標準的環流量定理（Kelvin's Theorem）指出：
  $\frac{d\Gamma}{dt} = \oint_C (\nu \nabla^2 \mathbf{v} - \frac{1}{\rho} \nabla p + \mathbf{g}) \cdot d\mathbf{r}$
  其中 $\nu$ 是運動黏滯係數。
  如果題目給的這個公式是正確的，它可能是在一個特定的控制體積或曲面上進行積分，並試圖推導出環流量的變化。

  讓我們假設題目公式的右側每一項代表了影響環流量變化的某種「力」或「通量」。

  • $\int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A}$: 質量通量隨時間的變化率在曲面上的積分。
  • $\int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A}$: 動量通量散度在曲面上的積分（這項通常與動量守恆有關）。
  • $\int_{S} \nabla p \cdot d\mathbf{A}$: 壓力梯度力在曲面上的積分。
  • $\int_{S} \rho \mathbf{g} \cdot d\mathbf{A}$: 重力在曲面上的積分。

  正壓流體的意義:
  在正壓流體中，$\nabla \rho$ 和 $\nabla p$ 是平行的。
  標準的環流量定理 (Kelvin's Theorem) 告訴我們，環流量的變化率 $\frac{d\Gamma}{dt}$ 可以由以下幾項引起：

  1. 黏滯力: $\oint_C \nu \nabla^2 \mathbf{v} \cdot d\mathbf{r}$ (如果存在黏滯)
  2. 非保守力: $\oint_C \nabla (\Phi) \cdot d\mathbf{r}$ (如重力，如果 $\Phi$ 是力勢，則此項為零)
  3. 壓力梯度與密度梯度不平行: $\oint_C \frac{1}{\rho} \nabla p \cdot d\mathbf{r}$ (此項在斜壓流體中是關鍵)

  對於一個無黏滯、無外力（或保守力）的 正壓流體，即 $\rho = \rho(p)$ 且 $\nabla \rho$ 與 $\nabla p$ 平行，並且沒有非保守力，那麼環流量是守恆的，即 $\frac{d\Gamma}{dt} = 0$。
  如果題目給的公式是正確的，並且我們考慮一個無黏滯、無外力的正壓流體，那麼：
  $0 = \int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A} + \int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A} + \int_{S} \nabla p \cdot d\mathbf{A} + \int_{S} \rho \mathbf{g} \cdot d\mathbf{A}$
  在這種情況下，這些積分項的總和必須為零。
  意義: 在正壓流體中，由於密度和壓力具有相同的空間分佈，壓力梯度和密度梯度不會產生額外的「扭轉」效應來改變環流量。如果流體又是無黏滯且受保守力作用，則環流量守恆。

  注意: 題目給的公式形式與標準環流量定理差異較大，我們在此基於標準定理進行解釋。若需嚴格使用題目公式，則需對其積分項做進一步推導和解釋，這可能超出了常規。

---

(2) 何謂斜壓流體？針對斜壓流體，在固定點來觀察，環流量會因為哪些條件（機制或物理過程）而改變？

• 斜壓流體 (Baroclinic Fluid):
  斜壓流體是指其密度 $\rho$ 不是壓強 $p$ 的單值函數，即 $\rho \ne \rho(p)$。這意味著在流場中，等密度線（isopycnals）與等壓強線（isobars）通常是不重合的。
  當等密度線與等壓強線不重合時，壓力梯度力 $\nabla p$ 和密度梯度力 $\nabla \rho$ 通常不平行。這種不平行性會產生一個額外的力，稱為「斜壓項」或「斜壓力」。

• 斜壓流體中環流量的改變:
  對於斜壓流體，環流量的改變主要來自於：
  a. 壓力梯度與密度梯度不平行:
  在標準環流量定理中，有一個項與壓力梯度和密度有關。對於斜壓流體，由於 $\nabla p$ 和 $\nabla \rho$ 不平行，這會導致產生額外的環流量。
  具體來說，考慮環流量 $\Gamma = \oint_C \mathbf{v} \cdot d\mathbf{r}$。
  在斜壓流體中，壓力梯度力 $\nabla p$ 和密度梯度力 $\nabla \rho$ 的作用會產生額外的扭轉效應。
  從動量方程來看，加速度 $\frac{D\mathbf{v}}{Dt}$ 由多項組成，包括壓力梯度、重力、黏滯力等。
  如果我們考慮環流量定理的一個形式：
  $\frac{d\Gamma}{dt} = \oint_C \frac{1}{\rho} \nabla p \cdot d\mathbf{r} - \oint_C \nabla \Phi \cdot d\mathbf{r} + \oint_C \nu \nabla^2 \mathbf{v} \cdot d\mathbf{r}$
  在斜壓流體中，$\nabla p$ 和 $\nabla \rho$ 不平行。這意味著 $\frac{1}{\rho} \nabla p$ 和 $\nabla \rho$（或重力 $\mathbf{g}$，如果 $\rho$ 是 $\mathbf{g}$ 的函數）之間存在夾角。
  這個不平行性是斜壓流體中環流量變化的根本原因。

  **從題目給的公式角度解釋**:
  題目給的公式是：
  $$ \frac{d\Gamma}{dt} = \int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A} + \int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A} + \int_{S} \nabla p \cdot d\mathbf{A} + \int_{S} \rho \mathbf{g} \cdot d\mathbf{A} $$
  在斜壓流體中，壓力梯度 $\nabla p$ 和密度 $\rho$ 的空間分佈不一致。
  *   $\int_{S} \nabla p \cdot d\mathbf{A}$: 壓力梯度力在曲面上的積分。
  *   $\int_{S} \rho \mathbf{g} \cdot d\mathbf{A}$: 重力在曲面上的積分。
  
  **關鍵點**: 斜壓性意味著，即使在沒有黏滯力或時間變化項的情況下，壓力梯度力 $\nabla p$ 和重力 $\rho \mathbf{g}$ 的組合作用，也可能因為 $\rho$ 和 $p$ 的空間分佈不一致而產生淨的扭轉效應，改變環流量。
  
  **機制**:
  1.  **壓力梯度與密度梯度不平行**: 這是斜壓性的定義，也是環流量改變的根本原因。
  2.  **時間變化**: $\int_{S} \frac{\partial (\rho \mathbf{v})}{\partial t} \cdot d\mathbf{A}$ (質量通量隨時間變化)。
  3.  **動量通量變化**: $\int_{S} \nabla \cdot (\rho \mathbf{v} \mathbf{v}) \cdot d\mathbf{A}$ (動量通量的淨散度)。
  4.  **壓力梯度力**: $\int_{S} \nabla p \cdot d\mathbf{A}$。
  

此題考驗對層流（Laminar Flow）與紊流（Turbulent Flow）基本概念、雷諾數（Reynolds Number）的應用，以及相似性原理（Similitude）的理解。

---

(a) 請說明 Laminar flow 及 Turbulent flow 的定義。

• Laminar Flow (層流):
  層流是一種流動狀態，其中流體質點沿著平滑、穩定的路徑（稱為流線，streamlines）運動。在層流中，相鄰的流體層之間只有分子擴散作用產生剪應力，沒有顯著的質點混合。流動是規整、有序且可預測的。在低流速、高黏度或小尺度通道中，流動傾向於呈現層流狀態。

  • 特徵:
    • 流線平滑、平行。
    • 流體質點沿著明確的路徑移動。
    • 缺乏顯著的質點混合（macroscopic mixing）。
    • 剪應力主要由分子黏滯引起。
    • 流動是穩定的（steady）或隨時間緩慢變化的。

• Turbulent Flow (紊流):
  紊流是一種高度無序、隨機且混亂的流動狀態，其特徵是流體質點的運動軌跡不規則，並伴隨有劇烈的質點混合和能量耗散。在紊流中，流體會形成大小不一的渦旋（eddies），這些渦旋在空間和時間上都在不斷變化。紊流通常發生在高流速、低黏度或大尺度流動中。

  • 特徵:
    • 流線不規則、隨機，不斷改變形狀。
    • 流體質點運動軌跡混亂，有顯著的橫向（垂直於主流方向）運動。
    • 存在顯著的質點混合（macroscopic mixing），導致動量、熱量和質量在流場中快速傳遞。
    • 剪應力除了分子黏滯外，還有由於渦旋運動引起的額外「紊流剪應力」（turbulent shear stress），通常遠大於分子黏滯剪應力。
    • 流動本質上是隨機和不穩定的（unsteady），即使平均速度場可能是穩定的。

---

(b) 請以 Reynolds number 區分 Turbulent flow 及 Laminar flow。

Reynolds Number (雷諾數, Re)
雷諾數是一個無量綱的參數，用於描述流體慣性力與黏滯力之間的相對大小。它對於預測流動狀態（層流或紊流）至關重要。
$Re = \frac{\text{慣性力}}{\text{黏滯力}} = \frac{\rho v L}{\mu} = \frac{v L}{\nu}$
其中：

• $\rho$: 流體密度。
• $v$: 流體的特徵速度（例如平均速度）。
• $L$: 流體的特徵長度（例如管道直徑、物體尺寸）。
• $\mu$: 流體的動力黏滯係數。
• $\nu = \mu/\rho$: 流體的運動黏滯係數。

區分標準:
層流和紊流之間的轉換通常發生在一個特定的雷諾數範圍內，這個範圍取決於具體的幾何形狀和流動條件。

• 管道內流動 (Pipe Flow):

  • $Re < 2300$: 通常是 層流。
  • $2300 < Re < 4000$: 過渡區（Transition Region），流動可能呈現層流和紊流的混合狀態，不穩定。
  • $Re > 4000$: 通常是 紊流。

• 外部流動（例如，平板上的邊界層）:

  • 層流邊界層在 $Re_x \approx 5 \times 10^5$ 左右開始轉捩為紊流邊界層（$Re_x = \frac{\rho v x}{\mu}$，其中 $x$ 是從平板前緣開始的距離）。

總結: 雷諾數提供了一個量化的指標來判斷流動是傾向於有序的層流還是混亂的紊流。較低的雷諾數意味著黏滯力佔主導，流動趨向於層流；較高的雷諾數意味著慣性力佔主導，流動趨向於紊流。

---

(c) 兩個流體系統具有動力相似性 (Dynamical Similarity) 的條件為何？

動力相似性是指兩個流體系統（一個是原型，另一個是模型）在運動和受力上是相似的。如果兩個系統具有動力相似性，那麼在模型上進行的實驗結果就可以用來預測原型的行為。

要使兩個流體系統（原型 P 和模型 M）具有動力相似性，必須滿足以下兩個條件：

此題考驗對向量微積分中梯度（gradient）概念的理解，以及其在流體力學中的物理意義。題目要求證明一個向量場在空間中微小位移下的變化，並解釋其中梯度項的含義。

---

證明 $d\mathbf{v} = \frac{\partial \mathbf{v}}{\partial x} dX + \frac{\partial \mathbf{v}}{\partial y} dY$

這實際上是向量場在笛卡爾座標系下的全微分（total differential）表達式。
假設我們有一個二維向量場 $\mathbf{v}(x,y) = \langle u(x,y), v(x,y) \rangle$。
我們想計算在點 $(X_0, Y_0)$ 附近，當位置從 $(X_0, Y_0)$ 變化到 $(X_0 + \delta X, Y_0 + \delta Y)$ 時，向量場 $\mathbf{v}$ 的變化量 $d\mathbf{v}$。
這裡，我們將 $x$ 視為 $X_0$ 的變量，$y$ 視為 $Y_0$ 的變量。
題目中給定的 $dX = \delta X$ 和 $dY = \delta Y$ 是位置的微小變化量。

根據向量微積分的鏈式法則，向量場 $\mathbf{v}$ 的全微分 $d\mathbf{v}$ 可以表示為：
$d\mathbf{v} = \left( \frac{\partial \mathbf{v}}{\partial x} \right) dx + \left( \frac{\partial \mathbf{v}}{\partial y} \right) dy$
這裡，我們將 $dx$ 替換為題目中的 $dX$（或 $\delta X$），將 $dy$ 替換為題目中的 $dY$（或 $\delta Y$）。

所以，
$d\mathbf{v} = \frac{\partial \mathbf{v}}{\partial x} dX + \frac{\partial \mathbf{v}}{\partial y} dY$
這與題目要求證明的形式一致。

更詳細的展開:
設 $\mathbf{v}(x,y) = \langle u(x,y), v(x,y) \rangle$。
那麼 $d\mathbf{v} = \langle du, dv \rangle$。
根據多元函數的全微分公式：
$du = \frac{\partial u}{\partial x} dx + \frac{\partial u}{\partial y} dy$
$dv = \frac{\partial v}{\partial x} dx + \frac{\partial v}{\partial y} dy$

將 $dx$ 替換為 $dX$，$dy$ 替換為 $dY$：
$du = \frac{\partial u}{\partial x} dX + \frac{\partial u}{\partial y} dY$
$dv = \frac{\partial v}{\partial x} dX + \frac{\partial v}{\partial y} dY$

所以，
$d\mathbf{v} = \langle du, dv \rangle = \left\langle \frac{\partial u}{\partial x} dX + \frac{\partial u}{\partial y} dY, \frac{\partial v}{\partial x} dX + \frac{\partial v}{\partial y} dY \right\rangle$

現在，我們來看 $\frac{\partial \mathbf{v}}{\partial x}$ 和 $\frac{\partial \mathbf{v}}{\partial y}$ 這兩項。
$\frac{\partial \mathbf{v}}{\partial x}$ 是向量場 $\mathbf{v}$ 對 $x$ 的偏導數。將 $\mathbf{v} = \langle u, v \rangle$ 代入：
$\frac{\partial \mathbf{v}}{\partial x} = \left\langle \frac{\partial u}{\partial x}, \frac{\partial v}{\partial x} \right\rangle$
同理，
$\frac{\partial \mathbf{v}}{\partial y} = \left\langle \frac{\partial u}{\partial y}, \frac{\partial v}{\partial y} \right\rangle$

將這些代回 $d\mathbf{v}$ 的表達式：
$d\mathbf{v} = \left\langle \frac{\partial u}{\partial x}, \frac{\partial v}{\partial x} \right\rangle dX + \left\langle \frac{\partial u}{\partial y}, \frac{\partial v}{\partial y} \right\rangle dY$
$d\mathbf{v} = \frac{\partial \mathbf{v}}{\partial x} dX + \frac{\partial \mathbf{v}}{\partial y} dY$
證明完畢。

---

說明 $\frac{\partial \mathbf{v}}{\partial x}$ 與 $\frac{\partial \mathbf{v}}{\partial y}$ 在數學上與物理上各代表什麼？

數學上 (Mathematically):

• $\frac{\partial \mathbf{v}}{\partial x}$:
  這是向量場 $\mathbf{v}(x,y)$ 對於空間變數 $x$ 的偏導數。它是一個向量，其 $i$-th 分量是 $\mathbf{v}$ 的 $i$-th 分量對 $x$ 的偏導數。
  若 $\mathbf{v} = \langle u, v \rangle$，則 $\frac{\partial \mathbf{v}}{\partial x} = \langle \frac{\partial u}{\partial x}, \frac{\partial v}{\partial x} \rangle$。
  它衡量了在 $y$ 座標固定時，向量場 $\mathbf{v}$ 沿著 $x$ 方向變化的「率」或「方向」。

• $\frac{\partial \mathbf{v}}{\partial y}$:
  這是向量場 $\mathbf{v}(x,y)$ 對於空間變數 $y$ 的偏導數。它是一個向量，其 $i$-th 分量是 $\mathbf{v}$ 的 $i$-th 分量對 $y$ 的偏導數。
  若 $\mathbf{v} = \langle u, v \rangle$，則 $\frac{\partial \mathbf{v}}{\partial y} = \langle \frac{\partial u}{\partial y}, \frac{\partial v}{\partial y} \rangle$。
  它衡量了在 $x$ 座標固定時，向量場 $\mathbf{v}$ 沿著 $y$ 方向變化的「率」或「方向」。

這兩個向量 $\frac{\partial \mathbf{v}}{\partial x}$ 和 $\frac{\partial \mathbf{v}}{\partial y}$ 共同構成了向量場 $\mathbf{v}$ 在空間中的「梯度」或「變形率」。

此題考驗對極座標下流體運動學的理解，特別是自由渦旋的特性，以及如何計算流函數、速度位、渦度和環流量。

1. 自由渦旋 (Free Vortex) 的定義與速度場
題目給定的自由渦旋在極座標 $(r, \theta)$ 下的速度場是：
$v_r = 0$ (徑向速度為零)
$v_\theta = \frac{\Gamma}{2\pi r}$ (切向速度與半徑 $r$ 成反比，$\Gamma$ 是常數)

這表示流體質點僅繞原點做圓周運動，速度大小隨半徑增大而減小。

• $\Gamma$ 是環流量 (Circulation)。
• 原點 $(r=0)$ 處速度趨於無限大，是一個奇異點 (singularity)。

---

2. 求流函數 (Stream Function, $\psi$)

在二維極座標下，流函數 $\psi$ 與速度分量的關係為：
$v_r = \frac{1}{r} \frac{\partial \psi}{\partial \theta}$
$v_\theta = -\frac{\partial \psi}{\partial r}$

我們已知 $v_r = 0$ 和 $v_\theta = \frac{\Gamma}{2\pi r}$。

• 由 $v_r = 0$：
  $\frac{1}{r} \frac{\partial \psi}{\partial \theta} = 0$
  這意味著 $\frac{\partial \psi}{\partial \theta} = 0$。
  因此，$\psi$ 不依賴於 $\theta$，只是一個 $r$ 的函數：$\psi = \psi(r)$。

• 由 $v_\theta = -\frac{\partial \psi}{\partial r}$：
  $\frac{\Gamma}{2\pi r} = -\frac{\partial \psi}{\partial r}$
  $\frac{\partial \psi}{\partial r} = -\frac{\Gamma}{2\pi r}$

  對上式進行積分以求解 $\psi(r)$:
  $\psi(r) = \int -\frac{\Gamma}{2\pi r} dr$
  $\psi(r) = -\frac{\Gamma}{2\pi} \ln(r) + C$
  其中 $C$ 是積分常數。在流函數的定義中，常數項可以任意選擇，通常設為零以簡化。
  所以，流函數為：
  $\psi(r) = -\frac{\Gamma}{2\pi} \ln(r)$

  （注意：有時流函數的定義符號可能不同，例如 $v_\theta = \frac{\partial \psi}{\partial r}$，此時 $\psi(r) = \frac{\Gamma}{2\pi} \ln(r)$。我們採用 $v_\theta = -\frac{\partial \psi}{\partial r}$ 的定義。）

---

3. 求速度位 (Velocity Potential, $\phi$)

速度位 $\phi$ 存在於無旋流 (irrotational flow) 中，其定義為 $\mathbf{v} = \nabla \phi$。
在二維極座標下，速度分量與速度位 $\phi(r, \theta)$ 的關係為：
$v_r = \frac{\partial \phi}{\partial r}$
$v_\theta = \frac{1}{r} \frac{\partial \phi}{\partial \theta}$

• 由 $v_r = 0$:
  $\frac{\partial \phi}{\partial r} = 0$
  這意味著 $\phi$ 不依賴於 $r$，只是一個 $\theta$ 的函數：$\phi = \phi(\theta)$。

• 由 $v_\theta = \frac{\Gamma}{2\pi r}$:
  $\frac{1}{r} \frac{\partial \phi}{\partial \theta} = \frac{\Gamma}{2\pi r}$
  $\frac{\partial \phi}{\partial \theta} = \frac{\Gamma}{2\pi}$

  對上式進行積分以求解 $\phi(\theta)$:
  $\phi(\theta) = \int \frac{\Gamma}{2\pi} d\theta$
  $\phi(\theta) = \frac{\Gamma}{2\pi} \theta + D$
  其中 $D$ 是積分常數，通常設為零。
  所以，速度位為：
  $\phi(\theta) = \frac{\Gamma}{2\pi} \theta$

  檢查是否為無旋流:
  自由渦旋的渦度 $\boldsymbol{\omega} = \nabla \times \mathbf{v}$。
  在極座標下，對於二維向量場 $\mathbf{v} = (v_r, v_\theta)$，渦度（垂直於 $r-\theta$ 平面的 $z$ 分量）為：
  $\omega_z = \frac{1}{r} \frac{\partial}{\partial r} (r v_\theta) - \frac{1}{r} \frac{\partial v_r}{\partial \theta}$
  代入 $v_r = 0$ 和 $v_\theta = \frac{\Gamma}{2\pi r}$: