定位格架下游湍流特性測(cè)量研究

曲文海，陳仕龍，黃慧劍，熊進(jìn)標(biāo)

(1.上海交通大學(xué) 核科學(xué)與工程學(xué)院，上海 200240；2.中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 成都 610213)

典型壓水堆核燃料組件中安裝了數(shù)個(gè)帶交混翼的定位格架，既能支撐棒束燃料組件，亦能增強(qiáng)子通道內(nèi)部和子通道之間的二次流和湍流脈動(dòng)。子通道程序亟需精確的子通道交混模型[1-2]。CFD模擬是燃料組件設(shè)計(jì)的重要方法之一，可用于研究格架交混效應(yīng)，但仍需大量實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證[3-5]。

Xiong等[6-7]使用激光多普勒測(cè)速儀測(cè)量了3×3和6×6棒束內(nèi)的三維流場(chǎng)，并討論了雷諾數(shù)效應(yīng)。Wang等[8-9]在水和氟化乙烯丙烯共聚物(FEP)折射率補(bǔ)償技術(shù)基礎(chǔ)上測(cè)量了光棒束內(nèi)軸向流場(chǎng)，并用CFD方法模擬棒束通道內(nèi)流場(chǎng)。以上研究揭示了光棒束通道內(nèi)的湍流特征。Nguyen等[10-12]使用基于水和FEP折射率補(bǔ)償技術(shù)的三維粒子成像測(cè)速儀(PIV)測(cè)量了5×5棒束通道內(nèi)定位格架下游的流場(chǎng)，獲得一個(gè)縱截面上三維速度場(chǎng)和湍流特征。Chang等[13-14]使用激光多普勒測(cè)速儀(LDV)測(cè)量了分流型交混翼和旋流型交混翼定位格架下游橫向速度場(chǎng)，獲得了交混翼產(chǎn)生的渦結(jié)構(gòu)特征和湍流脈動(dòng)特征。Li[15]使用水和FEP折射率補(bǔ)償技術(shù)輔助PIV測(cè)量定位格架下游的軸向和橫向流場(chǎng)，獲得了格架下游流場(chǎng)特征。Kim等[16]使用PIV方法和基于有機(jī)玻璃與碘化鈉水溶液的折射率補(bǔ)償技術(shù)，測(cè)量了4×6棒束通道內(nèi)的縱向速度場(chǎng)。為從光路上消除常規(guī)PIV測(cè)量中存在的視角誤差，Xiong等[17]使用遠(yuǎn)心鏡頭測(cè)量格架下游橫向流場(chǎng)。Qu等[18-20]對(duì)基于FEP和水折射率補(bǔ)償技術(shù)的二維PIV方法，進(jìn)行了系統(tǒng)的誤差和不確定度分析，在此基礎(chǔ)上獲得了定位格架下游的橫向流動(dòng)和湍流沿程變化特征以及雷諾數(shù)效應(yīng)，并開(kāi)發(fā)了滑移式壓降測(cè)量系統(tǒng)，獲得了子通道橫向壓降。Xiong等[21]基于流場(chǎng)和壓降實(shí)驗(yàn)數(shù)據(jù)，開(kāi)發(fā)了子通道交混模型和阻力模型。

在上述研究基礎(chǔ)上，本文使用基于FEP和水折射率補(bǔ)償技術(shù)的二維PIV方法，獲得不同雷諾數(shù)下棒束通道內(nèi)定位格架下游多個(gè)橫截面上的高精度精細(xì)化橫向流場(chǎng)實(shí)驗(yàn)數(shù)據(jù)，分析定位格架下游橫向湍流特征和雷諾數(shù)效應(yīng)。

1 實(shí)驗(yàn)裝置

本實(shí)驗(yàn)在上海交通大學(xué)中型熱工水力回路上開(kāi)展，該回路包括一個(gè)實(shí)驗(yàn)回路和輔助的冷卻回路。5×5棒束實(shí)驗(yàn)段和PIV布置如圖1所示，棒直徑D=9.5 mm，棒間距P=12.6 mm，子通道水力直徑Dh=11.78 mm。三維坐標(biāo)系定義如圖1c、d所示，坐標(biāo)原點(diǎn)位于格架條帶上平面的中心棒圓心，z軸與中心棒軸線重合，正方向沿著主流方向，x軸與y軸沿水平方向。棒束上安裝了1個(gè)交混翼傾角30°的定位格架，定位格架交混翼的設(shè)計(jì)參考了文獻(xiàn)[22]中的交混翼設(shè)計(jì)思想，尺寸如圖2所示。PIV系統(tǒng)包括1臺(tái)30 mJ高頻激光器(最高頻率為10 kHz)，1臺(tái)1 280×800像素高速相機(jī)(最高頻率為7 400 Hz，像素尺寸dr=20 μm)，1臺(tái)200 mm定焦鏡頭，1臺(tái)三維位移平臺(tái)(行程1 000 mm、1 000 mm、600 mm，精度0.01 mm)，1臺(tái)同步器和1套PIV圖像處理軟件。示蹤粒子為直徑dp=20 μm的聚苯乙烯基吡啶樹(shù)脂(PSP)粒子。實(shí)驗(yàn)在暗室中進(jìn)行，橫向速度測(cè)量時(shí)光學(xué)系統(tǒng)圖像放大倍數(shù)M=0.6，激光面厚度Δz=3 mm?？v向速度測(cè)量時(shí)光學(xué)系統(tǒng)圖像放大倍數(shù)M=0.4，激光面厚度Δz=1 mm。橫向測(cè)量面位于z=1Dh～28Dh范圍，覆蓋中心棒周圍的4個(gè)中心子通道，橫向空間分辨率為0.275 mm?？v向測(cè)量面位于y=-1.5P截面上，覆蓋z=1Dh～25Dh的范圍，空間分辨率為0.32 mm。實(shí)驗(yàn)段進(jìn)口壓力為0.15 MPa，溫度為25 ℃，4個(gè)試驗(yàn)工況的主流平均速度Wb分別為0.5、1、2和3 m/s，對(duì)應(yīng)子通道雷諾數(shù)分別為6 600、13 200、26 400、39 600。定義三維坐標(biāo)系原點(diǎn)于格架條帶頂部平面的中心棒圓心處。

a——橫截面測(cè)量PIV布置；b——縱截面測(cè)量PIV布置；c——縱截面位置；d——格架局部；e——棒束橫截面

a——格架尺寸；b——交混翼尺寸

2 實(shí)驗(yàn)結(jié)果分析

三維平均速度、脈動(dòng)速度均方根定義如下：

(1)

(2)

(3)

(4)

(5)

(6)

式中：U、V、W分別為x、y、z方向的平均速度；ui、vi、wi分別為x、y、z方向的瞬時(shí)速度；N為采樣數(shù)量；Urms、Vrms、Wrms分別為x、y、z方向的脈動(dòng)速度均方根。對(duì)平均速度和脈動(dòng)速度均方根以Wb進(jìn)行歸一化。

2.1 誤差和不確定度分析

系統(tǒng)誤差主要為視角誤差和當(dāng)?shù)厮俣忍荻日`差[18]，隨機(jī)誤差主要為粒子出平面運(yùn)動(dòng)誤差，湍流統(tǒng)計(jì)量還存在采樣數(shù)量誤差。糾正系統(tǒng)誤差后，湍流統(tǒng)計(jì)量不確定度由粒子出平面運(yùn)動(dòng)誤差和采樣數(shù)量誤差決定。平均速度和湍流脈動(dòng)速度均方根的測(cè)量結(jié)果最大不確定度小于1%的主流平均速度。

視角誤差εp的計(jì)算公式為：

(7)

式中：|Vout|為垂直于測(cè)量面的速度分量；z0為相機(jī)的物距。

粒子出平面運(yùn)動(dòng)誤差εout的計(jì)算公式為：

(8)

式中：Δt為激光脈沖間隔時(shí)間，即瞬時(shí)速度測(cè)量時(shí)間；Δz為激光厚度。

平均速度和脈動(dòng)速度均方根的采樣數(shù)量誤差εN，V、εN,Vrms的計(jì)算公式為：

(9)

(10)

為使采樣數(shù)量誤差降低到1%Wb以內(nèi)，采樣數(shù)量取為N=5 000。經(jīng)過(guò)系統(tǒng)誤差糾正，平均速度不確定度εV由采樣數(shù)量誤差決定，脈動(dòng)速度均方根不確定度εVrms由粒子出平面運(yùn)動(dòng)誤差和采樣數(shù)量誤差決定：

(11)

(12)

2.2 橫向流場(chǎng)分析

Re=39 600工況下，定位格架下游橫向流場(chǎng)如圖3所示。在z=1Dh處，1對(duì)分流型交混翼在子通道內(nèi)產(chǎn)生了1對(duì)同轉(zhuǎn)漩渦對(duì)和繞棒橫向速度，U峰值約為40%Wb，V峰值約為30%Wb。隨著流動(dòng)發(fā)展，子通道內(nèi)1對(duì)漩渦逐漸耗散擴(kuò)大融合，在z=5Dh處形成1個(gè)覆蓋子通道的大漩渦，U峰值約為25%Wb，V峰值約為28%Wb。在z=10Dh處，U峰值約為15%Wb，V峰值約為12%Wb。該漩渦在沿程發(fā)展中逐漸減弱消失。在z=28Dh處，U、V峰值約為5%Wb，橫向速度分布趨于光棒束充分發(fā)展。

a——z=1Dh；b——z=5Dh；c——z=10Dh；d——z=28Dh

縱截面y=-1.5P上的橫向速度如圖4a所示，相鄰兩個(gè)子通道內(nèi)正交布置的交混翼片產(chǎn)生方向相反的橫向速度。最大橫向速度出現(xiàn)棒間隙處，子通道中心附近橫向速度較小，橫向速度向下游逐漸衰減。在z=1Dh～5Dh區(qū)域，橫向速度大小顯示出子通道內(nèi)雙渦結(jié)構(gòu)轉(zhuǎn)變?yōu)閱螠u結(jié)構(gòu)的現(xiàn)象，隨后逐漸耗散。縱截面y=-1.5P上的橫向脈動(dòng)速度均方根如圖4b所示，在0～1Dh區(qū)域，橫向脈動(dòng)速度均方根分布與交混翼分布一致。在z=1Dh處，橫向脈動(dòng)速度均方根分布與二次流雙渦結(jié)構(gòu)一致，Urms有兩個(gè)峰值出現(xiàn)在子通道內(nèi)，約為32%Wb。在z=1Dh～5Dh區(qū)域，橫向脈動(dòng)速度均方根隨著子通道內(nèi)雙渦結(jié)構(gòu)轉(zhuǎn)變?yōu)閱螠u結(jié)構(gòu)而迅速耗散，在z=2Dh處，Urms峰值約為35%Wb；在z=5Dh處，Urms峰值約為15%Wb。隨后橫向脈動(dòng)速度均方根緩慢耗散。

圖4 Re=39 600時(shí)縱截面上的橫向速度(a)與橫向脈動(dòng)速度均方根(b)

橫向脈動(dòng)速度均方根分布如圖5、6所示，分布模式與渦結(jié)構(gòu)分布一致。在z=1Dh處，橫向脈動(dòng)速度均方根峰值出現(xiàn)在雙渦結(jié)構(gòu)處，因?yàn)闇u結(jié)構(gòu)內(nèi)強(qiáng)剪切流動(dòng)產(chǎn)生了強(qiáng)湍流脈動(dòng)，Urms峰值約為20%Wb，Vrms峰值約為25%Wb。隨著雙渦結(jié)構(gòu)耗散，橫向湍流脈動(dòng)急劇降低。在z=5Dh處，單渦結(jié)構(gòu)在子通道中心產(chǎn)生強(qiáng)湍流脈動(dòng)，并逐漸耗散，Urms峰值約為12.5%Wb，Vrms峰值約為12%Wb。橫向脈動(dòng)速度分量呈現(xiàn)出與交混翼分布一致的棋盤式分布。在z=28Dh處，Urms峰值約為4%Wb，Vrms峰值約為3%Wb，橫向脈動(dòng)速度趨于光棒束充分發(fā)展分布。

a——z=1Dh；b——z=5Dh；c——z=10Dh；d——z=28Dh

a——z=1Dh；b——z=5Dh；c——z=10Dh；d——z=28Dh

2.3 軸向流場(chǎng)分析

定位格架下游1Dh～12Dh的軸向平均速度如圖7a所示。在格架近場(chǎng)，交混翼顯著改變了軸向速度分布。在x∈[0,1]范圍子通道內(nèi)，軸向速度在交混翼根部受到指向x軸的交混翼片阻礙顯著減小，隨后沿著交混翼片方向逐漸增大，在子通道中心(x=0.5P)達(dá)到峰值。在x∈[-1,0]范圍子通道內(nèi)，由于交混翼片指向y軸方向，軸向速度未受到交混翼片的阻礙而在子通道中心達(dá)到峰值，軸向速度在棒間隙(x=0P)處受交混翼根部阻礙出現(xiàn)谷值。在z=1Dh處，軸向速度峰值位于子通道中心附近，約為117%Wb，最小值位于棒間隙附近，約為85%Wb；在z=2Dh處，軸向速度峰值位于子通道內(nèi)，約為108%Wb，最小值位于棒間隙附近，約為90%Wb；在z=5Dh處，軸向速度峰值位于子通道中心附近，約為110%Wb，谷值位于棒間隙附近，約為95%Wb；在z=10Dh處，軸向速度峰值位于子通道中心附近，約為105%Wb。

圖7 y=-1.5P、Re=39 600的歸一化軸向平均速度(a)和脈動(dòng)速度均方根(b)

定位格架下游1Dh～12Dh的軸向脈動(dòng)速度均方根分布如圖7b所示。在格架近場(chǎng)，軸向脈動(dòng)速度均方根分布與交混翼方向一致，并在棒間隙附近達(dá)到峰值，在子通道中心附近的值較小。隨著流動(dòng)發(fā)展，軸向脈動(dòng)速度均方根迅速衰減。在z=1Dh處，軸向脈動(dòng)速度均方根峰值在棒間隙處，約為20%Wb，谷值位于子通道中心，約為15%Wb。在z=1Dh～5Dh范圍，軸向脈動(dòng)速度均方根迅速耗散，隨后緩慢耗散。在z=2Dh處，軸向脈動(dòng)速度均方根峰值位于棒間隙處，約為18%Wb，谷值位于子通道內(nèi)，約為12%Wb；在z=5Dh處，軸向脈動(dòng)速度均方根峰值位于棒間隙處，約為15%Wb，谷值位于子通道中心，約為10%Wb；在z=10Dh處，軸向脈動(dòng)速度均方根峰值位于棒間隙處，約為15%Wb，谷值位于子通道中心，約8%Wb。

2.4 雷諾數(shù)效應(yīng)

雷諾數(shù)效應(yīng)如圖8、9所示，在格架近場(chǎng)，格架交混效應(yīng)強(qiáng)于雷諾數(shù)效應(yīng)，掩蓋了雷諾數(shù)對(duì)橫向速度的影響。隨著流動(dòng)發(fā)展，格架交混效應(yīng)衰減，雷諾數(shù)效應(yīng)逐漸表現(xiàn)出來(lái)?？傮w上隨雷諾數(shù)的增加，橫向平均速度增大，橫向脈動(dòng)速度均方根減小，并趨于雷諾數(shù)無(wú)關(guān)的自模化區(qū)。雷諾數(shù)效應(yīng)對(duì)脈動(dòng)速度均方根的作用更明顯。如圖8所示，在z=1Dh處，各雷諾數(shù)下橫向平均速度相差無(wú)幾；在z=5Dh處，橫向平均速度隨著雷諾數(shù)增大而增大，Re=6 600時(shí)最大橫向平均速度為20%Wb，而Re=39 600時(shí)最大橫向平均速度為23%Wb；在z=10Dh處，橫向平均速度隨著雷諾數(shù)增大而增大，Re=6 600時(shí)最大橫向平均速度為10%Wb，而Re=39 600時(shí)最大橫向平均速度為16%Wb；在z=28Dh處，橫向平均速度隨著雷諾數(shù)增大而增大，Re=6 600時(shí)最大橫向平均速度為5%Wb，Re=39 600時(shí)最大橫向平均速度為7.5%Wb。如圖9所示，在子通道中心線y=0.5P上，在z=1Dh、5Dh、10Dh處，橫向脈動(dòng)速度均方根無(wú)明顯的雷諾數(shù)效應(yīng)。在z=28Dh處，橫向脈動(dòng)速度均方根隨著雷諾數(shù)的增大而減小。

a——z=1Dh；b——z=5Dh；c——z=10Dh；d——z=28Dh

a——z=1Dh；b——z=5Dh；c——z=10Dh；d——z=28Dh

本文實(shí)驗(yàn)獲得了豐富的橫向流場(chǎng)和湍流脈動(dòng)場(chǎng)的實(shí)驗(yàn)數(shù)據(jù)，能夠解釋定位格強(qiáng)化子通道交混的機(jī)理，能夠用于開(kāi)發(fā)精細(xì)的子通道交混模型，為子通道分析程序提供基于實(shí)驗(yàn)數(shù)據(jù)庫(kù)的交混模型，提高子通道分析程序的預(yù)測(cè)能力和精度。

3 結(jié)論

本文使用基于水和FEP折射率補(bǔ)償技術(shù)的二維PIV方法，在4個(gè)不同雷諾數(shù)下測(cè)量帶交混翼格架的5×5棒束通道內(nèi)橫向流場(chǎng)，得到以下結(jié)論：1)定位格架下游二次流結(jié)構(gòu)經(jīng)歷了雙渦結(jié)構(gòu)向單渦結(jié)構(gòu)的轉(zhuǎn)變及單渦結(jié)構(gòu)衰減過(guò)程，湍流統(tǒng)計(jì)量沿程變化均受渦結(jié)構(gòu)演進(jìn)影響；2)橫向平均速度和脈動(dòng)速度均方根受雷諾數(shù)效應(yīng)和格架交混效應(yīng)影響，格架近場(chǎng)以格架交混效應(yīng)影響為主，格架遠(yuǎn)場(chǎng)以雷諾數(shù)效應(yīng)影響為主。