關(guān)善超,胡雪佳 ( 陜西師范大學(xué)物理學(xué)與信息技術(shù)學(xué)院，陜西 西安 710062)

沈壯志 (陜西師范大學(xué)物理學(xué)與信息技術(shù)學(xué)院 陜西省物理研究所,陜西省超聲學(xué)重點實驗室，陜西 西安 710062)

自激共振噴嘴與錐形噴嘴空化發(fā)生器的空化特性仿真研究

關(guān)善超,胡雪佳 ( 陜西師范大學(xué)物理學(xué)與信息技術(shù)學(xué)院，陜西 西安 710062)

沈壯志 (陜西師范大學(xué)物理學(xué)與信息技術(shù)學(xué)院 陜西省物理研究所,陜西省超聲學(xué)重點實驗室，陜西 西安 710062)

運用Fluent軟件、k-ε模型、氣液兩相流空化模型，模擬計算了相同噴嘴直徑下的自激共振噴嘴和錐形噴嘴空化發(fā)生器的空化特性，得出了自激共振噴嘴在同樣條件下出現(xiàn)空化現(xiàn)象所需時間比錐形噴嘴少，空化區(qū)域大，湍流動能強的結(jié)論，對設(shè)計空化發(fā)生器有一定的指導(dǎo)作用。

自激振動；空化；噴嘴

隨著科技的進步，生活中的污染尤其是水污染也越來越嚴重。我國是世界上最缺水的國家之一，怎樣方便且無二次污染的處理化學(xué)生活污水成了重要的研究課題，而空化現(xiàn)象由于其瞬間產(chǎn)生的高溫高壓能很好的分解大分子化合物且無二次污染，成為了人們重點關(guān)注的方法之一。產(chǎn)生空化的基本方法有超聲空化和水力空化，超聲空化已經(jīng)在加速化學(xué)反應(yīng)實現(xiàn)非均相反應(yīng)物間的均勻混合及超聲清洗、鉆孔方面做出了不少貢獻，但是超聲空化同時存在能量利用率低、空化效應(yīng)區(qū)域小的缺點[1]，而水力空化正好彌補了該缺點，在化工應(yīng)用方面優(yōu)勢明顯[2]。

近些年研究表明，空化發(fā)生器的結(jié)構(gòu)對空化場的分布有十分顯著的影響，漸變型空化發(fā)生器產(chǎn)生的空化區(qū)域范圍大，比突變性反應(yīng)發(fā)生器優(yōu)勢明顯，且最佳入口壓力為0.55MPa[3]。自激共振射流可以很好地產(chǎn)生空化效應(yīng)，且湍流動能更強[4]。王萍輝等的研究表明，諧振腔長度L與諧振腔噴嘴直徑d的比值L/d不同會影響出口速度及脈動壓力，而當(dāng)其比值在2.1附近時出口速度最大[5]。下面，筆者研究自激共振噴嘴與錐形噴嘴(自激共振噴嘴選用L/d=3的自激共振噴嘴，k-ε湍流方程和schnerr-sauer空化模型[6])在0.55MPa下的空化場及壓力場。

1 2種噴嘴模型的建立

圖1和圖2分別是錐形噴嘴空化發(fā)生器和自激共振噴嘴空化發(fā)生器2種噴嘴的平面示意圖，示意圖上所有尺寸的單位為毫米，D為入口直徑，d為噴嘴直徑。L為諧振腔長度，保證流場區(qū)域長度和尺寸一致。

圖1 錐形噴嘴空化發(fā)生器 圖2 自激共振噴嘴空化發(fā)生器

2 數(shù)值模擬與控制方程

Fluent模擬流體空化使用Gambit建立二維模型和劃分網(wǎng)格，網(wǎng)格數(shù)在10萬左右，充分保證了模擬的精確性。經(jīng)驗證即使網(wǎng)格數(shù)再增加，對模擬精確度影響很小而且影響計算速度。流域內(nèi)流體為水，其參數(shù)為fluent數(shù)據(jù)庫自帶的標(biāo)準(zhǔn)大氣壓下的水的參數(shù)。

1)模擬的動力學(xué)方程 流體為不可壓縮粘性流體，其運動規(guī)律受控于Navier-Stokes方程：

式中，u為流體流速；P為流體壓力；ρ為流體密度； ?為流體粘滯系數(shù)。

湍流方程采用標(biāo)準(zhǔn)k-ε模型，其中k表示湍流動能，ε表示湍流耗散項。其運輸方程為：

式中，Gk表示由層流速度梯度而產(chǎn)生的湍流動能；Gb是由浮力產(chǎn)生的湍流動能；YM是可壓縮湍流中過渡的擴散產(chǎn)生的波動；C1、C2、C3是常量；σk和σε是k方程和ε方程的湍流Prandtl數(shù)；Sk和Sε自己定義。Fluent默認值為：

C1ε=1.44C2ε=1.92Cμ=0.09σk=1.0σε=1.3

2)模擬的空化模型 采用Schnerr-Sauer空化模型，主相為水，第2相為水蒸氣，由于流體的高速運動，造成局部出現(xiàn)負壓區(qū)，當(dāng)負壓達到一定值(-9.78×104Pa)，則fluent默認出現(xiàn)空化現(xiàn)象?？栈鲌鍪芤韵路匠炭刂疲?/p>

式中，t表示時間；u表示速度；ρm、ρ?、ρ表示混合相密度、汽相密度和液相密度；δij表示克羅內(nèi)克數(shù)；α?表示汽相體積分數(shù)；μ和μt表示混合介質(zhì)動力粘度和湍流粘度；R表示相間質(zhì)量傳輸率。

3 數(shù)值模擬結(jié)果與分析

1)錐形噴嘴的壓力場和空化場分布結(jié)果 圖3顯示了錐形噴嘴在文丘里管中的空化過程中的壓力場變化，在t=0.0037s時開始出現(xiàn)空化現(xiàn)象，此時負壓區(qū)壓力為-9.78×104Pa，空化首先出現(xiàn)在噴嘴處，之后隨時間延長負壓區(qū)域逐漸擴大，在0.008s時達到穩(wěn)定，基本上遍布整個文丘里管。圖4是錐形噴嘴的空化云圖，同壓力云圖對照可知，空化場出現(xiàn)在流體的負壓區(qū)，最終基本遍布整個文丘里管。

圖3 錐形噴嘴的壓力云圖 圖4 錐形噴嘴的空化云圖

2)自激共振噴嘴的壓力場和空化場分布結(jié)果 圖5顯示了自激共振噴嘴在文丘里管中的空化發(fā)生過程，在t=0.00334s時達到空化條件，比錐形噴嘴空化開始早，但是最終達到穩(wěn)定空化時間比錐形噴嘴晚，負壓區(qū)域大。圖6顯示了自激共振噴嘴在文丘里管中的空化場分布，可以看出自激共振噴嘴空化現(xiàn)象出現(xiàn)時間早，且產(chǎn)生的空化場更大，比錐形噴嘴更有優(yōu)勢。

3)2種噴嘴湍流動能比較 由圖7和圖8可以清晰直觀的看到自激共振噴嘴的湍流動能最大值比錐形噴嘴強約55%，而且自激共振噴嘴在出口處湍流動能有一個增加過程，湍流動能在出口處有2處對稱增強區(qū)，因此也就解釋了為什么自激共振噴嘴的空化區(qū)域比較大。

圖5 自激共振噴嘴的壓力云圖 圖6 自激共振噴嘴的空化云圖

圖7 錐形噴嘴湍流動能曲線 圖8 自激共振噴嘴湍流動能曲線

4 結(jié) 論

1)自激共振噴嘴的空化現(xiàn)象出現(xiàn)時所需時間比錐形噴嘴少，達到穩(wěn)定的時間比錐形噴嘴長。

2)自激共振噴嘴的空化場范圍比錐形噴嘴空化場范圍大。

3)自激共振噴嘴的湍流動能比錐形噴嘴湍流動能強。

[1]沈壯志,柳楠.文丘里管反應(yīng)器空化泡的動力學(xué)特性[J].陜西師范大學(xué)學(xué)報(自然科學(xué)版),2011,40(1):23-28.

[2] 張曉冬,李志義,武君,等.水力空化對化學(xué)反應(yīng)的強化效應(yīng)[J].化工學(xué)報,2005,56(2):262-265.

[3] 關(guān)善超，沈壯志.水力空化發(fā)生器的結(jié)構(gòu)對空化場分布的影響[J].化學(xué)工程與裝備，2012(8)：52-55.

[4] 李根生,沈忠厚,周長山,等. 自振空化射流沖擊壓力脈動特性實驗研究[J]. 水動力學(xué)研究與進展(A輯),2003,18(5): 570-575.

[5] 王萍輝，朱單，馬飛.自激振動空化噴嘴的仿真研究[J].煤礦機械，2009,30(3):53-55.

[6] 劉厚林，劉東喜，王勇，等.3種空化模型在離心泵空化流計算中的應(yīng)用評價[J].農(nóng)業(yè)工程學(xué)報，2012,28(16)：54-59.

2012-10-24

國家自然科學(xué)基金項目(11174191)；陜西省自然科學(xué)基金項目(2011JM1002)；國家大學(xué)生創(chuàng)新訓(xùn)練計劃項目(CX12131)。

關(guān)善超(1989-)，男,現(xiàn)主要從事聲空化方面的學(xué)習(xí)工作。

O357.1

A

1673-1409(2013)01-0012-03

[編輯] 洪云飛