陳云龍，黃群英，朱志強(qiáng)，高 勝，吳慶生，F(xiàn)DS團(tuán)隊

（1．中國科學(xué)院核能安全技術(shù)研究所，安徽 合肥230031；

2．中國科學(xué)技術(shù)大學(xué)核科學(xué)與技術(shù)學(xué)院，安徽 合肥230027）

聚變能是目前認(rèn)識到的可最終解決能源危機(jī)和環(huán)境問題的重要途徑之一［1－4］。聚變堆液態(tài)金屬鋰鉛包層因其固有的特性和優(yōu)勢，被國際上普遍認(rèn)為是未來聚變電站最具發(fā)展前景的包層設(shè)計方案之一［5－9］。目前，國際熱核聚變實驗堆（ITER）成員方中，歐盟、美國和中國等均將液態(tài)金屬鋰鉛包層作為主要研究方案而重點(diǎn)發(fā)展［10－12］。

聚變堆中，包層所處的環(huán)境十分苛刻，液態(tài)金屬鋰鉛在包層和輔助系統(tǒng)中的流動將帶來極具挑戰(zhàn)的科學(xué)和技術(shù)問題，如：鋰鉛純化技術(shù)、傳熱和磁流體動力學(xué)（MHD）效應(yīng)等，對這些關(guān)鍵技術(shù)問題進(jìn)行研究，為實現(xiàn)聚變堆商業(yè)化具有重要意義。

液態(tài)鋰鉛純化技術(shù)是聚變堆鋰鉛包層關(guān)鍵技術(shù)之一，主要用以純化液態(tài)金屬，以避免雜質(zhì)長期在管道中堆積堵塞回路，影響回路中鋰鉛成分，并提高設(shè)備的使用壽命。冷阱作為國際上液態(tài)金屬回路的主要純化裝置之一，其根據(jù)液態(tài)金屬中雜質(zhì)溶解度隨溫度變化的原理，使雜質(zhì)結(jié)晶成核、沉淀析出進(jìn)而分離［13］。

FDS團(tuán)隊具有中國自主知識產(chǎn)權(quán)的DRAGON-Ⅳ鋰鉛實驗回路是國際上先進(jìn)的多功能強(qiáng)迫對流鋰鉛實驗回路之一（圖1）［14－16］，它包含多個子實驗系統(tǒng)如純化系統(tǒng)、TBM實驗段與高溫試驗段等，能開展多種實驗研究。本文擬以DRAGON-Ⅳ鋰鉛實驗回路冷阱為研究對象，利用數(shù)值模擬手段，討論分析單層過濾芯結(jié)構(gòu)參數(shù)（包括過濾芯直徑和過濾芯濾孔）及液態(tài)金屬流動速度對捕集效率的影響，為冷阱過濾芯的選取提供理論依據(jù)與參考。

圖1 DRAGON-Ⅵ強(qiáng)迫對流鋰鉛實驗回路Fig．1 DRAGON-Ⅵforced convection lithium lead test loop

1 冷阱過濾芯描述

DRAGON-Ⅳ冷阱選取SS316L作為結(jié)構(gòu)材料，其中過濾芯選取耐高溫鋰鉛腐蝕的鉬（Mo）絲。其中過濾芯的幾何參數(shù)是影響冷阱的捕集效率的重要因素，應(yīng)適當(dāng)選取過濾芯絲網(wǎng)，以免導(dǎo)致冷阱堵塞或雜質(zhì)捕集效率過低。

由于目前液態(tài)鋰鉛冷阱參數(shù)較少，因此設(shè)計參數(shù)主要參考技術(shù)相對成熟的鈉回路，鈉冷阱過濾芯濾孔孔徑一般在100μm以下［17］。參考工業(yè)絲網(wǎng)編織標(biāo)準(zhǔn)（如孔徑為100μm時，絲徑為65μm；孔徑為70μm 時，絲徑為45μm），為了分析過濾網(wǎng)單一幾何參數(shù)對捕集效率的影響，另外再分析計算孔徑為100μm，絲徑為45μm的絲網(wǎng)的捕集效率。

2 數(shù)值分析計算

2．1 液態(tài)鋰鉛及雜質(zhì)物性參數(shù)

由于鋰鉛熔點(diǎn)為235℃，且等溫冷阱［18］過濾區(qū)通過均勻布置過濾芯，可使雜質(zhì)均勻析出，提高過濾芯使用壽命。故本文假設(shè)冷阱過濾區(qū)溫度恒為250℃，鋰鉛主要物性參數(shù)隨溫度的變化可表述為式［19］（1）～式（3）。

密度，kg／m3：

其中R＝8．31J／（mol·K），T 為溫度，單位為K。

回路主體結(jié)構(gòu)材料與高溫液態(tài)鋰鉛長期接觸腐蝕導(dǎo)致雜質(zhì)的產(chǎn)生，一般情況下鋰鉛冷阱中雜質(zhì)顆粒尺寸主要集中在2～10μm［20］。

2．2 計算模型

2．2．1 數(shù)學(xué)模型

由于液態(tài)鋰鉛在冷阱中降溫后有固態(tài)雜質(zhì)析出，冷阱內(nèi)部流體運(yùn)動為兩相流?？紤]回路主體結(jié)構(gòu)在流動鋰鉛中的年腐蝕量為幾十微米［21］，故假設(shè)液態(tài)鋰鉛析出雜質(zhì)的體積含量小于10%，冷阱內(nèi)兩相流動為稀疏懸浮流，符合拉格朗日離散相模型（DPM模型）計算條件。

離散相模型中顆粒的軌跡通過積分拉氏坐標(biāo)系下的顆粒作用力微分方程來求解。顆粒在流場中的受力比較復(fù)雜，主要有粘性力、重力、浮力、曳力等。顆粒的作用力平衡方程在笛卡爾坐標(biāo)系下的形式（X方向）為［22］

式中：u是流體相速度；uP是顆粒速度；μ為流體動力黏度；ρ為流體密度；ρp為顆粒密度；dp為顆粒直徑；Re為相對雷諾數(shù)（顆粒雷諾數(shù)），其定義為式中：對于一定的球形顆粒，在一定的雷諾數(shù)范圍內(nèi)，a1、a2、a3為常數(shù)。

2．2．2 邊界條件

由于液態(tài)鋰鉛中雜質(zhì)降溫結(jié)晶析出需要時間，故冷阱中鋰鉛流速不能過大，DRAGON-Ⅳ回路冷阱中鋰鉛流速為0．01 m／s左右，本文考慮流速對雜質(zhì)捕集效率的影響，設(shè)鋰鉛流速為0．005 m／s、0．01 m／s、0．015 m／s、0．02 m／s；模型中垂直于過濾芯的四壁面設(shè)置為對稱邊界條件。

離散相邊界條件設(shè)置如下：離散相速度大小設(shè)為與鋰鉛流動速率相同，方向垂直于入口截面；顆粒沿噴射面均勻分布；離散相濃度：根據(jù)結(jié)構(gòu)材料在流動鋰鉛中年腐蝕速率估算，鋰鉛中離散相雜質(zhì)濃度約為1 000 ppm左右；離散相溫度與液相相同；固壁邊界：網(wǎng)絲壁面設(shè)置完全捕集。

3 雜質(zhì)濃度場分布

利用計算流體力學(xué)軟件Fluent計算在不同流速下不同顆粒尺寸雜質(zhì)（5μm、7μm、10μm）在上面三種規(guī)格的過濾網(wǎng)周圍的濃度場分布。

考慮到鋰鉛流動雷諾數(shù)較小（Re為0．84～6．8），采用層流計算模型。求解過程先采用一階迎風(fēng)格式，使計算快速收斂，再轉(zhuǎn)入二階迎風(fēng)格式計算，提高計算精度，得出連續(xù)相流場，再加入離散相顆粒，考慮離散相與連續(xù)相的相間耦合，交替計算連續(xù)相和離散相直到兩相的計算結(jié)果到達(dá)收斂標(biāo)準(zhǔn)。

以過濾芯直徑65μm，孔徑大小100μm的過濾網(wǎng)為例，流速為0．015 m／s時不同顆粒直徑（5μm／7μm／10μm）下的雜質(zhì)濃度場分布如圖2所示。

圖2 三種尺寸雜質(zhì)的濃度場分布（a）5μm，（b）7μm，（c）10μmFig．2 Three sizes of impurity concentration distribution（a）5μm，（b）7μm，（c）10μm

從濃度場分布圖可以看出：雜質(zhì)顆粒濃度相對較大的區(qū)域基本在過濾芯表面，而空間濃度為較小值，說明較大部分雜質(zhì)顆粒被過濾芯捕集。

4 捕集效率計算結(jié)果與分析

根據(jù)計算出的濃度場分布，可計算過濾芯對雜質(zhì)的捕集效率，考慮到計算時過濾芯模型四壁為對稱邊界條件，絲網(wǎng)的捕集效率ηST定義為

式中：C入口為圖2中入口雜質(zhì)濃度，C出口為出口雜質(zhì)濃度。

圖3～圖5為三種規(guī)格過濾芯在鋰鉛流速為0．005m／s、0．01m／s、0．015m／s、0．02m／s時，對不同顆粒尺寸雜質(zhì)的捕集效率。

圖3 絲網(wǎng)捕集效率與流速之間的關(guān)系（過濾芯直徑65μm、網(wǎng)孔孔徑100μm）Fig．3 Relations of filter trapping efficiency and flow velocity（Diameter 65μm，pore size 100μm filter）

圖4 絲網(wǎng)捕集效率與流速之間的關(guān)系（過濾芯直徑45μm、網(wǎng)孔孔徑70μm）Fig．4 Relations of filter trapping efficiency and flow velocity（Diameter 45μm，pore size 70μm filter）

圖5 絲網(wǎng)捕集效率與流速之間的關(guān)系（過濾芯直徑45μm、網(wǎng)孔孔徑100μm）Fig．5 Relations of filter trapping efficiency and flow velocity（Diameter 45μm，pore size 100μm filter）

4．1 流速對捕集效率的影響

從圖3～圖5可知過濾芯的捕集效率隨鋰鉛流速的增加而增大，且隨著鋰鉛流速增加到一定值，捕集效率趨于平緩，計算結(jié)果與過濾絲網(wǎng)對氣相中液滴捕集效率隨流速變化規(guī)律類似［23］。這可能是由于流速增大，單位時間經(jīng)過過濾芯的雜質(zhì)數(shù)目變多，導(dǎo)致過濾芯的捕集效率增大。另外鋰鉛流速過大，不利于雜質(zhì)結(jié)晶析出，因此實際設(shè)計冷阱時，需要考慮雜質(zhì)捕集效率與雜質(zhì)結(jié)晶析出速度的關(guān)系，以獲得鋰鉛最佳流速。

4．2過濾芯幾何參數(shù)對捕集效率的影響

比較圖3和圖5可知在過濾芯孔徑相同時，過濾芯直徑大，捕集效率較高，這是因為此時過濾芯的孔隙率小，過濾芯與雜質(zhì)的相對接觸面積大。同時隨著流速增加，兩種過濾芯捕集效率差距減少。

比較圖4和圖5可知：過濾芯直徑相同時，過濾芯孔徑較小，捕集效率較高，這也是因為過濾芯孔隙率小，過濾芯與雜質(zhì)的相對接觸面積大。

綜上所述，低鋰鉛流速時，增大流體速度和降低過濾芯孔隙率可以提高過濾芯對雜質(zhì)捕集效率，且隨著流速增大到一定數(shù)值，過濾芯對雜質(zhì)捕集效率趨于平緩。

5 結(jié)論

本文用數(shù)值模擬方法對冷阱過濾芯周圍兩相流進(jìn)行三維數(shù)值模擬，模型建立及邊界條件設(shè)置最大限度接近真實工況，計算表明：單層過濾絲網(wǎng)對鋰鉛中雜質(zhì)的捕集效率隨流體速度的增大而增大，隨過濾芯孔隙率的減小而增大，數(shù)值模擬結(jié)果可為冷阱過濾芯的優(yōu)化設(shè)計及液態(tài)金屬流速的確定提供理論依據(jù)與參考。

本文計算結(jié)果以鋰鉛溫度為250℃作前提，未來可繼續(xù)開展其他溫度下鋰鉛中雜質(zhì)捕集情況的數(shù)值計算和實驗研究，以進(jìn)一步對鋰鉛冷阱進(jìn)行設(shè)計優(yōu)化。

致謝

特別感謝FDS團(tuán)隊老師對本文工作予以的幫助。

［1］ Yican Wu．Conceptual design activities of FDS series fusion power planet in China［J］．Fusion Eng Des，2006，81（23-24）：2713-2718．

［2］ Alejaldre C，De Marco F，F(xiàn)inzi U，et al．Status report on fusion research［J］．Nucl Fusion，2005，45（10A）：1-28．

［3］ Yican Wu，F(xiàn)DS Team．Conceptual Design of the China Fusion Power plants in China［J］．Fusion Eng Des，2008，83：1683-1689．

［4］ Ilhli T，Basu T，Giancarli，et al．Review of blanket designs for advanced fusion reactors［J］．Fusion Eng Des，2008，83（7-9）：912-919．

［5］ Yican Wu．Design Status and development strategy of China liquid lithium-lead blankets and related material technology［J］．J Nucl Mater，2007，367：1410-1415．

［6］ Salavy J F，LV Boccaccini，R L?sser，et al．Overview of the last progresses for the European Test Blanket Modules projects［J］．Fusion Engineering and Design，2007，82：2105-2122．

［7］ Wu Y C，Qian J P，Yu J N．The Fusion-Driven Hybrid System and Its Material Selection［J］．Journal of Nuclear Materials，2002，307-311：1629-1636．

［8］ Ying A，Abdou M，Wong C，et al．An overview of US ITER test blanket module program ［J］． Fusion Engineering and Design，2006，81：433-411．

［9］ Yican Wu， FDS Team． Fusion-Based Hydrogen Production Reactor and Its Material Selection［J］．Journal of Nuclear Materials，2009，386-388：122-126．

［10］ Yican Wu，F(xiàn)DS Team．Conceptual Design and testing strategy of a dual function lithium-lead test blanket module in ITER and EAST［J］．Nucl Fusion，2007，47（11）：1533-1539．

［11］ Yican Wu，the FDS Team．Design Analysis of the China Dual-functional Lithium Lead （DFLL） Test Blanket Module in ITER［J］．Fusion Engineering and Design，2007，82：1893-1903．

［12］ Poitevin Y，Boccaccini L V，Cardella A，et al．The European breeding blankets development and the test strategy in ITER［J］．Fusion Eng Des，2005，75-79：741-749．

［13］ 邢朝青．鈉的凈化技術(shù)［J］．核動力工程，1989，10（1）：75-82．

［14］ 吳宜燦，黃群英，朱志強(qiáng)，等．中國系列液態(tài)鋰鉛實驗回路設(shè)計與研制進(jìn)展［J］．核能科學(xué)與工程，2009，29（2）：161-169．

［15］ Zhiqiang Zhu，Qunying Huang，Sheng Gao，et al．Design analysis of DRAGON-Ⅳ LiPb loop［J］．Fusion Eng Des，2011，86：2666-2669．

［16］ Yican Wu，Qunying Huang，Zhiqing Zhu，et al．R＆D of DRAGON Series Lithium-Lead Loops for Material＆Blanket Technology Testing in China［R］．ICFRM-15 meeting，South Carolina，America，October 16-22，2011．

［17］ 洪順章．鈉凈化裝置和鈉化學(xué)回路運(yùn)行經(jīng)驗［J］．原子能科學(xué)技術(shù)，1991，25（5）：90-96．

［18］ 趙福宇，付制浩，張貴勤．鈉冷阱過濾芯布置的最優(yōu)化分析［J］．核動力工程，1993，14（6）：533-537．

［19］ 汪衛(wèi)華，朱志強(qiáng)，李晉嶺，等．液態(tài)金屬回路DRAGON-Ⅱ鋰鉛合金流動與傳熱三維數(shù)值模擬［J］．核科學(xué)與工程，2010，30（2）：166-171．

［20］ Terlain A，F(xiàn)lament T，Dufrenoy T，et al．Mass transfer of 316L steel in Pb17Li CLIPPER loop after 12000 hour running［J］．J Nucl Mater，1992，191-194：984-987．

［21］ Chopra O K，Tortorelli P F．Compatibility of materials for use in liquid metal blanlets of fusion reactors［J］．J Nucl Mater，1984，122-123：1201-1212．

［22］ 單麗君，陳蕊娜．固液分離過濾器流場的數(shù)值模擬與分析［J］．大連交通大學(xué)學(xué)報，2009，30（2），6-9．

［23］ 張李偉．船用絲網(wǎng)除霧器分離效率計算［J］．船舶工程，2007，29（5），1-8．