摇摆条件下典型通道间湍流的流动传热特性

利用Fluent软件分析了摇摆条件对典型四棒束间的湍流流体流动和传热特性的影响机理。摇摆运动会对棒束间流体的流动传热特性产生一定影响,但不会对绝热通道与加热通道内流体流动相似性产生影响。而当摇摆幅度较大时,径向附加力会使通道横截面上的参数分布发生显著的变化,进而影响流体的流动与传热特性。在摇摆条件下,随着P/D（棒间距/棒直径）的逐渐减小,尤其是小于1.1时,典型棒束间流体的流动传热特性发生明显变化。     

均匀加热全长棒束过冷沸腾工况子通道参数场计算分析

采用壁面热分配模型对PSBT基准题中的5×5均匀加热全长棒束过冷沸腾传热进行了数值模拟研究,分析了均匀加热全长棒束通道中不同子通道和加热元件表面参数沿轴向的发展过程和径向的分布特性。研究发现,角通道和边通道是弱对流区域,其质量流速低于棒束平均值,但由于冷棒功率偏低,消除了流动不均衡性对传热效果的影响。在棒束径向方向,不同位置子通道间参数场存在差异,这是由于位于搅混格架横向导流对角方向的通道具有更有效的通道间对流效果,其传热效果更好。这种流动特性引起的参数差异在角通道中尤为显著。热棒表面过热度明显高于冷棒过热度,且位于非搅混格架横向导流方向的热棒具有更高的壁面过热度。     

7棒束紧密栅元流体流动传热数值研究

紧密栅元内的流体流动传热研究对高转化比反应堆燃料组件的优化有十分重要的意义。本文采用CFD方法对7棒束紧密栅元棒束通道内流体流动传热现象进行了数值模拟,并与7棒束紧密栅元内氟利昂流体传热的实验结果进行对比分析,详细分析了定位格架对棒束内流体传热流动的影响。结果表明：数值计算所得的非加热棒的壁面温度和实验吻合良好,定位格架的存在对其下游流体流动、棒束最高温度分布及交混系数有明显的影响,棒束某些位置因流动滞止导致温度大幅上升,在设计中应加以注意。     

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

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

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

事故条件下路基核反应堆以及受到海洋条件附加惯性力影响的浮动核反应堆一回路冷却剂会处于非稳定流动状态,进而改变冷却剂的流动和传热特性,影响反应堆的安全运行.本文应用锁相粒子图像测速(PIV)以及折射率匹配技术分别对脉动流条件下有无定位格架棒束通道内瞬时速度进行了测量.实验结果表明:对于不带定位格架的棒束通道,加速使得靠近...     

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

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

摇摆条件下子通道内湍流流体的数值分析

利用FLUENT软件分析了摇摆条件对典型四棒束间的湍流流体流动和传热特性的影响机理。摇摆运动会对棒束间流体的流动传热特性产生一定影响。RSM模型可以很好地描述摇摆条件下子通道内的参数分布。摇摆周期变化带来的径向附加力的变化不会对摩擦阻力系数、传热系数和Reynolds应力产生影响。在摇摆条件下,摩擦阻力系数、传热系数和Reynolds应力呈周期性变化,但最大摩擦阻力系数所在时刻并不固定,而最大传热系数却始终是在流速最大的时刻。     

六边形排列的7棒束通道内液态钠流动换热特性试验研究

以液态钠作为试验工质,对六边形排列的7棒束通道内液态钠流动换热特性进行了试验研究。试验流速为0~4 m·s^-1,热流密度为0~120 kW·m^-2,系统压力为1.5~200 kPa,对应的雷诺数和佩克莱数分别为4 000~60 000和0~340。深入分析了部分热工参数对7棒束通道内液态钠流动换热特性的影响,通过对7棒束通道内液态钠流动换热的试验数据的非线性拟合,得到适用于7棒束通道内液态钠流动换热的经验关系式。结果表明：拟合得到的摩擦系数关系式和努塞尔数关系式能准确地预测7棒束通道内的试验数据,其预测误差分别小于5%和6%。将获得的努塞尔数关系式与其他研究者的试验数据进行比较,与其他研究者985%的试验数据误差在30%以内,表明获得的关系式适用于7棒束通道内液态钠流动换热。     

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

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

高压条件下燃料组件内过冷沸腾过程传热不均匀性研究

基于流体动力学软件Fluent中的流体体积函数（VOF）两相流模型,通过编写用户自定义函数（UDF）程序添加控制方程源项,建立过冷沸腾模型,对压水堆带定位格架的5×5燃料组件棒束通道内的过冷沸腾现象进行数值模拟。根据模拟结果,从空泡份额、燃料棒周向传热方面对比分析各个子通道内传热特性。研究发现各子通道内空泡份额的分布不均匀性较大,同样加热条件下,边通道的沸腾程度高于角通道。此外,对棒束周向的传热特性进行了分析,燃料棒周向努塞尔数呈不均匀性分布,燃料棒0°、90°、180°、270°等方向附近的传热能力较强,其相应的横向速度较大,对应的沸腾程度较强。      

The flow and heat transfer characteristics of turbulent flow in typical rod bundles in ocean environment

The flow and heat transfer characteristic of turbulent flow in typical 4 and 7 rod bundles in ocean environment is investigated theoretically. In ocean environment, the periodic variation of secondary flow in 7 rod bundles is not obvious. Because of the velocity oscillation, there is a periodic heat accumulation on the tube wall. And the restriction of the channel wall on the rolling motion is considerable. In 7 rod bundles, because of the restriction of the channel wall, the effect of the additional force perpendicular to flowing direction is limited, and the turbulent flowing and heat transfer is mainly determined by the axial turbulent intensity and inlet velocity. However, in the 4 rod bundles, the restriction of the channel wall is small. The effect of the additional force perpendicular to flowing direction on the flowing and heat transfer is significant. And the additional force perpendicular to flowing direction can also affect the Reynolds stress.     

Numerical analysis of flowing characteristics of turbulent flow in rectangular channels in ocean environment

Because of the periodic effects of ocean waves, there are great discrepancies between the operational characteristics of nuclear power systems in ocean environment and that of land-based nuclear power systems. In some special operational status, like natural circulation, the additional forces due to ocean environment may impose so great disturbance on the coolant flow that theatres the safety operation of the systems. In the present paper, the turbulent flow in rectangular channels in ocean environments is investigated theoretically with CFD code FLUENT. The effects of several parameters on turbulent flow are analyzed. The effects of rolling motion includes two parts, the first part is the additional force parallel to flowing direction, which can affect on the pressure drop of the flow and change the flowing velocity, and the other part is the additional force perpendicular to flowing direction. In ocean environments, the flowing characteristics of turbulent flow are dominated by the additional force parallel to flowing direction. The effect of additional force perpendicular to flowing direction is very limited. In rolling and heaving motions, if the flowing velocity is the same, the flowing characteristics of turbulent flow are nearly the same, too. The bigger the Reynolds number is, the more serious the oscillation of turbulent kinetic energy and frictional resistance coefficient is, and the more the oscillation of turbulent flow is. The relationship between average frictional resistance coefficient and velocity oscillating amplitude is quadratic. And the oscillating amplitude of frictional resistance coefficient is in direct ratio with velocity oscillating amplitude.     

LES and URANS simulations of turbulent flow in 4 rod bundles in ocean environment

The flow and heat transfer of turbulent flow in typical 4 rod bundles in rolling motion is investigated with LES and URANS. The effect of rolling motion consists of two parts, the axial additional force which causes velocity oscillation and the radial additional force. The effect of rolling motion on the flowing similarity is considerable. The effect of radial additional force on the flow should not be neglected. In ocean environment, the effect of radial additional force on the flow should not be neglected. The average parameters are determined by the drive force and axial additional force, but the parameter profiles in the cross section are mainly determined by the radial additional force.     

Numerical investigation of supercritical water-cooled nuclear reactor in horizontal rod bundles

The commercial CFD code STAR-CD v4.02 is used as a numerical simulation tool for flows in the supercritical water-cooled nuclear reactor (SCWR). The basic heat transfer element in the reactor core can be considered as round rods and rod bundles. Reactors with vertical or horizontal flow in the core can be found. In vertically oriented core, symmetric characters of flow and heat transfer can be found and two-dimensional analyses are often performed. However, in horizontally oriented core the flow and heat transfer are fully three-dimensional due to the buoyancy effect. In this paper, horizontal rods and rod bundles at SCWR conditions are studied. Special STAR-CD subroutines were developed by the authors to correctly represent the dramatic change in physical properties of the supercritical water with temperature. In the rod bundle simulations, it is found that the geometry and orientation of the rod bundle have strong effects on the wall temperature distributions and heat transfers. In one orientation the square bundle has a higher wall temperature difference than other bundles. However, when the bundles are rotated by 90° the highest wall temperature difference is found in the hexagon bundle. Similar analysis could be useful in design and safety studies to obtain optimum fuel rod arrangement in a SCWR.     

Dynamic Model Representing Condensate Feedwater System of BWR Plants

In order to study the effect of burst temperature on the coolant flow channel restriction, burst tests of fuel bundles were performed. Each bundle consisted of 49 rods (7×7 rods), and bursts were conducted in flowing steam. Burst temperature was changed by changing the internal gas pressure in rods. After the burst, the ballooning behavior of each rod and the degree of coolant flow area restriction in the bundle were measured.

Maximum swelling of rod occurs when the burst temperature is around α and α+β phase boundary, and this phenomenon is almost the same as that in single rod burst tests. Maximum coolant flow area restriction is also observed in this condition.     

带交混翼矩形流道内超临界流动传热CFD研究

针对超临界水堆使用的棒束,取1根加热棒为中心的矩形区域作为研究对象,采用CFD中CFX10作为数值计算方法,研究了25MPa超临界压力下,有无交混翼时流体的流动和传热特点。计算结果表明：交混翼可增加流道内温度分布的均匀性,同时具有强化传热作用。     

The influence of ocean conditions on thermal-hydraulic characteristics of a passive residual heat removal system

There are many differences between the flow and heat transfer characteristics of nuclear reactors under ocean and land-based conditions for the effects of ocean waves. In this paper, thermal hydraulic characteristics of a passive residual heat removal system (PRHRS) for an integrated pressurized water reactor (IPWR) in ocean environment were investigated theoretically. A series of reasonable theoretical models for a PRHRS in an IPWR were established. These models mainly include the core, once-through steam generator, nitrogen pressurizer, main coolant pump, flow and heat transfer and ocean motion models. The flow and heat transfer models are suitable for the core with plate-type fuel element and the once-through steam generator with annular channel, respectively. A transient analysis code in FORTRAN 90 format has been developed to analyze the thermal–hydraulic characteristics of the PRHRS under ocean conditions. The code was implemented to analyze the effects of different ocean motions on the transient thermal-hydraulic characteristics of PRHRS. It is found that the oscillating amplitudes and periods of the system parameters are determined by those of the ocean motions. The effect of rolling motion is more obvious than that of pitching motion when the amplitudes and periods of rolling and pitching motions are the same. The obtained analysis results are significant to the improvement design of the PRHRS and the safety operation of the IPWR.     

CFD analyses in tight-lattice subchannels and seven-rod bundle geometries of a Super Fast Reactor

This paper presents CFD analyses in heat unsymmetric subchannels and heat symmetric seven-rod bundle geometries of a Super Fast Reactor (Super FR) fuel assembly using STAR-CD. The purpose of CFD analyses in heat unsymmetric subchannels is to evaluate the effect of the power differences on the heat transfer in subchannels of the Super Fast Reactor. For heat symmetric seven-rod bundles, the effects of the gap clearance between the fuel rod and the assembly wall and the displacement of the fuel rod on the circumferential temperature distributions and Maximum Cladding Surface Temperature (MCST) are analyzed. The results show that larger power difference between fuel rods gives larger circumferential temperature difference of the hottest fuel rods. Considering cross flow between edge and ordinary subchannels, 1 mm gap between the fuel rod and the assembly wall is better for small MCST although the circumferential temperature difference in edge subchannel is large. MCST increases exponentially with the displacement. The relative error of displacement should be less than 1% if the allowable increment of MCST due to displacement is less than 6 °C.     

超临界水四棒束传热数值分析

超临界水冷堆(SCWR)开发的关键是棒束内超临界水(SCW)的热工水力特性。本文针对超临界水四棒束流动传热实验进行CFD数值模拟,SSG湍流模型的计算结果与实验结果吻合良好。分析结果表明,流动方向对棒束截面内流量分布有显著影响。与下降流相比,尽管上升流时棒束间流动搅混较弱,但上升流时棒束截面流量及壁面周向温度分布更加均匀,加热棒壁面温度更低。可见,棒束横截面上的流量分布是影响加热棒壁面流动传热的主要因素。     

带格架四棒束超临界水流动传热数值分析

棒束内超临界水流动传热是超临界水堆堆芯热工水力研究的重要内容,但对其认识还十分有限。本文针对四棒束内超临界水的流动传热现象开展数值模拟,特别分析了定位格架对棒束通道内流动和传热的影响。结果表明,采用SSG湍流模型计算所得到的棒束壁面温度和实验结果吻合良好,定位格架的存在影响下游流体的速度分布,显著提高格架下游的传热特性,交混系数有大幅上升,使得加热棒周向壁面温度分布更加平均,最高温度出现位置发生改变。