钠冷快堆堆芯捕集器设计优化数值研究.pdf

第 卷第期原子能科学技术V o l ,N o 年月A t o m i cE n e r g yS c i e n c ea n dT e c h n o l o g yA p r 钠冷快堆堆芯捕集器设计优化数值研究曹胜,张斌,王文鹏,单建强,(西安交通大学 能源与动力工程学院,陕西 西安 ;西安交通大学 动力工程多相流国家重点实验室,陕西 西安 )摘要:堆芯捕集器是为有效分散堆芯熔融物并防止压力容器下封头与熔体大规模接触而被破坏的非能动预防和缓解装置.在钠冷快堆(S F R)中,堆芯捕集器的结构直接影响碎片床的堆积形状和分布,进而影响碎片床的再临界性和长期衰变热去除能力.本文针对堆芯捕集器的结构设计优化开展数值研究,重点关注其烟囱结构设计对碎片床形成和分布的影响机理及规律.基于无量纲刚度系数和无量纲阻尼系数改进离散元法(D EM),通过改变堆芯捕集器烟囱顶盖垂直投影边长、顶盖倾斜角度和烟囱间距,研究碎片颗粒的运动和碎片床的形成行为.结果表明,堆芯捕集器的烟囱顶盖垂直投影边长、烟囱顶盖倾角和烟囱间距对碎片床的堆积形状和分布均有重要影响,碎片颗粒的二次散射对于改善碎片床的均匀性至关重要.关键词:钠冷快堆;堆芯捕集器;离散元法;碎片床中图分类号:T L 文献标志码:A文章编号:()收稿日期:;修回日期:通信作者:张斌d o i:/y z k y o u x i a n N u m e r i c a l S t u d yo nD e s i g nO p t i m i z a t i o no fS o d i u m c o o l e dF a s tR e a c t o rC o r eC a t c h e rC AOS h e n g,Z HANGB i n,WANG W e n p e n g,S HANJ i a n q i a n g,(S c h o o l o fE n e r g ya n dP o w e rE n g i n e e r i n g,X ia nJ i a o t o n gU n i v e r s i t y,X ia n ,C h i n a;S t a t eK e yL a b o r a t o r yo fM u l t i p h a s eF l o wi nP o w e rE n g i n e e r i n g,X ia nJ i a o t o n gU n i v e r s i t y,X ia n ,C h i n a)A b s t r a c t:I nt h ee v e n t o f c o r ed i s a s s e m b l ya c c i d e n t(C D A)i nas o d i u m c o o l e d f a s t r e a c t o r(S F R),t h em o l t e n m a t e r i a lf r o mt h ec o r em i g r a t e sa n di n t e r a c t sw i t ht h el o w t e m p e r a t u r ec o o l a n t i nt h e l o w e rp l e n u m,r e s u l t i n gi nt h ef o r m a t i o no fc o r ed e b r i st h a te v e n t u a l l ys e t t l e st of o r m sad e b r i sb e d I m p r o p e rp l a c e m e n ta n di n a d e q u a t ec o o l i n go ft h ed e b r i sb e dc a nr e s u l t i nt h e f a i l u r eo f t h ep r e s s u r ev e s s e l T oe f f i c i e n t l yd i s p e r s e t h ec o r ed e b r i sa n dm i n i m i z ec o n t a c tb e t w e e nt h e l o w e rh e a do f t h ep r e s s u r ev e s s e l a n dt h em o l t e nm a t e r i a l,a c o r e c a t c h e rw a s s p e c i f i c a l l yd e s i g n e da s ap a s s i v ep r e v e n t i v e a n dm i t i g a t i n gd e v i c e T h e s t r u c t u r eo f t h e c o r e c a t c h e rd i r e c t l y i m p a c t s t h e s h a p e a n dd i s t r i b u t i o no f t h ed e b r i sb e d,w h i c hs u b s e q u e n t l ya f f e c t s i t sr e c r i t i c a l i t ya n dl o n g t e r md e c a yh e a tr e m o v a lc a p a b i l i t y T h i ss t u d ya i m st o p e r f o r m a n u m e r i c a li n v e s t i g a t i o nt oo p t i m i z e t h es t r u c t u r a ld e s i g no ft h ec o r ec a t c h e r,w i t has p e c i f i cf o c u so nt h em e c h a n i s m sa n d l a w so f t h ec h i m n e ys t r u c t u r eo f t h ec o r ec a t c h e ro nt h e f o r m a t i o na n dd i s t r i b u t i o no f t h ed e b r i sb e d T h e i m p r o v e dd i s c r e t e e l e m e n tm e t h o d(D EM),w h i c h i sb a s e do nd i m e n s i o n l e s ss t i f f n e s sa n dd i m e n s i o n l e s sd a m p i n gc o e f f i c i e n t s,w a se m p l o y e dt on u m e r i c a l l ys i m u l a t e t h e f o r m a t i o np r o c e s so f t h ed e b r i sb e do nt h ec o r ec a t c h e r T h i sm e t h o di n c o r p o r a t e ss o u r c et e r m si n t ot h em o m e n t u me q u a t i o no fd e b r i sp a r t i c l e st oa c c o u n tf o rt h ei n f l u e n c eo fc o o l a n to nt h e i r m o t i o n V a l i d a t i o ne x p e r i m e n t s w e r ec o n d u c t e dt ov e r i f yt h ec a p a b i l i t yo f t h en u m e r i c a l s i m u l a t i o na l g o r i t h mu s e di no p t i m i z i n gt h ed e s i g no f t h ec o r ec a t c h e r I nt h i sp a p e r,t h ee f f e c t so f t h r e e f a c t o r so np a r t i c l em o t i o na n dd e b r i sb e df o r m a t i o nm e c h a n i s m sw e r ei n v e s t i g a t e db yv a r y i n gt h ev e r t i c a lp r o j e c t i o ns i d e l e n g t ho f t h ec h i m n e yc o v e r,i n c l i n a t i o na n g l eo f t h ec o v e r,a n ds p a c i n gb e t w e e nc h i m n e y s T h er e s e a r c hf i n d i n g ss u g g e s tt h a tt h es e l e c t i o no fa p p r o p r i a t ev a l u e s f o r t h e s ep a r a m e t e r sh a s as u b s t a n t i a l i m p a c t o nt h e s h a p ea n dd i s t r i b u t i o no f t h ed e b r i sb e d T h e r e f o r e,i t i s r e c o mm e n d e d t oe s t a b l i s hr e a s o n a b l e c h i m n e yd e s i g np a r a m e t e r s t oa c h i e v eo p t i m a lp e r f o r m a n c e i np r a c t i c a la p p l i c a t i o n s F u r t h e r m o r e,t h i ss t u d yr e v e a l s t h a td u r i n gt h ef o r m a t i o no f t h ed e b r i sb e d,t h ed e b r i sp a r t i c l e sd e m o n s t r a t eas e c o n d a r ys c a t t e r i n ge f f e c t,w h i c hp l a y sav i t a l r o l e i ne n h a n c i n gt h eu n i f o r m i t yo f t h ed e b r i sb e d T h i ss t u d yw i l lp r o v i d ev a l u a b l ee n g i n e e r i n gr e f e r e n c e sf o ro p t i m i z i n gt h ed e s i g no f t h ec o r ec a t c h e r K e yw o r d s:s o d i u m c o o l e df a s t r e a c t o r;c o r ec a t c h e r;d i s c r e t ee l e m e n tm e t h o d(D EM);d e b r i sb e d钠冷快堆(S F R)堆芯解体事故(C D A)可能对堆内构件和压力容器造成严重损害,是放射性包容与纵深防御的重要安全隐患.当熔融物下降并与下腔室中较低温度的冷却剂接触后,会骤冷形成放射性碎片,最终沉降在反应堆下腔室形成碎片床.若碎片床没有被妥善地滞留和冷却,可能导致压力容器失效.因此,钠冷快堆中设计了堆芯捕集器作为预防与缓解装置,用于收集和容纳熔融物碎片,以尽可能地实现对碎片床的长期冷却,确保主容器的完整性.目前已有不少学者研究了碎片床的冷却特性和堆芯捕集器的性能.S u d h a等通过数值传热分析估算了钠冷快堆在无保护失流事故(U L O F A)和热阱丧失事故(P L OH S)下熔融物转移到堆芯捕集器的最短和最长时间,这些估算时间对于堆芯捕集器的性能设计具有重要的参考意义.D a v i d等分析了S F R池内自然对流的瞬态发展,并评估了不同堆芯熔化场景下堆芯捕集器底板的最高温度.R o y c h o w d h u r y等使用一维模型证明了印度原型快堆(P F B R)能够安全消除 个组件融化带来的放射性衰变热.V l a s i c h e v等使用S T O B E O代码对B N 捕集器上的碎片床开展了数值分析,并评估出碎片床内的最高温度低于燃料沸腾温度.在许多研究中碎片床的形状是任意选择的,如圆柱形、圆锥形、高斯形、堆状、圆锥状或同向堆积.然而,碎片床的形状和分布很重要,它是防止碎片床重返临界和衰变热长期冷却的关键影响因素.一些快堆在设计和建造阶段就依据反应堆特征对堆芯捕集器进行了优化.例如,英国D o u n r e a y快堆的下腔室中设计了燃料分散锥和防溅板,以避免堆芯熔融物结块,并将熔融物引导至 堆 芯 捕 集 器 ;英 国 商 业 示 范 快 堆(C D F R)采用多个垂直管道焊接成圆锥形管板收集堆芯碎片,并连接到一个中央烟囱以促进冷却剂的自然对流 ;俄罗斯B N 和B N 所配备的堆芯捕集器,采用个垂直引流管以增强自然对流;欧洲快堆(E P R)配备了箱式结构的堆芯捕集器,并采用双烟囱结构以更好地分散堆芯碎片,.此外,还有采用多托盘的堆芯捕集器,将堆芯碎片分开进行冷却 ,采原子能科学技术第 卷用多层托盘的堆芯捕集器设计,消除了冷却剂对堆芯捕集器部件长期侵蚀和腐蚀的隐患,.这些设计通过工程实践对堆芯捕集器的结构及组件进行了优化.烟囱作为堆芯捕集器的主要结构之一,不仅影响冷却剂的自然对流,而且对碎片床的堆积和分布具有重要影响,因此亟需进行详细的研究与设计优化.针对目前堆芯捕集器烟囱结构设计优化研究不足的现状,本文采用改进离散元法(D EM)作为数值模拟工具,以中国实验快堆(C E F R)为原型,对碎片床在堆芯捕集器上的形成过程进行数值模拟.通过改变堆芯捕集器托盘上的烟囱顶盖垂直投影边长、烟囱顶盖倾角以及烟囱间距,研究堆芯捕集器烟囱结构对碎片床堆积形状和分布的作用机理和影响规律,为碎片床的传热性能分析和长期可冷却性评估提供前提条件,并对钠冷快堆堆芯捕集器的设计及优化提供工程参考.数值模型本文采用D EM模拟堆芯捕集器上碎片床的形成过程.在D EM计算中,将熔融物碎片假设为粒子,其运动遵循牛顿定律.碎片颗粒的运动分为平移和旋转,并通过以下控制方程分别进行求解:miduidtFi()IididtTi()式中:mi和Ii为颗粒i的质量和转动惯量;ui和i为颗粒i的平移速度和旋转速度.两粒子碰撞时粒子i的受力如图所示,其中Fi和Ti分别为作用在粒子i上的力和力矩,由方程()和()求解.Fij(Fni jFsi j)Fb u o y,iFd r a g,iFg,i()Tij(RiFsi jTri j)()式中:Fni j和Fsi j为粒子i和j的法向接触力和切向接触力;Fb u o y,i为粒子i受到的浮力;Fd r a g,i为粒子i受到钠冷却剂的曳力;Fg,i为粒子i的重力;Ri为粒子i的等效半径;Tri j为颗粒i和j由于弹性滞后损失和黏性耗散导致的滚动摩擦力矩.图粒子i与粒子j碰撞的受力示意图F i g S c h e m a t i cd i a g r a mo f f o r c eo nc o l l i s i o nb e t w e e np a r t i c l eia n dp a r t i c l ej粒子间和粒子与壁面间碰撞产生的法向接触力Fn和切向接触力Fs采用弹簧阻尼模型计算,如式()、()所示,图为弹簧阻尼模型示意图.Fnknncndndt()Fs m i n(ksscsdsdt,Fn)s i g n(ksscsdsdt)()式中:为粒子接触点的相对位移;k为刚度系数;c为阻尼系数;为摩擦系数.图线性弹簧阻尼模型的粒子碰撞示意图F i g P a r t i c l ec o l l i s i o n i l l a s t r a t i o nu s i n g l i n e a rs p r i n g d a m p i n gm o d e l法向刚度系数kn和阻尼系数cn分别由式()和式()计算.kn m(t nt)()cnc nm kn()式中:t n为无量纲刚度修正系数;c n为无量纲第期曹胜等:钠冷快堆堆芯捕集器设计优化数值研究阻尼修正系数;t为数值模拟中的时间步长.数值模拟中,首先确定无量纲刚度系数、无量纲阻尼系数 和时间步 长的数值.在D EM中,时间步长的选取需要考虑计算效率和计算精度,并结合计算机的运行性能满足公式t(m/kn)/.通过单个颗粒自由落体后的回弹高度确定无量纲刚度系数t n的值.通过计算组颗粒沉降算例,限制粒子堆积的重叠度小于颗粒直径的,确定无量纲阻尼系数c n.使用t n、c n计算kn和cn后,通过式()、()计算切向刚度系数ks和阻尼系数cs.kskn()()cscn()()式中,为泊松比.在均匀等向性材料中,剪切模量G、杨氏模量E与满足GE/().本文采用t n和c n替代传统D EM计算中的kn、ks、cn和cs,不仅减少了需要确定的参数数量,也减少了确定参数取值带来的工作量.除颗粒间的碰撞外,颗粒也受到冷却剂的曳力Fd r a g、浮力Fb u o y以及颗粒自身的重力Fg,分别通过式()、()和()计算.Fd r a gCd r a g(d)du()Fb u o y(d)cg()Fg(d)cg()式中:d为颗粒的等效直径;d为颗粒的密度;c为冷却剂的密度;u为颗粒与冷却剂的相对速度;g为重力加速度;Cd r a g为曳力系数,是绕流雷诺数R e的函数.本文采用编程方式进行D EM求解,计算流程如图所示.在D EM计算中,采用显式方案更新粒子的速度和位置,计算数据通过开源软件T e c p l o t进行后处理,以可视化颗粒的运动和碎片床的形成过程.数值建模与验证 数值建模本文以C E F R堆芯捕集器为 研究对象.在C E F R堆芯捕集器设计中,考虑了两种事故假设:)假设转化区和活性区的所有燃料组件图D EM计算流程图F i g D EMc a l c u l a t i o nf l o wc h a r t熔化并迁移至堆芯捕集器;)假设/的燃料组件迁移到堆芯捕集器.事故参数如下:活性区燃料质量为 t,转化区燃料质量为 t,包壳与盒壁的质量为 t.假设a 整体结构;b 底部托盘俯视图图堆芯捕集器的建模图纸F i g M o d e l i n gs k e t c ho f c o r ec a t c h e r事故后混合的熔融物体积不变,约为 m.为确保安全裕度,假设事故形成的碎片床孔隙率为 ,堆芯捕集器的托盘设计容积为 m.综合考虑数值模拟的计算效率和验证实验的对照等因素后,本研究采用比例为的模型对堆芯捕集器进行建模和计算.图为堆芯捕集器的建模图纸,该装置由带有个烟囱的底部托盘和正上方的“漏斗”状释放装置组成,烟囱设置为个中央烟囱和个周边烟囱,且周边烟囱呈等距排布.考虑到堆芯捕集器具有旋转对称性,仅对堆芯捕集器的/进行建模和计原子能科学技术第 卷算,建模示意如图所示.堆芯捕集器结构由黑色粒子组成,熔融物碎片颗粒由绿色粒子构成,粒子直径为 mm.熔融物粒子静止于释放装置内,在模拟开始后依靠重力自由下落.数值模拟中熔融物粒子和堆芯捕集器壁面粒子的参数如下:密度均为 k g/m,泊松比均为 ,t n均为 ,c n均为 ,静摩擦系数均为 .图堆芯捕集器建模F i g S t r u c t u r a lm o d e l i n go f r e a c t o rc o r ec a t c h e r在钠冷快堆严重事故中,冷却剂通过阻力和浮力影响碎片的下落行为,本研究仅考虑冷却 剂 对 堆 芯 碎 片 的 单 向 作 用.通 过 在碎片颗粒 的 动 量 方 程 中 加 入 冷 却 剂 的 曳 力和浮力,考虑冷却剂对堆芯碎片的影响.曳力计算详见式(),其中曳力系数计算公式如式()所示.Cd r a gcd r a gdc(uduc)()式中:ud和uc分别为碎片颗粒和钠冷却剂的速度;cd r a g为阻力系数,与绕流雷诺数R e有关,通过式()计算.在本研究中,钠冷却剂视为大空间内流体,冷却剂密度为c,碎片颗粒的等效直径为d.cd r a gkR en()R ec(uduc)dc()由于堆芯碎片的密度远大于冷却剂的密度,因此认为碎片的运动处于平方阻力区,取k 、n.钠冷却剂的物性参数如下:液位为 m、密度为 k g/m、动力黏度为 Ns/m.数值验证本研究通过开展数值模拟验证实验,以证明所使用数值方法具备堆芯捕集器设计与优化的研究能力.图为验证实验装置示意图,该装置由装置主体、辅助支撑台架、数据采集系统和进排水系统组成.装置主体由透明有机玻璃制成的底部托盘和不锈钢漏斗释放装置组成,如图所示.底部烟囱由透明有机玻璃制成,为方便更换,与底部托盘采用分离式设计,通过磁铁和固定销固定.数据采集系统由两套装置组成,一套是高速摄像机,用于连续拍摄碎片的下落和迁移过程;另一套是三维激光扫描仪,在颗粒堆积稳定后对碎片床进行扫描,获取相关的堆积形状数据.图堆芯捕集器实验平台示意图F i g S c h e m a t i cd i a g r a mo f r e a c t o rc o r ec a t c h e re x p e r i m e n t a lp l a t f o r m图堆芯捕集器的装置主体F i g M a i nb o d yo fd e v i c e f o rc o r ec a t c h e rv a l i d a t i o ne x p e r i m e n t实验使用直径mm的氧化铝球形颗粒代替碎片颗粒,总体积为L.堆芯捕集器底部托盘内的中央烟囱与周围烟囱的距离设定为 mm,烟囱顶盖的垂直投影边长为 mm,顶盖倾角为 ,如图所示.实验开始前,将颗粒倒入释放装置中并抚平表面,向底部托盘内注入指定高度的水.待水面静止后,迅速释放漏斗内的颗粒.待颗粒堆积稳定后,通过排第期曹胜等:钠冷快堆堆芯捕集器设计优化数值研究水口缓慢排出托盘内的水.使用三维激光扫描仪对托盘内形成的碎片床进行扫描,得到碎片床堆积形状和颗粒分布数据.实验过程中高速摄像机拍摄到的照片如图所示,三维激光扫描仪扫描得到的颗粒分布如图 所示,其中缺失部分是烟囱占据的空间.图验证实验中烟囱的结构与参数F i g S t r u c t u r ea n dp a r a m e t e ro f c h i m n e yu s e df o rv a l i d a t i o ne x p e r i m e n t为定量比较数值模拟结果与实验结果之间的差异,选取个指标评价碎片床颗粒分布的均匀性,分别是碎片床的高度方差Zv a r、高度归一化方差Zn o r v a r和三维模化方差,计算公式如式()()所示,它们可以反映碎片床在高度方向和空间维度上的不均匀程度.Zv a rNNizi()Zn o r v a rNNi(zizz)()NNi(xixD)(yiyD)(zizH)()式中:N为碎片床的粒子数量;zi为粒子所处高度;z为碎片床内粒子的平均高度.图高速摄像机拍摄的不同时刻的粒子运动F i g P a r t i c l em o t i o na td i f f e r e n t t i m ec a p t u r e db yh i g h s p e e dc a m e r a图 碎片床稳定堆积后的分布F i g D i s t r i b u t i o nc l o u da t l a so fd e b r i sb e da f t e r s t a b l ea c c u m u l a t i o n将数值模拟获得的碎片床表面粒子坐标和三维 激 光 扫 描 仪 扫 描 的 粒 子 坐 标 按 照 式()、()计算,结果列于表.从表可看到,数值模拟结果与实验结果之间的相对误差最大值为 ,在可接受范围内,表明本研究所采用的数值模拟方法能够有效模拟堆芯捕集器上碎片床的形成过程.数值结果分析在数值模拟研究中,通过改变烟囱顶盖垂直投影边长、顶盖倾角和烟囱间距探究堆芯捕集器的烟囱设计对碎片床形状和分布的作用机理及影响规律,从而为堆芯捕集器的结构设计和性能优化提供建议和指导.表数值模拟和实验中颗粒床粒子群的高度方差和三维模化方差T a b l eH e i g h t v a r i a n c ea n dt h r e e d i m e n s i o n a ln o r m a l i z e dv a r i a n c eo fd e b r i sb e di nn u m e r i c a l s i m u l a t i o na n de x p e r i m e n t参数模拟值实验值相对误差/Zv a r,m 烟囱顶盖垂直投影边长在数值模拟中,保持中央烟囱与周围烟囱间距为 mm、烟囱顶盖倾角为 ,通过改变烟囱顶盖的垂直投影边长(mm)形成第组模拟工况(工况 a g),如图 所原子能科学技术第 卷示.在模拟完成后,得到碎片床粒子的高度方差和高度归一化方差,结果如图 所示.从图 可观察到,随着烟囱顶盖垂直投影边长的增加,碎片床的高度方差和高度归一化方差呈相同的变化趋势.当烟囱顶盖垂直投影边长为 mm时,碎片床的高度方差和高度归一化方差达到最 小,分别为 mm和 .这表明在该组工况中,烟囱顶盖垂直投影边长为 mm时,堆芯捕集器表现出最佳的碎片颗粒分散性能.图 烟囱结构参数示意图F i g S c h e m a t i c i l l u s t r a t i o no f c h i m n e yss t r u c t u r a lp a r a m e t e r对比发现,碎片颗粒在不同的烟囱顶盖上运动及后续行为具有差异.在工况 a中,没有观察到颗粒的二次散射,颗粒从烟囱顶盖滑落后于烟囱底部附近聚集.而在工况 g中,由于烟囱顶盖过大导致烟囱间缝隙较小,使得颗粒下落时在顶盖边缘处滞留并在烟囱缝隙下方发生明显堆积.在工况 c中,烟囱底部附近没有出现碎片床的凹陷,烟囱缝隙下的颗粒堆积也不明显,烟囱顶盖使部分颗粒在顶盖上发生二次散射,弹射到更远的位置,从而形成最均匀的碎片床分布.图 为工况 a g中碎片床在y z平面上的投影.由图 可观察到,随着烟囱顶盖垂直投影边长的增加,滞留在堆芯捕集 器 烟 囱 顶 盖 上 的 碎 片 颗 粒 数 量 也 增加.碎片颗粒滞留的原因是最后下落的颗粒不再受后续颗粒的撞击,导致摩擦力大于沿烟囱切线方向的重力分量.所有工况中碎片床的堆积角均约为 ,但碎片颗粒分布直径D(D(xy)/)/)不同,工况 d中的碎片床分布直径最大,为 mm,如图 所示.在设计和优化堆芯捕集器时,应确定合理的烟囱顶盖垂直投影边长,使堆芯捕集器具备更好的性能.烟囱顶盖倾角保持中央烟囱与周围烟囱间距为 mm,烟囱顶盖垂直投影边长为 m i n,设计烟囱顶盖倾角为 ,形成第组工况(工况 a g),如图 所示.通过数值模拟,得到碎片床的高度方差和高度归一化方差随堆芯捕集器烟囱顶 盖 倾 角 的 变 化 结 果,如 图 所 示.由图 可看出,随着烟囱顶盖倾角的增大,碎片床的高度方差和高度归一化方差呈先减小后增大的趋势.当烟囱顶盖倾角为 时,堆芯捕集器托盘上形成的碎片床的高度方差和高度归一化方差达到最小值,分别为 mm和 .这表明,烟囱顶盖倾角为 时,图 碎片床的高度方差和高度归一化方差F i g H e i g h tv a r i a n c ea n dh e i g h tn o r m a l i z e dv a r i a n c eo fd e b r i sb e d第期曹胜等:钠冷快堆堆芯捕集器设计优化数值研究图 不同烟囱盖垂直投影边长下碎片床的y z平面投影F i g y zp l a n ep r o j e c t i o no fd e b r i sb e du n d e rd i f f e r e n t c h i m n e yc a ps i d e l e n g t h s图 不同烟囱顶盖垂直投影边长下碎片床的最终分布直径F i g F i n a l d i s t r i b u t i o nd i a m e t e ro fd e b r i sb e du n d e rd i f f e r e n t c h i m n e yc a ps i d e l e n g t h s图 工况 a g的堆芯捕集器烟囱顶盖倾斜角度示意图F i g S c h e m a t i c i l l u s t r a t i o no f c h i m n e y sp a r a m e t e r s f o rc a s e s at o g堆芯捕集器上的碎片床分布最均匀.需要指出的是,工况 b碎片床的统计方差小于工况 c,这是由于在工况 b中少量颗粒滞留在释放装置中,因此在统计方差时这些粒子被排除在外.图 为下落颗粒在垂直于xy平面并与x轴成 夹角的投影.从图 可看出,除了工况 f和 g外,其他工况中碎片颗粒在烟囱顶盖上的滞留厚度均匀,这表明在下落过程中原子能科学技术第 卷图 碎片床高度方差和高度归一化方差随堆芯捕集器烟囱顶盖倾角的变化F i g V a r i a t i o n i nh e i g h tv a r i a n c ea n dh e i g h tn o r m a l i z e dv a r i a n c eo fd e b r i sb e dw i t hc h i m n e yc a p i n c l i n a t i o na n g l e图 模拟时间 s时不同烟囱顶盖倾角下颗粒的运动F i g M o t i o no fp a r t i c l e sa td i f f e r e n t c h i m n e yc a p i n c l i n a t i o na n g l e sa t s烟囱顶盖上的碎片颗粒数量处于动态平衡状态.在工况 a中,碎片颗粒几乎沿着烟囱顶盖向下滑落,没有观察到二次散射现象,且堆芯捕集器托盘上的碎片床分布直径也最小.在工况 e和 f中,可以观察到颗粒在堆芯捕集器烟囱顶盖上发生明显的二次散射.而在工况 f和 g中,烟囱顶盖上出现了明显的颗粒滞留现象,且二次散射的碎片颗粒数量也是该组模拟工况中最多的.堆芯捕集器的烟囱顶盖对碎片颗粒的运动起到了减缓或反弹的作用,导致颗粒可能出现两种不同的后续行为.第一种运动行为是碎片颗粒被滞留在烟囱顶盖表面,并沿顶盖滑落,最终堆积在烟囱底部附近.另一种运动行为是碎片颗粒的二次散射,颗粒在烟囱顶盖表面发生反弹,继而落在更远的位置.在碎片床的形成过程中,一些颗粒滞留在烟囱顶盖表面形成堆积层,而另一些颗粒则被堆积层反弹形成散射层.当烟囱顶盖的倾斜角度增加到一定程度时,落到烟囱顶盖上的颗粒与从烟囱顶盖上二次弹射的颗粒发生大量碰撞,导致颗粒的散射效应逐渐被抵消.需要指出的是,堆芯捕集器的最佳烟囱顶盖倾角与碎片颗粒的物性和机械性能相关.即在设计和优化具有相似结构的堆芯捕集器时,烟囱顶盖的倾斜角度应根据碎片颗粒的物性和机械性能来确定.通过优化堆芯捕集器烟囱顶盖倾角,可以使碎片床的形状和分布更加均匀,从而提升碎片床的热工水力性能.中央烟囱与周围烟囱的间距保持烟囱顶盖倾角为 ,烟囱顶盖垂直投影边长为 mm,设计中央烟囱与周围烟囱的距离为 mm,形成第组工况(工况 a l).图 为该组工况碎片床颗粒的高度第期曹胜等:钠冷快堆堆芯捕集器设计优化数值研究图 碎片床颗粒的高度方差和高度归一化方差随烟囱间距的变化F i g V a r i a t i o n i nh e i g h tv a r i a n c ea n dh e i g h tn o r m a l i z e dv a r i a n c eo fd e b r i sp a r t i c l e sw i t hd i s t a n c eb e t w e e nc e n t r a l c h i m n e ya n ds u r r o u n d i n gc h i m n e y s方差和高度归一化方差随堆芯捕集器烟囱间距的变化.结果显示,随着烟囱间距的增加,碎片床的高度方差和高度归一化方差呈整体下降趋势.当中央烟囱与周围烟囱间距为 mm时,碎片床的方差达到最小值,表明工况 l中碎片床分布最为均匀.通过观察碎片颗粒的运动轨迹发现,中央烟囱与周围烟囱的间距对碎片颗粒的运动具有复杂的影响.当周围烟囱的位置恰好位于粒子二次散射的落点时,它们对碎片床的均匀分布起到积极作用.在该组工况中,碎片床的统计方差随着烟囱间距的增加呈整体下降趋势,表明所选取的烟囱间距并未符合周围烟囱位于颗粒二次落点的情况.然而,碎片床的统计方差在烟囱间距为 mm的范围内出现了反复波动,由此推测粒子二次落点可能处于该范围内.但由于受计算工作量的限制,未对此猜想做进一步验证.在设计和优化钠冷快堆的堆芯捕集器时,需要谨慎考虑周围烟囱的数量和位置.若周围烟囱处于碎片颗粒二次散射的落点,可改善碎片床的均匀分布;反之,周围烟囱的存在会阻碍碎片颗粒的分散.烟囱间距越大,周围烟囱的阻碍作用越小,这可通过图 中碎片床的统计方差得以证实.结论与展望本工作采用改进的无量纲系数离散元法,对钠冷快堆碎片床形成过程和堆芯捕集器的优化设计进行了数值研究.主要探究了堆芯捕集器的烟囱结构对碎片床的堆积形状和碎片分布的影响机制.模拟结果表明,碎片颗粒在堆芯捕集器烟囱顶盖上的二次散射明显改善了碎片床的形状和分布.在设计和优化类似结构的堆芯捕集器时,可以选择广泛的结构参数进行初步实验或模拟.根据本工作探究的碎片床形状和分布规律,缩小参数取值范围,并最终确定最佳设计参数.需要注意的是,堆芯捕集器的烟囱顶盖形状和烟囱高度也可能对碎片床的形状和分布产生影响,这方面的研究将在后续工作中展开.参考文献:D AV I DDK,R AOPM,NA S H I N EBK,e t a l N u m e r i c a l i n v e s t i g a t i o no f t u r b u l e n tn a t u r a l c o n v e c t i o n i nt h e l o w e rp l e n u mo f s o d i u mc o o l e d f a s tr e a c t o ru n d e rc o r er e l o c a t i o ns c e n a r i oJN u c l e a rE n g i n e e r i n ga n dD e s i g n,:S U DHA