球团流与球流在球床内的混合卸料特性研究

原子能科学技术

AtomicEnercienceandTechnologySgy

Vol．５４,No．２Feb．２０２０

球团流与球流在球床内的混合卸料特性研究

∗２

,,葛　良１,桂　南１,杨星团１,屠基元１,姜胜耀１)皇家墨尔本理工大学工程学院,澳大利亚墨尔本　３２．０８３(清华大学核能与新能源技术研究院,北京　１１．０００８４;

摘要:为更好地了解球床中球团与单球的双组分体系的流动特性,本文采用离散单元法模拟分析了球床、、和１夹角.不同数量占比的球团和单球从球床底部的卸料口混合卸出,定量分析了球团６０°９０°１２０°８０°组合夹角和球团数量占比对球床双组分混合球流流动特性的影响.研究结果表明,球团的存在会降低卸料速率,定量分析该降低程度与球团的粘接夹角及球团数量占比的关系后,发现卸料速率与二者均存在一定负相关性,即卸料速率随球团夹角及球团数量占比的增加而减小.关键词:球床;球流;球团;双组分颗粒流;组合角;数量占比:/doi１０．７５３８zk．２０１９．ouxian．０１７５yy

()中图分类号:TL３３３　　　文献标志码:A　　　文章编号:１０００Ｇ６９３１２０２００２Ｇ０２３５Ｇ０７

中球团与单球二元混合颗粒材料的卸料特性.球团由３个粘接的单球通过刚性b形成ond组合而成,

BinarixtureDischarinharacteristicsyMggC

ofPebbleClusterFlowandPebbleFlowinPebbleBed

(１．InstituteouclearandNewEnerechnoloTsinhuaUniversitBeiin０００８４,China;fNgyTgy,gy,jg１

２．SchoolonineerinRMITUniversitelbourneVIC３０８３,Australia)fEgg,y,M１１∗１１２１

,,,GELianGUINanYANGXintuanTUJiuanJIANGShenaog,gygy

,

,

:AbstractInordertobetterunderstandtheflowfeaturesofbinarixturesinpebble　ym

,bedwhichiscomosedofpebbleclustersandpebblesdischarinharacteristicsofthepggc

,,６０°９０°１２０°and１８０°andbeinischaredfromthedischarinoleatthebottomofgdggghebblebedwithdifferentnumberfractions．Thedischarinharacteristicsofpebblepggc

binarixturematerialofebbleclusterandpebblesweresimulatedbEM．ThepebＧympyD

,bleclusterwascomosedofthreeparticlesconnectedbiidbondsforminnlesofpyrggag

clustersandpebblesofdifferentnumberfractionsandincludedanleswerequantified．g

,Theresultsshowthattheexistenceofpebbleclusterscanreducethedischarinateggrandthedereeofreductionisneativelorrelatedwiththeincreaseoftheincludedggycanlesandthenumberfractionofpebbleclusters．Whentheincludedanleandthegg

收稿日期:修回日期:２０１９Ｇ０３Ｇ１９;２０１９Ｇ０５Ｇ０９

)基金项目:国家自然科学基金资助项目(５１５７６２１１

,作者简介:葛　良(男,湖南岳阳人,博士研究生,核能科学与工程专业１９９０—):桂　南,∗通信作者:EＧmailuinan＠mail．tsinhua．edu．cngg

:///网络出版时间:网络出版地址:２０１９Ｇ０８Ｇ２８;htt∥kns．cnki．netkcmsdetail１１．２０４４．TL．２０１９０８２７．１０１１．００２．htmlp

２３６

原子能科学技术　　第５４卷

;includedanlenumberfractiong

,numberfractionofpebbleclustersincreasethedischarinatereduces．ggr

:p;p;bKeordsebblebedebbleflow;pebbleclusterinarixtureparticleflow;ymyw

　　世界上现在运行的模块化高温气冷实验堆

()是由清华大学研制的,于２HTRＧ１０００３年成功发电.得益于球床式堆芯的设计,高温气冷堆堆芯产生的“溢出热量”能自发向环境中排在冷却剂无法循环工作的情况下,堆芯最高温１]

,放[即使发生人为事故和遇到不可抗力天灾,的单元获得整个宏观系统的整体运动状态.D通过检索２个单球球EM在模拟过程中,心间距判断两者是否接触或碰撞,再使用合适的物理模型表征碰撞应变和能量的耗散,使用牛顿第二定律分析接触碰撞中每个单球的受力和运动学信息,从而更新下一时刻单球运动状度也不会超过安全限值,且可自发停止运行[

２

]因此球床式高温气冷堆安全性高,被认为是有.希望满足核能系统要求的堆型.在径的燃料球组成的重力驱动下缓慢的稠密颗粒HTRＧ１０堆芯球床内,

堆芯是由相同直流.本课题组已对等直径的球流运动特性开展了多年的研究,如建立了与实际堆芯尺寸比例

为基础上结合离散单元法１∶５的二维可视化球流实验装置[３Ｇ４

](计算模拟程序[５]

DEM)

开发了对,应在此的随机性和在宏观上具有不连续性等,发现球流运动在微观上具有[６]

验结果和模拟数据的对比验证,论证,了通过实

目前H球区T[７R]

Ｇ.１０的堆芯设计合理、

床内不存在滞留相比于实验,DEM具有能获得颗粒全部运动学信息的优势,因此被广泛应用于粉

体材料的研究[

８Ｇ９]

之前关于堆芯球床内球流运动的机理性研

.究主要是针对单组分球形燃料球,同时也因为球床堆内球流运动存在多体互相接触且具有随机性、间歇性和不连续性等复杂颗粒流动特性,假设极端异常事故发生时燃料球可能发生附着,球团的出现可能导致堆芯球排放卸料流动速度的降低.为进一步深入研究球流运动,对于球团和单球组成的双组分燃料球的流动问题仍有待研究,很多机理性问题尚不十分清楚.鉴于此,本文将研究球床内球团和单球组合球流的流动特性,为反应堆在极端情况下球流特性的预测和分析提供指导.

　计算模型和模拟条件

􀆰１　数学模型

DEM的基本原理[１０Ｇ１１]

就是将球流研究对象离散为的计算单元,通过计算分析每个

态和空间坐标,两个部分[１２]

,具体方程如下DEM基本方程包括平动和旋转

d:

midvti＝iC

j＋Fig

(Iidw􀰐jF１)其中:mdti＝i、Ii、vi和wi分别􀰐jMij(２

)为堆芯球i的质量、

转动惯量、速度及角速度;Fig

为重力;FiCj和Mij分别为堆芯球i受堆芯球j的接触力和转矩.接

触力可分解为法向接触力Fn和切向接触力FtFn＝knFξn－γṅξn(.

t＝ktttFξ－γ̇ξtt(３

４

))其中:kn和kt分别为法≤μFn向和切向的刚性系(数５);n和γt分别为法向和切向的阻尼系数;ξn和̇ξn分别为法向应变和速度;ξt和̇ξt分别为切向应

变和速度;μ为摩擦系数.

在课题组前期的研究[１３Ｇ１５]

中,已在多种球

流工况下验证了􀆰２　计算模型

DEM的准确性.

本文所采用的计算模型与实际实验反应堆球床几何尺寸相同,主要由球床圆HTRＧ１０柱体、锥形底座和卸料口.模拟计算所采用的主要参数列于表３部分组成,如图１所示芯球流是重力驱动下的稠密流,本文使用两种１.堆堆芯球,第１种堆芯球为与实际堆芯球等大,即直径dp由３个单球粘接成的球团＝６０mm的单球.第２种堆芯球为,又根据组成球团的单球粘接的角度不同以及球团所占球床内堆芯球总个数比不同,使用多球模型分别模拟了θ＝６０°、９０°、１２０°和１８０°４种不同的粘接夹角,以及在这些夹角下１６００、４５００∶１３５００、７２Rn００＝∶N５４c∶００Ns和＝９０１８００００∶∶

０

γ１１１２第２期　　葛　良等:球团流与球流在球床内的混合卸料特性研究

２３７

Ns为单球个数,９０００∶０即为球床中仅含球团不含单球.

等４种球团数量占比的工况,Nc为球团个数,

布于绿色单球体系内,局部近壁区可看到球团的直角构型形状.

表１　计算模型主要参数

Table１　Parameterusedinsimulation

参数堆芯直径堆芯高度堆芯锥角球直径de回复系数e杨氏模量E摩擦系数μ粘接夹角θ球团数量占比Rn

数值１􀆰８m３m３０°

６０mm０􀆰９７

７

１×１０Pa

１８００∶２１６００、４５００∶１３５００、()简记为１∶１２、１∶３、４∶３、９０００∶０

７２００∶５４００、９０００∶０

、、、６０°９０°１２０°１８０°

０􀆰３

ofhihtemeraturegasＧcooledreactorgp

Fi．１　Exerimentalinstallationgp

图１　高温气冷堆实验装置

整个卸载过程分为加载堆积和卸料两个步骤,初始状态是球团和单球在球床内自由落体随机堆积,此时卸料口关闭,加载结束到堆积状态稳定后再打开卸料口开始卸球.

口向周围扩散的一种依赖于时间的整体行为,局部交换较少,球床中的顶面始终保持接近水平.在卸球引入空隙和内壁约束的共同影响下,在球床的上部,球团和单球整体均匀地向下运动,在锥形底座附近,边缘的堆芯球朝中心汇聚,再经由卸料口排出.随着时间推进,球床内堆积的顶面逐渐降低.

２􀆰２　球团粘接夹角和数量占比对卸料速率的

影响

为更好地分析球团粘接夹角和数量占比对卸料速率的影响,定义卸球数量占比λt如下:

/λN０t＝Nt()６

从图２可看出,卸料的球流运动是由卸料

２　结果与讨论

２􀆰１　典型卸料过程

图２为粘接夹角θ＝９和Rn＝１∶１０°２时球团和单球混合卸料过程在t的单球＝５、１０、２０s分布.整个卸料过程中,单球和球团用不同颜色表示,绿色代表单球,其余颜色代表由３个单球组成的球团.初始堆积好后,各球团随机分

Fi．２　Dischareofpebbleforθ＝９０°andRn＝１∶１２gg

,)图２　堆芯球卸料过程(θ＝９０°Rn＝１∶１２

２３８

其中:Nt为t时刻已卸载的单球和球团数量之和;N０为初始堆积球床内单球和球团数量之

原子能科学技术　　第５４卷

和.不同粘接夹角及不同球团数量占比卸料过程的卸球数量占比计算结果如图３、４所示.

Fi．３　Percentofdischarinarticlesatdifferentincludedanlesgggpg

图３　不同粘接夹角卸料过程卸球数量占比的变化

Fi．４　Percentofdischarinarticlesatdifferentnumberfractionsgggp

图４　不同球团数量占比卸料过程卸球数量占比的变化

第２期　　葛　良等:球团流与球流在球床内的混合卸料特性研究

２３９

　　图３为球团夹角固定的情形下不同球团数量占比对卸球流率的影响,图中曲线上任意点的导数(斜率)等于该时刻瞬时卸球流率.从图３可清楚地看到,卸球数量占比λt随时间逐渐增大.

从图３可看出,在固定夹角θ下,随着球团

动、小粘接夹角的工况(Rn＝９０００∶０,θ＝６０°

)和θ＝９中,对０°λt随时间逐渐增加.特别地,

),于θ＝１且Rn＝９０图４８０°００∶０时的情况(d可能发生卸球数量占比短暂性持平的情形,表明此时卸料过程发生了间歇性阻塞.这种暂时性的阻塞在重力的驱动下随着时间推进会消失,表现为卸球数量占比恢复增加的情形.因此,在整个球床内不存在滞留区.

综上可见,球团的存在总是起到减缓卸料流率的作用,因此,球团的存在会降低卸料流动性.球团粘接夹角和数量占比对卸球特性的影响可分析.

２􀆰３　球床内堆芯球平均速度分析

为直观显示堆芯球流动速度特性,将球床根据径向距离和轴向高度划分成等大的网格,然后分别依次统计每一时刻内任意网格内球团和单球各自在重力方向的速度分量,最后对时间取均值得到卸载过程中重力方向的球流速度均值v计算公式如下:rＧz的空间分布,

１１)()r,z,t７z(j,􀰐vTNrＧzj＝１

、其中:轴向高度为z的网格内NrrＧz为径向距离为

NrＧz数量占比Rn的增加,这说λt的变化速率降低,明Rn的增加也会导致卸球流率降低.需指出的是,随着Rn的增大,不同粘接夹角的λt斜率降低卸球流率的作用越明显.这是因为球团和单球一样,运动也分为平移和旋转两个部分,但因球团并不具备单球中心对称的特性,球团的存在将会降低自身和周围球团/单球的旋转和

１６]

,平移[且球团的体积为单球的３倍,球团的存

的差异也逐渐增大,这说明随着球团数量增多,

/球团的挤压、在也会增加对周围单球碰撞、滑移和反弹等作用,加强相互之间的竞争和阻滞,因此球团的存在会降低整体的球流速度从而降低卸球流率.

角θ的变大,曲线的斜率逐渐减小,意味着增大

的球团的内夹角也起到了延迟或降低球床卸球

[７]

流率的作用.根据T的研究,当卸料口o等１

由图４可看出,在固定Rn下,随着粘接夹

vrＧz＝

直径大于球直径的５倍时,不会形成堵塞拱,在混合球流中,单球起到了润滑的作用,球流速度减缓,但不会出现堵塞现象.然而在纯球团流

球团或单球的数量;vz为对应球团或单球在重力方向上的速度分量,其计算结果如图５所示.

速度越大;在轴向高度上,距离卸料口越近,速

从图５可看出,在径向上,距离中心越近,

Fi．５　Satialdistributionofmeanvelocitfpebbleflowdurinischaringpyogdgg

图５　卸载过程中球流速度均值的空间分布

２４０

度越大.在球床圆柱本体空间内因距离卸料口较远而受到的影响较小,因此球流运动速度较小.在卸料出口附近,球流速度急剧增可发现,球团和单球之间的速度差异不大,这是因为在重力驱动下,球团和单球因为多体间的互相接触,存在驱动、阻滞和竞争等互相影响的机制,因此球流运动呈现整体性运动趋势,导致球团和多球有着相近的速度.且在球床的大部分空间内,重力方向上球流速度的空间分布与粘接角度和数量占比之间的关系不大,只有卸料口及附近的空间内,才出现了速度分化.因此,对于双组分的球床及不同的球团形状,球床本体内仍保留了较好的整体流流态特征.

为定量分析球团粘接夹角和数量占比对球、、需指出的是,通过分别对比图５acd可发现,、大,并在卸料口处达到最大值.对比图５ab

原子能科学技术　　第５４卷

、Rn＝１∶１２、１∶３、４∶３以及Rn＝４∶３,θ＝６０°

、、下重力方向的球流速度均值v９０°１２０°１８０°p,计算公式如下:

１１)()vvz,t８z(j,p＝

TNp􀰐１j＝

;其中:T为总的计算时间,sNp为颗粒数;t为

.瞬时时间,s的结果表明,球团和单球之间的平均速度存在差别,但绝对值的变化不大,球流呈整体运动趋势.对比单球和球团在相同夹角下不同数量占)比的情形(图６和在相同数量占比而不同夹a),角下的情形(图６可发现夹角的变化对整体b流动速度的影响更大,这意味着相比于球团数量占比的增大,球团夹角的增大会造成更大程度球流流动速度的降低.

速度均值vp计算结果如图６所示.定量

Np

,流速度的影响,沿轴向高度分别统计θ＝６０°

Fi．６　Distributionofdischarinelocitfpebbleflowwithheihtgggvyog

图６　球流速度沿高度的分布

３　结论

卸料进行了定量分析,讨论了球团粘接夹角θ和数量占比Rn对二元混合组分球流卸球流率的影响,得到如下结论.

)球团和单球的混合卸料过程为随时间１

变化的整体性流动,局部混合较少,球床圆柱本体内球流保持近似水平下降且堆芯的底部不存在滞留区.

)球团的存在会降低卸料速率,卸料速率２

降低的程度与粘接夹角及球团数量占比的增大呈正相关趋势.需注意的是,在θ＝１且R８０°n＝

本文采用DEM方法对球团和单球的混合

可能会出现间歇性的堵塞,但堵塞９０００∶０时,

只是暂时性的.

)球流运动的速度变化趋势是轴向越靠３

近卸料口,球流的速度越快;径向越靠近球床中心,球流的速度越快.因为球流呈整体性运动,同一工况下球团和单球在径向距离和轴向高度上的分布有一定区别,但相差不大,球流速度与粘接夹角和数量占比之间的关系也不大.参考文献:

[]１　KADAKAC,BERTEMV．AdvancedmoduＧ

[]laritesinfortheMITpebblebedreactorJ．ydg

第２期　　葛　良等:球团流与球流在球床内的混合卸料特性研究[]２　PebbleBedModularReactorLimited．PBMR

:B[R]．SouthAfricaPBMR,１９９９．,():NuclEnes２００６,２３６５Ｇ６５０２Ｇ５０９．gD

２４１

mechanicalinteractionsofgranularmediaina

,０/safetnalsisreort０１９２９Ｇ２０７４,Rev．yayp

[]hoerusinhediscreteelementmethodJ．ppgt[]T１０SUJIY,TANAKAT,ISHIDAT．Laraniangg

,PowderTechnolo２０１５,２７８:２８６Ｇ３０５．gy

[]３　YANGXT,HUWP,JIANGSY．ExerimenＧp

desinedpebbleＧbedhihＧtemeraturegasＧcooledggp

talinvestiationonfeasibilitftwoＧreionＧgyog[],:reactorJ．JNuclSciTechnol２００９,４６(４)

numericalsimulationofplulowofcohesionlessgf]articlesinahorizontalpie[J．PowderTechＧpp,():nolo１９９２,７１３２３９Ｇ２５０．gy

[]４　JIANGSY,YANGXT,TANGZW,etal．

ExerimentalandnumericalvalidationoftwoＧp

３７４Ｇ３８１．

[]C１１UNDALLPA,STRACKODL．Adiscrete

,Geotechniue１９７９,２９:４７Ｇ６５．q

[]numericalmodelforgranularassembliesJ．

[]C１２UNDALLPA．AcomutermodelforsimulaＧp

rceogrieo[nJＧ]designedpegbbDlebe２dr０１e２a,ct２o４rw６:i２t７hd７Ｇ２y８n５a．mic[５]　TtHORN．TNOucNlECn．Onest,herelationshipb

etweenehneergmyoodfultuhsecoofnsptiatrutiecnutlaptaerticmleedsiaJa．nJPdshuy

rsDfaceApplPhy

s,１９９３,２６(１０):１５８７Ｇ１５[]９１．:[６]　CHOIJ,KUDROLLIAgi,ROSALy

EＧSRRdriven,edtenalse．

gDriaffnuusliaorfnlaonwds[Jm]ix．PinhysnRgervaviLt(１７):１７４３０１．

ett,２００４,９２[７]　杨星团,刘志勇,胡文平,等运动的唯象学DEM模拟[J．H]．原TR子_１能０堆芯球流

科学技术,

２Y０A１３,４７(１２):２２３１Ｇ２２３７．ceotraelN．GgbDyyEXiSpMsngtheniuomanmul,eantiLooIlonoUgicfZhalpieybombngelgy

te,hfoldoHUwi[J]nHWen,．４７ATpt(１oRiＧn２m１g

i)c０

,[８]　C２２Ene３r１Ｇ２２c３i７en(icnCeahnidTneseec)h．

nolo,２０１３:mOetEhTodaZEnECdth．eCeaflfiebcrtatoiofnoparftitchleesdihsacp

reeteJel．ePmoewnt

derTechnology

,２０１６,２９７:５０Ｇ７０．[]Ｇ[９]　HstuÖdHNyoEntRheiD,nflWuIeRncTeoZfSp,arStiCclHesERhaEp

ReoV．nthA

etriaonRlScgpksroockymysgptroeessivelargescalemovementsinblockysmiu[mC]o∥PftrhoeIceednitnerg

nsaotifotnhaelISnotecirentaytiofnoＧr１３SocietyMf

eocparticlesiGrRhanzeaUINoiccks．Nndr,MeYecagiAhnoNacny

nGic,sFwiX,rdthoT１a９n,７cnte１eta．:International[]WUH,heml．Eifxfiencgtsoa

ndf

d[iJs]pe．rGsiroanonulfarMpebabtlteesir,nt２０１w６or,１８eg(i４o)n:７p６eb．blebed[１４]JmIeAntXaL,GUIN,YANGXT,etal．Exp

eriＧaflnodwiinlntesartmusiidltyobteannecdady[Jofvnal]．eyNrsiuysoclslEofvnwaelogDncedditys,e２nc０so１erre６pl,ea３bti０bo５len６２６Ｇ６３８．

:[１５]tLiIYono,nGpUeINbble,flYoAwNpGXT,etal．EffectofafcrticoＧ

[J]．AnnNuclEnergay

tt,e２r０n１i６n,p９e４b:bl３eb２Ｇ４e３dr．er

[１６]aGlR．AEnVtIaSnHgleNd,FRANKLINSV,HUDL,et

１７]２０１２,１０８(２０)gr:a２n０u８l０a０rm１．

edia[J]．Phy

sRevLett,[gT[rOJa]n．uKlPa,hrLy

spAiaIPcraAticleY,si,P２００nA２atKH,３１wK．Flowandjamof

５o:Ｇd１i７m４Ｇe１n８si０o．nalhopper