核反应堆热工水力课程设计.docx

1、核反应堆热工水力课程设计一、设计要求在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预 先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。目 前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：1.燃料元件芯块内最高应低丁其他相应燃耗下的熔化温度；2.燃料元件外表面不允许发生沸腾临界；3.必须保证正常运行工况下燃料元件和堆内构件得到充分冷却； 在事故工况下能提供足够的冷却剂以排除堆芯余热；4.在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。5.在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而 热通道（热管）则是堆芯中轴

2、向功率最高的通道，通过它确定堆芯功率的上 限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点 来确定DNBR二、设计任务某压水反应堆的冷却剂和慢化剂都是水， 用二氧化铀作燃料，Zr-4作燃料包 壳材料。燃料组件无盒壁，燃料元件为棒状，正方形排列，已知下列参数：系统压力P 15.8 MPa堆芯输出热功率Nt 1 820MW冷却剂总流量W 32500t / h反应堆进口温度t f i n 287 C堆芯高度L 3. 60m燃料组件数m 1 2 1燃料组件形式n0x n0 17 x 1 7每个组件燃料棒数n 265燃料包壳外径dcs 9. 5mm燃料包壳内径dci 8. 6mm燃料包

3、壳厚度a c 0. 57mm燃料芯块直径du 8. 19mm燃料棒间距（栅距）s 12.6mm两个组件间的水隙a 0. 8 mmUO犯块密度p UO2 95% 理论密度旁流系数1 5%燃料兀件发热占总发热份额Fa 97. 4%径向核热管因子 1. 33轴向核热管因子 1. 5201. 03热流量核热点因子 = 2. 022热流量工程热点因子焰升工程热点因子（未计入交混因子）1. 142交混因子水含升核热管因子 =堆芯进口局部阻力系数Ki n 堆芯出口局部阻力系数Kout堆芯定位格架阻力系数Kgr0. 951. 0850. 751. 01. 05若将堆芯自下而上分为3个控制体，其轴向归一化功率分

4、布见下表:自卜而上控制体号123456归一化功率分布0.481.0 21.5 01.560.960.4 8表 堆芯归一化功率分布（轴向等分 3个控制体）通过计算，得出：1.堆芯流体出口温度；2.燃料棒表面平均热流密度以及最大热流密度，平均线功率，最大线功率；3.热管内的流体温度（或焰）、包壳表面温度、芯块中心温度随轴向的分布;4.包壳表面最高温度，芯块中心最高温度；5.DNBR在轴向上的变化；6.计算堆芯压降三、设计正文（详细的计算过程、计算结果及分析）1.计算过程1.1堆芯流体出口温度（平均管）按流体平均温度以及压力由表中查得。假设 C,查表得经过输入所查程序不断迭代得如需更精确的值，可以继

5、续进行迭代计算1.2燃料表面平均热流密度CC误差小丁 0.5 C总式中总为堆芯燃料棒的总传热面积代入数据得燃料棒表面最大热流密度qmax代入数据得燃料棒平均线功率代入数据得燃料棒最大线功率代入数据得1.3平均管的情况平均管的流速V式中 堆芯内总流通面积为燃料组件内正方形排列时的每一排（列）的燃料元件数由压力以及流体的平均温度 查表得到：由1.1知 C ,查表得1.4为简化计算起见，假定热管内的流体流速 M和平均管的V相同 同样，热管四根燃料元件组成的单元通道内的流量代入数据得1.5热管中的计算(按一个单元通道计算)(1)热管中的流体温度 其中 取平均温度对应的参数值，需要进行迭代计算，下面给出

6、第一控制体出口处温度的算法假设- C，查表得 C，带入上式- C与假设误差较大，进行迭代，查表知 C- C误差 C,可以不再进行迭代，就取 - C同理由程序迭代可求得第二控制体出口处流体温度 一 C第三控制体出口处流体温度 一 C第四控制体出口处流体温度 一 C第五控制体出口处流体温度 一 C第六控制体出口处流体温度 C(2)第一个控制体出口处的包壳外壁温度式中：h(z)为单相水强迫对流换热系数 C所以式中流体的入、和Pr数根据流体的压力和温度由表查得。如果流体已经达到过冷沸腾，用 Jens-Lottes公式： 其中为气体的饱和温度，p的单位为MPa当 时，用前面的式子当 时，用 替换掉,可以

7、利用以下公式来求时， C代入数据得C 4 - CCCC 4 一第一控制体出口处-C,查表可得入 - -故- - A其余同理由程序计算得出结果如下 第二控制体出口处-C,查表可得入 C故一 一 一 C第三控制体出口处11CC故一 C第四控制体出口处-C,查表可得入 C 41CC故一 C第五控制体出口处C ,查表可得入 C 41CC故 A C第六控制体出口处C,查表可得入 C 41C C故 c(3)包壳内壁温度式中Zr-4的代入数据得:由丁 与平均温度有关，由程序迭代计算结果如下第一控制体出口处 一 C第控制体出口处 一 C第三控制体出口处 一 C第四控制体出口处 一 C第五控制体出口处 一 C第

8、六控制体出口处 C(4)燃料芯块外表面温度式中？是包壳与芯块间的气隙等效传热系数，这里取 ？ C代入数据得第一个控制体出口处第二个控制体出口处一 C第三个控制体出口处 C第四个控制体出口处 C第五个控制体出口处 C第六个控制体出口处C(5)燃料芯块中心温度用积分热导求解的方法，即其中代入数据得令由丁函数递增，可以通过二分法求解 f(x)的根得出二氧化铀中心温度通过编程可求得结果如下第一控制体出口处第二控制体出口处-C第三控制体出口处-C第四控制体出口处-C第五控制体出口处-C第六控制体出口处C1.6热管中的用w-3公式计算，同样对3个控制体都算式中：p为冷却剂工作压力（Pa）, G为冷却剂质量

9、流密度 ？，为冷却剂通道的当量直径（, ? 为冷却剂的饱和比焰（J/kg ） , ? 为控制体进口处冷却剂的比焰（J/kg ）,为计算点z处的平衡含气量， 为其绝对值。平衡含气量的计算式为其中？ 为汽化潜热（J/kg ）通过程序计算得出结果如下第一控制体出口处第二控制体出口处第三控制体出口处第四控制体出口处第五控制体出口处第六控制体出口处1.7DNBR勺计算第一控制体出口处同理可求其余段DNBR结果如下: 第二控制体出口处第三控制体出口处第四控制体出口处第五控制体出口处第六控制体出口处1.8计算热管中的压降单相流体的摩擦压降式中：使用公式编程分别计算六段控制体的摩擦压降 单相流体加速压降：同样

10、使用公式编程分别计算六段控制体的加速压降单相流体提升压降同上局部压降，出口：代入数据得局部压降，进口:代入数据得局部压降，定位格架出口压降代入数据得以上所使用的比热容和动力粘度都通过软件查询后输入代码中进行计算。 通过程序计算结果如下分段压降变化控制体123456提升压 降/ pa4.3974.3164.1733.9843.8133.708摩擦压 降/ pa1.0030.9350.4500.4450.4390.435加速压 降/ pa0.1170.2690.4530.5650.4120.228总压降为温度汇总表控制段控制体外流体出口温度/ C包壳外表面温度C包壳内表面温度芯块表面温度芯块中心温

11、度1291.69303.51313.13384.03568.792301.3326.15346.26492.83974.833314.86348.41377.54593.091383.574327.70348.43378.72602.891430.385334.89348.19366.89504.84958.776338.22347.90357.27426.25619.90临界热流密度和烧毁比汇总控制体123456临界热流密度5.2444.7313.8973.1092.6242.388烧毁比15.076.403.5852.7503.7226.8647单从表内数据看，DNBFW大丁 1,芯块中心

12、温度小丁二氧化铀的熔点 2878C, 所以理论上能够保证安全性。四、课程设计感想通过这次反应堆热工分析的课程设计，我加深了对反应堆内部传热的了解， 同时我发发现了自己的很多不足之处。拿到课程设计题目后，我首先根据题目中的问题去复习热工分析的知识， 去了解每个问题该怎么算。然后将计算思路和公式写出来，为下一步的编程做准备。 由丁大部分的求解都是非线性的方程，所以采用了迭代和二分法求解方程。当把整个问题的求解思路理活和确定计算过程中使用算法后， 整个问题就在水和水蒸气热物性如何导入上了。 在最开始，我使用通过一定的数据去拟合用到 的热物性在定压下与温度的关系。 但在拟合之后,通过得出的结果与接用水

13、和水 蒸气热物性查询软件得出的结果有一点差距，所以就放弃了。然后，我想将网上通过水和蒸汽热力性质公式 IAPWS_IF9m算的源码做成 接口，在计算过程中直接用它计算热物性。后来，发现太难，i(码有很多地方看 不懂，不知道如何接入。丁是最终也放弃了。后来就采用了最原始的方法，就是 在计算的过程中，通过水和水蒸气热物性查询软件查出结果，一步一步输进去。 最后使用程序计算时，在输入参数时耗时太久，我深深的感受到书到用处方恨少， 如果原来，把编程多学一点，将整个数据导入直接通过动态链接库完成， 就省事 太多。这次课设让我把原来学的编程与实际问题结合起来， 给我很大感触，原来敲着书上的例子，感觉漫无目

14、的。通过这次自己去思考如何解决一个问题， 让我理 解了拿到一个问题后，该如何处理。同时通过这次课设，我意识到反应堆热工分 析是个复杂而连续的过程，每个参数都受到大量的常数参数的影响也具有很多的 修正因子。我们不应该根据自己的常识来判断数据的变化情况， 相反地我们应该随时坚持以数据计算为引导，以实验作为验证。仔细客观认真地分析堆内数据的 变化，并且对堆内数据进行全程监控，防止堆内数据随时变化，对反应堆的危害 性。附录(设计流程图、程序)1、程序说明共7个程序由c语言编写，分别为计算流体出口温度，控制体出口流体温度,燃料包壳外壁温度，燃料包壳内壁温度，二氧化铀中心温度， qDNB压降。使用了迭代和

15、二分法求解部分非线性方程。(1)堆芯出口温度计算：此段根据任务书给出的基本参数和热量与流量之间关系，运用迭代的算法，求出堆芯的出口温度。(2)第一至第六控制体的各量计算：因为六个控制体的计算过程类似，这里只说明第一个控制体的计算过程。在现有的参数下，根据热流量与流量的关系和迭代算法，求出该控制体的出 口温度。通过流通截面积与湿周的关系求出栅元的当量直径。 再根据上面的温度，查出对应的热物性参数由雷诺数与努尔数的关系，解出控制体出口处的对流 换热系数。因为不知该处的流体状态，分别用单相强迫对流放热公式和詹斯-洛特斯 传热方程算出各自的膜温压，取较小的值加上出口处的流体温度即是包壳的 外表面温度。

16、由包壳的外表面的温度再根据圆管的传热方程运用迭代算法解出包壳内 表面的温度。芯块与包壳内表面之间的导热问题，根据间隙导热模型，即可 解出芯块表面的温度，根据内热源的导热模型，依据积分热导率与温度的对 应关系列出方程用二分法解出芯块中心的温度。接下来依据冷却剂的温度， 得出的控制体出口处的含汽量。进而依据 W-3公式求出该出的临界热流量 qDNB最后得出该出的烧毁比DNBR(3)热管的压降计算：热管的压降包括摩擦压降、提升压降、进出口局部压降、定位搁架出口压降。摩擦压降可由计算单相流的达西(Darcy)公式算得。提升压降可由根据位置 的变化算得，其中参数都取平均值。其余的压降根据形阻压降的基本公

17、式再 乘以相应的系数求得。最后各项相加得出热管的总压降。2、流程图2、程序/ 3.1.cpp : 定义控制台应用程序的入口点/流体出口温度迭代计算#include stdafx.h#include #include using namespacestd;double tfout( double cp, double atfout)double tfin = 287, Fa = 0.974,Nt = 1820,W=32500,plxs=0.05,tfout;tfout = tfin + 3600*Fa*Nt / (W*cp*(1 - plxs); return tfout;int main()d

18、ouble atfout,tfin=287,tav,cp,tfoutc=0;int panding = 1;printf(请输入假设出口温度n);cin atfout; |while (panding=1)tav = (atfout + tfin) / 2;printf(请输入 讣的 Cpn , tav);cin cp;tfoutc = tfout(cp, atfout);if (fabs(tfoutc - atfout) = 0.5)panding = 0;atfout = tfoutc;printf(出口温度为 %fn ,tfoutc);while ( true ) /3.3.1热管内流体

19、温度计算#include stdafx.h#include #include using namespacestd;double ffout( double tfin, double atfout)double qav, fnr=1.33, feh=1.03, fehm=0.95, dcs=9.5, gyh, cp, wh,l=3.60,tfout;bool panding = true ;cout 请输入平均热流密度( W/m*m qav;cout 请输入归一化参数 ” gyh;cout 请输入热管单元通道流量 Wh( t/h ) wh;while (panding)double tav =

20、 (atfout + tfin) / 2;cout 请输入 tav 下的 Cp(J/(kg. C) cp;tfout = tfin + 3.6*(qav*fnr*feh*fehm*3.14*dcs*gyh*l / 6) / (wh*cp*1000);if (fabs(tfout - atfout) = 0.5)panding = false ;atfout = tfout;)return tfout;)int main()(double tfin,atfout,tfout;cout 请输入流体进口温度” tfin;cout 请输入流体假设流体出口温度 atfout;tfout = ffout(

21、tfin, atfout);cout tfout;while ( true )()/ 3.3.02.cpp : 定义控制台应用程序的入口点/计算包壳外壁温度#include stdafx.h#include cmath#include iostreamusing namespacestd;int main()(double tfout6 = ( 291.69,301.38, 314.86, 327.70, 334.89, 338.22 ), gyhcs6=0.48,1.02,1.50,1.56,0.96,0.48;double pr6,drxs6,dlnd6,f16,f26,f6;int i;

22、for (i = 0; i 6; i+)cout 请输入 tfouti 下的普朗特数 pri;cout 请输入 tfouti 下的导热系数(W/(m? C) drxsi;cout ”请输入 tfouti 下的动力粘度 kg/(m ? s) dlndi;f1i = 2.25*pow(10, 4)*gyhcsi * pow(dlndi, 0.8)/(drxsi*pow(pri,0.4);f2i = 346.38 + 1.83*pow(gyhcsi, 0.25) - tfouti;if (f1if2i)fi = f2i;elsefi = f1i;cout e f1= f1i e f2= f2i f=

23、 fi + tfouti endl;)while ( true )( )/ 3.3.2.cpp : 定义控制台应用程序的入口点/计算内壁温度#include stdafx.h#includeiostream#includemath.husing namespacestd;double ftci( double tcs, double atci)(double frn = 1.33, feq = 1.03, gyhcs, dcs=9.5, dci=8.6, kc, tav,q,tci;bool panding = true ;cout 请输入归一化参数 gyhcs;cout 请输入平均线功率，单

24、位 w/m q;while (panding)(tav = (tcs + atci) / 2;kc = 0.00547*(1.8*tav + 32) + 13.8;tci = tcs + log(dcs/dci)*(q*frn*feq*gyhcs) / (2 * 3.14*kc);if (tci - atci = 0.5)panding = false ;atci = tci;return tci;int main()(double tcs, tci,atci;cout 请输入夕卜壁温度C tcs;cout 请输入内壁假设温度C atci;tci = ftci(tcs, atci);cout

25、tci;while ( true )(/ 3.3.3.cpp : 定义控制台应用程序的入口点/计算二氧化优中心温度#include stdafx.h#include iostream#include math.h#include using namespacestd;double f1( double t)double fout;fout = 38.24*log(t+402.55)/402.55)+4.788*pow(10, -13)*(pow(t + 273.15), 4)- pow(273.15, 4);return fout;double f2( double tci, double a

26、tu)double tu, fnr = 1.33, feq = 1.03, gyhcs, q, tmid, jieguo;bool panding = true ;cout 请输入归一化参数 gyhcs;cout 请输入线热流密度/w q;if (f1(atu) - f1(tci) - (q*fnr*feq*gyhcs) / (4 * 314) 0) tu = tci;panding = false ;while (panding)if (f1(atu) - f1(tci) - (q*fnr*feq*gyhcs) / (4 * 314) 0)panding = false ;tu = atu - 100;elseatu = atu + 100;while (atu-tu0.0000001) tmid = (atu + tu) / 2;jieguo = f1(tmid) - f1(tci) - (q*fnr*feq*gyhcs) / (4 * 314); if (jieguo = 0)tu = tmid;elseif (jieguo 0)atu = tmid;elsetu = tmid;return tu;int main()double tci, atu, tu;cout 请输入燃料芯块外壁温度/C tci;cout 请输入假设二氧化铀中心温度/C atu;tu=