CHF experiment and CFD analysis in a 2 × 3 rod bundle with mixing vane
Introduction
In a PWR nuclear fuel assembly, a spacer grid is used to maintain a gap between the fuel rods. Maintaining the gap is to thermal-hydraulically support uniform heat transfer and to thermal-hydraulically suppress the local occurrence of CHF. If a mixing vane is attached to the spacer grid, the functions of the spacer grid are further reinforced, and the CHF is greatly enhanced. For example, a commercial nuclear fuel assembly with a mixing vane can enhance the CHF by up to 10% and the economic benefits are great due to this enhancement. Despite the demerits of a mixing vane, such as the increase of the pressure drop and FIV and the debris problem, most nuclear fuel development includes the development of a mixing vane because the economical benefits are too great to ignore.
In the existing research, the enhancement of CHF using a mixing vane is one of many considerations in CHF studies on rod bundle geometry. ERPI, CE, and B&W correlations, the well-known correlations for rod bundles, include the effect of a spacer grid, a guide thimble, the non-heating area locations, and so on. In this case, the effect of the mixing vane is taken into consideration in the grid factor (Chang and Baek, 1997).
Fortini and Veloso (2002) and Lee (2000) have attempted to predict the CHF in rod bundles by applying the correction factor to the CHF correlation in a circular tube. They also included the effects of the mixing vanes as a correction factor. The reason that the studies have propriety is because the shapes of the existing mixing vanes are similar to each other and the CHF enhancements were not vastly different. Therefore, a correction factor that does not significantly affect the performance of the prediction should be used. Recently, Chen et al. (2004) and Akhtar et al. (2006) performed a CHF experiment using rod bundles with refrigerants: they showed that using a refrigerant enables proper results to be obtained through fluid-to-fluid modeling.
Some research has investigated the detailed mechanism of the CHF enhancement caused by the mixing vane. Generally, it is well-known that the CHF enhancement through the mixing vane is due to the generation of swirl and cross flows, and the increase of the turbulent intensity. Chung et al. (1996) summarized the effects as shown in Table 1. de Crecy (1994) conducted a radial non-uniform heating CHF experiment to investigate the change of location and the enhancement of the CHF value through the implementation of a mixing vane. According to Crecy's study, the CHF appeared at random radial locations and the CHF value increased in rod bundles with a mixing vane.
Mixing vanes in nuclear fuel assemblies are used to increase the thermal mixing and increase the DNB margin. To date, in most research on mixing vanes, the shape of the mixing vane is not important; rather the research only considers one parameter, such as the mixing vane factor, in the CHF correlation. However, due to the economical benefits gained by incorporating a mixing vane and the technical improvements in the CFD code, much more thermal-hydraulic research on mixing vanes has been undertaken recently.
Using the CFD analysis of a single phase flow, Ikeda et al. (2006) predicted the location of the DNB occurrence. They established that the location of the maximum local enthalpy in a single-phase analysis corresponded to the location of the DNB occurrence. However, the analysis could not directly predict the CHF value and could only predict the location of the DNB. Krepper et al. (2007) predicted various two-phase variables in the subcooled flow boiling using a wall boiling model. This study shows the potential of predicting the CHF through the CFD analysis using the wall boiling model in a two-phase flow. The study establishes that the CFD analysis results in proper flow characteristics regarding the CHF prediction. However, the study only shows the analysis results for a single-phase flow: the CFD analysis should also be applied to a two-phase flow to examine the detailed mechanism of the CHF enhancement
The primary objective of this paper is to find the representative mixing factor to mean the performance of the mixing vane on the CHF. To achieve the primary objective, a CFD analysis was performed under identical conditions to those of the CHF test. The two-phase flow parameters under the CHF condition were calculated using the CFD analysis. The CHF test was performed in a 2 × 3 rod bundle with a simple spacer grid and various mixing vanes at KAIST R-134a Loop. The R-134a was a working fluid and, for CFD analysis, CFX 11.0 with a beta version of the wall boiling model was used. By comparing the experimental and numerical results, the representative mixing factor related to the CHF enhancement was obtained. An additional objective of this study was to investigate the CHF characteristics of various mixing vanes and the effect of varying design parameters on the mixing vane. This study also proposes the design direction and provides useful information for designing mixing vanes.
Section snippets
Experimental works
Fig. 1 shows the test loop for the experiment. The loop is closed loop and R-134a was used as the working fluid. The detailed specifications of the loop can be found in the experimental apparatus section of Shin et al.’s study (Shin et al., 2005).
Fig. 2 shows the axial and cross sectional geometry and structure of the 2 × 3 rod bundle test section. The distance between heater rods was 12.65 mm; the heater rod diameter was 9.5 mm; the distance between the heater rod and the insulated wall was 1.575
Two-fluid model
In a two-phase flow, the vapor and liquid phases are treated as continua. The phases may interact with each other through interfacial forces as well as heat and mass transfer across the phase interfaces. The two-fluid model solves two sets of mass, momentum, and energy conservation equations, which are written for each phase as:
(a) Continuity equation
(b) Momentum equation
(c) Energy equation
Experimental results
In the CHF experiment with a 2 × 3 rod bundle and a simple spacer grid, the trends of the CHF regarding the major parameters were investigated. Fig. 9 shows the parametric trends of the CHF regarding pressure, mass flux, and inlet subcooling. The CHF value increased as the mass flux and inlet subcooling increased, and the pressure decreased. This trend is similar to that in a circular tube (Bartolomei and Chanturiya, 1967).
Fig. 10 presents a comparison of the experimental data in this study with
Conclusions
The major conclusions from the experimental work and numerical analyses are as follows:
- (1)
The CHF enhancement using mixing vanes was enlarged under relatively high mass flux and low pressure conditions. Under a high mass flux condition, the enlarged swirl and cross flows led to the reduction of the wall void fraction and enthalpy near the wall. Under a low pressure condition, the effect of the swirl and cross flows strengthened due to the high density ratio (liquid/vapor).
- (2)
Among the mixing vanes
Reference (18)
Modeling capability of R134a for a critical heat flux of water in a vertical 5 × 5 rod bundle geometry
Int. J. Heat Mass Transfer
(2006)Fluid-to-fluid modeling of critical heat flux in 4 × 4 rod bundles
Nucl. Eng. Des.
(2004)- et al.
Effects of the spacer and mixing vanes on critical heat flux for low-pressure water at low-velocities
Int. Commun. Heat Mass Transfer
(1996) The effect of grid assembly mixing vanes on critical heat flux values and azimuthal location in fuel assemblies
Nucl. Eng. Des.
(1994)Single-phase CFD applicability for estimating fluid hot-spot locations in a 5 × 5 fuel rod bundle
Nucl. Eng. Des.
(2006)Temperature and concentration profiles in fully turbulent boundary layers
Int. J. Heat Mass Transfer
(1981)A generalized correlation of critical heat flux for the forced convection boiling in vertical uniformly heated round tubes
Int. J. Heat Mass Transfer
(1978)CFD modelling of subcooled boiling-Concept, validation and application to fuel assembly design
Nucl. Eng. Des.
(2007)A critical heat flux approach for square rod bundles using the 1995 Groeneveld CHF table and bundle data of heat transfer research facility
Nucl. Eng. Des.
(2000)
Cited by (65)
Numerical investigation on critical heat flux for 5 × 5 fuel assemblies with a wide range of working conditions
2024, Progress in Nuclear EnergyA comprehensive review of numerical and experimental research on the thermal-hydraulics of two-phase flows in vertical rod bundles
2024, International Journal of Heat and Mass TransferNumerical analysis on subcooled boiling in PWR coolant channel based on a modified multi-scale interface model
2023, Applied Thermal EngineeringExperimental investigation on critical heat flux enhancement by mixing-vanes spacer grid in a vertical annular geometry
2023, Nuclear Engineering and DesignNumerical study on the effects of guide thimble on critical heat flux in a 5 × 5 fuel assembly
2022, Annals of Nuclear EnergyNumerical investigation of the critical heat flux in a 5 × 5 rod bundle with multi-grid
2022, Nuclear Engineering and TechnologyCitation Excerpt :According to the research results, the Eulerian two-fluid model has been able to predict the DNB phenomenon in a simple structure with a wide range of working conditions (pressure range: 0.1–16 MPa, mass flux range: 500–4000 kg/m2s). Some scholars attempted to apply the Eulerian two-fluid model to the calculation of flow boiling and CHF in complex structures (such as a bundle channel with mixing vanes) [23–26]. According to the literature summarized in this article, most of these studies were case by case calculations and only confirmatory calculations were carried out for a single working condition.