2cs

流体的kz、μ（z）和Pr数根据流体的压力好温度由表查得。（kλ传热系
数）
如果流体已经达到过冷沸腾，用Je
sLottes公式：
tw
ts

25
qFRNzFqE106
025

e
p62

fhz2twztfhztwztststfhz
5
f
25

qFRN
Z106
FqE
025
e
p62


tstfhz
当
fhz2
fhz2时，用前面的式子
当
fhz2
fhz2时，用fhz2代替fhz1
（3）第一个控制体出口处的包壳内壁温度
tcihztcshzchz

tcshz

qFRNFqE2kCz
z
l


dcsdci

式中：Zr4的kc00054718tc32138
tc

z


12
tcsztciz
4第一个控制体出口处的UO2芯块外表面温度
tuhz

tcihz

qFRNFqEdcidu
2
Zhg
（5）第一个控制体出口处的UO2芯块中心温度用积分热导求解的方法，即
Wm℃
kt0hz
0
u
t
dt
ktuhz
0
u
t
dtqFRNFqE4
Z
其他2个控制体的计算方法相同，重复上述过程即可。
q36热管中的DNBhz
用w3公式计算，同样对3个控制体都算
6
f37DNBR的计算
DNBRqDNBhz
qDNBhz
qzqFRNFqEz
38计算热管中的压降39单相流体的摩擦压降
式中：
LV2
pf

fDe
2
fLG2vDe2
f

fiso

wf




0316Re025
4


wf
06
单相流体加速压降：
pel
单相流体提升压降：

gL

i

out2
gL
局部压降，出口：
pout
kout
V2outout2

kout
G
v2out
2
pi
进口：

ki

i
Vi
2
2

ki

G2vi
2
pgr
定位格架出口压降：
kgr
grV22
G21kgr2
voutvi
2
其中，比容v按相应的流体压力和温度，由表查得。
7
f四、程序设计框图
开始
读输入参数
计算有关堆参数
估算控制体出口温度tf
计算控制体出口温度θ
计算该处含汽量
根据W3算临界热流
计算烧毁比
1θ0001
tf
是计算包壳外表面温
度包壳内表面温度
芯块表面温度
芯块中心温度
打印输出值
停机
8
否重估tf
f五、代码说明书
本代码主要由五个小部分组成。堆芯出口温度计算、堆热流量计算、堆平均参数计算、第一至第六控制体各量计算、热管的压降计算。
（1）堆芯出口温度计算：此段根据任务书给出的基本参数和热量与流量之间关系，运用迭代的算法，求出堆芯的出口温度。（2）堆热流量计算：先根据堆芯的输出功率和释热率以及总的传热面积，求出燃料元件表面平均热流量，再根据热管因子求出最大热流量。再求出平均线功率和最大线功率。（3）堆平均参数计算：根据基本的尺寸，求出堆体的流通截面积和一个栅元的流通截面积。然后再求出流经栅元的流量。依据上面的温度结果，查出热物性参数，再求出冷却剂的流速。（4）第一r