收藏 分销(赏)

干法后处理熔盐电解精炼过程数学模型研究.pdf

上传人:自信****多点 文档编号:3110573 上传时间:2024-06-18 格式:PDF 页数:10 大小:4.15MB
下载 相关 举报
干法后处理熔盐电解精炼过程数学模型研究.pdf_第1页
第1页 / 共10页
干法后处理熔盐电解精炼过程数学模型研究.pdf_第2页
第2页 / 共10页
干法后处理熔盐电解精炼过程数学模型研究.pdf_第3页
第3页 / 共10页
亲,该文档总共10页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

1、第58 卷第1期2024年1月原子能科学技术Atomic Energy Science and TechnologyVol.58,No.1Jan.2024干法后处理熔盐电解精炼过程数学模型研究王赛,林如山,李康祎,钟振亚,钱宾杰,张磊,唐洪彬(中国原子能科学研究院放射化学研究所,北京10 2 413)摘要:熔盐电解精炼是乏燃料干法后处理的核心工艺单元,通过数学模型探索高温熔盐电解精炼过程的化学与电化学变化,可为电解精炼工艺优化和设备设计提供参考依据。本文基于电化学热力学及物质传递公式建立了乏燃料熔盐电解精炼过程的数学模型,以铀环锆三元合金燃料为研究对象,计算了燃料中关键元素的电极电势、分电流及

2、物料分布随时间的变化。采用向后差分法对物料分布变化方程进行离散,通过文献实验数据对建立的数学模型进行了准确性验证。结果表明,模拟计算所得阴极沉积铀产品与实验数据的相对误差为2.8 0%,所建数学模型具有较好的拟合性。同时采用所建模型模拟计算了电流强度对乏燃料电解精炼过程的影响,结果表明电解速率与电流强度呈正比,不改变环铀锆的溶解和沉积顺序。关键词:乏燃料;干法后处理;熔盐电解精炼;数学模型;物料分布变化中图分类号:TL244doi:10.7538/yzk.2023.youxian.0346WANG Sai,LIN Rushan*,LI Kangyi,ZHONG Zhenya,QIAN Binj

3、ie,ZHANG Lei,TANG Hongbin(Department of Radiochemistry,China Institute of Atomic Energy,Beijing 102413,China)Abstract:Molten salt electrorefining is the core process unit of spent fuel dry reprocess-ing,and the chemical and electrochemical changes of high-temperature molten salt elec-trorefining pro

4、cess can be explored through mathematical model research,which canprovide a reference basis for electrorefining process optimization and equipment design.In this paper,based on electrochemical thermodynamics and material transport formu-las,a mathematical model of the electrorefining process of molt

5、en salt in spent fuel wasestablished,and the changes of electrode potential,partial current and material distribu-tion of key elements in spent fuel with time were calculated.The central differencemethod was used to discretize the material distribution change equation and obtain theoscillation resul

6、t.The backward difference method of discrete material distributionchange equation was used to verify the accuracy of the established mathematical modelthrough literature experimental data.The results show that the relative error between文献标志码:AMathematical Model of Dry ReprocessingElectrorefining of

7、Molten Salt文章编号:10 0 0-6 931(2 0 2 4)0 1-0 0 2 3-10收稿日期:2 0 2 3-0 5-11;修回日期:2 0 2 3-0 7-15*通信作者:林如山24simulated uranium product deposited at the cathode and the experimental data is2.80%,which shows that the mathematical model has good fit.During electrorefining,there are two mutations in the anode p

8、otential,indicating complete dissolution of Pu andcomplete dissolution of U separately.At the same time as the electrochemical dissolu-tion of Pu,a displacement reaction occurs with UCl,in the molten salt,so that theanode division current of Pu is larger than the input current.The order of dissoluti

9、on ofthe three elements at the anode is PuUZr.U elements are deposited in large quanti-ties at the cathode,and the deposition rate slows down after Zr dissolved at the anode.The concentration of elements in molten salt varies with the electrorefining process.Inthe simulated time range,the concentrat

10、ion of U element decreases to a certain extent atfirst,and the concentration of Pu element increases with time,and then the concentra-tion of the two elements remains almost unchanged.The molar concentration of Zrelements in molten salt remains traces.At the same time,the influence of current inten-

11、sity on the spent fuel electrorefining process is simulated and calculated.The resultsshow that with the increase of current intensity,the anode potential shows an overallincreasing trend,and the cathode potential shows an overall decreasing trend.The rateof anodic dissolution and cathodic depositio

12、n of spent fuel accelerates,and the partialcurrent of each element increases.The change rate of element concentration in moltensalt is accelerated.In actual working conditions,the current can increase within thepermissible range to improve cell efficiency.Key words:spent fuel;dry reprocessing;molten

13、 salt electrorefining;mathematicalmodel;varieties in mass distribution干法后处理是闭式循环快堆核能系统的关键环节-2 ,是高燃耗、高环含量、强放射性快堆乏燃料后处理的现实技术选择。熔盐电解干法后处理技术具有介质耐辐照、临界风险低、工艺流程短、废物量小等特点,适应性更高、处理对象更广,是最有应用前景的干法后处理技术3-41。美俄均已掌握快堆乏燃料干法后处理技术,并分别基于本国燃料循环策略建立了经工程规模热验证的干法后处理流程5-6 。近年来,我国也已围绕快堆闭式燃料循环,开展了熔盐电解干法后处理技术研究,并取得了较大进步7-8

14、 。熔盐电解精炼是乏燃料干法后处理的核心工艺单元,其反应过程复杂、腐蚀工况严苟、高温过程不可直接观察,同时涉及多相界面的流动及传质问题 1。通过理论模拟和计算分析,研究熔盐电解精炼反应器内的物理及化学变化规律,可为电解精炼工艺优化和设备设计提供参考依据。国外研究人员已利用先进计算机技术,开展了大量研究。如Hoover等10-11基于化学动力学、化学热力学和物质传递模型建立原子能科学技术第58 卷了Mark-IV电解精炼装置的模型,模拟了铀环锆三种元素在阳极的溶解情况,但该模型不能正确使用过电位。Zhang121基于扩散控制假设开发了一种电解精炼过程动力学模型,并用美国阿贡国家实验室(ANL)的

15、实验室数据验证了从液态镉阳极到固体阴极物料变化数据的准确性,该模型考虑了液态电极的体积变化使模拟过程更加精确。Yos等13 基于扩散层模型和Bulter-Volmer公式建立了电解精炼过程动力学模型,模拟所得物料变化与ANL的实验数据拟合良好,并进一步研究了模型结果对不同输人参数的敏感程度,提出不同输入参数应该考虑不同的扩散层厚度这一结论,该模型缺乏对于固体阳极的适应性研究。国外建立的液态镉阳极的电解精炼模型日趋成熟,对于固体电极的物料变化模拟也有一定进展,但基于物质传递原理对于固体电极的电极电势以及分电流的模拟报道较少。本文基于电化学热力学及物质传递公式建立乏燃料熔盐电解精炼过程的数学模型,

16、开展第1期王赛等:干法后处理熔盐电解精炼过程数学模型研究铀环锆三元合金燃料中关键元素的电极电势、分电流及物料分布计算,并采用有限差分法对物料分布变化方程进行离散,通过文献实验数据对建立的数学模型进行准确性验证。1数学模型构建乏燃料电解精炼的电化学行为主要包括乏燃料在阳极溶解、熔盐电解质中的物质传递以及阴极沉积过程,金属乏燃料电解精炼过程原理如图1所示。在实际工程应用中,通过控制电极电位使贵金属裂片产物和Zr保持不溶,形成阳极泥从阳极篮掉落进入Cd相。本文为探究电解精炼规律,在铀环溶解后继续考虑Zr的溶解。其中涉及到的电化学知识主要包括电池热力学、电极反应动力学以及由迁移和扩散引起的物质传递行为

17、14,并在此基础上建立数学模型。阳极吊篮+人人碱金属和碱土金属元素部分稀土FP部分稀土FP贵金属FP燃料合金中的ZT图1金属乏燃料电解精炼过程原理示意图Fig.1 Schematic diagram of electrorefiningprocess principle of metallic spent fuel1.1电化学热力学以U3+为例,其在阳极发生的反应为:UU3+3e根据热力学可逆性,通过Nernst方程可以得到电极电势E9和电极过程参与者浓度之间的关系:Eeq=EnFaR其中:E为标准电势,V;R为气体常数,取8.314J/(molK);T为开尔文温度,K;F为法拉第常数,取9

18、6 48 5C/mol;ao和aR分别为电解质中氧化物(O)和还原物(R)的活度。活度与电极状态有关,电极电势因组分不同有不25同的表现形式,本文考虑的两种电极均为固体电极。假设元素在固体阴极的活度为1,则电极上的Nernst方程可转化为:EW=E+RTIn onF参与阳极反应的各元素的界面电势差相同,等于阳极电势;参与阴极反应的各元素同理。控制方程可表示为:E.=Ei.=E;j=.=E.aE.=Ej.=.=E.其中:下标a和c分别表示阳极和阴极;i、j、n分别代表参与电极反应的元素。1.2物质传递及电流溶液中的物质传递表现形式为扩散、迁移及对流,假设在电极附近电活性物质的传递全部为扩散形式。

19、根据Fick定律,熔盐电解质与电极界面处元素的摩尔通量(f)可表示为:固体阴极液态镐阴极人人人人人人人UTRUJU,RECd(1)RT1nao(3)(4)(5)acf=-D;AaN其中:N为熔盐/电极界面法线;D;为元素i的扩散系数,cm/s;A 为界面面积,cm;c;为元素浓度,mol/cm。假设电极附近扩散层中离子仅在法线方向进行扩散传质。基于量纲分析法,用传质系数(K,cm/s)替代式(6)中扩散系数与法线位置偏微分离散形式之商,则在熔盐和固体电极的界面处Fick定律可简化为:f=-A.Km(c-c.)fme=A.Km(c-cms.)其中:ms表示介质在熔盐中;b代表本体溶液;K和K分别

20、为阳极侧与阴极侧元素在熔盐中的扩散系数;c为元素在熔盐中的本体浓度;cs.和cms.分别为元素在阳极侧和阴极侧与熔盐界面处的浓度。基于熔盐中的物质传递,带电荷的元素离子在阳极和阴极产生的分电流为:Ia=nFfms(2)I=nFfms控制方程为阳极和阴极处的总电流等于各元素分电流之和,阴极电流和阳极电流相等。其表达式为:NI.=I=I=I.=I(11)(6)(7)(8)(9)(10)N(15)261.3物料变化方程本体溶液中元素摩尔浓度随时间变化的微分方程为:dciLi.a-Ii.dtnFVms其中,Vms为熔盐体积,cm。物质在阳极溶解和阴极沉积的物质的量的变化为:dM=Ldtn;F其中,M,

21、为元素i在电极侧的物质的量,mol。2电解精炼过程模拟求解电解精炼过程数学模型的建立通过计算机编程实现。在程序开始前,首先通过假设各元素在阴极与熔盐界面处的浓度为0 mol/cm,计算阴极允许最大电流:NLe.mx=n,FKpcr.1式(14)等式右侧各参数均为已知。若施加的电流大于阴极允许最大电流,将发生预测之外的界面反应,应增大熔盐中元素的初始浓度或降低所施加的电流以满足该限制条件。判定满足电流限制条件后,熔盐性质、元素性质、电解精炼阴阳极尺寸以及物料阳极初始组成及熔盐中元素初始浓度均可通过已知实验-2.5a-2.6F八/审-2.7-2.850C1.00.50.02-0.5-1.050a-

22、电极电势;b-阳极元素物料变化;c-阳极元素分电流;d熔盐物料变化图2 中心差分法离散物料变化方程的振荡解Fig.2 Oscillation solution of discrete material change equation by central difference method原子能科学技术第58 卷数据或计算15-19 1得到。算法流程如下。1)在一个合适的范围内估计某时刻阴阳极电极电势,并基于该电势通过式(1)(2)、(6)(10)计算电极表面各元素的浓度以及电(12)极与熔盐界面的分电流,通过电流控制方程反向利用二分法拟合逼近阴阳极电极电势。2)通过有限差分法离散物料变化微分

23、方(13)程,得到下一时间节点各元素在熔盐中的浓度及摩尔质量数据。分别采用向后差分法和中心差分法离散方程。中心差分法离散熔盐中元素浓度变化:CZ.b.+1 c.b.-12t其中,t为时间步长,s。中心差分法离散阴阳极物质的量变化:M.a.+1 M.a.-(16)2t(14)2t模型所使用函数数值之间有间断,如某元素在阳极全部溶解后电流突变为OA。中心差分法在处理此种函数时,具有易振荡不收敛特性。在模拟过程中,发现无法得到收敛解,电极电势、分电流以及物料变化的振荡解如图2 所示。阳极0.10rb阴极0.05F0.00F200400时间/minUPuZr200400时间/minn;FVmsn;Fn

24、.FUPuZr6000d(Tu/louu)/卖1.00.5F0.06000(17)200400时间/minUPuZr200400时间/min600600第1期向后差分法精度为一阶精度,其离散方程如下。向后差分法离散熔盐中元素浓度变化:C2b.+1-cb.=Lia-li.eattn;FVms同样的方法离散阴阳极物质的量变化:Mi a.-+1-M.a(19)tn;Ftn;F3)时间步长增加1个单位长度,计算传质系数与阳极面积在当下时刻的值。传质系数与阳极面积随电解精炼过程均有不同程度减小,传质系数随时间的变化采用线性函数、阳极面积采用指数函数拟合11。返回步骤1。4)达到设置的反应时间时,停止程序

25、,并输出阴阳极电极电势、分电流及元素的物质的量在电极侧和熔盐中随时间变化的数据。程序操作流程如图3所示。3结果与讨论3.1电解精炼模型结果与验证本文采用日本中央电力研究所(CRIEPI)的实验研究数据2 0 1来验证模型的准确性。电解精炼器阳极采用9.6 g三元合金U-20%Pu-10%Zr(质量分数);阴极采用固体不锈钢阴极。容器中放入1kg熔融电解质,其中U、Pu、Zr的初始质量分数为0.7 3%、3.1%、0.0 0 3%。在7 53K下,采用恒电流0.4A模拟电解精炼过程。其余输人参数与实验条件保持一致。所需物性参数根据文献15-19 或计算得到,如表1所列,其中各元素标准电势以Cl2

26、/CI-为参比电极。模拟数据与实验数据的对比如图4所示。图4a模拟数据显示阳极电势在30 min附近发生突变,之后呈稳定趋势。模拟数据较实验数据低,这是由于忽略了实验过程存在的电池电阻。由图4b可见,在电解初期,U元素浓度下降,Pu元素浓度升高,之后二者浓度几乎保持不变。整个过程中Zr元素在熔盐中的摩尔浓度始终较低。与实验数据对比可发现,U元素摩尔浓度较真实值高,Pu元素较真实值低,这是因为熔盐中元素浓度对于标准电极电势非常敏感9,3%的变化会导致熔盐中U元素的量王赛等:干法后处理熔盐电解精炼过程数学模型研究参数poms,g/cm3Eu3+/u,VEPa+/Pu,VEzA+/z.VKicm/s

27、0.001.9KmsKPu3+,cm/sco.-mol/ml 2.41 9 10-l发生50%的偏差。通过电极电势相同以及阴阳极电流相等控制方程得到U元素阴极分电流,进一步根据物料变化方程得到阴极沉积的元素的物质的量。将同一时刻阴极沉积的U计算值0.0 2 7 8 mol与实验值0.0 2 8 6 mol对比,27输人:温度、电流、电极/电解槽几何形状、传质系数、燃料成分否(18)Ic,maxn,FKmcms?Ci,6工是计算:熔盐性质、电极性质、初始条件、表观电位(20)计算:阳极表面积及传质系数估计:阳极和阴极电位计算:电极表面浓度、各元素分电流、阳极电流、阴极电流否1,=1=1?是计算:

28、电极上的物质种类和元素浓度否-I+1输出结果并结束程序图3电解精炼程序流程图Fig.3Program flowchart of electrorefining process表17 53K下U、Pu、Zr 元素的物性参数Table 1 Physical property parametersof U,Pu and Zr at 753 K数值1.6289-2.5007z.mol/mL.ms-2.769 9Mu.a,mol-2.1596Mpu.a,molMzr.a.mol0.006.2Aa,cm20.0016Ac.cm2Vms,cm3一操作设定时间?是参数数值9.0959X10-48.8025X1

29、0-70.028 20.00800.01055.398547.7836613.919328计算得到相对误差为2.8 0%。模型充分考虑了U在阴极沉积的过程机理,具有较好的拟合效果。电流为0.4A时模拟得到的阴阳极电极电势、各元素分电流及物料随时间变化示于图5。由图5a可见,电解精炼器阳极电势存在两次突变。第1次突变表明阳极三元合金的Pu元素完全溶解,U元素开始溶解;第2 次突变表明阳极U元素全部溶解,Zr元素开始溶解。电极电势基于Nernst方程计算,模拟结果符合金属活动性顺序表。阴极电势随时间略有减小,初始时刻阴极电势大于阳极电势,之后小于阳极电势。图5b中,前30 minPu元素在阳极的分

30、电原子能科学技术第58 卷流从1.4A线性递减,元素U的分电流从一1A线性递增,意味着U在阳极原位沉积,该负值电流由置换反应产生;之后Pu元素的分电流突变为0 A,U元素的分电流突变为0.4A。这一时刻对应阳极电极电势的第1次突变。在439min时刻,U元素分电流突变为0 A,Zr 元素的分电流在此刻从0 A突变为0.4 A。该时刻对应阳极电势的第2 次突变。图5c显示,前30 min内,Pu首先在阳极溶解,同时U元素在阳极沉积,质量增加。因为Pu元素除在电流驱动下发生电极反应外,还同时与熔盐中的UCl3发生了置换反应,其化学反应方程式如下:a-2.02.2八/审-2.4-2.6-2.8E05

31、0100150200250300 350时间/mina-图4模拟数据与实验数据的对比Fig.4Comparison of simulated data and experimental data1.0元0.8(T/low)/卖-模拟数据0.6F实验数据0.40.210.0日11阳极电势;b一一熔盐中元素物料变化b050100150200250300350时间/min模拟实验UPuZa阳极-2.2F。阴极2-2.4-2.6-2.800100200300400500600时间/min0.4斤0.03oPuo/一明0.020.010.004口口口口口口口口0100200300400500600时间/

32、min一电极电势;b阳极元素分电流;c图5电流为0.4A时的模拟结果Fig.5SSimulation results at 0.4 A1.5bO-O-O1.0V/T0.5会一0.0-0.5-1.0元0100200300400500600时间/min0.040.3xUPuZr0.10.00100200300400500600时间/min阳极元素物料变化;d-阴极元素分电流;e阴极元素物料变化UPuZre0.030.020.010.000100200300400500600时间/minU.PuZr第1期Pu在固体电极的表观电位为一2.7 7 V,远小于U的表观电位一2.5V,此置换反应可以发生。3

33、0 min后,Pu完全溶解,U开始溶解,至439 min时完全溶解,溶解质量与时间呈线性关系。439 min后,惰性金属Zr开始溶解,600min后在阳极侧完全溶解。由图5de 可知,439 min前阴极侧元素分电流大部分由U元素产生。在0 439 min时段,U元素分电流逐渐增大,Zr元素分电流由于浓度降低而不断下降。439 min后,Zr元素由于在阳极溶解浓度升高,分电流逐渐增大,532min后大于U元素分电流,是电流的主要贡献者。U元素最先在阴极大量沉积,沉积量随时间呈线性变化。同时Zr元素在阴极低速沉积。439 min后,U沉积速率变慢,沉积曲线呈凸函数,Zr元素沉积速率加快。因此为控

34、制金属阴极得到较纯的铀枝晶,应严格控制通电时间,使阳极Zr元素不能溶解留在阳极篮中,同时熔盐中不能出现元素Zr。-2.1a-2.2-2.3-2.4-2.5-2.6-2.7-2.890Fig.6 Anode potential and cathode potential under different current intensities00.4A0.030.8 AA2A口4A0.020.01F0.0L0100200300400500600时间/minFig.7 Anode element material change under different current intensities王

35、赛等:干法后处理熔盐电解精炼过程数学模型研究UCl;+Pu PuCl;+U最小。-2.545b-2.550F00.4A-2.555个0.8 A42A04A阳极1100200300400500 600时间/min图6 不同电流强度下的阳极电势和阴极电势0.01200.4A80.8 A42Ao/一望64U200100200300400500600时间/min图7 不同电流强度下阳极元素物料变化29(21)3.2电流强度对电解精炼过程的影响在电解精炼设备放大过程中,电流强度也将随之增大至百安培级别。将电流强度设置为0.4、0.8、2、4A,时间截止点为三元合金全部溶解的前一时刻,分析不同电流强度对电

36、解精炼过程的影响,结果示于图6 11。在模拟过程中,阳极合金不断溶解,面积变化采用前文所述指数函数拟合,固体阴极面积不变。由图6 a可见,随着电流强度的增大,阳极电势整体呈增大趋势,包括Pu溶解时段、U溶解时段以及Zr溶解时段。同时电流增大会使2次电势突变提前,元素溶解的速率加快。因此在放大设备中增大电流有助于阳极乏燃料的溶解。同时不能超出电解精炼过程所允许通过的最大电流限制,若电流大于阴极允许最大电流,可能发生预料之外的化学反应。由图6 b可见,阴极电势整体呈降低趋势。Pu、Zr、U 阳极溶解过程阴极电势随时间逐渐线性降低,其中Pu元素斜率最大,U元素斜率2-2.5602.565-2.570

37、5-2.575口-2.5800口4APu00.40.8 A42A口4A阴极1100200300400 500600时间/min0.010F00.4A0.0080.8 A42A中0.0060.0040.002F0.0000100200300400500600时间/min口4AZr30543V/审210日0100200300400500600时间/minFig.8Partial current of anode element under different current intensities0.050.040.030.020.010.00元0100200300400500600时间/minF

38、ig.9Cathode element material change under different current intensities4.0P3.53.02.52.01.51.00.50.050100200300400.500600时间/minFig.10Partial current of cathode element under different current intensities00.4A0.21*0.8A2A口4A0.200.190.180100200300400500600时间/min图11不同电流强度下熔盐中的元素浓度变化Fig.11Element concentr

39、ation change in molten salt under different current intensities原子能科学技术第58 卷5500.4A00.4A0.8 A0.8 A42A44A3U2100100200300400500600时间/min图8 不同电流强度下阳极元素分电流1.00.6F0.200.4A0.8 A2A口4AU1图9 不同电流强度下阴极元素物料变化1.00.500.4A00.4A0.8 A0.8 A42A0.5F4A0.0U-0.5-1.05020406080100时间/min图10不同电流强度下阴极元素分电流(-0.92U0.9102040608010

40、0时间/min00.4A40.8 A42A42A4A台4A3V/0审2Pu0100200300400500600时间/min500.4A0.8 A42A口4A-0.2-0.6-1.05020406080100时间/min0.4F42A4APu00.4A0.8 A42A04APuZr00.4A0.8 A+A2Aouu/一明口4A32Pu10.30.20.1F0.00100200300400500600时间/min0.016O0.4A0.8 A0.0122A口4A0.080.040.000100200300400500600时间/minZr0100200300400500600时间/min00.4

41、A0.8 A42A口4AZTZr第1期王赛等:干法后处理熔盐电解精炼过程数学模型研究由图7 可见,随着电流强度的增大,各元素完全溶解的时间提前。电流扩大2 倍后,铀阳极完全溶解时间缩短近1/2。电解精炼初期,U元素在阳极沉积的速率不变,Pu元素在阳极溶解的速率随电流的增大逐渐加快,Zr元素溶解速率与电流强度呈正比。根据模拟结果,在实际电解精炼过程中可以适当提高电流强度,以提高乏燃料在阳极溶解的速率。由图8 可见,随着电流的增大,Pu溶解的分电流增大,U在阳极沉积的分电流增大;U溶解和Zr溶解的电流等于输人电流,U溶解结束后,Zr元素溶解电流突变为输人电流。由图9 可见,随着时间的增加,U元素大

42、量沉积;阳极Zr元素溶解后,Zr阴极沉积速率增大,电流强度增加,U元素沉积速率加快;Pu元素在阴极沉积的物质的量为0 mol。以上结果表明,增大电流能增大U元素在阴极的沉积速率,因此可在实际应用中增大电流。由图10 可见,Pu元素的分电流为0 A,U元素分电流占绝大部分。随着Zr元素在阳极的溶解,其分电流逐渐增加。电流增大导致Zr元素的分电流增大速率以及U元素的分电流减小速率加快。电流强度为0.4 A时,在电解精炼后期Zr元素分电流大于U元素分电流,电流强度为0.8、2、4A时U、Zr 元素分电流没有出现交点。由图11可见,U元素摩尔浓度在电解精炼初期减小,之后略有增大,在Zr开始溶解后继续降

43、低。随着电流强度的增加,变化速率加快。Pu元素的浓度在不同电流强度下增加至一定量后保持不变,电流越大,增大速率越大,所需时间越短。Zr元素的摩尔浓度先减小后增大,增大速率与电流强度呈正相关。4结论本文采用有限差分法及二分法建立了铀环锆三元合金燃料熔盐电解精炼数学模型,计算了电极电势、分电流以及物料分布随时间的变化,通过文献实验数据验证了程序的准确性。1)采用向后差分法离散物料变化微分方程建立的电解精炼数学模型可得到收敛结果,避免了中心差分法所建模型计算时的振荡现象。2)模拟计算所得阴极沉积铀产品与实验31数据的相对误差为2.8 0%,所建数学模型具有较好的拟合性。3)电解精炼过程中,环铀锆依次

44、在阳极溶解。若锆发生溶解,将在阴极优先沉积,进入铀产品。可根据对铀产品收率和净化的要求,合理控制电极电势。电流强度与电解速率呈正比关系,但不改变环铀锆的溶解和沉积顺序。参考文献:1徐快堆和我国核能的可持续发展J中国核电,2 0 0 9,2(2):10 6-110.XU Mi.Fast reactor and sustainable nuclearenergy development in ChinaJ.China NuclearPower,2009,2(2):106-110(in Chinese).2林如山,何辉,唐洪彬,等我国乏燃料干法后处理技术研究现状与发展J原子能科学技术,2020,54

45、(增刊):115-12 5.LIN Rushan,HE Hui,TANG Hongbin,et al.Progress and development of dry reprocessingtechnology of spent fuel in ChinaJ.Atomic En-ergy Science and Technology,2 0 2 0,54(Su p p l.):115-125(in Chinese).3刘学刚。乏燃料干法后处理技术研究进展J.核化学与放射化学,2 0 0 9,31(增刊):35-44.LIU Xuegang.Research on dry reprocessin

46、gtechnology of spent nuclear fuelJ.Journal ofNuclear and Radiochemistry,2009,31(Suppl.):35-44(in Chinese).4LI S X.Anodic process of electrorefining spentnuclear fuel in molten LiCl-KCL-UCl,/Cd systemJ.Proceedings of Electrochemical Society,2002.doi10.1149/200219.0541.51BENEDICT R W,SOLBRIG C,WESTPHA

47、LB,et al.Pyroprocessing progress _at IdahoNational LaboratoryC/GLOBAL 2007:Ad-vanced Nuclear Fuel Cycles and Systems.S.1.:s.n.,2007.6SKIBA O V.Technology of pyroelectrochemicalreprocessing and production of nuclear fuelCJ/lProc Int Conf on Future Nuclear Systems,GLOBAL93.S.1.J:s.n.,19 9 3.71LIU K,TA

48、NG H B,PANG J W,et al.Electro-chemical properties of uranium on the liquid galli-um electrode in LiCI-KCl eutecticJ.Journal ofthe Electrochemical Society,2 0 16,16 3(9):D554-D561.328LIU Y C,LIU Y L,ZHAO Y,et al.A simpleand effective separation of UO,and Lnz Os assis-ted by NH,Cl in LiCI-KCl eutectic

49、J.Journal ofNuclear Materials,2020,532:1-9.9NASH K L,LUMETTA G J.Advanced separa-tion techniques for nuclear fuel reprocessing andradioactive waste treatmentMJ.Berlin:Elsevi-er,2014.1o HOOVER R,PHONGIKAROON S,LI S,etal.A computational model of the Mark-IV elec-trorefiner:Phase I-Fuel basket-salt int

50、erfaceJ.Journal of Engineering for Gas Turbines andPower,2009,131(5):054503.11 HOOVER R,PHONGIKAROON S,SIMPSONM.Development of computational models for theMark-IV electrorefiner-effect of uranium,pluto-nium,and zirconium dissolution at the fuelbasket-salt interfaceJJ.Nuclear Technology,2010,171(3):2

展开阅读全文
相似文档                                   自信AI助手自信AI助手
猜你喜欢                                   自信AI导航自信AI导航
搜索标签

当前位置:首页 > 学术论文 > 论文指导/设计

移动网页_全站_页脚广告1

关于我们      便捷服务       自信AI       AI导航        获赠5币

©2010-2024 宁波自信网络信息技术有限公司  版权所有

客服电话:4008-655-100  投诉/维权电话:4009-655-100

gongan.png浙公网安备33021202000488号   

icp.png浙ICP备2021020529号-1  |  浙B2-20240490  

关注我们 :gzh.png    weibo.png    LOFTER.png 

客服