一种空间核电推进系统球床反应堆物理热工耦合分析方法技术方案

一种空间核电推进系统球床反应堆物理热工耦合分析方法，主要步骤如下：1、划分通道和控制体，输入球床堆的结构参数和边界条件，设定计算时间和时间步长。2、堆芯初始化计算，得到零时刻的初始值。3、采用吉尔算法，计算当前t时刻的堆芯功率分布、流场、温度场和压力场。4、根据步骤3的计算结果，采用吉尔算法进行下一时刻的计算，直到设定的总计算时间，计算停止。5、若计算结果收敛，终止计算；否则读取上一次停止时刻的计算结果作为初值，重新启动计算，直至结果收敛。本发明专利技术可以分析球床堆堆芯的稳态热工水力特性及事故工况下的瞬态响应特性，为核电推进球床堆的设计及安全特性分析提供建议与指导。

A physical and thermal coupling analysis method for pebble bed reactor of space nuclear propulsion system

【技术实现步骤摘要】
一种空间核电推进系统球床反应堆物理热工耦合分析方法
本专利技术涉及空间核动力
，具体涉及一种空间核电推进系统球床反应堆物理热工耦合分析方法。
技术介绍
随着空间安全和探索任务的拓展，航天器将面临更为复杂、严峻的空间环境。传统的太阳能或蓄电池航天器功率范围小、环境依赖性强，已无法满足未来的任务需求。对于深空探测、载人火星飞行和一些对机动性要求高的军事任务等，环境适应性强、大功率长寿命、安全可靠的空间核动力系统成为最优甚至是唯一的选择。在20世纪60年代，美国和前苏联就成功在航天任务中使用了空间核电源，此后，美国、俄罗斯、法国、德国和日本等航天大国都积极开展了将核能应用于太空的探索和研究。为了满足我国未来军事需求与国家战略需要，研究空间核动力技术刻不容缓。在空间核动力系统中，高温气冷球床反应堆具有功率密度高、结构紧凑、裂变产物完全包容等优点，可广泛应用于核热推进、核电推进及双模式推进，特别地，球床反应堆结合闭式布雷顿循环可满足空间兆瓦级的功率需求。针对用于核热推进的球床反应堆，新墨西哥大学的空间核动力研究所进行了一些堆芯物理计算和热工水力分析，而针对用于核电推进的球床反应堆，除提出相关概念设计外，并未进行深入研究。
技术实现思路
本专利技术公开一种空间核电推进系统球床反应堆物理热工耦合分析方法。该方法针对空间核动力系统中球床堆的结构特点，通过稳态热工水力计算，对堆芯的设计进行优化，对一些关键参数进行敏感性分析；通过瞬态计算，分析堆芯在一些典型事故工况下的响应特性，评估其安全性，为空间核电推进球床堆的设计及安全特性分析提供建议与指导。为了实现上述目的，本专利技术采取了以下技术方案予以实施：一种空间核电推进系统球床反应堆物理热工耦合分析方法，步骤如下：步骤1：将球床反应堆的堆芯根据需要沿轴向划分为若干个独立的通道，而后沿径向划分为若干控制体；输入球床反应堆的结构几何参数和边界条件，设定计算时间和时间步长；步骤2：根据球床反应堆的结构几何参数和边界条件，进行堆芯的稳态初始化计算，得到堆芯的热功率以及每个控制体的冷却剂质量流量、温度和压力，作为零时刻的初始值；堆芯的热功率：Pt＝Pn(1)第i个控制体的冷却剂质量流量：Mi＝Min(2)第i个控制体的冷却剂压力：Pi＝Pi+1+ΔPfi+ΔPai(3)第i个控制体的冷却剂温度：第i个控制体的燃料球表面温度：在堆芯的稳态初始化计算中，认为燃料球的各层材料温度等于燃料球的表面温度，即Tz1,i＝Tgm,i＝Tz2,i＝Thdc,i＝Tpyc,i＝Tu,i＝Ts,i(6)式(1)-(6)中：Pt——堆芯的热功率/WPn——堆芯的额定功率/WMi——第i个控制体的冷却剂质量流量/kg·s-1Min——入口冷却剂质量流量/kg·s-3Pi——第i个控制体的冷却剂压力/PaPi+1——第i+1个控制体的冷却剂压力/PaΔPfi——第i个控制体的摩擦压降/PaΔPai——第i个控制体的加速压降/PaTg,i——第i个控制体的冷却剂温度/TTg,i-1——第i-1个控制体的冷却剂温度/TAi——第i个控制体的流通面积/m2Ai-1——第i-1个控制体的流通面积/m2cp,g——冷却剂比容/J·kg-1·K-1ΔVi——第i个控制体的体积/m3q”'c——堆芯的功率密度/W·m-3hsg——对流换热系数/W·m-2·K-1asg——对流换热面积与体积比，Ts,i——第i个控制体的燃料球表面温度/TTz1,i——第i个控制体燃料球的ZrC涂层温度/TTgm,i——第i个控制体燃料球的活性区温度/TTz2,i——第i个控制体燃料球的TRISO颗粒ZrC涂层温度/TThdc,i——第i个控制体燃料球的TRISO颗粒致密石墨层温度/TTpyc,i——第i个控制体燃料球的TRISO颗粒疏松热解碳层温度/TTu,i——第i个控制体燃料球的TRISO颗粒UC-ZrC核心温度/T；步骤3：建立堆芯的中子物理模型、冷却剂热工水力模型、堆芯本体传热模型以及燃料球传热模型；堆芯中子物理模型，采用点堆动力学方程求解堆芯的裂变功率，方程考虑六组缓发中子和反应性反馈；式中：P(t)——t时刻堆芯裂变功率/Wt——时间/sρ(t)——t时刻总反应性/$βeff——缓发中子的总份额Λ——中子代时间/sλi——第i组缓发中子的衰变常数/s-1Ci(t)——t时刻第i组缓发中子先驱核的浓度/m-3βi——第i组缓发中子的份额冷却剂的热工水力模型：认为冷却剂在堆芯内的流动为不可压缩的径向一维流动，控制方程如下：质量控制方程：动量控制方程：能量控制方程：式中：ρg——冷却剂密度/kg·m-3r——径向坐标/mu——冷却剂径向速度/m·s-1P——冷却剂压力/PaS——动量源项cp,g——冷却剂比容/J·kg-1·K-1Tg,——冷却剂温度/TTs,——燃料球表面温度/Thsg——对流换热系数/W·m-2·K-1asg——对流换热面积与体积比，dp——燃料球直径/mε——球床孔隙率球床堆的堆芯本体传热模型，不考虑球床的轴向导热，球床的能量方程如下：式中：ρs——燃料球密度/kg·m-3cp,s——燃料球比容/J·kg-1·K-1λeff——球床的有效热导率/W·m-1·K-1q”'c——堆芯的功率密度/W·m-3燃料球的传热模型，燃料球包括活性区和ZrC涂层，活性区由众多TRISO型燃料颗粒弥散在石墨基质中组成；燃料球的ZrC涂层内没有内热源，其一维导热方程为：式中：ρz1——ZrC涂层密度/kg·m-3cz1——ZrC涂层比容/J·kg-1·K-1Tz1——ZrC涂层温度/Tλz1——ZrC涂层热导率/W·m-1·K-1燃料球的活性区看作是一种TRISO型燃料颗粒和石墨基质的均匀混合物，导热方程为：式中：ρgm——活性区密度/kg·m-3cgm——活性区比容/J·kg-1·K-1Tgm——活性区温度/Tλgm——活性区的有效热导率/W·m-1·K-1q”'gm——活性区的功率密度/W·m-3TRISO型燃料颗粒以UC-ZrC燃料为核心，由内向外依次有三层包覆层：疏松热解碳层、致密石墨层和外ZrC涂层；UC-ZrC核心的导热方程为：式中：ρ本文档来自技高网...

【技术保护点】
1.一种空间核电推进系统球床反应堆物理热工耦合分析方法，其特征在于：步骤如下：/n步骤1：将球床反应堆的堆芯根据需要沿轴向划分为若干个独立的通道，而后沿径向划分为若干控制体；输入球床反应堆的结构几何参数和边界条件，设定计算时间和时间步长；/n步骤2：根据球床反应堆的结构几何参数和边界条件，进行堆芯的稳态初始化计算，得到堆芯的热功率以及每个控制体的冷却剂质量流量、温度和压力，作为零时刻的初始值；/n堆芯的热功率：/nP

【技术特征摘要】
1.一种空间核电推进系统球床反应堆物理热工耦合分析方法，其特征在于：步骤如下：
步骤1：将球床反应堆的堆芯根据需要沿轴向划分为若干个独立的通道，而后沿径向划分为若干控制体；输入球床反应堆的结构几何参数和边界条件，设定计算时间和时间步长；
步骤2：根据球床反应堆的结构几何参数和边界条件，进行堆芯的稳态初始化计算，得到堆芯的热功率以及每个控制体的冷却剂质量流量、温度和压力，作为零时刻的初始值；
堆芯的热功率：
Pt＝Pn(1)
第i个控制体的冷却剂质量流量：
Mi＝Min(2)
第i个控制体的冷却剂压力：
Pi＝Pi+1+ΔPfi+ΔPai(3)
第i个控制体的冷却剂温度：



第i个控制体的燃料球表面温度：



在堆芯的稳态初始化计算中，认为燃料球的各层材料温度等于燃料球的表面温度，即
Tz1,i＝Tgm,i＝Tz2,i＝Thdc,i＝Tpyc,i＝Tu,i＝Ts,i(6)
式(1)-(6)中：
Pt——堆芯的热功率/W
Pn——堆芯的额定功率/W
Mi——第i个控制体的冷却剂质量流量/kg·s-1
Min——入口冷却剂质量流量/kg·s-3
Pi——第i个控制体的冷却剂压力/Pa
Pi+1——第i+1个控制体的冷却剂压力/Pa
ΔPfi——第i个控制体的摩擦压降/Pa
ΔPai——第i个控制体的加速压降/Pa
Tg,i——第i个控制体的冷却剂温度/T
Tg,i-1——第i-1个控制体的冷却剂温度/T
Ai——第i个控制体的流通面积/m2
Ai-1——第i-1个控制体的流通面积/m2
cp,g——冷却剂比容/J·kg-1·K-1
ΔVi——第i个控制体的体积/m3
q″′c——堆芯的功率密度/W·m-3
hsg——对流换热系数/W·m-2·K-1
asg——对流换热面积与体积比，
Ts,i——第i个控制体的燃料球表面温度/T
Tz1,i——第i个控制体燃料球的ZrC涂层温度/T
Tgm,i——第i个控制体燃料球的活性区温度/T
Tz2,i——第i个控制体燃料球的TRISO颗粒ZrC涂层温度/T
Thdc,i——第i个控制体燃料球的TRISO颗粒致密石墨层温度/T
Tpyc,i——第i个控制体燃料球的TRISO颗粒疏松热解碳层温度/T
Tu,i——第i个控制体燃料球的TRISO颗粒UC-ZrC核心温度/T
步骤3：建立堆芯的中子物理模型、冷却剂热工水力模型、堆芯本体传热模型以及燃料球传热模型；
堆芯中子物理模型，采用点堆动力学方程求解堆芯的裂变功率，方程考虑六组缓发中子和反应性反馈：






式中：
P(t)——t时刻堆芯裂变功率/W
t——时间/s
ρ(t)——t时刻总反应性/$
βeff——缓发中子的总份额
Λ——中子代时间/s
λi——第i组缓发中子的衰变常数/s-1
Ci(t)——t时刻第i组缓发中子先驱核的浓度/m-3
βi——第i组缓发中子的份额
冷却剂的热工水力模型，认为冷却剂在堆芯内的流动为不可压缩的径向一维流动，控制方程如下：
质量控制方程：



动量控制方程：



能量控制方程：



式中：
ρg——冷却剂密度/kg·m-3
r——径向坐标/m
u——冷却剂径向速度/m·s-1
P——冷却剂压力/Pa
S——动量源项
cp,g——冷却剂比容/J·kg-1·K-1
Tg,——冷却剂温度/T
Ts,——燃料球表面温度/T
hsg——对流换热系数/W·m-2·K-1
asg——对流换热面积与体积比，
dp——燃料球直径/m
ε——球床孔隙率
球床堆的堆芯本体传热模型，不考虑球床的轴向导热，球床的能量方程如下：



式中：
ρs——燃料球密度/kg·m-3
cp,s——燃料球比容/J·kg-1·K-1
λeff——球床的有效热导率/W·m-1·K-1
q″′c——堆芯的功率密度/W·m-3
燃料球的传热模型，燃料球包括活性区和ZrC涂层，活性区由众多TRISO型燃料颗粒弥散在石墨基质中组成；
燃料球的ZrC涂层内没有内热源，其一维导热方程为：



式中：
ρz1——ZrC涂层密度/kg·m-3
cz1——ZrC涂层比容/J·kg-1·K-1
Tz1——ZrC涂层温度/T
λz1——ZrC涂层热导率/W·m-1·K-1
燃料球的活性区看作是一种TRISO型燃料颗粒和石墨基质的均匀混合物，导热方程为：



式中：
ρgm——活性区密度/kg·m-3
cgm——活性区比容/J·kg-1·K-1
Tgm——活性区温度/T
λgm——活性区的有效热导率/W·m-1·K-1
q″′gm——活性区的功率密度/W·m-3
TRISO型燃料颗粒以UC-ZrC燃料为核心，由内向外依次有三层包覆层：疏松热解碳层、致密石墨层和外ZrC涂层；
UC-ZrC核心的导热方程为：



式中：
ρu——UC-ZrC核心密度/kg·m-3
cu——UC-ZrC核心比容/J·kg-1·K-1
Tu——UC-ZrC核心温度/T
λu——UC-ZrC核心热导率/W·m-1·K-1
q″′u——UC-ZrC核心的功率密度/W·m-3
疏松热解碳层、致密石墨层和外ZrC涂层内没有内热源，具有相同的导热方程：



式中：
ρ——包覆层密度/kg·m-3
c——包覆层比...

【专利技术属性】
技术研发人员：王成龙，张冉，秋穗正，田文喜，苏光辉，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61