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

核心觀念

本題考查流體力學中最基本的流動型態分類：層流（Laminar flow） 與 紊流／湍流（Turbulent flow） 的物理定義、流動特徵與判別準則。

• 層流（Laminar Flow）：流體質點呈現平滑、有序且分層滑動的流動狀態，層與層之間無宏觀的流體團塊混摻，質量、動量與能量的傳遞主要依賴分子擴散（Molecular diffusion）。
• 紊流（Turbulent Flow）：流體質點呈現高度無序、隨機脈動（Random fluctuations）與混亂漩渦（Eddies）的流動狀態，強烈的宏觀渦流混摻主導了動量、熱量與質量的傳遞。
• 流態判別準則：流態通常由**雷諾數（Reynolds number, $Re$）**決定，代表慣性力（Inertial force）與黏滯力（Viscous force）的相對比值：
  $Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$
  其中 $\rho$ 為流體密度，$U$ 為特徵流速，$L$ 為特徵長度，$\mu$ 為動力黏滯係數，$\nu$ 為運動黏滯係數。當 $Re$ 較小（黏滯力占優勢，能抑制流場擾動）時為層流；當 $Re$ 超過臨界雷諾數 $Re_{cr}$（慣性力占優勢，微小擾動被放大失穩）時流場過渡轉化為紊流。

---

解題方法

本題為名詞定義與觀念說明題，切入重點在於：

1. 精準定義：從質點運動模式、層狀結構及是否有橫向宏觀脈動切入。
2. 物理機制與傳遞特性對比：對比兩者的動量傳遞機制（分子擴散 vs. 渦流混合/雷諾應力）。
3. 無因次參數關聯：導入雷諾數 $Re$ 解釋外在流況如何決定流動型態。

---

詳細說明與分析

1. 層流（Laminar Flow）

• 定義：流體質點沿著平滑且平行的路徑（流線）呈層狀滑動的流態。各流層之間互不干擾交錯，不存在垂直於主流方向的宏觀混合或速度脈動。
• 主要特徵：
  1. 高度有序性：質點運動規則、速度場具確定性（Deterministic），可精確求解（如 Hagen-Poiseuille 流）。
  2. 傳遞機制：剪應力僅由分子黏滯性（Newton's law of viscosity）主導：
     $\tau = \mu \frac{\partial u}{\partial y}$
     動量、質量與熱量的橫向交換速率極慢，僅依賴分子擴散尺度。

核心觀念

1. 雷諾數（Reynolds number, $Re$）的定義與物理意義：
   雷諾數為流體力學中最重要的無因次參數，代表流體流動時「慣性力（Inertial force）」與「黏滯力（Viscous force）」的比值：
   $Re = \frac{\text{慣性力}}{\text{黏滯力}} = \frac{\rho V L}{\mu} = \frac{V L}{\nu}$
   其中：

   • $\rho$：流體密度（Fluid density, $\mathrm{kg/m^3}$）
   • $V$：特徵流速（Characteristic velocity, $\mathrm{m/s}$）
   • $L$：特徵長度（Characteristic length, $\mathrm{m}$）
   • $\mu$：動力黏度（Dynamic viscosity, $\mathrm{Pa \cdot s}$）
   • $\nu = \frac{\mu}{\rho}$：運動黏度（Kinematic viscosity, $\mathrm{m^2/s}$）

2. 流態穩定性機制：

   • 慣性力：傾向於放大流場中的隨機微小擾動，引發渦旋與流動失穩。
   • 黏滯力：流體分子間的內摩擦力，具耗散（dissipation）作用，傾向於抑制並平撫流場擾動。

---

解題方法

區分層流（Laminar flow）與紊流（Turbulent flow）的核心依據為臨界雷諾數（Critical Reynolds number, $Re_{cr}$），具體步驟與說明如下：

1. 物理機制的判定準則

• 層流（Laminar Flow）：
  當 $Re < Re_{cr}$ 時，黏滯力遠大於慣性力。流場中的微小擾動會被流體的黏滯剪應力迅速衰減並平撫，流體粒子呈平滑、有序的分層流動（Streamlines 平行且無垂直方向的宏觀混合）。
• 紊流（Turbulent Flow）：
  當 $Re > Re_{cr}$ 時，慣性力占絕對主導地位。微小擾動會被慣性力非線性放大，導致流動失穩，形成大小不一的渦旋（Eddies），流體質點呈現隨機、混沌且強烈的三維脈動與混合現象。

2. 結合流動幾何形狀與臨界值判斷

雷諾數的數值高低與所選取的「特徵長度 $L$」密切相關，不同流場幾何具有不同的臨界雷諾數作為區分基準：

核心觀念

本題評量流體力學中的相似性原理（Similitude / Similarity Requirements）。

在模型試驗（Model testing）或比較兩個流體流動系統（原型 Prototype 與模型 Model）時，流場的相似性可分為三個層次：

1. 幾何相似（Geometric Similarity）：兩系統之所有對應線性幾何尺寸比例恆定，且對應夾角相同。
2. 運動相似（Kinematic Similarity）：兩系統在幾何相似的前提下，所有對應點的速度向量方向相同、量值比例恆定（流線形態完全相同），亦即時間尺度與速度場成固定比例。
3. 動力相似（Dynamic / Dynamical Similarity）：兩系統在幾何相似與運動相似的前提下，所有對應點所承受的各種外力向量（如慣性力、黏滯力、壓力、重力、表面張力、彈性力等）方向相同、大小比例恆為定值。

若由統御方程式（Navier-Stokes 方程式）的無因次化觀點切入，達到動力相似的充要條件為：兩系統的無因次控制參數（相似參數 / 無因次群）完全相同。

---

解題方法

說明動力相似性的條件，應由「受力比例」與「無因次參數匹配」兩個層次完整論述：

1. 力比恆定條件（Force Ratio Condition）

兩系統達到動力相似的物理定義為：在對應的幾何點上，所有主導流動的力（Forces）之比值必須相同。
即對應的慣性力（$F_I$）、壓力（$F_p$）、黏性力（$F_v$）、重力（$F_g$）等外力之間，滿足：

$$\left(\frac{F_I}{F_v}\right)_m = \left(\frac{F_I}{F_v}\right)_p, \quad \left(\frac{F_I}{F_g}\right)_m = \left(\frac{F_I}{F_g}\right)_p, \quad \left(\frac{F_I}{F_p}\right)_m = \left(\frac{F_I}{F_p}\right)_p, \quad \dots$$

（式中下標 $m$ 代表模型 Model，$p$ 代表原型 Prototype）

2. 必要前提與無因次特徵參數相等的條件

要達到動力相似性，必須同時滿足下列條件：

• 前提條件：
  1. 兩系統必須具有幾何相似（Geometric similarity）：幾何形狀與邊界條件比例必須一致。
  2. 兩系統必須具有運動相似（Kinematic similarity）：速度場幾何形態一致（事實上，在純動力系統中，幾何相似加上力比恆定即可推導出運動相似）。
• 實務充要條件（無因次參數相符）：
  控制流動物理現象的所有關鍵無因次參數必須在兩系統中分別相等。常見重要無因次參數包括：
  1. 雷諾數相等（Reynolds Number, $Re$）：慣性力與黏滯力之比 $$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu} \implies Re_m = Re_p$$

核心觀念

本題旨在考驗流體力學中控制方程式的無因次化（Nondimensionalization），以及由此衍生出的**重要無因次參數（Dimensionless Parameters）**之物理意義與數學表示式。

1. 無因次化的目的：藉由選取流場具代表性的特徵尺度（長度 $L$、速度 $U$、特徵時間 $L/U$、特徵壓力差 $\Delta P$ 或特徵壓力 $P_0$），將各物理量正規化。這不僅能縮減獨立變數數量、利於動態相似性（Dynamic Similarity）與模型比例放大縮小實驗，更能透過各項係數評估各作用力相對大小，以利合理簡化方程式。
2. 無因次參數本質：通常表示為**慣性力（Inertial force）**與特定作用力（黏滯力、重力、壓力）的比值。

---

解題方法

1. 無因次變數設定與方程式轉換

根據題目設定，選取下列特徵尺度進行無因次化：

• 空間座標與微分算子：$\vec{x}^* = \frac{\vec{x}}{L} \implies \nabla^* = L\nabla$，$\nabla^{*2} = L^2\nabla^2$
• 速度場：$\vec{u}^* = \frac{\vec{u}}{U} \implies \vec{u} = U\vec{u}^*$
• 時間：$t^* = \frac{t}{L/U} \implies t = \frac{L}{U}t^*$
• 壓力：$P^* = \frac{P}{P_0}$（其中 $P_0$ 為特徵壓力）
• 重力單位向量：$\frac{\vec{g}}{g}$（$g = |\vec{g}|$ 為重力加速度純量）

原 Navier-Stokes 方程式（註：原式中黏滯項係數以動黏滯係數 $\nu = \mu/\rho$ 撰寫更為精確，依題意保留）：

$$\frac{d\vec{u}}{dt} = -\frac{1}{\rho}\nabla P + \vec{g} + \nu\nabla^2\vec{u}$$

將各無因次變數代入各項：

• 質點加速度項（慣性項）：

$$\frac{d\vec{u}}{dt} = \frac{U}{L/U}\frac{d\vec{u}^*}{dt^*} = \frac{U^2}{L}\frac{d\vec{u}^*}{dt^*}$$

• 壓力梯度項：

$$-\frac{1}{\rho}\nabla P = -\frac{1}{\rho}\frac{1}{L}\nabla^*(P_0 P^*) = -\frac{P_0}{\rho L}\nabla^* P^*$$

• 體積力（重力）項：

$$\vec{g} = g\frac{\vec{g}}{g}$$

• 黏滯擴散項：

$$\nu\nabla^2\vec{u} = \nu\left(\frac{1}{L^2}\nabla^{*2}\right)(U\vec{u}^*) = \frac{\nu U}{L^2}\nabla^{*2}\vec{u}^*$$

將各項組合後同除以慣性力特徵尺度 $\frac{U^2}{L}$：

$$\frac{d\vec{u}^*}{dt^*} = -\left(\frac{P_0}{\rho U^2}\right)\nabla^* P^* + \left(\frac{g L}{U^2}\right)\frac{\vec{g}}{g} + \left(\frac{\nu}{U L}\right)\nabla^{*2}\vec{u}^*$$

比對題目給定之方程式：

$$\frac{d\vec{u}^*}{dt^*} = -\frac{1}{E}\nabla^* P^* + \frac{1}{Fr^2}\frac{\vec{g}}{g} + \frac{1}{Re}\nabla^{*2}\vec{u}^*$$

即可直接比對出各無因次係數之數學定義。

---

各無因次參數之數學表示式與物理意義

1. Reynolds number（雷諾數，記作 $Re$）

• 數學表示式：

$$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$$

• 物理意義：
  • 慣性力（Inertial force）與黏滯力（Viscous force）之比：

$$Re = \frac{\text{Inertia force}}{\text{Viscous force}} \sim \frac{\rho U^2 / L}{\mu U / L^2} = \frac{\rho U L}{\mu}$$

• 應用意義：判斷流場為層流（Laminar flow）或紊流（Turbulent flow）的關鍵指標。

核心觀念

本題考查流體力學中的兩大基本流場假設與定義：

1. 無旋性流（Irrotational Flow）與位勢流（Potential Flow）：涉及流體微元的旋轉狀態（渦度，Vorticity）以及速度位函數（Velocity Potential, $\phi$）之存在條件。
2. 不可壓縮流體（Incompressible Fluid）：涉及流體微團在運動過程中的密度變化行為，以及連續方程式（Continuity Equation）在速度場散度（Divergence）上的簡化。

---

解題方法與觀念詳解

一、無旋性流（Irrotational Flow）與位勢流（Potential Flow）

1. 定義與物理意義：
   • 無旋性流：流場中各處流體微元（fluid element）自身的旋轉角速度向量 $\boldsymbol{\omega}$ 恆為零（即無自轉現象，純粹產生平移與剪切變形）。
   • 數學條件：定義渦度向量（Vorticity Vector）為速度場 $\mathbf{V}$ 的旋度（Curl）： $$\boldsymbol{\zeta} = \nabla \times \mathbf{V} = 2\boldsymbol{\omega} = 0$$
2. 位勢流之關聯：
   • 根據向量微積分定理，任何無旋向量場（旋度為零）皆可表示為某個純量函數的梯度（Gradient）。
   • 因此，存在一純量速度位函數（Velocity Potential Function）$\phi$，使得速度場可寫為： $$\mathbf{V} = \nabla \phi$$ （亦有部分教科書定義為 $\mathbf{V} = -\nabla \phi$）。
   • 滿足此條件的流場即稱為位勢流（Potential Flow）。故在流體力學中，「無旋流」與「位勢流」兩詞常等價互換。

---

二、不可壓縮流體（Incompressible Fluid）

1. 定義與物理意義：
   • 指流體微元在隨流運動過程中，其密度 $\rho$ 保持為常數，即流體的實質導數（Material Derivative）為零： $$\frac{D\rho}{Dt} = \frac{\partial \rho}{\partial t} + \mathbf{V} \cdot \nabla \rho = 0$$
   • 若流體處於均質（homogeneous）狀態，則全場之密度 $\rho = \text{常數}$。
2. 控制方程式（連續方程式）之簡化：
   • 通用質量守恆方程式（連續方程式）為：

核心觀念

本題考查流體力學中位流理論（Potential Flow Theory）的兩大基礎純量函數——**速度位（Velocity Potential, $\phi$）與流函數（Stream Function, $\psi$）**之定義、物理存在條件及其在二維笛卡兒座標系（Cartesian Coordinates）下與速度向量分量之關係：

1. 速度位（Velocity Potential, $\phi$）：

   • 存在條件：流場必須為無旋流（Irrotational Flow），即旋度為零： $$\nabla \times \vec{V} = 0$$
   • 根據向量微積分定理，旋度處處為零之向量場可表示為某一純量函數之梯度： $$\vec{V} = \nabla \phi$$ （註：部分教科書採用 $\vec{V} = -\nabla \phi$，但在流體力學中標準慣用為正號定義 $\vec{V} = \nabla \phi$）。

2. 流函數（Stream Function, $\psi$）：

   • 存在條件：流場必須滿足**不可壓縮流（Incompressible Flow）**之連續方程式（Continuity Equation）： $$\nabla \cdot \vec{V} = 0$$
   • 在二維平面流場中，連續方程式自然滿足流函數的全微分形式，保證質量守恆。

---

解題方法與完整推導

設二維平面卡氏座標系為 $(x, y)$，速度向量為：

$$\vec{V} = u(x, y)\,\hat{i} + v(x, y)\,\hat{j}$$

其中 $u$ 為 $x$ 方向速度分量，$v$ 為 $y$ 方向速度分量。

一、速度分量與速度位 $\phi$ 之關係推導

流體若為無旋流，渦度向量（Vorticity）$\vec{\omega}$ 為零：

$$\vec{\omega} = \nabla \times \vec{V} = \left( \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \right)\hat{k} = 0 \implies \frac{\partial v}{\partial x} = \frac{\partial u}{\partial y}$$

因為旋度為零，速度向量 $\vec{V}$ 必可寫為純量速度位 $\phi(x, y)$ 之梯度：

$$\vec{V} = \nabla \phi = \frac{\partial \phi}{\partial x}\hat{i} + \frac{\partial \phi}{\partial y}\hat{j}$$

將速度向量之卡氏分量與梯度分量對照比較，即得速度分量與速度位的關係式：

$$u = \frac{\partial \phi}{\partial x}$$ $$v = \frac{\partial \phi}{\partial y}$$

二、速度分量與流函數 $\psi$ 之關係推導

對於二維不可壓縮連續流體，二維連續方程式為：

$$\nabla \cdot \vec{V} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$

為使此偏微分方程式自動恆等成立，定義流函數 $\psi(x, y)$，使得：

核心觀念

本題考查**二維理想流體（二維、無黏性、不可壓縮且無旋性流動，即勢流 Potential Flow）的基本理論與流網（Flow Net）**的分析應用：

1. 無旋性（Irrotationality）：速度場滿足 $\nabla \times \vec{V} = 0$，保證存在標量速度勢函數（Velocity Potential）$\phi(x, y)$，使得速度向量可表示為勢函數的梯度： $$\vec{V} = \nabla \phi = \frac{\partial \phi}{\partial x}\hat{i} + \frac{\partial \phi}{\partial y}\hat{j}$$ 故速度分量為 $u = \frac{\partial \phi}{\partial x}$、$v = \frac{\partial \phi}{\partial y}$。
2. 不可壓縮連續方程式（Incompressibility）：二維不可壓縮流滿足 $\nabla \cdot \vec{V} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$，保證存在流線函數（Stream Function）$\psi(x, y)$，定義為： $$u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}$$
3. 等值線定義：
   • 流線（Streamline）：$\psi(x, y) = \text{常數}$，切線方向隨處與流速向量 $\vec{V}$ 平行。
   • 等勢線（Equipotential line）：$\phi(x, y) = \text{常數}$，法線方向隨處沿勢能梯度方向 $\nabla \phi$。

   符號說明：題目敘述中之「等 $\phi$ 線（streamline）」與「等 $\varphi$ 線（Equipotential line）」及後續「$\phi$ 及 $\psi$」，係流體力學文獻中常見的符號互換或印製筆誤。以下統一以國際及台灣工程力學標準通用符號：速度勢 $\phi$（等勢線 $\phi = \text{常數}$）與流函數 $\psi$（流線 $\psi = \text{常數}$）進行嚴謹論證。

---

解題方法

1. 物理觀點證明（Orthogonality from Physical Viewpoint）

• 流線的物理意義：
  依流線定義，流線上一點的切線方向即為該點的流速向量 $\vec{V} = (u, v)$ 方向，因此流線的切向單位向量為：

  $$\hat{t}_{\psi} = \frac{\vec{V}}{|\vec{V}|}$$

  流體粒子沿流線切線方向運動，法向沒有跨越流線的流動（流線不可穿透，$\vec{V} \cdot \hat{n}_{\psi} = 0$）。

• 等勢線的物理意義：
  由於 $\vec{V} = \nabla \phi$，在多元微積分與場論中，任意純量場 $\phi$ 的梯度向量 $\nabla \phi$ 必然處處正交（垂直）於其等值面（或二維下的等勢線 $\phi = \text{常數}$）。因此，等勢線在任意點的法線方向即為速度向量 $\vec{V}$ 的方向：

  $$\hat{n}_{\phi} = \frac{\nabla \phi}{|\nabla \phi|} = \frac{\vec{V}}{|\vec{V}|}$$

• 正交性結論：
  因為等勢線的法向向量 $\hat{n}_{\phi}$ 與流線的切向向量 $\hat{t}_{\psi}$ 同為 $\frac{\vec{V}}{|\vec{V}|}$（完全平行），這表示流線的切線方向正好垂直於等勢線的切線方向。

  因此，在流速不為零（$\vec{V} \neq 0$）的區域，等勢線必定處處與流線相互垂直。

---

2. 數學觀點證明（Orthogonality from Mathematical Viewpoint）

可由「曲線斜率乘積為 $-1$」與「法向量內積為 $0$」兩種方式嚴謹推導：

方法一：全微分與斜率法

• 對等勢線 $\phi(x, y) = C_1$ 取全微分：

  $$\mathrm{d}\phi = \frac{\partial \phi}{\partial x} \mathrm{d}x + \frac{\partial \phi}{\partial y} \mathrm{d}y = u\,\mathrm{d}x + v\,\mathrm{d}y = 0$$

  整理可得等勢線的幾何切線斜率：

  $$\left( \frac{\mathrm{d}y}{\mathrm{d}x} \right)_{\phi = \text{const}} = -\frac{u}{v} \quad (v \neq 0)$$

• 對流線 $\psi(x, y) = C_2$ 取全微分：

  $$\mathrm{d}\psi = \frac{\partial \psi}{\partial x} \mathrm{d}x + \frac{\partial \psi}{\partial y} \mathrm{d}y = -v\,\mathrm{d}x + u\,\mathrm{d}y = 0$$

  整理可得流線的幾何切線斜率：

  $$\left( \frac{\mathrm{d}y}{\mathrm{d}x} \right)_{\psi = \text{const}} = \frac{v}{u} \quad (u \neq 0)$$

• 計算兩切線斜率乘積：

  $$\left( \frac{\mathrm{d}y}{\mathrm{d}x} \right)_{\phi = \text{const}} \times \left( \frac{\mathrm{d}y}{\mathrm{d}x} \right)_{\psi = \text{const}} = \left( -\frac{u}{v} \right) \left( \frac{v}{u} \right) = -1$$

核心觀念

本題考查流體力學與熱力學中機械能（動能）方程式與熱能（內能）方程式的物理意義，以及兩種能量形式之間的能量轉換機制（Coupling / Exchange Mechanism）。

1. 機械能方程式（Kinetic Energy Equation）：由柯西動量方程式（Navier-Stokes 方程式）與速度向量 $\vec{u}$ 取內積（$\vec{u}\cdot\text{Momentum Eq.}$）純量化推導而得，描述單位體積流體微元動能的實質時變率（Substantial derivative）。
2. 內能方程式（Thermal/Internal Energy Equation）：由熱力學第一定律（總能量守恆方程式）減去機械能方程式而得，描述單位質量（或單位體積）內能 $e$ 的實質變化。
3. 能量轉換機制（Coupling Terms）：兩式中大小相等但正負號相反的交互項，正是動能與內能相互轉換的橋樑：
   • 壓力膨脹／壓縮功項（Reversible work）：$P(\nabla\cdot\vec{V})$
   • 黏滯耗散項（Viscous dissipation）：$\phi$

(註：題目的第二式開頭 $P\frac{de}{dt}$ 為考卷排版之筆誤，物理定義上應為密度 $\rho\frac{de}{dt}$，解說時依標準流體熱力學物理意義剖析。)

---

解題方法

本題為觀念推導與原理解釋題，解題切入點分為三大部分：

1. 逐項剖析式 (1) 各項物理意義（動能的產生、傳輸、轉換與耗散）。
2. 逐項剖析式 (2) 各項物理意義（內能的變化、熱傳導、做功轉換與黏滯加熱）。
3. 明確指認出連接兩方程式的轉換項（$P\nabla\cdot\vec{V}$ 與 $\phi$），並以可逆與不可逆性詳細說明轉換機制。

---

詳細解析

一、方程式 (1)：機械動能方程式之物理意義

$$\rho\frac{d}{dt}\left(\frac{1}{2}u_i^2\right) = \rho\vec{g}\cdot\vec{u} + \frac{\partial}{\partial x_j}(u_i\tau_{ij}) + P(\nabla\cdot\vec{V}) - \phi$$

此式描述流體微元之動能隨體積變化率，各項意義如下：

1. 左式 $\rho\frac{d}{dt}\left(\frac{1}{2}u_i^2\right)$：
   隨體流動的單位體積流體動能的實質時變率（Material derivative of kinetic energy per unit volume）。
2. 右式第 1 項 $\rho\vec{g}\cdot\vec{u}$：
   重力做功功率（Work done by body force per unit volume）。重力作用在流體上，將重力位能與動能進行可逆轉換。
3. 右式第 2 項 $\frac{\partial}{\partial x_j}(u_i\tau_{ij})$：
   黏滯表面應力功的空間發散項（傳輸項）。代表由表面黏滯剪應力所傳遞進出流體微元邊界的淨功率（Net rate of work done by surface viscous stress）。
4. 右式第 3 項 $P(\nabla\cdot\vec{V})$：
   壓力做功與流體膨脹／壓縮之耦合項。
   (註：在標準推導中，壓力梯度做功 $-\vec{u}\cdot\nabla P = -\nabla\cdot(P\vec{u}) + P(\nabla\cdot\vec{V})$，右邊此項代表動能與熱能轉換的壓力膨脹項。)
5. 右式第 4 項 $-\phi$：
   黏滯耗散率（Viscous Dissipation Rate）。黏滯應力將宏觀流動動能做功轉換為微觀分子熱運動，符號為負代表動能的永久損失。

---

二、方程式 (2)：熱能（內能）方程式之物理意義

$$\rho\frac{de}{dt} = -\nabla\cdot\vec{q} - P(\nabla\cdot\vec{V}) + \phi$$

此式描述流體微元之熱能（比內能 $e$）隨體積變化率，各項意義如下：

1. 左式 $\rho\frac{de}{dt}$：
   單位體積流體比內能（熱能）的實質時變率（Rate of change of internal

核心觀念

1. 平行流（Parallel Flow）定義：流場中流體質點的流線（Streamlines）皆互相平行且為直線。若選取笛卡兒座標系，使 $x$ 軸沿著流動方向，則速度向量僅有單一方向之速度分量（例如 $\vec{u} = (u, 0, 0)$）。
2. 連續方程式（質量守恆）與對流加速度消去：
   對於不可壓縮流體（或由單向流動幾何限制）：
   $\nabla \cdot \vec{u} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$
   當 $v = w = 0$ 時，必然有 $\frac{\partial u}{\partial x} = 0$，意味著流速 $u$ 沿流向不隨位置改變（充分發展流，Fully Developed Flow）。此特性使得實質導數（物質導數）中的非線性對流項（Convective Acceleration）完全消失：
   $(\vec{u} \cdot \nabla)\vec{u} = 0$
3. Navier-Stokes 方程式（動量守恆）線性化：
   非線性項消去後，原本難以求解的非線性偏微分方程退化為線性偏微分方程（或常微分方程），因而在數學上具有解析解（Exact Analytical Solution）。

---

解題方法

1. 平行流之數學定義與運動學分析

設速度向量 $\vec{u} = (u, v, w)$ 且流動平行於 $x$ 軸：
$v = 0, \quad w = 0, \quad \vec{u} = u(x, y, z, t)\hat{i}$

由連續方程式（不可壓縮流）：
$\nabla \cdot \vec{u} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = \frac{\partial u}{\partial x} + 0 + 0 = 0 \implies \frac{\partial u}{\partial x} = 0$
因此，速度分量 $u$ 僅為橫向座標及時間之函數：
$u = u(y, z, t)$

物質導數（Material Derivative）展開為：
$\frac{d\vec{u}}{dt} = \frac{\partial \vec{u}}{\partial t} + (\vec{u} \cdot \nabla)\vec{u}$
其中非線性對流加速度各項：
$(\vec{u} \cdot \nabla)\vec{u} = \left( u\frac{\partial}{\partial x} + v\frac{\partial}{\partial y} + w\frac{\partial}{\partial z} \right) (u\hat{i}) = u\frac{\partial u}{\partial x}\hat{i} = 0$
因此加速度項簡化為純局部加速度：
$\frac{d\vec{u}}{dt} = \frac{\partial u}{\partial t}\hat{i}$

此外，因 $\nabla \cdot \vec{u} = 0$，原方程式中之可壓縮黏滯耗散項亦自然消去：
$\frac{\mu}{3}\nabla(\nabla \cdot \vec{u}) = 0$
拉普拉斯算子項簡化為：
$\nabla^2\vec{u} = \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)\hat{i} = \left( \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)\hat{i}$

重力加速度取標準重力向量（設 $z$ 軸垂直向上）：
$\vec{g} = -g\hat{k} \quad \text{或一般形式 } \vec{g} = g_x\hat{i} + g_y\hat{j} + g_z\hat{k}$

---

2. 運動方程式之三個分量式

將上述條件代入題給之 Navier-Stokes 方程式：
$\frac{d\vec{u}}{dt} = -\frac{1}{\rho}\nabla P + \vec{g} + \mu\nabla^2\vec{u} + \frac{\mu}{3}\nabla(\nabla\cdot\vec{u})$
（註：若方程式中黏滯項係數以動力黏度 $\mu$ 表示，密度為 $\rho$，等號兩端單位維持每單位質量的力或加速度）

各空間方向之純量分量式分別為：

• $x$ 分量（流動方向）：
  $\frac{\partial u}{\partial t} = -\frac{1}{\rho}\frac{\partial P}{\partial x} + g_x + \mu\left( \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)$

• $y$ 分量（橫向方向）：
  $0 = -\frac{1}{\rho}\frac{\partial P}{\partial y} + g_y \implies \frac{\partial P}{\partial y} = \rho g_y$

• $z$ 分量（法向／垂直方向）：
  $0 = -\frac{1}{\rho}\frac{\partial P}{\partial z} + g_z \implies \frac{\partial P}{\partial z} = \rho g_z$

若流動為二維（例如兩無窮大平板間流動，與 $z$ 無關）且無沿流向的重力分量（$g_x = 0$，$z$ 朝上使 $g_z = -g$），則三個分量式退化為最常見的形式：
$\frac{\partial u}{\partial t} = -\frac{1}{\rho}\frac{\partial P}{\partial x} + \mu\frac{\partial^2 u}{\partial y^2}$
$\frac{\partial P}{\partial y} = 0$

核心觀念

本題評量流體力學中**連續方程式（質量守恆方程式）之不同座標觀點（歐拉觀點 vs. 拉格朗日觀點）的物理意義，以及物質導數（Material Derivative）與局部時變率（Local Time Derivative）**之間的關係：

1. 尤拉觀點（Eulerian specification / 空間固定控制體積）：
   $\frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\vec{V}) = 0$
   關注空間中固定位置（Fixed Control Volume）的物理量隨時間變化，以及通過控制面的質量淨通量。
2. 拉格朗日觀點（Lagrangian specification / 移動流體質點）：
   $\frac{d\rho}{dt} = -\rho\nabla\cdot\vec{V}$
   關注隨著流體運動的特定流體微元（Fluid Parcel / Material System），其體積變形（膨脹或壓縮）所導致的密度變化。
3. 物質導數關聯式：
   $\frac{d\rho}{dt} = \frac{\partial\rho}{\partial t} + \vec{V}\cdot\nabla\rho$
   其中 $\dfrac{\partial\rho}{\partial t}$ 為局部時變率（空間固定點），$\vec{V}\cdot\nabla\rho$ 為對流項（Convective term），$\dfrac{d\rho}{dt}$ 為沿質點路徑的全時變率。

---

解題方法與完整推導

第一部分：說明兩式所代表的質量保守意義（5 %）

1. 方程式一：$\dfrac{\partial\rho}{\partial t} + \nabla\cdot(\rho\vec{V}) = 0$

   • 觀點：歐拉觀點（Eulerian framework / 控制體積觀點）。
   • 物理推導：
     對空間中固定的任意控制體積 $V$（外法線向量為 $\vec{n}$，控制表面為 $A$），質量守恆定律要求「控制體積內流體質量的時變率」加上「通過控制表面淨流出的質量流率」必須為零：
     $\frac{d}{dt}\int_V \rho\, dV + \oint_A \rho(\vec{V}\cdot\vec{n})\, dA = 0$
     由於控制體積固定不隨時間改變，將時間微分移入積分號內變成偏微分，並利用散度定理（Divergence Theorem）將表面積分轉換為體積分：
     $\int_V \left[ \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\vec{V}) \right] dV = 0$
     因為控制體積 $V$ 為任意選取，被積函數處處為零，即得微分形式：
     $\frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\vec{V}) = 0$
   • 物理意義：
     在空間中任意固定的位置（點），該處流體密度隨時間的增加率（$\dfrac{\partial\rho}{\partial t}$），精確等於由周圍流入該點的質量淨通量率（即負的質量流散度 $-\nabla\cdot(\rho\vec{V})$）。質量既不憑空產生，也不憑空消失。

2. 方程式二：$\dfrac{d\rho}{dt} = -\rho\nabla\cdot\vec{V}$

   • 觀點：拉格朗日觀點（Lagrangian framework / 物質系統觀點）。
   • 物理推導：
     追蹤一個隨流體運動的特定微元（質點系統），其質量為 $\delta m = \rho\,\delta V = \text{常數}$。對時間取物質導數（全微分）：
     $\frac{d(\delta m)}{dt} = \frac{d(\rho\,\delta V)}{dt} = \delta V \frac{d\rho}{dt} + \rho \frac{d(\delta V)}{dt} = 0$
     兩邊同除以 $\delta V$：
     $\frac{d\rho}{dt} + \rho \left( \frac{1}{\delta V} \frac{d(\delta V)}{dt} \right) = 0$
     由運動學定義可知，單位體積的體積膨脹率（Volumetric strain rate）等於速度向量的散度：
     $\frac{1}{\delta V} \frac{d(\delta V)}{dt} = \nabla\cdot\vec{V}$
     代入整理即得：
     $\frac{d\rho}{dt} = -\rho\nabla\cdot\vec{V}$
   • 物理意義：
     跟隨特定流體質點運動時，該質點的密度時變率取決於流場的體積膨脹或壓縮。若流場呈發散流動（$\nabla\cdot\vec{V} > 0$），質點體積膨脹，導致密度隨時間降低（$\dfrac{d\rho}{dt} < 0$）；若流場呈聚合壓縮（$\nabla\cdot\vec{V} < 0$），質點體積縮小，密度隨時間上升。

---

核心觀念

流體運動學（Kinematics of Fluid Flow）中描述流體運動有兩種基本觀點：

1. 尤拉法（Eulerian specification）：以空間固定點（空間座標 $\vec{x}$ 與時間 $t$）為觀測視角，記錄流體質點流經該固定控制體或位置時的各項場變數（如速度場 $\vec{u}(\vec{x}, t)$、壓力場 $p(\vec{x}, t)$）。
2. 拉格朗日法（Lagrangian specification）：又譯「蘭格蘭基法」，以流體個別質點（Fluid particle）為追蹤對象，類似古典質點力學的觀點。透過「識別標籤（初始位置）」追蹤特定流體質點在隨後時間的空間軌跡與物理量變化。

---

解題方法

本題為流體運動學基礎名詞與數學表示式之定義說明題。在拉格朗日描述中，流場中的質點由其在基準時間 $t_0$（通常取 $t_0 = 0$）的初始空間位置向量標定。位置向量 $\vec{X}(\vec{r_0}, t) = \big(X_1(\vec{r_0}, t),\, X_2(\vec{r_0}, t),\, X_3(\vec{r_0}, t)\big)$ 代表的是特定質點在時間 $t$ 時於三維空間中的直角座標位置（軌跡方程式）。

逐項具體說明如下：

• $\vec{r_0}$（質點的識別標記／初始位置向量）：
  代表流體質點在初始時刻（通常為 $t = 0$）的空間位置向量，即 $\vec{r_0} = (a, b, c)$ 或 $(X_{1,0}, X_{2,0}, X_{3,0})$。在拉格朗日法中，$\vec{r_0}$ 作為質點的「特徵標籤（Labeling coordinates）」，用來唯一標記並識別某一特定的流體微元（Fluid particle）。

• $t$（時間變數）：
  代表自初始時刻起算的物理時間（Time），為自變數。

• $X_1$、$X_2$、$X_3$（特定質點於時刻 $t$ 的瞬時空間座標分量）：

核心觀念

本題旨在測驗流體力學與流體物理學之基石概念——連續體假設（Continuum Hypothesis）與其適用準則。核心考點包含：

1. 微觀離散與巨觀連續的橋樑：物質在微觀上是由不連續的分子與原子組成，且分子間存在間隙；但在巨觀尺度下，為避免追蹤天文數字般的粒子運動方程式，將流體視為連續分佈且無間隙的介質。
2. 代表性基本體積（Representative Elementary Volume, REV）：巨觀流體物理量（如密度、速度、壓力、溫度）是建立在「包含足夠多分子以消除統計漲落，但尺度又遠小於巨觀流動特徵長度」的微小體積極限上。
3. 連續體適用指標——克努森數（Knudsen Number, $Kn$）：分子平均自由徑（Mean free path, $\lambda$）與巨觀特徵長度（Characteristic length scale, $L$）的比值，為判定流體是否可視為連續體的無因次準則。

---

解題方法

本題為申論說明題，作答應從「物理定義與數學表示」切入說明連續體假設，接著透過「無因次參數（克努森數）」給出量化判定標準與大氣科學中的實際應用情境。

1. 什麼是連續體假設（The continuum hypothesis）？

• 物理意涵：
  連續體假設是一種唯象（phenomenological）的巨觀簡化模型。它假定流體由連續分佈於空間中的物質所構成，流體內部不存在任何真空隙縫或離散粒子結構。在此假設下，流場中的物理性質（如密度 $\rho$、速度向量 $\mathbf{V}$、壓力 $p$、溫度 $T$ 等）均可表示為空間座標 $\mathbf{x} = (x, y, z)$ 與時間 $t$ 的連續且處處可微之點函數（continuous and differentiable point functions），從而得以採用微積分與歐拉（Eulerian）場論建立微分形式的控制方程式（如 Navier-Stokes 方程式）。

• 數學定義與代表性體積（REV）：
  以流體在某空間點的「密度 $\rho$」為例，若在該點周圍取一包含質量 $\delta m$ 的微小體積 $\delta V$，其密度定義為：

  $$\rho = \lim_{\delta V \to \delta V^*} \frac{\delta m}{\delta V}$$

  式中 $\delta V^*$ 即為代表性基本體積（Representative Elementary Volume, REV）：

  • 若選取的體積過小（$\delta V < \delta V^*$），體積內的分子數過少，個別分子的進出與熱運動會引起劇烈的分子統計漲落（statistical fluctuations），導致平均比值不穩定。
  • 當體積選在 $\delta V \approx \delta V^*$ 時，體積內含有數以億計的分子，統計平均值趨於穩定，且此微小體積的尺度仍遠小於整個流場的幾何特徵尺寸。連續體假設即建立在此穩定極限的存在之上。

---

2. 在什麼條件下，我們可以假定流體為連續體？

假定流體為連續體的物理條件為：流體分子的平均碰撞距離遠小於流場的巨觀特徵尺度，使得分子碰撞足夠頻繁，足以在極小時間與空間內建立局部熱力平衡（local thermodynamic equilibrium）。

定量上，採用無因次參數**克努森數（Knudsen number, $Kn$）**做為評判準據：

$$Kn = \frac{\lambda}{L}$$

• $\lambda$ 為分子的平均自由徑（Mean free path），即分子在兩次連續碰撞之間所移動的平均距離。
• $L$ 為流場的巨觀特徵長度（Characteristic macroscopic length scale），可為幾何尺寸（如固體表面直徑、管道口徑）、流動變化尺度（如邊界層厚度 $\delta$、激波厚度），或大氣物理中的垂直大氣標高（Scale height）與擾動波長。

適用條件分區判定準則：

1. 連續體區域（Continuum regime）：$Kn < 0.01$