核反应堆热工分析09..
核反应堆热工分析09
![核反应堆热工分析09](https://img.taocdn.com/s3/m/81cc72593c1ec5da50e270d0.png)
燃料元件 的形状、 它在堆芯 内的布置 方式以及 栅距允许 变化的范 围
二回路对 一回路冷 却剂热工 参数的要 求
冷却剂流 过堆芯的 流程以及 堆芯进口 处冷却剂 流量的分 配情况
2
9.1 热工设计准则
压水动力堆设计中所规定的稳态热工设计准则:
燃料元件芯块内 最高温度低于其 相应燃耗下的熔 化温度
提高冷却剂的工作压力 提高冷却剂的流量 适当选定堆冷却剂的工作温度
从提高堆芯的功率密度来降低电能成本 从增加核燃料的燃耗深度来降低电能成本 从减少核电厂的厂用电来降低电能成本
13
9.3.2 单通道模型反应堆热工设计的一般步骤和方法
商定有关热工参数 确定燃料元件参数
根据热工设计准则中规定的内容进 行有关的计算
堆芯热管最大焓升 hh,max F = = 堆芯名义最大焓升 hn,max
E H
堆芯热点最大热流密度 qh,max F = = 堆芯名义最大热流密度 qn,max
E q
7
综合考虑核和工程两方面的因素后,定义热流密度热点因子Fq 和焓升热管因子 FH 为
q n,m ax q h,m ax qh,max Fq=F F = = q q n,m ax q
9
9.2.4 热管因子和热点因子的应用
乘积法:就是指把所有最不利的工程偏差都同时集中作用 在热管或热点上。 混和法:把工程误差分为非随机误差和随机误差两大类,先 分别计算各类误差,最后再把它们综合起来。
下面就根据前述工程热管因子及热点因子的计算方法,先分别计算 各分因子的值而后再综合成总的工程热管因子及工程热点因子
1 i ,l ,el 2
p
1 i ,l , a 2
p
核反应堆热工分析整理
![核反应堆热工分析整理](https://img.taocdn.com/s3/m/5e27b2272af90242a895e531.png)
二.简答题1.影响功率分布的因素有哪些,分别有什么影响?答:A.燃料布置:均匀装载:中心区域会出现较高的功率峰值,限制堆的总功率输出量,且平均燃耗也较低分区装载:中心功率水平降低,外区功率水平上升,整体功率分布得到展平,平均燃耗也提高了B.控制棒:用控制棒时堆芯功率峰值对平均功率之比可能高于未受扰动的堆芯的该比值。
径向:堆寿命初期,中央控制棒插入可使径向功率分布得到展平,即中央部分中子通量及功率水平下降了,外区中子通量及功率水平提高了轴向:插入控制棒给轴向功率分布带来不利影响(如压水堆中,寿期初堆顶部插入控制棒,中子通量向堆底部歪斜;寿期末抽出控制棒,中子通量向堆顶部歪斜)C.水隙及空泡:水隙:引起附加慢化作用,使该处中子通量上升,导致周围元件功率升高,从而增大功率分布不均匀度;空泡:使周围的堆芯反应性下降2.控制棒、慢化剂和结构材料中热量产生的来源?答:A.控制棒的热源:吸收堆芯的γ辐射以及控制棒本身吸收中子的(n, α)或(n, γ)反应B.慢化剂的热源:裂变中子的慢化、吸收裂变产物放出的β粒子的一部分能量、吸收各种γ射线的能量C.结构材料的热源:几乎完全是由于吸收来自堆芯的各种γ射线3.两种沸腾的定义及特点是什么?答:大容积沸腾:定义:指由浸没在具有自由表面原来静止的大容积液体内的受热面所产生的沸腾特点:液体流速很低,自然对流换热起主导作用流动沸腾:定义:指流体流经加热通道是发生的沸腾,亦称为对流沸腾特点:常发生在强迫对流工况下4.沸腾临界包括哪两类,它们的物理特点及发生的区域分别是什么?答:5.答:A.(选取合适的燃料材料和包壳材料,并限制元件芯块中心温度低于燃料的熔点,包壳表面热流密度低于临界热流密度以)保证燃料元件包壳在寿命期内完整性B.棒径选择除满足物理设计中的水铀比要求外,还须满足热工传热的要求C.(限制包壳外表面的最大温度和限制芯块----包壳交界面处的最大温度来)保证在整个寿命期不产生不良的物理化学作用D.满足结构方面要求并易于加工,工艺性能好E.经济性好,价廉6.UO 2陶瓷燃料的优缺点?答:优点:熔点高、高温和高辐照下几何形状比较稳定;高温水和液态钠中具有良好的耐腐蚀性能;与包壳材料锆合金、不锈钢的相容性好缺点:导热性能差;热梯度下的脆性7.单相流体的流动压降包括哪些,定义分别是什么?答:提升压降:指流体自截面1至截面2时由流体位能改变而引起的压力变化加速压降:指因流体速度发生改变而引起的压力变化摩擦压降:指流体沿等截面直通道流动时由沿程摩擦阻力的作用而引起的压力损失 形阻压降:指流体流经有急剧变化的固体边界时所出现的集中压力损失8.垂直加热通道中流型有哪些,分别有什么特点?答:泡状流:液相是连续相,汽相以汽泡的形式弥散在液相中,两相同时沿通道流动。
核反应堆热工分析
![核反应堆热工分析](https://img.taocdn.com/s3/m/288f7d539b89680203d825cf.png)
核科学与技术学院
2 堆芯功率的分布及其影响因素
轻水作慢化剂的堆芯中,水隙的存在引起附加慢化作 用,使该处的中子通量上升,提高水隙周围元件的功 率,增大了功率分布的不均匀程度
影
响
燃料布置
功
率
分
控制棒
布
的
因 水隙及空泡
素
克服办法:采用棒束型控制棒组件
核科学与技术学院
2 堆芯功率的分布及其影响因素
轻水作慢化剂的堆芯中,水隙的存在引起附加慢化作 用,使该处的中子通量上升,提高水隙周围元件的功 率,增大了功率分布的不均匀程度
整个堆芯的 热功率
计入位于堆 芯之外的反 射层、热屏 蔽等的释热 量
R f N f 正比 Nc 1.60211010 Fa Ef N f Vc
qv Fa E f N f
Nt Nc / Fa (qvVc )106 / Fa
106 E f N f Vc
核科学与技术学院
2 堆芯功率的分布及其影响因素
核科学与技术学院
1 核裂变产生能量及其分布
裂变碎片的动能 约占总能量的 84%
裂变能的绝大部 分在燃料元件内 转换为热能,少 量在慢化剂内释 放,通常取 97.4% 在 燃 料 元 件内转换为热能
核科学与技术学院
1 核裂变产生能量及其分布
不同核素所释放出来的裂变能量是有差异的,一般认为取
Ef 200MeV
核科学与技术学院
2 堆芯功率的分布及其影响因素
非均匀堆栅阵
用具有等效截面的圆来代替原来的正方形栅元 假设热中子仅在整个慢化剂内均匀产生
运用扩散理论,燃料元件内热中子通量分布的表达式:
AI0 (K0r)
若燃料棒表面处的热中子通量为 ,则在 s s,则:
核反应堆热工分析-课程设计
![核反应堆热工分析-课程设计](https://img.taocdn.com/s3/m/0f2557aadd3383c4bb4cd249.png)
%另一种带入方式%P=15.8;Nt=1.82*10^9;W=8916.66667;t_fin=287;L=3.66 ;m=121;n=265;n0=17;d_cs=0.0095 ;d_ci=0.00836 ; %%有问题d_u=0.00819 ; %%有问题s=0.0126 ;kexi=0.05;Fu=0.974;F_rn=1.35;F_zn=1.528;F_ln=1.11;F_h5e=0.95;F_qe=1.03;F_he=1.085;K_in=0.75;K_out=1.0;K_gr=1.05;f1=0.48;f2=1.02;f3=1.5;f4=1.56;f5=0.96;f6=0.48;hg=5678;%t_fout=input('请任意输入一堆芯出口温度(T=325):t_fout='); t_fout=325;r=1;kk=0;while r>=0.0001kk=kk+1;t_p=0.5*(t_fout+t_fin);if (280<t_p<=300)cp=19.35*t_p-350;elseif (300<t_p<=320)cp=34.43*t_p-4874;elseif (320<t_p<=340)cp=84.7*t_p-20960.4;endt1=t_fin+Fu*Nt/(W*(1-kexi)*cp);r=(t1-t_fout)/t_fout;t_fout=t1;endt_p=0.5*(t_fout+t_fin);% 计算燃料棒表面平均热流密度fz=m*n*pi*d_cs*L ; %计算堆心燃料棒的总传热面积q=Fu*Nt/fz ; %燃料棒表面平均热流密度F_qn=F_rn*F_zn*F_ln; % 有问题!!q_max=q*F_qn*F_qe;ql=q*pi*d_cs ;ql_max=ql*F_qn*F_qe;%计算平均管的流速dert=0.0008;Af=m*n*(s^2-pi*d_cs^2/4)+4*n0*s*dert*m/2; %堆心内流通面积if(300<t_p<=310)v=4.25*10^(-6)*t_p+1.0002*10^(-4);elseif(310<t_p<=320)v=5.112*10^(-6)*t_p-1.672*10^(-4);elseif(320<t_p<=330)v=6.444*10^(-6)*t_p-5.9344*10^(-4);elseif(330<t_p<=340)v=9.0976*10^(-5)*t_p-1.469128*10^(-3);endV=W*(1-kexi)*v/(Af);%fprintf('V=%.4f\n',V);%计算热管流量Ab=s^2-pi*d_cs^2/4;Wh=W*(1-kexi)*Ab/Af;%fprintf('Wh=%.4f\n', Wh);%热管中的计算De=4*(s^2-pi*d_cs^2/4)/(pi*d_cs);%t_fh1=input('请任意输入第一控制体热管流体温度(T=291):t_fh1=');t_fh1=291;r1=1;kk1=0;while r1>=0.0001kk1=kk1+1;t_fhp1=0.5*(t_fh1+t_fin);if (t_fhp1<=300)cp1=19.35*t_fhp1-350;elseif (t_fhp1<=320)cp1=34.43*t_fhp1-4874;elseif (t_fhp1<=340)cp1=84.7*t_fhp1-20960.4;endt11=t_fin+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp1*6)*f1;r1=(t11-t_fh1)/t_fh1;t_fh1=t11;endt_fhp1=0.5*(t_fh1+t_fin);%fprintf('t_fh1=%.4f\n',t_fh1);ts=346.19;%计算第一控制体出口处的包壳外壁温if(t_fhp1<=300)u1=(-0.40072*t_fhp1+209.504)/10^6;k1=(-1.42504*t_fhp1+994.816)/1000;pr1=(-0.40072*t_fhp1+209.504)*(0.01935*t_fhp1-0.35)/(-1.42504*t_fhp1+994.816); elseif(t_fhp1<=320)u1=(-1.044416*t_fhp1+402.6128)/10^6;k1=(-6.916*t_fhp1+2642.104)/1000;pr1=(-1.044416*t_fhp1+402.6128)*(0.03443*t_fhp1-4.874)/(-6.916*t_fhp1+2642.104); elseif(t_fhp1<=340)u1=(-1.044416*t_fhp1+402.6128)/10^6;k1=(-6.916*t_fhp1+2642.104)/1000;pr1=(-1.044416*t_fhp1+402.6128)*(0.0847*t_fhp1-20.9604)/(-6.916*t_fhp1+2642.104); endre1=Wh*De/(Ab*u1) ;h1=0.023*re1^0.8*pr1^0.4*k1/De ;t_csh11=t_fh1+q*F_rn*f1*F_qe/h1;t_csh12=ts+25*(q*F_rn*f1*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh11<=t_csh12)t_csh1=t_csh11;elseif(t_csh11>=t_csh12)t_csh1=t_csh12;end%fprintf('t_csh1=%.3f\n',t_csh1);%t_csh1=t_csh1;%计算包壳内壁温度%t_cih1=input('请任意输入一包壳内壁温度(T=315):t_cih1='); t_cih1=315;r12=1;kk11=0;while r12>=0.0001kk11=kk11+1;t_cp1=0.5*(t_cih1+t_csh1);k12=0.00547*(1.8*t_cp1+32)+13.8;t12=t_csh1+ql*F_rn*F_qe*f1*log(d_cs/d_ci)/(2*pi*k12);r12=(t12-t_cih1)/t_cih1;t_cih1=t12;end%fprintf('t_cih1=%.4f\n',t_cih1);%计算芯块表面温度t_uh1=t_cih1+ql*F_rn*F_qe*f1/(pi*0.5*(d_ci+d_u)*hg);%fprintf('t_uh1=%.4f\n',t_uh1);%计算芯块中心温度if(t_uh1<=400) %300 400ku1=0.051*t_uh1+6.02;elseif(t_uh1<=500) %400 500ku1=0.0451*t_uh1+8.38;elseif(t_uh1<=600) %500 600ku1=0.0404*t_uh1+10.73elseif(t_uh1<=700) %600 700ku1=0.0368*t_uh1+12.89;elseif(t_uh1<=800) %700 800ku1=0.0337*t_uh1+15.06;elseif(t_uh1<=900) %800 900ku1=0.0312*t_uh1+17.06;elseif(t_uh1<=1000) %900 1000ku1=0.0292*t_uh1+18.86;elseif(t_uh1<=1100) %1000 1100ku1=0.0255*t_uh1+22.56;elseif(t_uh1<=1200) %1100 1200ku1=0.028*t_uh1+19.81;endkuto1=ku1+ql*F_rn*F_qe*f1/(4*pi*100) ;if(kuto1<=34.93) %500 600to1=22.173*kuto1-185.8 ;elseif(kuto1<=38.65) %600 700to1=26.74*kuto1-335.02;elseif(kuto1<=42.02) %700 800to1=27.17*kuto1-446.88;elseif(kuto1<=45.14) %800 900to1=32.05*kuto1-546.79;elseif(kuto1<=48.06) %900 1000to1=34.25*kuto1-645.89;elseif(kuto1<=50.61) %1000 1100to1=39.22*kuto1-884.71;elseif(kuto1<=53.41)% 1100 1200to1=35.71*kuto1-707.5;elseif(kuto1<=55.84) %1200 1298to1=40.32*kuto1-953.63;elseif(kuto1<=58.4) %1298 1405to1=41.84*kuto1-1037;elseif(kuto1<=61.95) %1405 1560to1=43.668*kuto1-1145;elseif(kuto1<=66.87) %1560 1738to1=36.23*kuto1-682.25;elseif(kuto1<=68.86) %1378 1876to1=69.35*kuto1-2899.2;elseif(kuto1<=71.31) %1876 1990to1=46.53*kuto1-1328.1;elseif(kuto1<=74.88) %1990 2155to1=46.22*kuto1-1305.8;elseif(kuto1<=79.16) %2155 2343to1=43.93*kuto1-1134;end%fprintf('to1=%.3f\n',to1);%热管中的计算%t_fh2=input('请任意输入第二控制体热管流体温度(T=300):t_fh2='); t_fh2=300;r2=1;kk2=0;while r2>=0.0001kk2=kk2+1;t_fhp2=0.5*(t_fh2+t_fh1);if (t_fhp2<=300)cp2=19.35*t_fhp2-350;elseif (t_fhp2<=320)cp2=34.43*t_fhp2-4874;elseif (t_fhp2<=340)cp2=84.7*t_fhp2-20960.4;endt22=t_fh1+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp2*6)*f2;r2=(t22-t_fh2)/t_fh2;t_fh2=t22;endt_fhp2=0.5*(t_fh2+t_fh1);%fprintf('t_fh2=%.4f\n', t_fh2);%计算第二控制体出口处的包壳外壁温if(t_fhp2<=300)u2=(-0.40072*t_fhp2+209.504)/10^6;k2=(-1.42504*t_fhp2+994.816)/1000;pr2=(-0.40072*t_fhp2+209.504)*(0.01935*t_fhp2-0.35)/(-1.42504*t_fhp2+994.816); elseif(t_fhp2<=320)u2=(-1.044416*t_fhp2+402.6128)/10^6;k2=(-6.916*t_fhp2+2642.104)/1000;pr2=(-1.044416*t_fhp2+402.6128)*(0.03443*t_fhp2-4.874)/(-6.916*t_fhp2+2642.104); elseif(t_fhp2<=340)u2=(-1.044416*t_fhp2+402.6128)/10^6;k2=(-6.916*t_fhp2+2642.104)/1000;pr2=(-1.044416*t_fhp2+402.6128)*(0.0847*t_fhp2-20.9604)/(-6.916*t_fhp2+2642.104); endre2=Wh*De/(Ab*u2);h2=0.023*re2^0.8*pr2^0.4*k2/De;t_csh21=t_fh2+q*F_rn*f2*F_qe/h2;t_csh22=ts+25*(q*F_rn*f2*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh21<=t_csh22)t_csh2=t_csh21;elseif(t_csh21>=t_csh22)t_csh2=t_csh22;end%fprintf('t_csh2=%.4f\n',t_csh2);%计算第二控制体包壳内壁温度%t_cih2=input('请任意输入第二控制体包壳内壁温度(例如T=350):t_cih1='); t_cih2=350;r22=1;kk22=0;while r22>=0.0001kk22=kk22+1;t_cp2=0.5*(t_cih2+t_csh2);k22=0.00547*(1.8*t_cp2+32)+13.8;t22=t_csh2+ql*F_rn*F_qe*f2*log(d_cs/d_ci)/(2*pi*k22);r22=(t22-t_cih2)/t_cih2;t_cih2=t22;end%fprintf('t_cih2=%.4f\n',t_cih2);%计算芯块表面温度t_uh2=t_cih2+ql*F_rn*F_qe*f2/(pi*0.5*(d_ci+d_u)*hg) ;;%fprintf('t_uh2=%.4f\n',t_uh2);%计算芯块中心温度if(t_uh2<=400) %300 400ku2=0.051*t_uh2+6.02;elseif(t_uh2<=500) %400 500ku2=0.0451*t_uh2+8.38;elseif(t_uh2<=600) %500 600ku2=0.0404*t_uh2+10.73elseif(t_uh2<=700) %600 700ku2=0.0368*t_uh2+12.89;elseif(t_uh2<=800) %700 800ku2=0.0337*t_uh2+15.06;elseif(t_uh2<=900) %800 900ku2=0.0312*t_uh2+17.06;elseif(t_uh2<=1000) %900 1000ku2=0.0292*t_uh2+18.86;elseif(t_uh2<=1100) %1000 1100ku2=0.0255*t_uh2+22.56;elseif(t_uh2<=1200) %1100 1200ku2=0.028*t_uh2+19.81;endkuto2=ku2+ql*F_rn*F_qe*f2/(4*pi*100) ;if(kuto2<=38.65) %600 700to2=26.74*kuto2-335.02;elseif(kuto2<=42.02) %700 800to2=27.17*kuto2-446.88;elseif(kuto2<=45.14) %800 900to2=32.05*kuto2-546.79;elseif(kuto2<=48.06) %900 1000to2=34.25*kuto2-645.89;elseif(kuto2<=50.61) %1000 1100to2=39.22*kuto2-884.71;elseif(kuto2<=53.41)% 1100 1200to2=35.71*kuto2-707.5;elseif(kuto2<=55.84) %1200 1298to2=40.32*kuto2-953.63;elseif(kuto2<=58.4) %1298 1405to2=41.84*kuto2-1037;elseif(kuto2<=61.95) %1405 1560to2=43.668*kuto2-1145;elseif(kuto2<=66.87) %1560 1738to2=36.23*kuto2-682.25;elseif(kuto2<=68.86) %1378 1876to2=69.35*kuto2-2899.2;elseif(kuto2<=71.31) %1876 1990to2=46.53*kuto2-1328.1;elseif(kuto2<=74.88) %1990 2155to2=46.22*kuto2-1305.8;elseif(kuto2<=79.16) %2155 2343to2=43.93*kuto2-1134;end%fprintf('to2=%.4f\n',to2);%热管中的计算%t_fh3=input('请任意输入第三控制体热管流体温度(T=314):t_fh3='); t_fh3=314;r3=1;kk3=0;while r3>=0.0001kk3=kk3+1;t_fhp3=0.5*(t_fh3+t_fh2);if (t_fhp3<=300)cp3=19.35*t_fhp3-350;elseif (t_fhp3<=320)cp3=34.43*t_fhp3-4874;elseif (t_fhp3<=340)cp3=84.7*t_fhp3-20960.4;endt33=t_fh2+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp3*6)*f3;r3=(t33-t_fh3)/t_fh3;t_fh3=t33;endt_fhp3=0.5*(t_fh3+t_fh2);%fprintf('t_fh3=%.3f\n', t_fh3);%计算第三控制体出口处的包壳外壁温if(t_fhp3<=300)u3=(-0.40072*t_fhp3+209.504)/10^6;k3=(-1.42504*t_fhp3+994.816)/1000;pr3=(-0.40072*t_fhp3+209.504)*(0.01935*t_fhp3-0.35)/(-1.42504*t_fhp3+994.816); elseif(t_fhp3<=320)u3=(-1.044416*t_fhp3+402.6128)/10^6;k3=(-6.916*t_fhp3+2642.104)/1000;pr3=(-1.044416*t_fhp3+402.6128)*(0.03443*t_fhp3-4.874)/(-6.916*t_fhp3+2642.104); elseif(t_fhp3<=340)u3=(-1.044416*t_fhp3+402.6128)/10^6;k3=(-6.916*t_fhp3+2642.104)/1000;pr3=(-1.044416*t_fhp3+402.6128)*(0.0847*t_fhp3-20.9604)/(-6.916*t_fhp3+2642.104); endre3=Wh*De/(Ab*u3);h3=0.023*re3^0.8*pr3^0.4*k3/De;t_csh31=t_fh3+q*F_rn*f3*F_qe/h3;t_csh32=ts+25*(q*F_rn*f3*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh31<=t_csh32)t_csh3=t_csh31;elseif(t_csh31>=t_csh32)t_csh3=t_csh32;end%fprintf('t_csh3=%.3f\n',t_csh3);%计算第三控制体包壳内壁温度%t_cih3=input('请任意输入第三控制体包壳内壁温度(例如T=385):t_cih3=');t_cih3=385;r33=1;kk33=0;while r33>=0.0001kk33=kk33+1;t_cp3=0.5*(t_cih3+t_csh3);k33=0.00547*(1.8*t_cp3+32)+13.8;t33=t_csh3+ql*F_rn*F_qe*f3*log(d_cs/d_ci)/(2*pi*k33); r33=(t33-t_cih3)/t_cih3;t_cih3=t33;end%fprintf('t_cih3=%.3f\n',t_cih3);%计算芯块表面温度t_uh3=t_cih3+ql*F_rn*F_qe*f3/(pi*0.5*(d_ci+d_u)*hg); %fprintf('t_uh3=%.3f\n',t_uh3);%计算芯块中心温度if(t_uh3<=400) %300 400ku3=0.051*t_uh3+6.02;elseif(t_uh3<=500) %400 500ku3=0.0451*t_uh3+8.38;elseif(t_uh3<=600) %500 600ku3=0.0404*t_uh3+10.73 ;elseif(t_uh3<=700) %600 700ku3=0.0368*t_uh3+12.89;elseif(t_uh3<=800) %700 800ku3=0.0337*t_uh3+15.06;elseif(t_uh3<=900) %800 900ku3=0.0312*t_uh3+17.06;elseif(t_uh3<=1000) %900 1000ku3=0.0292*t_uh3+18.86;elseif(t_uh3<=1100) %1000 1100ku3=0.0255*t_uh3+22.56;elseif(t_uh3<=1200) %1100 1200ku3=0.028*t_uh3+19.81;endkuto3=ku3+ql*F_rn*F_qe*f3/(4*pi*100) ;if kuto3<=38.65 %600 700to3=26.74*kuto3-335.02;elseif kuto3<=42.02 %700 800to3=27.17*kuto3-446.88;elseif(kuto3<=45.14) %800 900to3=32.05*kuto3-546.79;elseif kuto3<=48.06 %900 1000to3=34.25*kuto3-645.89;elseif kuto3<=50.61 %1000 1100to3=39.22*kuto3-884.71;elseif kuto3<=53.41 % 1100 1200to3=35.71*kuto3-707.5;elseif kuto3<=55.84 %1200 1298to3=40.32*kuto3-953.63;elseif kuto3<=58.4 %1298 1405to3=41.84*kuto3-1037;elseif kuto3<=61.95 %1405 1560to3=43.668*kuto3-1145;elseif kuto3<=66.87 %1560 1738to3=36.23*kuto3-682.25;elseif kuto3<=68.86 %1378 1876to3=69.35*kuto3-2899.2;elseif kuto3<=71.31 %1876 1990to3=46.53*kuto3-1328.1;elseif kuto3<=74.88 %1990 2155to3=46.22*kuto3-1305.8;elseif kuto3<=79.16 %2155 2343to3=43.93*kuto3-1134;end%fprintf('to3=%.3f\n',to3);%热管中的计算De=4*(s^2-pi*d_cs^2/4)/(pi*d_cs);%t_fh4=input('请任意输入第四控制体热管流体温度(T=326):t_fh4='); t_fh4=326;r4=1;kk4=0;while r4>=0.0001kk4=kk4+1;t_fhp4=0.5*(t_fh4+t_fh3);if t_fhp4<=300cp4=19.35*t_fhp4-350;elseif t_fhp4<=320cp4=34.43*t_fhp4-4874;elseif t_fhp4<=340cp4=84.7*t_fhp4-20960.4;endt44=t_fh3+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp4*6)*f4;r4=(t44-t_fh4)/t_fh4;t_fh4=t44;endt_fhp4=0.5*(t_fh4+t_fh3);%fprintf('t_fh4=%.3f\n', t_fh4);ts=346.19;%计算第四控制体出口处的包壳外壁温if(t_fhp4<=300)u4=(-0.40072*t_fhp4+209.504)/10^6;k4=(-1.42504*t_fhp4+994.816)/1000;pr4=(-0.40072*t_fhp4+209.504)*(0.01935*t_fhp4-0.35)/(-1.42504*t_fhp4+994.816); elseif(t_fhp4<=320)u4=(-1.044416*t_fhp4+402.6128)/10^6;k4=(-6.916*t_fhp4+2642.104)/1000;pr4=(-1.044416*t_fhp4+402.6128)*(0.03443*t_fhp4-4.874)/(-6.916*t_fhp4+2642.104); elseif(t_fhp4<=340)u4=(-1.044416*t_fhp4+402.6128)/10^6;k4=(-6.916*t_fhp4+2642.104)/1000;pr4=(-1.044416*t_fhp4+402.6128)*(0.0847*t_fhp4-20.9604)/(-6.916*t_fhp4+2642.104); endre4=Wh*De/(Ab*u4);h4=0.023*re4^0.8*pr4^0.4*k4/De;t_csh41=t_fh4+q*F_rn*f4*F_qe/h4;t_csh42=ts+25*(q*F_rn*f4*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh41<=t_csh42)t_csh4=t_csh41;elseif(t_csh41>=t_csh42)t_csh4=t_csh42;end%fprintf('t_csh4=%.3f\n',t_csh4);%计算第四控制体包壳内壁温度%t_cih4=input('请任意输入第五控制体包壳内壁温度(T=386):t_cih4=');t_cih4=386;r44=1;kk44=0;while r44>=0.0001kk44=kk44+1;t_cp4=0.5*(t_cih4+t_csh4);k44=0.00547*(1.8*t_cp4+32)+13.8;t44=t_csh4+ql*F_rn*F_qe*f4*log(d_cs/d_ci)/(2*pi*k44); r44=(t44-t_cih4)/t_cih4;t_cih4=t44;end%fprintf('t_cih4=%.4f\n',t_cih4);%计算芯块表面温度t_uh4=t_cih4+ql*F_rn*F_qe*f4/(pi*0.5*(d_ci+d_u)*hg); %fprintf('t_uh4=%.3f\n',t_uh4);%计算芯块中心温度if t_uh4<=400 %300 400ku4=0.051*t_uh4+6.02;elseif t_uh4<=500 %400 500ku4=0.0451*t_uh4+8.38;elseif t_uh4<=600 %500 600ku4=0.0404*t_uh4+10.73elseif t_uh4<=700 %600 700ku4=0.0368*t_uh4+12.89;elseif t_uh4<=800 %700 800ku4=0.0337*t_uh4+15.06;elseif t_uh4<=900 %800 900ku4=0.0312*t_uh4+17.06;elseif t_uh4<=1000 %900 1000ku4=0.0292*t_uh4+18.86;elseif t_uh4<=1100 %1000 1100ku4=0.0255*t_uh4+22.56;elseif t_uh4<=1200 %1100 1200ku4=0.028*t_uh4+19.81;endkuto4=ku4+ql*F_rn*F_qe*f4/(4*pi*100);if kuto4<=38.65 %600 700to4=26.74*kuto4-335.02;elseif kuto4<=42.02 %700 800to4=27.17*kuto4-446.88;elseif kuto4<=45.14 %800 900to4=32.05*kuto4-546.79;elseif kuto4<=48.06 %900 1000to4=34.25*kuto4-645.89;elseif kuto4<=50.61 %1000 1100to4=39.22*kuto4-884.71;elseif kuto4<=53.41 % 1100 1200to4=35.71*kuto4-707.5;elseif kuto4<=55.84 %1200 1298to4=40.32*kuto4-953.63;elseif kuto4<=58.4 %1298 1405to4=41.84*kuto4-1037;elseif kuto4<=61.95 %1405 1560to4=43.668*kuto4-1145;elseif kuto4<=66.87 %1560 1738to4=36.23*kuto4-682.25;elseif kuto4<=68.86 %1738 1876to4=69.35*kuto4-2899.2;elseif kuto4<=71.31 %1876 1990to4=46.53*kuto4-1328.1;elseif kuto4<=74.88 %1990 2155to4=46.22*kuto4-1305.8;elseif kuto4<=79.16 %2155 2343to4=43.93*kuto4-1134;end%fprintf('to4=%.3f\n',to4);%热管中的计算%t_fh5=input('请任意输入第五控制体热管流体温度(T=332):t_fh5='); t_fh5=332;r5=1;kk5=0;while r5>=0.0001kk5=kk5+1;t_fhp5=0.5*(t_fh5+t_fh4);if t_fhp5<=300cp5=19.35*t_fhp5-350;elseif t_fhp5<=320cp5=34.43*t_fhp5-4874;elseif t_fhp5<=340cp5=84.7*t_fhp5-20960.4;endt55=t_fh4+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp5*6)*f5;r5=(t55-t_fh5)/t_fh5;endt_fhp5=0.5*(t_fh5+t_fh5);%fprintf('t_fh5=%.4f\n', t_fh5);%计算第五控制体出口处的包壳外壁温if(t_fhp5<=300)u5=(-0.40072*t_fhp5+209.504)/10^6;k5=(-1.42504*t_fhp5+994.816)/1000;pr5=(-0.40072*t_fhp5+209.504)*(0.01935*t_fhp5-0.35)/(-1.42504*t_fhp5+994.816); elseif(t_fhp5<=320)u5=(-1.044416*t_fhp5+402.6128)/10^6;k5=(-6.916*t_fhp5+2642.104)/1000;pr5=(-1.044416*t_fhp5+402.6128)*(0.03443*t_fhp5-4.874)/(-6.916*t_fhp5+2642.104); elseif(t_fhp5<=340)u5=(-1.044416*t_fhp5+402.6128)/10^6;k5=(-6.916*t_fhp5+2642.104)/1000;pr5=(-1.044416*t_fhp5+402.6128)*(0.0847*t_fhp5-20.9604)/(-6.916*t_fhp5+2642.104); endre5=Wh*De/(Ab*u5);h5=0.023*re5^0.8*pr5^0.4*k5/De;t_csh51=t_fh5+q*F_rn*f5*F_qe/h5;t_csh52=ts+25*(q*F_rn*f5*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh51<=t_csh52)t_csh5=t_csh51;elseif(t_csh51>=t_csh52)t_csh5=t_csh52;end%fprintf('t_csh5=%.4f\n',t_csh5);%计算第五控制体包壳内壁温度%t_cih5=input('请任意输入第五控制体包壳内壁温度(T=371):t_cih5=');t_cih5=371;r55=1;kk55=0;while r55>=0.00001kk55=kk55+1;t_cp5=0.5*(t_cih5+t_csh5);k55=0.00547*(1.8*t_cp5+32)+13.8;t55=t_csh5+ql*F_rn*F_qe*f5*log(d_cs/d_ci)/(2*pi*k55);r55=(t55-t_cih5)/t_cih5;end%fprintf('t_cih5=%.4f\n',t_cih5);%计算芯块表面温度t_uh5=t_cih5+ql*F_rn*F_qe*f5/(pi*0.5*(d_ci+d_u)*hg); %fprintf('t_uh5=%.4f\n',t_uh5);%计算芯块中心温度if t_uh5<=400 %300 400ku5=0.051*t_uh5+6.02;elseif t_uh5<=500 %400 500ku5=0.0451*t_uh5+8.38;elseif t_uh5<=600 %500 600ku5=0.0404*t_uh5+10.73;elseif t_uh5<=700 %600 700ku5=0.0368*t_uh5+12.89;elseif t_uh5<=800 %700 800ku5=0.0337*t_uh5+15.06;elseif t_uh5<=900 %800 900ku5=0.0312*t_uh5+17.06;elseif t_uh5<=1000 %900 1000ku5=0.0292*t_uh5+18.86;elseif t_uh5<=1100 %1000 1100ku5=0.0255*t_uh5+22.56;elseif t_uh5<=1200 %1100 1200ku5=0.028*t_uh5+19.81;endkuto5=ku5+ql*F_rn*F_qe*f5/(4*pi*100);if kuto5<=38.65 %600 700to5=26.74*kuto5-335.02;elseif kuto5<=42.02 %700 800to5=27.17*kuto5-446.88;elseif kuto5<=45.14 %800 900to5=32.05*kuto5-546.79;elseif kuto5<=48.06 %900 1000to5=34.25*kuto5-645.89;elseif kuto5<=50.61 %1000 1100to5=39.22*kuto5-884.71;elseif kuto5<=53.41 % 1100 1200to5=35.71*kuto5-707.5;elseif kuto5<=55.84 %1200 1298to5=40.32*kuto5-953.63;elseif kuto5<=58.4 %1298 1405to5=41.84*kuto5-1037;elseif kuto5<=61.95 %1405 1560to5=43.668*kuto5-1145;elseif kuto5<=66.87 %1560 1738to5=36.23*kuto5-682.25;elseif kuto5<=68.86 %1738 1876to5=69.35*kuto5-2899.2;elseif kuto5<=71.31 %1876 1990to5=46.53*kuto5-1328.1;elseif kuto5<=74.88 %1990 2155to5=46.22*kuto5-1305.8;elseif kuto5<=79.16 %2155 2343to5=43.93*kuto5-1134;end%fprintf('to5=%.3f\n',to5);%热管中的计算%t_fh6=input('请任意输入第六控制体热管流体温度(T=336):t_fh6='); t_fh6=336;r6=1;kk6=0;while r6>=0.0001kk6=kk6+1;t_fhp6=0.5*(t_fh6+t_fh5);if t_fhp6<=300cp6=19.35*t_fhp6-350;elseif t_fhp6<=320cp6=34.43*t_fhp6-4874;elseif t_fhp6<=340cp6=84.7*t_fhp6-20960.4;endt66=t_fh5+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp6*6)*f6;r6=(t66-t_fh6)/t_fh6;t_fh6=t66;endt_fhp6=0.5*(t_fh6+t_fh6);%fprintf('t_fh6=%.4f\n', t_fh6);%计算第六控制体出口处的包壳外壁温if(t_fhp6<=300)u6=(-0.40072*t_fhp6+209.504)/10^6;k6=(-1.42504*t_fhp6+994.816)/1000;pr6=(-0.40072*t_fhp6+209.504)*(0.01935*t_fhp6-0.35)/(-1.42504*t_fhp6+994.816); elseif(t_fhp6<=320)u6=(-1.044416*t_fhp6+402.6128)/10^6;k6=(-6.916*t_fhp6+2642.104)/1000;pr6=(-1.044416*t_fhp6+402.6128)*(0.03443*t_fhp6-4.874)/(-6.916*t_fhp6+2642.104); elseif(t_fhp6<=340)u6=(-1.044416*t_fhp6+402.6128)/10^6;k6=(-6.916*t_fhp6+2642.104)/1000;pr6=(-1.044416*t_fhp6+402.6128)*(0.0847*t_fhp6-20.9604)/(-6.916*t_fhp6+2642.104); endre6=Wh*De/(Ab*u6);h6=0.023*re6^0.8*pr6^0.4*k6/De;t_csh61=t_fh6+q*F_rn*f6*F_qe/h6;t_csh62=ts+25*(q*F_rn*f6*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh61<=t_csh62)t_csh6=t_csh61;elseif(t_csh61>=t_csh62)t_csh6=t_csh62;end%fprintf('t_csh6=%.4f\n',t_csh6);%计算第六控制体包壳内壁温度%t_cih6=input('请任意输入第六控制体包壳内壁温度(T=358):t_cih6=');t_cih6=358;r66=1;kk66=0;while r66>=0.00001kk66=kk66+1;t_cp6=0.5*(t_cih6+t_csh6);k66=0.00547*(1.8*t_cp6+32)+13.8;t66=t_csh6+ql*F_rn*F_qe*f6*log(d_cs/d_ci)/(2*pi*k66);r66=(t66-t_cih6)/t_cih6;t_cih6=t66;end%fprintf('t_cih6=%.4f\n',t_cih6);%计算芯块表面温度t_uh6=t_cih6+ql*F_rn*F_qe*f6/(pi*0.5*(d_ci+d_u)*hg); %fprintf('t_uh6=%.4f\n',t_uh6);%计算芯块中心温度if t_uh6<=400 %300 400ku6=0.051*t_uh6+6.02;elseif t_uh6<=500 %400 500ku6=0.0451*t_uh6+8.38;elseif t_uh6<=600 %500 600ku6=0.0404*t_uh6+10.73elseif t_uh6<=700 %600 700ku6=0.0368*t_uh6+12.89;elseif t_uh6<=800 %700 800ku6=0.0337*t_uh6+15.06;elseif t_uh6<=900 %800 900ku6=0.0312*t_uh6+17.06;elseif t_uh6<=1000 %900 1000ku6=0.0292*t_uh6+18.86;elseif t_uh5<=1100 %1000 1100ku6=0.0255*t_uh6+22.56;elseif t_uh6<=1200 %1100 1200ku6=0.028*t_uh6+19.81;endkuto6=ku6+ql*F_rn*F_qe*f6/(4*pi*100);if kuto6<=38.65 %600 700to6=26.74*kuto6-335.02;elseif kuto6<=42.02 %700 800to6=27.17*kuto6-446.88;elseif kuto6<=45.14 %800 900to6=32.05*kuto6-546.79;elseif kuto6<=48.06 %900 1000to6=34.25*kuto6-645.89;elseif kuto6<=50.61 %1000 1100to6=39.22*kuto6-884.71;elseif kuto6<=53.41 % 1100 1200to6=35.71*kuto6-707.5;elseif kuto6<=55.84 %1200 1298to6=40.32*kuto6-953.63;elseif kuto6<=58.4 %1298 1405to6=41.84*kuto6-1037;elseif kuto6<=61.95 %1405 1560to6=43.668*kuto6-1145;elseif kuto6<=66.87 %1560 1738to6=36.23*kuto6-682.25;elseif kuto6<=68.86 %1738 1876to6=69.35*kuto6-2899.2;elseif kuto6<=71.31 %1876 1990to6=46.53*kuto6-1328.1;elseif kuto6<=74.88 %1990 2155to6=46.22*kuto6-1305.8;elseif kuto6<=79.16 %2155 2343to6=43.93*kuto6-1134;end%fprintf('to6=%.4f\n',to6);%热管中的计算f0=1/L;t_fh0=t_fin;ts=346.19;%计算第一控制体出口处的包壳外壁温if(t_fh0<=300)u0=(-0.40072*t_fh0+209.504)/10^6;k0=(-1.42504*t_fh0+994.816)/1000;pr0=(-0.40072*t_fh0+209.504)*(0.01935*t_fh0-0.35)/(-1.42504*t_fh0+994.816); elseif(t_fh0<=320)u0=(-1.044416*t_fh0+402.6128)/10^6;k0=(-6.916*t_fh0+2642.104)/1000;pr0=(-1.044416*t_fh0+402.6128)*(0.03443*t_fh0-4.874)/(-6.916*t_fh0+2642.104); elseif(t_fh0<=340)u0=(-1.044416*t_fh0+402.6128)/10^6;k0=(-6.916*t_fh0+2642.104)/1000;pr0=(-1.044416*t_fh0+402.6128)*(0.0847*t_fh0-20.9604)/(-6.916*t_fh0+2642.104); endre0=Wh*De/(Ab*u0) ;h0=0.023*re0^0.8*pr0^0.4*k0/De ;t_csh01=t_fh0+q*F_rn*f0*F_qe/h0;t_csh02=ts+25*(q*F_rn*f0*F_qe/10^6)^0.25*exp(-P/6.2);if(t_csh01<=t_csh02)t_csh0=t_csh01;elseif(t_csh01>=t_csh02)t_csh0=t_csh02;end%fprintf('t_csh0=%.3f\n',t_csh0);%计算包壳内壁温度t_cih0=300;r0=1;kk0=0;while r0>=0.0001kk0=kk0+1;t_cp0=0.5*(t_cih0+t_csh0);k02=0.00547*(1.8*t_cp0+32)+13.8;t02=t_csh0+ql*F_rn*F_qe*f0*log(d_cs/d_ci)/(2*pi*k02); r0=(t02-t_cih0)/t_cih0;t_cih0=t02;end%计算芯块表面温度t_uh0=t_cih0+ql*F_rn*F_qe*f0/(pi*0.5*(d_ci+d_u)*hg); %计算芯块中心温度if(t_uh0<=400) %300 400ku0=0.051*t_uh0+6.02;elseif(t_uh0<=500) %400 500ku0=0.0451*t_uh0+8.38;elseif(t_uh0<=600) %500 600ku0=0.0404*t_uh0+10.73elseif(t_uh0<=700) %600 700ku0=0.0368*t_uh0+12.89;elseif(t_uh0<=800) %700 800ku0=0.0337*t_uh0+15.06;elseif(t_uh0<=900) %800 900ku0=0.0312*t_uh0+17.06;elseif(t_uh0<=1000) %900 1000ku0=0.0292*t_uh0+18.86;elseif(t_uh0<=1100) %1000 1100ku0=0.0255*t_uh0+22.56;elseif(t_uh0<=1200) %1100 1200ku0=0.028*t_uh0+19.81;endkuto0=ku0+ql*F_rn*F_qe*f0/(4*pi*100) ;if(kuto0<=26.42) %300 400to0=19.61*kuto0-118.04;elseif(kuto0<=34.93) %500 600to0=22.173*kuto0-185.8 ;elseif(kuto0<=38.65) %600 700to0=26.74*kuto0-335.02;elseif(kuto0<=42.02) %700 800to0=27.17*kuto0-446.88;elseif(kuto0<=45.14) %800 900to0=32.05*kuto0-546.79;elseif(kuto0<=48.06) %900 1000to0=34.25*kuto0-645.89;elseif(kuto0<=50.61) %1000 1100to0=39.22*kuto0-884.71;elseif(kuto0<=53.41)% 1100 1200to0=35.71*kuto0-707.5;elseif(kuto0<=55.84) %1200 1298to0=40.32*kuto0-953.63;elseif(kuto0<=58.4) %1298 1405to0=41.84*kuto0-1037;elseif(kuto0<=61.95) %1405 1560to0=43.668*kuto0-1145;elseif(kuto0<=66.87) %1560 1738to0=36.23*kuto0-682.25;elseif(kuto0<=68.86) %1378 1876to0=69.35*kuto0-2899.2;elseif(kuto0<=71.31) %1876 1990to0=46.53*kuto0-1328.1;elseif(kuto0<=74.88) %1990 2155to0=46.22*kuto0-1305.8;elseif(kuto0<=79.16) %2155 2343to0=43.93*kuto0-1134;end%fprintf('to0=%.3f\n',to0);fprintf('q=%.4f\nq_max=%.4f\nql=%.4f\nql_max=%.4f\n ',q,q_max,ql,ql_max); fprintf('t_fout=%.4f\n',t_fout)fprintf('V=%.4f\n',V);fprintf('Wh=%.4f\n', Wh);fprintf('t_fh0=%.4f\nt_csh0=%.4f\nt_cih0=%.4ft_uh0=%.4f\n5_to0=%.4f\n',t_fh0,t_csh0,t_cih0,t_ uh0,to0);fprintf('t_fh1=%.4f\n',t_fh1);fprintf('t_csh1=%.4f\n',t_csh1);fprintf('t_cih1=%.4f\n',t_cih1);fprintf('t_uh1=%.4f\n',t_uh1);fprintf('to1=%.4f\n',to1);fprintf('t_fh2=%.4f\n',t_fh2);fprintf('t_csh2=%.4f\n',t_csh2);fprintf('t_cih2=%.4f\n',t_cih2);fprintf('t_uh2=%.4f\n',t_uh2);fprintf('to2=%.4f\n',to2);fprintf('t_fh3=%.4f\n',t_fh3);fprintf('t_csh3=%.4f\n',t_csh3);fprintf('t_cih3=%.4f\n',t_cih3);fprintf('t_uh3=%.4f\n',t_uh3);fprintf('to3=%.4f\n',to3);fprintf('t_fh4=%.4f\n',t_fh4);fprintf('t_csh4=%.4f\n',t_csh4);fprintf('t_cih4=%.4f\n',t_cih4);fprintf('t_uh4=%.4f\n',t_uh4);fprintf('to4=%.4f\n',to4);fprintf('t_fh5=%.4f\n',t_fh5);fprintf('t_csh5=%.4f\n',t_csh5);fprintf('t_cih5=%.4f\n',t_cih5);fprintf('t_uh5=%.4f\n',t_uh5);fprintf('to5=%.3f\n',to5);fprintf('t_fh6=%.4f\n',t_fh6);fprintf('t_csh6=%.4f\n',t_csh6);fprintf('t_cih6=%.4f\n',t_cih6);fprintf('t_uh6=%.4f\n',t_uh6);fprintf('to6=%.4f\n',to6);Hfs=1642.8*10^3;Hgs=2589.98*10^3;Hfg=Hgs-Hfs;G=W*(1-kexi)/Af;g=9.8;if (t_fin<=290)H0=(5.146*t_fin-208.22)*10^3;elseif (t_fin<=300)elseif (t_fin<=310)H0=(5.59*t_fin-339.44)*10^3;elseif (t_fin<=320)H0=(5.932*t_fin-445.46)*10^3;elseif (t_fin<=330)H0=(6.422*t_fin-602.26)*10^3;elseif (t_fin<=340)H0=(7.23*t_fin-868.9)*10^3;elseif (t_fin<=350)H0=(104.63*t_fin-33984.9)*10^3;endXe0=(H0-Hfs)/Hfg;q_DNB0=3.154*10^6*((2.022-6.238*10^(-8)*P)+(0.1722-1.43*10^(-8)*P)*exp((18.177-5.987*1 0^(-7)*P)*Xe0))*((0.1484-1.596*Xe0+0.1729*Xe0*abs(Xe0))*0.2049*G/10^6+1.037)*(1.157-0. 869*Xe0)*(0.2664+0.8357*exp(-124*De))*(0.8258+0.341*10^(-6)*(Hfs-H0));DNBR0=q_DNB0/(q*F_rn*F_qe*f0);%H1if (t_fh1<=290)H1=(5.146*t_fh1-208.22)*10^3;elseif (t_fh1<=300)H1=(5.344*t_fh1-265.64)*10^3;elseif (t_fh1<=310)H1=(5.59*t_fh1-339.44)*10^3;elseif (t_fh1<=320)H1=(5.932*t_fh1-445.46)*10^3;elseif (t_fh1<=330)H1=(6.422*t_fh1-602.26)*10^3;elseif (t_fh1<=340)H1=(7.23*t_fh1-868.9)*10^3;elseif (t_fh1<=350)H1=(104.63*t_fh1-33984.9)*10^3;endXe1=(H1-Hfs)/Hfg;q_DNB1=3.154*10^6*((2.022-6.238*10^(-8)*P)+(0.1722-1.43*10^(-8)*P)*exp((18.177-5.987*1 0^(-7)*P)*Xe1))*((0.1484-1.596*Xe1+0.1729*Xe1*abs(Xe1))*0.2049*G/10^6+1.037)*(1.157-0. 869*Xe1)*(0.2664+0.8357*exp(-124*De))*(0.8258+0.341*10^(-6)*(Hfs-H0));DNBR1=q_DNB1/(q*F_rn*F_qe*f1);if (t_fh2<=290)elseif (t_fh2<=300)H2=(5.344*t_fh2-265.64)*10^3;elseif (t_fh2<=310)H2=(5.59*t_fh2-339.44)*10^3;elseif (t_fh2<=320)H2=(5.932*t_fh2-445.46)*10^3;elseif (t_fh2<=330)H2=(6.422*t_fh2-602.26)*10^3;elseif (t_fh2<=340)H2=(7.23*t_fh2-868.9)*10^3;elseif (t_fh2<=350)H2=(104.63*t_fh2-33984.9)*10^3;endXe2=(H2-Hfs)/Hfg;q_DNB2=3.154*10^6*((2.022-6.238*10^(-8)*P)+(0.1722-1.43*10^(-8)*P)*exp((18.177-5.987*1 0^(-7)*P)*Xe2))*((0.1484-1.596*Xe2+0.1729*Xe2*abs(Xe2))*0.2049*G/10^6+1.037)*(1.157-0. 869*Xe2)*(0.2664+0.8357*exp(-124*De))*(0.8258+0.341*10^(-6)*(Hfs-H1));DNBR2=q_DNB2/(q*F_rn*F_qe*f2);if (t_fh3<=290)H3=(5.146*t_fh3-208.22)*10^3;elseif (t_fh3<=300)H3=(5.344*t_fh3-265.64)*10^3;elseif (t_fh3<=310)H3=(5.59*t_fh3-339.44)*10^3;elseif (t_fh3<=320)H3=(5.932*t_fh3-445.46)*10^3;elseif (t_fh3<=330)H3=(6.422*t_fh3-602.26)*10^3;elseif (t_fh3<=340)H3=(7.23*t_fh3-868.9)*10^3;elseif (t_fh3<=350)H3=(104.63*t_fh3-33984.9)*10^3;endXe3=(H3-Hfs)/Hfg;q_DNB3=3.154*10^6*((2.022-6.238*10^(-8)*P)+(0.1722-1.43*10^(-8)*P)*exp((18.177-5.987*1 0^(-7)*P)*Xe3))*((0.1484-1.596*Xe3+0.1729*Xe3*abs(Xe3))*0.2049*G/10^6+1.037)*(1.157-0. 869*Xe3)*(0.2664+0.8357*exp(-124*De))*(0.8258+0.341*10^(-6)*(Hfs-H2));DNBR3=q_DNB3/(q*F_rn*F_qe*f3);if (t_fh4<=290)H4=(5.146*t_fh4-208.22)*10^3;elseif (t_fh4<=300)。
核反应堆热工分析课程设计报告书详细过程版本
![核反应堆热工分析课程设计报告书详细过程版本](https://img.taocdn.com/s3/m/9a07c68249649b6648d747d1.png)
课程设计报告( 20 13 -- 2014 年度第二学期)名称:核反应堆热工分析课程设计题目:利用单通道模型进行反应堆稳态热工设计院系:核科学与工程学院班级:实践核1101班学号:1111440306学生:佳指导教师:王胜飞设计周数:1周成绩:日期:2014 年 6 月19 日一、课程设计的目的与要求反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。
对于反应堆热工设计,尤其是对动力堆,最基本的要安全。
要求在整个寿期能够长期稳定运行,并能适应启动、功率调节和停堆等功率变化,要保证在一般事故工况下堆芯不会遭到破坏,甚至在最严重的工况下,也要保证堆芯的放射性物质不扩散到周围环境中去。
在进行反应堆热工设计之前,首先要了解并确定的前提为:(1)根据所设计堆的用途和特殊要求(如尺寸、重量等的限制)选定堆型,确定所用的核燃料、冷却剂、慢化剂和结构材料等的种类;(2)反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化围;(3)燃料元件的形状、它在堆芯的分布方式以及栅距允许变化的围;(4)二回路对一回路冷却剂热工参数的要求;(5)冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。
在设计反应堆冷却系统时,为了保证反应堆运行安全可靠,针对不同的堆型,预先规定了热工设计必须遵守的要求,这些要求通常就称为堆的热工设计准则。
目前压水动力堆设计中所规定的稳态热工设计准则,一般有以下几点:(1)燃料元件芯块最高应低于其他相应燃耗下的熔化温度;(2)燃料元件外表面不允许发生沸腾临界;(3)必须保证正常运行工况下燃料元件和堆构件得到充分冷却;在事故工况下能提供足够的冷却剂以排除堆芯余热;(4)在稳态额定工况和可预计的瞬态运行工况中,不发生流动不稳定性。
在热工设计中,通常是通过平均通道(平均管)可以估算堆芯的总功率,而热通道(热管)则是堆芯中轴向功率最高的通道,通过它确定堆芯功率的上限,热点是堆芯中温度最高的点,代表堆芯热量密度最大的点,通过这个点来确定DNBR。
大学精品课件:核反应堆热工分析(压水堆结构概述)
![大学精品课件:核反应堆热工分析(压水堆结构概述)](https://img.taocdn.com/s3/m/6bde4ac61a37f111f0855b34.png)
• 喷淋系统:位于稳压器顶部,包
括主喷淋和辅助喷淋,用于减缓系 统热冲击、水温均匀及化学浓度、 降低系统压力;
• 电加热器:直接浸没的直套管式
电加热器,用于升高压力;
• 安全阀组:安装于稳压器顶部,
由保护阀与隔离阀组成;
• 测量仪表:主要用于水位检测与
显示;
核科学与技术学院
反应堆冷却剂泵
• 水力机械部分:泵体、热屏组件、
路具有放射性,管板与U形管属于冷却 剂压力边界;
• 排污与给水:防止各种杂质高度浓
缩以及一回路向二回路泄漏,确保正 常工况与特殊工况的给水要求;
• 水位控制及相关测量:水位测量
及调整、给水流量、蒸汽流量、蒸汽 压力等信号测量;
核科学与技术学院
压力壳——Mn-Mo-Ni低合金碳钢; 燃料——二氧化铀; 包壳——锆-4合金(Zr-4); 控制棒——银-铟-镉合金/316,304不锈钢(Ag-In-
temperature
• 120–400 MWe
• 15–30 year core life
• Cartridge core for regional fuel processing
(LFR)
Benefits
• Proliferation resistance of long-life cartridge core
英国建造32MWe原型堆,1976-1988年,运行的AGR共有14座, 8.9GW,由于受到CO2与不锈钢元件包壳材料化学相容性的限制 (690℃ ),使出口温度难以进一步提高,再加上功率密度低、燃耗低的 限制,使其仍难以和压水堆在经济上竞争;
• 高温气冷堆:采用90%以上的浓缩铀,全陶瓷燃料元件及堆芯,采
反应堆热工水力特性分析与设计研究
![反应堆热工水力特性分析与设计研究](https://img.taocdn.com/s3/m/2ec41c2fa55177232f60ddccda38376bae1fe064.png)
反应堆热工水力特性分析与设计研究引言反应堆热工水力特性是设计与研究反应堆核心的重要方面。
反应堆是一种利用核能进行能量转化和控制的设备,因此对其热工水力特性的分析与设计至关重要。
本文将对反应堆热工水力特性的分析与设计进行研究,并探讨其在核能利用过程中的重要性。
1. 反应堆热工水力特性分析1.1 反应堆热工水力循环反应堆热工水力循环是反应堆系统中热能转移的重要环节。
通过循环系统,热能可以在核燃料与冷却剂之间进行传递。
热工水力循环的设计应考虑冷却剂的流动和热能转移效率,以满足反应堆的运行需求。
常见的热工水力循环包括单相流循环和两相流循环。
1.2 反应堆热工水力特性分析方法反应堆热工水力特性的分析通常通过数值模拟和实验方法进行。
数值模拟可以通过计算流体力学(CFD)等方法来模拟反应堆内部的流动和热传导过程,以获得反应堆的热工水力特性参数。
实验方法可以通过搭建实验装置来观测和测量反应堆内部的流动和温度分布情况,以验证数值模拟结果的准确性。
2. 反应堆热工水力特性设计研究2.1 热工水力特性参数设计在反应堆的设计过程中,重要的一步是确定热工水力特性参数。
这些参数包括热流密度、冷却剂流速、冷却剂温度等。
热工水力特性参数的选择将直接影响反应堆的工作性能和安全性。
因此,需要通过理论分析和实验研究来确定这些参数的合理取值。
2.2 热工水力特性优化设计反应堆的热工水力特性优化设计旨在提高反应堆的热能转移效率和热功率密度,以提高反应堆的运行效率和能源利用效率。
优化设计可以通过改变反应堆的几何形状、流动通道的设计和材料选择等方法来实现。
通过优化设计,可以使反应堆具有更好的热工水力特性,提高反应堆的运行稳定性和安全性。
3. 反应堆热工水力特性在核能利用中的重要性反应堆热工水力特性在核能利用中起到至关重要的作用。
合理设计和控制反应堆的热工水力特性可以提高核能的利用效率和安全性。
同时,热工水力特性的分析与设计研究还可以为核能发电领域的技术创新和发展提供科学依据。
核反应堆热工分析课设
![核反应堆热工分析课设](https://img.taocdn.com/s3/m/b67df9156c85ec3a87c2c560.png)
目录一、设计任务 (1)二、课程设计要求 (2)三、计算过程 (2)四、程序设计框图 (8)五、代码说明书 (9)六、热工设计准则和出错矫正 (10)七、重要的核心程序代码 (11)八、计算结果及分析 (17)一、设计任务某压水反应堆的冷却剂及慢化剂都是水,用二氧化铀作燃料,用Zr-4作包壳材料。
燃料组件无盒壁,燃料元件为棒状,正方形排列。
已知下列参数:系统压力 15.8MPa堆芯输出功率 1820MW冷却剂总流量 32100t/h反应堆进口温度287℃堆芯高度 3.66m燃料组件数 121燃料组件形式17×17每个组件燃料棒数 265燃料包壳直径 9.5mm燃料包壳内径 8.36mm燃料包壳厚度 0.57mm燃料芯块直径 8.19mm燃料棒间距(栅距) 12.6mm芯块密度 95%理论密度旁流系数 5%燃料元件发热占总发热的份额 97.4%径向核热管因子 1.35轴向核热管因子 1.528局部峰核热管因子 1.11交混因子 0.95热流量工程热点因子 1.03焓升工程热管因子 1.085堆芯入口局部阻力系数 0.75堆芯出口局部阻力系数 1.0堆芯定位隔架局部阻力系数 1.05若将堆芯自上而下划分为5个控制体,则其轴向归一化功率分布如下表:堆芯轴向归一化功率分布(轴向等分5个控制体)通过计算,得出1. 堆芯出口温度;2. 燃料棒表面平均热流及最大热流密度,平均线功率,最大线功率;3. 热管的焓,包壳表面温度,芯块中心温度随轴向的分布;4. 包壳表面最高温度,芯块中心最高温度;5. DNBR在轴向上的变化;6. 计算堆芯压降;二、课程设计要求1.设计时间为两周;2.独立编制程序计算;3.迭代误差为0.1%;4.计算机绘图;5.设计报告写作认真,条理清楚,页面整洁;6.设计报告中要附源程序。
三、计算过程目前,压水核反应堆的稳态热工设计准则有:(1)燃料元件芯块内最高温度应低于其相应燃耗下的熔化温度。
核反应堆热工分析
![核反应堆热工分析](https://img.taocdn.com/s3/m/00b5813e83c4bb4cf7ecd12c.png)
运用扩散理论,燃料元件内热中子通量分布的表达式:
AI0 ( K0r )
若燃料棒表面处的热中子通量为
s ,则:
,则在
s
处, r
R0
I0 ( K0r ) s I 0 ( K 0 R0 )
2
堆芯功率的分布及其影响因素
补偿棒
补偿棒是用于抵消寿期初大量的 剩余反应性的
2
堆芯功率的分布及其影响因素
轻水作慢化剂的堆芯中,水隙的存在引起附加慢化 作用,使该处的中子通量上升,提高水隙周围元件的 功率,增大了功率分布的不均匀程度
影 响 功 率 分 布 的 因 素
燃料布置 控制棒 水隙及空泡
克服办法:采用棒束型控制棒组件
2
1
核裂变产生能量及其分布
不同核素所释放出来的裂变能量是有差异的,一般认为 取 E f 200MeV
堆内热源及其分布还与时间有关,新装料、平衡运行和停堆后都不 相同 输出燃料元件内产生的热量的热工水力问题就成为反应堆设计的关 键
2
堆芯功率的分布及其影响因素
释热率
堆内热源的分 布函数和中子 单位体积的 通量的分布函 数相同 释热率
堆芯最大体积释热率
qv,max Fa E f N f 0
2
堆芯功率的分布及其影响因素
均匀裸堆中的中子通量分布
2
堆芯功率的分布及其影响因素
均匀装载燃料方案: 早期的压水堆采用此方案 优点:装卸料方便 缺点:功率分布过于不平均,平均燃耗低
燃料布置 控制棒 水隙及空泡
影 响 功 率 分 布 的 因 素
(r , z ) 0 J 0 (2.405
外推半径:R e
核反应堆热工分析
![核反应堆热工分析](https://img.taocdn.com/s3/m/87711b34bb1aa8114431b90d6c85ec3a87c28b8c.png)
核反应堆热工分析核反应堆热工分析是指对核反应堆的热力学和动力学性质进行评估和分析。
在核反应堆的设计和操作过程中,热工分析是至关重要的一步。
核反应堆是利用核反应的原理来产生大量能量的一种设备。
核反应堆通常用于发电、生产同位素、研究等领域。
核反应堆的基本工作原理是将一些放射性物质放入到反应堆中,在核反应的过程中释放出大量的热能,以此来驱动涡轮机发电。
核反应堆的热工分析主要关注的是反应堆内的热力学和动力学特性。
在核反应堆中,核燃料的裂变会产生大量的热,这些热需要通过反应堆内的冷却剂来传递到发电设备。
因此,热工分析的主要任务是评估反应堆中热量的产生和传递,以便在设计和操作过程中避免过热或过冷的问题。
核反应堆的热工分析可以通过多种方法进行,其中最常见的方法是数值模拟。
数值模拟是指使用计算机模拟反应堆内复杂的热力学和动力学过程,以便更好地理解反应堆内的热传递和热力学行为。
在数值模拟中,需要考虑的因素包括反应堆内燃料的构成和排布、冷却剂的流动和热传递、反应堆的几何形状等。
另外,通过实验来验证和修正数值模拟模型也是非常重要的一项工作。
实验可以获得反应堆内的温度、压力、流速等关键参数,以进一步改进数值模拟模型。
同时,在实际操作中对反应堆进行连续监测和评估也是必不可少的。
核反应堆热工分析的结果可以帮助反应堆设计师和操作员更好地了解反应堆内的热力学行为,并且可以预测反应堆在不同工作条件下的热传递行为。
这些分析结果可以用于优化反应堆的设计,提高反应堆的效率和安全性。
总之,核反应堆热工分析是核反应堆设计和操作中不可或缺的一步。
通过热工分析可以更好地了解反应堆内的热力学和动力学行为,预测反应堆在不同工作条件下的性能,提高反应堆的效率和安全性。
除了数值模拟和实验外,还有一些重要的因素需要考虑,这些因素包括:1. 燃料特性燃料的特性直接影响反应堆内的核反应过程,从而影响热传递效率。
例如,燃料的粒度和密度会影响其热传递特性,而燃料的化学成分和裂变产物的特性也会影响其热放射。
[精选]第七章核反应堆热工--资料
![[精选]第七章核反应堆热工--资料](https://img.taocdn.com/s3/m/73aee223844769eae009eda8.png)
2.2.2、控制棒对功率分布的影响(2)
2.2.3、水隙及空泡对功率分布的影响
2.3、燃料元件内的功率分布(1)
2.3、燃料元件内的功率分布(2)
2.4、核热管因子(1)
热管和热点的概念
2.4、核热管因子(2)
热管因子:
为了衡量各有关的热工参数的最大值偏离平均值 (或名义值)的程度,引入一个修正因子,这个修正因 子就称为热管因子。热管因子是用各有关的热工(或物 理)参数的最大值与平均值的比值来表示的。
3.1.1.2.3、沸腾放热(4)
过冷沸腾起始点的判据:
9 qONB 1.798 103 p 1.156 t w t s 5 t w t s 和系统压力为p时, 式中qONB 为在壁面过冷度 开始产生沸腾所需的热 流量
2.828 p 0.0234
沸腾临界:
1.1、核燃料(4)
固体核燃料: 陶瓷燃料:氧化物、碳化物、氮化物
氧化铀:特点(5点内容)(自修) 热物性(熔点、密度、热导率、比热)(自修) 钚、铀混合物:UO2+PuO2 UC+PuC UN+PuN
陶瓷型燃料颗粒均匀分布在非裂变材料的基体中。 基体材料:铝、不锈钢、锆合金、石墨等 缺点:基体材料所占百分比大,必须使用浓缩铀(加浓铀)
缺点:沸点低、存在沸腾临界、在高温下有腐蚀作用。
重水慢化堆采用重水作冷却剂的好处是可以减少核燃料的装载量或降低核燃料的浓 缩度。缺点是价格昂贵。
钠:钠作为冷却剂主要应用于快中子堆。
缺点:钠水剧烈反应、温度梯度质量迁移、金属的扩散结合、存在由反应性正空泡 效应引起的控制和安全问题。
氦气:氦气作为冷却剂主要应用于气冷堆。
核反应堆热工分析(热工部分)
![核反应堆热工分析(热工部分)](https://img.taocdn.com/s3/m/25cb5b44cf84b9d528ea7ae6.png)
影 响 功 率 分 布 的 因 素
燃料布置 控制棒 水隙及空泡
通常I区的燃料富集度是最低的,III区的燃料富集度最高
核科学与技术学院
二
2.堆芯功率的分布及其影响因素
控制棒一般均匀布置在高中子通量的区域,既提高 控制棒的效率,又有利于径向中子通量的展平
影 响 功 率 分 布 的 因 素
燃料布置 控制棒 水隙及空泡
核科学与技术学院
四
停堆后的功率
在反应堆停堆后,由于中子在很短一段时间内还会引起裂变,裂变产物的 衰变以及中子俘获产物的衰变还会持续很长时间,因而堆芯仍有一定的释 热率。这种现象称为停堆后的释热,与此相应的功率称为停堆后的剩余功 率。
核科学与技术学院
四
停堆后的功率
剩余中子引起的裂变 燃料棒内储存的显热
(r , z ) 0 J 0 (2.405
外推半径:R e
r z ) cos Re LRe
R R R 0.71tr 外推高度: LRe LR 2LR LR 1.42tr
堆芯的释热率分布
r z qv (r , z ) qv ,max J 0 (2.405 ) cos Re LRe
二
1.核裂变4% 裂变能的绝大部 分在燃料元件内 转换为热能,少 量在慢化剂内释 放 , 通 常 取 97.4% 在 燃 料 元 件内转换为热能
核科学与技术学院
二
1.核裂变产生能量及其分布
不同核素所释放出来的裂变能量是有差异的,一般认为取 E f 200 MeV
核科学与技术学院
反应堆热工水力学
堆内释热
核科学与技术学院
1 一
1.核反应堆热工分析的任务
安全:
稳定运行,能 适应瞬态稳态 变化,且保证 在一般事故工 况下堆芯不会 破坏,最严重 事故工况下也 要保证堆芯放 射性不泄漏
《核反应堆热工分析》复习资料大全
![《核反应堆热工分析》复习资料大全](https://img.taocdn.com/s3/m/d55c4d09bb1aa8114431b90d6c85ec3a87c28bde.png)
《核反应堆热工分析》复习资料大全1. 核反应堆分类:按中子能谱分快中子堆、热中子堆按冷却剂分轻水堆(压水堆,沸水堆)、重水堆、气冷堆、钠冷堆按用途分研究试验堆:研究中子特性、生产堆: 生产易裂变材料、动力堆:发电舰船推进动力2.各种反应堆的差不多特点:3.压水堆优缺点:4.沸水堆与压水堆相比有两个优点:第一是省掉了一个回路,因而不再需要昂贵的蒸汽发生器。
第二是工作压力能够降低。
为了获得与压水堆同样的蒸汽温度,沸水堆只需加压到约72个大气压,比压水堆低了一倍。
5.沸水堆的优缺点:6.重水堆优缺点:优点:●中子利用率高〔要紧由于D吸取中子截面远低于H〕●废料中含235U极低,废料易处理●可将238U 转换成易裂变材料238U + n →239Pu239Pu + n →A+B+n+Q(占能量一半)缺点:●重水初装量大,价格昂贵●燃耗线〔8000~10000兆瓦日/T〔铀〕为压水堆1/3〕●为减少一回路泄漏〔因补D2O昂贵〕对一回路设备要求高7.高温气冷堆的优缺点:优点:●高温,高效率〔750~850℃,热效率40%〕●高转换比,高热耗值〔由于堆芯中没有金属结构材料只有核燃料和石墨,而石墨吸取中子截面小。
转换比0.85,燃耗10万兆瓦日/T〔铀〕〕●安全性高〔反应堆负温度系数大,堆芯热容量大,温度上升缓慢,采取安全措施裕量大〕●环境污染小〔采纳氦气作冷却剂,一回路放射性剂量较低,由于热孝率高排出废热少〕●有综合利用的宽敞前景〔假如进一步提高氦气温度~900℃时可直截了当推动气轮机;~1000℃时可直截了当推动气轮机热热效率大于50%;~1000-1200℃时可直截了当用于炼铁、化工及煤的气化〕●高温氦气技术可为今后进展气冷堆和聚变堆制造条件8.钠冷快堆的优缺点:优点:●充分利用铀资源239Pu + n →A+B+2.6个n238U + 1.6个n →1.6个239Pu 〔消耗一个中子使1.6个238U 转换成239Pu 〕●堆芯无慢化材料、结构材料,冷却剂用量少●液态金属钠沸点为895℃堆出口温度可高于560 ℃缺点:●快中子裂变截面小,需用高浓铀〔达~33%〕●对冷却剂要求苛刻,既要传热好又不能慢化中子,Na是首选材料,Na是爽朗金属,遇水会发生剧烈化学反应,因此需要加隔水回路9.各种堆型的特点、典型运行参数第二章堆芯材料选择和热物性〔简答〕1.固体核燃料的5点性能要求:教材14页2.常见的核燃料:金属铀和铀合金、陶瓷燃料、弥散体燃料3.选择包壳材料,必须综合考虑的7个因素:包壳材料的选择•中子吸取截面要小•热导率要大•材料相容性要好•抗腐蚀性能 •材料的加工性能 •材料的机械性能 •材料的抗辐照性能只有专门少的材料适合制作燃料包壳,铝、镁、锆、不锈钢、镍基合金、石墨。
核反应堆热工分析答案2009
![核反应堆热工分析答案2009](https://img.taocdn.com/s3/m/f8fc98d36137ee06eef9180d.png)
一、 简答题 1. 压水堆的稳态热工设计准则有: (1) 、燃料元件芯块内最高温度应低于其熔化温度。1 分 (2) 、燃料元件外表面不允许发生沸腾临界。1 分 (3) 、在稳态额定工况和可预计的瞬态工况下,不发生流动不稳定性。2 分 (4) 、必须保证正常运行工况下燃料元件和堆内构件能得到充分冷却;在事故工况下能提 供足够的冷却剂以排出堆芯余热。2 分 2. 气液两相流的流量漂移静态不稳定性产生的原因是:气液两相流中,压降随流量变化 的趋势由两个因素决定:(1)、由于流量降低压降有下降趋势;(2)、由于产生沸 腾,汽水混合物体积膨胀流速增加,从而使压降反而随流量的减少而增大。当第二个 因数其主要作用时,压降随流量的变化趋势成 N 字形(如下图示),即在两相区域存 在负斜率区。在 2’、3’两点之间所包含的压降范围内一个压降可能对应三个不同的 流量。由于水动力曲线的这样的变化趋势,当提供一个外加驱动压头时通道中的流量 就可能在三个不同流量下中漂移,这就发生了流量漂移静态不稳定性。 4 分
当燃料的线功率密度提升 10%时,燃料的最高许用温度应提高到: o o 因为燃料芯块中心的最高温度为 2200 C,点处燃料外表面温度为 500 C。
to
由
tu
k (t )dt 4 则
u
ql
q l 4 (
to
2200
0
'
k u (t)dt
500
0
k u (t)dt) q l 564.82W / cm
4分 当反应堆功率增加时,堆芯热通道的实际热流密度曲线向上移仍为余弦函数分布,临界热流密 度线向下移。3 分 三、(10 分) 解:自然循环:自然循环是指在闭合回路内依靠热段(向上流)和冷段(向下流)中的流体密度差 所产生的驱动压头来实现的流动循环。4 分 驱动压头: 由于将反应堆看作加热点源, 蒸汽发生器看作冷却点源, 热源冷源之间的高度为 L。
核反应堆热工分析
![核反应堆热工分析](https://img.taocdn.com/s3/m/00b5813e83c4bb4cf7ecd12c.png)
影 响 功 率 分 布 的 因 素
燃料布置 控制棒 水隙及空泡
控制棒对轴向功率分布的影响
2
堆芯功率的分布及其影响因素
停堆棒通常在堆芯的外面,只有 在需要停堆的时候才迅速插入堆 芯
停堆棒
影 响 功 率 分 布 的 因 素 燃料布置 控制棒 水隙及空泡
分 类
调节棒
调节棒是用于反应堆正常运行时 功率的调节
2
堆芯功率的分布及其影响因素
燃料元件数很多的非均匀圆柱形堆芯的通量分布总趋势与均匀堆 的是一样的 非均匀堆中的燃料元件自屏效应,使得元件内的中子通量和它周 围慢化剂内的中子通量分布会有较大差异
2
堆芯功率的分布及其影响因素
非均匀堆栅阵
用具有等效截面的圆来代替原来的正方形栅元 假设热中子仅在整个慢化剂内均匀产生
2.19 103 c(1 ) exp(3.14 106 )
若是低富集度铀作为燃料的压水堆,可取c=0.6,a=0.2 上式忽略了其他俘获产物对衰变功率的贡献,通常间计算结果再乘以 系数1.1
第三章
研究目的:
堆的传热过程
在保证反应堆安全的前提下,尽可能地提高堆芯单位 体积的热功率、冷却剂的温度等,以提高核动力的经 济性
裂变率
热功率
热功率
计入位于堆 芯之外的反 射层、热屏 蔽等的释热 量
单位时间, 单位体积燃 料内,发生 的裂变次数
整个堆芯的 热功率
R f N f
正比
Nc 1.60211010 Fa E f N f Vc
Nt N c / Fa (qvVc )106 / Fa 106 E f N f Vc
与冷却剂之间直接接触时的热交换 换热过程所传递的热量可用牛顿冷却 定律来求得,即: Q hF f
哈工程核反应堆热工分析简答
![哈工程核反应堆热工分析简答](https://img.taocdn.com/s3/m/a4ec94d2d15abe23482f4d73.png)
1,堆内热源的由来和分布特点。
裂变(1)瞬发裂变碎片的动能(在燃料元件内);裂变中子的动能(大部分在慢化剂中);瞬发γ射线的能量(堆内各处)。
(2)缓发裂变产物衰变的β射线能(大部分在燃料元件内);裂变产物衰变的γ射线(堆内各处)过剩中子引起的(n,γ)反应瞬发和缓发来源是过剩中子引起的裂变反应加(n,γ)反应产物的β衰变和γ衰变能(堆内各处)。
2,影响堆芯功率分布的因素主要有哪些?(1)燃料布置,目前核电厂压水堆通常采用分压装载的方案。
优点,功率分布得到了展平;燃料的平均燃耗提高了。
(2)控制棒,均匀布置有利于径向中子通量的展平,但给轴向功率分布带来了不利影响。
(3)水隙及空泡,水隙引起的附加慢化作用使该处的中子通量上升因而使水隙周围远见的功率升高,从而增大了功率分布的不均匀程度。
克服方法,采用棒束控制棒组件。
空泡将会导致堆芯反应性下降,空泡的存在能减轻某些事故的严重性。
沸水堆,下部插入。
4,燃料的自屏效应。
5,堆内结构材料3,控制棒中的热源来源是什么?(1)吸收堆芯的r辐射。
(2)控制棒本身吸收的中子的(n,α)或者(n,γ)反应。
4,热中子反应堆中慢化剂中的热源来源是什么?(1)裂变中子的慢化。
(2)吸收裂变产物放出的β粒子的一部分能量。
(3)吸收各种γ射线的能量。
5,反应堆停堆后的功率由哪几部分组成?有何特点。
(1)燃料棒内储存的显热。
(2)剩余中子引起的裂变。
(3)裂变产物的衰变以及中子俘获产物的衰变。
特点,在反应堆停堆后,其功率并不是立刻降为零,而是按照一个负的周期迅速的衰减,周期的长短最终取决于寿命最长的放射缓发中子的裂变核群的半衰期。
6,铀-235作为燃料的压水堆,每次裂变释放出来的总能量约为多少?在大型压水堆的设计中,往往取燃料元件的释热量占堆总释热量的百分之几?200Mev97.4%7,与早期压水堆中采用的均匀装载方案相比,现代大型压水堆采用分区装载方案的优点是什么?(1)功率的分布得到了展平,这对提高整个反应堆的热功率都是有利的。
核反应堆热工分析课程设计
![核反应堆热工分析课程设计](https://img.taocdn.com/s3/m/35eacd39aaea998fcc220e54.png)
2.计算堆芯出口温度(平均管)℃出口温度℃进出口温度平均值℃313300 325.2113135306.1056567324.1126955305.5563478324.2170398305.6085199324.2071793计算堆芯出口温度(平均管)3.燃料元件表面热流密度3.燃料元件表面热流密度堆芯找总传热面积F总表面平均热流密度q平均燃料棒表面最大热流密度q max平均线功率ql平均最大线功率ql,max4.平均管情况堆芯总流通面积Af平均密度ρf平均管的流速V5.热管流量通道截面积Ab热管通道内流量Wh6.热管中的计算6.1热管中的流体温度与出口包壳外壁温度π的值取: 3.141592654设中间变量K,有Tfh=Tfin+K*i/Cp则K的值是:53.20947968入口的温度T0287控制体10.49出口温度℃进出口温度平均值℃305296 291.8677009289.4338505292.020934289.510467292.020*******.5102843控制体2 1.04出口温度℃进出口温度平均值℃300296.0102843 302.4264044297.2234865302.300162297.1603653302.3028651297.1617168控制体3 1.5出口温度℃进出口温度平均值℃315308.6514325 316.3405646309.3217149316.2884384309.2956518316.2904793309.2966722控制体4 1.52出口温度℃进出口温度平均值℃330323.1452396329.1687983322.7296388329.2146074322.7525433329.2120963322.7512878控制体50.96出口温度℃进出口温度平均值℃340334.6060482 336.3835807332.7978385336.561608332.8868521336.5531338332.8826151控制体60.49出口温度℃进出口温度平均值℃350343.2765669 339.6501479338.1016408340.0167062338.28492340.0054583338.279296340.005805338.27946946.2出口包壳内壁温度和UO2外表面温度,UO2中心温度。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
堆芯稳态热工分析
反应堆热工设计的任务:
就是要设计一个既安全可靠而又经济的堆芯输热系统。 在进行反应堆热工设计之前,由各专业讨论并初步确定的前提为:
根据所设 计堆的用 途和特殊 要求选定 堆型,确 定所用的 核燃料、 慢化剂、 冷却剂和 结构材料 的种类
反应堆的 热功率、 堆芯功率 分布不均 匀系数和 水铀比允 许的变化 范围
9
9.2.4 热管因子和热点因子的应用
乘积法:就是指把所有最不利的工程偏差都同时集中作用 在热管或热点上。 混和法:把工程误差分为非随机误差和随机误差两大类,先 分别计算各类误差,最后再把它们综合起来。
下面就根据前述工程热管因子及热点因子的计算方法,先分别计算 各分因子的值而后再综合成总的工程热管因子及工程热点因子
F
N H
hmax 热管最大焓升 = 堆芯平均管焓升 h
6
9.2.2 工程热管因子和工程热点因子
由于工程上不可避免的误差,会使堆芯内燃料元件的热流密度、 冷却剂流量、冷却剂焓升以及燃料元件的温度等偏离名义值。 定量分析由工程因素引起的热工参数偏离名义值的程度, 热流密度工程热点因子 FqE 和焓升工程热管因子 FE H
热管和热点对确定堆芯功率的输出量起着决定性的作用。
5
堆芯功率分布的不均匀程度常用热流密度核热点因FqN 来 表示,有时也称热流密度核热管因子。
堆芯最大热流密度 qmax N N F = FR FZ 堆芯平均热流密度 q
N q
热管和平均管中冷却剂焓升的比值,称为焓升核热管因子。 用 FN H 表示。
10
热流密度工程热点因子 FqE
FqE
qh,max qn,max
qn,max q qn,max
1
q 1 3( ) qn,max qn,max
qE
焓升工程热管因子
由燃料芯块直径、密度及裂变物质富集度的加工误差所引起的焓升工程热管因子
FE H ,1 1 3(
提高冷却剂的工作压力 提高冷却剂的流量 适当选定堆冷却剂的工作温度
燃料元件外表面 不允许发生沸腾 临界
必须保证正常运 行工况下燃料元 件和堆内构件能 得到充分冷却; 在事故工况下能 提供足够的冷却 剂以排出堆芯余 热
在稳态工况下和 可预计的瞬态运 行工况中,不发 生流动不稳定性
3Leabharlann 燃料元件外表面不允许发生沸腾临界,通常用临界热流密 度比(DNBR)来定量地表示这个限制条件。 DNBR的定义是:
[
2 (4 P dcs , n ) /( d cs , n ) 2 h ,min
2
d
2 cs , h ,max
) /( dcs ,h,max )
]
1 b 2 b
11
堆芯下腔室冷却剂流量分配不均匀的焓升工程热管分因子
F
E H ,3
Qn,max / Wh,min,3 Qn,max / W
W W h,min,3
Wh,min,3 Wh,min,4
考虑热管内冷却剂流量再分配时的焓升工程热管分因子
F
E H ,4
hh,max,4 hn,max,3
Qn,max / Wh,min,4 Qn,max / Wh,min,3
考虑相邻通道冷却剂间相互交混的焓升工程热管分因子
F
E H ,5
N q E q
hn,m ax hh,m ax hh,max FH=F F = = h hn,m ax h
N H E H
同时考虑核和工程两方面因素后,热管和热点的定义为: 热管是堆芯内具有最大焓升的冷却剂通道。
热点是燃料元件上限制堆芯功率输出的局部点。
8
9.2.3 降低热管因子和热点因子的途径
燃料元件 的形状、 它在堆芯 内的布置 方式以及 栅距允许 变化的范 围
二回路对 一回路冷 却剂热工 参数的要 求
冷却剂流 过堆芯的 流程以及 堆芯进口 处冷却剂 流量的分 配情况
2
9.1 热工设计准则
压水动力堆设计中所规定的稳态热工设计准则:
燃料元件芯块内 最高温度低于其 相应燃耗下的熔 化温度
热管因子和热点因子在反应堆设计时必须设法降低它们的数值。 要减小它们必须从核和工程两方面着手。
沿堆芯径向 装载不同富集度 的核燃料 在堆芯周围设 置反射层 在堆芯径向 不同位置布置一 定数量的控制棒 和可燃毒物棒
合理地控制有 关部件地加工及 安装误差 兼顾工程热管 因子和工程热点 因子数值地减少 和加工费用地增 加
hh,max,5 hn,max
最后综合各分因子求得的总焓升工程热管因子
E E E E E FE F F F F F H H ,1 H ,2 H ,3 H ,4 H ,5
12
9.3 单通道模型的堆芯稳态热工分析
9.3.1 核反应堆热工参数的选择
核电厂反应堆热工参数的选择
适当的选择反应堆热工参数,可以降低电能成本。主要有以下几条: 从提高动力循环热效率 t 来降低电能成本
E H ,1
H n,max
)
由燃料元件冷却剂通道尺寸误差引起焓升工程热管分因子
FE H ,2
W Wh,min,2
[ ADe
(
1b ) 2 b 1 b ) 2 b
]m ]h
[ ADe
(
P
Pn2
2 h ,min
4
2 dcs ,n
4
d cs ,h,max (4 P
DNBR 利用专门公式计算得到的堆内某处的临界热流密度 该处的实际热流密度
在整个堆芯内DNBR的最小值称为最小临界热流密度或最小偏离核 态沸腾比或最小DNBR比。
4
9.2 热管因子和热点因子
9.2.1 核热管因子和核热点因子
积分功率输出最大的冷却剂通道通常就称为热管或热通道; 某一燃料元件表面热流密度最大的点就称为热点。
堆芯热管最大焓升 hh,max F = = 堆芯名义最大焓升 hn,max
E H
堆芯热点最大热流密度 qh,max F = = 堆芯名义最大热流密度 qn,max
E q
7
综合考虑核和工程两方面的因素后,定义热流密度热点因子Fq 和焓升热管因子 FH 为
q n,m ax q h,m ax qh,max Fq=F F = = q q n,m ax q