本研究利用計算流體力學(Computational Fluid Dynamics, CFD)模擬探討薄膜接觸器應用於鹽水脫鹽與二氧化碳吸收之流力與熱質傳特性。薄膜接觸器使用多孔疏水性膜且模組型式包括空通道與具網狀間隔物通道。本研究使用FLUENT軟體之層流模式與VOF(Volume of Fluid)多相模式，薄膜為可滲透邊界條件。模式之驗證是利用文獻實驗數據且獲得相當接近之結果。
具間隔物模組之流力與熱質傳特性皆會受到間隔物配置之影響，而呈現上下震盪之變化，其壓降、熱傳係數與質傳係數均高於空通道。間隔物通道與空通道之濃度邊界層厚度分別約為1 mm與2-4 mm。使用具間隔物通道可大幅提升鹽水脫鹽系統與碳酸丙烯酯之二氧化碳物理性吸收系統之質傳通量。間隔物通道之水利損耗遠高於空通道，後者接近平板通道關聯式預測值。間隔物與空通道之熱傳係數與質傳係數均與文獻報導關聯式之預測值有很大偏差。

In this thesis, the fluid flow, heat and mass transfer characteristics in the membrane contactors for desalination and carbon dioxide physical absorption with propylene carbonate are studied using computational fluid dynamics (CFD) simulation. For desalination, the membrane modules are direct contact type. The membrane contactors simulated use hydrophobic membranes and either empty or spacer-filled channels. The laminar flow model and Volume of Fluid (VOF) multiphase model in FLUENT are employed with the membranes defined as permeable boundary conditions. The model simulation results are fairly close to the literature experimental data
The fluid flow, mass and heat transfer characteristics of spacer-filled channels show fluctuating patterns corresponding to the spacer configurations. The pressure drop, heat and mass transfer coefficients of spacer-filled channels are all higher than the empty channels. The concentration boundary layer thicknesses of spacer-filled and empty channels are approximately 1 mm and 2-4 mm, respectively. For both desalination and carbon dioxide absorption systems, the spacers significantly enhance the transmembrane mass fluxes. The hydraulic power consumptions of spacer-filled channels are much higher than the empty channels, which are close to the correlation for parallel-plate channels. The heat and mass transfer coefficients of both spacer-filled and empty channels are not close to the literature reported correlations.

目錄
中文摘要	I
英文摘要	II
目錄	III
圖目錄	VI
表目錄	XII
 第一章 緒論	1
1.1 前言	1
1.2 研究動機、範疇與方法	9
1.3 論文組織與架構	10
 第二章 文獻回顧	11
2.1	薄膜接觸器分離程序	11
2.1.1 薄膜蒸餾	11
2.1.2 薄膜吸收	12
2.2 紊流強化	15
2.2.1 粗糙面	15
2.2.2 間隔物強化	18
2.3 薄膜接觸器之計算流體力學模擬	21
 第三章 模擬系統	27
3.1 直接接觸式薄膜蒸餾平板模組	27
3.2 具間隔物平板模組	29
3.3 二氧化碳氣體吸收薄膜中空纖維模組	35
 第四章 薄膜熱質傳模式與物理性質	38
4.1 質傳模式	38
4.2 熱傳模式	43
4.3 鹽水系統物理性質	48
4.4 二氧化碳-碳酸丙烯酯系統物理性質	50
 第五章 計算流體力學模式建立	53
5.1 基本統制方程式	53
5.2 使用者定義函數	56
5.3 求解設定	57
5.3.1 離散方法	57
5.3.2 速度-壓力耦合方法	58
5.4 收斂準則	59
 第六章 網格無關化分析與模式驗證	60
6.1 鹽水系統直接接觸式平板模組	60
6.1.1 網格無關化分析	60
6.1.2 模式驗證	63
6.2 二氧化碳-碳酸丙烯酯系統中空纖維模組	68
6.2.1 網格無關化分析	68
6.2.2 模式驗證	71
 第七章 模擬結果與討論	74
7.1 鹽水脫鹽	78
7.1.1 基本個案內部特性	79
7.1.2 變因影響分析	97
7.1.3 關聯式比較	100
7.2 二氧化碳吸收	103
7.2.1 基本個案內部特性	104
7.2.2 變因影響分析	124
7.2.3 關聯式比較	127
 第八章 結論	131
 符號說明	134
 參考文獻	139

圖目錄
圖 1.1市售中空纖維型薄膜接觸器	1
圖 1.2薄膜蒸餾操作原理	4
圖 1.3薄膜蒸餾配置 (a)DCMD, (b)AGMD, (c)SGMD, (d)VMD (Meindersma et al, 2006)	5
圖 1.4二氧化碳薄膜吸收器	6
圖 1.5使用薄膜接觸器之二氧化碳吸收系統	6
圖 1.6薄膜層質傳機制	8
圖 1.7薄膜氣液接觸器之溫度極化與濃度極化	8
圖 2.1粗糙面高度與分離流程度	16
圖 2.2流體通道板 (a)開放式通道 (b)篩網式通道	20
圖 2.3網狀間隔物 (a)俯視圖 (b)側視圖	20
圖 2.4間隔物不同方位配置	23
圖 2.5間隔物幾何配置	23
圖 2.6間隔物之幾何參數(Meindersma et al., 2004)	24
圖 2.7狹窄通道模擬範圍與邊界定義	25
圖 2.8薄膜介面邊界之定義	26
圖 3.1直接接觸式薄膜蒸餾平板模組 (1)冷流體入口；(2)出口；(3)熱流體入口；(4)出口；(5)、(7)流體分離通道；(6)薄膜	27
圖 3.2直接接觸式薄膜蒸餾平板模組模擬系統	28
圖 3.3具間隔物模組	30
圖 3.4最小重複單元與關節點	31
圖 3.5間隔物接觸點處理	32
圖 3.6流體區塊內銳角處處理(水力角45o個案)	33
圖 3.7殼側自由表面	36
圖 3.8中空纖維模組模擬系統	36
圖 4.1薄膜內之質傳機制與阻力聯結	39
圖 4.2薄膜蒸餾與薄膜氣體吸收之熱傳機制與阻力聯結	44
圖 4.3複合膜結構與熱傳導阻力 (Martinez et al., 2008)	46
圖 5.1使用者定義函數(UDF)計算流程	57
圖 6.1鹽水系統直接接觸式平板模組尺寸	61
圖 6.2鹽水系統直接接觸式平板模組L2D3網格	62
圖 6.3鹽水系統直接接觸式平板模組驗證個案實驗數據(Martínez and Rodríguez-Maroto, 2007)	64
圖 6.4鹽水系統直接接觸式平板模組驗證結果	66
圖 6.5二氧化碳-碳酸丙烯酯系統中空纖維模組尺寸	68
圖 6.6二氧化碳-碳酸丙烯酯系統中空纖維模組L3D1網格	70
圖 6.7二氧化碳-碳酸丙烯酯系統中空纖維模組驗證個案實驗數據(Dindore et al, 2007)	71
圖 6.8二氧化碳-碳酸丙烯酯系統中空纖維模組模式驗證	73
圖 7.1 L8T90通道流體流動路徑示意圖	74
圖 7.2 L11T45通道流體流動路徑示意圖	75
圖 7.3鹽水系統L8T90通道基本個案速度等高分佈	80
圖 7.4鹽水系統L11T45通道基本個案速度等高分佈	81
圖 7.5鹽水系統L8T90與L11T45通道基本個案壓力等高分佈	82
圖 7.6鹽水系統L8T90與L11T45通道基本個案溫度等高分佈	83
圖 7.7鹽水系統L8T90與L11T45通道基本個案質量分率等高分佈	84
圖 7.8鹽水系統基本個案之質傳通量與膜面剪應力分佈	85
圖 7.9鹽水系統具基本個案之溫度與溫度極化係數分佈	88
圖 7.10鹽水系統基本個案之熱傳通量分佈	89
圖 7.11鹽水系統具基本個案之局部熱傳係數分佈	90
圖 7.12鹽水系統基本個案質傳通量	92
圖 7.13鹽水系統具基本個案熱側通道截線濃度分佈	94
圖 7.14鹽水系統基本個案之濃度與濃度極化係數分佈	95
圖 7.15鹽水系統之功率數隨雷諾數之變化(T1)	97
圖 7.16鹽水系統平均熱傳係數隨雷諾數之變化	98
圖 7.17鹽水系統之平均質傳通量隨進料速度之變化	99
圖 7.18鹽水系統摩擦因子模擬結果與關聯式之比較	100
圖 7.19鹽水系統熱傳係數模擬結果與關聯式之比較	101
圖 7.20鹽水系統雪耳伍德數模擬結果與關聯式之比較	102
圖 7.21二氧化碳-碳酸丙烯酯系統L8T90通道基本個案速度等高分佈	104
圖 7.22二氧化碳-碳酸丙烯酯系統L11T45通道基本個案速度等高分佈	105
圖 7.23二氧化碳-碳酸丙烯酯系統L8T90與L11T45通道基本個案壓力等高分佈	106
圖 7.24二氧化碳-碳酸丙烯酯系統L8T90與L11T45通道基本個案溫度等高分佈	107
圖 7.25二氧化碳-碳酸丙烯酯系統L8T90與L11T45通道基本個案質量分率等高分佈	108
圖 7.26二氧化碳-碳酸丙烯酯系統基本個案之質傳通量與膜面剪應力分佈	110
圖 7.27二氧化碳-碳酸丙烯酯系統基本個案之溫度分佈	112
圖 7.28二氧化碳-碳酸丙烯酯系統基本個案氣液側熱傳通量分佈	114
圖 7.29二氧化碳-碳酸丙烯酯系統基本個案之局部熱傳係數分佈	115
圖 7.30二氧化碳-碳酸丙烯酯系統基本個案質傳通量	117
圖 7.31二氧化碳-碳酸丙烯酯系統基本個案氣體通道截線濃度分佈	118
圖 7.32二氧化碳-碳酸丙烯酯系統基本個案液體通道截線濃度分佈	120
圖 7.33二氧化碳-碳酸丙烯酯系統基本個案之氣體側濃度與濃度極化係數分佈	122
圖 7.34二氧化碳-碳酸丙烯酯系統功率數隨雷諾數之變化	125
圖 7.35二氧化碳-碳酸丙烯酯系統平均熱傳係數隨雷諾數之變化	125
圖 7.36二氧化碳-碳酸丙烯酯系統之平均質傳通量隨進料速度之變化	126
圖 7.37二氧化碳-碳酸丙烯酯系統摩擦因子模擬結果與關聯式之比較	127
圖 7.38二氧化碳-碳酸丙烯酯系統氣體側納塞數模擬結果與關聯式之比較	128
圖 7.39二氧化碳-碳酸丙烯酯系統液體側納塞數模擬結果與關聯式之比較	129
圖 7.40二氧化碳-碳酸丙烯酯系統氣體側質傳係數模擬結果與關聯式之比較	130
圖 7.41二氧化碳-碳酸丙烯酯系統液體側雪耳伍德數模擬結果與關聯式之比較	130

表目錄
表 1.1利用薄膜接觸器之系統(Drioli et al., 2005)	2
表 1.2薄膜接觸器系統對應之傳統單元操作(Drioli et al., 2005)	2
表 1.3薄膜接觸器優缺點(Drioli et al., 2005)	3
表 3.1直接接觸式薄膜蒸餾薄膜與模擬通道規格	28
表 3.2具間隔物模組之規格	33
表 3.3具間隔物模組之模擬系統規格與薄膜特性	34
表 3.4中空纖維模組規格與薄膜特性	37
表 5.1求解設定	59
表 6.1鹽水系統直接接觸式平板模組網格無關化分析個案	61
表 6.2鹽水系統直接接觸式平板模組網格無關化分析結果	61
表 6.3純水系統直接接觸式平板模組驗證個案邊界條件設定	65
表 6.4鹽水系統直接接觸式平板模組驗證個案邊界條件設定	66
表 6.5純水系統直接接觸式平板模組透膜質傳通量比較	67
表 6.6鹽水系統直接接觸式平板模組透膜質傳通量比較	67
表 6.7二氧化碳-碳酸丙烯酯系統中空纖維模組網格無關化分析個案	69
表 6.8二氧化碳-碳酸丙烯酯系統中空纖維模組網格無關化分析結果	69
表 6.9二氧化碳-碳酸丙烯酯系統中空纖維模組驗證個案邊界條件設定	72
表 6.10二氧化碳-碳酸丙烯酯系統中空纖維模組透膜質傳通量比較	73
表 7.1層流摩擦因子關聯式	76
表 7.2具間隔物通道質傳係數關聯式	77
表 7.3鹽水系統直接接觸式模組之模擬個案	78
表 7.4二氧化碳-碳酸丙烯酯系統模擬個案	103

Ahmad, A. L., Lau, K. K., & Bakar, M. A. (2005). Impact of different spacer filament geometries on concentration polarization control in narrow membrane channel. Journal of Membrane Science, 262(1), 138-152.
Al-Sharif, S., Albeirutty, M., Cipollina, A., & Micale, G. (2013). Modelling flow and heat transfer in spacer-filled membrane distillation channels using open source CFD code. Desalination, 311, 103-112.
Amrei, S. M. H. H., Memardoost, S., & Dehkordi, A. M. (2014). Comprehensive modeling and CFD simulation of absorption of CO2 and H2S by MEA solution in hollow fiber membrane reactors. AIChE Journal, 60(2), 657-672.
Bird, R. B., Stewart, W. E., & Lightfoot, E. N. Transport Phenomena, Wiley, New York, 1960.
Reed, B. W., Semmens, M. J., & Cussler, E. L. (1995). Membrane contactors.Membrane Science and Technology, 2, 467-498.
Camacho, L. M., Dumée, L., Zhang, J., Li, J. D., Duke, M., Gomez, J., & Gray, S. (2013). Advances in membrane distillation for water desalination and purification applications. Water, 5(1), 94-196.
Cao, Z., Wiley, D. E., & Fane, A. G. (2001). CFD simulations of net-type turbulence promoters in a narrow channel. Journal of Membrane Science,185(2), 157-176.
 
Chernyshov, M. N., Meindersma, G. W., & De Haan, A. B. (2005). Comparison of spacers for temperature polarization reduction in air gap membrane distillation. Desalination, 183(1), 363-374.
Da Costa, A. R., Fane, A. G., Fell, C. J. D., & Franken, A. C. M. (1991). Optimal channel spacer design for ultrafiltration. Journal of Membrane Science,62(3), 275-291.
Da Costa, A. R., Fane, A. G., & Wiley, D. E. (1994). Spacer characterization and pressure drop modelling in spacer-filled channels for ultrafiltration. Journal of Membrane Science, 87(1), 79-98.
Dindore, V. Y., Brilman, D. W. F., Feron, P. H. M., & Versteeg, G. F. (2004). CO 2 absorption at elevated pressures using a hollow fiber membrane contactor. Journal of Membrane Science, 235(1), 99-109.
Dindore, V. Y., Brilman, D. W. F., & Versteeg, G. F. (2005). Modelling of cross-flow membrane contactors: physical mass transfer processes. Journal of Membrane Science, 251(1), 209-222.
Drioli, E., Criscuoli, A., & Curcio, E. (2011). Membrane Contactors: Fundamentals, Applications and Potentialities: Fundamentals, Applications and Potentialities (Vol. 11). Elsevier.
El-Bourawi, M. S., Ding, Z., Ma, R., & Khayet, M. (2006). A framework for better understanding membrane distillation separation process. Journal of Membrane Science, 285(1), 4-29.
Fane, A. G., Schofield, R. W., & Fell, C. J. D. (1987). The efficient use of energy in membrane distillation. Desalination, 64, 231-243.
 
Feron, P. H. M., & Jansen, A. E. (1997). The production of carbon dioxide from flue gas by membrane gas absorption. Energy conversion and management,38, S93-S98.
Fimbres-Weihs, G. A., & Wiley, D. E. (2010). Review of 3D CFD modeling of flow and mass transfer in narrow spacer-filled channels in membrane modules.Chemical Engineering and Processing: Process Intensification, 49(7), 759-781.
Glatzel, W. D., & Tomaz, M. C. (1966). The effect of turbulence promoters on heat transfer and pressure drop in plate heat exchanger. CSIR Special report CHEM, 53.
Goyal, N., Suman, S., & Gupta, S. K. (2015). Mathematical modeling of CO 2 separation from gaseous-mixture using a Hollow-Fiber Membrane Module: Physical mechanism and influence of partial-wetting. Journal of Membrane Science, 474, 64-82.
Ho, C. D., Chang, H., Chang, C. L., & Huang, C. H. (2013). Theoretical and experimental studies of flux enhancement with roughened surface in direct contact membrane distillation desalination. Journal of Membrane Science, 433, 160-166.
Hoff, K. A., & Svendsen, H. F. (2014). Membrane contactors for CO 2 absorption–Application, modeling and mass transfer effects. Chemical Engineering Science, 116, 331-341.
Khayet, M., Velázquez, A., & Mengual, J. I. (2004). Modelling mass transport through a porous partition: effect of pore size distribution. Journal of Non-Equilibrium Thermodynamics, 29(3), 279-299.
Koutsou, C. P., Yiantsios, S. G., & Karabelas, A. J. (2009). A numerical and experimental study of mass transfer in spacer-filled channels: effects of spacer geometrical characteristics and Schmidt number. Journal of Membrane Science, 326(1), 234-251.
Lawson, K. W., & Lloyd, D. R. (1997). Membrane distillation. Journal of Membrane Science, 124(1), 1-25.
Lewis, M. J. (1975). Optimising the thermohydraulic performance of rough surfaces. international Journal of Heat and Mass transfer, 18(11), 1243-1248.
Li, F., Meindersma, W., De Haan, A. B., & Reith, T. (2004). Experimental validation of CFD mass transfer simulations in flat channels with non-woven net spacers. Journal of Membrane Science, 232(1), 19-30.
Li, F., Meindersma, W., De Haan, A. B., & Reith, T. (2005). Novel spacers for mass transfer enhancement in membrane separations. Journal of Membrane Science, 253(1), 1-12.
Li, Y. L., & Tung, K. L. (2008). CFD simulation of fluid flow through spacer-filled membrane module: selecting suitable cell types for periodic boundary conditions. Desalination, 233(1), 351-358.
Martınez-Dıez, L., Vazquez-Gonzalez, M. I., & Florido-Dıaz, F. J. (1998). Study of membrane distillation using channel spacers. Journal of Membrane Science, 144(1), 45-56.
Martinez, L., & Rodríguez-Maroto, J. M. (2007). Effects of membrane and module design improvements on flux in direct contact membrane distillation.Desalination, 205(1), 97-103.
Martínez, L., & Rodríguez-Maroto, J. M. (2008). Membrane thickness reduction effects on direct contact membrane distillation performance. Journal of Membrane Science, 312(1), 143-156.
Mason, E. A., & Malinauskas, A. P. (1983). Gas transport in porous media: the dusty-gas model (Vol. 17). Elsevier Science Ltd.
Meindersma, G. W., Guijt, C. M., & De Haan, A. B. (2006). Desalination and water recycling by air gap membrane distillation. Desalination, 187(1), 291-301.
Phattaranawik, J., Jiraratananon, R., Fane, A. G., & Halim, C. (2001). Mass flux enhancement using spacer filled channels in direct contact membrane distillation. Journal of Membrane Science, 187(1), 193-201.
Phattaranawik, J., Jiraratananon, R., & Fane, A. G. (2003a). Heat transport and membrane distillation coefficients in direct contact membrane distillation.Journal of Membrane Science, 212(1), 177-193.
Phattaranawik, J., Jiraratananon, R., & Fane, A. G. (2003b). Effect of pore size distribution and air flux on mass transport in direct contact membrane distillation. Journal of Membrane Science, 215(1), 75-85.
Phattaranawik, J., Jiraratananon, R., & Fane, A. G. (2003c). Effects of net-type spacers on heat and mass transfer in direct contact membrane distillation and comparison with ultrafiltration studies. Journal of Membrane Science, 217(1), 193-206.
Razavi, S. M. R., Razavi, S. M. J., Miri, T., & Shirazian, S. (2013). CFD simulation of CO 2 capture from gas mixtures in nanoporous membranes by solution of 2-amino-2-methyl-1-propanol and piperazine. International Journal of Greenhouse Gas Control, 15, 142-149.
Schwinge, J., Wiley, D. E., Fane, A. G., & Guenther, R. (2000). Characterization of a zigzag spacer for ultrafiltration. Journal of Membrane Science, 172(1), 19-31.
Shakaib, M., Haq, M. E. U., Ahmad, I., & Yunus, R. M. (2010). Modeling the effect of spacer orientation on heat transfer in membrane distillation. World Academy of Science, Engineering and Technology, 72, 279-282.
Shirazian, S., Moghadassi, A., & Moradi, S. (2009). Numerical simulation of mass transfer in gas–liquid hollow fiber membrane contactors for laminar flow conditions. Simulation Modelling Practice and Theory, 17(4), 708-718.
Schock, G., & Miquel, A. (1987). Mass transfer and pressure loss in spiral wound modules. Desalination, 64, 339-352.
Van Medevoort, J., Jansen, A., Hanemaaijer, J. H., Dotremont, C., Nelemnas, B., Van Sonsbeek, E., & Biemans, R. (2008, November). Memstill: Seawater desalination a solution to water scarcity. In BMG-NMG Membrane Symposium, Belgian Membrane Group, Antwerp, Belgium.
Wang, M., Lawal, A., Stephenson, P., Sidders, J., & Ramshaw, C. (2011). Post-combustion CO 2 capture with chemical absorption: a state-of-the-art review. Chemical Engineering Research and Design, 89(9), 1609-1624.
Webb, R. L., Eckert, E. R. G., & Goldstein, R. (1971). Heat transfer and friction in tubes with repeated-rib roughness. International Journal of Heat and Mass Transfer, 14(4), 601-617.
Yang, X., Yu, H., Wang, R., & Fane, A. G. (2012). Analysis of the effect of turbulence promoters in hollow fiber membrane distillation modules by computational fluid dynamic (CFD) simulations. Journal of Membrane Science,415, 758-769.
Yu, H., Yang, X., Wang, R., & Fane, A. G. (2011). Numerical simulation of heat and mass transfer in direct membrane distillation in a hollow fiber module with laminar flow. Journal of Membrane Science, 384(1), 107-116.
Yu, H., Yang, X., Wang, R., & Fane, A. G. (2012). Analysis of heat and mass transfer by CFD for performance enhancement in direct contact membrane distillation. Journal of Membrane Science, 405, 38-47.
Zhang, H. Y., Wang, R., Liang, D. T., & Tay, J. H. (2006). Modeling and experimental study of CO 2 absorption in a hollow fiber membrane contactor.Journal of Membrane Science, 279(1), 301-310.
Zhang, H. Y., Wang, R., Liang, D. T., & Tay, J. H. (2008). Theoretical and experimental studies of membrane wetting in the membrane gas–liquid contacting process for CO 2 absorption. Journal of Membrane Science, 308(1), 162-170.
許建安. (2014). 薄膜蒸餾傳輸特性之計算流體力學研究. 淡江大學化學工程與材料工程學系碩士班學位論文, 1-154.