等位函数法应用在两液滴连续撞击壁面之三维数值计算.PDF

工程應用技術學刊 第二卷 第一期 等位函數法應用在兩液滴連續撞擊壁面之三維數值計算 劉昭忠 1* 劉源昌 2 余森年 3 1,2,3中州科技大學 電機與能源科技系 摘 要 本研究以三維數值方法計算兩液滴連續撞擊壁面的暫態過程。數值方法上以結構性多區 塊、交錯網格系統進行空間離散，並採用有限體積法來依序疊代求解連續、動量等守恆方程 式，而流場中的速度與壓力間藕合以 SIMPLER 處理。在兩液滴連續撞擊壁面過程中，所產 生的合併及擴展現象，乃採用等位函數法加以處理之。為了驗證三維數值程式的準確性，先 計算雷諾數 400及韋伯數 4.3 純水滴撞擊平板的暫態過程，並與實驗結果相比較，發現計算 結果非常準確。隨後，探討不同雷諾數的液滴 2 與固定雷諾數的液滴 1 兩者連續撞擊壁面的 過程，從數值計算過程中，不但可瞭解兩液滴撞擊的動態模式外，還能進一步獲知不同雷諾 數的液滴 2 對撞擊後兩液滴合併成單一液滴最大擴展距離之影響。 關鍵字：液滴、等位函數法、擴展距離 通訊作者 姓名 ：劉昭忠 E-mail ：c00jcl00@dragon.ccut.edu.tw 27 Journal of Applied Engineering Technology Vol.2 No.1 壹、緒論 自由液面乃是一移動邊界，其流場範圍或邊界形狀會隨時間變化，因此，不論從理論或 數值計算均很難處理。有鑑於此，遂利用 Osher 和 Sethian （1988,1994）提出一套無須簿記 （Bookkeeping ）而能追蹤移動邊界的準則，亦即Hamilton-Jacobi 等位函數（Level Set Function ） 表示法，此法可輕易解決邊界合併、破裂等拓樸性質改變的問題，推展至三維計算亦相當直 接，故發展出等位函數法數值模擬程式，且應用於噴霧系統中液滴撞擊模式的研究。 在真實噴霧系統中，許多霧化的細小液滴會連續撞擊壁面，因此，從過去迄今，已有許 多文獻加以研究（1983, 1991, 1993, 2000a, 2000b, 2007a, 2007b）。由於霧化的液滴是連續撞擊 壁面，而液滴連續撞擊壁面之行為模式，到目前為止，則是較缺乏清楚詳細的探討，本研究 以三維數值方法計算兩液滴連續撞擊壁面的暫態過程，以便對此撞擊模式有更深入的瞭解。 以數值方法模擬液滴連續撞擊壁面，是一個非常複雜且困難的問題。因為在極短時間內， 由於液滴的撞擊作用，使得液滴彼此間與壁面會產生一連串的變化，所以如何捕捉液滴自由 表面（Free Surface ）的變化，將是數值模擬的重點。 貳、文獻探討 Madejski （1983）為了克服液態金屬液滴撞擊平板的擴展過程，乃利用慣性力、黏滯力 及表面張力三者達平衡時的量階（ Order of Magnitude ）分析，發展出近似方法來捕捉液滴變 形擴展過程。Trapaga和 Szekely （1991）則使用數值模式，配合標誌法和流體體積法（Volume of Fluid ）的方式，來模擬平板上液滴的變形過程，並提出相關變數對此過程之影響。Fuki et al. （1993）使用有限元素法來模擬液滴垂直撞擊平板的過程。為了準確求得液滴複雜的變形行 為，配合特別的技巧來修正液滴變形後的計算網格，並提出雷諾數（Reynolds Number ）、韋 伯數（ Weber Number ）對油滴擴展的影響。 最近， Siddhartha et al. （2007a ）使用VOF 方法和利用靜態接觸角(static contact angle: SCA)及動態接觸角 (stat