核反应堆热工分析课程设计报告书详细过程版本

课程设计报告(2013--2014年度第二学期)名称：核反应堆热工分析课程设计题目：利用单通道模型进行反应堆稳态热工设计院系：核科学与工程学院班级：实践核1101班学号：1111440306学生姓名：蒋佳指导教师：王胜飞设计周数：1周成绩：日期：2014年6月19日课程热工课程设计报告1一、课程设计的目的与要求反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。对于反应堆热工设计，尤其是对动力堆，最基本的要求是安全。要求在整个寿期内能够长期稳定运行，并能适应启动、功率调节和停堆等功率变化，要保证在一般事故工况下堆芯不会遭到破坏，甚至在最严重的工况下，也要保证堆芯的放射性物质不扩散到周围环境中去。在进行反应堆热工设计之前，首先要了解并确定的前提为：（1）根据所设计堆的用途和特殊要求（如尺寸、重量等的限制）选定堆型，确定所用的核燃料、冷却剂、慢化剂和结构材料等的种类；（2）反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化范围；（3）燃料元件的形状、它在堆芯内的分布方式以及栅距允许变化的范围；（4）二回路对一回路冷却剂热工参数的要求；（5）冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。目前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：（1）燃料元件芯块内最高应低于其他相应燃耗下的熔化温度；（2）燃料元件外表面不允许发生沸腾临界；（3）必须保证正常运行工况下燃料元件和堆内构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；（4）在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而热通道（热管）则是堆芯中轴向功率最高的通道，通过它确定堆芯功率的上限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点来确定DNBR。热工课程设计主要是为了培养学生综合运用反应堆热工分析课程和其它先修课程的理论和实际知识，树立正确的设计思想，培养分析和解决实际问题的能力。通过本课程设计，达到以下目的：1、深入理解压水堆热工设计准则；2、深入理解单通道模型的基本概念、基本原理。包括了平均通道（平均管）、热通道（热管）、热点等在反应堆设计中的应用；3、掌握堆芯焓场的计算并求出体现在反应堆安全性的主要参数：烧毁比DNBR，最小烧毁比MDNBR，燃料元件中心温度及其最高温度，包壳表面温度及其最高温度等；4、求出体现反应堆先进性的主要参数：堆芯流量功率比，堆芯功率密度，燃料元件平均热流密度（热通量），最大热流密度，冷却剂平均流速，冷却剂出口温度等；5、掌握压降的计算；课程热工课程设计报告26、掌握单相及沸腾时的传热计算。7、理解单通道模型的编程方法。课程设计要求：1．设计时间为一周；2．独立编制程序计算；3．迭代误差为0.1%；4．计算机绘图；5．设计报告写作认真，条理清楚，页面整洁；6．设计报告中要附源程序。课程设计的考核方式：1、报告一份；2、计算程序及说明一份；3、答辩。二、设计任务（设计题目）探求某情况下压水堆核电站对应的热工参数。某压水反应堆的冷却剂和慢化剂都是水，用二氧化铀作燃料，Zr-4作燃料包壳材料。燃料组件无盒壁，燃料元件为棒状，正方形排列，已知参数如表一所示：若将堆芯自下而上分为3个控制体，其轴向归一化功率分布见下表：表一堆芯归一化功率分布（轴向等分3个控制体）自下而上控制体123归一化功率分布()z0.801.500.70表一某压水反应堆的热工参数参数名称参考值参数范围说明系统压力P15.5MPa14.5-16堆芯输出热功率Nt2895MW4000进、出口温度、流量、功率耦合，由3个量可推出另外一个冷却剂总流量W32100t/h*3反应堆进口温度fin288℃265-295反应堆出口温度fout310-330堆芯高度L3.66m3.5-4燃料组件数m121根据功率计算燃料组件形式n0xn017x1782,152,172等需要除去控制棒和中子通量管每个组件燃料棒数n265燃料包壳外径dcs9.5mm9-11mm三者耦合燃料包壳内径dci8.60mm8-9mm燃料包壳厚度δc0.57mm0.5-0.9mm课程热工课程设计报告3燃料芯块直径du8.19mm8-9mm燃料棒间距（栅距）s12.6mm12-14mm两个组件间的水隙δ0.8mm0.5-1mmUO2芯块密度2UO95%理论密度90-99%旁流系数5%4-7%燃料元件发热占总发热的份额Fa97.4%96-98%径向热管因子NRF1.351.2-1.6核热点因子由前3者乘积可以推出轴向热管因子NZF1.5281.5-2.0局部峰核热管因子NLF1.11热流量核热点因子NNNNqRZLFFFF2.29热流量工程热点因子EqF1.031.03-1.09焓升工程热管因子EHF（未计入交混因子）1.0851.05-1.16交混因子*EHmF0.95焓升核热管因子NNHRFF1.351.35-1.8堆芯入口局部阻力系数Kin0.75堆芯出口局部阻力系数Kout1.0堆芯定位隔架局部阻力系数Kgr1.05通过计算，得出：1.堆芯流体出口温度；2.燃料棒表面平均热流密度以及最大热流密度，平均线功率，最大线功率；3.热管内的流体温度（或焓）、包壳表面温度、芯块中心温度随轴向的分布；4.包壳表面最高温度，芯块中心最高温度；5.DNBR在轴向上的变化；6.计算堆芯压降。课程热工课程设计报告4三、设计正文（详细的计算过程、计算结果及分析）1.计算过程1.1堆芯流体出口温度（平均管）,,..(1).foutfinFaNtttWCp℃Cp按流体平均温度,,1()2ffinfoutttt以及压力由表中查得。1.2燃料表面平均热流密度q./qFaNtF总W/m2式中F总为堆芯燃料棒的总传热面积=....csFmndL总m2燃料棒表面最大热流密度qmaxmaxqq..NEqqFFw/m2燃料棒平均线功率lq.....cslcsqdLqqdLW/m燃料棒最大线功率,maxlq,max..NEllqqqqFFw/m1.3平均管的情况平均管的流速V(1)..ffefffWWVAAm/s式中，堆芯内总流通面积22000.(.)4(.)42fcsAmnnpdmnsn0为燃料组件内正方形排列时的每一排（列）的燃料元件数课程热工课程设计报告5f由压力以及流体的平均温度ft查表得到：1ffV1.4为简化计算起见，假定热管内的流体流速Vh和平均管的V相同。（实际上，应该按照压降相等来求。热管内的流体流速要小一些）。则Vh=V同样，热管四根燃料元件组成的单元通道内的流量(1)hhfWWAA224bcsASd1.5热管中的计算（按一个单元通道计算）（1）热管中的流体温度.,(),0..(.)()NEEzRHHmcsfhzfinzhqFFFdttzdWCp（2）第一个控制体出口处的包壳外壁温度1,(),(),()cshzfhzfhztt,()..()NEqqfhzqFFthz,()..().()NERqfhzqFzFthz式中：h(z)可以用0.80.4()0.023RePr()ehzDNukz来求。所以，0.80.4()()0.023RePr.ekzhzD式中：.Re.ehebGDWDA224()44csbecsSdADUd流体的k(z)、μ（z）和Pr数根据流体的压力好温度由表查得。（k=λ传热系数）如果流体已经达到过冷沸腾，用Jens-Lottes公式：课程热工课程设计报告60.256.26..().25.10pNERqwsqFzFtte2,(),(),fhzwfhzwssfhztzttzttt0.256.2,6..25.10pNERqsfhzqFZFett当22,,fhzfhz时，用前面的式子当22,,fhzfhz时，用21,,fhzfhz代替（3）第一个控制体出口处的包壳内壁温度,,,tcihzcshzchzt,...ln2NERqcscshzCciqFFzdtkzd式中：Zr-4的0.005471.83213.8ccktW/m.℃12ccsciztztzt(4)第一个控制体出口处的UO2芯块外表面温度,,..2NERquhzcihzciugqFFZttddh（5）第一个控制体出口处的UO2芯块中心温度用积分热导求解的方法，即0,,00..4hzuhzNEttRquuqFFZktdtktdt其他2个控制体的计算方法相同，重复上述过程即可。课程热工课程设计报告71.6热管中的,DNBhzq用w-3公式计算，同样对3个控制体都算1.7DNBR的计算,,...DNBhzDNBhzNERqqqDNBRqzqFFz1.8计算热管中的压降1.9单相流体的摩擦压降22..2..2feeLVpfDLGvfD式中：0.640.250.316..Renwwisoffff单相流体加速压降：单相流体提升压降：....2inoutelpgLgL局部压降，出口：2222outoutoutoutoutoutVGvpkk进口：2222ininininininVGvpkk定位格架出口压降：221222outingrgrgrgrGvvVpkk其中，比容v按相应的流体压力和温度，由表查得。课程热工课程设计报告82.计算结果1)流体堆芯出口温度outft,=344.0019℃；2)燃料棒表面平均热流密度q=8.0505e+005w；3)燃料棒表面最大热流密度maxq=1.7105e+006w；4)燃料棒平均线功率lq=2.4027e+004w/m；5)燃料棒最大线功率max,lq=5.1049e+004w/m；6)热管平均温度ft=316.0009℃；7)第一控制体出口流体温度hft,（L1）=299.9305℃；8)第一控制体出口处的包壳外壁温hcst,（L1）=330.1578℃；9)第一控制体出口处的包壳内壁温hcit,(L1)=338.8008℃；10)第一控制体出口处的芯块外表面温度hut,（L1）=517.2810℃；11)第一控制体出口处的芯块中心温度hot,(L1)=1.0629e+003℃；12)热管中的hDNBq,(L1)=5.5524e+006w；13)DNBR(L1)=6.200114)第二控制体出口流体温度hft,（L2）=320.6601℃；15)第二控制体出口处的包壳外壁温hcst,（L2）=348.6468℃；16)第二控制体出口处的包壳内壁温hcit,(L2)=364.3162℃；17)第二控制体出口处的芯块外表面温度hut,（L2）=698.9667℃；18)第二控制体出口处的芯块中心温度hot,(L2)=2.0864e+003℃；19)热管中的hDNBq,(L2)=5.0732e+006w；20)DNBR(L2)=3.021321)第三控制体出口流体温度hft,（L3）=329.4185℃；22)第三控制体出口处的包壳外壁温hcst,（L3）=349.2588℃；23)第三控制体出口处的包壳内壁温hcit,(L3)=356.6231℃；24)第三控制体出口处的芯块外表面温度hut,（L3）=512.7934℃；25)第三控制体出口处的芯块中心温度(ho