本研究主要是討論薄膜萃取系統並導入回流系統之解析解研究，可分為理論分析與實驗兩大部分。理論部份是針對共軛格拉茲問題的薄膜萃取系統進行理論分析，利用分離變數法、正交展開法及特徵函數以及數解展開來求得系統中的二維濃度分佈式。本文假設流體於兩平行板間的流體視為性質均勻之不可壓縮流體完全展開之速度分佈，忽略軸傳導，且流體物性不隨溫度改變，層流流動，進而去討論引入回流比、薄膜位置及水相溶液的體積流率的影響，其結果皆與二行程無迴流的狀況比較。結果顯示對於迴流比、薄膜位置及水相體積流率的增加，系統的萃取量也會隨之提升。
    本研究亦同時設計平板型薄膜萃取器實驗，並將其結果與二維理論值做比較，結果顯示二維質傳數學模型具有相當的準確性。

The mass transfer efficiency of parallel-plate membrane extractor module with recycle is investigated theoretically and experimentally in this study. The analytical solution is obtained by using an eigen-function expansion in terms of the power series. The experimental works in flat-plate module are conducted to confirm the accuracy of the theoretical results. The theoretical predictions were represented graphically with the volumetric flow rate, flow pattern, recycle rate, and subchannel thickness ratio as parameters compared with those obtained by the experimental runs. The influences of the design and operation parameters on the mass-transfer efficiency enhancement are also discussed.

中文摘要						I
英文摘要                                               II
目錄						III
圖目錄						VII
表目錄						XVII
第一章	緒論					1
1.1	前言					1
1.2	迴流效應對系統的影響			2
1.3	研究動機及目的				4
第二章	文獻回顧					5
2.1	文獻回顧					5
2.2	格拉茲問題				8
第三章	單行程液-液薄膜萃取器			11
3.1	順流式模型之理論分析			15
3.2	逆流式模型之理論分析			20
第四章	先順流後逆流式薄膜萃取器			24
4.1	雙行程先順後逆之無回流系統			24
4.2	雙行程先順後逆之具回流系統			31
4.2.1	左右等流出口迴流至入口模型	                  31
4.2.2	平板末端迴流至入口處模型	                  38
4.2.3	出口迴流至末端處模型			44
4.2.4	末端出口伴入口端迴流模型		         51
第五章	先逆流後順流式薄膜萃取器			58
5.1	雙行程先逆後順之無回流系統			58
5.2	雙行程先逆後順之具回流系統			65
5.2.1	左右等流出口迴流至入口模型	                  65
5.2.2	平板末端迴流至入口處模型	                  72
5.2.3	出口迴流至平板末端處模型	                  78
5.2.4	平板末端出口伴入口端迴流模型                  84
第六章	薄膜萃取實驗				91
6.1	實驗設備					91
6.2	藥品					92
6.3	實驗流程					94
6.4	實驗操作條件				95
6.5	注意事項					95
第七章	結果與討論			         97
7.1	先順流後逆流平板型薄膜萃取器		97
7.2	雙行程先順流後逆流式模型之結果與討論		110
7.2.1	左右等流出口迴流至入口			110
7.2.2	平板末端迴流至入口				120
7.2.3	出口迴流至平板末端				126
7.2.4	末端出口伴入口端迴流			132
7.2.5	雙行程先順流後逆流之四種迴流型式比較		132
7.3	雙行程先逆流後順流式模型之結果與討論		144
7.3.1	左右等流出口迴流至入口			144
7.3.2	平板末端迴流至入口				150
7.3.3	出口迴流至平板末端				156
7.3.4	末端出口伴入口端迴流			162
7.3.5	雙行程先逆流後順流之四種迴流型式比較		162
第八章	結論					175
第九章	未來研究方向				176
符號說明						177
參考文獻						179
附錄A	正交性質證明				184
附錄B	積分公式證明				188
附錄C	平板型系統展開係數				192
附錄D	平板系統平均濃度				203
圖目錄
圖3.1.1	順流式模型之薄膜萃取器示意圖		16
圖3.2.1	逆流式模型之薄膜萃取器示意圖		21
圖4.1.1	雙行程先順後逆無迴流之平板薄膜萃取器示意圖	25
圖4.2.1	雙行程先順後逆左右等流出口迴流至入口之薄膜萃取器示意圖						32
圖4.2.2	雙行程先順後逆平板末端迴流至入口處之薄膜萃取器示意圖						39
圖4.2.3	雙行程先順後逆平板出口迴流至平板末端之薄膜萃取器示意圖						45
圖4.2.4	雙行程先順後逆平板末端出口伴入口端迴流之薄膜萃取器示意圖						52
圖5.1.1	雙行程先逆流後順流無迴流之平板薄膜萃取器示意圖59
圖5.2.1	雙行程先逆流後順流左右等流出口迴流至入口之薄膜萃取器示意圖						66
圖5.2.2	雙行程先逆流後順流平板末端迴流至入口處之薄膜萃取器示意圖						73
圖5.2.3	雙行程平板出口迴流至平板末端之薄膜萃取器示意圖79
圖5.2.4	雙行程先逆流後順流平板末端出口伴入口端迴流之薄膜萃取器示意圖					85
圖6.3.1	平板型先順後逆之雙行程左右等流出口迴流至入口薄膜萃取系統實驗裝置圖					96
圖7.1.1	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同薄膜位置下，理論值與實驗值之無因次平均濃度與迴流比之關係						107
圖7.1.2	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同薄膜位置下，理論值與實驗值之總質傳速率與迴流比之關係						108
圖7.2.1	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在薄膜位置0.5萃餘相濃度變化分布圖(Qa=0.225cm3/s與Qb=0.25cm3/s)		                            112
圖7.2.2	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在薄膜位置0.5萃餘相濃度變化分布圖(Qa=0.45cm3/s與Qb=0.25cm3/s)					113
圖7.2.3	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					114
圖7.2.4	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					115
圖7.2.5	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					116
圖7.2.6	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					117
圖7.2.7	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖          118
圖7.2.8	先順流後逆流之雙行程左右等流出口迴流至入口系統中，在不同薄膜位置下，平均謝塢數對格拉茲數Gzb作圖          119
圖7.2.9	雙行程先順流後逆流平板末端迴流至入口系統中，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)				          121
圖7.2.10	雙行程先順流後逆流平板末端迴流至入口系統中，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)		                            122
圖7.2.11	雙行程先順流後逆流平板末端迴流至入口系統中，在不同的薄膜位置下，萃取速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					123
圖7.2.12	雙行程先順流後逆流平板末端迴流至入口系統中，在不同的薄膜位置下，萃取速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					124
圖7.2.13	雙行程先順流後逆流平板末端迴流至入口系統中，在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖	         125
圖7.2.14	雙行程先順流後逆流出口迴流至平板末端系統中，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					127
圖7.2.15	雙行程先順流後逆流出口迴流至平板末端系統中，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					128
圖7.2.16	雙行程先順流後逆流出口迴流至平板末端系統中，在不同的薄膜位置，萃取速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					129
圖7.2.17	雙行程先順流後逆流出口迴流至平板末端系統中，在不同的薄膜位置，萃取速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					130
圖7.2.18	雙行程先順流後逆流出口迴流至平板末端系統中，在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖	         131
圖7.2.19	雙行程先順流後逆流平板末端出口伴入口端迴流系統中，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					133
圖7.2.20	雙行程先順流後逆流平板末端出口伴入口端迴流系統中，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					134
圖7.2.21	雙行程先順流後逆流平板末端出口伴入口端迴流系統中，在不同的薄膜位置下，萃取速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					135
圖7.2.22	雙行程先順流後逆流平板末端出口伴入口端迴流系統中，在不同的薄膜位置下，萃取速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					136
圖7.2.23	雙行程先順流後逆流平板末端出口伴入口端迴流系統中，在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖		137
圖7.2.24	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.5與Qa=0.225cm3/s				138
圖7.2.25	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.5與Qa=0.45cm3/s				139
圖7.2.26	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.25與Qa=0.225cm3/s				140
圖7.2.27	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.25與Qa=0.45cm3/s				141
圖7.2.28	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.75與Qa=0.225cm3/s				142
圖7.2.29	四種迴流形式先順流後逆流之總質傳速率與迴流比之比較Delta=0.75與Qa=0.45cm3/s				143
圖7.3.1	雙行程先逆流後順流雙行程左右等流出口迴流至入口，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					145
圖7.3.2	雙行程先逆流後順流雙行程等流出口迴流至入口，在不同的薄膜位置下，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					146
圖7.3.3	雙行程先逆流後順流雙行程左右等流出口迴流至入口，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					147
圖7.3.4	雙行程先逆流後順流雙行程左右等流出口迴流至入口，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					148
圖7.3.5	雙行程先逆流後順流左右等流出口迴流至入口在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖			149
圖7.3.6	雙行程先逆流後順流平板末端迴流至入口處，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					151
圖7.3.7	雙行程先逆流後順流平板末端迴流至入口處，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					152
圖7.3.8	雙行程先逆流後順流平板末端迴流至入口處，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)						153
圖7.3.9	雙行程先逆流後順流平板末端迴流至入口處，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)                                          154
圖7.3.10	雙行程先逆流後順流平板末端迴流至入口處，在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖		         155
圖7.3.11	雙行程先逆流後順流出口迴流至平板末端，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)						157
圖7.3.12	雙行程先逆流後順流出口迴流至平板末端，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)						158
圖7.3.13	雙行程先逆流後順流出口迴流至平板末端，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)						159
圖7.3.14	雙行程先逆流後順流出口迴流至平板末端，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)						160
圖7.3.15	雙行程先逆流後順流出口迴流至平板末端在不同薄膜位置下，平均謝塢數對格拉茲數作Gza圖			161
圖7.3.16	雙行程先逆流後順流平板末端出口伴入口端迴流，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					164
圖7.3.17	雙行程先逆流後順流平板末端出口伴入口端迴流，在不同的薄膜位置，總質傳速率與迴流比的關係(Qa=0.45cm3/s與Qb=0.25cm3/s)					165
圖7.3.18	雙行程先逆流後順流平板末端出口伴入口端迴流，在不同的薄膜位置，萃取效率與迴流比的關係(Qa=0.225cm3/s與Qb=0.25cm3/s)					166
圖7.3.19	雙行程先逆流後順流平板末端出口伴入口端迴流，在不同的薄膜位置，萃取效率與迴流比的關係
(Qa=0.45cm3/s與Qb=0.25cm3/s)			167
圖7.3.20	雙行程先逆流後順流平板末端出口伴入口端迴流在不同薄膜位置下，平均謝塢數對格拉茲數Gza作圖	         168
圖7.3.21	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.5與Qa=0.225cm3/s)			         169
圖7.3.22	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.5與Qa=0.45cm3/s)				170
圖7.3.23	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.25與Qa=0.225cm3/s)				171
圖7.3.24	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.25與Qa=0.45cm3/s)				172
圖7.3.25	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.75與Qa=0.225cm3/s)				173
圖7.3.26	四種迴流形式先逆流後順流之總質傳速率與迴流比之比較(Delta=0.5與Qa=0.45cm3/s)				174
 
表目錄
表7.1.1	順流系統在薄膜位置為0.5之特徵值		99
表7.1.2	順流系統在薄膜位置為0.25之特徵值		100
表7.1.3	順流系統在薄膜位置為0.75之特徵值		101
表7.1.4	逆流系統在薄膜位置為0.5之特徵值		102
表7.1.5	逆流系統在薄膜位置為0.25之特徵值		103
表7.1.6	逆流系統在薄膜位置為0.75之特徵值		104
表7.1.7	先順流後逆流式薄膜位置為0.5之萃取與萃餘相第一行程係數						105
表7.1.8	先順流後逆流式薄膜位置為0.5之萃取與萃餘相第一行程係數						106
表7.1.9	實驗值與理論值之平均濃度與總質傳速率之誤差	109

1.M. C. Porter, “Handbook of Industrial Membrane Technology”, Noyes Publications, New Jersey, pp. 175, 1990.
2.郭文正和曾添文, 薄膜分離, 高立圖書公司, pp. 1-18 , 1988
3.W. S. W. Ho and K. K. Sirkar, “Membrane Handbook”, Van Nostrand Reinhold, New York, 1992.
4.P. M. Bungy, H. K. Lonsdale, and M. N. de Pinho, “Synthetic Membranes: Science, Engineering, and Applications”, D. Reidel, Publishing Company, pp.3-8, 1986.
5.Warren L. McCabe, Julian C. Smith, Peter Harriott, Unit Operations of Chemical Engineering, McGraw-Hill, New York (1993).
6.S. W. Tsai, and H. M. Yeh , “A study of the Separation Efficiency in Horizontal Thermal Diffusion Columns with External Refluxes”, Can. J. Chem. Eng., 63 (1985), 406.
7.R. Marquart , and H. Blenke, “Circulation of Moderately to Highly Viscous Newtonian and Non-Newtonian Liquids in Propeller-Pumped Circulation Loop Reactor”, Int. Chem. Eng., 20 (1980), 368.
8.R. Marquart, “Circulation of High-Viscosity Newtonian and Non-Newtonian Liquids in Jet Loop Reactor”, Int. Chem. Eng., 20 (1981), 399.
9.G. Dussap, and J. B. Gros, “Energy Consumption and Interfacial Mass Transfer Area in an Air-Lift Fermentor”, Chem. Eng. J., 25, (1982) 151.
10.Siegel, M. H. J. C. Merchuk, and K. Schugerl, “Air-Lift Reactor Analysis: Interrelationships Between Riser , Downcomer, and Gas-Liquid Separator Behavior, Including Gas Recirculation Effects”, AiChE J., 32 (1986), 1585.
11.R. K. Shah., and A. L. London, Laminar Flow Forced Convection in Ducts, Academic Press, New York, 1987, 196.
12.V-D. Dang., and M. Steinberg, Convective Diffusion with homogeneous and heterogeneous reaction in a tube, J. Phys. Chem., 84 (1980) 214.
13.T. L. Perelman, On conjugated problems of heat transfer, Int. J. Heat Mass Transfer, 3 (1961) 293.
14.R. J. Nunge, and W. N. Gill, Analysis of Heat or Mass Transfer in Some Countercurrent Flows, Int. J. Heat Mass Transfer, 8 (1965) 873.
15.D. Murkerjee, and E. J. David, Direct-contact heat transfer immiscible fluid layers in laminar flow, AIChE J., 18 (1972) 94.
16.S. S. Kim, and D. O. Cooney, Improved theory for hollow-fiber enzyme reactor, Chem. Eng. Sci., 31 (1976) 289.
17.E. J. David, and S. Venkatesh, The solution of conjugated multiphase heat and mass transfer problems, Chem. Eng. Sci., 34 (1979) 775.
18.E. Papoutsakis, and D. Ramkrishna, Conjugated Graetz problems. I:general formalism and a class of solid-fluid problem, Chem. Eng. Sci., 36 (1981) 1381.
19.E. Papoutsakis, D. Ramkrishna, Conjugated Graetz problems. II: fluid-fluid problems, Chem. Eng. Sci., 36 (1981) 1393.
20.X. Yin, H. H. Bau, The conjugated Graetz problem with axial conduction, Trans. ASME, 118 (1996) 482.
21.H. M. Yeh, S. W. Tsai and C. L. Chiang, “Recycle Effects on Heat and Mass Transfer through a Parallel-Plate Channel”, AIChE J., 33 (1987) ,1743.
22.H. M. Yeh, Yu-Shu Hsu, Analysis of membrane extraction through rectangular mass Exchangers, Chem. Eng. Sci., 54 (1999) 897.
23.G. M. Brown, Heat or Mass Transfer in a Fluid in Laminar Flow in a Circular or Flat Conduit, AIChE J., 6 (1960) 179.
24.T. L. Perelman, On Conjugated Problems of Heat Transfer, Int. J. Heat Mass Transfer, 3 (1961) 293.
25.A. P. Hatton and A. Quarmby, Heat Transfer in the Thermal Entry Length with Laminar Flow in an Annulus, Int. J. Heat Mass Transfer, 5 (1962) 973.
26.R. J. Nunge and W. N. Gill, Analysis of Heat or Mass Transfer in Some Countercurrent Flows, Int. J. Heat Mass Transfer, 8 (1965) 73.
27.R. J. Nunge and W. N. Gill, An Analytical Study of Laminar Counterflow Double-Pipe Heat Exchangers. AIChE J., 12 (1966) 279.
28.C. J. Hsu, Heat Transfer in a Round Tube with Sinusoidal Wall Heat Flux Distribution, AIChE J., 11 (1965) 690.
29.E. J. Davis, Exact Solutions for a Class of Heat and Mass Transfer Problems, Can. J. Chem. Eng., 51 (1973) 562.
30.H. M. Yeh, T. W. Chang and S. W. Tsai, A Study of the Graetz Problems in Concentric-Tube Continuous-Contact Countercurrent Separation Process with Recycles at Both Ends, Separation Science and Technology, 21 (1986) 403.
31.M. A. Ebadian and H. Y. Zhang, An Exact Solution of Extend Graetz Problems with Axial Heat Conduction, Int. J. Heat Mass Transfer, 32 (1989) 1709.
32.X. Yin and H. H. Bau, The Conjugated Graetz Problem with Axial Conduction, Journal of Heat Transfer, 118 (1996) 428.
33.C. D. Ho, H. M. Yeh and W. S. Sheu, An analytical study of heat and mass transfer through a parallel-plate channel with recycly, Int. J. Heat Mass Transfer, 44 (1998) 2589.
34.A. Kiani, R. R. Bhave and K. K. Sirkar, “Solvent Extraction with Immobilized Interfaces in a Microporous Hydrophobic Membrane”, J. Membrane Sci., 20 (1984)123.
35.N. A. D’ Elia, L. Dahuron and E. L. Cussler, “Liquid-Liquid Extraction with Microporous Hollow Fibers”, J. Membrane Sci., 29(1986) 309.
36.R. Prasad, A. Kian, R. R. Bhave and K. K. Sirkar, “Further Studies on Solvent Extraction with Immobilized Interfaces In a Microporous Hydrophilic Membrane”, J. Membrane Sci., 26 (1986)79.
37.R. Prasad, and K. K. Sirkar, “Solvent Extraction with Microporous Hydrophilic and Composite Membrane”, AChE J., 33 (1987) 1057.
38.R. Prasad, and K. K. Sirkar, “Microporous Membrane Solvent Extraction”, Separation Science and Technology, 22 (1987) 619.
39.P. R. Alexander and R. W. Callahan, “Liquid-Liquid Extraction and Stripping of Gold with Microporous Hollow Fibers”, J. Membrane Sci., 35(1987) 57.
40.A. Sengupta, R. Basu and K. K. Sirkar, “Separation of Solute from Aqueous by Contained Liquid Membranes”, AIChE J., 34 (1988)1698.
41.S. Majumdar, A. K. Guha and K. K. Sirkar, “A New Liquid Membrane Technique for Gas Separation”, AIChE, J., 34 (1988) 1135.
42.C. D. Ho, H. M. Yeh and R. C. Wang, Heat-transfer Enhancement in Double-Pass Flat-Plate Solar Air Heaters with Recycle, Energy, 30 (2005) 2796.