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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

圖4 噴霧場(chǎng)液滴粒徑分布

3.全流場(chǎng)數(shù)值模擬與實(shí)驗(yàn)結(jié)果對(duì)比

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

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

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

4.結(jié)論

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