蒸汽发生器一次侧流阻数值模拟研究.docVIP

蒸汽发生器一次侧流阻数值模拟研究 　　【摘 要】采用CFD的方法，分别采用k-ε湍流模型和SST湍流模型对蒸汽发生器一次侧流动阻力进行数值模拟，获得了蒸汽发生器一次侧入口段、U形管段及出口段的压降，并将CFD计算结果与按经验公式计算结果进行对比分析，提出了对经验公式修正建议。 　　【关键词】蒸汽发生器；一次侧流动阻力；数值模拟 　　0 前言 　　蒸汽发生器是压水堆核电站的关键设备，是核电厂一、二次侧的重要枢纽。在压水堆核电厂中蒸汽发生器一次侧阻力约占反应堆冷却剂系统总流动阻力的40%左右[1]，因此，准确的计算蒸汽发生器一次侧阻力对反应堆冷却剂系统的水力计算和系统布置具有重要意义。 　　在蒸汽发生器的工程设计中，通常采用专用热工水力程序进行蒸汽发生器的水力计算，如岭澳核电站55/19型蒸汽发生器采用GVSPA程序进行蒸汽发生器一次侧流阻计算，秦山二期核电站60F蒸汽发生器采用GENF程序进行一次侧流阻计算。在上述专用程序中，流动阻力采用经验公式计算，利用经验公式手册查得对应结构的流动阻力关系式，通过简化完成压降计算。这种计算方法需对应于特定结构，针对性较强。目前随着CFD技术的成熟，有学者采用CFD手段开展了蒸汽发生器一次侧流动特性研究[2]。 　　本文针对AP1000核电站蒸汽发生器的结构特点，采用经验公式方法和CFD方法对蒸汽发生器一次侧流阻进行了计算，并对两种方法的计算结果进行了分析和讨论。 　　1 经验公式法计算 　　针对AP1000蒸汽发生器的结构特点，蒸汽发生器一次侧流动阻力可分为5个部分进行计算，如图1所示。其中包括一次侧入口接管突扩局部阻力K1、一次侧出口接管突缩局部阻力K5、传热管入口突缩局部阻力K2、传热管出口突扩局部阻力K4和传热管阻力K3。 　　蒸汽发生器流阻计算相关参数可由蒸汽发生器稳态热工水力计算获得。 　　采用经验公式计算时，选用不同的摩擦阻力系数计算公式，摩擦阻力系数略有差别，但最大差值为3.1%，在工程计算误差可接受范围内。 　　2 数值模拟计算 　　2.1 几何模型 　　由于蒸汽发生器传热管数量众多，且传热管长度较长，横向和纵向尺度差别大，如果对蒸汽发生器一次侧整体进行数值模拟，则会由于网格数目过于巨大，难以实现，因此本文将蒸汽发生器分为进口部分、传热管部分和出口部分，分别进行模拟，最后将结果相加获得总压降。 　　图2为进口部分流体区域，与进口部分相连接的传热管截取高度为150mm，内径为15.4mm，这样高度大约为内径的10倍，可以保证流体充分发展而不受进口条件的影响。出口部分流体域的确定采用类似方法，传热管截取高度为150mm。 　　2.2 网格划分 　　本文利用Ansys中的ICEM CFD进行网格划分，对于变径面以及弧面进行了局部网格加密，对于管束区的网格通过网格拉伸完成，对和流体流向不垂直的壁面添加了附面层网格。进口部分总网格数为4986万，出口部分总网格数为2400万。 　　工程实际中传热管的数量巨大，难以按照实际情况进行模拟，传热管沿程压降利用最短一根传热管和最长的一根传热管进行数值模拟，通过求取二者压降的平均值作为U型传热管的压降损失。最短的U型管的网格总数为274300，最长的U型管的网格总数为386960。 　　2.3 边界条件 　　蒸汽发生器发生器一次侧进口部分水的温度和压力基本相同，且下封头不存在热交换，因此可采用绝热模型，流体区域内水的物性参数根据进口部分水的物性参数确定。进口流速为v=23m/s，温度T=321℃。 　　蒸汽发生器发生器一次侧出口部分同样可采用绝热模型，流体区域内水的物性参数根据出口处水的物性参数确定。 　　蒸汽发生器传热管需要与二次侧进行热交换，U形换热管数量巨大，导致换热量较大，但是单根传热管内冷却剂的温度并没有发生很大改变，因此传热管内冷却剂的物性参数根据进出口处水的平均温度给定，根据传热管内平均流速给定入口流速。 　　2.4 计算模型 　　本文将蒸汽发生器一次侧分为3部分进行数值模拟，单独考虑进口部分或出口部分，流体温差不大，且为单相液体流，水的密度近似为常数。因此，蒸汽发生器一次各部分侧流场可视为三维定常不可压粘性等温流场。湍流模型分别选用k-ε湍流和SST湍流模型。本计算中的压力为相对压力。 　　2.5 计算结果及分析 　　2.5.1 局部流动特性 　　从计算结果发现，采用k-ε和SST这两种湍流模型对计算结果影响不大，两种计算模型的总压降差值仅约为3.2%，属于工程计算误差可接受范围，下文以k-ε模型的结果为例进行说明。 　　蒸汽发生器一次侧进口部分总压力分布和速度矢量图如图3和图4所示。由图3的进口总压力分布可以看出下封头的进口接管至传热管的压降较大，说明进口部分阻力较大