课程设计报告( 2013 -- 2014 年度第二学期)名称：核反应堆热工分析课程设计题目：利用单通道模型进行反应堆稳态热工设计院系：核学院班级：核电1101班学号：1111440113学生：漆圣培指导教师：向斌设计周数：一周成绩：日期：2014 年06 月29 日一、课程设计的目的与要求反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。
对于反应堆热工设计，尤其是对动力堆，最基本的要安全。
要求在整个寿期能够长期稳定运行，并能适应启动、功率调节和停堆等功率变化，要保证在一般事故工况下堆芯不会遭到破坏，甚至在最严重的工况下，也要保证堆芯的放射性物质不扩散到周围环境中去。
在进行反应堆热工设计之前，首先要了解并确定的前提为：（1）根据所设计堆的用途和特殊要求（如尺寸、重量等的限制）选定堆型，确定所用的核燃料、冷却剂、慢化剂和结构材料等的种类；（2）反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化围；（3）燃料元件的形状、它在堆芯的分布方式以及栅距允许变化的围；（4）二回路对一回路冷却剂热工参数的要求；（5）冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。
在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。
目前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：（1）燃料元件芯块最高应低于其他相应燃耗下的熔化温度；（2）燃料元件外表面不允许发生沸腾临界；（3）必须保证正常运行工况下燃料元件和堆构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；（4）在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。
在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而热通道（热管）则是堆芯中轴向功率最高的通道，通过它确定堆芯功率的上限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点来确定DNBR。
热工课程设计主要是为了培养学生综合运用反应堆热工分析课程和其它先修课程的理论和实际知识，树立正确的设计思想，培养分析和解决实际问题的能力。
通过本课程设计，达到以下目的：1、深入理解压水堆热工设计准则；2、深入理解单通道模型的基本概念、基本原理。
包括了平均通道（平均管）、热通道（热管）、热点等在反应堆设计中的应用；3、掌握堆芯焓场的计算并求出体现在反应堆安全性的主要参数：烧毁比DNBR，最小烧毁比MDNBR，燃料元件中心温度及其最高温度，包壳表面温度及其最高温度等；4、求出体现反应堆先进性的主要参数：堆芯流量功率比，堆芯功率密度，燃料元件平均热流密度（热通量），最大热流密度，冷却剂平均流速，冷却剂出口温度等；5、掌握压降的计算；6、掌握单相及沸腾时的传热计算。
7、理解单通道模型的编程方法。
课程设计的考核方式：1、报告一份；2、计算程序及说明一份；3、答辩。
二、设计任务（设计题目）2.2 已知压水反应堆的热功率3500tN Mw=；燃料元件包壳外径9.5csd mm=，包壳径8.6cid mm=，芯块直径8.19ud mm=；燃料组件采用17x17正方形排列，共157组燃料组件；每个组件有24个控制棒套管和一个中子通量测量管；燃料棒中心间栅距P＝13mm,组件间水隙1wmmδ=。
系统工作压力p＝15.51MPa，冷却剂平均温度305Rt C=︒，堆芯冷却剂平均温升27.4t C∆=︒；冷却剂旁流系数6%ζ=；冷却剂设计总流量14314Kg/s，2.55qF=，1.65NHF∆=；DNBR=1.3；又设燃料元件释热份额占总释热量的97.4%；堆芯高度取L＝4.2672 m；并近似认为燃料元件表面最大热流密度、元件表面最高温度和元件中心最高温度都发生在元件半高度处；已知元件包壳的热导率0.00547(1.832)13.8[/()]c csk t W m C=++︒。
试用单通道模型求燃料元件中心温度。
三．设计正文冷却剂有效流量焓升工程热管因子取1.085（未计入交混因子）1）确定燃料元件的实际最大热流密度m axq因为压水堆的安全限值首先是燃料元件表面的最小DNBR，其次才是燃料元件的中心温度，故实际的m axq值由热点处的DNBq值除以DNBR而得。
堆芯燃料棒数目41448)1241717(157=--⨯⨯=N平均热流密度2max/65.155.2646.0mMWFqqq=⨯==--2max/14.23.165.1mMWDNBRqqDNB=⨯==--2）确定燃料元件表面平均热流密度3）堆芯等效直径式中T为正方形组件每边边长（m）。
因为组件无盒壁，组件间水隙，相邻组件的燃料元件棒中心距为，故得236/646.0414482672.4105.9%4.97103500mMWLNdFNqcsut=⨯⨯⨯⨯⨯⨯==--ππ将带入到中，得4)热管版高度处水的比焓堆芯平均管焓升堆芯热管最大焓升假设冷却剂温度变化是线性的，则借助水和蒸汽计算程序，查得当，p=15.51MPa时，故借助水和蒸汽计算程序，在工作压力下5)热管半高处冷却剂流速热管冷却剂流速（或质量流速）的精确计算可按教科书中介绍的方法求解，也可按热管与平均管压降相等的原则进行迭代求解。
作为例子，为简化计算，取热管半高处冷却剂流速近似等于平均管半高处的流速，则式中：为堆芯燃料元件周围的冷却剂总有效流通截面积（）；为冷却剂平均温度下的比容（）；借助水和蒸汽计算程序，可得时，比容。
应由两部分组成：一部分是组件燃料元件棒之间冷却剂的流通面积；另一部分是组件间水隙的横截面积，因为流过这个水隙的冷却剂是冷却燃料组件最外面一排燃料元件的，所以它也属于有效冷却剂的流通面积。
因此有式中为组件间的水隙宽度为0.001mm。
于是6）计算热管半高度处燃料元件表面与冷却剂间的对流换热系数h式中为一个燃料元件栅元中冷却剂通道的当量直径（m），热工参数，在给定的热工参数下，运动粘度，而在给定的热工参数下，水的热导率故7）计算燃料元件表面最高温度式中：为冷却剂的饱和温度（），时，对应饱和温度为由于，故8）燃料元件包壳表面最高温度式中：为包壳厚度（m），为包壳材料导热率，近似认为燃料元件表面最高温度和表面最高温度都发生在元件半高度处因而故得9）燃料芯块表面最高温度式中是包壳与芯块间的气隙等效传热系数，这里取，则于是10）查课本表3-7得查课本表3-7可得四、课程设计总结或结论堆芯部的燃料元件的温度设计终于完成了，经历了很多，也学到了很多。
首先，经过此次的课程设计，掌握了核反应堆堆工设计的基本设计步骤，并在设计中对上课时所学的积分热导率，反应堆和工程热管热点因子等知识从了解到熟练运用；其次学会了office软件中公式编辑器的使用，可以对较难的公式进行编辑，快捷美观；另外，温习了以前学习的MALAB的编程知识，在编程的过程和编程的能力有一定的提升。
尤其在matlab 编程过程中，反复修改程序容，再次体会细心的重要性，感叹编程的不容易。
理论与实践的结合是实践环节的核心，热工作为一名核电学子不但要学，更是要学好的专业课，所以在本次的课程设计，不仅学会了课程的知识，而且也学会了使用自己学过的学习工具应用与本门课程，希望在学习过程中各门课也多把学科之间的学习工具应用于其中，更快更好地学习好相关的课程！！附录（设计流程图、程序、表格、数据等）Nt=input('请输入堆芯功率MW：');d1=input('请输入燃料元件包壳外径mm：');d2=input('请输入燃料元件包壳径mm：');du=input('请输入燃料元件芯块直径mm：');n1=input('请输入燃料组件边长排列数：');n2=input('请输入燃料组件数目：');x1=input('请输入水隙宽度mm：');x2=input('请输入燃料元件栅距mm：');p0=input('请输入压强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*(T^2)/(pi/4));%热管版高度处水的比焓hf=1266.97e3+0.5*(Nt*1e6*fh*1.085)/(w1*(1-w0/100));tl=323.21;at=n2*264*(xp^2-pi/4*(dcs)^2)+n2*4*n1*xp*xw;v=(w1*(1-w0/100))*1.3764e-3/at;%计算热管半高度处燃料元件表面与冷却剂间的对流换热系数De=4*(xp^2-pi/4*dcs)/(pi*dcs);Re=v*De/0.1179e-6;Nu=0.023*Re^0.8*(0.9718)^0.4;hl=Nu*0.514/De;%计算燃料元件表面最高温度ol=q*fq/hl;ojl=345+25*(qmax/1e6)^0.25*exp(-p0/6.2)-tl; if ol<=ojltcs=tl+ol;elsetcs=tl+ojl;end%燃料元件包壳表面最高温度xc=(dcs-dci)/2;kc=0.00547*(1.8*tcs+32)+13.8;oc=q1*fq*xc/(kc*pi*(dcs+dci)/2);tci=tcs+oc;%燃料芯块表面最高温度hg=5678;og=q1*fq/(hg*pi*(dci+dcs)/2);tu=tci+og;%计算燃料芯块的中心最高温度ou=q1*fq/(4*pi)*0.01;t=sym('t');ku=38.24/(t+402.55)+4.788e-13*(t+273.15)^3; fto=double(int(ku,t,0,tu)+ou);q=int(ku);p=q-subs(q,0);to=solve(p-fto);if double(to)>=2200disp('错误');elsedisp('燃料芯块的中心最高温度t(o,max)='); disp(to);end。