1 / 38

週期性干擾 T 型微流道流體混合之模擬分析

週期性干擾 T 型微流道流體混合之模擬分析. 邱煥釗 1,2 、 陳志敏 1 1 國立中興大學 機械工程學系 2 勤益科技大學機械工程學系 中國機械工程學會第二十六屆全國學術研討會 成功大學 台南 November 21&22, 2009. Numerical Study of Fluid Mixing in a T-type Micromixer Subject to Periodic Pressure Disturbances. 大綱. 研究目的與方法 統御方程式 數值方法 模擬結果分析 結論. 研究目的與方法. 目的：

cleavant-mauricio
Download Presentation

週期性干擾 T 型微流道流體混合之模擬分析

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 週期性干擾T型微流道流體混合之模擬分析 • 邱煥釗1,2、陳志敏1 • 1國立中興大學 機械工程學系 • 2勤益科技大學機械工程學系 • 中國機械工程學會第二十六屆全國學術研討會 成功大學 台南 November 21&22, 2009 Numerical Study of Fluid Mixing in a T-type Micromixer Subject to Periodic Pressure Disturbances

  2. 大綱 • 研究目的與方法 • 統御方程式 • 數值方法 • 模擬結果分析 • 結論

  3. 研究目的與方法 目的： • 模擬微T型混合器加入週期性壓力干擾方式之三維動態流場現象與混合效果 方法： • 利用CFD軟體Fluent 6.3作為數值模擬工具 • 分析T型混合器的動態濃度場與流場

  4. 分析假設 • 兩流體接觸後不會產生化學反應 • 兩種工作流體同為水，物理性質相同 • 不可壓縮流且具有固定黏度 • Re<100,為層流流動

  5. 統御方程式 • 連續方程式 continuity equation (1) • 動量方程式 Navier-Stokes equations (2) • 濃度分布的擴散-對流方程式 Fick’s law (3)

  6. 數值方法—模型建立 模型建立： 由Gambit 2.3所建構而成 網格(Mesh)建立： 均勻網格 (Uniform Mesh) (Elements：Hex/Wedge，Type：Cooper) 網格數目：控制在10萬左右

  7. 網格獨立性測試 • 以寬深比AR=10(流道寬度:500μm) • 無壓力振盪 • 邊界條件: 入口流速Vin=0.028 m/s 出口壓力 Pout=0 • 流體性質: μ=1.03e-3 kg/m．s D=1e-8 m2/s

  8. 不同網格數出口斷面中心線速度分佈

  9. 數值方法—邊界條件 入口條件： 均勻的入口速度，左右入口設定流體A與流體B。 出口條件： 設定為大氣壓力，即錶壓力為0 邊界條件： 壁面設定為無滑動條件(no slip condition) 左右振盪流道之振盪壓力設定為：

  10. 濃度場之邊界條件 • 左側流道入口與左側振盪流道之濃度0，右側入口流道與右側振盪流道之濃度為1 • 壁面為濃度通量為0(zero diffusive flux) 即 ， 為壁面之單位法向量。

  11. 混合器流道結構示意圖 邊界條件 Vin=constant AR=b/h Pout=0

  12. 數值模擬步驟 使用軟體: Fluent 6.3 Gambit 2.3 Tecplot 360

  13. 模擬參數 流體性質:黏滯係數 μ=1.03e-3 kg/m．s 擴散係數 D=1e-8 m2/s 幾何尺寸:出口長度:18.25mm 微流道斷面:2.5×10-2 mm2 流道寬深比AR(b/h):1.6~10 入口流速範圍:0.028~0.084 m/s (Re=2.94~8.82) 干擾壓力振盪的相位差範圍:0~π 振盪頻率: 0~50Hz

  14. C=0.5 C=1 C=0 量化混合程度-定義混合效率 出口處斷面之混合效率公式： (Xia et al., 2005) Ci：格點i上的濃度 C∞：完全混合的濃度 n：格點數

  15. 無壓力振盪干擾時的混合情形 不同入口流速下,不同寬深比 出口處數值混合效率的比較 AR=6.4, 不同入口流速，出口 剖面之濃度分佈截面 誤差0.5～5 %

  16. 壓力振盪干擾=0時的混合情形 AR= 6.4，Vin=0.056 m/s，出口處混合效率 隨時間變化的情形。

  17. 流場與濃度場分析 • 流道中間剖面及其混合流道在各位置之流線圖。 • 出口及中央剖面之濃度分佈圖。 AR=6.4，Vin=0.056 m/s，=0，在頻率為10 Hz，t=1.0 s時

  18. 流線圖

  19. 濃度分佈圖 (a)出口剖面 (b)中央剖面

  20. 壓力振盪干擾=0～π之混合效益 AR=6.4，Vin=0.056 m/s，=π，在不同之振盪頻率下，混合效率隨時間變化的情形。 AR=6.4，Vin=0.056 m/s， =0～π混合效率隨頻率變化情形。

  21. 流場與濃度場分析(2) AR=6.4，Vin=0.056 m/s，=π，在頻率 為10 Hz，t=1.0 s時 • 流道中間剖面及其混合流道在各位置之流線圖。 • 出口及中央剖面之濃度分佈圖。

  22. 流線圖(2)

  23. 濃度分佈圖(2) (a) 出口剖面 (b)中央剖面

  24. 入口流速對混合效益之影響 AR=6.4，=π，不同入口流速下之混合效率隨頻率變化的情形。

  25. 流道寬深比對混合效益的影響 Vin=0.056 m/s，=π時，不同寬深比AR， 混合效率隨頻率變化的情形。

  26. 流場與濃度場分析(3) AR=1.6，Vin=0.056 m/s，=π，在頻率 為10 Hz，t=1.0 s時 • 流道中間剖面及其混合流道在各位置之流線圖。 • 出口及中央剖面之濃度分佈圖。

  27. 流線圖(3)

  28. 濃度分佈圖(3) (a) 出口剖面 (b)中央剖面

  29. AR=1.6，=π，入口流速對混合效益之影響 。 η=96%

  30. AR=1.6，Vin=0.028 m/s，=π，在20 Hz，t=1.0 s時之濃度分佈。 η=96% (a) 出口剖面 (b)中央剖面

  31. 結論 • 利用兩側振盪流道加入壓力干擾，確實可提高混合效率約2～5倍 • 兩側振盪流道導入無相位差的壓力干擾，可使混合效率提高約2倍，且混合效率並不會隨頻率變化。 • 兩側導入具有相位差之壓力干擾，則會有較佳的混合效果，並且混合效率會隨振動頻率先升高達到最佳值後會隨頻率升高而下降的趨勢 。 • 在=π時，流體界面會有最大彎曲，偏折的現象而具有最大之接觸面積，在以擴散主宰的低雷諾數流動範圍，會有最佳的混合效果。

  32. 4. 流道寬比深AR愈低，混合效率愈好。在AR=1.6時分別在軸向及橫向剖面都會產生渦流，使流體界面有摺疊(folding)，扭曲(distortion)的現象，更能幫助兩流體的混合。 5. 在較低的入口流速，有較長的擴散時間，可提高混合的效率。 6. 本研究中當流道寬深比AR=1.6，入口流速Vin=0.028 m/s，兩側壓力振盪=π，兩流體的混合效率可達到η=96%。

  33. 報告完畢 謝謝大家

  34. 不同網格數出口斷面混合效率的比較

  35. 混合流道不同位置剖面流場觀察區

More Related