压水反应堆的热功率

1、 课程设计报告( 2013 - 2014 年度第 二 学期)名 称： 核反应堆热工分析课程设计 题 目：利用单通道模型进行反应堆稳态热工设计院 系： 核学院 班 级： 核电1101班 学 号： 1111440113 学生姓名： 漆圣培 指导教师： 李向斌 设计周数： 一周 成 绩： 日期： 2014 年 06 月 29 日一、课程设计的目的与要求反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。对于反应堆热工设计，尤其是对动力堆，最基本的要求是安全。要求在整个寿期内能够长期稳定运行，并能适应启动、功率调节和停堆等功率变化，要保证在一般事故工况下堆芯不会遭到破坏，甚至在最严重的工

2、况下，也要保证堆芯的放射性物质不扩散到周围环境中去。在进行反应堆热工设计之前，首先要了解并确定的前提为：（1）根据所设计堆的用途和特殊要求（如尺寸、重量等的限制）选定堆型，确定所用的核燃料、冷却剂、慢化剂和结构材料等的种类；（2）反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化范围；（3）燃料元件的形状、它在堆芯内的分布方式以及栅距允许变化的范围；（4）二回路对一回路冷却剂热工参数的要求；（5）冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。

3、目前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：（1）燃料元件芯块内最高应低于其他相应燃耗下的熔化温度；（2）燃料元件外表面不允许发生沸腾临界；（3）必须保证正常运行工况下燃料元件和堆内构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；（4）在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而热通道（热管）则是堆芯中轴向功率最高的通道，通过它确定堆芯功率的上限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点来确定dnbr。热工课程设计主要是为了培养学生综合运用反应堆热工分析课程

4、和其它先修课程的理论和实际知识，树立正确的设计思想，培养分析和解决实际问题的能力。通过本课程设计，达到以下目的：1、深入理解压水堆热工设计准则；2、深入理解单通道模型的基本概念、基本原理。包括了平均通道（平均管）、热通道（热管）、热点等在反应堆设计中的应用；3、掌握堆芯焓场的计算并求出体现在反应堆安全性的主要参数：烧毁比dnbr，最小烧毁比mdnbr，燃料元件中心温度及其最高温度，包壳表面温度及其最高温度等；4、求出体现反应堆先进性的主要参数：堆芯流量功率比，堆芯功率密度，燃料元件平均热流密度（热通量），最大热流密度，冷却剂平均流速，冷却剂出口温度等；5、掌握压降的计算；6、掌握单相及沸腾时的

5、传热计算。7、理解单通道模型的编程方法。课程设计的考核方式：1、 报告一份；2、计算程序及说明一份；3、答辩。2、 设计任务（设计题目）2.2 已知压水反应堆的热功率；燃料元件包壳外径，包壳内径，芯块直径；燃料组件采用17x17正方形排列，共157组燃料组件；每个组件内有24个控制棒套管和一个中子通量测量管；燃料棒中心间栅距p13mm,组件间水隙。系统工作压力p15.51mpa，冷却剂平均温度，堆芯冷却剂平均温升；冷却剂旁流系数；冷却剂设计总流量14314kg/s，；dnbr=1.3；又设燃料元件内释热份额占总释热量的97.4%；堆芯高度取l4.2672 m；并近似认为燃料元件表面最大热流密度

6、、元件表面最高温度和元件中心最高温度都发生在元件半高度处；已知元件包壳的热导率。试用单通道模型求燃料元件中心温度。三设计正文冷却剂有效流量wef=wt1-=143141-6%=13455kg/s焓升工程热管因子fhe取1.085（未计入交混因子）1）确定燃料元件的实际最大热流密度因为压水堆的安全限值首先是燃料元件表面的最小dnbr，其次才是燃料元件的中心温度，故实际的值由热点处的值除以dnbr而得。堆芯燃料棒数目 平均热流密度 2）确定燃料元件表面平均热流密度qq=0.646mw/m2q1=qdcs=0.6463.140.0095=0.0193mw/mq1,max=q1fq=0.01932.5

7、5=0.0491mw/m3）堆芯等效直径defdef=121t2/4式中t为正方形组件每边边长（m）。因为组件无盒壁，组件间水隙w=1mm，相邻组件的燃料元件棒中心距为13+1=14mm，故得t2=(171310-3+20.510-3)2=0.049m2将t2带入到def中，得def=2.75m4)热管版高度处水的比焓hl2hl2=hf,in+hmax2堆芯平均管焓升h=ntwef=350010613455=260kj/kg堆芯热管最大焓升hh,max=hfhnfhe=2601031.651.085=465.5kj/kg假设冷却剂温度变化是线性的，则tf,in=286.7借助水和蒸汽计算程序，

8、查得当tf,in=286.7，p=15.51mpa时，hf,in=1266.96 kj/kg故hl2=hf,in+hmax2=1266.96+465.52=1499.71kj/kg借助水和蒸汽计算程序，在工作压力下tl2=327.39 5)热管半高处冷却剂流速热管内冷却剂流速（或质量流速）的精确计算可按教科书中介绍的方法求解，也可按热管与平均管压降相等的原则进行迭代求解。作为例子，为简化计算，取热管半高处冷却剂流速近似等于平均管半高处的流速，则v=wefvfat式中：at为堆芯燃料元件周围的冷却剂总有效流通截面积（m2）；vf为冷却剂平均温度下的比容（m3/kg）；借助水和蒸汽计算程序，可得t

9、r=305，p=15.51mpa时，比容vf=1.396710-3 m3/kg。at应由两部分组成：一部分是组件内燃料元件棒之间冷却剂的流通面积；另一部分是组件间水隙的横截面积，因为流过这个水隙的冷却剂是冷却燃料组件最外面一排燃料元件的，所以它也属于有效冷却剂的流通面积。因此有at=np2-4dcs2+n4(17p)式中为组件间的水隙宽度为0.001mm。at=1572641310-32-40.00952+1574171310-3110-3=4.20m2于是v=wefvfat=134551.396710-34.20=4.4744m/s6）计算热管半高度处燃料元件表面与冷却剂间的对流换热系数hr

10、el2=vdev式中de为一个燃料元件栅元中冷却剂通道的当量直径（m），de=4abal=4p2-4dcs2dcs=41310-32-40.009520.0095=0.0132m热工参数p=15.51mpa， tfl2=327.39 在给定的热工参数下，运动粘度=0.120910-6m2/s，pr=0.8751 rel2=vde=4.47440.01320.120910-6=4.88105而nul2=0.023re0.8pr0.4=0.0234.881050.80.87510.4=774.97在给定的热工参数下，水的热导率=552.810-3w/(m)故hl2=nude=0.5528774.9

11、70.0132=32.45103w/m27）计算燃料元件表面最高温度tcs,maxtcs,max=tfl2+fl2fl2=fl2 当fl2f,jl2f,jl2 当f,jl2fl2fl2=qfqhl2=0.6461062.5532.45103=50.76f,jl2=ts+25qmax1060.25exp-p6.2-tfl2式中：ts为冷却剂的饱和温度（），p=15.51mpa时，对应饱和温度为344.84f,jl2=344.84+251.651061060.25exp-15.516.2-327.39=19.77由于f,jl2fl2，故tcs,max=tfl2+f,jl2=327.39+19.77

12、=347.168）燃料元件包壳内表面最高温度tci,maxtci,max=tcs,max+cl2cl2=q1fqccdcs+dci2式中：c为包壳厚度（m），c=dcs-dci2=9.5-8.62=0.4510-3mc为包壳材料导热率w/m，近似认为燃料元件表面最高温度和内表面最高温度都发生在元件半高度处c=0.005471.8tcs,max+32+13.8=0.005471.8347.16+32+13.8=17.39w/m因而c=0.01931062.550.4510-317.393.149.5+8.6210-3=44.81故得tci,max=347.16+44.81=391.979）燃料芯

13、块表面最高温度tu,maxtu,max=tci,max+gl2gl2=q1fqhgdci+du2式中hg是包壳与芯块间的气隙等效传热系数，这里取hg=5678w/m2，则g=19.31032.5556783.148.6+8.19210-3=328.8于是tu,max=391.97+328.8=720.7710）计算燃料芯块的中心最高温度to,max0to,maxudt=0tu,maxudt+q1fq4=0916.43udt+19.31032.5543.1410-2查课本表3-7得0673.38udt=45.140to,maxudt=38.65+39.18=77.83查课本表3-7可得to,ma

14、x=21552200四、课程设计总结或结论堆芯内部的燃料元件的温度设计终于完成了，经历了很多，也学到了很多。首先，经过此次的课程设计，掌握了核反应堆堆工设计的基本设计步骤，并在设计中对上课时所学的积分热导率，反应堆和工程热管热点因子等知识从了解到熟练运用；其次学会了office软件中公式编辑器的使用，可以对较难的公式进行编辑，快捷美观；另外，温习了以前学习的malab的编程知识，在编程的过程和编程的能力有一定的提升。尤其在matlab编程过程中，反复修改程序内容，再次体会细心的重要性，感叹编程的不容易。理论与实践的结合是实践环节的核心，热工作为一名核电学子不但要学，更是要学好的专业课，所以在本

15、次的课程设计，不仅学会了课程的知识，而且也学会了使用自己学过的学习工具应用与本门课程，希望在学习过程中各门课也多把学科之间的学习工具应用于其中，更快更好地学习好相关的课程！附录（设计流程图、程序、表格、数据等）nt=input(请输入堆芯功率mw：);d1=input(请输入燃料元件包壳外径mm：);d2=input(请输入燃料元件包壳内径mm：);du=input(请输入燃料元件芯块直径mm：);n1=input(请输入燃料组件边长排列数：);n2=input(请输入燃料组件数目：);x1=input(请输入水隙宽度mm：);x2=input(请输入燃料元件栅距mm：);p0=input(请

16、输入压强mpa：);t0=input(请输入冷却剂平均温度：);dt=input(请输入堆芯冷却剂温升：);w0=input(请输入冷却剂旁流系数（%）:);w1=input(请输入堆芯体积流量kg/s：);fq=input(请输入热流密度热点因子：);fh=input(请输入焓升热管因子：);fu=input(燃料元件内的释热量占堆总热功率的份额(%)：);dnbr=input(请输入最小热流密度(dnbr)：);l=input(请输入燃料原件高度m：);dcs=d1/1000;dci=d2/1000;xw=x1/1000;xp=x2/1000;%确定燃料元件的实际最大热流密度n=157*(17*17-25)q= nt*fu/(3.14*dcs*l*n)qmax=q*fq%确定燃料元件表面平均热流密度q1=q*pi*dcs;q1max=q1*fq;%堆芯等效直径t=n1*xp+xw;def=sqrt(121*(t2)/(pi/4);%热管版高度处水的比焓hf=1266.97e3+0.5*(nt*1e6*fh*1.085)/(w1*(1-w0/100);tl=323.21;at=n2*264*(xp2-pi/4*(dcs)2)+n2*4*n1*xp*xw;v=(w1*(1-w0