一种锂冷核反应堆与斯特林发电机耦合传热特性分析方法与流程

1.本发明涉及核反应堆领域换热技术，具体涉及一种锂冷核反应堆与斯特林发 电机耦合传热特性分析方法。


背景技术：

2.目前，动态能量转换如斯特林循环、布雷顿循环等，因其转换效率高已成为 大功率空间核反应堆设计采用的主要能量转换方式。斯特林循环结合液态金属锂 冷却反应堆可实现百千瓦级的功率输出，是大型空间堆可行设计方案之一。为了 确定复杂宇宙环境下，锂冷反应堆耦合斯特林发电机瞬态传热特性的影响，建立 堆芯热工水力模型、斯特林发电机系统换热模型，从而为更加全面、有效地评估 复杂宇宙环境下斯特林发电机耦合锂冷反应堆系统安全运行特性提供依据。
3.

技术实现要素：

4.为了克服上述现有技术存在的问题，本发明的目的在于提供一种锂冷核反应 堆与斯特林发电机耦合传热特性分析方法，可以准确的反映出系统的运行特性， 可以实现不同结构与功率的空间锂冷核电源系统，包括堆芯、斯特林发电机的计 算，降低了对空间锂冷核电源的结构与参数的要求，有效增加了本方法对不同问 题的适应性。
5.为了实现上述目的，本发明采取了以下技术方案：
6.一种锂冷核反应堆与斯特林发电机耦合传热特性分析方法，该方法包括以下 步骤：
7.步骤1：输入锂冷核反应堆以及斯特林发电机的结构、几何参数，确定堆芯 燃料的功率分布和冷却剂的温度，输入反应堆堆芯流量及堆芯斯特林发电机冷端 温度，设定计算时间；
8.步骤2：对核反应堆堆芯沿径向及轴向划分若干控制体，初始时刻设定各控 制体温度、压力；
9.步骤3：建立关于堆芯控制体的非线性微分方程，通过吉尔算法得到当前时 刻堆芯控制体温度；
10.锂冷核反应堆堆芯由燃料元件和液态锂冷却剂构成，燃料元件主要由燃料区 域、裂变气体间隙、燃料包壳构成；
11.首先采用考虑六组缓发中子的点堆模型动态方程来求解反应堆裂变功率：
[0012][0013][0014]
式中：
[0015]
p(t)——t时刻反应堆裂变功率/w；
[0016]
t——计算时间/s；
[0017]
λ——中子代时间/s；
[0018]
β——总的有效缓发中子份额；
[0019]
β
i
——第i组缓发中子份额；
[0020]
λ
i
——第i组缓发中子的衰变常数/s
‑1；
[0021]
c
i
(t)——t时刻第i组缓发中子的浓度/m
‑3；
[0022]
n
c
——缓发中子组数；
[0023]
ρ(t)——总的反应性/$；
[0024]
以上计算点堆模型动态方程得到反应堆裂变功率，将反应堆裂变功率折算为 控制体内热源：
[0025][0026]
q
v,i
——燃料控制体i的内热源/w
·
m
‑3；
[0027]
p(t)——t时刻反应堆裂变功率/w；
[0028]
λ
p,i
——输入第i个控制体的功率分布系数
[0029]
v
i
——控制体i的体积/m3；
[0030]
燃料芯块及包壳、气隙区域温度变化率计算如下：
[0031][0032][0033]
式中：
[0034]
ρ
u
——燃料芯块的密度/kg
·
m
‑3；
[0035]
c
u
——燃料芯块的比热/j
·
kg
‑1·
k
‑1；
[0036]
t
u
——燃料芯块的温度/k；
[0037]
λ
u
——燃料芯块的导热系数/w
·
m
‑1·
k
‑1；
[0038]
r——燃料芯块的半径/m；
[0039]
q
v
——燃料控制体的热源密度/w
·
m
‑3；
[0040]
ρ
i
——区域i的密度/kg
·
m
‑3；
[0041]
c
i
——区域i的比热/j
·
kg
‑1·
k
‑1；
[0042]
t
i
——区域i的温度/k；
[0043]
λ
i
——区域i的导热系数/w
·
m
‑1·
k
‑1；
[0044]
r
i
——区域i的半径/m；
[0045]
g1——包壳区域
[0046]
g2——气隙区域
[0047]
假设液态锂冷却剂在堆芯内的流动为不可压缩的一维单向流动，控制方程如 下：
[0048]
质量控制方程：
[0049][0050]
动量控制方程：
[0051][0052]
能量控制方程：
[0053][0054]
式中：
[0055]
ρ——冷却剂密度/kg
·
m
‑3[0056]
t——时间/s
[0057]
w——冷却剂流量/kg
·
s
‑1[0058]
p——冷却剂压力/pa
[0059]
f——摩擦系数
[0060]
a——冷却剂流通面积/m2[0061]
h——冷却剂比焓/j
·
kg
‑1[0062]
z——冷却剂轴向高度/m
[0063]
u——冷却剂加热周长/m
[0064]
d——冷却剂水力直径/m
[0065]
方程组(6)
‑
(7)获得堆芯流量随时间及空间的变化率，求解方程(8)获 得堆芯流体温度随时间的变化；利用吉尔算法求解上述锂冷核反应堆堆芯传热方 程组(4)、(5)、(8)获得芯块、包壳、气隙、堆芯流体的温度分布；
[0066]
步骤4：计算当前时刻斯特林发电机换热量；
[0067]
堆芯流体从堆芯出口流经斯特林发电机热端后,重新流入堆芯进口，斯特林发 电机热端温度为堆芯进出口平均温度；斯特林发电机模型包括斯特林发电机冷热 端壁面、气体工质、回热器部分；对于斯特林发电机冷热端包壳有如下热平衡方 程：
[0068][0069][0070]
式中：
[0071]
t
h
——堆芯出口流体温度
[0072]
t
c
——斯特林冷端温度
[0073]
π
in
——热端换热周长
[0074]
π
in
——冷端换热周长
[0075]
ρ
in
——热端包壳材料密度/kg
·
m
‑3[0076]
ρ
out
——冷端包壳材料密度/kg
·
m
‑3[0077]
t
in
——热端包壳内壁温度/k
[0078]
t
out
——冷端包壳内壁温度/k
[0079]
t
f
——为气体温度/k
[0080]
a
in
——热端截面面积/m2[0081]
a
out
——热端截面面积/m2[0082]
λ
in
——热端材料热导率/w
·
m
‑1·
k
‑1[0083]
λ
out
——冷端材料热导率/w
·
m
‑1·
k
‑1[0084]
δ
in
——热端包壳壁面厚度/m
[0085]
δ
out
——冷端包壳壁面厚度/m
[0086]
c
in
——热端包壳材料比热容/j
·
kg
‑1·
k
‑1[0087]
c
out
——冷端包壳材料比热容/j
·
kg
‑1·
k
‑1[0088]
η
st
——斯特林发电机效率
[0089]
加热器、冷却器与回热器为有限温差传热，具有相似控制方程，其循环换热 量如下：
[0090]
q＝ha
w
(t
w
‑
t
g
)(1/n)
ꢀꢀ
(11)
[0091]
式中：
[0092]
q——循环换热量/w
[0093]
h——工质与壁面平均换热系数/w
·
m
‑2·
k
‑1[0094]
t
w
——内壁面温度/k
[0095]
t
g
——气体温度/k
[0096]
n——机轴转速/r
·
s
‑1[0097]
对于实际斯特林发电机回热器，有一定回热损失，定义有效度ε，ε为实际 循环回热量与理想循环回热量之比；回热器部分的冷热流体的热容量流率相同， 若忽略回热器导热热阻，则有：
[0098][0099]
式中：
[0100]
a
wg
——换热面积/m2[0101]
c
p
——工质热容/j
·
kg
‑1·
k
‑1[0102]
a——工质自由流动面积/m2[0103]
st——斯坦顿数，对于工质氦气，st＝0.46re
‑
0.4
pr
‑1。
[0104]
回热器的轴向热损失可由下式计算：
[0105][0106]
式中：
[0107]
q
loss
——轴向热损失/w
[0108]
λ
r
——回热器外壳处导热系数/w
·
m
‑1·
k
‑1[0109]
a
r
——回热器横截面积/m2[0110]
l
r
——为回热器长度/m
[0111]
t
h
——回热器热侧平均温度/k
[0112]
t
c
——回热器冷侧平均温度/k
[0113]
n——机轴转速/r
·
s
‑1[0114]
通过式(11)
‑
(13)的求解，获得当前时刻斯特林发电机的实际换热量：
[0115]
q
ac
＝q
·
ε
‑
q
loss
ꢀꢀꢀ
(14)
[0116]
式中：
[0117]
ε——有效度
[0118]
q——循环换热量/w
[0119]
q
loss
——轴向热损失/w
[0120]
q
ac
——实际换热量/w
[0121]
步骤5：根据所有已知条件，利用吉尔算法求解下一时刻的堆芯温度、压力 等参数，循环计算直到达到设定时间。与现有技术相比，本发明有如下突出特点：
[0122]
对锂冷核反应堆与斯特林发电机耦合的系统进行研究，可以准确计算系统整 体温度分布，可以实现不同结构与功率的锂冷核反应堆与斯特林发电机的系统安 全分析，包括了堆芯、斯特林发电机的传热计算，降低了对锂冷核反应堆的结构 与参数的要求，有效增加了本方法对不同问题的适应性。本方法能计算通用的锂 冷核反应堆与斯特林发电机耦合的瞬态传热特性，为此类系统瞬态运行时运行策 略、电气系统控制方案等提供了研究方法。
附图说明
[0123]
图1为本发明方法流程图。
具体实施方式
[0124]
下面结合附图和具体实施方式对本发明作进一步详细说明：
[0125]
本发明一种锂冷核反应堆与斯特林发电机耦合传热特性分析方法，采用吉尔 算法求解锂冷核反应堆堆芯温度分布，通过迭代求解每一时刻斯特林发电机换热 量。如图1所示，该方法具体流程包括以下步骤：
[0126]
步骤1：输入锂冷核反应堆以及斯特林发电机的结构、几何参数，确定堆芯 燃料的功率分布和冷却剂的温度，输入反应堆堆芯流量及堆芯斯特林发电机冷端 温度，设定计算时间；
[0127]
锂冷反应堆为液态锂用作冷却剂的反应堆，斯特林发电机可实现将反应堆热 能转换为电能的动态能量转换。通过设置步骤1锂冷核反应堆及斯特林发电机的 几何、结构参数，并通过步骤2
‑
5能实现对不同几何结构的锂冷核反应堆与斯特 林发电机的耦合传热计算。
[0128]
步骤2：对核反应堆堆芯沿径向及轴向划分若干控制体，初始时刻设定各控 制体温度、压力；
[0129]
控制体数目可任意设置，提高控制体数目可增加计算精度，初始时刻堆芯各 控制体温度压力与初始输入冷却剂参数一致。
[0130]
步骤3：建立关于堆芯控制体的非线性微分方程，通过吉尔算法得到当前时 刻堆
芯控制体温度；
[0131]
锂冷核反应堆堆芯由燃料元件和液态锂冷却剂构成，燃料元件主要由燃料区 域、裂变气体间隙、燃料包壳构成；对堆芯的传热计算流程包括确定燃料元件产 生功率、燃料元件导热、液态锂冷却剂流动传热过程。
[0132]
首先采用考虑六组缓发中子的点堆模型动态方程来求解反应堆裂变功率：
[0133][0134][0135]
式中：
[0136]
p(t)——t时刻反应堆裂变功率/w；
[0137]
t——计算时间/s；
[0138]
λ——中子代时间/s；
[0139]
β——总的有效缓发中子份额；
[0140]
β
i
——第i组缓发中子份额；
[0141]
λ
i
——第i组缓发中子的衰变常数/s
‑1；
[0142]
c
i
(t)——t时刻第i组缓发中子的浓度/m
‑3；
[0143]
n
c
——缓发中子组数；
[0144]
ρ(t)——总的反应性/$；
[0145]
以上计算点堆模型动态方程得到堆芯裂变功率，将裂变功率折算为控制体内 热源：
[0146][0147]
q
v,i
——燃料控制体i的内热源/w
·
m
‑3；
[0148]
p(t)——t时刻反应堆裂变功率/w；
[0149]
λ
p,i
——输入第i个控制体的功率分布系数
[0150]
v
i
——控制体i的体积/m3；
[0151]
燃料芯块及包壳、气隙区域温度变化率计算如下：
[0152][0153][0154]
式中：
[0155]
ρ
u
——燃料芯块的密度/kg
·
m
‑3；
[0156]
c
u
——燃料芯块的比热/j
·
kg
‑1·
k
‑1；
[0157]
t
u
——燃料芯块的温度/k；
[0158]
λ
u
——燃料芯块的导热系数/w
·
m
‑1·
k
‑1；
[0159]
r——燃料芯块的半径/m；
[0160]
q
v
——燃料控制体的热源密度/w
·
m
‑3；
[0161]
ρ
i
——区域i的密度/kg
·
m
‑3；
[0162]
c
i
——区域i的比热/j
·
kg
‑1·
k
‑1；
[0163]
t
i
——区域i的温度/k；
[0164]
λ
i
——区域i的导热系数/w
·
m
‑1·
k
‑1；
[0165]
r
i
——区域i的半径/m；
[0166]
g1——包壳区域
[0167]
g2——气隙区域
[0168]
假设液态锂冷却剂在堆芯内的流动为不可压缩的一维单向流动，控制方程如 下：
[0169]
质量控制方程：
[0170][0171]
动量控制方程：
[0172][0173]
能量控制方程：
[0174][0175]
式中：
[0176]
ρ——冷却剂密度/kg
·
m
‑3[0177]
t——时间/s
[0178]
w——冷却剂流量/kg
·
s
‑1[0179]
p——冷却剂压力/pa
[0180]
f——摩擦系数
[0181]
a——冷却剂流通面积/m2[0182]
h——冷却剂比焓/j
·
kg
‑1[0183]
z——冷却剂轴向高度/m
[0184]
u——冷却剂加热周长/m
[0185]
d——冷却剂水力直径/m
[0186]
方程组(6)
‑
(7)获得堆芯流量随时间及空间的变化率，求解方程(8)获 得堆芯流体温度随时间的变化；利用吉尔算法求解上述锂冷核反应堆堆芯传热方 程组(4)、(5)、(8)获得芯块、包壳、气隙、堆芯流体的温度分布；
[0187]
步骤4：计算当前时刻斯特林发电机换热量；
[0188]
堆芯流体从堆芯出口流经斯特林发电机热端后,重新流入堆芯进口，斯特林发 电机热端温度为堆芯进出口平均温度；斯特林发电机模型包括斯特林发电机冷热 端壁面、气体工质、回热器部分；对于斯特林发电机冷热端包壳有如下热平衡方 程：
[0189]
[0223]
a——工质自由流动面积/m2[0224]
st——斯坦顿数，对于工质氦气，st＝0.46re
‑
0.4
pr
‑1。
[0225]
回热器的轴向热损失可由下式计算：
[0226][0227]
式中：
[0228]
q
loss
——轴向热损失/w
[0229]
λ
r
——回热器外壳处导热系数/w
·
m
‑1·
k
‑1[0230]
a
r
——回热器横截面积/m2[0231]
l
r
——为回热器长度/m
[0232]
t
h
——回热器热侧平均温度/k
[0233]
t
c
——回热器冷侧平均温度/k
[0234]
n——机轴转速/r
·
s
‑1[0235]
通过式(11)
‑
(13)的求解，获得当前时刻斯特林发电机的实际换热量：
[0236]
q
ac
＝q
·
ε
‑
q
loss
ꢀꢀ
(28)
[0237]
式中：
[0238]
ε——有效度
[0239]
q——循环换热量/w
[0240]
q
loss
——轴向热损失/w
[0241]
q
ac
——实际换热量/w
[0242]
步骤5：根据所有已知条件，利用吉尔算法求解下一时刻的堆芯温度、压力 等参数，循环计算直到达到设定时间。
[0243]
步骤1
‑
5能准确计算从堆芯至斯特林发电机的传热过程，并获得堆芯与斯特 林发电机实际换热量，能计算得出不同时刻下的堆芯或斯特林发电机各控制体的 实时温度。