反应堆压力容器下降段水-蒸汽CCFL实验与模型研究

为探究反应堆压力容器下降段在喷放末期冷段安注过程中的水-蒸汽逆流特性,建立下降段逆向流动限制(CCFL)模型,开展了基于压力容器模化本体的下降段CCFL实验研究以及建模分析。通过实验研究获得了不同入口安注水流量、安注水过冷度、堆芯蒸汽流量等条件下的下降段环腔内的安注特性数据,并基于实验数据进行了CCFL建模分析。结果表明,开始发生CCFL的蒸汽无量纲流速与入口安注水无量纲流速呈现正相关,基于无量纲流速建立的模型斜率与入口安注水无量纲流速呈现高度指数关联。本文建立了适用于从不发生CCFL至不完全CCFL,再到完全CCFL的下降段水-蒸汽气液逆流全过程预测模型。     

CAP1400非能动堆芯冷却整体试验关键现象分析

在CAP1400小破口失水事故中,非能动堆芯冷却系统所有设备均投入,显现复杂且独特的物理现象,为验证设计和程序以及识别重要现象,开展了CAP1400非能动堆芯冷却系统性能试验。整体试验台架主要特征为：1∶3高度比、最大工作压力9.2 MPa、等压模拟。试验结果表明,在小破口失水事故中堆芯不裸露,试验过程中发现了冷管段分层流产生机制、安注箱排空后氮气注射及其影响以及压力容器下降段流体温度不均匀性等关键物理现象。最后利用RELAP5程序对整体试验关键现象进行了分析和验证。     

“华龙一号”压力容器直接安注比例模化可视化试验研究

压力容器直接注入（DVI）技术可降低安注系统容量的要求,同时优化支持系统的配置,可以在保证安全的前提下简化安注系统设计、提高电厂经济性。本试验以“华龙一号”采用DVI技术优化的安注系统为研究对象,在模化比例为1:5的可视化模拟体上采用有色试剂跟踪法,以空气代替蒸汽,观察大破口失水事故（大LOCA）再淹没阶段DVI安注旁流特性,以及不同DVI管嘴结构对安注旁流的影响。发现大LOCA再淹没阶段安注旁流主要体现为直接旁流,在不带导流管嘴的4个DVI安注口的安注系统设计中,大LOCA再淹没初期的安注旁流约为7%,再淹没后期约为4%;带导流管嘴的DVI安注,再淹没初期的安注旁流约为2%,导流管嘴可有效降低安注旁流份额。本研究为“华龙一号”安全系统持续优化创新提供了重要参考。     

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

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

基于DVI管失水事故试验的CATHARE程序模拟评价

采用CATHARE程序对直接注入(DVI)管失水事故(LOCA)试验进行了数值模拟。研究发现:DVI管LOCA中系统卸压、非能动安注、堆芯冷却等主要过程和物理现象得到了较好的模拟。一回路系统压力、堆芯补水箱(CMT)安注流量、安注箱(ACC)安注流量、内置换料水箱(IRWST)安注流量以及堆芯流体温度等参数的计算结果和试验数据符合较好。研究结果表明,CATHARE程序可以用于失水事故下非能动安注系统瞬态特性模拟分析。     

基于ECC安注热混合试验的LOCUST 1.2分析与验证

在失水事故（LOCA）工况下安注系统投入使用时,蒸汽与安注冷却剂会发生流体热力学混合,热混合过程中冷腿段的冷却是直接影响堆芯再淹没与否的重要因素。中国广核集团有限公司自主研发了一款两相流热工水力系统分析软件LOCUST,可用于压水堆核电厂事故工况的分析计算。基于西安交通大学堆芯应急冷却系统（ECCS-XJTU）试验台架进行的堆芯应急冷却（ECC）安注热混合试验,本文使用LOCUST软件对ECC热混合试验进行了几何建模及计算分析。ECC热混合试验工况主要为不同流量下主管纯蒸汽与安注管过冷水的混合,蒸汽流量为25~125 kg/h,过冷水流量为100~500 kg/h。模拟计算结果和试验结果的对比分析表明：试验段出口质量流量计算值的最大相对误差在13.8%以内,混合后温度计算值的最大相对误差在8%以内,LOCUST在计算高温蒸汽和过冷水混合时的计算结果相对保守,总体上验证了LOCUST在LOCA下两相热混合安注计算的可靠性和准确性。     

气-液喷射器内两相流流型分析

对水平安装气-液喷射器内两相流流型进行了分析研究。根据圆管内气-液两相流流型转变的经验准则式,结合气-液喷射器性能方程,得出了气-液喷射器内环状流和雾状流的流型区间。以空气-水喷射器为例,给出了喷嘴出口当地最大马赫数与引射流体入口／出口压比的关系式,证明了当工作气体（空气）在喷嘴中作超临界膨胀时,喷射器内流型至少为环状流。     

无内热源有序饱和多孔介质通道内气液两相流阻力特性研究

对无内热源有序饱和多孔介质内蒸汽-水两相流阻力特性进行了实验研究,在多孔介质通道内蒸汽-水两相受力分析的基础上,结合实验数据,得到了多孔介质内空泡份额及气液两相相间作用关系式,通过分析热工水力特征参数和多孔介质几何特征参数对两相流阻力特性的影响,得到了多孔介质内蒸汽-水两相流阻力关系式。结果表明,本文提出的两相流阻力关系式计算结果与实验结果符合良好,且优于其他关系式。本研究结果为进一步开展含内热源多孔介质内气液两相流阻力及传热特性研究提供了实验技术及理论依据。     

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

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

溶液堆内气-液两相流流动及换热特性数值研究

在溶液堆台架模型数值模拟研究的基础上,对实际堆结构的堆芯内气-液两相流流动及冷却盘管与堆内溶液间的换热特性进行了数值模拟研究.采用欧拉两相流模型描述堆芯内气-液两相流流动,MUSIG(MUltiple-SIze-Group)模型描述堆芯内气泡尺度分布和相互作用,流-固耦合模型描述溶液与盘管间换热.数值计算得到了堆芯内的温度、速度、气泡组分等分布及冷却盘管的换热效率.数值计算结果表明:在有气泡扰动时堆内温度分布比没有气泡时均匀,冷却盘管可将堆内产生的85%热量带出,与试验测量结果一致.额定功率时,不同气体产生量对于冷却盘管换热影响的研究表明,随着堆内气泡产生量的增加,溶液与冷却盘管之间的换热得到强化.     

An Experimental Study of Thermal-hydraulic Phenomena in the Downcomer with a Direct Vessel Injection System of APR1400 during the LBLOCA Reflood Phase

The direct vessel injection (DVI) mode is adopted as a safety injection system in the place of a conventional cold leg injection (CLI) mode in the Advanced Power Reactor 1400MW (APR1400). It is expected that “sweep-out” and “direct ECC (Emergency Core Cooling) water bypass” are two most important bypass mechanisms of ECC water injected through the DVI lines during the LBLOCA reflood phase. Using the test facility of plane-channel type scaled down to 1/7 ratio of the prototype reactor (APR1400), we carry out the following tests for the investigation of the two mechanisms: water film spreading test, sweep-out test, and direct ECC water bypass test. From the water film spreading test, it was found that the curvature effect is negligible and the present modified linear scaling law is more appropriate than the linear scaling law. In the sweep-out test, the continuous onset is used to analyze the water height in the downcomer and the amount of ECC water bypass by sweep-out is compared with the previous correlations. The direct ECC water bypass test is performed in order to understand the ECC water film behavior in the downcomer.     

Effect of the yaw injection angle on the ECC bypass in comparison with the horizontal DVI

The comparison tests for the direct emergency core cooling (ECC) bypass fraction were experimentally performed with a typical direct vessel injection (DVI) nozzle and an ECC column nozzle having a yaw injection angle to the gravity axis. The ECC yaw injection nozzle is newly introduced to make an ECC water column in the downcomer region. The yaw injection angle of the ECC water relative to the gravity axis is varied from 0 to (±)90° stepped by 45°. The tests are performed in the air–water separate effect test facility (direct injection visualization and analysis (DIVA)), which is a 1/7.07 linearly scaled-down model of the APR1400 nuclear reactor. The test results show that (1) if the ECC water column is injected into the wake region which is induced by the hot leg blunt body in the downcomer annulus, the ECC bypass fraction is greatly reduced compared with the typical horizontal ECC injection which makes ECC film on the downcomer wall. At the same time, the ECC penetration toward the lower downcomer region becomes larger than those of a typical horizontal type of direct vessel injection on the downcomer wall vertically. (2) If the ECC water column is injected near the broken cold leg, the ECC water is directly bypassed. Thus, the ECC penetration fraction is greatly reduced compared with a typical film type of the horizontal ECC injection. (3) In order to minimize the ECC bypass fraction, the ECC water should be injected toward the wake region of the hot leg blunt bodies.     

Experimental validation of the modified linear scaling methodology for scaling ECC bypass phenomena in DVI downcomer

In the direct vessel injection (DVI) system downcomer, the direct emergency core coolant (ECC) bypass is activated during the reflood phase of a large-break loss-of-coolant accident (LBLOCA) by the interaction between the downward-flowing liquid-film and the transverse gas flow. Direct ECC bypass is reportedly the major bypass mechanism of ECC, and various experiments have been performed to obtain detailed information about the ECC bypass in a DVI downcomer. These lead to a proposed new scaling methodology, named ‘modified linear scaling’, which is expected to preserve the phase distribution in the downcomer and the ECC bypass phenomena. In the present study, modified linear scaling was experimentally validated in air–water tests comprising Test 21-D of the upper plenum test facility (UPTF). The counterpart tests of UPTF Test 21-D were performed in 1/7.3 and 1/4.0 scale models of a UPTF downcomer, and the test results were compared with the experimental data of UPTF Test 21-D. Two important parameters of direct ECC bypass – the normalized liquid-spreading width on the downcomer wall and the direct ECC bypass fraction, which is the fraction of input water bypassed out the broken cold-leg – were considered in the validation. The comparison revealed that the scaling parameters of direct ECC bypass are well preserved in the prototype and reduced models, from which we conclude that the modified linear scaling methodology is appropriate for designing a reduced test facility and for a scaling analysis of direct ECC bypass in the reflood phase of an LBLOCA.     

Air-water multidimensional impingement flow on vertical flat wall

The emergency core cooling (ECC) water is supplied from the direct vessel injection (DVI) system in the Advanced Power Reactor 1400 MWe (APR1400) during a postulated large-break loss-of-coolant accident (LBLOCA). The velocity of ECC water exceeds 10 m/s in the early high pressure phase of LBLOCA and then is decreased to 2-3 m/s in the late phase of reflood. During the injection the flow behavior exhibits a complex mode involving impingement, bypass, entrainment, sweepout and condensation in the reactor downcomer. There is currently no model to accurately simulate the local and complicated flow behavior in the APR1400 downcomer during a LBLOCA. This study is aimed at developing models for the water film flow and deformation, both of which are expected to sizably affect the other multidimensional flow characteristics in the downcomer. Experimental studies are conducted to benchmark the predictive model by furnishing the boundary conditions for the analysis resorting to the Accelerated Liquid Phase Hydrodynamics Apparatus (ALPHA) and the Kinetic Aerodynamic Physics Parallelepiped Apparatus (KAPPA). The Poisson equation and potential theory are applied to formulate the behavior of the water film and air flow. In both the experimental and numerical studies, the temperature-dependent thermodynamic properties and the reactor vessel curvature are neglected to render the problem at hand tractable. The model is found to reasonably describe the downward film flow behavior. The water film is developed in proportion to the initial injection velocity of the ECC water. The downward velocity of water film is increased with the heights of injection. Regarding the film deformation the calculated results tend to deviate from the experimental data as the injected air velocity is increased. The disagreement is attributed to limitations inherent in the two-dimensional treatment and point source approach.     

Investigation of Sweepout Mechanism and Critical Void Height in Annular Downcomer

The APR1400 (Advanced Power Reactor 1,400 MWe) has adopted the direct vessel injection (DVI) in lieu of the conventional cold leg injection for its emergency core cooling system (ECCS). In this reactor, sweepout from the water surface by gas (vapor or air) flow plays an important role in analyzing the mass and momentum transfer in the reactor downcomer of multidimensional geometry during a loss-of-coolant accident (LOCA) by decreasing the water level in the downcomer. The core water level will tend to decrease rapidly if a considerable amount of the entrained water stream and droplets bypasses through the break. The amount of entrained water is mostly determined by the interacting gas flow rate, the geometric condition, and the interfacial area between the gas and the water. The sweepout is observed to take place in three rather distinct steps: the beginning of undulation, the full wave and the wave peak (droplet separation). In view of these observations we investigated the relation between the gas flow rate and the amount of bypass as a function of time. The current experimental results shed light on the flow mechanism and the semi-empirical relations for the three-dimensional sweepout in a large-diameter annulus such as the reactor downcomer. A physico-numerical model is being developed to predict the multidimensional bypass flow rate resulting from the sweepout and entrainment in the downcomer.     

Direct vessel inclined injection system for reduction of emergency core coolant direct bypass in advanced reactors

Multidimensional thermal hydraulics in the APR1400 (Advanced Power Reactor 1400 MWe) downcomer during a large-break loss-of-coolant accident (LBLOCA) plays a pivotal role in determining the capability of the safety injection system. APR1400 adopts the direct vessel injection (DVI) method for more effective core penetration of the emergency core cooling (ECC) water than the cold leg injection (CLI) method in the OPR1000 (Optimized Power Reactor 1000 MWe). The DVI method turned out to be prone to occasionally lack in efficacious delivery of ECC to the reactor core during the reflood phase of a LBLOCA, however. This study intends to demonstrate a direct vessel inclined injection (DVII) method, one of various ideas with which to maximize the ECC core penetration and to minimize the direct bypass through the break during the reflood phase of a LBLOCA. The 1/7 scaled down THETA (Transient Hydrodynamics Engineering Test Apparatus) tests show that a vertical inclined nozzle angle of the DVII system increases the downward momentum of the injected ECC water by reducing the degree of impingement on the reactor downcomer, whereby lessening the extent of the direct bypass through the break. The proposed method may be combined with other innovative measures with which to ensure an enough thermal margin in the core during the course of a LBLOCA in APR1400.     

Advanced DVI for ECC direct bypass mitigation

An ECC direct bypass fraction during a late reflood phase of a LBLOCA is strongly dependent on the characteristics of the cross flow and the geometrical configuration of a DVI in the downcomer of a pressurized light water reactor. The important design parameters of a DVI are the elevation, the azimuthal angle, and the separator to prevent a steam-water interaction. An ECC sub-channel to separate or to isolate an ECC water from a high-speed cross flow is one of the important design features to mitigate the ECC bypass phenomena. A dual core barrel cylinder as an ECC flow separator is located between a reactor vessel and a core barrel outer wall in the downcomer annulus. A new narrow gap between the core barrel and the additional dual core barrel plays the role of a downward ECC flow channel or an ECC flow separator in a high-speed cross flow field of the downcomer annulus. The flow zone around a broken cold leg in the downcomer annulus has the role of a high ECC direct bypass due to a strong suction force while the wake zone of a hot leg has the role of an ECC penetration. Thus, the relative azimuthal angle of the DVI nozzle from the broken cold leg is an important design parameter. A large azimuthal angle from a cold leg to a hot leg needs to avoid a high suction flow zone when an ECC water is being injected. The other enhancing mechanism of an ECC penetration is a grooved core barrel which has small rectangular-shaped grooves vertically arranged on the core barrel wall of the reactor vessel downcomer annulus. These grooves have the role for a generation of a vortex induced by a high-speed cross flow. Since the stagnant flow in a lateral direction and rotational vortex provides the pulling force of an ECC drop or film to flow down into the lower downcomer annulus by gravity, the ECC direct bypass fraction is reduced when compared to the current design of a smoothed wall. An open channel of grooves generates a stagnant vortex, while a closed channel of grooves creates an isolated ECC downward flow channel from a high-speed lateral flow. In this study, new design concepts for a dual core barrel cylinder, grooved core barrel, and a reallocation of the DVI azimuthal angle are proposed and tested by using an air-water 1/5 scaled air-water test facility. The ECC direct bypass reduction performances of the new design concepts have been compared with that of the standard type of a DVI injection. The azimuthal angle of the DVI nozzle from a broken cold leg varies from −15° to +52° toward a hot leg. The test results show that the azimuthal injection angle is an effective parameter to reduce the ECC direct bypass fraction. The elevation of the DVI nozzle is also an important parameter to reduce the ECC direct bypass fraction. The most effective design for reducing the ECC direct bypass fraction is a dual core barrel. The reduction fraction when compared to the standard DVI is about −30% for the dual core barrel while it is −15% for the grooved core barrel.     

Numerical analysis and visualization experiment on behavior of borated water during MSLB with RCP running mode in an advanced reactor

The core bypass phenomenon of borated water injected through direct vessel injection (DVI) nozzles in APR1400 (Advanced Power Reactor 1400MWe) during main steam line break (MSLB) accidents with a reactor coolant pump (RCP) running mode has been simulated using a two-channel and one-dimensional system analysis model code (MARS), and a three-dimensional computational fluid dynamics (CFD) code (FLUENT). A visualization experiment has also been performed using a scaled-down model of the APR1400. The MARS analysis has predicted a serious core bypass phenomenon of borated water, while the CFD analysis has shown results opposite to the MARS results. The CFD analysis has shown that the flow pattern in the downcomer is fully three-dimensional and that vortex flow structures are formed near the cold legs so that the borated water might pass without difficulty into the high flow region of the cold legs and flow well into the lower downcomer. The visualization experiment has shown that the borated water flows well to the lower plenum, as in the CFD analysis. Both the CFD analysis and visualization experiment have proved that a serious core bypass phenomenon of borated water might not happen in the APR1400. These results are quite different from those predicted by MARS.     

Scaling for the ECC bypass phenomena during the LBLOCA reflood phase

As one of the advanced design features of the APR1400 (Advanced Power Reactor), a direct vessel injection (DVI) system is adopted instead of the conventional cold leg injection (CLI) system. It is known that the DVI system greatly enhances the reliability of the emergency core cooling (ECC) system. However, there is still a dispute on its performance in terms of water delivery to the reactor core during the reflood period of a large-break loss-of-coolant accident (LOCA). Thus, experimental validation is underway. In this paper, a new scaling method, using the time and velocity reduced “modified linear scaling law”, is suggested for the design of a scaled-down experimental facility to investigate the direct ECC bypass phenomena in the PWR downcomer.