The structure of turbulent flow in triangular array rod bundles

A wind tunnel study of fully developed uniform-density turbulent flow through triangular array rod bundles is described. Measurements were made for three tube spacings (

Full-size image (<1K)

) over a Reynolds number range of 12 000–84 000. The data include friction factors, local wall shear stresses, and the distributions of mean axial velocity, Reynolds stresses and eddy diffusivities. The secondary flow pattern is from the available evidence.     

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.     

Experimental investigation of turbulent flow through spacer grids in fuel rod bundles

This paper contains experimental data of pressure, velocity and turbulence intensity in a 24-rod fuel bundle with spacer grids. Detailed pressure measurements inside the spacer grid have been obtained by use of a sliding pressure-sensing rod. Laser Doppler Velocimetry technique was used to measure the local axial velocity and its fluctuating component upstream and downstream of the spacer grid in sub-channels with different blockage ratios. The measurements show a changing pattern in function of radial position in the cross-section of the fuel bundle. For sub-channels close to the box wall, the turbulence intensity suddenly increases just downstream of the spacer and then gradually decays. In inner sub-channels, however, the turbulence intensity downstream of the spacer decreases below its upstream value and then gradually increases until it reaches the maximum value at approximately two spacer heights. The present study reveals that spacer effects, such as local pressure distribution and turbulence intensity enhancement, not only depend exclusively on the local geometry details, but also on the location in the cross-section of the rod bundle.     

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.     

Two-phase turbulent mixing and buoyancy drift in rod bundles

This paper describes the development of generalized relationships for single- and two-phase intersubchannel turbulent mixing in vertical and horizontal flows, and lateral buoyancy drift in horizontal flows.The relationships for turbulent mixing, together with a recommended one for void drift, have been implemented in a subchannel thermalhydraulics code, and assessed using a range of data on enthalpy migration in vertical steam–water flows under BWR and PWR diabatic conditions. The intent of this assessment was to optimize these relationships to give the best agreement with the enthalpy migration data for vertical flows. The optimized turbulent mixing relationships were then used as a basis to benchmark a proposed buoyancy drift model to give the best predictions of void and enthalpy migration data in horizontal flows typical of PHWR CANDU¹ reactor operation under normal and off-normal conditions.Overall, the optimized turbulent mixing and buoyancy drift relationships have been found to predict the available data quite well, and generally better and more consistently than currently used models. This is expected to result in more accurate calculations of subchannel distributions of phasic flows, and hence, in improved predictions of critical heat flux (CHF).     

Measurements of turbulent velocity and temperature in axial flow through a heated rod bundle

Fully developed turbulent air flow in a heated 37-rod bundle with a pitch to diameter ratio of 1.12 has been investigated. Measurements were performed with a hot-wire probe with x-wires and a temperature wire. Besides the distributions of the mean velocity, mean fluid temperature, the wall shear stress and wall temperature, the turbulent quantities such as the turbulent kinetic energy, the Reynolds-stresses and the turbulent heat fluxes were measured and compared with data from isothermal flow and heated flow in pipes.     

A scale analysis of the turbulent mixing rate for various Prandtl number flow fields in rod bundles

The turbulent mixing rate is a very important variable in the thermal–hydraulic design of nuclear reactors. In this study, the turbulent mixing rate for the flow through rod bundles is estimated with the scale analysis on the flow pulsation generated by periodic vortices that is pointed out as a main cause of the mixing in rod bundles. Based upon the assumption that turbulent mixing is composed of molecular motion, isotropic turbulent motion (turbulent motion without the flow pulsation), and flow pulsation, the scale relation is derived as a function of P/D, Re, and Pr. The derived scale relation is compared with the published experimental results and shows good agreement. Since the scale relation is applicable to various Prandtl number fluid flows, it is expected to be useful for the thermal–hydraulic analysis of liquid metal coolant reactors as well as moderate Prandtl number coolant reactors.     

Annular flow experiments in rod bundles with spacers

Measurements of pressure drop, film flow rate and film thickness have been made on a special test section modelling the four subchannels between ten fuel rods in a LWR geometry. Tests were carried out with air/water at 2.5–3-bar (a). Measurements were made with and without spacer grids inserted in the channel. In the former case, the measurements were made upstream and downstream of the spacers. In the case without spacers, the results show similar trends to data from tubular geometries with a similar hydraulic diameter. The two spacers tested both showed similar pressure losses. There was an increase in film flow rate and film thickness downstream of the ULTRAFLOW spacer. For a conventional egg-crate spacer there was a small change, with in some cases a decrease in film flow rate.     

Critical heat flux prediction in vertical bottom-closed rod bundles

This paper presents the analysis of experimental data and calculational relationships for heat the transfer crisis in LWR rod bundle with closed bottom. A new relationship for critical heat flux prediction in the rod bundle with closed bottom based on the improved drift model is described. The comparison of critical heat flux values given by different correlations (including Groeneveld's algorithm used in RELAP5/MOD3.1 Code) and those obtained from the tests in the wide range of regime and geometric parameters is presented.     

Analytical prediction of friction factors and Nusselt numbers of turbulent forced convection in rod bundles with smooth and rough surfaces

A simple analytical method was developed for the prediction of the friction factor, f, of fully developed turbulent flow and the Nusselt number, Nu, of fully developed turbulent forced convection in rod bundles arranged in square or hexagonal arrays. The friction factor equation for smooth rod bundles was presented in a form similar to the friction factor equation for turbulent flow in a circular pipe. An explicit equation for the Nusselt number of turbulent forced convection in rod bundles with smooth surface was developed. In addition, we extended the analysis to rod bundles with rough surface and provided a method for the prediction of the friction factor and the Nusselt number. The method was based on the law of the wall for velocity and the law of the wall for the temperature, which were integrated over the entire flow area to yield algebraic equations for the prediction of f and Nu. The present method is applicable to infinite rod bundles in square and hexagonal arrays with low pitch to rod diameter ratio, P/D<1.2.     

An analytical method of evaluating the effect of rod displacement on circumferential temperature and heat flux distributions in turbulent flow through unbaffled rod bundles

An analytical method of evaluating the circumferential variations of temperature and heat flux fields inside and around a displaced fuel rod in triangular rod bundles in turbulent flow is presented with illustrative examples. The analysis consists mainly of the derivation of the simultaneous solutions of a set of heat conduction equations for fuel, cladding and coolant under the assumption of fully developed flow and heat transfer conditions. The local coolant velocity distribution, which is necessary for deriving the temperature field in coolant, is determined by solving the Navier-Stokes equation and the turbulent mixing of coolant is taken into consideration. The results show how the circumferential variations in the temperature and heat flux fields on the outer surface of the cladding increase the lower the ratio and the larger the fuel rod displacement due to thermal conduction and peripheral coolant flow velocity distribution.     

Analytical prediction of turbulent friction factor for a rod bundle

An analytical calculation has been performed to predict the turbulent friction factor in a rod bundle. For each subchannel constituting a rod bundle, the geometry parameters are analytically derived by integrating the law of the wall over each subchannel with the consideration of a local shear stress distribution. The correlation equations for a local shear stress distribution are supplied from a numerical simulation for each subchannel. The explicit effect of a subchannel shape on the geometry parameter and the friction factor is reported. The friction factor of a corner subchannel converges to a constant value, while the friction factor of a central subchannel steadily increases with a rod distance ratio. The analysis for a rod bundle shows that the friction factor of a rod bundle is largely affected by the characteristics of each subchannel constituting a rod bundle. The present analytic calculations well predict the experimental results from the literature with rod bundles in circular, hexagonal, and square channels.     

CFD application to prediction of void distribution in two-phase bubbly flows in rod bundles

This paper is concerned with the prediction of the void fraction distribution in two-phase bubbly flows in fuel rod bundles. Special attention has been devoted to the phenomena which govern the void fraction distribution in the lateral direction of a channel. A two-fluid model of two-phase flow has been formulated and implemented into a commercial computational fluid dynamics (CFD) code. The model has been used for the prediction of the void distribution in three different channels: a circular channel (inside diameter (ID), 34.5 mm) with a single heated rod of 13.9 mm outside diameter (OD), and circular channels (ID, 71 mm) with six heated rods (13.8 and 13.9 mm OD each). The predicted axial and lateral avoid fraction distributions in subcooled and bulk boiling regions have been area averaged in three lateral zones and compared with experimental data: in all cases, satisfactory agreement between the predictions and measurements has been obtained.     

Critical heat flux prediction in rod bundles under upward low mass flux densities

The analysis of experimental data and results of calculations for heat transfer crisis in heated channels under low upward coolant mass flux densities is presented. This analysis allows the determination of the basic features of the boiling crisis phenomenon. It is shown that the methods currently used for critical heat flux (CHF) prediction have insufficient accuracy in the given range of parameters. A new relationship for the CHF calculation is presented. It should be used for the water–water energy reactor (WWER) and uran–graphite channel reactor—Chernobyl-type (RBMK) rod bundles, and is verified by the test data. The comparison of results obtained by a new CHF correlation and the relationship used in RELAP5/MOD3.1 Code is presented. It is shown that the latter overpredicts the CHF values at atmospheric pressure and for x_cr>0.4 and does not provide conservative estimations for the RBMK fuel bundles.     

The structure of turbulent flow through a wall subchannel of a rod bundle

An experimental investigation was performed to establish reliable information on the transport properties of turbulent flow through subchannels of rod bundles. Detailed data were measured of the distributions of the time-mean velocity, the turbulence intensities in all directions and hence, the kinetic energy of turbulence, of the shear stresses in the directions normal and parallel to the walls and of the wall shear stresses for a wall subchannel of a rod bundle of four parallel rods. The pitch to diameter ratio of the rods equal to the wall to diameter ratio was 1.07, the Reynolds number of this investigation was Re = 8.7 × 10⁴.On the basis of the data measured the eddy viscosities in the directions normal and parallel to the walls were calculated. Thus, detailed data of the eddy viscosities in direction parallel to the walls in rod bundles were obtained for the first time. The experimental results were compared with predictions by the VELASCO code. There are considerable differences between calculated and measured data of the time-mean velocity and the wall shear stresses. Attempts to adjust the VELASCO code against the measurements were not successful. The reasons of the discrepancies are discussed.     

Developed single phase turbulent flow through a square-pitch rod cluster

An experimental study of developed single phase turbulent flow through a square pitched array of rod bundles is described. Measurements were made at two spacings (p/d = 1.194, 1.107) of the mean velocity distribution and wall shear stress variation, together with the six terms of the symmetrical Reynolds stress tensor. The departure of the turbulent flow structure from axisymmetric pipe flow, particularly in the rod gap region, was found to depend strongly on the (p/d) ratio.     

A theoretically based critical heat flux prediction for rod bundles at PWR conditions

A theoretically based procedure developed for round tubes has been applied to the prediction of DNB heat fluxes in rod bundles at PWR conditions. State-of-the-art subchannel analysis procedures were used to determine local flows and enthalpies. Very good comparison between DNB predictions and experimental observations are found for rod bundles which both uniform and non-uniform axial heat fluxes.     

Turbulent momentum transport in rod bundles

An experimental investigation was performed with air to obtain detailed information on the velocity and turbulence distribution for parallel turbulent flows through subchannels of rod bundles. Experimental results were obtained for wall and corner subchannels of rod bundles of four parallel rods. The pitch-to-diameter ratios were varied between 1.07 and 1.4. The Reynolds numbers ranged from 6 × 10⁴ to 2 × 10⁵, depending on the rod bundle arrangement.On the basis of the data measured, the eddy viscosities in the directions normal and parallel to the wall were calculated. The experimental results of the velocity and wall shear stress distributions are compared with the predictions by the VELASCO code. There are considerable differences between computed and experimental results especially for low pitch-to-diameter ratios. The reasons for the discrepancies are discussed together with the results of attempts to adjust the VELASCO code against the experimental data.     

Large eddy simulation of turbulent flow surrounding two simulated CANDU fuel bundles

Three-dimensional turbulent flow in one and a half simulated 43-element CANDU^®1 fuel bundle at the inlet of a fuel channel is solved using large eddy simulation. Wake generated after the endplates and the flow development in the inlet bundle are investigated based on the simulation results. Spatial distribution and frequency spectra of fluid force components are also examined. The simulation results are compared to experimental data available in the literature as well as the measurements conducted by the authors. The current investigation provides a basic understanding of the fluid excitation in a simulated 43-element fuel bundle. The results may be used in a flow-induced vibration analysis for fuel bundles.