第八章 煙流擴散
8.2 理論背景
8.2.3 大氣邊界層內擴散之物理模擬試驗
8.2.3.1 對模擬參數之考慮
在邊界層風洞內，應用縮尺模型進行物理模擬實驗，一般言之，
即使是使用大型之邊界層風洞，前二節所述之五個無因次參數亦無法 同時都能得到滿足。所幸，視問題的不同，在某些情況之下，可以放 鬆一些參數的相似性要求。
(1)羅斯培數， ，代表科氏力(Coriolis force)對大氣運動的影響 程度，在風洞實驗中對科氏力模擬有極大的困難，無法達到在原型與 限制，Halitsky[18]指出在方形模型周遭之流場，當雷諾數達 11,000 時，所量測到的污染濃度值已不再隨雷諾數而變化。美國環保局 比之乘積。普朗特數為動量擴散(momentum diffusivity)與熱擴散
(thermal diffusivity)之比，休密特數則為動量擴散與質量擴散(mass diffusivity)之比。可見普朗特數與休密特數皆為流體本身之性質而非 流場之性質。如果以空氣為介質進行擴散之物理模擬實驗，在原型與 模型之間其普朗特數及休密特數（對幾乎任何釋入空氣中之污染氣 體）幾乎相等。所以在衡量培克烈數與雷諾-休密特數是否相當時，
重點不在普朗特數與休密特數，而在雷諾數；如果一個流場之雷諾數 足夠高，則其間污染物質之傳輸主要是由紊流中之大型結構所帶動，
分子傳輸部份之貢獻極為微小。所以美國環保署之模擬準則[16]中指 出，在原型與模型間培克烈數或雷諾-休密特數相等之要求並不重 要，只要雷諾數足夠高，流場顯示出其對雷諾數之獨立性即可。
(1)福祿數， ，為慣性力與浮升力比值之平方根，與李查遜數 (Richardson number)之平方根成倒數關係。在模擬大氣擴散的實驗中 福祿數應該是一個最重要的參數。通常有二個福祿數必須考慮，即大 氣氣流之福祿數以及煙囪排氣之福祿數。在實驗中為達到與原型相似 的較小福祿數值，必須將風洞之風速調低，但為了達到模型與原型間 雷諾數相似的要求又必須儘可能將風速提高。此一互相矛盾兩極化的 要求，其解決方法是在滿足福祿數相似之要求下，儘可能提高風洞的 運轉風速以保證流場特性的雷諾數獨立性。
Fr
8.2.3.2 對煙囪廢氣昇流之考慮
如前所述模擬邊界層內之擴散現象，必須對煙流之上升加以考 慮。有關對模擬廢氣煙流上昇之討論可於 Strom & Halitsky[19]，
Cermak[15]，Isyumov & Tanaka[20]，以及 Snyder[16]等之論文中見到。
廢氣昇流之相關變數有：
Hs：煙囪高度
對由動量效應(momentum effect)及浮昇力效應(buoyancy effect) 所主宰之煙流，其相關的模擬參數如下：
式中 及 分別為煙氣昇流現象之動力尺度(momentum length) 與浮力尺度(buoyancy length)。煙流上昇之弧形軌跡及高度可由l 與l 表示。在模擬煙囪廢氣昇流之實驗中，在原型與模型間如何保持(5.9)
(c)W ， 以及W 三參數在原型與模型間必須 相似。
s /U ρ /s ρa s/(gD∆ρ/ρa)1/2
(2) 無煙氣下沖之問題 — 在地形地物等之氣動力影響之下之煙流擴 散模擬相似要求。
(a)煙氣的雷諾數儘可能提高，最好能超過 15,000，如果必須將煙 氣之雷諾數減小至小於 2,300，則可能需將煙囪內壁加粗以保 捲增常數(entrainment coefficients)， ， 。
h β1 β2
lm及 可分別視為因為動量及浮升力所造成煙流曲線的初始半 (Snyder,1981)，提出如下之建議：
(1) 保証煙氣釋流之紊流特性
本實驗風速量測所使用之儀器包括有：
1.皮托管(pitot-static tube)：用於風洞內風速之監控（如圖八）。
2.煙線(smoke wire)：用於風速在 1.2 米/秒以下風速之量測與監控。
3.恒溫式熱膜風速儀(Constant-temperature hot-film anemometer)：用 於風速剖面及紊流特性茲之量測(包括軸向及垂向之平均流速U 及
，以及紊流強度
W u′2 U及 w′2 U)。整個儀器系統包括:
(i) IFA-300 智慧型風速儀(anemometer)，
(ii) TSI Model 1243 邊界層〝 X 〞型探針(〝 X 〞probe)，
(iii) TSI Model 1125 探針率定器（probe calibration）。〝 X 〞型探 針之率定是根據製造廠商TSI 所設定之標準程序。
由 風 速 儀 所 輸 出 之 電 子 訊 號 經 過 類 比/ 數 位 轉 換 器 （ A/D converter）之二個頻道，由電腦系統之讀入進行計算、分析及紀錄存 檔。
8.3.2 煙氣之模擬
初步考慮只用中性浮力之煙流以抬升源之方式釋入邊界層內。模 擬之煙氣將以 1﹪之甲烷與空氣混和。
8.3.3 釋放系統
示蹤氣體（tracer）包括由甲烷、空氣所組合而成之混合氣體。
混 合 氣 體 之 百 分 比 是 經 由 兩 架 單 管 之 浮 子 流 量 計 （single tube flowmeter）所控制。
非標準狀況及非空氣之流體其相等於標準狀況下之空氣流量可
以下式計算得到。
控制之中心數據採集、處理、分析及結果展示系統。
筆者自行設計製造之氣體採樣系統是一架先進之採樣系統，它可 同時採集 48 個樣本，所採樣本不為外界環境所污染，表現穩定，能 在相當寬之濃度範圍內採集到重複性高之樣本，對採樣器設有標準之 檢定程序，以確保採樣器之性能表現。
單管（single column）氣相層析儀包括有：
1.火焰離子化偵測器(flame-ionization detector) 2.自動溫度控制器(automatic temperature controler) 3.差分電表(differential electrometer)
4.流動系統（flow system）- 包括壓力調控及壓力表（pressure regulator and pressure gauge）
氣相層析儀使用前之率定包括有經由標準濃度之混合氣體，以及 實驗室所自行混和之不等濃度之氣體進行系統之線性 - 非線性反應 之率定，以及零點偏移量（zero shift）之確定。
所 能 量 測 到 之 濃 度 下 限 ， 由 儀 器 之 靈 敏 度 及 背 景 濃 度
（background concentration）所定。在整個濃度之測試過程中，背景 濃度之量測都包括在內，每一測點之濃度值都經由背景濃度加以修
8.4 煙流擴散實驗結果
進行擴散實驗之前先對流場作一初步的測試，其空洞水平風壓分 佈穩定，上下震盪的幅度約為百分之一（圖42），下游約 25 米處（第 二迴轉盤）所量得的邊界厚度約 30cm 厚（圖 43）。為模擬大氣邊界 條件，於測試區上游安置粗糙元（如圖40 及圖 41）以增加邊界層厚 度，經量測後邊界層厚度增加為 170cm 厚。在迴轉盤中心處設置一 模擬煙囪排放擴散氣體，於排放源下游 1m、2m、3m 及 3.5m 處採集 氣體，經由氣體分析儀可測得其濃度，圖45 及圖 46 分別為不同距離 下水平及垂直方向濃度分佈。
第九章 教育訓練
本計畫亦將協助完成協助內政部建研所培訓風洞操作與維護所 需人力，以落實培訓風洞操作與維護所需人力的目的；而教育訓練包 含了儀器設備的教育訓練及性能驗證實驗的教育訓練。
中鋼公司、仕大電機及豪頓公司三家設備廠商已於今年5 月中旬 風洞驗收時完成了相關的教育訓練，而量測儀器部分有皮托公司及鑫 群科技兩家廠商分別針對多孔式壓力掃描閥及示波器、數位電表與電 源供應器等相關儀器進行操作示範（如圖47）。
風洞基本性能測試實驗於 93 年 12 月 23 日進行相關教育訓練，
內容包含有皮托管與熱線測速儀介紹、基本性能實驗項目、基本性能 實驗結果及實際示範皮托管及熱線測速儀操作（如圖 48）。
第十章 結語
歷經三年時程，很高興能順利完成此風洞的籌建工作。此風洞的 硬體驗收在今年五月底完成，隨即進行一系列實驗測試項目，其目的 為驗證本風洞的性能。整體而言，所得的測試結果顯示風洞測試區流 場品質正在逐項驗證之中，另外完成了橋樑斷面測試項目，與煙流擴 散測試項目。由這些實驗結果進一步描繪出此風洞已具備大型建築物 模型測試能力相信假以時日，用心經營，這個實驗室能為我國的建築 界產官學研究做很多的工作，同時也可藉此實驗室的成果促進國際科 技交流。
另一個角度而言，由於本風洞剛籌建完成，其各項性能有待進一 步確認，而測試的技巧也有待發展建立，這些工作均需要人力與時間 的投入，期此風洞的建立能吸引更多人投入，所累積的成果，才是實 驗室最重要的資產。
在此特別感謝內政部建築研究所的同仁、成大航太中心及相關單 位同仁的協助，使得在籌建過程工作進行順利。
參考文獻
Prandtl L (1927), “U¨ ber den Reibungswiderstand stro¨mender Luft.” Ergebnisse AVA Go¨ttingen III: 1–5
Cermak, J. E., “Applications of fluid mechanics to wind engineering – a freeman scholar lecture,” Trans. of the ASME, J. of Fluids Eng., vol. 97, pp. 9-38, 1975.
Cermak, J. E., “Aerodynamics of building,” Annual Review of Fluid Mechanics,” vol. 8, pp. 75-106, 1976.
Houghton, E. L., & Carruther, N. B. Wind forces on building
and structures, John Willy & Sons, New York, 1976.
8. Obasaju, E. D., Ermshaus, R. and Naudasher, E. (1989), "Vortex induced stream wise oscillations of square-section cylinder in a uniform stream", Journal of Fluid Mech., Vol. 213, pp.171-189.
9. 鄭啟明，林堉溢，葉博文，”以斷面模型實驗探討斜張橋之顫 振特性”，結構工程，第十六卷，第四期，第 85~98 頁, 2001。
10. Cheng, C.M., Lin, Y.Y., ”The Effects of Deck Shape and
Oncoming Turbulence on the Bridge Aerodynamics” Presented at The Inaugural Meeting of Wind Engineering Research Center at Tokyo Institute of Polytechnics, Atsugi/Japan, 2001.
11. Simiu, E. and Scanlan, R. H. " Wind Effects on Structures ",
John Wiley & Sons.,N.Y.（1986）
12. 林堉溢，鄭啟明，蔡子文，1997，"斜張橋受風載重之顫振 及亂流效應分析"，結構工程，第十二卷，第一期，
pp.107-117。
13. Scanlan, R. H. and Tonko, J. J. "Airfoil and Bridge Deck Flutter Derivative", Journal of Engineering Mechanics Division, ASCE, Vol.97, pp.1717-1737（1971）
14. Kubo, Y., Miyazaki, M. and Kato, K. (1989), "Effect of end plate and blockage of structure members on drag force", Journal
of Wind Eng. and Industrial Aerodynamics, 32, pp.329-342.
15. Cermak, J.E., “Application of fluid mechanics to wind
engineering,” A Freeman Scholar Lecture. ASME J. of Fluids Engng., 97, no.1, pp.9-38, 1975.
16. Snyder,W.H., “Guideline for fluid modeling of atmospheric diffusion,” EPA-600/8-81-009, 1981.
17. Lumley, J.L. & Panofsky, H.A., “The structure of Atmospheric Turbulence,” Interscience, N.Y. , p.239, 1964.
18. Halitsky, J. “Validation of scaling procedures for wind tunnel model testing of diffusion near building,” Report No.TR-69-8, Geophysical Sciences Laboratory, New York University, 90 p, 1969.
19. Strom, G.H. & Halitsky, J., “Important considerations in the use of the wind tunnel for pollution studies of power plants,” Paper
No.54SA-41 ASME Semi-Annual Meeting, Pittsburgh, PA, 1954.
20. Isyumov, N. & Tanaka, H. “Wind tunnel modeling of stack gas dispersion difficulties and approximations,” Wind engineering (ed. By Cermak, J.E.), Pergamon Press, Oxford, 2 , pp.987-1001, 1979.
21. 張能復，「台北市內湖垃圾焚化廠廢氣排放對廠址附近環境 空氣品質的影響」，研究報告 424p, 1988.
22. Briggs, G.A., “Plume rise predictions,” ATDL No. 75/15, Atmos. Turb. and Diff. Lab., NOAA Environ Res. Lab., Oak Ridge, TN, 53p, 1975
23. Hoult, D. P. “Simulation of buoyant pollutants in the atmospheric boundary layer flow studies in air and water pollution,” ASME, N. Y., pp. 61-69, 1973.
附 圖
圖1. 93 年 5 月 26 日風洞驗收
圖2 風扇測試現況暨風扇試車檢討會議
薄膜式壓力傳送器 皮扥管
圖 3 皮扥管及薄膜式壓力傳送器圖
(a)一維量測之熱線探針
(b)二維量測之熱線探針
(c)三維量測之熱線探針 圖 4 不同型式之熱線探針示意圖
圖5 三孔的 cobra probe 圖
圖 6 風扇運轉時間與其溫度、第一測試區入口速度之關係圖
0
Runner speed (rpm)
V elo city ( m /s )
First test section Second test section
y = 0.0937x - 0.1418
y = 0.0585x - 0.1775
圖8 風扇轉速與速度之關係圖
圖 9 第一測試區下游 2.5m 之平均速度分布圖
圖10 第一測試區下游 2.5m 之二維平均速度剖面圖
圖11 第一測試區下游 15m 之平均速度分布圖
圖 12 第一測試區下游 15m 之二維平均速度剖面圖
圖 13 第一測試區下游 25.5m 之平均速度分布圖
圖14 第一測試區下游 25.5m 之平均速度剖面圖
圖 15 第一測試區下游 25.5m 之紊流強度分布圖
圖 16 第一測試區下游 25.5m 之紊流強度剖面圖
圖17 第一測試區下游 3m 邊界層厚度
圖18 第一測試區下游 15m 邊界層厚度