連續(xù)液相碎化過程數(shù)值模擬研究

＊蘭治科 李勇 蘇光輝 昝元鋒

（1.中國核動力研究設(shè)計院 四川 610041 2.西安交通大學(xué)核科學(xué)與技術(shù)學(xué)院 陜西 710049）

液體霧化過程包含了能量和質(zhì)量轉(zhuǎn)化、連續(xù)液體顆?；任锢磉^程，在現(xiàn)代社會被廣泛應(yīng)用各個領(lǐng)域，如農(nóng)業(yè)噴霧、噴涂、發(fā)動機(jī)、航空航天、核能等。根據(jù)用途不同，對霧化性能的要求也不盡相同。因此獲得特定噴霧頭的相關(guān)性能參數(shù)顯得尤為重要，特別是霧化后的液滴粒徑、流量密度分布和霧化錐角。

由于噴霧機(jī)理的復(fù)雜性，到目前為止還沒有能夠完整、準(zhǔn)確描述噴霧頭性能的好方法。目前的研究方法主要有三種，分別為實驗研究、理論計算和計算流體力學(xué)（CFD）模擬。其中理論計算主要是建立在半經(jīng)驗公式的基礎(chǔ)上，而關(guān)于實驗研究，高昂的實驗成本以及測量精度，仍是主要的制約因素。特別是應(yīng)用于某些特殊領(lǐng)域的大型噴霧頭，如穩(wěn)壓器、安全殼中的噴霧設(shè)備，其共同特點是：噴霧頭尺寸和噴霧流量大，液滴粒徑分布范圍較廣，普通粒徑測量儀器不適用于該類型噴霧頭測量。此外，對于噴霧頭內(nèi)部流場以及出口處液膜破碎過程等微觀現(xiàn)象研究，目前尚未找到適合的實驗方法。

隨著計算機(jī)性能提高，CFD被廣泛應(yīng)用于霧化性能各項參數(shù)研究。采用數(shù)值模擬方法研究霧化現(xiàn)象有助于更加深入認(rèn)識霧化機(jī)理和指導(dǎo)相關(guān)實驗，同時可以彌補(bǔ)噴霧實驗不足，提供噴霧實驗無法得到的流場分析結(jié)果。

然而，對于液體霧化數(shù)值模擬，CFD技術(shù)仍然面臨著較大挑戰(zhàn)。主要是因為霧化過程不僅牽涉到從單相液體到氣液共存兩相流動的過渡，同時，對于液相來說，還存在從連續(xù)相的液體到離散相液滴的轉(zhuǎn)變。目前幾乎所有的霧化場模擬均使用Eulerain-Lagrangian兩相流模型[1]。這種模型通常假定噴口附近的液滴粒徑和速度分布作為已知的邊界條件已經(jīng)被確定，而獲得這些邊界條件的方法又往往依賴實驗或者經(jīng)驗。因此，從嚴(yán)格意義上說，這種數(shù)值方法并不能實現(xiàn)噴霧過程的真正模擬。而Eulerain-Eulerain方法，相對于網(wǎng)格尺寸而言，僅局限于對兩種連續(xù)相共存的流場進(jìn)行計算，并不能獲得完全離散化的液滴[2]。

基于以上原因，本文根據(jù)最大不穩(wěn)定波增長率液膜破碎理論，在連續(xù)液相出口位置構(gòu)建初級破碎模型-PBM，與噴霧頭內(nèi)部流場和霧化場的數(shù)值模擬進(jìn)行耦合，提出了一種全流場的數(shù)值模擬方法，實現(xiàn)從噴霧頭入口到霧化場的全過程數(shù)值模擬。最后，通過已有的實驗數(shù)據(jù)對以上方法進(jìn)行了驗證，計算結(jié)果符合較好，證明了全流場數(shù)值模擬方法的有效性。

1.全流場模型構(gòu)建

目前，在已有的兩相流模擬中，Eulerain-Eulerain法將兩相視為能夠互相穿插，但不相溶的連續(xù)介質(zhì)。該方法使用空間和時間上為連續(xù)函數(shù)的相體積分?jǐn)?shù)來模擬兩相在空間的分布。從噴霧頭入口全液相流動到出口連續(xù)液膜與氣相并存，符合Eulerain-Eulerain方法兩相模擬的要求。

Eulerain-Lagrangian方法用于噴霧場模擬，將噴霧空間的氣體介質(zhì)作為連續(xù)介質(zhì)，為其建立Eulerain坐標(biāo)系；而將液滴顆粒當(dāng)作離散介質(zhì)，為其建立Lagrangian坐標(biāo)系。該方法假定液相以離散液滴的形式進(jìn)入連續(xù)氣相中，通過質(zhì)量交換、能量交換和動量交換進(jìn)行兩個相間的耦合。因此Eulerain-Lagrangian方法僅用于噴霧頭出口連續(xù)液膜破碎后的兩相流場模擬。

本文使用Eulerain-Eulerain法中VOF模型捕捉兩相之間界面，模擬從單相液體進(jìn)入噴霧頭，直到出口位置形成液膜的內(nèi)部流場；使用Eulerain-Lagrangian方法中DPM模型模擬離散液滴在氣相中的碰撞、破碎以及再分布過程。本文采用最大表面擾動波增長率破碎理論，構(gòu)建了適合于大型噴霧頭的初級破碎模型（PBM）。通過在FLUENT求解器上使用UDF，編譯相關(guān)程序與以上兩個模型耦合，形成VOF+PBM+DPM的完整流場的數(shù)值模擬方法，模型結(jié)構(gòu)見圖1所示。

圖1 連續(xù)液相碎化過程數(shù)值模擬模型結(jié)構(gòu)圖

2.全流場數(shù)值模擬過程及結(jié)果分析

(1)噴霧頭內(nèi)部流場數(shù)值模擬

VOF方法是在1981年由Hirt和Nichols[3]首先提出，它使用一個標(biāo)量場函數(shù)γ體積份數(shù)來表征第二相流體在計算網(wǎng)格中的分布，其大小為第二相流體在網(wǎng)格中所占體積百分比。因此，VOF又被稱為流體體積法。在本文霧化模擬中，環(huán)境介質(zhì)空氣作為第一相流體，霧化介質(zhì)水作為第二相流體，所以體積分?jǐn)?shù)γ代表網(wǎng)格控制體中霧化介質(zhì)水的體積份額。

根據(jù)體積分?jǐn)?shù)的標(biāo)量場γ，全計算域流體物性密度ρ以及粘性μ的表達(dá)式如下：

其中ρ1、ρg分別為液相和氣相所對應(yīng)的密度；μ1、μg分別為液相和氣相所對應(yīng)的粘性應(yīng)力。因此，在液相、氣相以及兩相混合區(qū)域均使用式（2）中流體物性，可以使在流體力學(xué)相關(guān)控制方程中按單流體模型僅使用一套方程組。

體積分?jǐn)?shù)函數(shù)γ滿足如下輸運方程：

噴霧頭內(nèi)部結(jié)構(gòu)和流道對霧化特性起決定作用，因此在計算過程中需要對網(wǎng)格質(zhì)量進(jìn)行獨立性驗證。本文定義了4種不同的網(wǎng)格方案，對網(wǎng)格逐漸加密對流量特性及霧化錐角的影響進(jìn)行了比較。當(dāng)網(wǎng)格精度進(jìn)一步提高，以上特性參數(shù)不再明顯變化時即認(rèn)為網(wǎng)格質(zhì)量滿足獨立性條件。

本文為螺旋槽式旋流噴頭提供連續(xù)液相入口，在大多數(shù)情況下均表現(xiàn)為強(qiáng)烈的旋流特征，為了避免標(biāo)準(zhǔn)κ-ε模型的缺陷，湍流模型考慮κ-ε模型的修正方案：RNGκ-ε模型。FLUENT求解器三維旋流計算中通過旋流因子as對RNG κ-ε模型中的μt進(jìn)行修正。計算前對as=0.070，0.075和0.080三種不同旋流因子時流動特性計算值與實驗值進(jìn)行了比較，認(rèn)為as=0.080作為湍流模型旋流因子比較適合該類型噴霧頭的旋流強(qiáng)度。

連續(xù)液相破碎數(shù)值模擬，使用VOF模型從噴霧頭入口為時間和迭代開始起始點。圖2為噴霧頭內(nèi)部流場瞬態(tài)圖。噴霧頭在流動擾動過程中，在噴霧頭混合腔接近出口位置出現(xiàn)較為明顯的旋流、直流分層流動。其中直流部分位于正中心豎直向下，而旋流則沿四周壁面繞流。而在旋流作用下，混合腔內(nèi)部易形成負(fù)壓區(qū)，容易發(fā)生大氣壓力下環(huán)境介質(zhì)中空氣倒流現(xiàn)象。而直流與旋流交界處則成為相對較為“脆弱”的位置。

圖2 連續(xù)液相碎化為液帶過程圖

(2)噴霧頭出口初級破碎模擬

在噴霧頭出口建立液膜射流的幾何模型，假定液膜以錐形射進(jìn)入環(huán)境空氣。液膜與空氣相互作用，在液膜表面上形成擾動波。

通常認(rèn)為，擾動波滿足如下關(guān)系式：

假設(shè)擾動波是導(dǎo)致液膜破碎[4]的主要原因，當(dāng)擾動波增長率達(dá)到最大值時，液膜發(fā)生破裂，形成初始液滴。通過色散方程對式（4）進(jìn)行求解，最終得到最大增長率對應(yīng)下液膜破碎破碎模型PBM。

使用UDF中的DEFINE_EXECUTE_AT_END宏，在瞬態(tài)計算的每個時間步長統(tǒng)計噴霧頭孔口所在面每個網(wǎng)格單元的流動信息。通過采集兩相VOF模型在噴霧出口的計算結(jié)果作為邊界條件輸入霧化初級破碎模型。對旋流噴霧頭，假定孔口液膜以環(huán)狀形式存在，通過液相截面積相等的方法獲得PBM初始液膜厚度。

通過PBM模型獲得破碎距離后，在FLUENT求解器中通過UDF功能統(tǒng)計某一時間段內(nèi)在該距離橫截面上瞬態(tài)流量分布，獲得該截面上霧化場流量密度分布。UDF由編輯宏DEFINE_EXECUTE_AT_END程序以及定義F_UDMI存儲瞬態(tài)數(shù)據(jù)完成，在橫截面每個face類型的網(wǎng)格單元中累積該時間步下經(jīng)過該網(wǎng)格的質(zhì)量流量，最終得到某一時間段內(nèi)截面上的流量分布如圖3所示。

圖3 液膜破碎截面流量密度分布

(3)噴霧場數(shù)值模擬

在DPM模型中，通過積分拉式坐標(biāo)下的顆粒作用力微分方程來求解離散相顆粒的軌道。顆粒運動方程：

其中u、up分別為流體和顆粒的速度；ρ、ρp分別為流體和顆粒的密度；FD(u-up)為顆粒的單位質(zhì)量曳力。此外，液滴碰撞模型基于O'Rourke理論；液滴破碎模型適用TAB模型。

為了降低計算成本，本文構(gòu)建2-D幾何模型對噴霧場豎直截面內(nèi)離散液滴流動特性進(jìn)行數(shù)值模擬。該幾何模型在原有VOF模型的正下方，選取過噴霧頭軸心的豎直截面作為計算區(qū)域，其大小為4000mm×2200mm。

通過DPM模型可得到噴霧場溫度、速度等各項參數(shù)分布。圖4為霧化場液滴粒徑分布計算結(jié)果。從圖中可以看出，噴霧場中較大粒徑液滴主要位于霧化場外側(cè)。原因在于噴霧頭霧化過程以旋流特征為主，即大部分霧化液體以旋流方式沿霧化錐角外側(cè)方向進(jìn)入噴霧場。其次，霧化場粒徑分布圖中左半側(cè)液滴粒徑均值明顯大于右半側(cè)，且顆粒分布較密。主要是因為在離散相邊界條件上，噴霧頭霧化液體并非對稱分布，大部分液體分布在邊界條件的左側(cè)區(qū)域。

圖4 噴霧場液滴粒徑分布

3.全流場數(shù)值模擬與實驗結(jié)果對比

針對本文研究對象，在數(shù)值模擬工況下開展了相應(yīng)實驗研究，具體實驗研究內(nèi)容和方法詳見文獻(xiàn)[5]，此處不再贅述。

圖5給出了通過全流場數(shù)值模擬方法計算得到的噴霧場液滴粒徑概率密度與實驗數(shù)據(jù)的比較。圖中橫坐標(biāo)為液滴粒徑尺寸，曲線為計算得到粒徑分布，離散點為實驗測試得到的實際粒徑概率分布，兩者分布趨勢符合較好，表明了全流場數(shù)值模擬的有效性。

圖5 液滴粒徑概率密度計算值與實驗數(shù)據(jù)比較

4.結(jié)論

本文針對連續(xù)液相，提出一種全流場液相破碎數(shù)值模擬方法。通過構(gòu)建初級破碎模型PBM，與用于內(nèi)部流場計算的VOF模型和霧化場計算的DPM模型進(jìn)行耦合，實現(xiàn)了從連續(xù)液相到霧化顆粒的噴霧全過程數(shù)值模擬。最后與實驗結(jié)果對比，證明了VOF+PBM+DPM全流場模型的有效性。解決了目前霧化數(shù)值模擬中，Eulerain-Eulerain和Eulerain-Lagrangian方法僅能模擬局部過程的問題，為連續(xù)液相破碎及霧化場研究提供了一套完善的數(shù)值計算方法。