核反应堆热工水力课程设计报告书.doc
一、设计要求在设计反应堆冷却系统时.为了保证反应堆运行安全可靠.针对不同的堆型.预先规定了热工设计必须遵守的要求.这些要求通常就称为堆的热工设计准则。目前压水动力堆设计中所规定的稳态热工设计准则.一般有以下几点：1. 燃料元件芯块内最高应低于其他相应燃耗下的熔化温度；2. 燃料元件外表面不允许发生沸腾临界；3. 必须保证正常运行工况下燃料元件和堆内构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；4. 在稳态额定工况和可预计的瞬态运行工况中.不发生流动不稳定性。5. 在热工设计中.通常是通过平均通道平均管可以估算堆芯的总功率.而热通道热管则是堆芯中轴向功率最高的通道.通过它确定堆芯功率的上限.热点是堆芯中温度最高的点.代表堆芯热量密度最大的点.通过这个点来确定DNBR。二、设计任务某压水反应堆的冷却剂和慢化剂都是水.用二氧化铀作燃料.Zr-4作燃料包壳材料。燃料组件无盒壁.燃料元件为棒状.正方形排列.已知下列参数：系统压力P 15.8MPa堆芯输出热功率Nt 1820MW冷却剂总流量W 32500t/h反应堆进口温度tfin 287堆芯高度L3.60m燃料组件数m121燃料组件形式n0×n0 17×17每个组件燃料棒数n 265燃料包壳外径dcs 9.5mm燃料包壳内径dci 8.6mm燃料包壳厚度c 0.57mm燃料芯块直径du 8.19mm燃料棒间距栅距s 12.6mm两个组件间的水隙 0.8mmUO2芯块密度UO2 95%理论密度旁流系数 5%燃料元件发热占总发热份额Fa 97.4%径向核热管因子 1.33轴向核热管因子1.520热流量核热点因子= 2.022热流量工程热点因子1.03焓升工程热点因子<未计入交混因子> 1.142交混因子0.95焓升核热管因子= 1.085堆芯进口局部阻力系数Kin 0.75堆芯出口局部阻力系数Kout 1.0堆芯定位格架阻力系数Kgr 1.05若将堆芯自下而上分为3个控制体.其轴向归一化功率分布见下表： 表 堆芯归一化功率分布轴向等分3个控制体自下而上控制体号123456归一化功率分布0.481.021.501.560.960.48通过计算.得出：1. 堆芯流体出口温度；2. 燃料棒表面平均热流密度以及最大热流密度.平均线功率.最大线功率；3. 热管内的流体温度或焓、包壳表面温度、芯块中心温度随轴向的分布；4. 包壳表面最高温度.芯块中心最高温度；5. DNBR 在轴向上的变化；6. 计算堆芯压降三、 设计正文详细的计算过程、计算结果及分析1.计算过程1.1堆芯流体出口温度平均管按流体平均温度 以及压力由表中查得。假设.查表得经过输入所查程序不断迭代得误差小于0.5。如需更精确的值.可以继续进行迭代计算。1.2燃料表面平均热流密度式中为堆芯燃料棒的总传热面积代入数据得燃料棒表面最大热流密度qmax代入数据得燃料棒平均线功率代入数据得燃料棒最大线功率代入数据得1.3平均管的情况 平均管的流速V式中堆芯内总流通面积为燃料组件内正方形排列时的每一排列的燃料元件数由压力以及流体的平均温度查表得到: 由1.1知.查表得1.4为简化计算起见.假定热管内的流体流速Vh和平均管的V相同。同样.热管四根燃料元件组成的单元通道内的流量代入数据得1.5热管中的计算按一个单元通道计算1热管中的流体温度其中取平均温度对应的参数值.需要进行迭代计算.下面给出第一控制体出口处温度的算法假设,查表得,带入上式与假设误差较大.进行迭代.查表知误差.可以不再进行迭代.就取同理由程序迭代可求得第二控制体出口处流体温度第三控制体出口处流体温度第四控制体出口处流体温度第五控制体出口处流体温度第六控制体出口处流体温度<2>第一个控制体出口处的包壳外壁温度式中：h<z>为单相水强迫对流换热系数.可以利用以下公式来求所以式中流体的、和Pr数根据流体的压力和温度由表查得。如果流体已经达到过冷沸腾.用Jens-Lottes公式：其中为气体的饱和温度.p的单位为MPa.当时.用前面的式子当时.用替换掉代入数据得第一控制体出口处.查表可得故 其余同理由程序计算得出结果如下第二控制体出口处.查表可得故 第三控制体出口处.查表可得故 第四控制体出口处.查表可得故 第五控制体出口处.查表可得故 第六控制体出口处.查表可得故 3包壳内壁温度式中Zr-4的代入数据得：由于与平均温度有关.由程序迭代计算结果如下第一控制体出口处第二控制体出口处 第三控制体出口处 第四控制体出口处 第五控制体出口处 第六控制体出口处 <4>燃料芯块外表面温度式中是包壳与芯块间的气隙等效传热系数.这里取代入数据得第一个控制体出口处第二个控制体出口处第三个控制体出口处第四个控制体出口处第五个控制体出口处第六个控制体出口处<5>燃料芯块中心温度用积分热导求解的方法.即其中代入数据得令由于函数递增.可以通过二分法求解f<x>的根得出二氧化铀中心温度通过编程可求得结果如下第一控制体出口处第二控制体出口处第三控制体出口处第四控制体出口处第五控制体出口处第六控制体出口处1.6热管中的用w-3公式计算.同样对3个控制体都算式中：p为冷却剂工作压力Pa.G为冷却剂质量流密度,为冷却剂通道的当量直径m,为冷却剂的饱和比焓J/kg,为控制体进口处冷却剂的比焓J/kg,为计算点z处的平衡含气量.为其绝对值。平衡含气量的计算式为其中为汽化潜热J/kg。通过程序计算得出结果如下第一控制体出口处第二控制体出口处第三控制体出口处第四控制体出口处第五控制体出口处第六控制体出口处1.7DNBR的计算第一控制体出口处同理可求其余段DNBR.结果如下：第二控制体出口处第三控制体出口处第四控制体出口处第五控制体出口处第六控制体出口处1.8计算热管中的压降单相流体的摩擦压降 式中：使用公式编程分别计算六段控制体的摩擦压降单相流体加速压降：同样使用公式编程分别计算六段控制体的加速压降单相流体提升压降同上局部压降.出口：代入数据得局部压降.进口： 代入数据得局部压降.定位格架出口压降代入数据得以上所使用的比热容和动力粘度都通过软件查询后输入代码中进行计算。通过程序计算结果如下分段压降变化控制体123456提升压降/pa4.3974.3164.1733.9843.8133.708摩擦压降/pa1.0030.9350.4500.4450.4390.435加速压降/pa0.1170.2690.4530.5650.4120.228总压降为2.计算总结温度汇总表控制段控制体外流体出口温度/包壳外表面温度包壳内表面温度芯块表面温度芯块中心温度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.0776.403.5852.7503.7226.864单从表内数据看.DNBR都大于1.芯块中心温度小于二氧化铀的熔点2878.所以理论上能够保证安全性。四、课程设计感想通过这次反应堆热工分析的课程设计.我加深了对反应堆内部传热的了解.同时我发发现了自己的很多不足之处。拿到课程设计题目后.我首先根据题目中的问题去复习热工分析的知识.去了解每个问题该怎么算。然后将计算思路和公式写出来.为下一步的编程做准备。由于大部分的求解都是非线性的方程.所以采用了迭代和二分法求解方程。当把整个问题的求解思路理清和确定计算过程中使用算法后.整个问题就在水和水蒸气热物性如何导入上了。在最开始.我使用通过一定的数据去拟合用到的热物性在定压下与温度的关系。但在拟合之后.通过得出的结果与接用水和水蒸气热物性查询软件得出的结果有一点差距.所以就放弃了。然后.我想将网上通过水和蒸汽热力性质公式IAPWS_IF97计算的源码做成接口.在计算过程中直接用它计算热物性。后来.发现太难.源码有很多地方看不懂.不知道如何接入。于是最终也放弃了。后来就采用了最原始的方法.就是在计算的过程中.通过水和水蒸气热物性查询软件查出结果.一步一步输进去。最后使用程序计算时.在输入参数时耗时太久.我深深的感受到书到用处方恨少.如果原来.把编程多学一点.将整个数据导入直接通过动态链接库完成.就省事太多。这次课设让我把原来学的编程与实际问题结合起来.给我很大感触.原来敲着书上的例子.感觉漫无目的。通过这次自己去思考如何解决一个问题.让我理解了拿到一个问题后.该如何处理。同时通过这次课设.我意识到反应堆热工分析是个复杂而连续的过程.每个参数都受到大量的常数参数的影响也具有很多的修正因子。我们不应该根据自己的常识来判断数据的变化情况.相反地我们应该随时坚持以数据计算为引导.以实验作为验证。仔细客观认真地分析堆内数据的变化.并且对堆内数据进行全程监控.防止堆内数据随时变化.对反应堆的危害性。附录设计流程图、程序1、 程序说明共7个程序由c语言编写.分别为计算流体出口温度.控制体出口流体温度.燃料包壳外壁温度.燃料包壳内壁温度.二氧化铀中心温度.qDNB,压降。使用了迭代和二分法求解部分非线性方程。1堆芯出口温度计算：此段根据任务书给出的基本参数和热量与流量之间关系.运用迭代的算法.求出堆芯的出口温度。2第一至第六控制体的各量计算：因为六个控制体的计算过程类似.这里只说明第一个控制体的计算过程。在现有的参数下.根据热流量与流量的关系和迭代算法.求出该控制体的出口温度。通过流通截面积与湿周的关系求出栅元的当量直径。再根据上面的温度.查出对应的热物性参数由雷诺数与努尔数的关系.解出控制体出口处的对流换热系数。因为不知该处的流体状态.分别用单相强迫对流放热公式和詹斯-洛特斯传热方程算出各自的膜温压.取较小的值加上出口处的流体温度即是包壳的外表面温度。由包壳的外表面的温度再根据圆管的传热方程运用迭代算法解出包壳内表面的温度。芯块与包壳内表面之间的导热问题.根据间隙导热模型.即可解出芯块表面的温度.根据内热源的导热模型.依据积分热导率与温度的对应关系列出方程用二分法解出芯块中心的温度。接下来依据冷却剂的温度.得出的控制体出口处的含汽量。进而依据W-3公式求出该出的临界热流量qDNB.最后得出该出的烧毁比DNBR。3热管的压降计算：热管的压降包括摩擦压降、提升压降、进出口局部压降、定位搁架出口压降。摩擦压降可由计算单相流的达西Darcy公式算得。提升压降可由根据位置的变化算得.其中参数都取平均值。其余的压降根据形阻压降的基本公式再乘以相应的系数求得。最后各项相加得出热管的总压降。2、流程图开始读输入参数计算有关堆参数估算控制体出口温度tf计算控制体出口温度|<0.001重估tf否计算该处含汽量是计算包壳外表面温度根据W-3算临界热流包壳内表面温度芯块表面温度计算烧毁比芯块中心温度打印输出值停机2、程序/ 3.1.cpp : 定义控制台应用程序的入口点。/流体出口温度迭代计算#include"stdafx.h"#include<math.h>#include<iostream>usingnamespace std;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<>double atfout,tfin=287,tav,cp,tfoutc=0;int panding = 1;printf<"请输入假设出口温度n">cin >> atfout;while<panding>=1>tav = <atfout + tfin> / 2;printf<"请输入%f下的Cpn", tav>cin >> cp;tfoutc = tfout<cp, atfout>if <fabs<tfoutc - atfout> <= 0.5>panding = 0;atfout = tfoutc;printf<"出口温度为%fn",tfoutc>while <true>热管内流体温度计算#include"stdafx.h"#include<math.h>#include<iostream>usingnamespace std;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 " << endl;cin >> qav;cout << "请输入归一化参数 " << endl;cin >> gyh;cout << "请输入热管单元通道流量Wht/h" << endl;cin >> wh;while <panding>double tav = <atfout + tfin> / 2;cout << "请输入" << tav << "下的Cp<J/<kg.>>" << endl;cin >> 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 << "请输入流体进口温度" << endl;cin >> tfin;cout << "请输入流体假设流体出口温度" << endl;cin >> atfout;tfout = ffout<tfin, atfout>cout << tfout;while <true>/ 3.3.02.cpp : 定义控制台应用程序的入口点。/计算包壳外壁温度#include"stdafx.h"#include"cmath"#include"iostream"usingnamespace std;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;for <i = 0; i < 6; i+>cout << "请输入" << tfouti << "下的普朗特数" << endl;cin >> pri;cout << "请输入" << tfouti << "下的导热系数W/<m>" << endl;cin >> drxsi;cout << "请输入" << tfouti << "下的动力粘度kg/<ms>" << endl;cin >> 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 <f1i>f2i>fi = f2i;elsefi = f1i;cout << "f1=" << f1i << "f2=" << f2i <<"f="<< fi + tfouti << endl;while <true>/ 3.3.2.cpp : 定义控制台应用程序的入口点。/计算内壁温度#include"stdafx.h"#include"iostream"#include"math.h"usingnamespace std;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 << "请输入归一化参数" << endl;cin >> gyhcs;cout << "请输入平均线功率.单位w/m" << endl;cin >> 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 << "请输入外壁温度" << endl;cin >> tcs;cout << "请输入内壁假设温度"<<endl;cin >> atci;tci = ftci<tcs, atci>cout << tci;while <true>/ 3.3.3.cpp : 定义控制台应用程序的入口点。/计算二氧化优中心温度#include"stdafx.h"#include"iostream"#include"math.h"#include<cmath>usingnamespace std;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 atu>double tu, fnr = 1.33, feq = 1.03, gyhcs, q, tmid, jieguo;bool panding = true;cout << "请输入归一化参数" << endl;cin >> gyhcs;cout << "请输入线热流密度/w" << endl;cin >> 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-tu>0.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 << "请输入燃料芯块外壁温度/" << endl;cin >> tci;cout << "请输入假设二氧化铀中心温度/" << endl;cin >> atu;tu=f2<tci, atu>cout << tu;while <true>/ 3.5.cpp : 定义控制台应用程序的入口点。/计算DNBR#include"stdafx.h"#include"iostream"#include"cmath"#define pi 3.14#define P 15800000usingnamespace std;double fqDNB<double hfout, double hfin, double hfs, double hfg, double g, double dlzj>double x, a1, a2, a3, a, b, c, d, e, qDNB;x = <hfout - hfs> / hfg;a1 = 2.022 - 6.238*pow<10, -8>*P;a2 = 0.1722 - 1.43*pow<10, -8>*P;a3 = exp<<18.177 - 5.987*pow<10, -7>*P>*x>a = 3.154 * pow<10, 6> * <a1 + a2*a3>b = <0.1484 - 1.596*x + 0.1729*x*fabs<x>>*0.2049*g / 1000000 + 1.037;c = 1.157 - 0.869*x;d = 0.2664 + 0.8357*exp<-124 * dlzj>e = 0.8258 + 0.341*pow<10, -6>*<hfs - hfin>qDNB = a*b*c*d*e;return qDNB;int main<>double hfout6, hfin6, gyhcs6,qDNB6,DNBR6, hfs, hfg, dcs = 0.0095, g, s = 0.0126, q, frn=1.33, fqe=1.03,dlzj;int i;dlzj = 4 * <s*s - pi*dcs*dcs / 4> / <pi*dcs>for <i = 0; i < 6; i+>cout << "请输入第" << i+1 << "段控制体出口冷却剂比焓J/kg" << endl;cin >> hfouti;cout << "请输入第" << i +1<< "段进口温度下的冷却剂比焓J/kg" << endl;cin >> hfini;cout << "请输入第" << i+1 << "段控制体归一化参数" << endl;cin >> gyhcsi;cout << "请输入15.8MPa下的冷却剂的饱和比焓J/kg" << endl;cin >> hfs;cout << "请输入15.8MPa下的冷却剂的汽化潜热J/kg" << endl;cin >> hfg;cout << "请输入热管中冷却剂质量流量密度kg/<m*m*h>" << endl;cin >> g;cout << "请输入线热流密度<W/m*m*m>" << endl;cin >> q;for <i = 0; i <6; i+>qDNBi = fqDNB<hfouti, hfini, hfs, hfg, g, dlzj>DNBRi = qDNBi/ <q*fqe*frn*gyhcsi>cout <<"第"<<i+1<<"段控制体qDNB为"<< qDNBi << endl;cout << "第" << i + 1 << "段控制体DNBR为" << DNBRi << endl;while <true>/ 3.6.cpp : 定义控制台应用程序的入口点。/计算堆芯压降#include"stdafx.h"#include"iostream"#include"cmath"usingnamespace std;constdouble g= 9.8;constdouble Ko = 1.0;constdouble Ki = 0.75;constdouble Kgr = 1.05;constdouble L = 3.60;constdouble V = 3.85;constdouble G = 9.88e+6/3600;constdouble dlzj = 1.178e-2;constdouble tfout6 = 291.69, 301.38, 314.86, 327.70, 334.89, 338.22 ;constdouble tfin6 = 287, 291.69, 301.38, 314.86, 327.70, 334.89 ;constdouble tcs6 = 303.51, 326.15, 348.41, 348.43, 348.19, 347.90 ;double fpel<doublemidu> /计算提升压降returnmidu*L*g/6;double fpf<doublef, doubledlzj, doublemidu> /计算摩擦压降returnf*L*midu*V*V / <2 * dlzj*6>double fpa< doublevi, doublevo> /计算加速压降return G*G*<vo - vi>double fpout<doublevo> /计算出口局部压降return Ko*G*G*vo / 2;double fpin< doublevi> /计算进口局部压降return Ko*G*G*vi / 2;double fpgj<doublevi,doublevo>/计算格架出口压降return Kgr*G*G*<vo + vi> / 4;double fmxs<doublere, doubleuw, doubleuf> /计算摩擦系数return 0.184*pow<uw / uf, 0.6> / pow<re, 0.2>double fre<doubledlzj, doublemidu, doubleu> /计算雷诺数returndlzj*midu*V / u;int main<>double pel6, pf6, pa6, pout, pin, pgj, mxs6, re6, midu6, p6, pt = 0;int i; double birong6 = 0.001337049414, 0.001362247294, 0.001409055895, 0.001475639644, 0.001541851039, 0.001585457905 ;double vi6 = 0.001329349786, 0.001344993020, 0.001380774324, 0.001440939652, 0.001515991823, 0.001570780783 ;double vo6 = 0.001344993020, 0.001380774324, 0.001440939652, 0.001515991823, 0.001570780783, 0.001601107797 ;double uw6 = 8.724724e-005, 7.803018e-005, 2.328906e-005, 2.328904e-005, 2.328945e-005, 2.329043e-005 ;double uf6 = 9.277614e-005, 8.996754e-005, 8.543740e-005, 8.011150e-005, 7.572227e-005, 7.319442e-005 ;for <i = 0; i < 6; i+>midui = 1 / birongi;rei = midui*V*dlzj/ufi;mxsi = 0.184*pow<uwi / ufi, 0.6>/pow<rei,0.2>pout = fpout<vo5>pin = fpin<vi0>pgj = fpgj<vi0, vo5>for <i = 0; i < 6; i+>peli = fpel<midui>pfi = mxsi * L*midui * V / <dlzj * 6 * 2>pai = fpa<vii,voi>pi = peli + pfi + pai;pt = pt + pi;pt = pt + +pout + pin + pgj;for <i = 0; i < 6; i+>cout <<i+1<<"段提升"<< peli << endl;cout << i+1<< "段摩擦"<<pfi << endl;cout << i+1 << "段加速"<< pai << endl;cout <<"出口"<< pout << endl;cout <<"入口"<< pin << endl;cout <<"格架"<< pgj << endl;cout << pt << endl;while <true>
