压水反应堆的热功率.docx
《压水反应堆的热功率.docx》由会员分享,可在线阅读,更多相关《压水反应堆的热功率.docx(14页珍藏版)》请在冰豆网上搜索。
压水反应堆的热功率
华址电力*孑
课程设计报告
(2013--2014年度第二学期)
名
称:
核反应堆热工分析课程设计
题
目:
利用单通道模型进行反应堆稳态热工设计
院
系:
核学院
班
级:
核电1101班
学
号:
13
学生姓名:
漆圣培
指导教师:
李向斌
设计
倜数:
一周
成绩:
、课程设计的目的与要求
对于反应堆
反应堆热工设计的任务就是要设计一个既安全可靠又经济的堆芯输热系统。
热工设计,尤其是对动力堆,最基本的要求是安全。
要求在整个寿期内能够长期稳定运行,并能适应启动、功率调节和停堆等功率变化,要保证在一般事故工况下堆芯不会遭到破坏,甚至在最严重的工况下,也要保证堆芯的放射性物质不扩散到周围环境中去。
在进行反应堆热工设计之前,首先要了解并确定的前提为:
根据所设计堆的用途和特殊要求(如尺寸、重量等的限制)选定堆型,确定所用
反应堆的热功率、堆芯功率分布不均匀系数和水铀比允许的变化范围;
燃料元件的形状、它在堆芯内的分布方式以及栅距允许变化的范围;
冷却剂流过堆芯的流程以及堆芯进口处冷却剂流量的分配情况。
在设计反应堆冷却系统时,为了保证反应堆运行安全可靠,针对不同的堆型,预先规定了热工设计必须遵守的要求,这些要求通常就称为堆的热工设计准则。
目前压水动力堆设计中所规定的稳态热工设计准则,一般有以下几点:
1)燃料元件芯块内最高应低于其他相应燃耗下的熔化温度;2)燃料元件外表面不允许发生沸腾临界;3)必须保证正常运行工况下燃料元件和堆内构件得到充分冷却;在事故工况下能提
供足够的冷却剂以排除堆芯余热;
4)在稳态额定工况和可预计的瞬态运行工况中,不发生流动不稳定性。
在热工设计中,通常是通过平均通道(平均管)可以估算堆芯的总功率,而热通道(热
管)则是堆芯中轴向功率最高的通道,通过它确定堆芯功率的上限,热点是堆芯中温度最高
热工课程设计主要是为了培养学生综合运用反应堆热工分析课程和其它先修课程的理论和实际知识,树立正确的设计思想,培养分析和解决实际问题的能力。
通过本课程设计,达到以下目的:
1、深入理解压水堆热工设计准则;
2、深入理解单通道模型的基本概念、基本原理。
包括了平均通道(平均管)、热通道(热管)、热点等在反应堆设计中的应用;
DNBR最小烧
3、掌握堆芯焓场的计算并求出体现在反应堆安全性的主要参数:
烧毁比毁比MDNB,燃料元件中心温度及其最高温度,包壳表面温度及其最高温度等;
4、求出体现反应堆先进性的主要参数:
堆芯流量功率比,堆芯功率密度,燃料元件平均热流密度(热通量),最大热流密度,冷却剂平均流速,冷却剂出口温度等;
5、掌握压降的计算;
6、掌握单相及沸腾时的传热计算。
7、理解单通道模型的编程方法。
课程设计的考核方式:
1、报告一份;2、计算程序及说明一份;3、答辩。
二、设计任务(设计题目)
已知压水反应堆的热功率
Nt3500MW;燃料元件包壳外径dcs9.5mm
包壳内径dci8.6mm,芯块直径du8.19mm;燃料组件采用17x17正方形排列,
共157组燃料组件;每个组件内有24个控制棒套管和一个中子通量测量管;燃料棒中
心间栅距P=13mm组件间水隙w1mm。
系统工作压力p=,冷却剂平均温度
6%;冷却
tR305C,堆芯冷却剂平均温升t27.4C;冷却剂旁流系数
剂设计总流量14314Kg/s,Fq2.55,FH1.65;dnbr=又设燃料元件内释热
份额占总释热量的%堆芯高度取L=m;并近似认为燃料元件表面最大热流密度、元件
表面最高温度和元件中心最高温度都发生在元件半高度处;已知元件包壳的热导率
kc0.00547(1.8tcs32)13.8[W/(mC)]。
试用单通道模型求燃料元件中心
温度。
3.设计正文
冷却剂有效流量
训卄二WJIO=14314Xa-旳I3455kg/S焓升工程热管因子F為取(未计入交混因子)
1)确定燃料元件的实际最大热流密度qmax
度,故实际的qmax值由热点处的qDNB值除以DNBR而得。
平均热流密度
qmaxDNBR1.651.32.14MW
2)确定燃料元件表面平均热流密度q;q三0.64册伽?
[qi=qnd„=0646X3.14X0.0095=O.OmWI/n,[qi.niax=qiFq=0-0193X2.55=0.0491Mff/in3)堆芯等效直径D硏
(121X
D屏二J式中T为正方形组件每边边长(m)。
因为组件无盒壁,组件间水隙*■二1"",相邻组件的燃料元件棒中心距为13+1=IE,故得
F=(17X13X+2X0,5X10・-6049in^
D吋•2.75<11
堆芯平均管焓升
_Nt3500X1少
少二==13455=力恥他堆芯热管最大焓升
k%押=总h•F^h'=260X10^X1.65X1085=465.SkJAg
借助水和蒸汽计算程序,查得
当切油二286.7C,p=时,Mm=1266.96kj/l倨
故
S+=1266.96+=曲.71kJAgj
借助水和蒸汽计算程序,在工作压力下
t目=327.3?
P5)热管半高处冷却剂流速
热管内冷却剂流速(或质量流速)的精确计算可按教科书中介绍的方法求解,也可按热
管与平均管压降相等的原则进行迭代求解。
作为例子,为简化计算,取热管半高处冷却剂流速近似等于平均管半高处的流速,则
V=
式中:
At为堆芯燃料元件周围的冷却剂总有效流通截面积(为冷却剂平均温度下的
借助水和蒸汽计算程序,可得tfi=305p二1551服i时,比容
Vf=1.3967X107
另一部分是组
At应由两部分组成:
一部分是组件内燃料元件棒之间冷却剂的流通面积;
件间水隙的横截面积,因为流过这个水隙的冷却剂是冷却燃料组件最外面一排燃料元件的,所以它也属于有效冷却剂的流通面积。
因此有
NH:
吧Fr»[4X(17P)6]
式中b为组件间的水隙宽度为。
At=157X264X(13X10一沪—£X(0.0095)^十157L-*J
X(4X[17X13X10■X1X110■=4.20m"
6)计算热管半高度处燃料元件表面与冷却剂间的对流换热系数
=32739P
7)计算燃料兀件表面最咼温度tciJHaK
当乂申(dW比叭』可
当®』扣
•"彳—一32739=1977C
fL\/I65X101“
e.j©=34464.25X(—-—j
由于®・J(m<也&£),故
8)燃料元件包壳内表面最高温度
ZL
ttl.Hdh-十ASc|2
QlFq&C
2~
为燃料元件表面最高温度和内表面最高温度都发生在元件半高度处
Ke=0.00547(1.etesHx卜32)'13.S=0.00547X(1.8X347.16+32)+13.9=17.3911N/(m-C)]
因而
9)燃料芯块表面最高温度
tu.nav.*tcLiiax*△日J?
eg)二%
=328.St:
19.3X10^X2.55
\3r=
理Z6.6+8,19
5678X314X\X10"
于是
I饥肿M=3?
1-?
7328.8=720.77^C
10)计算燃料芯块的中心S高温度环
「血和仲mFqr"站19.3X心X2.55.
4X314
K=IK卜=—二IK+_X10■
查课本表3-7得A冷能K血=A514
;
4忤
K,dt=38.65+39.18=77.83
0
查课本表3-7可得[to.max=21拧C<22O0t
四、课程设计总结或结论
堆芯内部的燃料元件的温度设计终于完成了,经历了很多,也学到了很多。
首先,经过此次的课程设计,掌握了核反应堆堆工设计的基本设计步骤,并在设计中对
其次学会
上课时所学的积分热导率,反应堆和工程热管热点因子等知识从了解到熟练运用;
了office软件中公式编辑器的使用,可以对较难的公式进行编辑,快捷美观;另外,温习
了以前学习的MALAB的编程知识,在编程的过程和编程的能力有一定的提升。
尤其在matlab编程过程中,反复修改程序内容,再次体会细心的重要性,感叹编程的不容易。
理论与实践的结合是实践环节的核心,热工作为一名核电学子不但要学,更是要学好的
专业课,所以在本次的课程设计,不仅学会了课程的知识,而且也学会了使用自己学过的学习工具应用与本门课程,希望在学习过程中各门课也多把学科之间的学习工具应用于其中,更快更好地学习好相关的课程!
!
附录(设计流程图、程序、表格、数据等)
Nt=input('
请输入堆芯功率MW');
d1=input('
请输入燃料元件包壳外径
mm
');
d2=input('
请输入燃料元件包壳内径
mm
');
du=input('
请输入燃料元件芯块直径
mm
');
n仁input('
请输入燃料组件边长排列数:
');
n2=input('
请输入燃料组件数目:
');
x1=input('
请输入水隙宽度mm');
x2=input('
请输入燃料元件栅距mm
');
p0=input('
请输入压强MPa:
');
t0=input('
请输入冷却剂平均温度:
');
dt=input('
请输入堆芯冷却剂温升:
');
w0=input('
请输入冷却剂旁流系数(
%):
');
w1=input('
请输入堆芯体积流量kg/s:
');
fq=input('
请输入热流密度热点因子:
');
fh=input('
请输入焓升热管因子:
');
fu=input('
燃料元件内的释热量占堆总热功率的份额(%):
');
DNBR=input('请输入最小热流密度(DNBR):
');
L=input('请输入燃料原件高度m:
');
dcs=d1/1000;
dci=d2/1000;
xw=x1/1000;
xp=x2/1000;
%确定燃料元件的实际最大热流密度
N=157*(17*17-25)
q=Nt*fu/*dcs*L*N)
qmax=q*fq
%确定燃料元件表面平均热流密度
q1=q*pi*dcs;
q1max=q1*fq;
%堆芯等效直径
T=n1*xp+xw;
Def=sqrt(121*(TA2)/(pi/4));
%热管版高度处水的比焓
hf=+*(Nt*1e6*fh*/(w1*(1-w0/100));
tl=;
at=n2*264*(x卩人2-pi/4*(dcs)A2)+n2*4*n1*xp*xw;
v=(w1*(1-w0/100))*at;
%计算热管半高度处燃料元件表面与冷却剂间的对流换热系数
De=4*(xpA2-pi/4*dcs)/(pi*dcs);
Re=v*De/;
Nu=*ReA*A;
hl=Nu*De;
%计算燃料元件表面最高温度
ol=q*fq/hl;
ojl=345+25*(qmax/1e6)A*exp(-pO/-tl;
ifol<=ojl
tcs=tl+ol;
else
tcs=tl+ojl;
end
%燃料元件包壳内表面最高温度
xc=(dcs-dci)/2;
kc=**tcs+32)+;
oc=q1*fq*xc/(kc*pi*(dcs+dci)/2);
tci=tcs+oc;
%燃料芯块表面最高温度
hg=5678;
og=q1*fq/(hg*pi*(dci+dcs)/2);
tu=tci+og;
%计算燃料芯块的中心最高温度
ou=q1*fq/(4*pi)*;
t=sym('t');
ku=(t++*(t+A3;
fto=double(int(ku,t,O,tu)+ou);
q=int(ku);
p=q-subs(q,O);
to=solve(p-fto);
ifdouble(to)>=22OO
disp('
错误');
else
disp(to);
end