稳压器波动管流场及应力数值模拟.docVIP

精品论文 参考文献 稳压器波动管流场及应力数值模拟 （山东核电有限公司 山东烟台 265116） 摘要:在核电厂中,稳压器波动管是调整和维持反应堆一回路压力的重要管道.稳压器中的液相水与主管道中的冷却剂存在较大的温差,当发生正波动或负波动时,就会出现热分层现象.这种温度在空间上的非均匀分布和时间上的非定常性会带来波动管热变形在径向和整个管道长度上的不均匀,波动管管道产生静态热应力或振荡热应力,进而导致波动管产生整体上的弯曲变形、截面变形和局部的热应力集中,产生非预期的大位移和支撑载荷,进而影响波动管的运行可靠性和寿期,因此对热分层导致的波动管热应力分布特性进行分析研究是十分必要的. 本文根据实际建立了波动管的物理计算模型，运用Fluent和Ansys软件对典型工况进行了数值模拟，分析了其热分层现象和热应力分布现象，找到了在一定条件下的热分层规律和热应力分布和形变规律，为压水堆核电厂的运行和改进提供参考。 关键词:稳压器波动管;热分层;应力分析;形变 1．引言 一回路系统在稳态运行或变工况运行时,各种扰动因素或负荷变化会使冷却剂温度发生变化,并且当一回路上冲下泄出现不平衡（或泄漏）时,又会使冷却剂容积发生变化.在封闭的一回路系统中任何温度和容积的变化都会影响到压力,如果压力过高,会危及设备的安全;压力低于额定压力过多时,则反应堆堆芯内产生大量沸腾,有导致燃料熔化的危险.这些压力都是通过稳压器进行调节和控制的.而波动管是连接稳压器和一回路某一环路热管段的唯一管路,在压力控制中起着传输通道的作用. 当主管道压力偏高时,主冷却剂从主管道热段经波动管流向稳压器;当主管道压力偏低时,主冷却剂从稳压器经波动管流向主管道热段.由于主冷却剂与稳压器存在较大的温差,不同温度下的流体比重不同,当波动管流量较小时在水平管段,两种 不同温度流层之间缺少混合,来自主回路的冷却剂温度低,密度高,来自稳压器的冷却剂温度高,密度低,热的介质由于热浮力的影响,将向上运动占据管线的上方,这时两层流体之间被一薄的边界层分开,这种现象叫做热分层.与此同时,在冷热流体的交界面,由于一股缓慢流动流体的干扰,会形成象表面波一样的震荡. 热分层在管道中会带来两个方面的主要影响,即导致波动管产生整体上的弯曲变形、截面变形和局部的热应力集中,产生非预期的大位移和支撑载;由反复出现的热分层和热振荡引起的总体应力和局部应力使管道的疲劳寿命进一步降低.因此美国核管会于1988年两次颁发公告(NRC, Bulletin No．88-08＆88-II)[1,2]要求对核电站中有温差存在的管线系统进行热分层研究,计算热应力,进而对管线的安全设计作出更合理的评估.国内关于稳压器波动管热分层的研究基本上处于空 白.虽然国外有关这方面的报道较多,但具体的研究方法很少报道. 本文目利用商业软件Fluent 和 Ansys 采用数值模拟的方法对稳压器波动管及其连接部件内冷却剂进行流场计算,对热分层现象进行了分析,并对波动管应力即形变进行数值模拟.得到了不同参数条件下稳压器波动管内的流场分布特性,稳压器波动管内及管壁的温度场分布特性及稳压器波动管的应力分布和形变特性.将有助于了解不同参数和运行工况条件下波动管热分层引起的静态或动态热应力特性,为新型核反应堆设计中波动管的改进和反应堆运行提供参考. 2计算模型 2.1流体瞬态时均方程 连续方程(质量方程) 2.2 RNG k-epsilon;模型[3] 求解流动及换热问题时,采用目前最广泛应用的k-epsilon;模型,但是由于标准的k-epsilon;模型用于强旋流或带有弯曲壁面的流动时会出现一定的失真,我们选用了RNG k-epsilon;模型,该模型可以更好地处理高应变率及流线弯曲程度较大的流动. 湍动能k和湍动能耗散epsilon;的输运方程为: 2.3热分层发生条件与热浮力计算模型热分层发生条件 热分层现象的发生取决于浮力与流体惯性力的比值.浮力越大或惯性力越小,热分层现象就越容易出现.为了表示这种关系,定义无量纲Richardson数[4] 由上式可以看出,波动管发生热分层现象的3个主要条件是:①有一段水平或近似水平的管段;②稳压器和主管道热段内冷却剂有显著温差;③波动管内流体流动速度应较低.[5] 对于包含浮升力的计算,动量方程增加一个源项.这个附加的源项通常有两种计算方法——全浮力模型和波斯纳斯克假设模型.波斯纳斯克假设模型由于附加的