核電作為一種高效的清潔能源受到了能源行業(yè)的青睞，并逐漸成為一種重要的能源生產(chǎn)方式。壓力容器中的堆芯是熱能產(chǎn)生的部位，也是反應(yīng)堆的關(guān)鍵部件之一，其熱工水力特性直接關(guān)系到燃料元件的冷卻效率及安全問題，因此重點(diǎn)研究方向之一。在極端情況下，燃料組件子通道可能因?yàn)榘羰[脹、異物堆積等發(fā)生通道堵流現(xiàn)象，發(fā)生傳熱惡化，甚至超溫、堆芯融化等危險(xiǎn)工況，因此有必要開展相關(guān)研究。采用實(shí)驗(yàn)與數(shù)值模擬手段測(cè)量、模擬堵流情況下的子通道流場(chǎng)分布、獲取子通道流量分布特性可以充分發(fā)揮實(shí)驗(yàn)與模擬優(yōu)勢(shì)。本文利用激光粒子測(cè)速方法測(cè)量了5×5棒束堵塞情況下子通道的流場(chǎng)分布特性以及子通道流量分布特性，并基于國(guó)內(nèi)首款核工業(yè)多質(zhì)多相流體動(dòng)力學(xué)計(jì)算分析軟件“麒麟”軟件和商用CFD軟件進(jìn)行了堵流現(xiàn)象的數(shù)值模擬，結(jié)合實(shí)驗(yàn)結(jié)果分析棒束子通道堵流特性。

1 實(shí)驗(yàn)裝置簡(jiǎn)介
本次研究首先采用可視化試驗(yàn)方法對(duì)堵流情況下5×5棒束組件子通道流場(chǎng)進(jìn)行了測(cè)量，實(shí)驗(yàn)回路如圖1所示，主要包括實(shí)驗(yàn)本體、水箱、主泵、流量計(jì)、閥門和測(cè)控儀表等。

圖1 棒束組件內(nèi)流場(chǎng)分布試驗(yàn)回路流程圖

圖2a所示為實(shí)驗(yàn)本體結(jié)構(gòu)示意圖，本體主體為四面可視的矩形通道結(jié)構(gòu)，進(jìn)出口有穩(wěn)流裝置，內(nèi)部放置5×5棒束組件。組件由支撐段和測(cè)量段組成，支撐段為硬質(zhì)材料，測(cè)量段為薄壁FEP管，以最大程度接近水的折射率。堵流件安裝在FEP管外，堵流區(qū)域如圖2c所示為3×3棒束區(qū)域，中間通道堵流比例最大，5×5棒束截面的總堵流比例為72%。堵流件安裝位置距入口大于40D，距出口距離也在20D以上，以最大程度減弱入口和出口效應(yīng)。測(cè)量區(qū)域約為20D。

圖2 實(shí)驗(yàn)本體及棒束組件結(jié)構(gòu)示意圖

如圖2c所示，棒束組件由5×5棒組成，測(cè)量段為外徑9.5mm，中心距為12.6mm的FEP管。PIV面激光以XZ平面照射測(cè)量區(qū)域，高速相機(jī)則沿Y方向拍攝。采用跨幀模式進(jìn)行圖像采集，然后使用專業(yè)軟件進(jìn)行速度場(chǎng)計(jì)算，即獲得了測(cè)量區(qū)域的流場(chǎng)分布。相機(jī)拍攝幀率為15Hz,每次針對(duì)同一工況拍攝200組以上照片，然后求取平均流場(chǎng)。

2 數(shù)值模擬方法
數(shù)值模擬基于核工業(yè)多質(zhì)多相流體動(dòng)力學(xué)計(jì)算分析軟件“麒麟”進(jìn)行，該軟件采用浸入界面方法進(jìn)行網(wǎng)格劃分，基于NS方程求解流場(chǎng)、溫度場(chǎng)、壓力場(chǎng)等流體力學(xué)問題。該軟件支持常見CAD軟件三維幾何模型的導(dǎo)入然后進(jìn)行IST方法的笛卡爾網(wǎng)格劃分，之后進(jìn)行邊界條件設(shè)置、模型設(shè)置、求解設(shè)置等，最后進(jìn)行迭代求解和后處理。該流場(chǎng)計(jì)算軟件可以基于Windows和Linux系統(tǒng)進(jìn)行，二者求解方式?jīng)]有任何區(qū)別，本次求解基于Windows系統(tǒng)進(jìn)行。
圖3所示為本次模擬的三維幾何和網(wǎng)格劃分示意圖。圖3a所示為浸入界面法的基本原理，即在網(wǎng)格劃分時(shí)不區(qū)分流體與固體域，使用流固界面符號(hào)距離函數(shù)Φs來(lái)數(shù)值地區(qū)分流體與固體域，在求解NS方程時(shí)，通過(guò)引入Heaviside函數(shù)

自動(dòng)地區(qū)別流體與固體域。該函數(shù)在流體內(nèi)等于1，在固體內(nèi)等于0，在界面上等于0.5。帶堵流件的棒束幾何由三維CAD軟件根據(jù)實(shí)驗(yàn)裝置尺寸測(cè)量結(jié)果建立，然后導(dǎo)入“麒麟”軟件中進(jìn)行網(wǎng)格劃分，圖3b固體域即是繪制的棒束三維幾何。劃分網(wǎng)格時(shí)直接根據(jù)試驗(yàn)裝置矩形流道大小繪制矩形流體域，即圖3c中背景網(wǎng)格，這樣三維幾何就“浸沒”在背景網(wǎng)格中（流體域）。之后基于重疊網(wǎng)格技術(shù)在堵流件區(qū)域添加加密的子網(wǎng)格對(duì)堵流塊進(jìn)行加密，以提高堵流件區(qū)域流場(chǎng)解析精度。本次模擬主要針對(duì)堵流件上下游總計(jì)約9D的范圍，為了避免入口效應(yīng)，也繪制了堵流件上游10D的網(wǎng)格。網(wǎng)格無(wú)關(guān)性分析表明使用總量約50萬(wàn)的網(wǎng)格可以實(shí)現(xiàn)網(wǎng)格無(wú)關(guān)。

圖3 采用浸入界面法時(shí)網(wǎng)格劃分示意圖

本次模擬針對(duì)實(shí)驗(yàn)工況，因此所涉及的模型為冷態(tài)流動(dòng)，模擬為穩(wěn)態(tài)模擬，主要求解X、Y、Z三個(gè)方向的速度以及壓力。模擬邊界條件為流量入口，流量、壓力、溫度根據(jù)實(shí)驗(yàn)測(cè)量結(jié)果設(shè)置，模擬中迭代計(jì)算到殘差小于1e-5或者不再減小為止。

3 結(jié)果及分析
3.1 流場(chǎng)測(cè)量結(jié)果
本次實(shí)驗(yàn)中維持進(jìn)口流量穩(wěn)定在約14t/h，對(duì)應(yīng)的水平截面平均質(zhì)量流速約1566 kg/m2/s。實(shí)驗(yàn)中測(cè)量了堵流件上下游約7D范圍內(nèi)棒束中心豎直截面、子通道豎直中心截面的流場(chǎng)。測(cè)量得到的典型豎直截面流場(chǎng)圖如所示，圖中示出了C1和C2通道在不安裝和安裝堵流件時(shí)的上下游流場(chǎng)測(cè)量結(jié)果。在沒有堵流件的正常狀態(tài)下，堵流件所在位置上下游流場(chǎng)基本都以本子通道的豎直流動(dòng)為主，流場(chǎng)未發(fā)生偏斜，不同子通道流速分布相似。在堵流情況下，C1和C2都出現(xiàn)了流場(chǎng)偏斜和流速不均的情況。對(duì)照?qǐng)D1可知C1子通道豎直截面僅有一半流場(chǎng)被堵流件遮擋，因此可見其流場(chǎng)所受影響較小，僅堵塊下游少量流場(chǎng)受影響。對(duì)于C2，被堵塞的通道區(qū)域較寬，因此可以發(fā)現(xiàn)橫向位置1和2所在區(qū)域堵流件下游出現(xiàn)了流速很低的區(qū)域，未堵塞子通道流速則出現(xiàn)明顯增大，表明堵塞對(duì)流場(chǎng)影響很大。

圖4 不同狀態(tài)下堵流件上下游豎直截面流場(chǎng)典型結(jié)果

為了直觀地分析堵流對(duì)子通道流速的影響，這里使用速度不均勻指數(shù)來(lái)量化。該指數(shù)計(jì)算為求某豎直截面高度Z處絕對(duì)速度值與截面平均速度的均方根。圖5所示為根據(jù)實(shí)驗(yàn)測(cè)量結(jié)果得到的通道C1、C2、C3、C4在不同高度處的速度不均勻指數(shù)，其中曲線所在位置為堵流件，軸向位置0為堵流件軸向中心位置。可以發(fā)現(xiàn)通道C1和C4的速度不均勻指數(shù)變化較小，僅在虛線所示堵流件區(qū)域出現(xiàn)了較小的增加，并很快恢復(fù)正常，表明其受堵流件的影響較小，冷卻劑可以正常流動(dòng)。而通道C2和C3速度不均勻指數(shù)在堵流件處出現(xiàn)了較大的增加，而后緩慢減小；并且速度不均勻指數(shù)在堵流件上游1D處已經(jīng)開始增加，表明其對(duì)流場(chǎng)的影響在上游1D處已經(jīng)開始。雖然在堵流件下游約3D以后速度不均勻指數(shù)下降較緩，但仍未恢復(fù)到堵流前水平。上述結(jié)果表明堵流件對(duì)流速產(chǎn)生了較大影響，實(shí)驗(yàn)結(jié)果表明在堵流通道，在截面平均質(zhì)量流速約1500kg/m2/s下，堵流通道下游2D距離處內(nèi)堵塞子通道平均流速僅為截面平均流速的60%；5D距離處恢復(fù)到截面平均流速約85%。

圖5 速度不均勻指數(shù)隨Z變化曲線

3.2 模擬結(jié)果及分析
使用前述數(shù)值模擬方法進(jìn)行數(shù)值迭代求解可以得到5×5棒在堵流情況下的模擬結(jié)果，本次研究也使用商業(yè)CFD軟件進(jìn)行了模擬，二者網(wǎng)格總數(shù)、湍流模型、邊界條件等相同，以對(duì)比數(shù)值分析軟件與商業(yè)CFD軟件。圖6所示為C1、C2豎直截面的流場(chǎng)圖，圖中灰色框部分為堵流件所在位置。對(duì)于C1，子通道堵塞比例為10%~20%，小于C2通道的20~40%，因此流場(chǎng)受到的影響較小。C2通道在堵流件區(qū)域因子通道堵塞出現(xiàn)了明顯的流速降低，在堵流件下游形成了低速區(qū)，這與試驗(yàn)測(cè)量結(jié)果吻合。比較“麒麟”軟件和商用CFD軟件結(jié)果，可以發(fā)現(xiàn)二者流場(chǎng)分布十分接近。

圖6 模擬所得堵流情況下豎直截面流場(chǎng)典型結(jié)果

為了更加直觀地對(duì)比“麒麟”軟件和商用CFD軟件差別，提取了堵流件下游不同距離處試驗(yàn)與模擬結(jié)果繪制成圖7所示曲線圖。圖中分別對(duì)比了堵流件下游0.5D和2D距離處實(shí)驗(yàn)與模擬結(jié)果，可以看出在兩種距離處實(shí)驗(yàn)與模擬值的匹配性很好，同時(shí)也可以看出“麒麟”軟件的模擬結(jié)果與商用CFD軟件的結(jié)果十分接近，表明其流場(chǎng)模擬能力與商用CFD軟件接近。

圖8計(jì)算了模擬結(jié)果與實(shí)驗(yàn)結(jié)果的均方根誤差，圖中橫坐標(biāo)表示下游計(jì)算位置與堵流件位置的距離。可以看出，在堵流件2D范圍內(nèi)模擬值與實(shí)驗(yàn)值的均方根誤差比更遠(yuǎn)距離的明顯大，表明模擬與試驗(yàn)偏差比其他區(qū)域大，這可能是采用RANS湍流模型的導(dǎo)致堵流件緊鄰流場(chǎng)的強(qiáng)湍流區(qū)域解析能力不足導(dǎo)致的。整體而言，“麒麟”軟件的均方根誤與商業(yè)CFD軟件在走勢(shì)和值上十分接近。

圖8 模擬均方根誤差對(duì)比圖

4結(jié)論

本文基于棒束子通道流場(chǎng)可視化實(shí)驗(yàn)裝置和5×5棒束子通道堵流實(shí)驗(yàn)本體獲得了子通道堵流72%比例下的上下游流場(chǎng)分布，并使用“麒麟”軟件對(duì)堵流情況下的子通道流場(chǎng)進(jìn)行了數(shù)值模擬分析。本文獲得的主要結(jié)論為：子通道堵流約72%情況下，堵塞子通道上游約1D范圍流場(chǎng)開始受到明顯影響；堵流通道下游2D距離處內(nèi)堵塞子通道平均流速僅為截面平均流速的60%；5D距離處恢復(fù)到截面平均流速約85%。基于浸入界面方法的“麒麟”軟件能夠?qū)Χ铝髑闆r下的棒束子通道流場(chǎng)進(jìn)行合理模擬，其模擬結(jié)果與商用CFD軟件相當(dāng)。

“麒麟KILI”是中國(guó)核動(dòng)力研究設(shè)計(jì)院研發(fā)的國(guó)內(nèi)首款核工業(yè)多質(zhì)多相流體動(dòng)力學(xué)計(jì)算分析軟件，該軟件填補(bǔ)了國(guó)內(nèi)核工業(yè)自主研發(fā)計(jì)算流體動(dòng)力學(xué)軟件的空白，實(shí)現(xiàn)了關(guān)鍵軟件技術(shù)的自主可控。上海積鼎信息科技有限公司為麒麟KILI軟件開發(fā)唯一合作單位。

審核編輯 黃宇