第二章 文獻回顧
第一節 相關研究
風洞依照其送風方式可以分為吹送式及吸入式兩種,而本所風洞設 計之初即規劃為兩用式可滿足不同試驗需求。本所風洞在驗收時有針對 相關風場性能進行量測驗證是否合於標準,所得結果於所內92 年所內研 究報告中有相關資料呈現。
目前國內外著名風洞會對外宣告其測試區尺寸、風速範圍及紊流強 度等項目,本研究將其規格作一表格化(如表 2-1 所示)整理以利比較,由 表上可以知道所內風洞相較於其他大型風洞其測試區或測試風速範圍都 是相當優越,目前欠缺為性能之認證,本研究便以此為出發點將進行相 關TAF 認證規劃。
3
表 2-1 國內外低速風洞設備整理
Colorado State University,USA
M
第二章 文獻回顧 University of
Western
Meteorological Wind Tunnel, Environmental Protection
Agency, EPA, U.S. 2.13.718.3 1.5~8 0.5 邊界層可達2m Industrial Aerodynamics WT,
Department of Engineering Science, Oxford University
4×2×14 0~30 < 0.6 建築及工業用途
New Boundary Layer WT, Danish Maritime Institute,
Denmark 2.6×1.8×20.8 0~26 < 0.5 NUS-HDB WT, Nation
University of Singapore 2.851.8~2.319 0~15 < 1.0 模擬大氣邊界層
上海同濟大學土
表 2-1 國內外低速風洞設備整理(續) Boundary Wind Tunnel,
Building Research Establishment, Gt.
Britain
218 0~20 0.6 建築物空氣 動力學
New large WT, Politecnico di Milano, University, Milan, Italy
[15]
13.8×3.8×30 0~15 < 2 氣動力分析
4×3.8×5 0~50 < 0.2 風工程 研究 (資料來源:本研究整理)
第二章 文獻回顧
橫向流動角度變化spanwise variation of flow angledM
d x l
馬赫數梯度 Mach gradient '
nF n
干擾頻率含量 disturbance frequency contentT
T
D 露點溫度餘度 dew point temperature marginp
'密度擾動以及
p
'的空間相關性Density fluctuation and integral spatial correlation of
p
'w
:向上流速度分量、U
:自由流速度、c
:翼弦長、x
:流動方向座標、l
:平均自由徑、u '
:速度擾動分量、'
p
:靜壓擾動量的方均根值、q
:自由流動壓、n
:縮減頻率 (資料來源:本研究整理)TAF 認證是一個國家級實驗室永續經營重要的環節,並建立實驗室 儀器設備相關校正記錄、標訂作業程序及人員訓練等相關項目,其中試 驗數據誤差分析可因為長時間不間斷紀錄,藉以掌握實驗室試驗設備目 前狀況,可以維持實驗室一定水準試驗準確性及客戶對實驗室之信賴,
以下就量測不確定度及項目規劃進行說明:
(1)量測不確定度
試驗在分析數據上都會考量不確定度因素,是因為量測結果經常不是 一個固定的數據,也就是說量測的數據會有一定程度的變化量(有時候會 增加,有時候也會減少),如果經過長期性的量測就可以發現有這些變化量 的存在,因此我們將量測數據區分為兩種分量來表示:量測數據 =基準量
+變化量。量測不確定度是表示不能確定的程度,如果量測數據沒有含不 確定度時,則該數據就不能夠充分表達它的真實性,所以在國際上各個國 家的量測標準實驗室之校正或檢測報告, 都應具有量測不確定度值才能 符合科學與工商業界的要求。
在1980年國際度量衡局(BIPM)召集了國際上各個國家的量測標準實 驗室之專家,共同研討提出一種評估量測不確定度的建議大綱,1981年10 月國際度量衡委員會( CIPM )正式採用。
傳統方式的量測不確定度評估是將量測不確定度的成份區分為隨機 不 確 定 度 (random uncertainty , 簡 稱 Ur) 和 系 統 不 確 定 度 ( system uncertainty, 簡稱 Us),
然後再以適當的方式加以合併而得到量測系統的總不確定度。
第二章 文獻回顧
2
第二章 文獻回顧
第三節 有限高圓柱之尾流分析
Cheol-Woo [31] 等人以熱線測速儀量測有限高圓柱之尾流速度,以 及利用流場可視化描繪出圓柱尾流的狀態。實驗的規格與細節如表 2-3 所示。
表 2-3 熱線測速儀量測實驗規格
風洞測試段尺寸 (m)
長 6
寬 0.72 高 0.6 自由流速度 (m/s) 10
自由流紊流強度 (%) < 0.08 邊界層厚度 (mm) 4.1
圓柱直徑 (mm) 30
雷諾數 20000
展弦比 5、6、10、13、17.3(2D) (資料來源: 本研究整理)
由2-2 及 2-3 可看到有限高圓柱尾流支流場可視化結果。而由圖三可 看到紊流強度隨著下游距離變化也會和展弦比有關。
圖 2-2 尾流流場可視化之側視圖。
展弦比為5。
(資料來源: Cheol-Woo [31])
圖 2-3 尾流流場可視化之後視圖。
左邊流場截面位於圓柱後方X/D=0,右邊流場截面位於圓柱後方 X/D=0.6。
(資料來源: Cheol-Woo [31])
第二章 文獻回顧
圖 2-4 紊流強度隨下游距離變化示意圖。
量測高度為Z/L=0.5。
(資料來源: Cheol-Woo [31])
另外Cheol-Woo [32] 等人在不同邊界層厚度下以同一模型進行尾流 紊流強度量測以及流場可視化實驗,邊界層速度剖面圖則以大氣邊界層 指數律作為指標:
0
U z z
nU h
其中
U
0是自由流速度,h
為圓柱高度。實驗中則使用兩種不同的n,以n=0.14 來模擬平原上方的大氣邊界層,n=0.22 則模擬郊區上方的大氣 邊界層。圖2-5 顯示這次實驗所使用的自由流速度剖面圖及紊流強度剖面 圖情況。
圖 2-5 左邊為不同狀況下的自由流速度剖面圖,右邊為不同狀況下 的紊流強度剖面圖。
(資料來源: Cheol-Woo [32])
由圖2-6 可以看到在兩種不同的情況下,圓柱尾流流場可視化的側視 圖。圖六則為在n=0.14 之大氣邊界層下從不同下游距離處觀察圓柱尾流 流場可視化。圖 2-7 則為在 n=0.14 之大氣邊界層下不同高度流場可視化 的俯視圖。由圖 2-8 可以看到不同情況下的速度剖面圖和紊流強度剖面 圖,圖2-9 則為不同情況下紊流強度分布的示意圖。
第二章 文獻回顧
圖 2-6 有限高圓柱的尾流流場可視化側視圖。
展弦比為6,Y/D=0。
(資料來源: Cheol-Woo [32])
圖 2-7 有限高圓柱不同下游距離的尾流流場可視化之後視圖。
n=0.14
(資料來源: Cheol-Woo [32])
圖 2-8 有限高圓柱在不同高度的流場可視化俯視圖。
n=0.14
(資料來源: Cheol-Woo [32])
圖 2-9 在不同情況下,下游距離 X/D=5 時的速度剖面圖及紊流強度 剖面圖。
量測截面高度Z/D=0.5。
(資料來源: Cheol-Woo [32])
第二章 文獻回顧
圖 2-10 不同情況下的紊流強度分布圖。
量測截面高度Z/D=0.5。
(資料來源: Cheol-Woo [32])
M. S. Adaramola [33] 等人以二維熱膜測速儀量測有限高圓柱之尾流 速度,以及尾流的紊流強度分布。實驗的規格與細節如表2-4 所示。
表 2-4 二維熱膜測速儀量測實驗規格
風洞測試段尺寸 (m)
長 1.96
寬 1.13
高 0.91
自由流速度 (m/s) 30
自由流紊流強度 (%) < 0.5
邊界層厚度 (mm) 94.5
圓柱直徑 (mm) 31.5
雷諾數 60000
展弦比 3、5、7、9
(資料來源: 本研究整理)
圖 2-11、2-12 呈現不同展弦比情況下尾流的流線、垂直方向速度輪 廓圖,圖 2-13、2-14 則呈現不同展弦比情況下尾流的流線、垂直方向紊 流強度的輪廓圖。圖2-15 則為位於不同下游截面,不同展弦比情況下的 紊流強度、雷諾應力分布圖。
圖 2-11 不同展弦比有限高圓柱尾流之流線方向速度輪廓圖。
量測截面在下游距離X/D=6。
(資料來源: M. S. Adaramola [33])
圖 2-12 不同展弦比有限高圓柱尾流之垂直方向速度輪廓圖。
量測截面在下游距離X/D=6,虛線代表負(向下)速度。
(資料來源: M. S. Adaramola [33])
第二章 文獻回顧
圖 2-13 不同展弦比有限高圓柱尾流之流線方向紊流強度分布圖。
量測截面在下游距離X/D=6。
(資料來源: M. S. Adaramola [33])
圖 2-14 不同展弦比有限高圓柱尾流之垂直方向紊流強度分布圖。
量測截面在下游距離X/D=6。
(資料來源: M. S. Adaramola [33])
圖 2-15 不同下游處、不同展弦比的尾流輪廓圖。
(a)圖為流線方向紊流強度,(b)圖為垂直方向紊流強度,(c)圖為雷諾 應力。
(資料來源: M. S. Adaramola [33])
D. Sumner [34] 等人以一維熱線測速儀、一維熱膜測速儀量測有限高 圓柱之尾流速度,且以六力平衡儀測量阻力。實驗的規格與細節如表2-5 所示。
第二章 文獻回顧
表 2-5 以六力平衡儀測量阻力實驗規格
風洞測試段尺寸 (m)
長 1.96 寬 1.13 高 0.91 自由流速度 (m/s) 30
自由流紊流強度 (%) < 0.5
邊界層厚度 (mm) 82
圓柱直徑 (mm) 31.5
雷諾數 60000
展弦比 3、5、7、9
表 2-5 呈現不同展弦比、不同雷諾數以及不同邊界層對直徑比值情況 下的有限高圓柱的阻力係數。
(資料來源: 本研究整理)
表 2-
6 阻力係數比較雷諾數
D
展弦比備註 3 5 7 9
6 10
4 2.6 0.61 0.74 0.78 0.81 [34]7 10
4 - 0.63 0.67 [35]6 10
4 2.0 0.68 [36]3 10
4 3.3 0.68 [36]4 4
0.3 10 ~ 6 10
2.5 0.64 0.69 [37]1.3 10
4 0.15 0.75 0.76 0.87 1 [38]4 4
2.5 10 ~ 4.7 10
0.13 0.85 [39](資料來源: 本研究整理)
圖2-16 呈現不同展弦比有限高圓柱尾流中的流線方向渦度輪廓圖,
而圖2-17 呈現在不同下游距離尾流中流線方向渦度的發展。
圖 2-16 流線方向的無因次化渦度輪廓圖。
量測截面位於下游處X/D=6,實線為逆時針渦度,虛線為順時針渦度。
(資料來源: D. Sumner [34])
第二章 文獻回顧
圖 2-17 同一個展弦比(AR=9)、不同下游處的流線方向渦度輪廓圖。
(a)圖為 X/D=6,(b)圖為 X/D=8,(c)圖為 X/D=10,實線為逆時針渦 度,虛線為順時針渦度。
(資料來源: D. Sumner [34])
圖 2-18 同一個展弦比(AR=5)、不同下游處的流線方向渦度輪廓圖。
(a)圖為 X/D=4,(b)圖為 X/D6,(c)圖為 X/D=8,(d)圖為 X/D=10,
實線為逆時針渦度,虛線為順時針渦度。
(資料來源: D. Sumner [34])
第二章 文獻回顧
第四節 圓柱流場
Bearman [40]量測了許多雷諾數情況下的圓柱的表面壓力,以及使用 熱線測速儀量測頻率,呈現頻率與雷諾數的關係。圓柱直徑為0.178 公尺 而長度為2.13 公尺,垂直立於風洞。如圖 2-19 呈現 Z 方向的基部壓力係 數變化,可以發現隨著風速越大、雷諾數越高,基部壓力受到邊界的影 響越大。圖2-20 則呈現基部壓力係數對雷諾數的關係,基部壓力係數在 臨界雷諾數(
R ~ 10
5)這個範圍內變化相當劇烈。圖 2-21 為無因次化頻率 隨著雷諾數的變化,可發現頻率在某一小段雷諾數內有相當大的變化,而在
R 1 10
5附近沒有太大的變化。圖 2-19 沿著 Z 方向的基部壓力係數變化。
Z 為從圓柱根部算起的距離,L 為圓柱長度,
代表R 2 10
5,代 表R 3.68 10
5, 代表R 4 10
5。(資料來源: Bearman [40] )
圖 2-20 基部壓力係數對雷諾數的變化
(資料來源: Bearman [40])圖 2-21 無因次化頻率(St)隨著雷諾數的變化
(資料來源: Bearman [40])第三章 研究方法
第三章 研究方法
第一節 試驗設施與量測設備
本研究之風洞實驗在本所台南歸仁風洞實驗室進行,該風洞本體為 垂直向的封閉式循環風洞系統,需要時可轉為開放式,總長度為77.9m,
最大寬度為9.12m,最大高度為 15.9m。本風洞有兩個測試區,第一測試 區截面積寬4 米,高度為 2.6 米,內部有可調式上蓋板,可依據所需不同 邊界層厚度變化進行調整,最高可以調整至 3 米高。第一測試區內設置 有 2 個迴轉盤,可進行不同風攻角模型之氣動力試驗,其角度旋轉最小 精度可達1 度。第二測試區截面積寬度為 6 米,高度為 2.6 米,如同第一 測試區設置有可調式上蓋板,可依據所需不同邊界層厚度變化進行調 整,最高可以調整至 3 米高,本試驗區域最主要可進行橋樑或其他大跨 度模型之氣動力試驗,測試區內設置有一迴轉盤可進行模型風向角之調 整。
圖 3-1 本所風洞設計圖
(資料來源:內政部建築研究所研究報告,2003 年)
本研究團隊將針對第一測試區基礎相關風場資料進行量測比對,所 選擇量測位置將由該測試區之入口第一迴轉盤中心,主要進行風速及紊 流強度分布資料擷取,本研究風洞實驗所使用之量測儀器介紹如下:
(1) 皮托管
本實驗中採用皮托管進行來流平均風速之量測,由皮托管所量測到 的壓力差值,利用伯努利方程式(Bernoulli equation),即依據後式計算出 相應之風速。
air
U p
2
(3-1) 同時將所得之風速為參考風速進行無因次化雷諾數的計算,
UD
D
Re
(3-2)
第三章 研究方法
圖 3-
2 皮托管(資料來源:本研究整理)
(2) 壓力轉換器
入口風速量測是利用皮拖管經由壓力管將壓力差傳遞至儀控室,然 後透過薄膜式壓力轉換器(pressure transducer,VALIDYNE DP103-18,參 見圖3-5)將壓力資料以數位形式輸出。量取兩端壓力管輸入之壓差,以 類比電壓訊號輸出,在電磁感應線圈造成不同電感,利用電橋變化出不 同電壓,搭配訊號放大器CD-23 調整電壓訊號輸出,將輸出電壓設定在 -10V~+10V 之間。壓力轉換器若與皮托管(pitot tube)連接,經率定後可 以量測流場平均速度。
薄膜式壓力轉換器率定應配合壓力轉換器內部的壓電膜片的受壓範
圍,依照其膜片可承受範圍,利用壓力校正器(DPI 610)連接兩條短油管 傳輸壓力給薄膜式轉換器之動壓與靜壓。壓力由小至大,直到可承受之 最大壓力,透過資料擷取系統將所測之電壓值轉換存檔後,其迴歸率定 曲線呈線性型態。
圖 3-3 薄膜式壓力轉換計
(資料來源:本研究整理)
圖 3-4 Validyne 壓力計訊號放大器 CD-23
(資料來源:本研究整理)
第三章 研究方法 (3) 資料擷取系統
實驗所量得之類比訊號係經由NI CompacDAQ-9172 擷取後作類比數 位(analog-digital)轉換。本系統最高可連結 8 個模組,本系統最高可連結 八個模組,時間解析度(Timing resolution)為 50ns,取樣頻率(Sampling rate) 最高可達3.2MHz。配合四個 NI 9215 模組,每個模組共有四個訊號輸入
實驗所量得之類比訊號係經由NI CompacDAQ-9172 擷取後作類比數 位(analog-digital)轉換。本系統最高可連結 8 個模組,本系統最高可連結 八個模組,時間解析度(Timing resolution)為 50ns,取樣頻率(Sampling rate) 最高可達3.2MHz。配合四個 NI 9215 模組,每個模組共有四個訊號輸入