脉动流下棒束通道内流场与湍流特性的PIV实验研究

事故条件下路基核反应堆以及受到海洋条件附加惯性力影响的浮动核反应堆一回路冷却剂会处于非稳定流动状态,进而改变冷却剂的流动和传热特性,影响反应堆的安全运行。本文应用锁相粒子图像测速（PIV）以及折射率匹配技术分别对脉动流条件下有无定位格架棒束通道内瞬时速度进行了测量。实验结果表明：对于不带定位格架的棒束通道,加速使得靠近通道壁面附近流体速度变大,而靠近中心区域流体速度变小。此外湍流强度分量随流体加速而逐渐变小,随流体减速而逐渐增加。对于流向压力梯度驱动的周期性脉动流,横向脉动速度均方根分量u′滞后于流向脉动速度均方根分量v′,且二者都滞后于流速的变化;对于带定位格架的棒束通道,带有搅浑翼的定位格架强烈的交混作用极大地减小了流体加速度对棒束通道内速度分布和湍流强度带来的影响。实验结果有助于更加清晰地揭示脉动流在棒束通道中的作用机理。     

脉动流下定位格架下游时均流场分布特性研究

研究流量波动下棒束通道内定位格架下游瞬时流场演变特性对于揭示海洋条件下燃料组件内流动换热机理具有重要意义。本文应用粒子图像测速（PIV）技术获得了脉动流下棒束通道内定位格架下游时空演变流场结构,分析了脉动参数（脉动周期和脉动振幅）对定位格架下游速度分布和湍流特性的影响。结果表明,脉动流下定位格架下游时均速度与定常流动下时均速度差异较小,且基本不随脉动振幅和脉动周期变化而变化;脉动流下的定位格架下游横向速度和轴向速度均方根与定常流动下的速度均方根存在明显差异,且随脉动参数变化呈现出不同的变化趋势。本文研究结果有助于揭示燃料组件在非稳态条件下瞬态特性,并为燃料组件的设计和优化奠定基础。      

矩形通道内加减速条件下流态转捩特性研究

针对恒压差变化率下矩形通道内的流态转捩特性进行理论和实验研究。通过理论研究,得到恒压差变化率下的层流速度分布。以层流速度分布为基础,应用能量梯度法,对加速和减速情况下的转捩雷诺数(Re)进行计算,并分析影响转捩Re的因素。分析结果表明:加速过程中转捩Re低于稳态情况下的转捩Re和减速过程中的转捩Re高于稳态情况下的转捩Re;初始稳态压差越接近稳态转捩时的压差,转捩Re偏离稳态转捩Re就越小;压差变化率对转捩Re几乎没有影响。通过加减速实验研究得到的转捩Re变化规律和理论结果基本一致,但在压差变化率影响方面存在差异。     

不同棒束结构稠密栅元通道内的湍流CFD研究

采用URANS（UnsteadyReynoldsAveragedNavierStokes）方法对不同棒束结构稠密栅元通道（P/D=1.001～1.2）内的湍流流动进行CFD模拟。研究分析了不同Re（Re=5000~215000）的湍流流动的主流速度、壁面剪应力、湍动能等参数。研究表明：在较稠密的棒束（P/D<1.1）通道内,P/D的变化对子通道内主流速度和剪应力分布均有较大影响。本文的模拟结果也验证了在达到临界P/D前（即使δ/D<0.011）,交混因子Y和δ/D成反比关系。对于固定的棒束结构（P/D=1.062）,当Re达到一定值（Re=9600）时,子通道内主流速度和剪应力分布对Re的变化不敏感。     

棒束内湍流流动特性的PIV与CFD研究

在压水堆燃料组件的定位格架下游,局部扰动沿流动方向逐渐衰减,流场最终趋于稳定。光滑棒束区冷却剂的湍流流动和交混特性是影响反应堆经济性和安全性的重要因素,有必要进行深入研究。本文采用粒子图像测速（PIV）与数值模拟(CFD)相结合的方法,对3×3小规模棒束内水的流动特性进行研究,得到了一阶平均流速以及二阶湍流统计信息。结果表明,中心子通道的速度明显高于棒间隙区,但轴向均方根速度呈现出相反的变化趋势。在相邻子通道横向速度梯度的作用下,棒束内出现了大尺度的流量脉动现象,且脉动波长随雷诺数的增加而增大。此外,实验得到的湍流交混系数较压水堆采用的Castellana公式预测值偏高10%左右,这一偏差随雷诺数的增加有减小的趋势。     

5×5棒束通道内流动转捩特性研究

棒束通道的特殊结构导致其内部流动转捩情况较为复杂,探究其内部流动转捩规律具有重要意义。本文针对棒束通道内的流动转捩特性开展实验与CFD模拟研究,通过实验获得了棒束通道内沿程阻力系数的变化规律;采用不同湍流模型进行了数值模拟。结果表明,SST k-ω模型能较好地反映实验结果。进一步对比了不同雷诺数工况下通道内不同位置的沿程阻力系数与湍流强度,发现对于不同子通道,中心子通道湍流强度与沿程阻力系数高于边角子通道;对于同一子通道,子通道中心处湍流强度与壁面切应力高于子通道边缘处。这一结果说明,受壁面影响,棒束内湍流强度、壁面切应力、阻力特性具有不均匀性,这些空间上的不均匀性相互作用会引起总体上棒束转捩点不明显。     

含绕丝棒束低流速时阻力特性的实验研究

钠冷快堆含绕丝棒束组件低流速时的摩擦阻力特性,是设计钠冷快堆非能动事故余热排出系统的重要参数。使用水作为测量介质,采用改进的液柱测压方法,准确测量了一个有机玻璃材质37棒绕丝棒束组件雷诺数Re=200～1 100时的摩擦阻力特性。结果表明,被测量组件从层流向湍流转捩的临界Re为370。组件的入口段长度随着流速的增加而增大,随着流动开始向湍流转变而减小。现有的绕丝棒束摩擦阻力系数经验公式中,UCTD公式与实验值最接近：在层流区,UCTD公式高估了约7%,在层流向湍流转换区,UCTD公式计算值与实验值基本一致。同时,UCTD和CTS公式还高估了组件从层流向湍流转捩的临界Re。     

高宽比对矩形窄通道内单相水流动特性的影响机理

通过理论和实验研究了矩形窄通道内单相水的流动特性。结果表明：矩形通道截面高宽比对其流动阻力有重要影响,层流区和紊流光滑管区摩擦阻力均随高宽比的减小而增大。可视化结果表明：矩形窄通道临界Re为2 500~2 600,高宽比对窄通道内流态转捩无明显影响。层流区高宽比通过改变截面湿周及速度分布两方面影响摩阻特性;紊流区高宽比通过改变黏性底层紊流时均速度分布及截面湿周两方面影响摩阻特性。     

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

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

加减速流条件下棒束通道阻力特性研究

压水反应堆的启动、停堆、事故等多种工况会使反应堆堆芯内部冷却剂产生加减速流,影响反应堆的热工水力特性。在雷诺数(Re)700相似文献   

绕丝棒束组件低流速时摩擦阻力实验研究

钠冷快堆在事故停堆余热排放期间,堆芯组件内钠流为自然循环流动,流速很低,因此准确确定绕丝棒束组件低流速时的摩擦阻力系数对钠冷快堆非能动余热排出系统的设计具有重要意义。本文以水为流动介质,准确测量了37棒和19棒绕丝棒束组件在低流速（Re<1 000）时的摩擦阻力系数。实验结果表明,随着流速的降低,绕丝棒束组件的摩擦阻力系数迅速升高,流动从层流向过渡流转变时,摩擦阻力系数有明显跃升。将实验测量值与绕丝棒束摩擦阻力系数经验公式的计算结果进行比较,发现在低流速时,经验公式计算结果较实验测量值明显偏小,同时经验公式计算的绕丝棒束层流向过渡流转变的临界Re较实验值偏大。     

棒束通道内混合对流传热实验研究

以自然循环下堆芯内可能会发生的低流量传热为研究背景,对5×5棒束通道内的混合对流传热现象进行了实验研究。实验压力为6 MPa, 质量流量为25~150 kg/(m2·s),热流密度为25~300 kW/m2,实验雷诺数Re为1000~30000,浮升力参数Bo*为2×10-7~3×10-3。实验发现,随着Bo*的增大,棒束通道内传热产生先弱化后强化的趋势。浮升力对棒束通道内传热造成影响的起始点为Bo*=3.5×10-6,当Re >15000时,浮升力依然可对传热造成弱化现象。基于实验数据,提出了适用于棒束通道的混合对流经验关系式。      

快堆带绕丝棒束组件低雷诺数下的水力特性分析

低雷诺数(Re)流动存在于正常运行或事故停堆工况的各类组件中,对于快堆的安全运行具有重要意义。利用CFX程序对低Re下的中国实验快堆不同类型的带绕丝棒束组件的水力特性进行了分析。结果表明,通过利用1个螺距的带绕丝棒束组件计算得到的低Re下的水力特性与实验结果以及Engel关系式符合较好。通过利用4个螺距的带绕丝棒束组件计算结果表明,绕丝产生的横向流动使组件6个壁面上压力分布有所不同,但在流动充分发展时,每个面轴线方向的压降按螺距均匀分布,从而进行带绕丝棒束组件水力特性测量时,需在组件同一面上按照整数倍螺距来布置测点,才能避免由于横向流动对测量带来的影响。     

棒束通道内过冷沸腾起始点的实验研究

本文以去离子水为实验介质,在进口温度80～100 ℃、质量流速0～100 kg/(m²•s)、热流密度0～80 kW/m²的条件下对棒束通道内的过冷沸腾起始点（ONB）进行了实验研究。分析了部分热工参数和棒束特殊的几何结构对ONB的影响,通过引入雷诺数,对棒束通道内ONB的数据进行非线性回归分析,得到适用于棒束通道ONB的经验关系式。结果表明：新拟合得到的关系式能较准确地预测棒束通道内ONB的热流密度,其预测值的相对误差为14.75%。     

Numerical simulation of turbulent flow in a 37-rod bundle

The unsteady Reynolds averaged Navier–Stokes equations, combined with a Reynolds stress model, were solved numerically to determine fully developed isothermal turbulent flow in a 60° sector of a 37-rod bundle. It was found that this flow contained large-scale coherent structures, which affected strongly the local velocity fluctuations, especially near the gaps between rods or between rods and the surrounding wall. The time-averaged mean velocity and Reynolds stresses were in good agreement with experimental results in a similar channel. Coherent velocity fluctuations at different locations throughout the entire rod bundle were strongly correlated with each other.     

Mixed convective flow penetration in vertical channels

The primary purpose of the study is to investigate the factors relevant to the decay heat removal system in pool-type liquid metal reactors which are designed to remove decay heat in a passive way utilizing natural circulation. The reactor geometry is simulated by a vertical rod bundle channel connected to an upper plenum. Penetration of cold fluid from the upper plenum into the rod bundle channel is investigated experimentally and analytically with water as a working fluid. Three correlations to predict the onset of penetration, the penetration depth, and the ratio of penetrating to forced flowrate were developed. The correlations were found to agree well with experimental results for the range of Reynolds number in which experimental data were obtained.     

Flow Transition Characteristics in 5×5 Rod Bundle Channel

The special geometric structure of the rod bundle channel can induce complicated flow transition of the coolant, and investigation on the flow transition rules is sufficiently important. In the current study, experimental and numerical study on the flow transition characteristics in the 5×5 rod bundle channel was carried out. Experiments were performed to obtain the variation characteristics of the resistance coefficient and CFD simulation was performed using different turbulence models in ANSYS Fluent. The results show that the simulation with SST k-ω turbulence model agrees well with the experimental data. The simulated turbulence intensity and resistance coefficient at different measurement locations and in different flow conditions were compared. For different subchannels, the turbulence intensity and the resistance coefficient are higher in the center subchannel than those in the edge subchannel. For the same subchannel, the turbulence intensity and the shear stress in the subchannel center are higher than those in the subchannel edge. This result indicates that the turbulence intensity, shear stress and resistance coefficient in the rod bundle are not uniform due to the influence of the wall surface. This non-uniform spatial interaction makes the transition point obscure.     

Analysis of flow distribution a PWR fuel rod bundle model containng A 90% blockage

An experimental investigation, covering a Reynolds number range from 2 × 10³ to 3.5 × 10⁴, was conducted to study the velocity and turbulence intensity distributions due to the presence of a blockage in an unheated 7 × 7 rod bundle. The blockage configuration, consisting of a 4 × 4 rod array, created a maximum flow area reduction of 90% in the central nine subchannels. The blockage sleeve length was 38.3 × rod diameter and the 90% blockage zone length extended for 16.4 × rod diameter. The results showed that upstream of the blockage, the flow was not influenced by the blockage until it reached the location where the inlet taper section of the swelling started. At the downstream end, the flow disturbance was extensive and persisted over a distance of about 83 rod diameters. Compared to the downstream velocity profiles, the turbulence intensity measurements however showed a faster recovery from the blockage influence. At the higher Reynolds number, velocity profiles calculated using the COBRA subchannel computer code compared consistently with the experimental data. The general flow behaviour of the various subchannels was reasonably well predicted. However, at low Reynolds number, due mainly to the frictional form loss calculation scheme in COBRA and uncertainty in the flow transition, the flow diversion due to the blockage to the surrounding unblocked subchannels was overestimated. The influence of the degree of recovery from the rod swelling on the flow was also studied using COBRA.