108 年 國立臺灣大學生物機電系碩士班丁組《單元操作與輸送現象》
第 1 題20 分
Consider the fluidic lab-on-a-chip dye laser system. The liquid used in the laser is ethanol with a viscosity at room temperature of mPa·s. The dimensions of the rectangular channel are length mm, width µm, and height µm. For proper functioning the flow rate in the channel must be µL/h. Derive and calculate the pressure needed to run this device properly.
登入後即可作答並保存紀錄。
核心觀念
本題考查矩形微流道中的:
- 穩態、不可壓縮、層流流動
- 壓力驅動流
- Navier–Stokes 方程在微流道中的簡化
- 矩形流道的流量—壓降關係
由於流道寬度遠大於高度:
可將流動近似為兩平行板間的流動,主要速度梯度位於高度方向。此時流量與壓降的關係為
其中 為流道兩端的壓力差。
解題方法與推導
取流道長度方向為 ,高度方向為 ,並假設流動充分發展,因此速度僅隨 改變,即 。
穩態層流的動量方程可簡化為
由於壓力沿流道均勻下降,
因此
兩次積分後,配合壁面無滑移條件
得到速度分布
流量為速度對截面的積分:
代入速度分布:
整理得所需壓力差:
單位換算
黏度:
流道長度:
流道寬度:
流道高度:
流量:
壓力差計算
代入公式:
因為
第 2 題10 分
Consider a pressure-driven flow in a straight channel of length with a circular cross-section of radius given the Navier boundary condition with a non-zero slip length at the channel wall. This so-called Navier boundary condition reads
where is the slip length or Navier length defined as the distance behind the boundary where the tangent of the velocity field intersects the x axis. Note that for we recover the usual no-slip boundary condition.
(a) Derive the expressions for the slip length dependent velocity profile and for the hydraulic resistance . (10%)
(b) Derive the expression for the hydraulic resistance based on the velocity profile derived above. (10%)
(c) Calculate the relative change in the hydraulic resistance when changing from 0 nm to 50 nm for the two values of the radius µm and µm. (10%)
登入後即可作答並保存紀錄。
核心觀念
本題考查圓管內穩態、不可壓縮、充分發展之壓力驅動流,以及 Navier 滑移邊界條件。
主要公式如下:
-
壓力梯度提供流動驅動力:
-
圓柱座標下的軸向動量方程式:
-
體積流率與液壓阻力:
題目中的邊界條件寫成 ,但圓管中心線 不可能是管壁,因此應解讀為管壁 的 Navier 邊界條件。以徑向座標 向外增加表示,正確形式為
負號是因為速度由中心向管壁遞減,故 。
解題方法
先由圓管內的軸向動量方程式求得速度分布,再利用速度分布對截面積分求體積流率 ,最後由 得到液壓阻力。
這種方法可直接呈現滑移長度 如何改變速度剖面與流量。
(a) 速度分布與液壓阻力
1. 求速度分布
由
並令 ,可得
兩邊乘以 :
積分一次:
為使 處速度梯度保持有限,必須有 ,因此
再積分一次:
在管壁使用 Navier 滑移條件:
由
可得
另一方面,
故
因此速度分布為
代入 :
當 時,
即為通常的 Hagen–Poiseuille 拋物線速度分布。
管壁速度不再為零,而是
中心線最大速度為
2. 求液壓阻力
先求體積流率:
代入速度分布:
整理得
各項積分為
因此
又因為 ,
所以液壓阻力為
亦可寫成
(b) 由速度分布推導液壓阻力
由上題速度分布,
體積流率為
對圓形截面 ,因此
第 3 題20 分
Saturated steam at 130°C is flowing inside a ³⁄-in steel pipe having an ID of 0.824 in and an OD 1.05 in. The pipe is insulated with 1.5 in of insulation on the outside. The convective coefficient for the inside steam surface of the pipe is estimated as W/m²·K, and the convective coefficient on the outside of the lagging is estimated as W/m²·K. The mean thermal conductivity of the metal is W/m·K and W/m·K for the insulation. Calculate the heat loss for 1 m of pipe using resistance if the surrounding air is at 26 °C.
登入後即可作答並保存紀錄。
本題考驗複合壁結構的熱傳導問題,特別是圓柱形結構的熱阻計算。題目要求計算 1 米長管道的總熱損失,需考慮內部對流、鋼管導熱、保溫層導熱以及外部對流。
核心觀念:
- 熱阻概念 (Thermal Resistance): 在穩態熱傳導中,熱阻類似於電阻,可用來計算通過不同介質的熱流。總熱阻是各個熱阻的串聯總和。
- 圓柱座標下的熱傳導: 對於圓柱形壁面,熱傳導方程的解涉及對數函數。
- 對流熱傳阻:
- 導熱熱阻(圓柱):
解題步驟:
-
轉換單位並確定幾何參數:
所有尺寸都應轉換為國際單位制 (SI units)。- 內部蒸汽溫度 C
- 外部環境溫度 C
- 管道長度 m
- 鋼管內徑 in
- 鋼管外徑 in
- 保溫層外徑 in (注意 1.5 in 是保溫層的厚度)
轉換英吋 (in) 為公尺 (m):1 in = 0.0254 m
- 鋼管內徑 m m 內半徑 m
- 鋼管外徑 m m 鋼管外半徑 m
- 保溫層外徑 in
- 保溫層外徑 m m 保溫層外半徑 m
熱傳係數和導熱係數:
- 內部對流係數 W/m²·K
- 外部對流係數 W/m²·K
- 鋼管導熱係數 W/m·K
- 保溫層導熱係數 W/m·K
-
計算各部分的熱阻:
熱量從內部蒸汽流向外部空氣,經過三個熱傳阻串聯:內部對流熱阻、鋼管導熱熱阻、保溫層導熱熱阻、外部對流熱阻。-
內部對流熱阻 :
作用在內表面積 上。
m²
K/W -
鋼管導熱熱阻 :
作用在鋼管壁上,半徑從 到 。
K/W -
保溫層導熱熱阻 :
作用在保溫層上,半徑從 到 。
-
第 4 題10 分
Consider a microchannel perfusion bioreactor system within which the oxygen concentration can be controlled by adjusting the flow rate of media through the device. Nutrient media flows through the bioreactor with an average velocity equal to . On the lower surface, there is a monolayer of attached cells that consume the oxygen in the flowing stream.
The entering oxygen concentration is . Assuming steady-state, please
(a) derive the shell balance and the resulting differential equation for oxygen in a shell volume of the fluid with width from to and from to . (10%)
(b) Specify the boundary conditions and solve the differential equation by separation of variables. (20%)
登入後即可作答並保存紀錄。
核心觀念
本題考查微流道內氧氣的穩態傳遞,主要使用:
- 對流傳遞:
- Fick 第一定律:
- 穩態殼層衡算
- 對流—擴散方程式
- 分離變數法與本徵函數展開
設定微流道長度方向為 ,垂直方向為 。細胞附著於下表面 ,流體上表面為 ,流體平均速度為 ,氧氣擴散係數為 。細胞表面氧氣消耗通量以一階形式表示為
其中 為細胞層的表面消耗傳質係數。題目未提供 、 與細胞消耗速率,因此以下採用此標準假設。
(a) 殼層衡算與微分方程式
取流體內一微小殼層,其範圍為
流道垂直紙面方向寬度為 。
在穩態、無體積反應的情況下,殼層衡算為:
方向的對流傳遞
流體沿 方向流動,因此對流流入與流出量分別為
其差值為
方向的擴散傳遞
由 Fick 第一定律, 方向擴散通量為
因此兩側擴散通量的淨流入量為
方向的擴散傳遞
方向擴散通量為
因此上下表面的淨流入量為
完整殼層衡算
合併各項:
利用 Taylor 展開並除以 ,取 ,得到
即穩態二維對流—擴散方程式:
微流道通常具有高 Peclet 數,軸向擴散相較於軸向對流可忽略,因此採用常見的邊界層近似:
細胞表面的消耗邊界條件
在下表面 ,氧氣由流體向下擴散至細胞。由於 軸正方向向上,流向細胞的氧氣通量大小為 ,故有
在上表面 ,若上方為不可透過氧氣的壁面,則
(b) 邊界條件與分離變數解
邊界條件
入口氧氣濃度均勻為 :
細胞表面:
上壁面無氧氣通量:
因此要求解的方程式為
分離變數
令
代入控制方程式:
整理為
其中 為分離常數。於是得到兩個常微分方程式:
以及