濕式排煙脫硫法(FGD)因脫硫效率高、系統穩定，以及成本相對低廉等原因被廣泛的運用於鋼鐵業。本研究針對煙氣脫硫進行模擬，FGD模擬主要要考慮的有氣、液兩相流之流場以及化學反應兩個部分。在流場部分，本研究採用計算流體力學(CFD)方法，先建構二維脫硫塔的CFD模型，並以Eulerian-Eulerian Model進行氣液雙相的模擬，縮短模擬時間發現問題並解嘗試各種方法解決，再進行網格數較大的三維流場模擬，此外，本研究有針對加入脫硫塔篩板進行詳細的模擬與討論；在化學反應的研究，則以文獻整理出FGD相關的化學反應路徑，再進行模擬化學物質隨時間的濃度變化測試。最後流場與化學反應做耦合模擬脫硫系統的運作，即在計算流體力學的守恆方程式中加入化學反應，模擬的結果再與實驗脫硫數據比較。

This research focuses on the 3D dynamic behaviors of a wet flue-gas desulfurization (WFGD) by numerical simulation. Our simulation utilizes the Eulerian-Eulerian model on multiphase flow and couple with chemical reaction. Then the gas and liquid hydrodynamics inside a Sieve Tray Tower has been simulated using a commercial computational fluid dynamics (CFD) code, while a model that treats the absorption process of SO2 has been developed and implemented in the software through dedicated modules. The simulation has been validated by absorption of SO2 measurements of the experimental data, and further analyses have been conducted through the volume of fraction and velocity contours. Our goal is to analyze the result of simulation on every sieve tray and capture more complicated flow structures.
According to the results, the dynamic behaviors of the flow field and liquid-gas ratio can affect the sulphur dioxide removal efficiency significantly. The result shows that chemical reaction couple with Eulerian-Eulerian model can help this research to have better understandings of the complexity of the FGD and prepare for the further optimization.

中文摘要 i
Abstract ii
目錄 iii
圖目錄 iv
表目錄 vi
符號說明 vii
第一章、緒論 1
1.1 脫硫塔背景之簡介 1
1.2 影響脫硫塔之重要因素 3
1.2.1 流場動態行為 3
1.2.2 篩板幾何 4
1.2.3 化學反應 5
1.3 研究目的 6
第二章、研究方法 8
2.1 模擬FGD 之數值方法與流程 8
2.1.1 流體統御方程 9
2.1.2 化學反應機制 12
2.1.3 流場幾何、邊界條件及網格設定 14
第三章、遇到之困難與討論 20
3.1 Mixture model 與Eulerian-Eulerian Model 比較 21
3.2 下邊界修正 21
3.3 Porous media 與篩板對流場的影響之比較 22
第四章、模擬結果 26
4.1 模擬與實驗之驗證 26
4.2 篩板的效果分析 27
4.3 化學反應與流場耦合 28
第五章、結論 37
未來工作 37
參考文獻 37

[1] Health and Environmental Impacts of SO2. 美國國家環保局. [2008-12-18]
[2] Biondo, S.J. and Marten, J.C., A History of Flue Gas Desulphurization Systems Since 1850, Journal of the Air Pollution Control Association, Vol. 27, No. 10, pp 948–961, October 1977.
[3] AS Mujumdar（2007）。工業烘乾手冊。CRC出版社。頁。710，國際標準書號 1-57444-668-1。
[4] Holleman, A. F.; Wiberg, E. Inorganic Chemistry. San Diego: Academic Press. 2001.
[5] http://viewonline.in/scrubbing.html 2014/5/14
[6] Japanese Advanced Environment Equipment http://www.gec.jp/JSIM_DATA/AIR/AIR_2/html/Doc_065.html
[7] Arun Kumar Sharma, D.S. N. Prasad, Shveta Acharya, and Rashmi Sharma. Utility and Application of FGD System (Flue Gas Desulphurization) in Chemical and Environmental Engineering. International Journal of Chemical Engineering and Applications, Vol. 3, No. 2, April 2012.
[8] A. STERGARŠEK1, M. GERBEC1, R. KOCJANČIČ1, and I. HRASTEL2. OPTIMISED INTEGRATED FGD PROCESS
[9] Juncong Sai†, Shaohua Wu, Rui Xu, Rui Sun, Yan Zhao and Yukun Qin
Korean J. Chem. Eng., 24(3), 481-488 (2007)
[10] 王福軍. 計算流體動力學分析-CFD軟體原理與應用. 北京: 清華大學出版社, 2005, 1
[11] 周力行. 湍流氣粒多相流動和燃燒的理論與數值模擬. 北京: 科學出版社, 1994, 153 – 200
[12] 歐陽潔, 李靜海. 類比氣固多相流動非均勻結構的顆粒運動分解軌道模型. 中國科學(B輯), 1999, 29(1): 29 – 38
[13] Griffiths W D, Boysan F. Computational fluid dynamics (CFD) and empirical modeling of the performance of a number of cyclone samplers. Journal of Aerosol Science, 1996, 27(2): 281 – 304
[14] Haim M, Weiss Y, Kalman H, et al. The effect of the inlet conditions on the numerical solutions of particle-gas flows. Advanced Powder Technology, 2003, 14(1): 87 – 110
[15] Zhang Z, Xie Z L. Numerical simulation of fluid-solid two-phase flows. Journal of Chemical Industry and Engineering, 2001, 52(1): 1 - 12
[16] Eric P, Bo L. Fundamentals of turbulent gas-solid flows applied to circulating fluidized bed combustion. Progress in Energy and Combustion Science, 1998, 24(4): 259 – 296
[17] Thibaut Neveux, Yann Le Moullec. Wet Industrial Flue Gas Desulfurization Unit Model Development and Validation on Industrial Data, Fluid Mechanics for Energy and Environment (MFEE) Department, Research and Devlopment, EDF, 6 Quai Watier, 78400 Chatou, France.
[18] N. Fueyo, A. Gomez, and J.F. Gonzalez. A Comprehensive Mathematical Model of Flue-gas Desulfurization.
[19] Luca M, Modeling of the fluid dynamics and SO2 absorption in a gas-liquid reactor. Chemical Engineering Journal 162(2010)217-226.
[20] Xiang G, Wang H, Z. Y. LUO, K.F. CEN. CFD simulation with enhancement factor of sulfur dioxide absorption in the spray scrubber. Journal of Zhejiang University SCIENCE A ISSN 1673-565X (Print); ISSN 1862-1775
[21] Luca Marocco, Fabio Inzoli .Multiphase Euler–Lagrange CFD simulation applied to Wet Flue Gas Desulphurisation technology. International Journal of Multiphase Flow 35 (2009) 185–194
[22] Fang Zeng, LianQing Yin and Li Chen. Numerical Simulation and Optimized Design of the Wet Flue Gas Desulphurization Spray Tower.
[23] 過小玲, 金保升, 孫志翱. 裝有氣流分布板的脫硫噴林塔流場數值模擬.
[24] Teleken, J. G；Werle, L. O. CFD simulation of multiphase flow in a sieve tray of a distillation column. Brazilian journal of petroleum and gas
[25] R. KRISHNA (FELLOW), J. M. VAN BATEN, J. ELLENBERGER , A. P. HIGLER and R. TAYLOR. CFD SIMULATIONS OF SIEVE TRAY HYDRODYNAMICS. Trans IChemE, Vol 77, Part A, October 1999
[26] H. Darcy, Les Fontaines Publiques de la Ville de Dijon, Dalmont, Paris (1856).
[27] Abbas Firoozabadi , L.K.Thomas, Bert Todd. High-Velocity Flow in Porous Media.
[28] 景启国, 徐康富, 马永亮, 童 玫. 環境汙染治理技術與設備. 第6卷第11期2005年11月
[29] 邱滔, 羅光輝, 陆咏凯. 環境科學與技術. 第34卷第9期2011年9月
[30] Reynolds, Osborne, 1895: "On the Dynamical Theory of Incompressible Viscous Fluids and the Determination of the Criterion." Philosophical Transactions of the Royal Society of London. A, v. 186, pp. 123-164.
[31] Eymard, R. Gallouët, T. R. Herbin, R. (2000) The finite volume method Handbook of Numerical Analysis, Vol. VII, 2000, p. 713–1020.
[32] ANSYS FLUENT Theory Guide. Release 14.0 (2011).
[33] Courant, Richard; Isaacson, E; Rees, M. (1952). "On the Solution of Nonlinear Hyperbolic Differential Equations by Finite Differences"
[34] Computer Methods in Applied Mechanics and Engineering 19 (1): 59–98,
[35] Fundamentals of FLUID MECHANICS. Sixth edition.
[36]on Kármán, Th. (1930), "Mechanische Ähnlichkeit und Turbulenz", Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Fachgruppe 1 (Mathematik) 5: 58–76 (also as: “Mechanical Similitude and Turbulence”, Tech. Mem. NACA, no. 611, 1931).
[37]Schiller, L., Naumann, Z., 1935. Z Ver. Deutsch. Ing. 77, 318.
[38] 鍾泰 南京理工大學學報 第23卷第1期1999年2月
[39] Solving the Cartesian Cut cell Interpolation Problem with a Tetrahedral Mesh.
B. Eng, Carleton University(1999)
[40] dropletsize search on web 2014/4/24
file:///C:/Users/Asus/Downloads/droplets_v6.pdf