NB_T 20462-2017压水堆乏燃料干法贮存设施热工分析.pdf

ICS 27. 120. 20F 65备案号：59655—2017NB中华人民共和国能源行业标准NB/T20462—2017压水堆乏燃料干法贮存设施热工分析Thermal analysis of spent fuel dry storage installationforpressurizedwater reactor2017-04-01发布2017-10-01实施国家能源局发布 NB/T 204622017A.5物理模型A.5.1计算域A5.1.1以立式混凝土筒仓忙存容器为例，可将混凝土筒仓和金属忙建之间的空气流道（以下简称空气域，见图A.3）和金属贮罐内的密闭氧气域（以下简称氨气域，见图A.4）设置为独立的流体域。A5.1.2空气域和氢气域为浮力驱动的自然对流，由于温度变化较大，建议采用理想气体状态方程描述空气和氮气的热力学状态。热工分析时应选取合理的参考压力和参考密度，以减少数值计算的误差。A5.1.3对于燃料组件区域，采用多孔介质进行等效处理时，需要设置合理的孔隐率、渗透率和阻力系数,用于计算气流经乏燃斜组件时的压降AP。图A.3空气域的自然对流图A.4氮气域的自然循环 NB/T 20462—2017A. 5. 2流动状态由于乏燃料组件衰变余热，空气域和氮气域的壁面上存在较大温差，从而产生密度驱动的自然对流。选取热工分析的端流模型前，需要对计算域内的流动状态进行判定。通常可采用格拉晓夫数Gr判断自然对流的流动状态。Gr = gα,4773v2...(A.)式中：g-重力加速度：(,体积热膨账系数：A7—温度差：特征长度：V——运动黏度。对于竖圆柱形状的空气域和氨气域：—Gr为1×10~3×10°时，空气域和氨气域内的流动为层流：Gr为3×10～2×10时，空气域和氢气域内的流动为过渡流：Gr2X10时，空气域和氢气域内的流动为端流。A. 5. 3满流模型确定空气域和氮气域的流动状态后，应选用合理的流体分析模型。空气域可采用k-模型，氢气域可采用层流模型或k-t模型。A.5.4辐射模型A.5.4.1在乏燃料干法贮存设施可能出现的温度范围内，空气、氮气发射和吸收辐射的能力非常弱，可近似为热辐射的透明体，热工分析中可以忽略空气和氨气的辐射和吸收特性，并在此基础上选用合适的辐射换热模型。A.5.4.2采用FLUENT程序计算辐射传热时可采用分散坐标模型（DiscreteOrdinates，简称DO模型）或分散传输辐射模型（DiscreteTransferRadiationModel，简称DTRM模型）；使用CFX程序计算辐射传热时可采用分散传输模型（DiscreteTransferModel,简称DTM模型）。A.5.4.3对于燃料组件区域，由于辐射传热已通过等效导热系数入r加以考虑，因此不再考虑该区域内的辐射传热。A.5.5衰变热源热工分析时，乏燃料的衰变热作为热源设置在乏燃料组件区域内，由于乏燃料组件燃耗深度的不同。每个乏燃料组件的衰变热功率可能存在差异，应根据实际装载方案确定每个乏燃料组件的衰变热功率。衰变热应设置在燃料组件的活性段区域（即燃料芯块区域），通常活性段中间区域的燃耗最大，设置衰变热时应考虑沿燃料组件轴向的功率分布。A5.6边界条件以立式混凝土筒仓贮存容器为例，热工分析时，需设置以下的边界条件：通风孔进口，采用压力进口边界，进口压力为环境压力，进口温度为环境温度；通风孔出口，采用压力出口边界，出口压力为环境压力8 NB/T 20462—2017混凝土筒仓顶部和侧面，综合考虑自然对流换热和辐射换热：对于自然对流换热，可采用如下的关联式：Nu = C(GrPr)*..(A.2)式中;Nu-努塞尔数：格拉晓夫数；Pr-普朗特数。对于几种典型的表面形状及其布置情况，由实验确定的常数C和n值如下表A.1所示，表A1大空间自然对流关联式的常数C和n加热表面系数C及指数形状与位置流动情况示意流态Gr数适用范围c层流0.591/410′~3×10°竖平板及聚园柱过渡流0.02920.393×10°~2×1010油流0. 111/3 2×10*层流0.481/410°~5.76×108横园柱过渡流0.0445L05.76×10*~4.65×10*0.101/34.65×10a对于混凝土简仓底部可采用绝热边界进行保守考虑，也可设置对流换热系数，对流换热系数可按照3m厚度土壤的热传导进行等效计算，A.6求解控制A. 6.1离散格式为保证计算精度，选择离散格式时应尽量避免使用一阶离散格式，建议采用二阶及以上精度的离散格式，如使用FLUENT时，压力项可选择二阶精度（SecondOrder）、动量项和端动能项可选择二阶迎风格式(SecondOrderUpwind)或更高精度的设置。使用CFX时，对流项和端流数值项可采用高精度