翁宗賢臺灣大學:應用力學研究所莊峻維Zhuang, Jun-WeiJun-WeiZhuang2010-06-022018-06-292010-06-022018-06-292008U0001-3107200808495800http://ntur.lib.ntu.edu.tw//handle/246246/184814摘要於兩種或兩種以上微量流體之快速混合，在微型系統與裝置裡是一項極為重要的主題，也是影響系統功能的重要因素。論文提出一新型被動式微混合器之設計，先以熱流數值摹擬軟體FLUENT分析微流道內的流場狀態、混合指標及壓力係數，進而改變流道結構、尺寸及擋塊角度，比較流場特性，設計出符合低壓降、高混合效果之微混合器。據數值摹擬結果，本文研製平面式微混合器，以矽晶圓作為基材，利用電感耦合電漿蝕刻機(ICP)蝕刻出母模，進而以聚二甲基矽氧烷(polydimethysiloxane, PDMS)翻製出模型，最後以氧電漿接合上玻璃上蓋，完成微混合器製作。實驗量測方法則利用螢光染料及食用色素分別加入甲醇與去離子水做為工作流體，經由影像分析軟體將CCD攝影機所擷取到之混合影像灰階化，再分析比較各種流道之混合效果，目的在設計出結構簡單、低壓降、高混合效能之微混合器。其中3-U turns & 2 streamlined obstacles per turn with 5 degree rotated之設計在Re = 30時具有0.79之混合指標，與蛇型流道相同，但僅需要蛇型流道壓降值之 。因此，本文所研製的微混合器將可應用於微型化直接甲醇燃料電池(Micro Direct Methanol Fuel Cell, μDMFC)有實質幫助。Abstractt is an important topic for fast mixing two or more kinds of different μfluids in the micro systems and devices. It is also a significant factor that affected the functions of the systems.his thesis proposes a novel design of passive mixers. First, we analysis the flow fields, mixing indexes, and pressure coefficients by the commercial computational fluid dynamics software, FLUENT. In addition, we design a μmixer which meets the demands of low pressure drop, and high mixing performance by changing the μchannel structure size, and the orientation of the obstacle.ccording to the results of the simulations, the fabrication of the planar μmixers was based on the silicon wafer. After etching by the inductive coupling plasma etcher (ICP), we successfully fabricated our molds. Then we poured the PDMS on the molds to make our mixers and lifted the mixers from the molds. Finally, by using the oxygen plasma, we bonded the mixers with glasses to complete the chips. We use the pure methanol with fluorescent dye and D.I. water with food coloring as the working fluids. By using the charge-coupled device to capture the images, and employing the image analyzing software to quantity the concentration distribution, we can compare all the designs to find the simple structured, lowest pressure drop, and highest mixing performance mixer. Among our designs, 3-U turns & 2 streamlined obstacles per turn with 5 degree rotated design can meet the demand. At Re = 30, it has 0.79 of the mixing index, same as the serpentine design. However, it just needs three-quarter of the pressure coefficient. And we hope our design will essentially help the development of the μ-DMFC.目錄要 IBSTRACT II錄 III目錄 V目錄 VI號表 XIV一章 緒論 1-1 研究動機 1-2 文獻回顧 3-2-1 主動式微混合器 3-2-2 被動式微混合器 4-3 本文內容 8二章 微混合器之理論、設計與數值模擬 9-1 統御方程式 9-2 模型之設計、建立 10-3 數值方法與模擬 12-4 混合指標 13-5 壓力係數 14三章 數值模擬結果與分析 16-1 基本假設條件、初始條件及邊界參數設定 16-2 網格獨立性測試 17-3 混合指標與壓力係數之結果分析 17四章 微混合器之製作與量測 22-1 晶片製造流程 22-2 晶片之製作程序 22-2-1 光罩之製作 23-2-2 晶圓之清潔 23-2-3 晶圓母模之製作 24-2-4 PDMS轉印製程 26-2-5 PDMS晶片接合 27-3 流道製作結果 28-4 混合指標之量測 28五章 結果與討論 30六章 結論與未來展望 32-1 結論 32-2 未來展望 32考文獻 33目錄1 RCA清洗流程 37目錄1-1 (a)~(c)壓力場擾動 (d)電場擾動 (e)介電泳擾動 (f)、(g)電動力擾動 381-2 Tsai等人[10]利用熱氣泡產生熱擾動，增加混合效果 381-3 (a)典型T字型混合器 (b)Y字型混合器 (c)平行層疊式示意圖 (d) 水力聚焦型混合器 391-4 串聯層疊式混合器之概念示意圖 391-5 Melin等人[19]利用表面張力及幾何配置所設計之混合器 391-6 Miyake等人[20]所設計之注射式混合器 401-7 Wong等人[22]在壁面上佈入阻塊以提升混合效果 401-8 Lin等人[23]在流道內佈入圓柱型阻塊以增加混合效果 411-9 Hong等人[24] 運用Coanda效應去製造流體回流以增加混合效果 411-10 Liu等人[25] 設計一以混沌對流方式增加混合的3D蛇型流道(上)，並與方波型流道(中)及長直型流道(下)比較之 421-11 Stroock等人[27] 以人字型的凹槽設計產生混沌對流效應，增加混合效果 422-1 流道混合腔體設計圖 432-2 流道尺寸、阻塊尺寸設計圖 432-3 3-U turns & 2 streamlined obstacles per turn with 5 degree rotated所切之局部網格 442-4 數值疊代求解流程圖 453-1 網格測試之混合指標對網格數關係 463-2 網格測試之壓力係數對網格數關係 463-3 網格測試之混合效率對網格數關係 473-4 Y-tpye 在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 483-5 Y-tpye 在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 483-6 Y-tpye 在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 493-7 Y-tpye 在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 493-8 Y-tpye 在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 503-9 2-U turns & 2-streamlined obstacles per turn在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 503-10 2-U turns & 2-streamlined obstacles per turn在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 513-11 2-U turns & 2-streamlined obstacles per turn在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 513-12 2-U turns & 2-streamlined obstacles per turn在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 523-13 2-U turns & 2-streamlined obstacles per turn在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 523-14 2-U turns & 2-streamlined obstacles per turn (open)在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 533-15 2-U turns & 2-streamlined obstacles per turn (open)在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 533-16 2-U turns & 2-streamlined obstacles per turn (open)在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 543-17 2-U turns & 2-streamlined obstacles per turn (open)在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 543-18 2-U turns & 2-streamlined obstacles per turn (open)在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 553-19 2-U turns serpentine在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 553-20 2-U turns serpentine在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 563-21 2-U turns serpentine在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 563-22 2-U turns serpentine在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 573-23 2-U turns serpentine在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 573-24 2-U turns & 3-streamlined obstacles per turn在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 583-25 2-U turns & 3-streamlined obstacles per turn在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 583-26 2-U turns & 3-streamlined obstacles per turn在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 593-27 2-U turns & 3-streamlined obstacles per turn在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 593-28 2-U turns & 3-streamlined obstacles per turn在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 603-29 2-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 603-30 2-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 613-31 2-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 613-32 2-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 623-33 2-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 623-34 2-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 633-35 2-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 633-36 2-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 643-37 2-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 643-38 2-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 653-39 3-U turns & 2-streamlined obstacles per turn在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 653-40 3-U turns & 2-streamlined obstacles per turn在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 663-41 3-U turns & 2-streamlined obstacles per turn在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 663-42 3-U turns & 2-streamlined obstacles per turn在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 673-43 3-U turns & 2-streamlined obstacles per turn在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 673-44 3-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 683-45 3-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 683-46 3-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 693-47 3-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 693-48 3-U turns & 2-streamlined obstacles per turn (5 degree rotated)在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 703-49 3-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 703-50 3-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 713-51 3-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 713-52 3-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 723-53 3-U turns & 2-streamlined obstacles per turn (10 degree rotated)在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 723-54 3-U turns serpentine在Re=10時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 733-55 3-U turns serpentine在Re=30時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 733-56 3-U turns serpentine在Re=50時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 743-57 3-U turns serpentine在Re=70時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 743-58 3-U turns serpentine在Re=90時之流道混合狀況，左圖為流道平面，右圖為混合位置截面混合濃度分佈狀況 753-59 case1各流道混合指標對雷諾數之比較圖 763-60 case1各流道壓力係數對雷諾數之比較圖 763-61 case1各流道混合效率對雷諾數之比較圖 773-62 case2各流道混合指標對雷諾數之比較圖 773-63 case2各流道壓力係數對雷諾數之比較圖 783-64 case2各流道混合效率對雷諾數之比較圖 783-65 case3各流道混合指標對雷諾數之比較圖 793-66 case3各流道壓力係數對雷諾數之比較圖 793-67 case3各流道混合指標除以壓力係數對雷諾數之比較圖 804-1 微流道製程流程圖 814-2 微混合器之光罩設計 814-3 蝕刻製程示意圖 824-4 蝕刻成功之母模 834-5 PDMS流道量測之深度 834-6 氧電漿產生器 844-7 製作完成之實驗晶片 844-8 針筒式注入泵浦 845-1 Y-type流道混合影像，左側為甲醇，右側為水 855-2 2-U turns serpentine流道混合影像，左側為甲醇，右側為水 865-3 3-U turns & 2 streamlined obstacles per turn with 5 degree rotated流道混合影像，左側為甲醇，右側為水 875-4 Y-type實驗及模擬混合指標比較圖 885-5 2-U turns serpentine實驗及模擬混合指標比較圖 885-6 3-U turns & 2 streamlined obstacles per turn with 5 degree rotated實驗及模擬混合指標比較圖 89application/pdf4623052 bytesapplication/pdfen-US微流體微混合器計算流力μDMFCμfluidsμmixerCFDthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/184814/1/ntu-97-R95543067-1.pdf