华蛀也力*孑课程设计报告（20 13 -- 2014 年度第二学期）名称：核反应堆热工分析课程设计题目：利用单通道模型进行反应堆稳态热工设计院系：_____________________ 班级：实践核1101班 _______________________学号：1111440306 _______________________学生姓名：_____________________指导教师：_____________________设计周数：________________________成绩：________________________日期：2014 年6月19日、课程设计的目的与要求反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。
对于反应堆热工设计，尤其是对动力堆，最基本的要求是安全。
要求在整个寿期内能够长期稳定运行，并能适应启动、功率调节和停堆等功率变化，要保证在一般事故工况下堆芯不会遭到破坏，甚至在最严重的工况下，也要保证堆芯的放射性物质不扩散到周围环境中去。
在进行反应堆热工设计之前，首先要了解并确定的前提为：（1）根据所设计堆的用途和特殊要求（如尺寸、重量等的限制）选定堆型，确定所用的核燃料、冷却剂、慢化剂和结构材料等的种类；（2 ）反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化范围；（3 ）燃料元件的形状、它在堆芯内的分布方式以及栅距允许变化的范围；（4）二回路对一回路冷却剂热工参数的要求；（5 ）冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。
在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。
目前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：（1 ）燃料元件芯块内最高应低于其他相应燃耗下的熔化温度；（2）燃料元件外表面不允许发生沸腾临界；（3）必须保证正常运行工况下燃料元件和堆内构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；（4 ）在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。
在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而热通道（热管）则是堆芯中轴向功率最高的通道，通过它确定堆芯功率的上限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点来确定 DNBR。
热工课程设计主要是为了培养学生综合运用反应堆热工分析课程和其它先修课程的理论和实际知识, 树立正确的设计思想，培养分析和解决实际问题的能力。
通过本课程设计，达到以下目的：1、深入理解压水堆热工设计准则;2、深入理解单通道模型的基本概念、基本原理。
包括了平均通道（平均管）、热通道（热管）、热点等在反应堆设计中的应用；3、掌握堆芯焓场的计算并求出体现在反应堆安全性的主要参数：烧毁比DNBR最小烧毁比 MDNBR 燃料元件中心温度及其最高温度，包壳表面温度及其最高温度等；4、求出体现反应堆先进性的主要参数：堆芯流量功率比，堆芯功率密度，燃料元件平均热流密度（热通量），最大热流密度，冷却剂平均流速，冷却剂出口温度等；5、掌握压降的计算;6、掌握单相及沸腾时的传热计算。
7、理解单通道模型的编程方法。
课程设计要求：1•设计时间为一周；2.独立编制程序计算；3.迭代误差为0.1%;4.计算机绘图；5.设计报告写作认真，条理清楚，页面整洁;6.设计报告中要附源程序。
课程设计的考核方式：1、报告一份；2、计算程序及说明一份；3、答辩。
二、设计任务（设计题目）探求某情况下压水堆核电站对应的热工参数。
某压水反应堆的冷却剂和慢化剂都是水，用二氧化铀作燃料，Zr-4作燃料包壳材料。
燃料组件无盒壁，燃料元件为棒状，正方形排列，已知参数如表一所示：若将堆芯自下而上分为3个控制体，其轴向归一化功率分布见下表：表一堆芯归一化功率分布（轴向等分3个控制体）表一某压水反应堆的热工参数1.堆芯流体出口温度；2•燃料棒表面平均热流密度以及最大热流密度，平均线功率，最大线功率；3.热管内的流体温度（或焓）、包壳表面温度、芯块中心温度随轴向的分布;4.包壳表面最高温度，芯块中心最高温度；5.DNBR在轴向上的变化；6.计算堆芯压降。
三、设计正文（详细的计算过程、计算结果及分析） 1•计算过程1.1堆芯流体出口温度（平均管）1tf ' （tf,intf ,out ）CP 按流体平均温度2以及压力由表中查得。
1.2燃料表面平均热流密度q式中F 总为堆芯燃料棒的总传热面积燃料棒表面最大热流密度q max燃料棒平均线功率q iW/m1.3平均管的情况平均管的流速VV fW( 1-)n o 为燃料组件内正方形排列时的每一排（列）的燃料元件数式中，堆芯内总流通面积A 』.（rv 叽p 2蔦d ； m 4(n °.s)-t f ,outtf ,inFa .Nt W .(1 - ).Cpq = Fa.Nt / F 总W/m 2F 总=mn . d es L .m 2N Eqm ax~q.F q Fqw/m 2q 「d es LV .d es燃料棒最大线功率qi,maxqi, mN Eq i.Fq. Fqw/mm/s—— 1 f二——°由压力以及流体的平均温度t f查表得到：V f1.4为简化计算起见，假定热管内的流体流速 V h 和平均管的V 相同。
降相等来求。
热管内的流体流速要小一些)。
则V h =V 同样，热管四根燃料元件组成的单元通道内的流量(2)第一个控制体出口处的包壳外壁温度if h Z二 t f ,h zu =h⑺De= o.o23Re 0.8 Pr 0.4式中：h(z)可以用k(z)来求。
式中:(实际上，应该按照压：d 2d cs 4 cs1.5热管中的计算(按一个单元通道计算)(1) 热管中的流体温度.q”FR. F 「H F.：H口 …dA b =S 2tf,h Z 厂t f in ,Z。
⑵d Z"t f ,h ZN Eq.F q .F q +)h(z)所以,h(z) =O.O23Re 0.Pr也D eR e=G.D e W h D e[1 D e4( S 2- ■ : 2dCS ) 4:d U(z )和Pr 数根据流体的压力好温度由表查得。
(k=X流体的k(z)、如果流体已经达到过冷沸腾，用 Jen s-Lottes 公式:传热系数)式中：Zr-4 的 k c"00547「8 t c3213.81t c2 c S z ci t z(4)第一个控制体出口处的UO 2芯块外表面温度q.F R NF qE ® (Z) t = t R3_u，hz恥zd . d.h g 2 g(5)第一个控制体出口处的UO 2芯块中心温度 用积分热导求解的方法，即N E ： IRFq . ' Z丿心t dt —k u t dt ------- -/兀其他2个控制体的计算方法相同，重复上述过程即可。
Ff, h( z)△ 9仁h z2t w - t s f, A ef，hZ2时, =25q.F, . z F qE106h zt - ts0.25t f/S!q.F 「¥(Z)F q E'0.25用前面的式子f ,h z2时，用 当(3)第一个控制体出口处的包壳内壁温度t ,c stc j (h )z106■-■-f,hz 2代替宀f,h z二 tcs,h z2 二 k C z一上e®pt s-t f,hzd csd ciW/m. C711.6热管中的qD N B ，hZ用w-3公式计算，同样对3个控制体都算1.7DNBR 的计算1.8计算热管中的压降 1.9单相流体的摩擦压降D e G 2vD e式中:单相流体加速压降:进口:其中，比容V 按相应的流体压力和温度，由表查得。
—2定位格架出口压降:grgrgrVG 2- V outkgrV inDNBRqDNB,h zqDNB,h z「qF N .F q E . zf. r0. 3 16 i Re 0.250 6w < f 丿单相流体提升压降：.g 丄：?in+ P out2.g.L 局部压降，出口:P outkoutV o./ 'out V outoutG 2v °ut2P ink inin Vink inG 2V・iniso2•计算结果1)流体堆芯出口温度t f,out = 344.0019 C；2)燃料棒表面平均热流密度q=8.0505e+005w;3)燃料棒表面最大热流密度q max= 1.7105e+006w ；4)燃料棒平均线功率G= 2.4027e+004w/m ；5)燃料棒最大线功率q i,max = 5.1049e+004w/m；6)热管平均温度t f =316.0009 C；7)第一控制体出口流体温度t f,h(L1) =299.9305 C；8)第一控制体出口处的包壳外壁温t cs,h(L1) =330.1578 C；9)第一控制体出口处的包壳内壁温t ci,h(L1)= 338.8008C；10)第一控制体出口处的芯块外表面温度t u,h( L1) = 517.2810 C；11)第一控制体出口处的芯块中心温度t o,h(L1)= 1.0629e+003C；12)热管中的q DNB,h(L1)=5.5524e+006w ；13)DNBR(L1)=6.200114)第二控制体出口流体温度t f,h(L2) =320.6601 T；15)第二控制体出口处的包壳外壁温t cs,h (L2) =348.6468 C；16)第二控制体出口处的包壳内壁温t ci,h(L2)=364.3162 C；17)第二控制体出口处的芯块外表面温度t u,h( L2) =698.9667 C；18)第二控制体出口处的芯块中心温度t o,h (L2)=2.0864e+003 C；19)热管中的q DNB,h(L2)=5.0732e+006w ；20)DNBR(L2)=3.021321)第三控制体出口流体温度t f,h (L3) =329.4185C；22)第三控制体出口处的包壳外壁温t cs,h (L3) =349.2588 C；23)第三控制体出口处的包壳内壁温t ci,h (L3)=356.6231 C；24)第三控制体出口处的芯块外表面温度t u,h( L3) =512.7934 C；25)第三控制体出口处的芯块中心温度(t o,h L3)=1.0407e+003 C；26)热管中的q DNB,h(L3)=4.3073e+006w ；27)DNBR(L3)=5.496928）单相流体的摩擦压降P f =2.3180e+004Pa29）单相流体加速压降：P a=0Pa30）单相流体提升压降P ei =2.4742e+004Pa31）堆芯出口局部压降:P out =6.9301e+003Pa32）堆芯进口局部压降P n =4.2641e+003Pa33）定位格架出口压降：Fg r =6.6232e+003 Pa34）总的压降P=6.5739e+004Pa3.计算结果分析计算结果误差分析：由于采用的是W-3公式，并且基本没引入对应的修正因子。