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

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

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

3、19mm燃料棒间距（栅距）s 12.6mm两个组件间的水隙 0.8mmUO2芯块密度UO2 95%理论密度旁流系数 5%燃料元件发热占总发热份额Fa 97.4%径向核热管因子FRN 1.33轴向核热管因子FZN 1.520热流量核热点因子FqN=FRN FZN 2.022热流量工程热点因子FqE 1.03焓升工程热点因子FHE(未计入交混因子) 1.142交混因子FH·mE 0.95焓升核热管因子FHN=FRN 1.085堆芯进口局部阻力系数Kin 0.75堆芯出口局部阻力系数Kout 1.0堆芯定位格架阻力系数Kgr 1.05 若将堆芯自下而上分为3个控制体，其轴向归一化功率分布见

4、下表： 表 堆芯归一化功率分布（轴向等分3个控制体）自下而上控制体号123456归一化功率分布0.481.021.501.560.960.48通过计算，得出：1. 堆芯流体出口温度；2. 燃料棒表面平均热流密度以及最大热流密度，平均线功率，最大线功率；3. 热管内的流体温度（或焓）、包壳表面温度、芯块中心温度随轴向的分布；4. 包壳表面最高温度，芯块中心最高温度；5. DNBR 在轴向上的变化；6. 计算堆芯压降三、 设计正文（详细的计算过程、计算结果及分析）1.计算过程1.1堆芯流体出口温度（平均管）tf,out=tf,in+FaNtW(1-)Cp Cp按流体平均温度tf=12(tf,in+

5、tf,out) 以及压力由表中查得。假设tf,out=330，查表得Cp=5.610kJ/(kg)经过输入所查Cp程序不断迭代得tf,out=323.9误差小于0.5。如需更精确的值，可以继续进行迭代计算。1.2燃料表面平均热流密度qq=FaNt/F总 式中F总为堆芯燃料棒的总传热面积F总=mndcsL 代入数据得F总=121×265××9.5×10-3×3.60=3443.40 m2q=97.4%×1820×1063443.40=5.29×105 W/m2燃料棒表面最大热流密度qmax qmax=qFqNFqE

6、代入数据得 qmax=5.29×105×2.022×1.03=1.10×106 W/m2燃料棒平均线功率ql ql=qdcsLL=qdcs 代入数据得ql=5.29×105××9.5×10-3=1.57×104 W/m燃料棒最大线功率ql,maxql,max=qlFqNFqE 代入数据得ql,max=qlFqNFqE=1.57×104×2.022×1.03=3.26×104 W/m1.3平均管的情况 平均管的流速VV=W(1-)Atf 式中At堆芯内总流通面积At=

7、mn0×n0s2-4dcs2+m4n0s×2 n0为燃料组件内正方形排列时的每一排（列）的燃料元件数f由压力以及流体的平均温度tf查表得到: f=1vf 由1.1知tf=323.9+2872=305.5，查表得vf= 0.001397680614m3/kgf=10.001397680614=715.471 kg/m3At=121×17×17×12.6×10-32-4×9.5×10-32+121×4×17×12.6×10-3×0.8×10-32=3.11 m

8、2V=32500×(1-5%)3.11×715.471×3.6=3.85 m/s 1.4为简化计算起见，假定热管内的流体流速Vh和平均管的V相同。同样，热管四根燃料元件组成的单元通道内的流量Wb=W(1-)AtAb Ab=s2-4dcs2 代入数据得Ab=12.6×10-32-49.5×10-32=0.88×10-4 m2Wb=32500×（1-5%）3.11×0.88×10-4=0.87 t/h1.5热管中的计算（按一个单元通道计算）（1）热管中的流体温度tf(z)=tf,in+qFRNFHEFHmEd

9、csWbCp0z(z)dztf(z)=287+5.29×105×1.35×1.142×0.95×9.5×10-30.87×10003600Cp0zzdz=287+9.56×104Cp0zzdz其中Cp取平均温度对应的参数值，需要进行迭代计算，下面给出第一控制体出口处温度的算法假设tfL6=300,查表得Cp=5.3348 kJ/(kg),带入上式tfL6=287+9.56×1045.3348×1000×0.8×3.606=291.6 与假设误差较大，进行迭代，查表知Cp=5.

10、222 kJ/(kg)tfL6=287+9.56×1045.222×1000×0.8×3.663=291.69 误差291.60-291.69<0.5，可以不再进行迭代，就取tfL6=291.69 同理由程序迭代可求得第二控制体出口处流体温度tf2L6=301.38 第三控制体出口处流体温度tf3L6=314.86 第四控制体出口处流体温度tf4L6=327.70 第五控制体出口处流体温度tf5L6=334.89第六控制体出口处流体温度tfL=338.22 (2)第一个控制体出口处的包壳外壁温度tcs(z)=tf(z)+f1(z)=tf(z)+qF

11、RNzFqEhz式中：h(z)为单相水强迫对流换热系数W/(m2)，可以利用以下公式来求Nu=hzD=0.023Re0.8Pr0.4所以 hz=0.023Re0.8Pr0.4D式中Re=GD=WbAbDD=4AbU=4(s2-4dcs2)dcs流体的、和Pr数根据流体的压力和温度由表查得。如果流体已经达到过冷沸腾，用Jens-Lottes公式：f2(z)=ts+25qFRNFqEz1060.25e-p6.2-tf(z)其中ts为气体的饱和温度，p的单位为MPa，p=15.8MPa时，ts=346.38当f2f1时，用前面的式子当f2<f1时，用f2替换掉f1代入数据得D=4×1

12、2.6×10-32-4×9.5×10-32×9.5×10-3=11.78×10-3 mRe=0.863.60.88×10-4×11.78×10-3=32.35hz=0.023×32.350.8Pr0.411.78×10-3=31.51Pr0.40.8f1z=5.29×105×1.33z1.03hz=2.25×1040.8zPr0.4f2z=346.38+255.29×105×1.33×1.10z1060.25e-15.86.2

13、-tfz=346.38+1.83z0.25-tfz第一控制体出口处tfL6=291.69 ，查表可得=0.57799W/(m) L6=0.48 =9.1858×10-5 kg/(ms) Pr=0.8371f1L3=2.25×104×9.1858×10-50.8×0.80.57799×0.83710.4=11.82f2L3=346.38+1.82×0.480.25-303.91=56.21 故 tcsL6=tfL6+f1L6=291.69+11.82=303.51 其余同理由程序计算得出结果如下第二控制体出口处tf2L6=30

14、1.38 ，查表可得=0.5601W/(m) 2L6=1.02 =8.8×10-5 kg/(ms) Pr=0.8621f12L6=24.79 f22L6=46.84 故 tcs2L6=tf2L6+f12L6=301.38+24.79=326.15 第三控制体出口处tf3L6=314.86，查表可得=0.5333 W/(m) L=1.50 =8.275×10-5 kg/(ms) Pr=0.9157f1L=35.55 f2L=33.55 故 tcs3L6=tfL+f2L=348.41 第四控制体出口处tf4L6=327.70，查表可得=0.5044 W/(m) L=1.56 =

15、7.735×10-5 kg/(ms) Pr=1.0051f1L=35.67 f2L=20.73 故 tcs4L6=tfL+f2L=348.43 第五控制体出口处tf4L6=334.89，查表可得=0.4868W/(m) L=0.96 =7.401×10-5 kg/(ms) Pr=1.0873f1L=21.28 f2L=13.30故 tcs5L6=tfL+f2L=348.19 第六控制体出口处tfL=338.22，查表可得=0.4782 W/(m) L=0.48 =7.235×10-5 kg/(ms) Pr=0.9157f1L=11.39 f2L=9.68故 tcs

16、L=tfL+f2L=347.90 （3）包壳内壁温度tci(z)=tcs(z)+qlFRNFqE(z)2kc(z)lndcsdci式中Zr-4的kc=0.00547×1.8×tc+32+13.8 W/mtcz=12tcs(z)+tci(z)代入数据得：tciz=tcsz+1.57×104×1.33×1.03z2kczln9.58.36由于kc与平均温度有关，由程序迭代计算结果如下第一控制体出口处 tci1L6=313.13 第二控制体出口处 tci2L6=346.26 第三控制体出口处 tci3L6=377.54 第四控制体出口处 tci4L6

17、=378.72第五控制体出口处 tci5L6=366.89 第六控制体出口处 tciL=357.27 (4)燃料芯块外表面温度tu(z)=tci(z)+qlFRNFqE(z)dci+du2hg式中hg是包壳与芯块间的气隙等效传热系数，这里取hg=5678 W/(m2)代入数据得tuz=tciz+1.57×104×1.33×1.03(z)8.6+8.192×5678×10-3=tciz+143.70(z)第一个控制体出口处tuL6=313.13+143.70×0.48=384.03 第二个控制体出口处tu2L6=346.26+143.7

18、0×1.02=492.83第三个控制体出口处tu3L6=377.54+143.70×1.50=593.09第四个控制体出口处tu4L6=378.72+143.70×1.56=602.89 第五个控制体出口处tu5L6=366.89+143.70×0.96=504.84 第六个控制体出口处tuL=357.27+143.70×0.48=426.25(5)燃料芯块中心温度用积分热导求解的方法，即0t0zkutdt=0tuzkutdt+qlFRNFqEz4其中 ku=38.24t+402.55+4.788×(t+273.15)3代入数据得0t

19、0zkutdt=0tuzkutdt+1.57×104×1.33×1.03z4=0tuzkutdt+17.12(z) W/cm0tkutdt=0t38.24t+402.55+4.788×t+273.153dt=38.24lnt+402.55402.55+4.788×10-134×(t+273.15)4-273.154令 fx=0t0zkutdt-0tuzkutdt+17.12(z)由于函数递增，可以通过二分法求解f(x)的根得出二氧化铀中心温度通过编程可求得结果如下第一控制体出口处 t0L6=568.79 第二控制体出口处 t02L6=

20、974.83 第三控制体出口处 t03L6=1383.57 第四控制体出口处 t04L6=1430.38第五控制体出口处 t05L6=958.77第六控制体出口处 t0L=619.90 1.6热管中的qDNB(z)用w-3公式计算，同样对3个控制体都算qDNB=3.154×1062.022-6.238×10-8p+(0.1722-1.43×10-8p)×exp18.177-5.987×10-7pxe0.1484-1.596xe+0.1729xexe×0.2049G106+1.037(1.157-0.869xe)0.2664+0.835

21、7exp(-124D)×0.8258+0.341×10-6(hfs-hf,in)式中：p为冷却剂工作压力（Pa），G为冷却剂质量流密度kg/(m2h),D为冷却剂通道的当量直径（m）,hfs为冷却剂的饱和比焓（J/kg）,hf,in为控制体进口处冷却剂的比焓（J/kg）,xe为计算点z处的平衡含气量，xe为其绝对值。平衡含气量xe的计算式为xe=hf(z)-hfshfg其中hfg为汽化潜热（J/kg）。hfs=1641.4kJ/kg hfg=945kJ/kg p=15.8MPa D=11.78×10-3mG=WbAb=0.87×10000.88×

22、;10-4=9.88×106kg/(m2h)通过程序计算得出结果如下第一控制体出口处qDNBL6=5.244×106W/m2第二控制体出口处qDNB2L6=4.731×106W/m2第三控制体出口处qDNB3L6=3.897×106W/m2第四控制体出口处qDNB4L6=3.109×106W/m2第五控制体出口处qDNB5L6=2.624×106W/m2第六控制体出口处qDNBL=2.388×106W/m21.7DNBR的计算DNBR(z)=qDNB(z)qFRNFqE(z)第一控制体出口处DNBRL6=5.244×

23、;1065.29×105×1.33×1.03×0.48=15.077同理可求其余段DNBR，结果如下：第二控制体出口处DNBR2L6=6.40第三控制体出口处DNBR3L6=3.585第四控制体出口处DNBR4L6=2.750第五控制体出口处DNBR5L6=3.772第六控制体出口处DNBRL=6.8641.8计算热管中的压降单相流体的摩擦压降 pf=fLDV22 式中： f=fisowfn=0.184Re0.2(wf)0.6使用公式编程分别计算六段控制体的摩擦压降单相流体加速压降：pa=G2vout-vin同样使用公式编程分别计算六段控制体的加速压降单

24、相流体提升压降pel=gL同上局部压降，出口：pout=koutoutVout22=koutG2vout2代入数据得pout=1.0×（9.88×1063600）2×0.001667672=6.03×103 Pa局部压降，进口： pin=kininVin22=kinG2vin2代入数据得pin=0.75×（9.88×1063600）2×0.001329352=5.0×103 Pa局部压降，定位格架出口压降pgr=kgrgrV22=kgrG212vout+vin2代入数据得pgr=1.05×（9.88

25、15;1063600）2×0.00166767+0.001329354=5.79×103 Pa以上所使用的比热容和动力粘度都通过软件查询后输入代码中进行计算。通过程序计算结果如下分段压降变化控制体123456提升压降/103pa4.3974.3164.1733.9843.8133.708摩擦压降/103pa1.0030.9350.4500.4450.4390.435加速压降/103pa0.1170.2690.4530.5650.4120.228总压降为p=pout+pgr+pin+pel+pf+pa=4.6981×104Pa2.计算总结温度汇总表控制段控制体外流体

26、出口温度/包壳外表面温度包壳内表面温度 芯块表面温度 芯块中心温度 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临界热流密度×106W/m25.2444.7313.8973.1092.6242.388烧毁比15.077

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

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

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

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

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

32、据形阻压降的基本公式再乘以相应的系数求得。最后各项相加得出热管的总压降。2、流程图开始读输入参数计算有关堆参数估算控制体出口温度tf计算控制体出口温度|1-tf|<0.001重估tf否计算该处含汽量是计算包壳外表面温度根据W-3算临界热流包壳内表面温度芯块表面温度计算烧毁比芯块中心温度打印输出值停机2、程序/ 3.1.cpp : 定义控制台应用程序的入口点。/流体出口温度迭代计算#include "stdafx.h"#include <math.h>#include<iostream>using namespace std;double tfo

33、ut(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

34、+ 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)/3.3.1热管内流体温度计算#include "stdafx.h"#include <math.h>#include<iostream>using

35、namespace 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 <

36、;< "请输入热管单元通道流量Wh（t/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 (fa

37、bs(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;wh

38、ile (true)/ 3.3.02.cpp : 定义控制台应用程序的入口点。/计算包壳外壁温度#include "stdafx.h"#include "cmath"#include "iostream"using namespace 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;i

39、nt i;for (i = 0; i < 6; i+)cout << "请输入" << tfouti << "下的普朗特数" << endl;cin >> pri;cout << "请输入" << tfouti << "下的导热系数（W/(m)）" << endl;cin >> drxsi;cout << "请输入" << tfouti <&

40、lt; "下的动力粘度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 <&

41、lt;"f="<< fi + tfouti << endl;while (true)/ 3.3.2.cpp : 定义控制台应用程序的入口点。/计算内壁温度#include "stdafx.h"#include"iostream"#include"math.h"using namespace std;double ftci(double tcs,double atci)double frn = 1.33, feq = 1.03, gyhcs, dcs=9.5, dci=8.6, kc, tav

42、,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

43、(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 : 定义控制台应用程序的入口

44、点。/计算二氧化优中心温度#include "stdafx.h"#include"iostream"#include"math.h"#include<cmath>using namespace 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, doubl

45、e 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 = fa

46、lse;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;el

47、setu = 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"

48、;#include"iostream"#include"cmath"#define pi 3.14#define P 15800000using namespace 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.1

49、722 - 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 qD

50、NB;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 &g

51、t;> hfouti;cout << "请输入第" << i +1<< "段进口温度下的冷却剂比焓（J/kg）" << endl;cin >> hfini;cout << "请输入第" << i+1 << "段控制体归一化参数" << endl;cin >> gyhcsi;cout << "请输入15.8MPa下的冷却剂的饱和比焓（J/kg）" <<

52、; 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(

53、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 : 定义控制台应用程序的

54、入口点。/计算堆芯压降#include "stdafx.h"#include"iostream"#include"cmath"using namespace std;const double g= 9.8;const double Ko = 1.0;const double Ki = 0.75;const double Kgr = 1.05;const double L = 3.60;const double V = 3.85;const double G = 9.88e+6/3600;const double dlzj = 1.17

55、8e-2;const double tfout6 = 291.69, 301.38, 314.86, 327.70, 334.89, 338.22 ;const double tfin6 = 287, 291.69, 301.38, 314.86, 327.70, 334.89 ;const double tcs6 = 303.51, 326.15, 348.41, 348.43, 348.19, 347.90 ;double fpel(double midu) /计算提升压降return midu*L*g/6;double fpf(double f, double dlzj, double midu) /计算摩擦压降return