許志紅，王 喆，楊 萍，史國寶

（上海核工程研究設計院，上海 200233）

?

燃料組件邊角柵元偏離泡核沸騰比分析

許志紅，王 喆，楊 萍，史國寶

（上海核工程研究設計院，上海 200233）

摘要：燃料組件邊角柵元的阻力系數(shù)大于典型柵元和冷壁柵元。本文基于標準子通道模型，通過調整子通道模擬區(qū)域的位置，修改了通道、燃料棒以及格架阻力等相關參數(shù)，直接模擬燃料組件邊角柵元的綜合效應，進行偏離泡核沸騰比（Departure from Nuclear Boiling Ratio，簡稱DNBR）分析。分析結果表明：邊角柵元的偏離泡核沸騰比能被典型和冷壁柵元的計算結果所包絡。本文同時對燃料組件間隙和邊角柵元的偏離泡核沸騰比的影響進行了進一步分析。

關鍵詞：燃料組件；邊角柵元；偏離泡核沸騰比

反應堆熱工設計的目的主要是為了提供與堆芯產熱量相匹配的恰當?shù)膫鳠崮芰?，保持堆芯的完整性。其中偏離泡核沸騰（Departure from Nuclear Boiling，簡稱DNB）設計基準是熱工水力設計中的一項重要準則。即在95%的置信度水平下，在I、II類工況下堆芯極限燃料棒至少有95%的概率不發(fā)生偏離泡核沸騰。如果在事故過程中堆芯的最小偏離泡核沸騰比大于安全分析限值，則認為此類事故滿足偏離泡核沸騰設計準則的要求［1］。

目前的反應堆熱工水力分析一般采用參數(shù)統(tǒng)計學方法開展偏離泡核沸騰比計算［1 -5］，相關熱工水力參數(shù)取名義值。有多種因素會影響偏離泡核沸騰比的計算結果，其中在子通道分析中的格架阻力系數(shù)作為燃料組件形阻的最重要部分對計算結果影響較大。多種類型的水冷堆均開展過考慮格架阻力系數(shù)的熱工水力分析［6 -10］。為了分析燃料組件邊角柵元阻力系數(shù)的影響，本文進行了子通道分析。

燃料組件邊角柵元（兩個組件之間的燃料棒和間隙形成的子通道）的阻力系數(shù)與典型柵元（燃料組件中四根燃料棒圍成的子通道）和冷壁柵元（燃料組件中燃料棒和導向管圍成的子通道）的阻力系數(shù)是有差別的。燃料組件邊角柵元的阻力系數(shù)大于典型柵元和冷壁柵元的阻力系數(shù)。阻力系數(shù)大是不利的，它可能會導致邊角柵元的最小偏離泡核沸騰比小于典型柵元和冷壁柵元的偏離泡核沸騰比。盡管燃料組件邊角柵元的阻力系數(shù)大于典型柵元和冷壁柵元的阻力系數(shù)，這會導致流動阻力增大，由此降低通道中的流量。但是，邊角柵元中無導向管，這減少了冷壁效應，而且邊角柵元流通面積大于典型柵元和冷壁柵元，綜合效應可能使得流經邊角柵元的冷卻劑流量相比典型柵元和冷壁柵元要大。

這是一個共性的問題，此前并未進行過詳細分析。一般燃料組件格架的阻力系數(shù)都來自于試驗數(shù)據(jù)，這些數(shù)據(jù)是基于整個格架的，試驗獲取的數(shù)據(jù)并不區(qū)分格架的所在的柵元。通常的子通道分析并未將邊、角柵元與典型、冷壁柵元進行區(qū)分，一般進行子通道分析時邊、角柵元仍然按照是否包含導向管而采用典型或冷壁柵元的阻力系數(shù)。

AP1000核電廠基于試驗獲取的格架損失系數(shù)，通過經驗關系式，推算出更為精細的邊、角柵元和典型、冷壁柵元子通道阻力系數(shù)。AP1000核電廠典型的燃料組件包含有8個中間格架和4個攪混（Intermediate Flaw Mixing，簡稱IFM）格架。中間格架有兩種典型柵元，區(qū)別是柵元中彈簧和剛凸的數(shù)目不同（部分典型柵元包含有雙剛凸條帶）；攪混格架中僅有剛凸，而沒有彈簧，因此攪混格架的典型柵元只有一種。中間格架有四種冷壁柵元，區(qū)別也是柵元中彈簧和剛凸的數(shù)目不同。這主要是由于柵元相對于導向管的位置不同以及雙剛凸條帶位置不同引起的；攪混格架的冷壁柵元只有一種。中間格架有三種邊柵元，一種是因為包含有雙剛凸條帶；另外兩種是由于外條帶導向翼結構不同導致的柵元結構不同。攪混格架有兩種邊柵元，區(qū)別是外條帶上的頂部導向翼分為短導向翼和長導向翼兩種。中間格架和攪混格架的角柵元都只有一種。根據(jù)典型、冷壁、邊、角柵元各自特點，結合中間和攪混格架同一類型柵元的種類進行加權計算，獲取各種類型柵元的阻力系數(shù)。典型、冷壁柵元與邊、角柵元的結構和堵塞面積有所不同，邊角柵元位于燃料組件的外圍，其格架堵塞比（格架處柵元流通面積/無格架處的柵元流通面積）較小，格架堵塞阻力系數(shù)較大，且無內部攪混翼，剛凸的個數(shù)少于典型柵元和冷壁柵元。因此燃料組件邊角柵元的阻力系數(shù)大于典型柵元和冷壁柵元的阻力系數(shù)。

因為邊角柵元阻力系數(shù)較大，且子通道分析中未做專門考慮，該問題在采用美國西屋電氣公司AP1000技術的三門和海陽核電廠的安全分析報告的審查中被安審當局要求做補充回答。西屋公司首先定性地回復認為當前子通道分析中假設的保守的徑向功率分布足以包絡邊、角柵元的較大的阻力系數(shù)影響，隨后采用子通道全堆芯模型進行了分析，證明邊角柵元的偏離泡核沸騰比計算結果能夠被典型和冷壁柵元包絡。西屋定性的回答沒有量化評估邊角柵元阻力系數(shù)的影響，西屋在采用全堆芯模型分析時選取的一種功率分布不能被證明是足夠保守和能夠包絡峰值功率出現(xiàn)在邊角柵元的情況，因此其回答和計算分析結果都不具有說服力。

失流事故是一個典型的偏離泡核沸騰比極限事故，對此，有過許多的相關研究［11 -15］。本文采用VIPRE-W程序，基于安全分析采用的標準子通道模型，調整子通道模擬區(qū)域的位置，修改通道、燃料棒以及格架阻力相關參數(shù)，對邊角柵元的綜合效應進行了直接模擬；同時針對安全分析采用的保守徑向功率分布進行了處理，確保邊、角柵元保守地處于徑向功率峰值附近，在此基礎上選取正常運行和失流事故進行了分析，以確定保守假設條件下邊角柵元偏離泡核沸騰比仍然能被典型柵元和冷壁柵元的計算結果所包絡。VIPRE-W程序是基于美國電力研究院（Etectric Power Research Institute，簡稱EPRI）開發(fā)的VIPRE -01程序的美國西屋公司改進版本，其穩(wěn)態(tài)模型可用于堆芯偏離泡核沸騰比及堆芯通道間橫向流的計算。

1 邊角柵元子通道模型簡介

如圖1所示，標準的子通道模型劃分是針對堆芯的1/8部分進行的。首先，選取中心組件的1/8位置作為熱組件，將其劃分為11個通道和13根棒；然后將其余燃料組件集總表示為通道和棒，依次增加所代表的組件數(shù)目，分別為12號通道（14號棒）、13號通道（15號棒）、14號通道（16號棒）。與圖1對應的通道及燃料棒份額數(shù)據(jù)詳見表1。

圖1 標準模型子通道劃分Fig. 1 Sub-channel Scheme of Standard Model

如圖2所示，邊角柵元子通道模型的劃分也是針對堆芯的1/8部分進行的。首先，選取相鄰組件（標準的子通道模型的12號通道和13號通道之間）的位置作為1/8熱組件，將其劃分為11個通道和16根棒（因為無導向管，該部分燃料棒增加3根）；然后，將其余的燃料組件集總表示為通道和棒，依次增加所代表的組件數(shù)目，分別為12號通道（17號棒）、13號通道（18號棒）、14號通道（19號棒）。

表1 標準子通道模型參數(shù)Table 1 Parameters of Standard Sub -channel Model

圖2 邊角柵元模型子通道劃分Fig. 2 Sub-channel Scheme of Side and Corner Cell Model

2 邊角柵元偏離泡核沸騰比分析建模

在邊角柵元的偏離泡核沸騰比分析中，首先需要針對子通道位置進行調整，重新計算相關通道的流通面積等參數(shù)；同時根據(jù)燃料組件邊角位置無導向管布置，調整對應的燃料棒編號和數(shù)目，以及各燃料棒在相應子通道中的份額；然后根據(jù)新的燃料棒數(shù)目，采用保守的徑向功率分布，重新確定各燃料棒的功率份額；最后根據(jù)邊角柵元的特殊性，對應設定其格架的阻力系數(shù)。

2. 1 通道參數(shù)

如圖2所示，邊角柵元子通道模型由于位置特殊，橫跨了堆芯中的4個組件。其中燃料棒3、5、6圍成的子通道3和燃料棒6、9、10圍成的子通道6，由于一邊是燃料棒中心距，另一邊是燃料組件間隙，二者長度不相等，使得通道3和6不是規(guī)則的等邊三角形形狀，但通道3和6在對稱上是相同的。

燃料組件之間的間隙寬度為燃料組件邊長減去16個燃料棒中心距。根據(jù)燃料組件間隙寬度和燃料棒中心距可得通道3的銳角角度。根據(jù)相關棒在其他位置相對應的份額，分別求得3號棒、6號棒、10號棒在各自對應通道中的份額。在計算各子通道參數(shù)的基礎上，將通道1至11（燃料棒1至16）按照組件中心線位置劃分為三角區(qū)域上下兩部分。邊角柵元子通道模型的通道12是在標準的子通道模型通道12的基礎上，增加1/8組件（標準通道1 -11號通道所在范圍），減去三角區(qū)域的上半部份。邊角柵元子通道模型的通道13是在標準的子通道模型通道13的基礎上，減去三角區(qū)域下半部份。邊角柵元子通道模型10個子通道中已不包含導向管份額。

2. 2 燃料棒參數(shù)

邊角柵元子通道模型總的燃料棒數(shù)目與標準的子通道模型相同，只是因為1/8熱組件的位置特殊布置，使得與子通道相連的燃料棒3、6、10在各自對應的通道中的份額比較特殊。同時由于1/8熱組件在通道12和13之間，使得通道12 （17號棒）和通道13（18號棒）的燃料棒份額計算比較復雜。

基于標準的子通道模型的徑向功率分布，如圖1所示，圍繞熱通道（通道5）的燃料棒徑向功率因子都較高。由于邊角柵元子通道模型選取的組件相鄰位置作為1/8熱組件，子通道1到10相連無導向管。因此，相比標準的子通道模型，3個導向管的位置新增加了3根燃料棒。為了與保守的徑向功率分布保持一致，將這三根燃料棒（3號、10號、12號）的徑向功率因子設置為與外圍通道中集總的19號棒一致。這樣既保持了原有的分布不變，使得熱通道（角柵元）周圍仍然有高功率份額的棒，計算結果仍然具有保守性，又使得總的徑向功率分布容易歸一。將燃料棒3號、10號、12號以及集總的19號份額求和，歸一處理后獲得3、10、12、19號棒的功率因子，詳見表2。

2. 3 格架阻力參數(shù)

邊角柵元格架阻力系數(shù)與典型和冷壁柵元的比較詳見表3。邊角柵元子通道模型選取的組件相鄰位置作為1/8熱組件，子通道3、4、6、8表示邊柵元，子通道5表示角柵元。在子通道程序輸入卡中增加4個阻力系數(shù)，然后在格架阻力部分中分別定義邊角柵元的阻力系數(shù)。

表3 邊角柵元阻力系數(shù)Table 3 Loss Coefficient of Side and Corner Sub-channels

3 邊角柵元偏離泡核沸騰比計算結果分析

選取正常運行和全失流事故工況，根據(jù)第2節(jié)計算的相關參數(shù)，考慮邊角柵元的綜合效應，進行計算分析。在正常運行和全部失流事故工況［2］下，典型柵元和冷壁柵元的最小偏離泡核沸騰比以及邊角柵元的偏離泡核沸騰比比較詳見表4。

計算結果表明：對典型柵元和冷壁柵元計算的最小偏離泡核沸騰比都能包絡邊角柵元。

表4 邊角柵元偏離泡核沸騰比分析計算結果Table 4 Results of side and corner sub-channel DNBR analysis

針對典型柵元、冷壁柵元和邊角柵元，分別抽取計算的子通道冷卻劑流量和空泡份額進行比較，結果表明，考慮邊角柵元的綜合效應建立的邊角柵元子通道模型計算所得，角柵元（子通道5）比標準子通道模型計算的典型柵元（子通道5）質量流量明顯增大，通道中同一軸向位置的空泡份額明顯減小。邊柵元（子通道3）和標準子通道模型計算的冷壁柵元（子通道3）相比也有同樣的結果。流量增大，含汽率減小，表明盡管邊角柵元阻力系數(shù)較大，且徑向功率分布仍然保守的情況下，其流通面積較大、相鄰間隙增加，且邊角柵元中無導向管，減少了冷壁效應，這些綜合效應使得邊角柵元的偏離泡核沸騰比不如典型柵元和冷壁柵元極限。

4 邊角柵元偏離泡核沸騰比分析展望

以上分析都是基于通用的子通道分析方法，采用燃料組件幾何參數(shù)的名義值進行計算。由于邊角柵元涉及燃料組件的間距問題，本文拓展考慮燃料組件間距的多種可能性進行進一步分析。另外選取了燃料組件安裝布置過程中可能出現(xiàn)的典型偏差進行分析。計算結果詳見表5。

在組件外圍棒中心距名義值（13.46 mm）的基礎上，分別考慮下管座名義間隙和最小間隙之差（0.31mm），以及中間格架間隙（1.22mm）進行分析。不利的情況是：分別假設下管座為最小間隙，使得組件外圍棒中心距減小0.31mm，變?yōu)?3.16mm；以及假設中間格架碰在一起，間隙為零，此時組件外圍棒中心距變?yōu)?2.24mm。有利的情況是：假設下管座取最大間隙，使得組件外圍棒中心距名義值增加0.31mm，變?yōu)?3.77mm。

分析結果表明，即使考慮到下管座間隙最小，計算所得邊角柵元的偏離泡核沸騰比小于外圍棒中心距為名義值的結果，但仍然大于標準模型計算的典型柵元和冷壁柵元；若考慮下管座間隙最大，則所得邊角柵元偏離泡核沸騰比裕量有進一步增加。只有考慮中間格架間隙為零時，邊角柵元偏離泡核沸騰比會小于標準模型計算的典型柵元和冷壁柵元，因為此時燃料組件外圍棒中心距已經小于燃料組件中燃料棒的中心距（12. 60mm），這表明此時邊角柵元的流通面積已小于典型柵元和冷壁柵元，在阻力系數(shù)較大，徑向功率分布仍然保守的情況下，計算所得偏離泡核沸騰比肯定更小。

這一分析的目的只是為了進一步研究邊角柵元綜合效應中組件外圍棒中心距（或燃料組件之間的間隙）影響。盡管計算中假設的中間格架間隙為零這一計算結果表明在這種情況下邊角柵元的偏離泡核沸騰比小于標準模型計算的典型柵元和冷壁柵元。但這只是計算中一種特定的極限假設，即使真實的燃料組件中出現(xiàn)了中間格架相互碰撞，間隙為零的情況，也不會同時出現(xiàn)整個燃料組件自下至上都緊貼在一起的情況。而且如果組件間隙過小，水隙減少，會引起中子慢化減弱，相應的棒功率可能也會降低。本分析中同時假設邊角柵元周圍燃料棒具有最高的徑向功率份額，這種組合假設是極不可能出現(xiàn)的情況。

本文計算是基于子通道分析采用的典型徑向功率分布，這種分布是工程經驗結合計算分析獲得的，被證明是保守的。至于組件間隙變化對這一分布的影響，本文未做考慮。這一影響可能需要在未來實現(xiàn)全堆芯三維核熱耦合計算，針對每根燃料棒和每個柵元進行精細模擬，再詳細考慮子通道建模時組件間隙的不利偏差，同時在物理分析程序中也考慮組件間隙的不利偏差時予以解決。

表5 不同相鄰組件外圍棒中心距計算偏離泡核沸騰比比較Table 5 Compare of DNBR with different pitch distance

參考文獻

［1］楊萍，賈紅軼，王喆. RTDP方法在大型先進壓水堆熱工設計中的應用初步研究［J］.原子能科學技術，2012，47 （7）：1182 -1186.

［2］周洲，傅先剛. RTDP方法用于大亞灣核電站熱工水力設計的初步研究［C］.全國反應堆熱工流體會議，1999.

［3］張志達.統(tǒng)計DNBR的先進方法及應用研究［C］.全國反應堆熱工流體會議，1999.

［4］董博，匡波，朱學農.非參數(shù)統(tǒng)計方法在DNBR限值計算中的應用［J］.核技術，2013，36（3）：1 -6.

［5］王煦嘉.應用抽樣統(tǒng)計方法計算DNBR限值［J］.原子能科學技術，2012，46（4）：448 -452.

［6］李經緯.秦山核電二期工程反應堆熱工水力設計［J］.核動力工程，1999，20（4）：308 -321.

［7］陳軍，周有新，李石磊，等.跨間交混格架對EPR堆芯DNBR裕量的影響［J］.原子能科學技術，2013，47（2）：249 -253.

［8］陳毓湘，周全福.秦山核電廠反應堆熱工水力設計［J］.核科學與工程，1985，5（4）：45 -51.

［9］丁振鑫. 5MW低功率堆穩(wěn)態(tài)熱工水力計算及分析［J］.核動力工程，1992，13（4）：12 -17.

［10］王榮忠. 10MW研究堆堆芯熱工水力設計［J］.核動力工程，1998，19（4）：293 -296.

［11］郝老迷，謝輝.堆芯冷卻劑流量下降DNBR計算［J］.原子能科學技術，1993，27（5）：422 -425.

［12］郝老迷.秦山核電廠主泵軸卡死事故堆芯DNBR計算［J］，原子能科學技術，1993，27（4）：309 -314.

［13］許克，俞爾俊.秦山核電廠失流事故的研究［C］.北京：中國原子能科學研究院年報，1987.

［14］李運文，郝老迷.秦山核電廠失流事故分析計算［C］.第五屆全國反應堆熱工流體學術交流會論文集，1995.

［15］靖劍平，張春明，孫微，等. AP1000冷卻劑強迫流動全部喪失事故瞬態(tài)特性分析［J］.核科學與工程，2014，34（4）：541 -546.

DNBR Analysis of Fuel Assembly Side and Corner Channels

XU Zhihong，WANG Zhe，YANG Ping，SHI Guobao
（Shanghai Nuclear Engineering Research & Design Institute，Shanghai 200233，China）

Abstract：The loss coefficient of side and corner channels of fuel assembly is greater than that of the typical and thimble channels. Based on the standard sub-channel model，DNBR analysis of the side and corner channels is conducted by adjusting the sub-channel modeling position，and updating the parameters of related channels，rods and grid loss coefficient，as well as directly simulating the side and corner channels. The result indicates that DNBR of the side and corner channels can be enveloped by that of the typical and thimble channels. This paper further analyzes the impact of fuel assembly pitch toleranceon the side and corner channel DNBR.

Key words：fuel assembly；side and corner channel；DNBR

中圖分類號：TL429

文章標志碼：A

文章編號：1672-5360（2016）02-0084-05

收稿日期：2016-02-05 修回日期：2016-03-21

基金項目：上?？茖W技術委員會資助課題，項目編號 13DZ2250200

作者簡介：許志紅（1983 -），男，湖南岳陽人，博士，現(xiàn)主要從事反應堆熱工水力設計及安全分析工作