1、週期性干擾週期性干擾T T型微流道流體混合型微流道流體混合之模擬分析之模擬分析 Numerical Study of Fluid Mixing in a T-type Micromixer Subject to Periodic Pressure Disturbances 邱煥釗邱煥釗1,2、陳志敏陳志敏1 11 1國立中興大學國立中興大學 機械工程學系機械工程學系2勤益科技大學機械工程學系勤益科技大學機械工程學系 中國機械工程學會第二十六屆全國學術研討會成功大學 台南November 21&22,2009大綱大綱n研究目的與方法n統御方程式 n數值方法n模擬結果分析n結論研究目的與方法研究目
2、的與方法目的：n模擬微T型混合器加入週期性壓力干擾方式之三維動態流場現象與混合效果 方法：n利用CFD軟體Fluent 6.3作為數值模擬工具n分析T型混合器的動態濃度場與流場分析假設分析假設n兩流體接觸後不會產生化學反應n兩種工作流體同為水，物理性質相同 n不可壓縮流且具有固定黏度 nRe100,為層流流動統御方程式統御方程式n連續方程式 continuity equation (1)n動量方程式 Navier-Stokes equations (2)n濃度分布的擴散-對流方程式 Ficks law (3)0V2DVpVDt 2cVcDct 數值方法數值方法模型建立模型建立模型建立：由Gam
3、bit 2.3所建構而成網格(Mesh)建立：均勻網格(Uniform Mesh)(Elements：Hex/Wedge，Type：Cooper)網格數目：控制在10萬左右網格獨立性測試網格獨立性測試n以寬深比 AR=10(流道寬度:500m)n無壓力振盪n邊界條件:入口流速Vin=0.028 m/s 出口壓力 Pout=0 n流體性質:=1.03e-3 kg/ms D=1e-8 m2/s不同網格數出口斷面中心線速度不同網格數出口斷面中心線速度分佈分佈-3-2-101234x 10-400.010.020.030.040.050.060.070.080.09x(mm)velocity(m/s)
4、出 口 斷 面 中 心 線 速 度 分 佈mesh no.=38750mesh no.=77500mesh no.=125850mesh no.=310000數值方法數值方法邊界條件邊界條件入口條件：均勻的入口速度，左右入口設定流體A與流體B。出口條件：設定為大氣壓力，即錶壓力為0 邊界條件：壁面設定為無滑動條件(no slip condition)左右振盪流道之振盪壓力設定為：01 sin(2)inppft濃度場之邊界條件濃度場之邊界條件 n左側流道入口與左側振盪流道之濃度0，右側入口流道與右側振盪流道之濃度為1 n壁面為濃度通量為0(zero diffusive flux)即 ，為壁面之單
5、位法向量。0cnn混合器流道結構混合器流道結構示意圖示意圖01 sin(2)inppftPout=0Vin=constant邊界條件邊界條件AR=b/h數值模擬步驟使用軟體:Fluent 6.3Gambit 2.3Tecplot 360模擬參數模擬參數流體性質:黏滯係數=1.03e-3 kg/ms 擴散係數 D=1e-8 m2/s幾何尺寸:出口長度:18.25mm 微流道斷面:2.510-2 mm2 流道寬深比AR(b/h):1.610入口流速範圍:0.0280.084 m/s(Re=2.948.82)干擾壓力振盪的相位差範圍:0振盪頻率:050Hz量化混合程度量化混合程度-定義混合效率定義混
6、合效率 211()(1)100%niiCCnC出口處斷面之混合效率公式：(Xia et al.,2005)Ci：格點i上的濃度C：完全混合的濃度n：格點數C=0C=1C=0.5無壓力振盪無壓力振盪干擾干擾時的混合情形時的混合情形 不同入口流速下,不同寬深比出口處數值混合效率的比較Vin(m/s)0.0280.0560.084混合效率(%)(%)AR=1.6=1.634.121.116.2AR=6.4=6.421.714.611.6AR=10=1019.312.810.0AR=10=10(實驗值實驗值)18.312.510 AR=6.4,不同入口流速，出口剖面之濃度分佈截面 誤差0.55%壓力振
7、盪干擾壓力振盪干擾=0=0時的混合情形時的混合情形 00.511.52051015202530354030HzTime(sec)Mmixing efficiency(%)2Hz4Hz50HzAR=6.4，Vin=0.056 m/s，出口處混合效率隨時間變化的情形。流場與濃度場分析AR=6.4，Vin=0.056 m/s，=0，在頻率為10 Hz，t=1.0 s時n流道中間剖面及其混合流道在各位置之流線圖。n出口及中央剖面之濃度分佈圖。流線圖(a)出口剖面(b)中央剖面濃度分佈圖壓力振盪干擾壓力振盪干擾=0之混合效益之混合效益00.511.522.501020304050607080902HzT
8、ime(sec)Mmixing efficiency(%)4Hz8Hz10Hz20Hz30Hz40Hz50HzAR=6.4，Vin=0.056 m/s，=0混合效率隨頻率變化情形。AR=6.4，Vin=0.056 m/s，=，在不同之振盪頻率下，混合效率隨時間變化的情形。010203040501020304050607080Freguency(Hz)Mixing efficiency(%)=0=0.2=0.4=0.5=0.8=流場流場與與濃度場分析濃度場分析(2)(2)AR=6.4，Vin=0.056 m/s，=，在頻率為10 Hz，t=1.0 s時n流道中間剖面及其混合流道在各位置之流線圖。
9、n出口及中央剖面之濃度分佈圖。流線圖流線圖(2)(2)濃度分佈圖濃度分佈圖(2)(2)(a)出口剖面(b)中央剖面入口流速對混合效益之影響入口流速對混合效益之影響 AR=6.4，=，不同入口流速下之混合效率隨頻率變化的情形。01020304050102030405060708090Freguency(Hz)Mixing efficiency(%)Vin=0.028(m/s)Vin=0.056(m/s)Vin=0.084(m/s)流道寬深比對混合效益的影響流道寬深比對混合效益的影響 Vin=0.056 m/s，=時，不同寬深比AR，混合效率隨頻率變化的情形。0102030405010203040
10、50607080Freguency(Hz)Mixing efficiency(%)AR=1.6AR=6.4AR=10流場與濃度場分析流場與濃度場分析(3)(3)AR=1.6，Vin=0.056 m/s，=，在頻率 為10 Hz，t=1.0 s時n流道中間剖面及其混合流道在各位置之流線圖。n出口及中央剖面之濃度分佈圖。流線圖流線圖(3)(3)濃度分佈圖濃度分佈圖(3)(3)(a)出口剖面(b)中央剖面AR=1.6，=，入口流速對混合效入口流速對混合效益之影響益之影響。=96%010203040502030405060708090100Freguency(Hz)Mixing efficiency(
11、%)Vin=0.028(m/s)Vin=0.056(m/s)AR=1.6，Vin=0.028 m/s，=，在在20 Hz，t=1.0 s時之濃度分佈。時之濃度分佈。(a)出口剖面(b)中央剖面=96%結論1.利用兩側振盪流道加入壓力干擾，確實可提高混合效率約25倍2.兩側振盪流道導入無相位差的壓力干擾，可使混合效率提高約2倍，且混合效率並不會隨頻率變化。3.兩側導入具有相位差之壓力干擾，則會有較佳的混合效果，並且混合效率會隨振動頻率先升高達到最佳值後會隨頻率升高而下降的趨勢。4.在=時，流體界面會有最大彎曲，偏折的現象而具有最大之接觸面積，在以擴散主宰的低雷諾數流動範圍，會有最佳的混合效果。4
12、.流道寬比深AR愈低，混合效率愈好。在AR=1.6 時分別在軸向及橫向剖面都會產生渦流，使流體界面有摺疊(folding)，扭曲(distortion)的現象，更能幫助兩流體的混合。5.在較低的入口流速，有較長的擴散時間，可提高混合的效率。6.本研究中當流道寬深比AR=1.6，入口流速Vin=0.028 m/s，兩側壓力振盪=，兩流體的混合效率可達到=96%。報告完畢謝謝大家不同網格數出口斷面混合效率的不同網格數出口斷面混合效率的比較比較網格數混合效率(%)3875019.227750019.2512585019.313100019.39混合流道不同位置剖面流場觀察區混合流道不同位置剖面流場觀察區
