上旋式無(wú)窗散裂靶件自由界面數(shù)值研究

曾 月，顧漢洋，胡 晨

（上海交通大學(xué)核科學(xué)與工程學(xué)院，上海 200240）

上旋式無(wú)窗散裂靶件自由界面數(shù)值研究

曾 月，顧漢洋，胡 晨

（上海交通大學(xué)核科學(xué)與工程學(xué)院，上海 200240）

無(wú)窗散裂靶件是加速器驅(qū)動(dòng)的次臨界系統(tǒng)的重要組成部分。為了使靶件維持穩(wěn)定的自由界面并具有良好的熱移出功能，本文提出一種上旋式無(wú)窗靶件設(shè)計(jì)，以水為工質(zhì)，采用ANSYS Fluent軟件對(duì)該靶件設(shè)計(jì)進(jìn)行數(shù)值模擬。結(jié)果表明，該結(jié)構(gòu)散裂靶件能較好地維持穩(wěn)定自由界面，不會(huì)出現(xiàn)漩渦滯止區(qū)，并得到不同入口流速對(duì)自由界面形狀、流場(chǎng)特征及壓降的影響。

加速器驅(qū)動(dòng)的次臨界系統(tǒng)；無(wú)窗靶件；自由界面；VOF

加速器驅(qū)動(dòng)的次臨界系統(tǒng)（ADS）為國(guó)際公認(rèn)的從根本上處理嬗變長(zhǎng)壽命核廢物的新型核系統(tǒng)。散裂靶件是ADS的核心部件，該部件通過(guò)散裂反應(yīng)產(chǎn)生中子源，同時(shí)釋放大量熱量。為了在產(chǎn)生穩(wěn)定中子源的同時(shí)將散裂過(guò)程中產(chǎn)生的熱量有效帶出，ADS散裂靶件的發(fā)展經(jīng)歷了固態(tài)到液態(tài)、有窗到無(wú)窗的過(guò)程。無(wú)窗靶件無(wú)窗體材料壽命限制、散裂反應(yīng)熱沉積快速輸出等特點(diǎn)使之成為目前研究的熱點(diǎn)［1］。液態(tài)無(wú)窗靶件研究的重點(diǎn)包括以下兩點(diǎn)：一是如何在特定流動(dòng)邊界下維持穩(wěn)定光滑的自由界面，防止液滴飛濺污染質(zhì)子束入射管；二是如何盡可能減小與自由界面相耦合的非穩(wěn)態(tài)流場(chǎng)中滯止漩渦的尺度，增強(qiáng)流場(chǎng)的熱移出能力。國(guó)際上關(guān)于ADS的學(xué)術(shù)交流以及科技研究一直很活躍，美國(guó)、日本、俄羅斯和歐洲都有相應(yīng)的研究計(jì)劃，如ATW計(jì)劃、OMEGA計(jì)劃及SPIN計(jì)劃，其他的一些國(guó)家也參與到國(guó)際合作中，如韓國(guó)的KOMAC-HYPER、比利時(shí)的MYRRHA等［2］。多數(shù)無(wú)窗靶件結(jié)構(gòu)均與歐洲XT-ADS設(shè)計(jì)相似，上海交通大學(xué)也進(jìn)行了類似結(jié)構(gòu)的模擬與實(shí)驗(yàn)［3］。上海交通大學(xué)的實(shí)驗(yàn)和模擬結(jié)果顯示，錐形下落式無(wú)窗靶件的流場(chǎng)特點(diǎn)決定其存在自由界面穩(wěn)定性與減少滯止漩渦的尺度之間的矛盾［4］。本文提出一種上旋式結(jié)構(gòu)的無(wú)窗靶件，并采用ANSYS Fluent軟件對(duì)其自由界面進(jìn)行數(shù)值分析，驗(yàn)證其設(shè)計(jì)的可行性。

1 結(jié)構(gòu)設(shè)計(jì)

圖1為上旋式無(wú)窗散裂靶件結(jié)構(gòu)設(shè)計(jì)圖。上旋式無(wú)窗散裂靶件通道包括入口葉輪、圓柱直管段和喇叭口漸擴(kuò)段三部分。液態(tài)工質(zhì)從下部入口流入流道，經(jīng)由旋葉產(chǎn)生周向速度，從而旋轉(zhuǎn)向上，流過(guò)圓柱段至漸擴(kuò)段，從喇叭口流出，在喇叭段形成錐形自由液面。自由界面的位置可通過(guò)入口流速及下部葉輪葉片角度等來(lái)控制。

圖1 上旋式無(wú)窗散裂靶件結(jié)構(gòu)設(shè)計(jì)圖Fig.1 Design diagram of upward rotation windowless spallation target

2 計(jì)算模型

參照上旋式無(wú)窗靶件結(jié)構(gòu)設(shè)計(jì)圖建立數(shù)值模擬的幾何結(jié)構(gòu)模型及網(wǎng)格模型，如圖2所示。網(wǎng)格采用結(jié)構(gòu)化六面體網(wǎng)格。加密后網(wǎng)格單元總數(shù)為515萬(wàn)。數(shù)值模擬計(jì)算使用ANSYS Fluent軟件，采用k-ε模型對(duì)湍流特性進(jìn)行瞬態(tài)求解，自由界面重構(gòu)使用VOF方法，速度和壓力方程組的求解使用SIMPLE算法。邊界條件設(shè)置下部入口為速度入口，頂部出口為壓力出口。流體工質(zhì)為常溫、常壓下的水，相關(guān)資料［5］顯示，20℃的水與400℃的液態(tài)鉛鉍合金具有相似的水力學(xué)特性，故本文用水模擬液態(tài)鉛鉍合金進(jìn)行計(jì)算。計(jì)算主要研究入口速度和葉片傾斜角度對(duì)自由界面的影響，出口壓力設(shè)為大氣壓力。每個(gè)工況出口和入口的流量（min－mout）／min＜0.1%作為計(jì)算收斂標(biāo)準(zhǔn)。相關(guān)資料［6］顯示，液態(tài)無(wú)窗靶件的最大流速不能超過(guò)2m／s，所以最大入口流速為2m／s［6］。模擬計(jì)算中采用的葉片傾角為45°，入口流速分別為0.2、0.4、0.6、0.8、1.0、1.2、1.4、1.6、1.8、2.0m／s。

圖2 上旋式ADS液態(tài)無(wú)窗靶件計(jì)算網(wǎng)格模型Fig.2 Computing mesh model of upward rotation ADS liquid windowless target

3 計(jì)算結(jié)果分析

3.1 流動(dòng)結(jié)構(gòu)分析

圖3示出入口速度為1.2m／s時(shí)的流體旋流跡線，流體由下部入口進(jìn)入，經(jīng)過(guò)旋葉獲得切向速度旋轉(zhuǎn)向上，從周圍出口向下流出。為研究流體沿著高度方向軸向移動(dòng)的過(guò)程，在不同高度位置上，取6個(gè)橫截面，作出其速度云圖，如圖4所示。由圖4可看出，流體進(jìn)入流道后，經(jīng)過(guò)旋葉產(chǎn)生切向速度，形成旋轉(zhuǎn)上升流，在離心力的作用下，流體向管道周向移動(dòng)，在出口喇叭口區(qū)域形成自由界面。

圖3 流體旋流跡線Fig.3 Fluid flow trace

圖4 不同高度橫截面上速度云圖Fig.4 Velocity contours of cross section on different heights

3.2 自由液面形狀

計(jì)算葉輪傾角為45°時(shí)不同流速下的流動(dòng)情況，當(dāng)計(jì)算穩(wěn)定時(shí)，截取過(guò)中軸線的平面，可看到部分工況的自由液面分布如圖5所示。可見不同的流速下都能形成清晰的自由界面。圖6示出了不同流速下自由液面界面的變化曲線。計(jì)算結(jié)果顯示，不同入口流速下均能形成較為穩(wěn)定的自由液面，并隨入口流速的增大，流體流出的離心速度越大，導(dǎo)致自由液面中心更低，更接近圓錐形，在高速工況下計(jì)算結(jié)果表明自由液面會(huì)出現(xiàn)波動(dòng)?？偟膩?lái)講，由于離心力的作用，這種上旋自由界面較之下落自由界面更加穩(wěn)定，不存在漩渦滯止區(qū)，自由界面中心不會(huì)出現(xiàn)強(qiáng)烈的氣液攪混。

圖5 不同入口流速下的自由界面Fig.5 Free surfaces under different inlet velocities

圖6 入口流速與自由液面高度曲線Fig.6 Curve of inlet velocity and height of free surface

3.3 自由液面區(qū)流場(chǎng)特性

自由液面區(qū)流場(chǎng)特性是影響自由靶件沉積熱輸出的關(guān)鍵因素。故取不同入口流速下的自由液面區(qū)速度云圖，如圖7所示。

流體自下而上旋轉(zhuǎn)上升至自由界面，形成錐形面后從圓周邊緣旋轉(zhuǎn)溢出。由上述速度分布云圖可看出，無(wú)論入口速度較大或較小，自由液面下流體都保持一種穩(wěn)定狀態(tài)，不存在回流和漩渦滯止區(qū)。而且入口速度越小，自由界面下速度越小，流體的流動(dòng)越平穩(wěn)。

3.4 壓降與流速的關(guān)系

進(jìn)出口壓降與流速關(guān)系如圖8所示?？煽吹?，隨入口流速的增大，壓降大致呈線性變化，入口流速越大，進(jìn)出口壓降越大。但總的壓降不是很大，不需考慮壓力變化對(duì)材料及工質(zhì)造成的影響，且不需在高壓下工作，不存在工程實(shí)施困難。

圖7 不同入口流速下自由液面區(qū)速度云圖Fig.7 Velocity contours of free surface under different inlet velocities

圖8 進(jìn)出口壓降與進(jìn)口流速關(guān)系Fig.8 Inlet and outlet pressure drops under different inlet velocities

4 結(jié)論

本文提出了一種加速器驅(qū)動(dòng)的次臨界系統(tǒng)上旋式液態(tài)無(wú)窗散裂靶件并對(duì)其自由界面和流場(chǎng)進(jìn)行了數(shù)值模擬計(jì)算。通過(guò)對(duì)自由液面及自由界面下流場(chǎng)速度的對(duì)比分析，得出以下結(jié)論：

1）上旋式無(wú)窗散裂靶件幾乎不會(huì)出現(xiàn)自由界面的巨大波動(dòng)、撕裂及自由液面下的漩渦回旋等情況，解決了自由界面穩(wěn)定性與界面區(qū)出現(xiàn)流場(chǎng)渦流滯止之間的矛盾；

2）隨入口流速的變大，散裂靶件的自由液面的下凹越明顯；

3）入口流速越大，壓降越大，但總的壓降不是很大，不存在工程實(shí)施困難。

［1］ ROELOFS F，de JAGER B，CLASS A，et al.European research on HLM thermal-hydraulics for ADS applications［J］.Journal of Nuclear Materials，2008，376：401-404.

［2］ de BRUYN D.Accelerator driven systems for transmutation：Main design achievements of the XT-ADS and EFIT systems within the FP6IPEUROTANS integrated project［C］∥Proceeding of ICAPP10.CA：［s.n.］，2010：1 808-1 816.

［3］ CHENG X，PETTAN C，KNEBEL J U，et al.Experimental and numerical studies on thermalhydraulics of spallation targets［C］∥Emerging Nuclear and Transmutation Systems：Core Physics and Engineering Aspects.USA：［s.n.］，2000：308-325.

［4］ SU G Y，GU H Y，CHENG X.Experimental and numerical studies on free surface flow of windowless target［J］.Annals of Nuclear Energy，2012，43：142-149.

［5］ 蘇冠宇.ADS無(wú)窗散裂靶件模擬實(shí)驗(yàn)與數(shù)值研究［D］.上海：上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院，2011.

［6］ CLASS A G，ANGELI D，BATTA A.XT-ADS windowless spallation target thermohydraulic design ＆experimental setup［J］.Journal of Nuclear Materials，2011，415：378-384.

Numerical Investigation on Free Surface of Upward Rotation Windowless Spallation Target

ZENG Yue，GU Han-yang，HU Chen
（School of Nuclear Science and Engineering，Shanghai Jiao Tong University，Shanghai 200240，China）

The windowless spallation target is one of most important components of the accelerator driven sub-critical system.An upward rotation windowless spallation target which chooses water as the working medium was proposed.The target can keep the interface steady and move out the heat in time.The numerical simulation was proceeded by ANSYS Fluent.Through the study it is proved that a steady surface can be kept with this kind of target.There is no whirlpool exists and the velocity below interface as same as the pressure drop under various inlet velocities can be found out.

accelerator driven sub-critical system；windowless target；free surface；VOF

TL33

：A

：1000-6931（2015）01-0096-04

10.7538／yzk.2015.49.01.0096

2013-10-29；

2013-12-02

國(guó)家自然科學(xué)基金資助項(xiàng)目（91026020）

曾 月（1990—），女，四川成都人，碩士研究生，核能科學(xué)與工程專業(yè)