脉动流下棒束通道内相位差及瞬态流场研究

事故条件及海洋条件下反应堆处于非稳态工况,堆芯燃料组件内热工水力行为具有瞬变及多因素耦合特性,对反应堆的安全提出更高挑战,因此有必要对燃料组件内瞬态特性进行研究。本文通过测量棒状燃料组件内压降和流量之间延迟时间开展棒束通道脉动流条件下相位差研究,对比了相位差在不同振幅、不同流动状态下的变化特性,并分析了定位格架对脉动流相位差的作用特点。另外,基于粒子图像测速（PIV）技术开展了脉动流条件下棒束通道内流场分布特性研究,对比了相同流量条件下稳态工况与瞬态工况下流场分布差异,分析了主流具备不同加速度时棒束通道内流场分布特征。实验结果表明：定位格架可减小脉动流下棒束通道内相位差;棒束通道内流场演化滞后于主流量变化。实验结果有助于揭示燃料组件在非稳态条件下瞬态特性,并为燃料组件的设计和优化奠定基础。     

Impact of gas in sodium flow on the temperature variation in an LMFBR rod bundle

The impact of gas in sodium flow on the temperature variation in an LMFBR rod bundle was studied in two types of experiments: (1) The gas fraction of the subchannels as well as the gas bubble spectra across the outlet of an unheated 61-rod bundle with wire spacers were measured in water/air flow. The distributions of the gas fractions at the inlet of the bundle were performed under uniform and non-uniform conditions. The results show that the distribution of the averaged gas fractions between the individual subchannels at the outlet of the bundle was almost the same as the distribution at the inlet. The measured bubble spectra show a dependency existing between the bubble frequencies, the bubble lengths, and the gas fraction in a subchannel. (2) A model for computing the transient temperature distributions within a heated rod was supported by experiments performed in a sodium/argon flow. For slug flow conditions a comparison indicates that the measured variations of wall temperatures can be well interpreted as being functions of the bubble contact time, rod power, and gas fraction in the flow.     

自然循环系统摇摆条件下棒束通道内传热特性研究

棒束燃料元件子通道间流体存在搅混与横向二次流,流动及阻力特性相较矩形通道、圆管等简单通道更为复杂。核动力舰船、船舶、小型浮动核电站等会受到海浪影响,经常处于倾斜、摇摆、垂荡等瞬变运动下。目前的相关研究多集中在低压工况的研究领域,高温高压自然循环运动条件下的研究较少。本文采用实验研究方法,对自然循环系统摇摆条件下棒束通道内流动传热特性进行了研究,获得了过冷沸腾和饱和沸腾两种条件下摇摆角度和摇摆周期对棒束壁面温度变化和传热系数的影响,并获得了摇摆周期内棒束通道内的传热系数计算关系式。结果表明,饱和沸腾传热系数变化比过冷沸腾的剧烈;在本文实验工况范围内,棒表面传热系数波动幅值随着摇摆幅度的增大而增大;摇摆条件下棒束通道过冷沸腾和饱和沸腾工况时均传热系数基本不变。     

Development of a Transient Boiling Transition Analysis Method Based on a Film Flow Model

A new single-channel, transient boiling transition (BT) prediction method based on a film flow model has been developed for a core thermal-hydraulic code. This method could predict onset and location of dryout and rewetting under transient conditions mechanically based on the dryout criterion and with consideration of the spacer effect. The developed method was applied to analysis of steady-state and transient BT experiments using BWR fuel bundle mockups for verification. Comparisons between calculated results and experimental data showed that the developed method tended to predict occurrence of rewetting earlier, however, onset time of BT and maximum rod surface temperature were well predicted within 0.6 s and 20°C, respectively. Moreover, it was confirmed that consideration of the spacer effect on liquid film flow rate on the rod surface was required to predict dryout phenomena accurately under transient conditions.     

Slug to annular flow transition during boiloff in a rod bundle under high-pressure conditions

High-pressure boiloff experiments in a wide range of bundle powers by using the Two-Phase Flow Test Facility (TPTF) were conducted. Two kinds of boiloff patterns were observed in these experiments. One is the boiloff pattern in a low bundle power, in which the dryout points of rods locate at a certain elevation in the bundle because the mixture level controls the dryout points. The other is the boiloff pattern in a high bundle power, in which the clear mixture level can not be observed and the dryout points of rods locate in a wide range of vertical directions. The vertical scatter of the dryout points is considered to be due to the break of the thin water film on the heater rods under the annular flow pattern.A simple model to predict the slug to annular flow transition in the rod bundle is proposed. In the model, the slug to annular flow transition takes place when the interferences of the water films on the neighboring rods cease. The model appeares to give good predictions of the previous flow transition experiment conducted in a rod bundle. The slug-annular transition below the dryout points was predicted with the present model in the high power boiloff experiments of TPTF. No slug-annular transition below the dryout points is predicted with the present model in the low power boiloff experiments.     

加/减速流动下棒束通道内速度分布和湍流特性研究

事故工况及海洋条件下反应堆处于非稳态工况,堆芯燃料组件内热工水力行为复杂多变,对反应堆安全提出了更高挑战,因此有必要对非稳态下燃料组件内流动换热特性开展研究。基于粒子图像测速（PIV）技术,结合远心镜头和脉冲控制器,实现对燃料组件内复杂流场的高时空分辨率、长时间的连续测量,获得了流量波动下燃料组件内时空演变的流场结构,分析了棒束通道内速度分布、湍流强度、雷诺应力等瞬时流场信息的空间演变特性。以定常流动下流场分布特性为基准,对比分析了加速度对燃料组件内空间流场分布的贡献特点。实验结果表明：加速流动提高了棒束通道内流层之间的速度梯度,抑制了横向速度和湍流强度;减速流动减弱了棒束通道内流层之间的速度梯度,提高了横向速度和湍流强度。实验结果有助于揭示燃料组件在非稳态条件下的瞬态特性,并为燃料组件的设计和优化奠定基础。     

Comparisons of the RALIZA-2/02 two-phase flow model with experimental data

The RALIZA-2 computer program was designed for thermal-hydraulic analysis of flow channel and fuel element of PWR/BWR at steady-state and transient conditions. A nonhomogeneous, nonequilibrium model of a two-phase flow and a two-dimensional heat conduction model of fuel pin are used in the program. A fully implicit integration scheme for both models is used. The steady-state constitutive correlations set is used. The void fraction, pre- and post-DNB heat transfer mechanism are compared with data. Also a loss of flow experiment was calculated and compared with nuclear heated rod bundle experimental data for typical PWRs. A very good agreement was obtained.     

数值堆热工流体程序CVR-PACA验证及典型应用

谱元方法是一种高精度的数值计算方法,采用该方法开发了数值堆高精度热工水力并行CFD计算程序CVR-PACA。应用CVR-PACA对单棒光棒通道湍流流场、3×3光棒棒束湍流流场、Matis-H压水堆棒束通道基准题、19棒带绕丝组件通道湍流流场进行了仿真计算。通过与实验测量值对比,研究定量验证了大涡模拟(LES)模型及非稳态雷诺时均(URANS)模型对各类棒束通道流场预测的准确性。算例在建模过程中采用网格分裂技术实现了复杂几何的纯六面体网格划分,用于支撑谱元方法计算。研究较为全面地积累了高精度谱元方法模拟流场流动及换热的建模经验,获取了各类棒束通道内丰富的流动和换热细节,获得的建模经验能更加精准有力地指导相关设计的优化改进。     

冷却剂与含格架燃料棒束多场耦合分析

燃料棒束作为压水堆燃料组件的组成部分,其热工和结构特性直接关系到反应堆的安全。本文利用ANSYS WORKBENCH软件分析了冷却剂在5×5含定位格架燃料棒束通道内流动的分布,采用冷却剂与燃料棒束多场耦合的方式研究了燃料棒束的流动传热特性和结构形变特性。结果表明：定位格架扰动冷却剂形成横向二次流并在下游棒束间形成绕流;多场耦合条件下二次流峰值速度和平均速度均小于单流场的;二次流与燃料棒的热应力使棒束发生形变,功率和流动分布的不均匀导致形变在轴向和径向的不均匀;相较于无格架情况,定位格架的存在使冷却剂的搅混流动更加明显,冷却剂对燃料棒冲击增大;在有、无定位格架两种情况下棒束形变均很小,可保持原本结构的稳定。     

基于丝网传感器的5×5棒束通道空泡分布测量研究

为研究压水反应堆燃料组件棒束通道内的两相分布规律,设计并制造了适用于棒束通道的丝网传感器模块,开展了5×5棒束通道内空气-水泡状流的空泡分布测量实验,分析了棒束通道内空泡份额的分布规律及气泡尺寸对空泡分布的影响。实验结果表明,发生横升力方向反转的小气泡在壁面附近聚集、大尺寸气泡则聚集在子通道中心;常温常压下发生横升力方向反转的临界气泡直径在4~6 mm之间,证明了横升力模型在棒束通道中的适用性。      

Prediction of Transient Boiling Transition in BWR Fuel Rod Bundles by Subchannel Analysis Program MENUETT

Transient CHF (critical heat flux) tests of a 4 X 4 rod bundle were analyzed by the subchannel analysis program MENUETT. MENUETT is based on a non-equilibrium, five equation, two-phase flow model and is available both for steady state and transient analyses. Turbulent mixing and void drift effects are taken into account to calculate cross flows in fuel rod bundles. The tendency of calculated subchannel mass fluxes and qualities agreed with experimental data. By using a critical quality correlation obtained from steady state CHF data, the position of the earliest boiling transition could be predicted regardless of non-uniform axial heat flux distributions. This transition occurrence time was predicted within a difference of 0.1～0.3 s from the experimental time. MENUETT applicability was confirmed for transient calculations predicting thermalhydraulic behavior in bundles.     

Experimental and numerical simulation of air-water two-phase flow in the rod bundle with grid spacer

1 Introduction Grid spacer is the key part of reactor fuel assem-bly. The presence of spacers in fuel assemblies affectsvarious thermal-hydraulic characteristics of the reactorcore. The grid spacer with fine performance can im-prove thermal-hydraulic performance of the core fuelassembly and enhance the critical heat flux withouttoo much augment of the pressure loss. As a result,the implementation of grid spacer with high thermalperformance provides more thermal margin, then in-creases s…     

Estimation of fuel and clad temperature of a research reactor during dry period of de-fuelling operation

Heat transfer and fluid flow studies related to spent fuel bundle of a research reactor in fuelling machine has been carried out. When the fuel is in reactor core, the heat generated in the fuel bundle is removed by heavy water under normal reactor operation. However, during the de-fuelling operation, the fuel bundle is exposed to air for some period called dry period. During this period, the decay heat from fuel bundle has to be removed by air flow. This flow of air is induced by natural convection only. In this period, the temperatures of fuel and clad rise. If clad temperature rises beyond a certain limit, structural failure may occur. This failure can result into release of fission products from fuel rod. Hence the temperature of clad has to be within specified limit under all conditions. The objective of this study is to estimate the clad temperature rise during the dry period.In the CFD simulation, the turbulent natural convection flow over fuel and radiation heat transfer are accounted. Standard k-? model for turbulence, Boussinesq approximation for computing the natural convection flow and IMMERSOL model for radiation are used.The steady state and transient CFD simulation of flow and heat is performed, using the CFD code PHOENICS. The steady state analysis provides the maximum temperature the clad will attain if fuel bundle is left exposed to air for sufficiently long time. For safe operation, the clad temperature should be limited to a specified value. From steady state CFD analysis, it is found that steady state clad temperature for various decay powers is higher than the limiting value. Hence transient analysis is also performed. In the transient analysis, the variation of clad temperature with time is predicted for various decay powers. Safe dry time, i.e. the time required for clad to reach the limiting value, is predicted for various decay powers. Determination of safe dry time helps in deciding the time available to the operator to drop the bundle in light water pool for storage. The analysis is found useful in optimizing the de-fuelling process.     

The Effect of Mixing-vane Arrangements in a Subchannel Turbulent Flow

Large eddy simulation (LES) of developed turbulent flows in a rod bundle was carried out for four spacer designs. The mixing-vanes attached at the spacer were inclined at 30° or 20° they were arranged to promote the swirling or convective flow. These arrangements are possible elements to compose an actual rod bundle. Our LES technique with a consistent higher-order immersed boundary method and a one-equation dynamic sub-grid scale model contributed to an efficient treatment of the complex wall configurations of rods and spacers. The computational results reasonably reproduced experimental results for the drag coefficient and the decay rate of swirling flow. The profiles of the axial velocities and the turbulence intensities indicated reasonable trend for the turbulent flow in the rod bundle. The effect of mixing-vane arrangement on the lateral flows was successfully clarified: the cross flow took the longer way on the rod surface than the swirling flow and then was more significantly influenced by momentum diffusion at the no-slip wall. Therefore, the largely inclined mixing-vanes promoted the cross flow only in the neighborhood of the spacer; the swirling flow inside a subchannel could reach farther downstream than the cross flow.     

数值反应堆堆芯通道级三维热工水力程序CorTAF开发及初步验证

堆芯是核动力系统的核心部件,其完整性是反应堆安全运行的重要前提。传统核反应堆堆芯热工水力分析方法无法满足未来先进核动力系统的高精度模拟需求。本文依托开源CFD平台OpenFOAM,针对压水堆堆芯棒束结构特点建立了冷却剂流动换热模型、燃料棒导热模型和耦合换热模型,开发了一套基于有限体积法的压水堆全堆芯通道级热工水力特性分析程序CorTAF。选取GE3×3、Weiss和PNL2×6燃料组件流动换热实验开展模型验证,计算结果与实验数据基本符合,表明该程序适用于棒束燃料组件内冷却剂流动换热特性预测。本工作对压水堆堆芯安全分析工具开发具有参考和借鉴意义。     

基于MPS算法的2×2棒束结构熔化行为数值模拟

以典型压水堆燃料组件2×2棒束结构为研究对象,建立了含定位格架和不含定位格架的棒束三维模型,基于半隐式运动粒子（MPS）算法对严重事故背景下棒束结构的熔化行为进行了数值模拟,分析了定位格架对棒束熔化过程中流道堵塞进程的影响。结果表明:MPS算法能够较好地模拟棒束结构熔化行为,定位格架会加快堆芯的熔化进程和冷却流道的堵塞速度,本文研究结果有利于严重事故下堆芯熔化模型的优化改进。      

带MVG和MSMG的5×5全长棒束单相流场和温度行为数值分析

定位格架作为燃料组件的关键部件之一,直接影响到燃料组件的热工性能。本文对带结构格架(MVG)和跨间搅混格架(MSMG)的5×5全长加热棒束单相流场和温度场采用计算流体力学(CFD)程序进行数值分析研究,获得该特征棒束组件出口二次流场以及温度场分布特性。研究表明,定位格架下游流场受定位格架和距离的影响,定位格架上游流场对下游二次流几乎无影响,定位格架导致流体强烈的横向二次流,增强了流体和加热棒之间的换热能力,使得棒束子通道截面流体温度更加均匀。与5×5全长棒束出口子通道温度的实验数据对比分析表明,获得的计算模型可以较好地分析该型棒束组件结构温度场行为。     

Single-phase convective heat transfer in rod bundles

The convective heat transfer for turbulent flow through rod bundles representative of nuclear fuel rods used in pressurized water reactors is examined. The rod bundles consist of a square array of parallel rods that are held on a constant pitch by support grids spaced axially along the rod bundle. Split-vane pair support grids, which create swirling flow in the rod bundle, as well as disc and standard support grids are investigated. Single-phase convective heat transfer coefficients are measured for flow downstream of support grids in a rod bundle. The rods are heated using direct resistance heating, and a bulk axial flow of air is used to cool the rods in the rod bundle. Air is used as the working fluid instead of water to reduce the power required to heat the rod bundle. Results indicate heat transfer enhancement for up to 10 hydraulic diameters downstream of the support grids. A general correlation is developed to predict the heat transfer development downstream of support grids. In addition, circumferential variations in heat transfer coefficients result in hot streaks that develop on the rods downstream of split-vane pair support grids.     

Cooling mechanism during transient reflooding of a reactor fuel bundle after loss of coolant

Cooling efficiency during transient reflooding under loss of normal coolant conditions has been examined with a 7 × 7 simulated fuel rod bundle and jet pump bypass. The bundle contains 49 electrically heated rods with 3600 mm heated length and a pseudo cosine axial power distribution. Water is injected into the lower plenum and the superheated bundle is reflooded from the bottom with some flow diverted to the simulated jet pump bypass. The results show that effective cooling can be maintained.     

5×5棒束通道中简单支撑格架对流动过冷沸腾传热的影响分析

采用壁面热分配模型(即RPI模型)对PSBT基准题中的5×5均匀加热全长棒束过冷沸腾传热进行了数值模拟研究。重点分析了加热段末端搅混格架（MVG）下游简单支撑格架（SSG）对棒束通道内流动过冷沸腾传热特性的影响。在水力特性方面,研究发现SSG的形阻压降约为MVG的53%,且对棒束通道内的横向流动具有显著抑制作用。为反映SSG对搅混过程的影响,采用子通道平均横流速度比沿轴向的发展过程对其进行了分析。分析发现,在SSG附近横流速度比迅速衰减,衰减后的横流速度比与光棒束时的大小相当。由于SSG对横流过程的破坏,改变了发热表面的传热特性,在其下游气相迅速包覆加热表面,蒸发热流密度较无SSG情况偏高5%,加热段末端空泡份额偏高0.006,壁面过热度偏高0.3 ℃。