申请/专利权人：西安交通大学

申请日：2021-12-30

公开（公告）日：2024-04-09

公开（公告）号：CN114330169B

主分类号：G06F30/28

分类号：G06F30/28;G06F113/08;G06F119/08;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.04.09#授权;2022.04.29#实质审查的生效;2022.04.12#公开

摘要：本发明公开了一种小型模块化超级安全气冷堆物理热工耦合分析方法，步骤如下：1、划分控制体，输入小型模块化超级安全气冷堆的结构参数和边界条件，给定不同控制体的群常数，设定计算时间和时间步长。2、堆芯初始化计算，得到零时刻的初始值。3、采用PIMPLE算法，计算得到当前t时刻的堆芯功率分布、流场、温度场和压力场。4、在求解每个方程前，计算初始残差，若小于用户设定的公差，则方程不进行求解。当所有方程的初始残差小于公差时，采用PIMPLE进行下一时刻的计算，直到设定的总计算时间，计算停止。本发明可以快速分析小型模块化超级安全气冷堆堆芯的稳态热工水力特性及事故工况下的瞬态响应特性，为气冷堆的设计及安全特性分析提供建议与指导。

主权项：1.一种小型模块化超级安全气冷堆物理热工耦合分析方法，其特征在于：步骤如下：步骤1：将小型模块化超级安全气冷堆的堆芯根据不同组件以及不同功能划分为若干控制体；输入小型模块化超级安全气冷堆的结构几何参数和边界条件，给定不同控制体的群常数，设定计算时间和时间步长；步骤2：根据小型模块化超级安全气冷堆的结构几何参数和边界条件，进行堆芯的稳态初始化计算，得到堆芯的热功率以及每个控制体的冷却剂质量流量、温度和压力，作为零时刻的初始值；步骤3：建立堆芯的冷却剂热工水力模型、中子物理模型以及燃料棒导热模块；堆芯的冷却剂热工水力模型，为实现对小型模块化超级安全气冷堆堆芯的快速模拟，采用多孔介质热工水力模型，求解氦气在堆芯内的可压缩流动；多孔介质热工水力模型的控制方程如下：质量守恒方程： 式中：γ——孔隙率ρ——冷却剂密度kg·m-3t——时间su——冷却剂流速m·s-1动量守恒方程： 式中：μT——湍流运动粘度Pa·sp——冷却剂压力Papi——冷却剂与亚尺度结构界面间的压力PaFg——重力引起的单位体积阻力N·m-3Fss——亚尺度结构相互作用引起的的单位体积阻力N·m-3能量守恒方程： 式中：e——冷却剂总能量J·kg-1kT——湍流热导率W·m·K-1Qss——亚尺度结构对流体的导热量W·m-3假设p＝pi，并定义则上三式可简化为：质量守恒方程： 动量守恒方程： 能量守恒方程： 此时，孔隙率γ的导数仅出现在动量与能量方程中等式右侧的最后一项；在核工程领域的应用当中，孔隙率在多孔介质区域内不发生变化，因此将其导数项忽略不计；动量守恒方程中，Fss代表各控制体对流体的阻力，与速度成正比； 式中：fD——达西摩擦因数 ——常数项Dh——水力直径能量守恒方程中，Qss代表流体与控制体之间的热交换； 式中：AV——流体与亚尺度结构间的单位体积面积m2·m-3h——换热系数W·m-2·K-1Tss——亚尺度结构温度T堆芯的中子物理模型，采用多群扩散中子动力学方程求解堆芯的裂变功率，对中子能量进行分群，根据控制体内的中子数守恒原理，即某区域内中子通量密度的变化率等于该区域内中子的产生率减去消失率，得到如下的多群扩散中子动力学方程；其中群数与群常数由用户指定，方程考虑缓发中子和反应性反馈；多群扩散中子动力学方程如下： 其中，源项Sn,i,Sd,j和SS,i分别为： 式中： ——能群i的中子通量m-2·s-1vi——能群i的平均中子速度m·s-1D——中子扩散系数mβeff,t——总缓发中子占比χf,i——瞬发中子产额χd,i——缓发中子产额keff——有效增值因子λk——第k组缓发中子衰变常数s-1燃料棒导热模块，采用有限体积法求解得到石墨基质、燃料温度；取小型模块化超级安全气冷堆的一个最小单元，将其近似看作多层圆柱体导热，由外到内依次是燃料层、石墨层、氦气层；假设每个网格内存在一个燃料棒单元，使用有限体积法计算得到燃料温度、石墨基质温度；首先读取堆芯的冷却剂热工水力模型、中子物理模型以及燃料棒导热模块参数，包括燃料层和石墨层的内径、外径、密度、定压比热容、热导率、分层数；径向温度分布公式： 式中：q——功率密度W·m-3r——半径mk——热导率W·m-1·K-1C1、C2——常数项多层圆柱导热微分方程： 式中：cp——定压比热容kJ·kg-1·K-1由径向温度分布公式计算得到多层圆柱体各层的温度Trad用于后续计算；对多层圆柱体微分方程使用有限体积法进行离散，建立线性方程M·Trad＝S进行温度分布迭代；其中，矩阵M为： 矩阵S为：S＝qdA+Told·XdA16S＝Told·XdA+2πrh最内侧点17B、C、X均由边界条件计算得到： X＝ρcpdt20式中：Told——各节点上一次迭代的温度TΔr——各节点间距m采用PIMPLE算法，依次求解质量守恒方程、k-ε方程、动量守恒方程、燃料棒导热模块中方程、能量守恒方程、多群扩散中子动力学方程，得到当前t时刻的堆芯功率分布、流场、温度场和压力场；步骤4：在求解每个方程前，计算初始残差，若小于用户设定的公差，则方程不进行求解；当所有方程的初始残差小于公差时，采用PIMPLE算法进行下一时刻的计算，直到设定的总计算时间，计算停止。

全文数据：

权利要求：

百度查询： 西安交通大学 一种小型模块化超级安全气冷堆物理热工耦合分析方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD