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

、设计要求在设计反应堆冷却系统时，为了保证反应堆运行安全可靠，针对不同的堆型，预先规定了热工设计必须遵守的要求，这些要求通常就称为堆的热工设计准则。目前压水动力堆设计中所规定的稳态热工设计准则，一般有以下几点：燃料元件芯块内最高应低于其他相应燃耗下的熔化温度；燃料元件外表面不允许发生沸腾临界；必须保证正常运行工况下燃料元件和堆内构件得到充分冷却；在事故工况下能提供足够的冷却剂以排除堆芯余热；在稳态额定工况和可预计的瞬态运行工况中，不发生流动不稳定性。在热工设计中，通常是通过平均通道（平均管）可以估算堆芯的总功率，而热通道（热管）则是堆芯中轴向功率最高的通道，通过它确定堆芯功率的上限，热点是堆芯中温度最高的点，代表堆芯热量密度最大的点，通过这个点来确定DNBR。二、设计任务某压水反应堆的冷却剂和慢化剂都是水，用二氧化铀作燃料，Zr-4作燃料包壳材料。燃料组件无盒壁，燃料元件为棒状，正方形排列，已知下列参数：系统压力P15.8MPa堆芯输出热功率Nt1820MW冷却剂总流量W32500t/h反应堆进口温度tfin287°C堆芯高度L3.60m燃料组件数m121燃料组件形式n0Xn017X17每个组件燃料棒数n265燃料包壳外径des9.5mm燃料包壳内径dei8.6mm燃料包壳厚度6e0.57mm燃料芯块直径du8.19mm燃料棒间距（栅距）s12.6mm两个组件间的水隙60.8mmU02芯块密度PU0295%理论密度旁流系数Z5%燃料元件发热占总发热份额Fa97.4%径向核热管因子硝1.33轴向核热管因子F£1.520热流量核热点因子即曙僧2.022热流量工程热点因子Ffq1.03焓升工程热点因子厲未计入交混因子）1.142交混因子FEHm0.95焓升核热管因子臨FN1.085堆芯进口局部阻力系数Kin0.75堆芯出口局部阻力系数Kout1.0堆芯定位格架阻力系数Kgr1.05若将堆芯自下而上分为3个控制体，其轴向归一化功率分布见下表：表堆芯归一化功率分布（轴向等分3个控制体）自下而上控制体号123456归一化功率分布0.481.021.501.560.960.48通过计算，得出：堆芯流体出口温度；燃料棒表面平均热流密度以及最大热流密度，平均线功率，最大线功率；热管内的流体温度（或焓）、包壳表面温度、芯块中心温度随轴向的分布包壳表面最高温度，芯块中心最高温度；DNBR在轴向上的变化；计算堆芯压降三、 设计正文（详细的计算过程、计算结果及分析）1.计算过程1.1堆芯流体出口温度（平均管）F•Nt=t+ atf,out f,in W•（1—（）•G'“ 卩q按流体平均温度牛=2（?m+tf,out）以及压力由表中查得。假设tfout=330T,查表得q=5.610kJ/（kg・。0经过输入所查乔程序不断迭代得niit=323.9弋误差小于0.5。「P J,OUt如需更精确的值，可以继续进行迭代计算。1.2燃料表面平均热流密度@q=F•N/F“at'总式中F为堆芯燃料棒的总传热面积总F=m・n・兀•d■厶总 CS代入数据得F=121x265xux9.5x10-3x3.60=3443.40m2总

q=97.4%X1820X10q=97.4%X1820X1063443.40=5.29X105"/m2燃料棒表面最大热流密度qmaxq=q•Fn•FeLmaxLqq代入数据得q=5.29X105x2.022X1.03=1.10X106^/m2

丄max /燃料棒平均线功率&0=可兀"des』=@•兀•d5Zes代入数据得qt=5.29X105x冗X9.5X103=1.57X104"/m燃料棒最大线功率qlmaxa=7•Fn•Fe1l,maxllqq代入数据得a=7•Fn•Fe=1.57X104x2.022X1.03=3.26X104"/mTOC\o"1-5"\h\zLl,max Llqq 11.3平均管的情况平均管的流速V式中4堆芯内总流通面积「 兀 18-=m•(n0Xn0)[S2 畀es2]+m4(n0•sX2)_%为燃料组件内正方形排列时的每一排(列)的燃料元件数石由压力以及流体的平均温度5查表得到：TOC\o"1-5"\h\z1由i.i知仁=323.9+287=305.5T，查表得耳=0.001397680614m3/gJ 2 J1作=0.001397680614=715.471kg/m3「 冗 1坷=121X(17X17)X[(12.6X103)2_x(9.5X103)2+121X(4X17X12.6X103X0.8+121X(4X17X12.6X103X0.8X103)=3.11m2V「补00：(鳥5%「=3.85M3.11X715.471X3.61.4为简化计算起见，假定热管内的流体流速Vh和平均管的V相同。同样，热管四根燃料元件组成的单元通道内的流量代入数据得n每=(12.6X10—3)2—才(9.5X10—3)2=0.88X10—4m232500X(1—5%)亿= — X0.88X10—4=0.87t/hb 3.111.5热管中的计算(按一个单元通道计算)(1)热管中的流体温度q•Fn•Fe•Fe兀drzt(z)=t+R纱，Em csI0(z)dzTOC\o"1-5"\h\z/ /，in •C Jbp 0X105x1.35X1.142X0.95兀X9.5X10—3(z$(z)=287+ 0.87X1000厂 J0(z)血3600 •Cp 09.56X104cz

=287+ Ig(z)dzp丿0其中J取平均温度对应的参数值，需要进行迭代计算,下面给出第一控制体出口处温度的算法假设。G)=300OC，查表得C=5.3348灯/(耳•弋),带入上式>6 P9.56X104287+9.56X104287+5.3348X1000X0.8X3.60~~6-=291.6C与假设误差较大，进行迭代，查表知J=5.222k//(®•C)L(6)=287+9.56X10L(6)=287+9.56X1045.222X1000X0.8X3.66~3~=291.69C误差|291.60—291.691<0.5C,可以不再进行迭代，就取"G)=291.69CJ6同理由程序迭代可求得第二控制体出口处流体温度"(互)=301.38CJ6第三控制体出口处流体温度$仔）=314.86T第四控制体出口处流体温度$（严）=327.70T第五控制体出口处流体温度$（学）=334.89T第六控制体出口处流体温度$（厶）=338.22T2）第一个控制体出口处的包壳外壁温度式中：h（z）为单相水强迫对流换热系数［W/（m2-T）］,可以利用以下公式来求h（z）DN= '丿£=0.023^^0-8-Pr0.4所以式中h(z)=所以式中h(z)=0.023Re0-8Pr0-4・_LG・DWDRe= £=b•_£D=4%=4©-4吧)£“兀s流体的入、口和Pr数根据流体的压力和温度由表查得。如果流体已经达到过冷沸腾，用Jens-Lottes公式：106些2(z)=J+25(吟橙◎).总-片(z)106其中。为气体的饱和温度，p的单位为MPa,p=15.8MPa时，ts=346.38T当铃2二岂1时，用前面的式子当△牛2<"门时，用△牛2替换掉"门代入数据得=11.78x10-3m4x((12.6x10-3)2一4X(9.5x10-3=11.78x10-3m兀x9.5x10-3

0.86■36"Re0.86■36"Re=0.88x10411.78x103 32.35Xh(z)=0.023x几=31.51^^1TOC\o"1-5"\h\z仏 11.78X103 “0.8X105x1.33•0(z)•1.03 u0.8m(z)h(z)△备⑵= / =225x104h(z)106许(沪346.38+25(5.2955電3"10106=346.38+1.83•g(z)o.25 $(z)第一控制体出口处rG)=291.69T，查表可得入=0.57799W/(m•T)申G)=0.48J6 6|i=9.1858x105g/(m・s) Pr=0.8371L气1(L气1(3)=2.25x104x(9.1858x105)0.8x0.80.57799x0..83710.4=11.82T△冬2(|)=346.38+1.82x0.480.25 303.91=56.21T故厶 厶 Ltcs（召）=tf（召）+△蛰（召）=291.69+11.82=303.51T其余同理由程序计算得出结果如下第二控制体出口处r(垒)=301.38T，查表可得入=0.5601W/(m•T)申(垒)=1.02J6 6|i=8.8x105k^/(m•s)Pr=0.86212厶△冬］(石)=24.79T2厶A0 (丁)=46.84T/26故2厶 2厶 2厶t(―)=t(―)+A0(―)=301.38+24.79=326.15Tcs6 /6 J16第三控制体出口处t(理)=314.86T，查表可得入=0.5333W/(m•T)申(厶)=1.50f6卩=8.275x105k^/(m•s)Pr=0.9157轧(L)=35.55T△0(厶)=33.55TJ2

3Lt—=tL+LGL=348.41°Ccs6f f2第四控制体出口处t如=327.70C，查表可得入=0.5044W/m•C申厶=1.56f6|i=7.735x105kg/m•sPr=1.0051厶=35.67C厶=20.73C/2故4厶t=tL+△0L=348.43Ccs6f/2第五控制体出口处t如=334.89C，查表可得入=0.4868W/m•C申厶=0.96J6|i=7.401x105kg/m•sPr=1.0873L=21.28CL=13.30C5LIs~6下L+△作2L=348.19C第六控制体出口处tfL=338.22C,查表可得入=0.4782W/m•C申厶=0.48|i=7.235x105kg/m•sPr=0.9157L6L=11.39Cf1△0厶=9.68CJ2=t=tfL+△3f2L=347.90CJ233）包壳内壁温度式中Zr-4的仪tciz=t式中Zr-4的仪tciz=tosZ=0.00547dIn严

a.Cl2兀kZc1.8xf+32+13.8^/m•C12[tcsZ+tciZ]代入数据得：代入数据得：tzci=ttzci=tcs*1.57x104x1.33x倔阳加黑8.362nkzc第一个控制体出口处第一个控制体出口处第一控制体出口处乙（严）=313.13T第二控制体出口处tci（垒）=346.26T第三控制体出口处tci（半）=377.54T第四控制体出口处乙（严）=378.72T第五控制体出口处J（严）=366.89T第六控制体出口处tci（L）=357.27T4） 燃料芯块外表面温度式中鶴是包壳与芯块间的气隙等效传热系数，这里取鶴=5678"/（尬2•T）代入数据得t（-）=313.13+143.70X0.48=384.03°Cu\6'第二个控制体出口处=346.26+143.70X1.02=492.83C第三个控制体出口处t=377.54+143.70X1.50=593.09C第四个控制体出口处第五个控制体出口处t=366.89+143.70X0.96=504.84C

第六个控制体出口处t(L)=357.27+143.70X0.48=426.25°Cu5）燃料芯块中心温度用积分热导求解的方法，即「叫（咖」叫（酗+祝警⑵/ u / u 4〃00其中k=38-24+4.788X（t+273.15）3ut+402.55代入数据得g必=gdt+1.57X104x1.33X1.03g(z)g必=gdt+4兀rt(z)Iuk^(t)dt+17.12g(z)"/cm0ft ft38.24TOC\o"1-5"\h\zI乞(±)dt=It+40255+4.788X(t+273.15)3dt丿0 丿0 .t+402.55 4.788X1013=38.24加“cll+ X(t+273.15)4 273.154402.55 4令f(x)=R(z)k(t)dtR(z)花(r)dt+17.12卩(z)J0 u J0 u由于函数递增，可以通过二分法求解f(x)的根得出二氧化铀中心温度通过编程可求得结果如下第一控制体出口处tG)=568.79C06第二控制体出口处t（丝）=974.83C06第三控制体出口处r（理）=1383.57C06第四控制体出口处r（丝）=1430.38C06第五控制体出口处t（匹）=958.77C06第六控制体出口处q（厶）=619.90°C1.6热管中的⑵用w-3公式计算，同样对3个控制体都算qDNB=3.154X106{(2.022-6.238X10-8p)+(0.1722一1.43X10-8p)Xexp[(18.177—5.987X10—7卩)化]}[(0.1484—1.596化+0.1729%e|xe|)X。.囂“+1.037](1.157—0.869耳)[0.2664+0.8357exp(—124DjX［0.8258+0.341X10—6（件$—件曲）］式中：p为冷却剂工作压力（Pa）,G为冷却剂质量流密度［g/（m2・h）］,q为冷却剂通道的当量直径（m）,怯为冷却剂的饱和比焓（J/kg）片加为控制体进口处冷却剂的比焓（J/kg），耳为计算点z处的平衡含气量，|%|为其绝对值。平衡含气量墓的计算式为件（Z）—件S其中冷9为汽化潜热(J/kg)。h=1641.4k//k^ h=945kJ/kg p=15.8MPa/s jgD=11.78X10—3m8W, 0.87X1000G=r0.88X10—4 9.88X106S/(2h)通过程序计算得出结果如下第一控制体出口处厶你丽(6)=5.244X106W/m2第二控制体出口处2厶你曲(石)=4.731X106W/m2第三控制体出口处3厶你丽("6)=3.897X106W/m2第四控制体出口处

第五控制体出口处第六控制体出口处第五控制体出口处第六控制体出口处q =3.109x106W/m2^DNB65Lq——=2.624x106"/m2^DNB6qL=2.388x106"/m21.7DNBR的计算DNBRzqz5NB

qFNFE^zRq第一控制体出口处5.244x106 1.7DNBR的计算DNBRzqz5NB

qFNFE^zRq第一控制体出口处5.244x106 =15.077x105x1.33x1.03x0.48同理可求其余段DNBR,结果如下：第二控制体出口处2LDNBR==6.406第三控制体出口处3LDNBR丁=3.5856第四控制体出口处4LDNBR丁=2.7506第五控制体出口处5LDNBR==3.7726DNBR第六控制体出口处DNBRL=6.8641.8计算热管中的压降单相流体的摩擦压降P卩2~2式中：f_f %:几。石0.184心0.2石W0.6使用公式编程分别计算六段控制体的摩擦压降单相流体加速压降：=G2(Vout同样使用公式编程分别计算六段控制体的加速压降单相流体提升压降叽=D•9•厶同上局部压降，出口：△Pout=△Pout=kout0V2厂outout=上outG2Vout~2~代入数据得△p=1.0x△p=1.0xout(9.88x106x0.00166767=6.03x103Pa局部压降，进口：△Pin局部压降，进口：△PinG2V.=k inin-2~代入数据得9-88x109-88x106)△咕=0.75x(16002x0-00132935=5.0x103Pg局部压降，定位格架出口压降△卩=k局部压降，定位格架出口压降△卩=k厂gr gr=kgr1G2(v+V.)2outin2代入数据得△p=1.05x△p=1.05xgr9.88x106"3600 )2X0-00166767+0-0O132935=5.79x103Pa以上所使用的比热容和动力粘度都通过软件查询后输入代码中进行计算。通过程序计算结果如下分段压降变化控制体123456提升压降/103pa4.3974.3164.1733.9843.8133.708摩擦压降/103pa1.0030.9350.4500.4450.4390.435加速压降〃103pa0.1170.2690.4530.5650.4120.228总压降为Yp=Ap+Ap+An + + =4.6981x104PaJ丄 ’outLgr1mLelL/ 1a2.计算总结温度汇总表控制段控制体外流体出口温度/C包壳外表面温度C包壳内表面温度芯块表面温度芯块中心温度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临界热流密度x106"/m25.2444.7313.8973.1092.6242.388烧毁比15.0776.403.5852.7503.7226.864单从表内数据看，DNBR都大于1,芯块中心温度小于二氧化铀的熔点2878°C,所以理论上能够保证安全性。四、课程设计感想通过这次反应堆热工分析的课程设计，我加深了对反应堆内部传热的了解，同时我发发现了自己的很多不足之处。拿到课程设计题目后，我首先根据题目中的问题去复习热工分析的知识，去了解每个问题该怎么算。然后将计算思路和公式写出来，为下一步的编程做准备。由于大部分的求解都是非线性的方程，所以采用了迭代和二分法求解方程。当把整个问题的求解思路理清和确定计算过程中使用算法后，整个问题就在水和水蒸气热物性如何导入上了。在最开始，我使用通过一定的数据去拟合用到的热物性在定压下与温度的关系。但在拟合之后，通过得出的结果与接用水和水蒸气热物性查询软件得出的结果有一点差距，所以就放弃了。然后，我想将网上通过水和蒸汽热力性质公式IAPWS_IF97计算的源码做成接口，在计算过程中直接用它计算热物性。后来，发现太难，源码有很多地方看不懂，不知道如何接入。于是最终也放弃了。后来就采用了最原始的方法，就是在计算的过程中，通过水和水蒸气热物性查询软件查出结果，一步一步输进去。最后使用程序计算时，在输入参数时耗时太久，我深深的感受到书到用处方恨少，如果原来，把编程多学一点，将整个数据导入直接通过动态链接库完成，就省事太多。这次课设让我把原来学的编程与实际问题结合起来，给我很大感触，原来敲着书上的例子，感觉漫无目的。通过这次自己去思考如何解决一个问题，让我理解了拿到一个问题后，该如何处理。同时通过这次课设，我意识到反应堆热工分析是个复杂而连续的过程，每个参数都受到大量的常数参数的影响也具有很多的修正因子。我们不应该根据自己的常识来判断数据的变化情况，相反地我们应该随时坚持以数据计算为引导，以实验作为验证。仔细客观认真地分析堆内数据的变化，并且对堆内数据进行全程监控，防止堆内数据随时变化，对反应堆的危害性。附录（设计流程图、程序）1、程序说明共7个程序由c语言编写，分别为计算流体出口温度，控制体出口流体温度，燃料包壳外壁温度，燃料包壳内壁温度，二氧化铀中心温度，qDNB，压降。使用了迭代和二分法求解部分非线性方程。（1） 堆芯出口温度计算：此段根据任务书给出的基本参数和热量与流量之间关系，运用迭代的算法，求出堆芯的出口温度。（2） 第一至第六控制体的各量计算：因为六个控制体的计算过程类似，这里只说明第一个控制体的计算过程。在现有的参数下，根据热流量与流量的关系和迭代算法，求出该控制体的出口温度。通过流通截面积与湿周的关系求出栅元的当量直径。再根据上面的温度，查出对应的热物性参数由雷诺数与努尔数的关系，解出控制体出口处的对流换热系数。因为不知该处的流体状态，分别用单相强迫对流放热公式和詹斯-洛特斯传热方程算出各自的膜温压，取较小的值加上出口处的流体温度即是包壳的外表面温度。由包壳的外表面的温度再根据圆管的传热方程运用迭代算法解出包壳内表面的温度。芯块与包壳内表面之间的导热问题，根据间隙导热模型，即可解出芯块表面的温度，根据内热源的导热模型，依据积分热导率与温度的对应关系列出方程用二分法解出芯块中心的温度。接下来依据冷却剂的温度，得出的控制体出口处的含汽量。进而依据W-3公式求出该出的临界热流量qDNB，最后得出该出的烧毁比DNBR。

3）热管的压降计算：热管的压降包括摩擦压降、提升压降、进出口局部压降、定位搁架出口压降。摩擦压降可由计算单相流的达西（Darcy）公式算得。提升压降可由根据位置的变化算得，其中参数都取平均值。其余的压降根据形阻压降的基本公式再乘以相应的系数求得。最后各项相加得出热管的总压降。2、流程图2、程序//3.1.cpp:定义控制台应用程序的入口点。//流体出口温度迭代计算#include"stdafx.h"#include〈math.h〉#include〈iostream〉usingnamespacestd;doubletfout(doublecp,doubleatfout){doubletfin=287,Fa=0.974,Nt=1820,W=32500,plxs=0.05,tfout;tfout=tfin+3600*Fa*Nt/(W*cp*(1—plxs));returntfout;}intmain(){doubleatfout,tfin=287,tav,cp,tfoutc=0;intpanding=1;printf("请输入假设出口温度\n");cin>>atfout;while(panding〉=1){tav=(atfout+tfin)/2;printf("请输入％f下的Cp\n",tav);cin>>cp;tfoutc=tfout(cp,atfout);if(fabs(tfoutc—atfout)〈二0.5)panding=0;atfout=tfoutc;}printf("出口温度为%f\n",tfoutc);while(true){}//3.3.1热管内流体温度计算#include"stdafx.h"#include〈math.h〉#include〈iostream〉usingnamespacestd;doubleffout(doubletfin,doubleatfout){doubleqav,fnr=1.33,feh=1.03,fehm=0.95,dcs=9.5,gyh,cp,wh,l=3.60,tfout;boolpanding=true;cout〈〈"请输入平均热流密度(W/m*m) "〈〈endl;cin>>qav;cout〈〈"请输入归一化参数"<<endl;cin>>gyh;cout〈〈"请输入热管单元通道流量Wh(t/h)"<<endl;cin>>wh;while(panding){doubletav=(atfout+tfin)/2;cout〈〈"请输入"〈〈tav〈〈"下的Cp(J/(kg.°C))"〈〈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;}returntfout;}intmain(){doubletfin,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"usingnamespacestd;intmain(){doubletfout[6]={291.69,301.38,314.86,327.70,334.89,338.22},gyhcs[6]={0.48,1.02,1.50,1.56,0.96,0.48};doublepr[6],drxs[6],dlnd[6],f1[6],f2[6],f[6];inti;for(i=0;i<6;i++){cout<<"请输入"〈〈tfout[i]<<"下的普朗特数"〈〈endl;cin>>pr[i];cout〈〈"请输入"〈〈tfout[i]〈〈"下的导热系数(W/(m・°C))"〈〈endl;cin>>drxs[i];cout〈〈"请输入"〈〈tfout[i]〈〈"下的动力粘度kg/(m・s)"〈〈endl;cin>>dlnd[i];f1[i]=2.25*pow(10,4)*gyhcs[i]*pow(dlnd[i],0.8)/(drxs[i]*pow(pr[i],0.4));f2[i]=346.38+1.83*pow(gyhcs[i],0.25)—tfout[i];if(f1[i]〉f2[i])f[i]=f2[i];elsef[i]=f1[i];cout〈〈"0f1="〈〈f1[i]〈〈"0f2="〈〈f2[i]〈〈"f="〈〈f[i]+tfout[i]〈〈endl;}while(true){}}//3.3.2.cpp:定义控制台应用程序的入口点。//计算内壁温度#include"stdafx.h"#include"iostream"#include"math.h"usingnamespacestd;doubleftci(doubletcs,doubleatci){doublefrn=1.33,feq=1.03,gyhcs,dcs=9.5,dci=8.6,kc,tav,q,tci;boolpanding=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;}returntci;}intmain(){doubletcs,tci,atci;cout<<"请输入外壁温度°C"<<endl;cin>>tcs;cout<<"请输入内壁假设温度C〃〈〈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〉usingnamespacestd;doublef1(doublet){doublefout;fout=38.24*log((t+402.55)/402.55)+4.788*pow(10,-13)*(pow((t+273.15),4)—pow(273.15,4));returnfout;}doublef2(doubletci,doubleatu){doubletu,fnr=1.33,feq=1.03,gyhcs,q,tmid,jieguo;boolpanding=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;}returntu;}intmain(){doubletci,atu,tu;cout<<"请输入燃料芯块外壁温度/°C〃〈〈endl;cin>>tci;cout<<"请输入假设二氧化铀中心温度/C"<<endl;cin>>atu;tu=f2(tci,atu);cout〈〈tu;while(true){}//3.5.cpp:定义控制台应用程序的入口点。//计算DNBR#include"stdafx.h"#include"iostream"#include"cmath"#definepi3.14#defineP15800000usingnamespacestd;doublefqDNB(doublehfout,doublehfin,doublehfs,doublehfg,doubleg,doubledlzj){doublex,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(T24*dlzj);e=0.8258+0.341*pow(10,-6)*(hfs-hfin);qDNB=a*b*c*d*e;returnqDNB;}intmain(){doublehfout[6],hfin[6],gyhcs[6],qDNB[6],DNBR[6],hfs,hfg,dcs=0.0095,g,s=0.0126,q,frn=1.33,fqe=1.03,dlzj;inti;dlzj=4*(s*s-pi*dcs*dcs/4)/(pi*dcs);for(i=0;i<6;i++){cout<<"请输入第"〈〈i+1<<"段控制体出口冷却剂比焓(J/kg)"<<endl;cin>>hfout[i];cout〈〈"请输入第"〈〈i+1〈〈"段进口温度下的冷却剂比焓(J/kg)"〈〈endl;cin>>hfin[i];cout〈〈"请输入第"〈〈i+1〈〈"段控制体归一化参数"〈〈endl;cin>>gyhcs[i];}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++){qDNB[i]=fqDNB(hfout[i],hfin[i],hfs,hfg,g,dlzj);DNBR[i]=qDNB[i]/(q*fqe*frn*gyhcs[i]);cout<<"第"〈〈i+1〈〈"段控制体qDNB为"〈〈qDNB[i]<<endl;cout<<"第"<<i+1<<"段控制体DNBR为"<<DNBR[i]<<endl;}while(true){}}//3.6.cpp:定义控制台应用程序的入口点。//计算堆芯压降#include"stdafx.h"#include"iostream"#include"cmath"usingnamespacestd;constdoubleg=9.8;constdoubleKo=1.0;constdoubleKi=0.75;constdoubleKgr=1.05;constdoubleL=3.60;constdoubleV=3.85;constdoubleG=9.88e+6/3600;constdoubledlzj=1.178e-2;constdoubletfout[6]={291.69,301.38,314.86,327.70,334.89,338.22};constdoubletfin[6]={287,291.69,301.38,314.86,327.70,334.89};constdoubletcs[6]={303.51,326.15,348.41,348.43,348.19,347.90};doublefpel(doublemidu) //计算提升压降{returnmidu*L*g/6;}doublefpf(doublef,doubledlzj,doublemidu)//计算摩擦压降{returnf*L*midu*V*V/(2*dlzj*6);}double邙玄(doublevi,doublevo) //计算加速压