初旭陽，胡鵬華，王 攀，陳 剛

(核工業(yè)北京化工冶金研究院，北京 101149)

在地浸采鈾生產(chǎn)中，通常利用蒸發(fā)池來處理從井場多抽出的放射性廢液，以減少液態(tài)的、易擴散的放射性物質的體積；蒸發(fā)濃縮后的廢液，通常采用混凝土固化掩埋再處理。但蒸發(fā)池占地面積大，會對生態(tài)環(huán)境和地浸礦山擴大產(chǎn)能造成一定影響。

為了加快蒸發(fā)速率、減小蒸發(fā)池面積、提高蒸發(fā)效率，可采用多種不同的輔助蒸發(fā)方式，包括多級加熱[1]、浮動鋁片[2-3]、人工濕地[4]、機械噴霧[5-7]、輔助風速蒸發(fā)[8]等。

輔助風速蒸發(fā)技術利用泵輸送廢水，使掛簾形狀的纖維材料被廢水沾濕，依靠掛在支撐桿上的掛簾增加蒸發(fā)面積[9]，提高蒸發(fā)效率。除了蒸發(fā)面積影響蒸發(fā)效率外；蒸發(fā)面上的空氣流動也是一個主要因素，存在蒸發(fā)面上空氣流速越高，蒸發(fā)效率越高的正相關關系。內部掛簾折疊的密集程度會影響掛簾表面的空氣流速，因此掛簾如何布置是影響裝置整體蒸發(fā)效率的另一個因素。

為了使風能輔助蒸發(fā)裝置有足夠的表面空氣流速，提高裝置的蒸發(fā)效率，研究了簾間距、簾開角、簾懸掛方式、風向對于風能輔助蒸發(fā)裝置空氣流場的影響。

1 模型建立和邊界條件

1.1 模型建立與網(wǎng)格劃分

人力計算并不能很好地得到計算流體力學精確解[10-11]，通常依靠計算軟件來解決工程實際問題。ANSYS Fluent是國際通用的流體模擬軟件，在各類流動問題中的應用廣泛[12]，適用于此次的簾結構空氣流場模擬。

本次模擬使用ANSYS SpaceClaim DesignModeler進行三維模型的建立。在符合相似準則[13]的前提下，為了節(jié)省計算量，避免過于精細的結構帶來大量的網(wǎng)格數(shù)，對于W型懸掛方式，模型采用2個簾間隔來計算；對于Ⅰ型懸掛方式，模型采用3個簾間隔來計算。不同懸掛方式模型如圖1所示。簾厚度為10 mm，長和高均為2 m。設置長方體外流場，上下距離為500 mm，左右、前后距離為5 000 mm。

I型懸掛時，模擬簾間距10、30、50、100 mm的情況；W型懸掛時，模擬支撐管徑DN15、DN25、DN40(即對應簾間距23、32、48 mm)的空氣流場情況。簾開角模型是在W型懸掛、DN40支撐管徑模型的基礎上進行修改，僅模擬15°簾開角情況。

風向模擬角度為25°、45°、60°，風向模擬建模使用通風較充分的簾開角模型，如圖2所示。

圖2 風向模擬中的風向改變示意圖

由于各類懸掛方式有各種不規(guī)則的結構，不方便采用ICEM四面體網(wǎng)格進行手動多次劃分[14]，故上述模型統(tǒng)一采用自動六面體網(wǎng)格劃分，但在網(wǎng)格加密部分稍作區(qū)別。以單片懸掛模型的網(wǎng)格劃分為例(圖3)，所有網(wǎng)格劃分均保證簾材料部分的六面體網(wǎng)格質量Elements Quality>0.9，偏斜度Skewness<0.7。

1.2 邊界條件

考慮試驗場地的實地風速，設置進風口的邊界條件為速度入口，風速為9 m/s，風向垂直于進風口。在簾間距最小L=10 mm時(即雷諾數(shù)最小時)，雷諾數(shù)Re=ρvL/μ>4 000；其余簾間距下雷諾數(shù)更大，即入風流均設置為湍流。設置出口壓力(Gauge pressure)為0 Pa，進出口流體均為空氣，所有壁面均為無滑移固壁。在簾材料選材尚未確定的情況下，模擬時簾材料暫默認為鋁。計算采用RNGk-ε的數(shù)學模型。

圖3 單片懸掛模型網(wǎng)格劃分示意圖

2 簾間距與懸掛方式模擬研究

在簾開角為0°下，進行簾間距與懸掛方式模擬研究，模擬I型懸掛和W型懸掛方式下，簾間距對風能輔助蒸發(fā)裝置空氣流場的影響。對比I型懸掛和W型懸掛的空氣流場，分析懸掛方式對空氣流場的影響。

2.1 I型懸掛的簾間距模擬

簾間距對于空氣流場的影響主要反映為簾間風速的變化。對I型懸掛方式建模計算，取ZX方向(與簾材料垂直方向)的截面觀察流場情況。100 mm簾間距下，在ZX方向Y=1 m截面處的速度云圖如圖4所示。為了研究具體簾間距對截面上的速度影響，對不同簾間距的ZX方向橫截面圖進行了放大處理，如圖5所示。10 mm簾間距ZY方向橫截面速度放大云圖如圖6所示。

圖4 ZX方向橫截面速度總覽

(a)100 mm;(b)50 m;(c)30 mm;(d)10 mm;圖5 不同簾間距的ZX方向橫截面速度放大云圖

圖6 10 mm簾間距ZY方向橫截面速度放大云圖

由圖5可見，除10 mm簾間距外，不同簾間距的空氣流場均呈近壁面的速度分層現(xiàn)象，出現(xiàn)了湍流邊界層。這和Ludwig Prandtl提出的流動邊界層、附面層一致[15]，同時簾中心都有比較穩(wěn)定的空氣流速。但從圖5(d)和圖6可以看到，在10 mm的簾間距情況下，簾與簾中心存在大量流速為0的區(qū)域(深藍色)，此時裝置整體蒸發(fā)效率受流速過低影響。

為研究不同簾間距對表面風速的影響，取這4種情況下的簾間隙中心流速進行對比，結果見表1。

表1 I型懸掛方式不同簾間距下的簾間隙中心流速

由表1可知，對于I型懸掛方式，除簾間距10 mm以外，簾間隙中心流速總體隨簾間距的增大而提高；但流速上升有限，簾間距為100 mm和30 mm時，簾間隙中心流速在數(shù)值上只相差0.8 m/s。簾間距過小(10 mm)，將會導致流體不能充分流過簾的間隙。

2.2 W型懸掛的簾間距模擬

對W型懸掛方式建模計算，同樣取ZX橫截面，DN40管徑模擬橫截面如圖7所示。可以看出，與I型模擬時類似，也存在邊界層現(xiàn)象；在省略管徑的間隙中心部位(掛簾彎曲處)的風速比正中心的間隙中心部位的風速稍高。

圖7 DN40管徑模擬橫截面速度云圖

取不同簾間距(管徑)下Y軸中心處的簾間隙中心流速進行對比，結果見表2?？梢钥闯?，簾間距對于簾間隙中心流速的影響并不明顯。

表2 W型懸掛方式不同簾間距下的簾間隙中心流速

2.3 小結

結合I型和W型的模擬結果，簾間距大于23 mm的情況下，簾間距對于表面流速的影響不大；但在簾間距小于23 mm的情況下(本研究僅模擬了10 mm的情況)，簾間距將可能對表面流速造成較大影響。

對比2種懸掛方式，發(fā)現(xiàn)W型懸掛支撐管附近會產(chǎn)生一個風速稍高的區(qū)域，可能會導致裝置有些許受力不均；但考慮W型懸掛方式比I型懸掛方式的操作要簡易很多，這一影響可忽略。

綜上考慮，裝置宜采用W型懸掛方式，并保證簾間距大于23 mm(使用管徑大于DN15的支撐管)。

3 簾開角模型模擬研究

簾開角是指簾與鉛垂面的夾角，數(shù)值上為兩簾之間夾角的1/2。在簾開角15°時，Y軸中心處ZX截面和X軸中心處ZY截面的速度云圖如圖8所示。

圖8 簾開角15°時橫截面風速云圖

從圖8(a)可看出，簾開角為15°時，在兩簾之間存在主流區(qū)，風速為9 m/s左右；靠近簾的部分，邊界層十分明顯，具有充分的空氣流場。相比DN40的中心流速，在有簾開角的情況下，風速僅提升了0.3 m/s左右。簾開角的設置會使掛簾的表面風速有微弱提升；但在實際裝置設計中，設置簾開角有助于設置落水孔，對裝置布液有較大益處。

4 風向模型模擬研究

選取之前效果最佳的DN40(48 mm)簾間距、15°簾開角作為模型，研究風向與簾材料夾角為25°、45°、60°時空氣流場的變化。這3種風向情況下，Y軸中心ZX截面的空氣流速如圖9所示，25°風向時的流線圖如圖10所示。3種不同風向下的簾間隙中心流速和繞流最大風速見表3。

表3 風向變化對空氣流場變化數(shù)據(jù)

從試驗結果看出，風向25°時，裝置附近未產(chǎn)生太多渦流，對裝置幾乎沒有影響(圖9)；但在風速較高的地方(圖10紅色區(qū)域)產(chǎn)生了繞流，可能會導致裝置不穩(wěn)定。在45°和60°風向情況下，簾間隙風速將會受到很大影響，且風向角度越大風速越小；而且45°、60°風向的繞流速度過大(表3)，將可能導致裝置傾翻。

(a)25°;(b)45°;(c)60°。圖9 不同風向下的ZX截面空氣流速云圖

圖10 風向25°時的空氣流線圖

5 模擬結果驗證試驗

利用圖11的實驗室裝置進行試驗，因實驗室風扇最高風速約為4.3 m/s，無法達到現(xiàn)場9 m/s的風速；因此改變模擬進口風速為4.30 m/s來驗證模擬結果的準確性。經(jīng)計算，風速為4.30 m/s時，風流依然為湍流。

利用風速儀測量進風口和出風口15 cm位置的風速，作為掛簾間隙空氣流速進行比對分析。進口和出口分別選擇3個高度測量，試驗結果見表4?？梢钥闯觯诓煌c位，風速降低的趨勢和幅度具有一致性。與模擬結果對比，模擬的進、出口平均風速分別為4.30 m/s和3.68 m/s，與實驗組的4.20 m/s和3.38 m/s誤差不大。由此可見，使用數(shù)值模擬方法研究風能輔助蒸發(fā)裝置布簾結構對空氣流場的影響情況是可行的，通過模擬得出的有關簾間距、簾開角、懸掛方式、風向對空氣流場影響的結論可信。

圖11 風能輔助蒸發(fā)實驗裝置

表4 風能輔助蒸發(fā)實測結果與數(shù)值模擬結果對比

6 結論

在研究條件下，風向角度是影響風能輔助蒸發(fā)裝置蒸發(fā)效率和穩(wěn)定性的關鍵因素，風向與簾面平行時蒸發(fā)效率最佳。由于邊界層效應，當簾間距過小時，簾間距對蒸發(fā)效率有影響；當簾間距大于23 mm時，簾間距對蒸發(fā)效率的影響不顯著。簾開角對裝置蒸發(fā)效率影響不大，但有助于裝置布液。懸掛方式對空氣流場無太大影響，但考慮到施工簡易程度，推薦選用W型懸掛方式。

建議在風能輔助蒸發(fā)裝置設計中，簾間距大于23 mm(DN15)，采用W型懸掛方式，合理設置簾開角；并設計裝置能夠根據(jù)氣象參數(shù)變化自動旋轉，使風向平行吹入簾與簾的間隙。