反应堆热工分析课程设计

合集下载

反应堆课程设计

反应堆课程设计

反应堆课程设计一、课程目标知识目标:1. 了解核能反应堆的基本原理,掌握其组成部分及功能;2. 掌握核反应堆的热能转换过程,理解热效率的计算方法;3. 了解核反应堆的安全特性,掌握核安全的基本知识。

技能目标:1. 能够运用所学知识,分析核反应堆的运行原理,并进行简单的热能转换计算;2. 能够通过实例,解释核反应堆的安全措施,评估核事故的风险;3. 能够运用团队合作,设计并展示一个核反应堆模型。

情感态度价值观目标:1. 培养学生对核能的客观认识,提高对能源问题的关注和责任感;2. 增强学生的环保意识,认识到核能在可持续发展中的重要性;3. 培养学生严谨的科学态度,提高对核安全的重视。

本课程针对高年级学生,结合学科特点,注重理论联系实际,提高学生的分析、计算和创新能力。

在教学过程中,关注学生的个体差异,鼓励积极参与,培养团队合作精神。

通过本课程的学习,使学生能够掌握核反应堆的基本知识,具备初步的核能应用能力,并形成正确的能源观和安全观。

二、教学内容1. 核反应堆基本原理:包括核裂变与聚变的概念、链式反应的原理、中子慢化与扩散过程。

- 教材章节:第三章 核反应堆原理2. 核反应堆的组成部分及功能:重点介绍燃料元件、慢化剂、冷却剂、控制棒等。

- 教材章节:第四章 核反应堆的组成部分3. 热能转换过程:讲解核反应堆热能的产生、传递与转换,以及热效率的计算方法。

- 教材章节:第五章 热能转换4. 核反应堆安全特性:阐述核反应堆的安全措施、事故类型及预防措施。

- 教材章节:第六章 核反应堆安全5. 核反应堆模型设计与展示:结合所学知识,进行团队合作,设计并展示核反应堆模型。

- 教材章节:第七章 核反应堆设计与实践本教学内容根据课程目标,科学、系统地安排了核反应堆的基础知识、关键技术和安全特性。

在教学过程中,教师应按照教学大纲,逐步引导学生掌握核反应堆的相关知识,注重理论与实践相结合,提高学生的实际操作能力。

同时,鼓励学生进行团队合作,培养学生的创新意识和实践能力。

核反应堆工程课程设计

核反应堆工程课程设计

成绩评定标准
等级 评定标准 1. 全面完成设计任务,设计内容正确,设计图纸质量高 优秀 2. 设计说明书内容正确,文字精练、流畅、工整 3. 设计过程或答辩过程中,能准确回答与设计内容有关的问题 4. 工作态度认真、严谨、独立工作能力强,模范遵守纪律 1. 全面完成设计任务,设计内容正确 良好 2. 设计说明书内容正确,表达清楚,书写认真 3. 设计过程或答辩过程中,能较好的回答与设计内容有关的问题 4. 工作态度认真,独立完成设计任务,遵守纪律 1. 全面完成设计任务,设计内容正确,设计图纸基本正确 中等 2. 设计说明书内容正确,表达清楚,书写认真 3. 设计过程或答辩过程中,经提示基本回答相关的问题 4. 工作态度认真,具有一定的独立工作能力,遵守纪律 1. 能完成主要设计任务,质量较差或有较大错误,经启发能予纠正 及格 2. 设计说明书内容有个别错误,书写较草 3. 设计过程或答辩时,对有些问题的回答出现概念性的错误 4. 工作态度一般或独立工作能力较差,基本能遵守纪律 1. 未完成设计任务,或设计质量差并不加以改正 不及 格 2. 设计说明书内容有较大错误,或书写非常不认真 3. 设计过程或答辩时,回答问题出现严重的概念错误,或答不出 4. 工作态度不认真,或独立工作能力差,或不遵守纪律
完成日期 学 生


1. 学院、专业、年级均填全称,如:动力工程学院、核工程与核技术、2007。 2. 本表除签名外均可采用计算机打印。本表不够,可另附页,但应在页脚添加页 码。 3. 本课程设计是动力工程学院核能系新开的课程设计,相应任务书需要在该课程 教学和实践中不断完善,因此指导书为修订版。
图目录
图 1 图 2 图 3 图 4 图 5 图 6 图 7 图 8 图 9 图 10 图 12 MCNP5 中对燃料栅元的建模.....................................................................3 keff 与水铀比的关系曲线图..........................................................................3 轴向归一化功率分布,kcode=5000 时...................................................... 4 kcode=50000 时的轴向归一化中子通量密度分布.................................... 6 全燃料组件范围内的中子能谱................................................................... 7 燃料组件内各栅元的相对功率分布........................................................... 8 燃料组件内各栅元的相对功率分布(3D 图像)..................................... 8 可燃毒物布置方案....................................................................................... 9 热管内温度分布......................................................................................... 25 热管内温度分布(不含芯块)............................................................... 25 堆芯内压降沿轴向变化示意图............................................................... 34

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

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

六.计算结果分析:计算结果误差分析:由于采用的是W-3公式,且该设计中的给出参数与该公式的适用范围有些偏差,但是其算出的结果还是能客观反映出热管中各量的变化趋势的。

热管的焓、包壳表面温度、芯块中心温度随轴向的分布如下:控制体为6个:表1 各温度的汇总表各种温度控制体流体出口温度单位(℃)流体出口比焓(kJ/kg)出口处的包壳外壁温度单位℃出口处的包壳内壁温度单位℃出口处的uo2芯块外表面温度单位℃燃料芯块的中心最高温度单位℃堆芯高度L/m第一控制体291.54 1292.1 303.25 303.95 372.25 550 0.61 第二控制体301.29 1343.9 325.71 327.21 472.35 953 1.22 第三控制体315.38 1424.5 348.32 350.42 563.86 1411 1.83 第四控制体330.13 1517.2 348.34 350.44 572.41 1469 2.44 第五控制体339.21 1582.1 348.11 349.41 486.01 939 3.05 第六控制体343.75 1618.8 347.83 348.43 416.73 605 3.66表2 临界热流与烧毁比的汇总表DNBR 控制体DNBR临界热流密度qDNB10^6 单位W/m2第一控制体15.6 5.3 第二控制体 6.5 4.7 第三控制体 3.7 3.9 第四控制体 2.7 3 第五控制体 3.5 2.4 第六控制体 6 22602803003203403600.611.221.832.443.053.66堆芯高度L(m)流体出口温度(℃)图1 流体出口温度(单位℃)分析:由图可知,流体出口温度随着堆芯高度由下到上逐渐上升,到最后一个控制体的末尾,也就是堆芯出口处,达到最大值。

200400600800100012001400160018000.611.221.832.443.053.66堆芯高度L/m流体出口比焓(k J /k g )图2 流体出口比焓(kJ/kg )分析:由图可知,流体出口比焓和流体出口温度一样随着堆芯高度由下到上逐渐上升,到最后一个控制体的末尾,也就是堆芯出口处,达到最大值。

反应堆热工课程设计研究论文[5篇]

反应堆热工课程设计研究论文[5篇]

反应堆热工课程设计研究论文[5篇]第一篇:反应堆热工课程设计研究论文摘要:任务驱动法是一种以任务为主线,以学生为主体的教学方法。

本文结合反应堆热工课程设计的教学实践,探讨任务驱动法在课程设计中的应用。

并以“压水堆堆芯单通道稳态热工分析”为例,通过科学设计课程任务、合理设置实施方案和有效优化评价体系三个过程,阐述任务驱动法在反应堆热工课程设计中的应用。

结果表明不仅有利于提高学生的实践能力,还能培养学生分析问题和解决问题的能力,充分发掘学生的创造潜能。

关键词:任务驱动法;反应堆热工;课程设计引言近年来,核电技术的迅速拓展及深入对高校核电技术人才培养提出了更高的要求。

培养学生实践创新能力已成高校人才培养的主要目标。

课程设计是提高学生实践创新能力的重要途径。

我校针对核工程与核技术专业反应堆工程方向开设反应堆热工课程设计,该课程以我校核能与核技术工程国家级虚拟仿真实验教学中心为依托,是我校培养核专业学生实践能力和创新精神的重要手段。

传统反应堆课程设计教学是以教师为主的手把手式的教学模式,不能充分激发学生的创造潜能。

采用任务驱动教学法,明确“任务”的目标性和教学情境的创建,使学生带着真实任务在探索中学习,从而培养出独立探索、勇于开拓进取的自学能力。

一、反应堆热工课程设计的任务“反应堆热工课程设计”是我校在《反应堆热工》理论课程学习后开展的一门独立课程。

该课程的任务主要包含以下几点:其一,以反应堆相关理论知识为基础,建立符合课程设计任务要求的热工模型;其二,强调利用数学思维解决问题,引导学生学习先进数值模拟方法,其三,以我校国家级核能与核技术虚拟仿真实验教学中心为依托,利用各种编程语言开发求解程序,完成题目内容、分析对比、撰写报告。

以“压水堆堆芯单通道稳态热工分析”为例,其主要任务包含:利用压水堆堆芯单通道相关理论知识,建立堆芯单通道稳态模型,采用数值计算方法如四阶龙格库塔方法求解模型方程组,最后开展计算分析并撰写报告。

反应堆热工分析课程设计全自动版

反应堆热工分析课程设计全自动版

压降的相 关计算 进口温度 tfin 出口温度tfout
287 #NAME?
进口密度ρin #NAME?
出口密度ρout #NAME?
tf #NAME? #NAME?
#NAME? #NAME?
#NAME? #NAME?
#NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
#NAME?
#NAME?
Pr #NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
出口比容υout #NAME?
#NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
ΔPin #NAME?
ΔPout #NAME?
ΔPgr #NAME?
动力粘度 μf1 μf2 μf3 μf4 μf5 μf6
15104.90559
hfs #NAME?
hgs #NAME?
G 2734.147628
32148.52582
K
Cp
Re
h
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
#NAME?
0.9
冷却剂总流量W 8916.666667 两个组件间的水隙
δ 0.0008 FΔHN=FRZ
1.35
堆芯高度L 3.66
ρUO2 95% Kin 0.75
入口温度 假设出口温度
tfin tfout

核反应堆热工分析-课程设计

核反应堆热工分析-课程设计

%另一种带入方式%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)。

核反应堆热工设计 6页word文档

核反应堆热工设计 6页word文档

课程设计报告( 2019 -- 2019 年度第二学期)名称:核反应堆热工分析题目:单通道模型反应堆稳态热工设计院系:核科学与工程学院班级:学号:学生姓名:指导教师:设计周数: 2成绩:日期:2019年 6 月26 日一、课程设计的目的与要求该课程设计的主要目的为:培养学生综合运用反应堆热工分析课程和其它先修课程的理论和实际知识,树立正确的设计思想,培养分析和解决实际问题的能力。

该课程设计的基本要求为:在堆型和为进行热工设计所必需的条件已经确定的前提下,利用单通道模型进行反应堆稳态热工设计,并对热工设计准则进行验证;二、设计正文已知压水反应堆的热功率N t=2895Mw;燃料元件包壳外径d cs=9.5mm,包壳内径d ci=8.36mm,芯块直径d u=8.19mm;燃料组件采用17*17正方形排列,共157组燃料组件;每个组件内有24个控制棒套管和一个中子通量测量管;燃料棒中心间栅距P=12.6mm,组件间水隙δw=1mm。

系统工作压力p=15.5Mpa,冷却剂平均温度t R=310℃,堆芯冷却剂平均温升∆t=34.6℃;冷却剂旁流系数ζ=6.5%;冷却剂设计总流量71370m3/h;DNBR=2.08;又设燃料元件内释热份额占总释热量的97.4%;堆芯高度取L=3.66m;并近似认为燃料元件表面最大热流密度、元件表面最高温度和元件中心最高温度都发生在元件半高度处;已知元件包壳的热导率k c=0.00547(1.8t cs+32)+13.8[W/(m·℃)]。

适用单通道模型求燃料元件中心温度。

(大亚湾)求解步骤:由题知F u=97.4% , W ef=1-6.5%=93.5%。

取,F∆H N=1.51, F∆H E=1.07,;F Z N=1.54;F q E=1.03 F q=2.39 F q N=F qF q E=2.32相应温压下,水的c p=5.816kJ/(kg.C)1、确定燃料元件的实际最大热流密度q max因为压水堆的安全限值首先是燃料元件表面的最小DNBR,其次才是燃料元件的中心温度,故q max值由热点处的q DNB值除以DNBR而得。

核反应堆初步设计课程设计说明书

核反应堆初步设计课程设计说明书

z
表达式................................................................................... 10
2.3.2 热管的温度场表达式...........................................................................................11 2.3.3 包壳表面温度分布函数...................................................................................... 12 2.3.4 燃料中心温度的分布函数.................................................................................. 13 第 3 章 反应堆稳态工况水力计算...............................................................................................15 3.1 引言.....................................................................................................................................15 3.2 一回路内的流动压降和水泵功率................................................................................. 15 3.3 反应堆自然循环能力分析.................................................................................................17 第 4 章 反应堆压力容器设计与部件图绘制................................................................................ 17 4.2 薄壁容器的应力分析和强度设计................................................................................... 18 4.3 管道的强度分析.................................................................................................................19 第 5 章 计算程序与计算结果.........................................................................................................20 5.1 Matlab 程序如下................................................................................................................ 20 5.2 程序计算结果.....................................................................................................................28 参考文献........................................................................................................................................... 33 附录 1:反应堆稳态热工计算表....................................................................................................34 附录 2 燃料温度校核计算表..........................................................................................................38 附录 3 反应堆稳态工况水力计算表..............................................................................................41 附录 4 反应堆压力容器强度计算表..............................................................................................41
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
4
当用单通道模型分析计算时,热管冷却剂焓升为:
ΔH h ,max
∫ =
L
0
E qFΔN H FΔHϕ ( z ) dz
W

∫ =
L
0
qFRN FΔE Hϕ ( z ) dz W
此时,燃料元件表面最大热流量为:
qh ,max = q FqN FqE
二、设计目的
通过本课程设计,达到以下目的: 1.深入理解压水堆热工设计准则; 2.深入理解单通道模型的基本概念、基本原理。包括了解平均通道(平均管) 、热通道 (热管) 、热点等在反应堆热工设计中的应用; 3.掌握堆芯焓场的计算并求出体现反应堆安全性的主要参数;烧毁比 DNBR,最小烧 毁比 MDNBR,燃烧元件中心温度 t0 及其最高温度 t0,max ,包壳表面温度 tcs 及其最高温度
FqN = FRN ⋅ FZN =
堆芯最大热流量 堆芯平均热流量
N
在局部峰(控制棒、空泡、水隙所致)明显的堆内还必须考虑局部峰核热管因子 FL 的 影响。因而,由总的通量分布(径向、轴向和局部峰)引起的热流量核热点因子 Fq 为
N
FqN = FRN ⋅ FZN ⋅ FLN
在实际计算中,还应考虑到在堆芯核设计中如应用 R-Z 坐标计算时方位角影响,以及核计 算不准确性所造成的误差,故 Fq 为
N
3
表示。在反应堆发展早期,人为地把热点位于热管内,故 Fq 也称为热流量核热管因子。 核热管因子的定义是:
N
FRN =径向核热管因子=
热管的平均热流量 堆芯平均管的平均热流量 热管的最大热流量 热管的平均热流量
N
FZN =轴向核热管因子=
故,因中子通量分布(径向和轴向)引起的热流量核热点因子 Fq 可表示为:
考虑到堆芯内组件布置的对称性,对称安排后应取 n=121。 5.求堆芯等效直径 Def
2 Def =
4nT 2
π
式中:T 为正方形组件每边长,m。 设燃料组件无盒壁,考虑到装卸料的要求,组件间的水隙取为 0.8mm,即相 邻组件的燃料棒中心距为
T 2 = (15 × 13.3 × 10−3 + 0.8 × 10−3 ) 2 = 0.04m 2
o o
6
外径 d cs 为 10mm,包壳内径 d ci 为 8.7mm,芯块直径 du 为 8.53mm;栅距 P 为 13.3mm。燃 料元件采用正方形栅格排列,每个燃料组件中燃料元件排列的数目为 15 × 15 ,而其中有 20 根是控制棒套管,一根是中子通道测管,因而每个燃料组件内实际有 204 根燃料元件棒。燃 料元件内的释热量占堆芯总功率的 97.4%,可忽略冷却剂的释热量。旁流系数为 9%。设热 点位于热管内, FR = FΔH = 1.8 ,热管轴向归一化功率分布的最大值 FZ = 1.78 ,热流量
3.求全堆燃料棒根数 N 为:
N=
取堆芯高度为 2.9m,有
N t Fu π d cs Lq
N=
953 × 106 × 0.974 = 2470 3.14 × 0.01× 2.9 × 0.42 × 106
4.确定堆芯燃料元件组件数 n 因为每个组件内燃料棒根数为 204,则
n=
N 24270 = 119 204 204
N N N
焓升工程热管因子 FΔH = 1.08 。 要求燃料元件中心最高温度不超 工程热点因子 Fq = 1.04 ,
E
E
过 2200 C 。试用单通道模型对堆芯进行热工分析。 求解步骤: 1.确定燃料元件的实际最大热流密度 qmax
o
5
压水动力堆的安全限制首先是燃料元件表面的 MDNBR,其次才是燃料元件的中心温 度。因此, qmax 值由热点处的 qDNB 值除以 DNBR 而得。 取 qDNB = 2.8MW / m , DNBR = 2.0
N
FqN = FRN ⋅ FZN ⋅ FLN ⋅ FUN ⋅ FθN
式中 FU 为核计算误差修正系数, Fθ 为方位角修正系数。即热流量核热点因子为:
N N
FqN =
N N
q max N N N N N = FR ⋅ FZ ⋅ FL ⋅ FU ⋅ Fθ q
一般将 FL 及 Fθ 归并在 FR 中。 焓升核热管因子为 FΔH =
《核反应堆热工分析》课程设计



编 写:苏光辉 郭玉君 审 校:贾斗南
ቤተ መጻሕፍቲ ባይዱ
西安交通大学 核科学与技术学院 二零零六年十一月
1
一、概述
在利用核能作为能源的核动力装置中, 反应堆及堆芯输热系统时将核能转变为其他形式 能量的一个中间枢纽。因此,反应堆热工设计在整个核反应堆设计中占有极其重要的地位。 反应堆热工设计的任务实设计出一个安全、可靠、经济的堆芯输热系统。做好反应堆热工设 计并非易事,他涉及的面广,与许多方面的设计有关,如反应堆本体结构、堆物理、堆控制、 堆材料及二回路系统等等。 在进行堆芯输热系统的设计时,必须保证反应堆安全、可靠。为此,针对不同的堆型, 预先制定了热工设计必须遵守的要求,在反应堆的整个运行寿期内,不论是处于稳态工况, 还是处于预期的事故工况, 反应堆的热工参数都必须满足这些要求。 这些要求就是通常所讲 的反应堆的热工设计准则。这些准则非常重要,它不断是热工设计的依据,也是安全保护系 统设计的原始条件,而且也是制定运行规程的出发点。反应堆热工设计准则的内容,不仅随 反应堆的堆型有关, 而且随科学技术的发展, 堆设计与运行经验的积累以及堆用材料性能和 加工工艺等的改进而发生变化。 目前,压水动力堆的稳态热工设计准则有: 1. 燃料元件外表面不允许发生沸腾临界。 通常用 DNBR 来定量地表示这个限制条件。DNBR 是根据堆内某处燃料元件周围的冷 却剂状态使用专门的计算公式而得到的临界热流密度与该处燃料元件表面的实际热流密度 的比值。DNBR 随堆芯通道的长度是变化的,在整个堆芯内,DNBR 的最小值称为最小 DNBR, 用 MDNBR 或 DNBRmin 表示。 为了确保燃料元件不烧毁, 当计算的最大热功率下, MDNBR 不应低于某一规定值。如果计算热流密度的公式没有误差,则当 MDNBR=1 时, 表示燃料元件表面要发生沸腾临界。若该公式存在误差,则 MDNBR 就要大于 1。例如, W-3 公式的误差为 23%,所以当使用 W-3 公式计算 DNBR 时,就要求 MDNBR ≥ 1.3。 2. 燃料元件芯块最高温度应低于熔化温度。 目前,压水堆大多采用 UO2 作为燃料。二氧化铀的熔点约为 2805 ±15 C ,经辐照后, 其熔点会有所降低。燃耗每增加 10 兆瓦 i 日/吨铀,其熔点下降 32 C 。在通常所达到的燃 耗深度下,熔点将降至 2650 C 左右。在稳态热工设计中,一般将燃料元件中心最高温度限 制在 2200~2450 C 之间。 3. 在稳态额定工况下,要求在计算的最大热功率下,不发生流动不稳定性。 对于压水堆, 只要在堆芯最热通道出口附近冷却剂中的含气量不大于某一数值, 就不会 发生流动不稳定性。 由于气冷堆不存在压水堆那样的燃料元件表面的沸腾临界问题,所以气冷堆的热工设 计准则与压水堆有所不同。 它主要要求燃料元件表面最高温度、 中心最高温度以及燃料元件 和结构部件的最大热应力不超过允许值。 对于用水作冷却剂的生产堆, 一般就把燃料元件包壳与水发生加速腐蚀时的包壳表面温
2
qmax =
qDNB 2.8 = = 1.40 MW / m 2 DNBR 2.0

2.确定燃料元件表面平均热流密度 q
q=


qmax 1.40 = = 0.42 MW / m 2 N E Fq Fq 1.80 × 1.78 × 1.04
− −
q l = q π d cs = 0.42 × 3.14 × 0.01 = 0.013MW / m 2 ql ,max = q l FqN FqE = 0.044 MW / m
2
4
度作为其设计准则之一。 这是因为新的裂变燃料的生产量与堆的热功率成正比, 要增加新的 裂变燃料的生产量, 就必须尽可能提高反应堆的热功率; 而功率的提高会使燃料的表面温度 也跟着升高。当这个温度高到某个数值时,包壳就会被水加速腐蚀,从而影响燃料元件的寿 命。 但包壳加速腐蚀并不会立即影响到反应堆的安全, 但它会影响反应堆燃料的换料周期和 运行周期,从而影响新的裂变燃料的生产量。因此,把引进包壳加速腐蚀的温度确定为生产 堆的一条热工设计准则。 重水反应堆堆芯的性能要求和安全要求之一就是正常运行、 运行瞬态和中等频率事故过 程中防止燃料发生由热工或水力导致的失效。 为了满足这个准则, 规定了下列堆芯热工水力 设计基准: (1)临界功率比(CPR) ,在正常运行和缓慢失去反应性控制条件下,保证具有堆芯最 小临界功率比的燃料元件包壳表面不发生烧干(或称干涸) 。使用临界含汽量——沸腾长度 ( X C − BL )关系式以及考虑了关系式的不确定性后的临界功率比控制在 1.30 以上就能满 足这个设计标准。这个 30%的设计临界功率比裕量考虑了如下影响:换料时功率波动、临 界功率比分析中的不确定性、仪表误差以及通道功率测量不确定性等。 (2) 水力学稳定性 在正常运行工况和超功率瞬态过程中不应导致燃料管道组件内流动 不稳定性。 (3) 燃料管道组件流量 在所有四台主泵都运行的情况下, 每个燃料管道组件中的实际 流量必须大于和等于设计流量,但小于上限流量。 每根燃料管道组件中的设计流量以燃料管道组件中的功率为基础的,以便在电厂寿期 末的正常满功率运行条件下在出口集管处获得大约 4%的含汽量。 每根燃料管道的组件上限流量是疲劳试验的结果为基础的。上限流量被定义为这样一 个流量,即在这个流量下,流动所引起的燃料棒束的振动水平,棒束的磨损以及压力管的磨 损都是可以接受的,已有的试验结果表明:在单相流时,流量可上升至 30kg/s,在两相流时 流量上升至 24kg/s,不会引起过分的流致振动和棒束以及压力管的磨损。 ( 4 )燃料温度 在正常运行和中等频率事故下燃料的峰值温度小于它的熔化温度(~ 2840℃) 。同时也要求保证燃料元件包壳和燃料棒束的完整性。 以上这些设计基准要求与压水动力堆的设计准则相当的。 在反应堆内,即使燃料元件的形状、尺寸、密度和裂变物质的浓缩度都相同,堆芯内的中子 通量分布也是不均匀的,再加上堆芯内存在控制棒,水隙、空泡及反射层的影响,中子通量 的分布更是不均匀的。从而,堆芯内的热功率分布也是不均匀的。而燃料元件在加工、安装 及运行中的各类工程因素也能造成实际值与设计值之间产生偏差。 为了表示有关的热工参数 的最大值偏离平均值(或名义值)的程度,引入了热管因子的概念。分核热管因子和工程热 管因子。 堆芯功率分布(有时称堆芯功率整体分布)的不均匀程度常用热流量核热点因子 Fq 来
相关文档
最新文档