承压热冲击下反应堆压力容器中流体的瞬态混合特性

热冲击下,反应堆压力容器中的热工水力特性是一个与反应堆安全密切相关的课题。本文在1/10的模型体上进行了高温高压下安全注水时流体的瞬态混合特性实验,得到了在有回路流动和无回路流动时以及不同的环腔流体温度下的混合特征。结果表明,环腔无流动时,随着安全注水流速的提高,混合函数下降得更快,幅度更大;环路有流动时,混合函数变化缓慢：当环腔内的流体温度达到一定的数值后,压力容器部分区域的混合函数发生明显变化。     

PTS瞬态流动与换热数值模拟研究

为了研究压水堆因“直接安注”冷水注入压力容器下降环腔而导致的承压热冲击（PTS）热工水力问题,基于1:10比例模型,应用计算流体力学软件FLUENT5.4进行了紊流流动换热的数值模拟分析,同时进行了常压瞬态传热实验研究。针对下降环腔折算流速0.5 m/s,安注流速10m/s的典型工况,研究了安注水开启后下降环腔内的瞬态流动换热特性,数值模拟与实验结果吻合良好。考察了压力容器安注接管出口区环形焊缝区及堆芯段筒体中子强辐照区所承受的热冲击状况,基于稳态流动研究了下降环腔内流体混合特性及流动机理,为热冲击分析提供参考。     

承压热冲击下AP1000压力容器直接安注瞬态数值模拟研究

基于计算流体动力学(CFD)分析方法,采用流固共轭传热方式,对非能动堆芯冷却系统(PXS)的堆芯补水箱(CMT)热态功能试验、CMT注入同时自动减压系统(ADS)动作、蓄压安注箱(ACC)安注后CMT再注入以及常规余热排出系统运行等4种工况下反应堆压力容器(RPV)环腔内流动传热状态进行瞬态数值模拟,研究RPV壁面温度瞬态变化以及环腔下降段内流体的混合特性。结果表明:4种工况下直接安注(DVI)接管管嘴与RPV内壁面相交斜面处冷却水混合剧烈,冷段是否有流体注入环腔对其内流体温度分布变化影响巨大,且DVI接管管嘴局部区域将发生较大的温度变化。     

稳压器波动接管内套管的温度保护作用

在西安交通大学动力工程多相流国家重点实验室水力试验台上进行了5：1 稳压器波动接管的瞬态温度模拟试验，实验研究了稳压器波动接管在密封套管和非密封套管两种结构形式，流体在正波动和负波动两种流动模式下，Re在4500－750000时，各测试点的瞬态温度变化，实验结果表明：密封套管具有良好的温度保护作用，能够有效保护波动接管免受热冲击和热疲劳作用，非密封套管的保护作用在低流速时微弱，高流速时基本消失。     

直接安注反应堆压力容器下降环腔内射流传热试验研究

失水事故工况 (LOCA)下反应堆下降环腔内的流动和传热研究 ,对反应堆压力容器 (RPV)的安全具有重要的意义。通过对一种直接安注的反应堆压力容器内流动和传热的研究 ,将流动分为横穿射流和冲击射流 ,比较了在两种射流下下降环腔内流动和传热的特点 ,分析了流速比和对流换热系数及温度的关系 ,当流速比在 1～ 1 0时 ,流动属于横穿射流 ,对流换热主要由环腔流速决定 ;流速比大于 1 0后 ,属于冲击射流 ,环腔内对流换热主要决定于安注流速 ,此时局部对流换热能力随安注流速的增加而增加     

AP1000小破口事故下非能动氮气安注箱的瞬态特性研究

根据AP1000非能动氮气安全注入水箱的结构和工作原理建立了热工水力模型并开发了计算分析程序TACAP。利用TACAP计算得到了AP1000非能动氮气安全注入水箱在两种小破口失水事故（包括25.4 cm等效直径冷管破口和5.08 cm等效直径冷管破口）下的瞬态特性,得到了箱内水位及注入流量等关键参数的瞬态变化。计算结果表明：安注箱在小破口失水事故后能提供高效的安全注入,对一回路快速地进行冷却和降压,有效地缓解事故后果。TACAP计算结果与西屋公司NOTRUMP程序计算结果基本一致,表明了TACAP程序的适用性和正确性。     

一回路一台泵停运-单环路余热排出工况下CEFR钠池三维瞬态热工特性数值模拟

一回路一台泵停运-单环路余热排出是池式钠冷快堆的设计基准事故之一,有必要对该工况下钠池内的热工特性进行分析。由于钠池整体尺寸大,难以开展实验研究,通常采用数值模拟的方法进行研究。因此,本研究基于计算流体动力学（CFD）方法,开展了该工况下CEFR钠池三维瞬态数值模拟,得到在一回路泵惰转、返流和非对称余热排出作用下钠池内三维瞬态流动、温度分布以及堆芯出口温度、中间热交换器（IHX）进出口温度等关键参数。计算结果表明,故障环路中泵、IHX存在返流现象。在900 s内,堆芯出口温度降至394.9℃。正常环路IHX出口温度在400 s左右达到最大值360.5℃,随后逐渐降低。故障环路IHX出口温度先下降后上升,900 s时接近364.3℃。具有余热排出的环路具有事故缓解能力,钠池整体温度没有明显升高。研究结果能够为一回路一台泵停运-单环路余热排出事故下池式钠冷快堆安全分析提供参考。     

基于主曲线方法研究加载速率对RPV材料断裂韧性的影响

利用3组不同材料预制的裂纹夏比试样（PCCv)研究了不同温度和不同加载速率对反应堆压力容器材料断裂韧性的影响，对采集到的实验数据用ASTM E1921—97标淮，计算出材料度的断裂韧性值和参考温度(T0)。从稳态到瞬态加载条件下的参考温度幅值用主曲线方法确定。研究结果表明，T0依赖于加载速率，并随加载速率的增大而增加，同时当稳态加载时的T0较小时，瞬态加载时的T0增值（△T0)较大。     

熔盐堆单参数扰动瞬态特性分析

熔盐堆（MSR）作为一种新型的反应堆,其热工水力特性与其他堆型有很大差异,扰动瞬态分析有助于从根本上了解其安全特性和运行状态。为了研究MSR的运行瞬态特性,本研究以液态燃料MSR为研究对象,利用经过修改的RELAP5/ MOD4.0程序进行了稳态运行工况下的扰动瞬态分析。干扰变量包括反应性引入、一回路熔盐质量流量、二回路质量流量、空气散热器质量流量、空气散热器入口空气温度。分析了主要运行参数,如功率、堆芯进出口温度、二回路进出口温度、特征时间等。结果表明MSR在各种扰动瞬态下的最终状态都趋于稳定,而不存在严重的瞬态变化,这是对其固有稳定性特性的直观表征。根据功率和温度等变量在扰动下的变化,提出了功率和不同回路温度的控制方法。      

承压热冲击下压水堆压力容器壁面换热特性的数值模拟研究

为了研究压水堆因安注冷水直接注入反应堆压力容器下降环腔而导致的承压热冲击（PTS）热工水力问题，基于1∶10比例模型，应用计算流体力学商用软件FLUENT5．4进行了紊流流动换热的数值模拟分析，同时进行了常压传热实验研究。针对下降环腔折算流速0．5m/s，安注流速10m/s的典型工况，研究了压力容器下降环腔的壁面换热特性。通过分析下降环腔内的流动及混合特性，从流动机理上解释了压力容器内壁上准重接触点附近换热强烈的现象，并指出壁面换热强弱与近壁流体紊流脉动动能密切相关，为热冲击分析提供参考。     

600MW压水堆安注箱设计研究

本文用美国核管会热工水力程序TRACE和图形化建模软件SNAP,建立了600 MW两环路压水堆一回路和二回路热工水力系统分析模型,并对安注箱的各设计方案进行大破口失水事故(LBLOCA)模拟计算,通过对比各设计方案在LBLOCA事故下计算出的峰值包壳温度,研究安注箱在大破口失水事故工况下的安注性能,最后给出了优化的设计方案,并提出了可行的设计改进建议。研究结果表明,上腔室和下降段同时注入的方式较冷段注入和下降段注入更有效,且恰当地选取初始安注箱压力,可有效降低峰值包壳温度,提高LOCA裕量。     

基于RELAP5的LOCA喷放阶段下降段内CCFL特性研究

反应堆失水事故（LOCA）后下降段通道内形成的两相逆流状态极有可能引发汽-液逆向流动限制（CCFL）,不利于应急冷却水顺利进入堆芯,极大影响了核反应堆系统的安全性能。本研究基于RELAP5程序采用Wallis溢流关系式对UPFT实验装置进行建模并计算LOCA喷放阶段的下降段注水行为;通过对比下腔室蓄水量、下降段内压力及破口处蒸汽流量瞬态变化以验证模型的有效性,并对下降段通道内汽相速度场、液相体积分数分布特性进行分析。结果表明,由于下降段通道结构的三维特征引起的流动不均匀性影响了汽-液CCFL特性,随着蒸汽流量增大,在破口环路与下降段连接区域的压力梯度与向上流速度梯度越大,较少节点的划分方法很难真实反映下降段通道局部区域内汽-液溢流关系;在靠近破口的环路内注入的冷却水更难到达下腔室,而在远离破口环路的冷却水容易进入到下腔室;过热的蒸汽在流动过程中被冷却水冷却发生凝结现象,导致出口蒸汽流量小于进口蒸汽流量,且随着进口蒸汽流量的增大,凝结效应则随之减小。本研究所建立的模型与方法能够适用于LOCA喷放阶段下降段通道内的汽-液CCFL预测。      

Boron dilution transients during natural circulation flow in PWR—Experiments and CFD simulations

Coolant mixing in the cold leg, downcomer and the lower plenum of pressurized water reactors is an important phenomenon mitigating the reactivity insertion into the core. Therefore, mixing of the de-borated slugs with the ambient coolant in the reactor pressure vessel was investigated at the four loops 1:5 scaled Rossendorf coolant mixing model (ROCOM) mixing test facility. In particular thermal hydraulics analyses have shown, that weakly borated condensate can accumulate in the pump loop seal of those loops, which do not receive a safety injection. After refilling of the primary circuit, natural circulation in the stagnant loops can re-establish simultaneously and the de-borated slugs are shifted towards the reactor pressure vessel (RPV).In the ROCOM experiments, the length of the flow ramp and the initial density difference between the slugs and the ambient coolant was varied. From the test matrix experiments with 0 resp. 2% density difference between the de-borated slugs and the ambient coolant were used to validate the CFD software ANSYS CFX. To model the effects of turbulence on the mean flow a higher order Reynolds stress turbulence model was employed and a mesh consisting of 6.4 million hybrid elements was utilized. Only the experiments and CFD calculations with modeled density differences show stratification in the downcomer. Depending on the degree of density differences the less dense slugs flow around the core barrel at the top of the downcomer. At the opposite side, the lower borated coolant is entrained by the colder safety injection water and transported to the core. The validation proves that ANSYS CFX is able to simulate appropriately the flow field and mixing effects of coolant with different densities.     

Numerical boron mixing studies for Loviisa nuclear power plant

A program is in the process of studying numerically boron mixing in the downcomer of Loviisa NPP (VVER-440). Mixing during the transport of a diluted slug from the loop to the core might serve as an inherent protection mechanism against severe reactivity accidents in inhomogenous boron dilution scenarios for PWRs. The commercial general purpose Computational Fluid Dynamics (CFD) code PHOENICS is used for solving the governing fluid flow equations in the downcomer geometry of VVER-440. So far numerical analyses have been performed for steady state operation conditions and two different pump driven transients. The steady state analyses focused on model development and validation against existing experimental data. The two pump driven transient scenarios reported are based on slug transport during the start of the sixth and first loop, respectively. The results from the two transients show that mixing is case and plant specific; the high and open downcomer geometry of VVER-440 seems to be advantageous from mixing point of view. In addition the analyzing work for the ‘first pump start' scenario brought up some considerations about flow distribution in the existing experimental facilities.     

CAP1400反应堆压力容器下降段内气-液逆向流试验研究

为研究真实工况下CAP1400反应堆压力容器下降段气-液逆向流现象,以CAP1400为原型,搭建压力容器下降段高度和直径比为1:1、60°切片的试验台架。试验工质为空气和水,试验研究了不同安注（DVI）供水量、不同气量的气-液两相流动和应急堆芯冷却剂（ECC）旁通现象。试验结果表明,DVI供水量相同时,随着供气量的增加,气-液逆向流现象明显,当质量流速达到4 kg/s及以上时,安注水不能全部进入堆芯;Kutateladze经验关系式和UPTF经验关系式都与试验结果存在较大偏差,不适用于CAP1400压力容器下降段试验;基于试验数据,拟合了新的经验关系式,且通过比较有无DVI挡块的试验数据,验证了DVI挡块可以降低ECC旁通水量,增强安注能力。      

Investigation of fluid mixing with direct vessel injection

The object of this work is to investigate fluid mixing phenomena as they related to pressurized thermal shock (PTS) in a pressurized water reactor vessel downcomer during transient cooldown with direct vessel injection (DVI), using test models. The test model designs were based on ABB Combustion Engineering (CE) System 80+ reactor geometry. A cold-leg, small-break loss-of-coolant accident (LOCA) and a main steam line break were selected as the potential PTS events for the ABB-CE System 80+. This work consists of two parts. The first part provides the visualization tests of the fluid mixing between DVI fluids and existing coolant in the downcomer region, and the second part presents the results of thermal mixing tests with DVI in the other test model. Flow visualization tests with DVI have clarified the physical interaction between DVI fluid and primary coolant during transient cooldown. A significant temperature drop was observed in the downcomer during the tests of a small-break LOCA. The measured transient temperature profiles compare well with the predictions from the REMIX code for a small-break LOCA, and with the calculations from the COMMIX-1B code for a stream line break event.     

高压安注条件下冷管段和环腔流体混合特性的数值分析

压水堆高压安注条件下冷热流体混合会导致承压热冲击现象,影响压力容器的使用寿命。本文基于ROCOM实验装置的实验数据,使用CFD方法对高压安注条件下有密度差的冷热流体混合现象进行了模拟,并对模拟结果进行了验证与分析。结果表明,在冷管段和下降段环腔中流体混合的主导因素分别为强迫流动混合和浮升力驱动混合。在仅有1条冷管段注入的情况下,进入下腔室的流体会再次回流至环腔,从而对冷却剂的混合特性产生影响。     

Multidimensional two-phase flow regime distribution in a PWR downcomer during an LBLOCA refill phase

The multidimensional countercurrent two-phase flow regimes that occur in a pressurized-water reactor (PWR) vessel downcomer during the refill phase of a large-break loss-of-coolant accident are studied using a transparent 1/10 scale model of a PWR vessel. The various flow regimes and their distribution in the downcomer have been identified and mapped for a range of air-water flooding experiments. The two-phase flow patterns that are identified in the downcomer included various types of film flows, droplet flows, countercurrent churn flows and cocurrent flows depending on the flooding condition. Through observation of the two-phase flow dynamics it was deduced that the physical mechanisms associated with the flooding processes could be separated into a liquid entrainment process and a film flow reversal process. In addition to the above exercise, the effect of non-uniform injection of water into the downcomer via different combinations of cold leg was studied similarly by determining flooding curves and flow pattern maps. It was found that differences in the flooding characteristic were noticeable for various water inlet configurations when compared with the uniform injection case. The differences could be explained qualitatively in terms of the flooding mechanisms identified previously by examining the flow patterns in the downcomer for the non-uniform injection tests.     

Boron mixing in complex geometries: flow structure details

Quantitative measurements of scalar dilution transients within a scaled model of a PWR downcomer are presented. The injection transients correspond to the case of a single pump start-up injecting into an initially stagnant system. Detailed results of the flow path and velocity field are presented for a nominal case where a neutrally buoyant front of dilute water is injected into the downcomer. Additional tests are made for a finite-sized slug of dilute water with a volume approximately equal to that of the downcomer and the lower plenum. The results of the tests indicate that the flow is dominated near the injection leg by an impinging jet, which is bifurcated into two primary streams by the cylindrical geometry of the downcomer. The two streams exit the downcomer on the opposite side of the injection leg. The splitting of the streams creates a bell-shaped recirculation region in the lower downcomer beneath the injection leg, which is the last volume to be displaced by the injected fluid flowing up from the lower plenum. A third plume sometimes forms at the boundary between the streams and the recirculation region, introducing mixed fluid into parts of this relatively stagnant area at earlier times than is typically observed. The results are discussed in the context of computational fluid dynamics codes that are starting to be used to simulate the mixing phenomena encountered in such complex systems.