行政院國家科學委員會專題研究計畫 成果報告
聲波引致懸浮微粒聚結的理論研究及數值計算(2/2)
計畫類別： 個別型計畫 計畫編號： NSC91-2212-E-002-074- 執行期間： 91 年 08 月 01 日至 92 年 07 月 31 日 執行單位： 國立臺灣大學應用力學研究所 計畫主持人： 李雨 計畫參與人員： 劉義賢，羅凱升 報告類型： 完整報告 處理方式： 本計畫可公開查詢中 華 民 國 92 年 10 月 31 日
行政院國家科學委員會補助專題研究計畫
Ç
成 果 報 告
□期中進度報告
聲波引致懸浮微粒聚結的理論研究及數值計算(2/2)
Theoretical and numerical study of acoustic coagulation of suspended particles
計畫類別：
Ç
個別型計畫 □ 整合型計畫
計畫編號：NSC 91－2212－E－002－074
執行期間： 91 年 8 月 1 日至 92 年 7 月 31 日
計畫主持人：
李雨
共同主持人：
計畫參與人員：
劉義賢、羅凱升
成果報告類型(依經費核定清單規定繳交)：□精簡報告
Ç
完整報告
本成果報告包括以下應繳交之附件：
□赴國外出差或研習心得報告一份
□赴大陸地區出差或研習心得報告一份
□出席國際學術會議心得報告及發表之論文各一份
□國際合作研究計畫國外研究報告書一份
處理方式：除產學合作研究計畫、提升產業技術及人才培育研究計畫、
列管計畫及下列情形者外，得立即公開查詢
□涉及專利或其他智慧財產權，□一年□二年後可公開查詢
執行單位：
國立臺灣大學應用力學研究所
中 華 民 國 92 年 10 月 30 日
中文摘要
本計劃延續上年度計劃，以數值計算方法研究二相懸浮流中微粒的聲波聚結問題。本 文完成之工作有二：(1)就旅波(traveling wave)中微粒聚結的數值計算方面，本文解決了上 年度計劃中採用流力作用(hydrodynamic interaction)模式及聲波尾跡(acoustic wake)模式計 算時質量不守恆的問題。雖然本文的計算結果較前人者更接近實驗結果，但如同上年度計 劃中所採用的正動力作用(orthokinetic interaction)模式的模擬結果，我們需要採用一經驗的 碰撞效率係數方可使數值計算與實驗結果相吻合。這顯示目前文獻中現有的碰撞頻率函數 理論模式並不完備，有進一步研究的必要。(2)就駐波(traveling wave)中微粒聚結的數值計 算方面，本文對微粒顆數濃度隨體積分佈的廣義動力方程式，提出一項可對總顆數濃度作 分析模擬的簡化理論計算模式；假設微粒顆數濃度隨體積分佈呈對數常態分佈(log-normal distribution)，則我們可將廣義動力方程式此一微積分方程式(integral-differential equation)簡 化為一偏微分方程式，而可大幅降低數值模擬所需的時間。此一簡化理論計算模式所獲結 果與實驗結果吻合，而可應用於聲波集塵器的參數分析及設計。 關鍵詞：聲波微粒聚結，旅波與駐波，理論分析，數值模擬。 Abstract
This second year project studies numerically the acoustic coagulation of small particles in a two-phase suspension flow. Two topics have been completed in the present work. (1) For the acoustic coagulation of particles in the traveling wave, we have solved the false mass conservation problem for the simulation using both the hydrodynamic interaction and the acoustic wake model in the last year. Although the agreement between the present calculations and the experiments is better than the calculations in the literatures, we still need to incorporate an empirical collision efficiency factor in the calculation, as that in the last year using orthokinetic interaction model. This indicates that the available models for the collision frequency function in the literatures have certain drawbacks, and further theoretical analyses are required. (2) For the acoustic coagulation of particles in the standing wave, we proposed a simplified theoretical model for the general dynamic equation of the continuous distribution function of particles. The log-normal distribution of particles is assumed, and a resulting partial differential equation is obtained from the integral-differential general dynamic equation. Numerical results of such partial differential equation agree nicely with the experiments. Thus the proposed simplified model can be applied for the parameter analyses and the design of the acoustic agglomerator.
Keywords：acoustic coagulation of small particles, traveling and standing wave, theoretical analysis, numerical calculation.
前言 在環境保護及工業應用上，常需在流體中移除微米及次微米微粒。但此等微粒因體積 細小，慣性效應不彰；而其中0.1-1μm 範圍的較大微粒，擴散機制也相對薄弱，因此要在 流體中移除微米及次微米微粒實有先天上的困難。在眾多前人所提出的方法中，有一項有 趣而有效的方法，乃是利用聲波在二相懸浮流中傳播，引致微粒聚結成較大顆微粒，再用 傳統方法將聚結後的大顆微粒移除。早期相關回顧性文獻可參考Mednikov (1965)；近日者 則可參考Hoffmann (1993)及黃家增(2000)。 研究目的 本文延續去年計劃(見李雨，2002)，主要工作有二：一為繼續對旅波中的聚結模擬進行 研究，以了解目前文獻中現有理論模式的可用性。另一則為對駐波中的聚結發展一項簡化 的理論模式，並進行數值研究。後者成果可直接應用到聲波聚結器的設計工作上。 物理模式與數值計算 (i)旅波中的分析 本項工作的相關理論推導、數值計算方法、及詳細結果可參考羅凱升(2003)。在旅波中 流體速度只為時間的函數，令 為體積為 v 的微粒數目濃度，其隨時間 t 的演變滿足聚 結方程式(見 Friedlander, 1977) ( ) n v 0 ( ; ) 1 ( , ) ( ; ) ( ; ) 2 v dn v t v v v n v t n v v t dv dt =
∫
β − − n v t( ; )0β( , ) ( ; )v v n v t dv (1) ∞ −∫
其中β( , )v v 代表體積為 v 的微粒和體積為 的微粒的碰撞頻率函數，v β( ,v v −v 代表體積為) v 以外的懸浮微粒和其他”適當”體積的微粒碰撞、而產生體積為 v 的懸浮微粒的碰撞頻率函 數。碰撞頻率函數方面，視不同機制有不同的表示式。就布朗擴散、正動力作用(orthokinetic interaction)、流力作用(hydrodynamic interaction)及聲波尾跡等模式(acoustic wake model)等的 碰撞頻率函數β β βBr, ok, hi及βho，我們在上年度報告中已加以詳述(請參考李雨，2002)。在 本節中，我們將取以下兩種狀況進行模擬：β β= song ≡Eokβok+Ehiβhi (1a)
及 β β= hoffman ≡Ehoβho, (1b) 其中 及 為相對應模式的碰撞效率係數，在本文中為可調整的常數。其中(1a)式為
Song et al. (1994)所提出者，Ezekoye and Wibowo (1999)曾就 , ok hi E E Eho 0.5 ok hi E =E = 進行數值模擬。 聲波尾跡模式為Hoffman (1997)所提出者，他採用再填充機制(refill mechanism)將正動力作 用及流力作用以較含乎物理的方式結合，而非(1a)式的線性組合。在(1a)及(1b)式中， 2 2 0 2 | | ( , ) ( ) 2( ) 1 ( ) 1 ( ) v w ok ok v w vw v w v w U v w r r u r r 2 ω τ τ β β π ωτ ωτ − = = + = + + + , (2a) ( , ) ( )2 , (2b) hi v w rv rw Q β =β =π + 及 0 2 12 ( ) ( , ) 1 1 2 v v w w c ho ho ok U r l r l T u w L β β β π ⎛ + ⎞ = = ⎜ + ⎟ ⎝ − , (2c) ⎠ 其中 及rv rw分別為體積為 v 及 w 微粒的半徑，τv及τw分別為體積為 v 及 w 微粒的鬆馳時間 (relaxation time)，U0及ω 分別為旅波的速度振幅和角頻率；在史篤氏阻力的條件下，鬆馳
) 時間 2 2/(9 p p pr g τ = ρ ρη ，下標 p 可代表(2a)式中的 v 或 w，η為載運流體的運動黏滯係數， p ρ 及ρg分別為微粒及流體的密度。又(2b)式中，按 Dianov et al. (1968)的分析，
(
2 2 2 2)
2 2 0 0 0 2 2 2 3 6 ( ) ( 2 v v w w v v w w v w 1 v v w w d d U U U Q r l r l r l r l r r r l r l L L η π πη π π ⎡ ⎤ = ⎢ + + + ⎥− + − + ⎣ ⎦ 9 ) 6 Ld 2 0 2 3 ( )( 8 d v w v w w v v w U l l l q l q r r L ω + − − ) 其中 , 1 p g2 g p h l µ µ + = 2 2 1 p g g p p p h h q µ µ µ + + = , , 2 9 0 p p p p g p R r U h πωτ ω πρ ρ = = 0 p p U r R η = , , 1 1 2 2 P p τ ω µ + = 2 2 1 p p g τ ω ωτ µ + = , 及 為兩微粒間的距離，在本文中取 ，其中 N 為微粒的瞬時總顆數濃度。上述 各式中下標 p 可代表(2a)至(2c)式中的下標 v 或 w。另(2c)式中， 為聲波週期，L 為聚結體 積的長度，定義為兩倍微粒 v 和微粒 w 間相對速率(按正動力作用)和 的乘積。 d L Ld =N−1/ 3 c T c T (1)式中所示的微-積分方程式並無分析解，本研究中我們基本上採 Temkin(1994)的方法 進行計算，但在計算過程中加入了適當的判斷條件，並不對尚未聚結出的微粒種類進行計 算，因而大幅節省了計算時間(見黃國榮，2002)。在分析過程中，首先我們將(1)式無因次 化，再將無因次方程式的積分項離散化，也就是將(1)式所示的連續微粒顆數濃度分佈改採 離散形式分佈，最後可得一組非線性常微分聯立方程式組，在給定始起狀態的微粒濃度分 佈下，我們採用四階朗次-庫他(Runge-Kutta) 算則進行求解。在本年度作中，我們修正了去 年的一項程式錯誤，使模擬過程中的微粒總體積守恆。 (ii)駐波中的分析 在駐波中流體速度為時間及空間的函數，微粒數目濃度n v x t( ; , )滿足(
)
(
)
0 ( ; , ) 1 ( ; , ), ( ; , ) ( , ) ( ; , ) ( ; , ) 2 v n v x t n v x t D n v x t v v v n v x t n v v x t dv t β ∂ + ∇ ⋅ = ∇ ⋅ ∇ + − − ∂ V∫
)
⋅ , (3)(
0 ( ; , ) ( , ) ( ; , ) ( ; , ) , n v x t β v v n v x t dv n v x t ∞ −∫
− ∇ c 如將駐波前進方向定為x 軸，i 為其上的單位向量，則流場速度 0 2 sin sin u U kx ωt = = V i i (3a) 飄移速度(drifting velocity)按 Grinshpun et al. (2000)的分析可表示為(
)
0 2 sin 2 2 4 1 kx c U ε α , α = − + c i = i (3b) 擴散係數按Friedlander (1997)第二章為 1 2 1 2exp 3 3 3 KT KT l d D C A A A d d d l πµ πµ , ⎧ ⎡ ⎛ ⎞⎤⎫ = = ⎨ + ⎢ + ⎜− ⎟⎥⎬ ⎝ ⎠ ⎣ ⎦ ⎩ ⎭ (3c) 其中2U k0, 及ω分別為駐波的速度振幅、波數及角頻率，K 為波士曼常數，T 為絕對温度，(
)
0 0 2U /Cs 2U / /k ε = = ω 為速度振幅與聲速C 的比值，s α ωτ= ，而τ ρ= pd2/(18 )µ 為微粒 的鬆馳時間，µ為載運流體的黏滯係數，ρp為微粒的密度，l 為氣體分子的平均自由徑；而 碰撞頻率函數β則可採用(2a)-(2c)所述各項模式。(3)式為一統御n v x t( ; , )的微積分方程式，如對 v 此一自變數離散化，可由(3)式導出一組聯立偏微分方程式組，但因為微粒聚結的關 係，n 隨 v 變化的分佈範圍會隨著時間的推演而快速增加，也就是聯立的偏微分方程式數目 會快達增加(一般案例可達數萬條)，而非目前一般電腦所能處理；再者，文獻中的實驗結果 只提供平均粒徑隨時間的演變(見 Riera & Gallego, 1986)。因此本研究將不對(2)式直接進行 模擬，而是配合適當的假設，由(2)式推導出總數目濃度 0 ( , ) ( ; , ) , N x t n v x t dv ∞ = ∫ (4) 的統御方程式，再對所得的一條偏微分方程式進行數值求解。對(3)式就 v 變數積分，並利 用(4)式的定義、及(3a)-(3c)式，可得
( )
2 2 0 0 0 1 ( ; , ) ( , ; , ) ( ; , ) ( ; , ) 2 v uN N Dn v x t dv v v v x t n v x t n v v x t dvdv t x x β ∞ ∞ ∂ ∂ + = ∂ + − − ∂ ∂ ∂∫
∫ ∫
(
)
0 0 0 ( , ) ( ; , ) ( ; , )v v n v x t n v x t dvdv cn v x t dv; , , x β ∞ ∞ ∂ ∞ − − ∂∫ ∫
∫
(5) 因為本文所欲模擬之實驗案例其初始粒徑為0.1µm 的量階，布朗擴散引致的微聚結效應相對 重要，因此(5)式中的碰撞頻率函數我們取 β =Ehoβho +βbr, (5a) 其中βho如(2c)式所示，我們採用 Hoffman(1997)的聲波尾跡模式來考慮聲波聚結，Eho為碰 撞效率，為一可調參數；βbr為布朗擴散引致的微聚結效應，按Friedlander (1977)，可表示 為(
)
1/ 3 1/ 3(
1/ 3 1/ 32
1
1
,
3
br i j i j i jkT
v v
v
v
C
v
v
β
µ
⎛
⎞
=
⎜
⎜
+
⎟
⎟
+
⎝
⎠
)
,
(5b) 其中參數 C 如(3c)式所示。按文獻中的實驗結果，n 隨 v 的變化大約呈對數常態分佈，因此我 們假設(
)
2 2 ln / ( ; , ) 1 exp ( , ) 2 (3ln ) 18ln g v v v v n v x t N x t v π σ σ ⎡ ⎤ = ⎢− ⎥ ⎢ ⎥ ⎣ ⎦ , (6)其中v 為微粒的幾何平均體積g (geometry mean volume)，σv為標準偏差(standard deviation)。
利用(6)式，我們可處理(5)式中的積分項，而導得一條統御 N(x,t) 的偏微分方程式，可以無 因次方式表示如下：
(
)
2(
)
2 2 2 1 1 2 2 3 4 sin sin 2 sin 2 N x N x N N I t A I A I N I t ε x x ε ∂ ∂ ∂ + = ∂ + ⎛ − ⎞ + ⎜ ⎟ ∂ ∂ ∂ ⎝ ⎠ ∂x (7) 其中 0 0 0 , , , , , 0 , g g v N v v N x kx t t v v v N ω v v = = = = = = v 0 N ,d0及v0分別為微粒的初始總數目濃度、最小直徑及最小體積， 2 1 0 3 k KT A d ω πµ = 及 2 0 0 0 2 2 4 N U d A ω = 為無因次參數， 1 1( , ,g v), 2 2( ,0 , ,g v), 3 3( ,0 , ,g v), 4 4( ,0 , ,g ), I = I Kn v σ I =I α Kn v σ I =I α Kn v σ I =I α Kn v σv 為含相關無因次參數α0,Kn v, g及σv的定積分，Kn=d0/l為富森數(Knudsen number)，及 2 0 0 pd0 /(18 ) α =ωτ =ωρ µ 。我們將利用IMSL 副程式來計算I I I 及1, ,2 3 I ，利用有限差分法解4 (7)式，再考慮微粒的總質量守恆，以所算出的瞬時總數目濃度來計算瞬時的微粒幾何平均 體積及直徑。配合實驗給予適當的初始條件，及採用週期性邊界條件，我們將在本文中模擬Riera & Gallego (1986)的結果，以驗證我們上述所發展的模式。 結果與討論 (i)旅波中的狀况 在本文的計算中，初始懸浮微粒均採對數常態分佈，其總顆數濃度 ，幾何算術平均 體積 0 N g V ，標準差 σ。在上年度計劃中，我們對正動力作用的情況下進行了詳細模擬。在本 年度工作中，我們採(1a)及(1b)進行計算，並與前人實驗結果進行相互驗證。在圖一及二中， 我們分別以(1a)及(1b)式模擬 Brandt et al. (1936)的實驗，結果顯示我們可籍調整碰撞效率係 數來使計算與實驗結果相吻合，且在模擬過程中可保持微粒總體積守恆。在圖三及四中， 我們將計算所得與Magill et al. (1991)的實驗作比較，結果顯示因參數不同，我們需要選擇 不同於圖一及二的碰撞效率係數方可使計算與實驗大致相符；此項結果雖屬差強人意，但 還是優於Ezekoye & Wibowo (1999)的計算(即圖二中的 MAEROS simulation)。雖然我們進 行了更多的模擬，但因參數眾多，我們並不能歸納出 及 隨各參數變化的關係；事 實上目前碰撞頻率函數的物理模式並不完備，有待進一步作深入的理論探討。雖然聲波聚 結數值模擬的定量結果的正確性須引入可調經驗常數 或 ，但對聲波聚結各操控參 數的定性行為，我們仍可藉着數值模擬來獲致全面的了解，詳細結果可參考羅凱升(2003)。 , ok hi E E Eho , ok hi E E Eho (ii)駐波中的狀况
Riera & Gallego (1986)的實驗在目前現有文獻中相當完整，其相關實驗參數及結果如表 一所示，本文將針對該表中所示各案例進行模擬，以驗證本計劃中我們所發展的理論模式。 在圖五、六及七所顯示的不同案例中，如選取適當的 值(見各圖的說明)，本文的模擬結 果與實驗結果相符；且所耗用電腦計算量不大，因此本計劃所提出的理論模式可作為聲波 集塵器的一項初步不錯的研發及設計工具。圖五(a)、六及七乃比照 Riera & Gallego (1986) 的方式，以瞬時平均粒徑與初始粒徑比值 G 隨時間的演變來表示聲波聚結的效率，如改以 總微粒顆數來表示，微粒總顆數在聲波的作用下快速下降；如以圖五(b)為例，總顆數在一 秒後即下降四個量階，在二秒後下降至五個量階，在四秒後則下降六個量階。也就是說一 般而言，最重要的聚結發生在聲波作用的初期，這是因為聚結後微粒顆數濃度迅速下降、 而造成後續聚結效應變弱，所以理論模式首重聚結初期的正確性。圖六顯示聚結效果隨聲 波強度的增加而增加，圖七顯示聚結效果隨初始粒徑的下降而上升，後者之原因與初始總 顆數濃度隨初始粒徑的下降而上升(見表一)有密切的關係；因顆數濃度愈大，微粒碰撞聚結 的機率也愈大。圖八顯示聚結效果隨初始顆數濃度的增加而快速增加。更詳細的結果可參 考劉義賢(2003)。 ho E 參考文獻
1. Brandt, O., Freund, H., and Hiedemann, E., Kolloid zeitschrift, Vol. 77, pp. 103-115, 1936. 2. Dianov, D. B., A. A. Podol’skii and V. I. Turubarov, Soviet Physics – Acoustics, V.13, No.3,
pp.314-319 (1968).
3. Ezekoye, O. A. and Y. W. Wibowo, J. Aerosol Sci., V.30, No.6, 1117-1138, (1999).
4. Friedlander, S. K., Smoke, dust and haze – Fundamentals of aerosol behavior, John Wiley & Sons (1977).
5. Hoffmann, T. L., Visualization of particle interaction and agglomeration in an acoustic field, PhD Dissertation, Pennsylvania State Univ. (1993).
6. Hoffmann, T. L., J. Aerosol Sci., 28, 919-936 (1997).
7. Grinshpun S. A., Y. N. Redcoborody, S. G. Kravchuk, V. I. Zadorozhnii and V. I. Zhdanow, International J. Multiphase Flow, 26, 1305-1324 (2000).
8. Magill, J. S., Ph. Caperan, J. Somers, K. Richter, G. Rodriguze-Corral, E. Riera-Franco de Sarabia and J. A. Gallego-Juarez, J. Aerosol Sci., 22, Suppl.I, pp.S27-S30, (1991)
9. Mednikov, E. P., Acoustic coagulation and precipitation of aerosols, Consultants Bureau, New York (1965).
10. Riera-Franco De Sarabia, E. and J. A. Gallego-Juarez, J. Sound Vib., 110, 413-427 (1986). 11. Song, L., G. H. Koopmann and T. L. Hoffmann, ASME J. Vib. Acoustics, 116, 208-214
(1994).
12. Temkin, S., Phys. Fluids, 6, 2294 -2303 (1994).
13. 李雨, “聲波引致懸浮微粒聚結的理論研究及數值計算(1/2),” 國科會研究計劃報告, NSC 90-2212-E-002-199 (2002) 14. 黃家增, “聲波引致氣體中懸浮微粒的聚集研究,” 國立臺灣大學應用力學研究所碩士論 文 (2000). 15. 黃國榮, “聲波及布朗擴散對次微米微粒聚結效應的研究,” 國立臺灣大學應用力學研究 所碩士論文(2002). 16. 羅凱升, “兩項微粒聲波聚集模式的數值研究,” 國立臺灣大學應用力學研究所碩士論文 (2003). 17. 劉義賢, “駐波聲場中次微米微粒的聚集研究,” 國立臺灣大學應用力學研究所碩士論文 (2003). (b) (a) 圖一 初始懸浮微粒為對數常態分佈，N0=4.64×1012顆/m3，d0=0.766µm，ε=10-4 (誤差参數)， Vg=3.59×10-18m3，σ=1.2，聲波強度 138dB，聲波頻率 10KHz；在(1a)式所示之正動力加流 力作用下， (a)不同碰撞效率之無因次總顆數濃度n / N0隨時間的演變，(b)無因次微粒總體 積V/V0隨時間的演變。N0 及V0 為初始的總顆數濃度及微粒總體積。 (b) (a) 圖二 模擬條件與圖一同，但採用(1b)式所示之聲波尾跡效模式進行計算。
(b) (a) 圖三 初始懸浮微粒為對數常態分佈，N0=1012顆/m3，d0=0.323µm，ε=10-5，Vg=2.68×10-19m3， σ=1.35，聲波強度 155dB，聲波頻率 9.4KHz，ρg=1115kg/m3；在(1a)式所示之正動力加流力 模式作用下，(a)不同碰撞效率之無因次總顆數濃度n / N0隨時間的演變，(b)無因次微粒總體 積V/V0隨時間的演變。N0 及V0 為初始的總顆數濃度及微粒總體積。 (a) 圖四 模擬條件與圖三同，但採用(1b)式所示之聲波尾跡效模式進行計算。 (b) (b) (a)
圖五(a) 表一中 CASE 1 數值計算與實驗結果的比較。(a) 瞬時平均粒徑與初始粒徑比值 G 隨時間的演變，可用以表示聲波聚結效果的大小。(b)無因次總顆數濃度 隨時間的演 變。其中 為瞬時的總顆數濃度， 為初始總顆數濃度。本例中 ，聲波頻率 20.4KHz，。 0 / N N N N0
E
ho=
0.1
圖六 表一中 CASE 8、CASE 9、CASE10 數值計算與實驗結果的比較。本圖顯示不同聲波強 度的效應。本例中
E
ho=
0.15
，聲波頻率 20.4KHz。圖七 表一中 CASE 2、CASE 7、CASE12 數值計算與實驗結果的比較。本圖顯示不同初始粒 徑的效應。本例中聲波頻率 20.4KHz，
E
ho=
0.1(R0 =0.1µm)、0.15(R0 =0.22µm)、 0.45(R0 =0.43µm)。圖八 初始總顆數濃度N0對聚合趨勢的影響。本例中微粒初始平均半徑R0 =0.22µm，標準 差(σ )1.68，聲波頻率 20.4KHz、及聲波強度( I )160dB。
表一 Riera & Gallego (1986)實驗案例的參數與結果。他們以三種不同初始微粒數目濃度及 五種不同聲波強度組合成 15 組實驗。 0( 為初始微粒平均半徑，vg( )為初始平均體 積，c( )初始單位體積微粒重量， ( )初始微粒總數目濃度， 3 m 3 / g m N0 顆/ m3 σ 為分佈的標準平 均差，I( )為聲波強度，dB 為以分貝表示的聲波強度，G 為時間 0 至 5 秒瞬時平均微 粒粒徑和初始平均粒徑的比值，其變化表示聲波聚集的趨勢。 2 / W cm ) R µm 計劃成果自評 本報告內容與原計劃書所述相符，達成預期目標。相關成果具一定的學術價值，現正 撰寫論文，擬投稿於學術期刊(Journal of aerosol science)及會議(中華民國燃燒學會第十四屆 年會)。所發展的理論模式可助研發聲波集塵器，對 IGCC 等相關先進燃煤發電系統的研發 可能會有所貢獻。
