Thermal-hydraulic performance analysis of a subchannel with square and triangle fuel rod arrangements using the entropy generation approach

The present paper discusses entropy generation in fully developed turbulent flows through a subchannel,arranged in square and triangle arrays. Entropy generation is due to contribution of both heat transfer and pressure drop. Our main objective is to study the effect of key parameters such as spacer grid, fuel rod power distribution,Reynolds number Re, dimensionless heat power ω, lengthto-fuel-diameter ratio λ, and pitch-to-diameter ratio ξ on subchannel entropy generation. The analysis explicitly shows the contribution of heat transfer and pressure drop to the total entropy generation. An analytical formulation is introduced to total entropy generation for situations with uniform and sinusoidal rod power distribution. It is concluded that power distribution affects entropy generation.A smoother power profile leads to less entropy generation.The entropy generation of square rod array bundles is more efficient than that of triangular rod arrays, and spacer grids generate more entropy.     

A study of the effect of rod-bowing on critical heat flux

An experimental study was carried out to determine the effect of rod-bowing on critical heat flux, using an electrically-heated rod cluster. In this experiment, rod-bow was set to occur in the severest subchannel and axially at the middle between the last two spacers, with uniform axial heat flux. The minimum gap between the outer and inner rods was reduced variously to 1.6 mm, 1.0 mm and zero from the nominal value of 2.1 mm. Other experimental conditions were as follows: pressure 7 MPa; mass velocity 640–2600 kg/m² sec; inlet subcooling 40–560 kJ/kg.Experimental results show only a slight rod-bowing effect, if any, compared with normal spacing, as confirmed by analysis of three-dimensional heat conduction around the rod-bowing area and by the local steam quality deviations calculated by subchannel analyses.     

Local heat transfer and hot-spot factors in wire-wrap tube bundle

Heat transfer coefficients and hot-spot factors have been determined from measured local temperatures and calculated local mass flux in seven adjacent tubes and associated subchannels of a 61 wire-wrap tube bundle characteristic of the blanket of a GCFR (Gas Cooled Fast Reactor). The bundle consisted of 2.11 cm OD stainless steel tubes on a triangular array with a pitch/diameter ratio of P/D = 1.05. The helical wire of 0.1067 cm in diameter was coiled on the tube with a respective initial orientation of 0–120–240°C and 30.48 cm helical pitch. The experiment used water at atmospheric pressure and temperature as coolant. The resulting dimensionless correlation for heat transfer is applicable to gases and all non-metal fluids in one phase flow when the fluid properties at subchannel bulk temperature are used. This correlation is based on local subchannel mass flux and is applicable to all wire-wrap configurations. Local subchannel mass fluxes were determined with a computer program COBRA IV and used to correlate the average Nusselt number for each subchannel in terms of local Reynolds number and fluid Prandtl number. The differences of up to 19% between that correlation and the one presented in earlier work are discussed in the text. The hot-spot factors on the convective heat transfer coefficient for tubes and subchannels are given as a function of Reynolds number based on a bundle average mass flux and a local subchannel hydraulic diameter. These factors are specific to the bundle configuration and are also dependent on the wire-wrap configuration.     

Power up-grading study for the first Egyptian research reactor

In the present work, power up-grading study is performed, for the first Egyptian Research Reactor (ET-RR-1), using the present fuel basket with 4×4 fuel rods, (17.5 mm pitch), and a proposed fuel basket with 5×5 fuel rods, (14.0 mm pitch), without violating the thermal hydraulic safety criteria. These safety criteria are; fuel centerline temperature (fuel melting), clad surface temperature (surface boiling), outlet coolant temperature, and maximum heat flux (critical heat flux ratio). Different thermal reactor powers (2–10 MW) and different core coolant flow rates (450, 900, 1350 m³ h⁻¹) are considered. The thermal hydraulic analysis was performed using the subchannel code COBRA-IIIC for the estimation of temperatures, coolant velocities and critical heat flux. The neutronic calculations were performed using WIMS-D4 code with 5 — group neutron cross section library. These cross sections were adapted to use in the two-dimensional (2-D) diffusion code DIXY for core calculations. The study concluded that ET-RR-1 power can be upgraded safely up to 4 MW with the present 4×4-fuel basket and with the proposed 5×5-fuel basket up to 5 MW with the present coolant flow rate (900 m³ h⁻¹). With the two fuel arrays, the reactor power can be upgraded to 6 MW with coolant flow rate of 1350 m³ h⁻¹ without violating the safety criterion. It is also concluded that, loading the ET-RR-1 core with the proposed fuel basket (5×5) increases the excess reactivity of the reactor core than the present 4×4 fuel matrix with equal U-235 mass load and gave better fuel economy of fuel utilization.     

Natural circulation capability of Pb-Bi cooled fast reactor: PEACER

First-principle calculations were performed to analyze the natural circulation heat removal from the core of a liquid metal reactor (LMR). The lead-bismuth (Pb-Bi) was chosen as the primary coolant for the LMR system. From the single channel analysis the temperature and the pressure drop are calculated along the fuel assembly. The total pressure drop of the core is less than 100kPa due to the large pitch-to-diameter ratio and the small height of the fuel pin. The natural circulation potential is a key characteristics of the LMR design. The steady-state momentum and energy equations are solved along the primary coolant path. The calculations are divided into two parts: an analytical model and a one-dimensional lumped parameter flow loop model. Results of the analytical model indicate that the elevation difference of 4.5m between thermal centers of the core and the steam generators could remove as much as 10% of the nominal operating reactor power. The flow loop model yielded the total pressure drop and the natural circulation heat removal capacity.     

Matrix calculation of radiant heat transfer in LWR fuel bundles under accident conditions

A computer code was developed for calculating the radiant heat transfer in a LWR fuel bundle under accident conditions. The calculation method is a modular one: a fuel bundle or its part is divided into unit cells, each of which is composed of a coolant subchannel surrounded by several segments of solid or imaginary faces. The view factor matrix in each cell is expanded over the whole bundle using the concept of ‘boundary face’ between cells, and the resultant heat transfer equations are simultaneously solved for solid wall temperatures. The geometrical flexibility of this method is suitable for treating various simulation experiments for accidents. The method is also effective for repeated calculations of the radiant heat transfer reflecting state or material property changes when analyzing fuel rod behaviour under accident conditions.     

Analysis of radiant heat transfer in a BWR fuel assembly

A computer code ‘CIDER’ was developed which analyzes radiant heat transfer in a BWR fuel rod bundle under loss of coolant conditions. In the code, (1) a channel box and fuel rods are considered to be gray bodies, (2) reflection and absorption of radiation beams in the atmosphere is neglected, (3) a fuel rod is approximated by a regular polygonal rod, and (4) radiant heat flux is calculated considering circumferential temperature distribution on each fuel rod surface, which is determined from radial and circumferential heat conduction calculations in a fuel rod. It was found that the conventional model with uniform cladding temperature overestimated heat flux about 30% in a typical situation, or correspondingly underestimated the temperature rises.     

Development of a bundle correction method and its application to predicting CHF in rod bundles

A bundle correction method, based on the conservation laws of mass, energy, and momentum in an open subchannel, is proposed for the prediction of the critical heat flux (CHF) in rod bundles from round tube CHF correlations without detailed subchannel analysis. It takes into account the effects of the enthalpy and mass velocity distributions at subchannel level using the first derivatives of CHF with respect to the independent parameters. Three different CHF correlations for tubes (Groeneveld's CHF table, Katto correlation, and Biasi correlation) have been examined with uniformly heated bundle CHF data collected from various sources. A limited number of CHF data from a non-uniformly heated rod bundle are also evaluated with the aid of Tong's F-factor. The proposed method shows satisfactory CHF predictions for rod bundles both uniform and non-uniform power distributions.     

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.     

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.     

Prediction of the single-phase turbulent mixing rate between two parallel subchannels using a subchannel geometry factor

This paper presents a simple method for predicting the single-phase turbulent mixing rate between adjacent subchannels in nuclear fuel bundles. In this method, the mixing rate is computed as the sum of the two components of turbulent diffusion and convective transfer. Of these, the turbulent diffusion component is calculated using a newly defined subchannel geometry factor F^* and the mean turbulent diffusivity for each subchannel which is computed from Elder's equation. The convective transfer component is evaluated from a mixing Stanton number correlation obtained empirically in this study. In order to confirm the validity of the proposed method, experimental data on turbulent mixing rate were obtained using a tracer technique under adiabatic conditions with three test channels, each consisting of two subchannels. The range of Reynolds number covered was 5000–66 000. From comparisons of the predicted turbulent mixing rates with the experimental data of other investigators as well as the authors, it has been confirmed that the proposed method can predict the data in a range of gap clearance to rod diameter ratio of 0.02–0.4 within about ±25% for square array bundles and about ±35% for triangular array bundles.     

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.     

A compilation of rod array critical heat flux data sources and information

The geometric characteristics of rod array test sections employed in critical heat flux (CHF) tests with water coolant, and the ranges of the operating parameters for the tests, are presented for 126 test sections. The corresponding 4277 CHF data points have been stored on a magnetic tape for ease of reference and analysis. A versatile computer program associated with the data library has been used to determine the distributions of the data with respect to geometric and operating parameters. The dependence of CHF on operating parameters and the importance of subchannel conditions are shown through the use of some of the data. Tables are given for CHF data with a Freon coolant, for CHF data from test sections which only simulate a rod array, and for CHF data for transient situations.     

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.     

A method of evaluating local turbulent diffusivity of heat in a gridded fuel subassembly

In evaluating the turbulent diffusivity of heat associated with the coolant flow past a grid spacer within an FBR fuel subassembly, a heat diffusion technique is usually employed. However, measurement of subchannel bulk coolant temperature using thermocouples usually involves difficulty due to a steep and non-linear temperature gradient in the subchannels adjacent to a heater pin.A series solution of the heat conduction equation for the coolant flow in subchannels past a grid spacer and a heated section of a dummy fuel pin was derived under a slug flow approximation where the boundary conditions on dummy fuel pins were satisfied by means of the point-matching technique. The solution may be utilized in analyzing the turbulent diffusivity of heat within subchannel coolant flow as a function of distance from a grid spacer based on the measured temperature distribution on the wall of dummy fuel pins, which may be obtained without affecting the subchannel coolant temperature.In an illustrative example, the turbulent diffusivity of heat was most exaggerated at about 50 mm beyond a grid spacer and was approximately five times larger than the corresponding diffusivity without a grid spacer.     

Thermal-hydraulic modeling and analysis of a tank in pool reactor for normal operation and loss of flow transient

A transient thermal-hydraulic model entitled Tank in Pool Reactor Thermal-Hydraulic Analysis (TPRTHA) has been developed to simulate the steady-state operation and loss of flow transient for a tank in pool type research reactor. The model solves the momentum equation, energy equation and general conduction equation in cylindrical coordinates in order to predict the coolant velocity, coolant temperature and fuel rod temperature distribution respectively. The analytical solution is utilized for steady-state calculation while the finite difference technique with implicit scheme is adopted for transient calculation. The model divides the active core into a specified axial regions and the fuel rod into a specified radial zones, then a nodal calculation is performed for both average and hottest rods with a chopped cosine shaped heat generation flux. The model also predicts the heat flux leading to onset of nucleate boiling and the critical heat flux as well as the flow inversion phenomenon. The model is used to simulate a 2 MW reactor with downward flow direction and different types of fuel bundles of different powers and different flow rates. The best-estimate thermal-hydraulic safety margins are determined and the model results are analyzed and discussed.     

Modelling of azimuthal effects arising from interaction between clad deformation and heat transfer under LB LOCA conditions

The study of clad ballooning and the subchannel blockage arising from it is an important part of LWR safety research. A large amount of experimental data has been obtained over the years which has shown the importance of thermohydraulic effects on the clad ballooning process. Thus the distribution of clad strain, both azimuthally and axially within a fuel bundle are critically dependent on the detailed temperature distribution. Azimuthal clad temperature variations arise because of variations around the rod of pellet-clad gap, neutron flux, coolant temperature and surface heat transfer. They determine the pattern of straining but are also strongly influenced by it.A simple 3D computer model RODSWELL is being developed specifically to investigate the interactions between deformation and heat transfer. It is being used to identify the accident conditions potentially capable of producing significant subchannel blockage and to quantify the various parameters needed in the design of meaningful simulation experiments.     

A subchannel analysis code MATRA-LMR for wire wrapped LMR subassembly

In sodium cooled liquid metal reactors design limits are imposed on the maximum temperatures of the cladding and fuel pins. Thus an accurate prediction of the core coolant/fuel temperature distribution is essential to LMR core thermal hydraulic design. The detailed subchannel thermal hydraulic analysis code MATRA-LMR is being developed for LMFBR core design and analysis based on COBRA-IV-I and MATRA. The major modifications and improvements implemented in MATRA-LMR are as follows: sodium property calculation subprogram, sodium coolant heat transfer correlations, and most recent pressure drop correlations. To assess the development status of this code, benchmark calculations were performed with the ORNL 19 pin tests and EBR-II seven-assembly SLTHEN calculation results. The calculation results of MATRA-LMR were compared to the measurements and to the SABRE4 and SLTHEN code calculation results, respectively. Finally, the major technical results of the conceptual design for the KALIMER U-10%Zr binary alloy fueled core have been compared with the calculations of the MATRA-LMR, SABRE4 and SLTHEN codes.