CFD analysis of a turbulent jet behavior induced by a steam jet discharge through a single hole in a subcooled water tank

A computational fluid dynamics (CFD) analysis for a turbulent jet flow induced by a steam jet discharged into a subcooled water pool was performed for 10 s of transients to investigate whether the currently available CFD codes can be suitably used as a tool to investigate the applicability of the existing semi-analytical correlations to a condensing jet-induced turbulent jet and to analyze the thermal-hydraulic behavior, such as global circulation and local hot spot, in a condensation pool for advanced light water reactors. As for the numerical experiment, a series of sensitivity calculations was conducted systematically to elucidate the major factors which can cause different analysis results by varying the mesh distributions, numerical models for a convection term and an eddy viscosity term. The effect of a difference in the velocity and the temperature distribution in a region between the sparger and the pool wall has not been observed in the afore-mentioned sensitivity calculations. The comparison of the CFD results with the test data shows that the CFD analysis does not accurately simulate the local phenomenon of a turbulent jet existing downstream of a steam jet. It was found that the value of the turbulent intensity at the inlet of the turbulent jet region is the most important factor because it can determine the boundary of a turbulent jet through a momentum diffusion process in a radial direction. The comparison of the CFD results with the test data shows that the CFD analysis can accurately simulate the local phenomenon of a turbulent jet existing downstream of a steam jet only when the CFD analysis reflects the physics of a turbulent jet.     

Debris transport evaluation during the blow-down phase of a LOCA using computational fluid dynamics

The performance of the emergency recirculation water sump under the influence of debris accumulation following a loss-of-coolant accident (LOCA) has long been of safety concern. Debris generation and transport during a LOCA are significantly influenced by the characteristics of the ejected coolant flow. One-dimensional analyses previously have been attempted to evaluate the debris transport during the blow-down phase but the transport evaluation still has large uncertainties. In this work, a computational fluid dynamics (CFD) analysis was utilized to evaluate small and fine debris transport during the blow-down phase of a pressurized water reactor, OPR1000. The coolant ejected from the ruptured hot-leg was assumed to expand in an isenthalpic process. The transport of small and fine debris was assumed to be dominated by water-borne transport, and the transport fractions for the upper and lower parts of the containment were quantified based on the CFD analysis. It was estimated that 73% of small and fine debris is transported to the upper part of the containment. This value is close to the values estimated by nuclear regulatory bodies of The United States and Korea using one-dimensional models while it shows a large discrepancy from the value suggested in the NEI 04-07 baseline analysis.     

RANS modeling for flow in nuclear fuel bundle in pressurized water reactors (PWR)

This paper presents use of Reynolds-averaged Navier-Stokes (RANS) based turbulence model for single-phase CFD analysis of flow in pressurized water reactor (PWR) assemblies. An open source code called OpenFoam was used for computational fluid dynamics (CFD) study using computational meshes generated using Shari Harpoon. The PWR assembly design used in this analysis represents a 5 × 5 pin design including structural grid equipped with mixing vanes. The design specifications used in this study were obtained from the experimental setup at Texas A&M University and the results obtained are used to validate the CFD software, algorithm, and the turbulence model used in this analysis.     

Development of a numerical analysis methodology for the multi-dimensional and multi-phase phenomena of a sodium–water reaction in an SFR steam generator

A numerical approach on the exothermic sodium–water reaction (SWR) in a SFR steam generator is carried out by using a commercial computational fluid dynamics (CFD) code. The applicability of the analysis models and the physical limitations of some codes was investigated to select the most powerful CFD code to simulate a chemically reacting flow with various phases and components. In order to model the phenomena, among the several chemical reaction models studied, the eddy dissipation model (EDM) was employed because the EDM is the proper one when the reaction rate is sufficiently high when compared to the flow transport time. Based on the basic understandings for the characteristics of the SWR phenomena and the capabilities of the CFD codes, the numerical analysis methodology for a SWR was developed and transient analyses up to 0.05 s and 0.1 s with a time step of 0.0001–0.0005 s were carried out with a consideration of the geometric effect. The vapor mass flow rate and the corresponding hydrogen production rate were also calculated and compared with the conventional one-dimensional analysis results. As a result, it was found that the multi-dimensional approach underestimates the hydrogen production rate by 17% when compared to the theoretical values, and the difference is mainly caused by a multi-dimensional effect of the chemical reaction. The analysis performed in this study presents detailed information on each phase and the components of the SWR process and it also reflects the realistic SWR phenomena well. In order to confirm the applicability of the methodology, a multi-dimensional analysis was also carried out for the 49 tube bundle condition, and it was found that the results of the analysis were satisfactory.     

Simulation of hydrogen mitigation in catalytic recombiner. Part-II: Formulation of a CFD model

This paper aims at formulation of a model compatible with CFD code to simulate hydrogen distribution and mitigation using a Passive Catalytic Recombiner in the Nuclear power plant containments. The catalytic recombiner is much smaller in size compared to the containment compartments. In order to fully resolve the recombination processes during the containment simulations, it requires the geometric details of the recombiner to be modelled and a very fine mesh size inside the recombiner channels. This component when integrated with containment mixing calculations would result in a large number of mesh elements which may take large computational times to solve the problem. This paper describes a method to resolve this simulation difficulty. In this exercise, the catalytic recombiner alone was first modelled in detail using the best suited option to describe the reaction rate ( [Prabhudharwadkar et al., 2005] and [Prabhudharwadkar et al., 2011]). A detailed parametric study was conducted, from which correlations for the heat of reaction (hence the rate of reaction) and the heat transfer coefficient were obtained. These correlations were then used to model the recombiner channels as single computational cells providing necessary volumetric sources/sinks to the energy and species transport equations. This avoids full resolution of these channels, thereby allowing larger mesh size in the recombiners. The above mentioned method was successfully validated using both steady state and transient test problems and the results indicate very satisfactory modelling of the component.     

Non-intrusive experimental investigation of flow behavior inside a 5 × 5 rod bundle with spacer grids using PIV and MIR

The validity of the simulation results from computational fluid dynamics (CFD) is still under scrutiny. Some existing CFD closure models for complex flow produce results that are generally recognized as being inaccurate. Development of improved models for complex flow simulation requires an improved understanding of the detailed flow structure evolution with dynamic interaction of the flow multi-scales. Thus, the goal of this work is to contribute to a better understanding of presupposed and existent events that could affect the safety of nuclear power plants. The fundamental phenomena of fluid flow in rod bundles with spacer grids can be elucidated by using state-of-the-art measurement techniques. This study aims to develop an experimental data base with high spatial and temporal resolution of fluid flow velocity inside a 5 × 5 rod bundles with spacer grids. The full-field detailed data base is intended to validate CFD codes at various temporal-spatial scales. Measurements are carried out using dynamic particle image velocimetry (DPIV) technique inside an optically transparent rod bundle utilizing the matching index of refraction (MIR) approach. This work presents full field velocity vectors and turbulence statistics for a rod bundle under single phase flow conditions.     

Investigation of appropriate mesh size and solid angle number for CFD simulating the characteristics of pool fires_with experiments assessment

In order to implement NFPA 805 in the performance-based fire design for nuclear power plants (NPPs), a reliable computational fluid dynamics (CFD) fire model is needed. However, numerical treatments including mesh size and number of solid angles significantly influence the accuracy of CFD predicted results of the thermal-hydraulic behaviors involved in fire scenarios. Therefore, the majority of this paper is to investigate appropriate mesh size and solid angle number used for CFD simulating the characteristics of pool fires. Based on the present sensitivity studies of different mesh sizes and solid angle numbers, appropriate numerical treatment could be selected by comparing the predicted results of flame height and radiative heat flux for the pool fires. Two experiments with pan sizes of 20 and 30 cm, respectively, are also conducted to assess the CFD predicted results obtained using these selected mesh size and number of solid angles. Good agreement between the experiments and predictions clearly shows that the optima mesh resolution for the flame height and radiative heat flux is at the normalized mesh size (R*, ratio of the mesh size to the flame characteristic diameter) = 0.077 and 500 solid angles for thermal radiation are sufficient to reasonably capture the radiative characteristics from the pool fires with the size less than 30 cm. Herein, solid angle means included angle between rays from thermal radiator.     

A second order turbulence model based on a Reynolds stress approach for two-phase boiling flow and application to fuel assembly analysis

High-thermal performance PWR (pressurized water reactor) spacer grids require both low-pressure loss and high critical heat flux (CHF) properties. Numerical investigations on the effect of angles and position of mixing vanes and to understand in more details the main physical phenomena (wall boiling, entrainment of bubbles in the wakes, recondensation) are required.In the field of fuel assembly analysis or design by means of CFD codes, the overwhelming majority of the studies are carried out using two-equation Eddy Viscosity Models (EVM), especially the standard K-? model, while the use of Reynolds Stress Transport Models (RSTM) remains exceptional.The simulation of swirling flow generated by the mixing vanes plays an important role for the prediction of the CHF for the fuel assemblies. For this reason, according to [14] and [Mimouni et al., 2009b], rotation effects and RSTM model are more specifically addressed in the paper.Before comparing performance of EVM and RSTM models on fuel assembly geometry, we performed calculations with simpler geometries, the DEBORA case and the ASU-annular channel case. ASU-annular channel case has already been addressed in [14] and [Mimouni et al., 2009b].Then, a geometry closer to actual fuel assemblies is considered. It consists of a rectangular test section in which a 2 × 2 rod bundle equipped with a simple spacer grid with mixing vanes is inserted. The influence of the turbulence model on target variables linked to CHF limitation will be discussed. Moreover, the sensitivity to the mesh refinement will be particularly examined. The study of this case is a further step towards the modelling of the two-phase boiling flow in real-life grids and rod bundles.     

Study on Solid-Liquid Two-Phase Flow on PWR Sump Clogging Issue

It has been a concern that sump screen clogging would occur in pressurized water reactors (PWRs) in the case of a loss-of-coolant accident (LOCA), because two-phase jet flow would strip off thermal insulation from the piping and wash down the broken and fragmented debris to sump screens. It is necessary for the evaluation of the effectiveness of sump screens to estimate the amount of transported debris from a break position to sumps. In general, conservative logic trees have been used to determine debris transport rates. Realistic debris transport evaluation is useful for considering measures and rational decision making in licensing. The purpose of this study is to develop a debris transport evaluation model and to apply the model to this issue. We developed a solid-liquid multiphase model that is capable of simulating debris transport, settling, and resuspension. The model is able to treat solid particles of different sizes, which are smaller than uniform-sized liquid particles. This approach contributes to reducing the calculation cost in a large-scale simulation. The model and a turbulence model were implemented into a code based on the moving particle semi-implicit (MPS) method. Several open-channel hydraulic experiments with fibrous debris were conducted. The code named SANSUI 2.0 was validated by the comparison of the analytical results with experiments. This method was applied to the debris transport analysis of a full-scale PWR containment vessel floor, and the debris transport behavior was evaluated.     

Numerically simulating the thermal-hydraulic characteristics within the fuel rod bundle using CFD methodology

With the dramatic progress in the computer processing power, computational fluid dynamics (CFD) methodology can be applied in investigating the detailed knowledge of thermal-hydraulic characteristics in the rod bundle, especially with the spacer grid. These localized information, including flow, turbulence, and heat transfer characteristics, etc., can assist in the design and the improvement of rod bundles for nuclear power plants. In this paper, a three-dimensional (3D) CFD model with the Reynolds stresses turbulence model is proposed to simulate these characteristics within the rod bundle and subsequently to investigate the effects of different types of grid on the turbulent mixing and heat transfer enhancement. Two types of grid designs are used herein, including the standard grid and split-vane pair one, respectively. Based on the CFD simulations, the secondary flow can be reasonably captured in the rod bundle with the grid. The split-vane pair grid would enhance both the flow mixing and the heat transfer capability more than the standard grid does, as clearly shown in the simulation results. In addition, compared with the results of experiment and correlation, the present predicted result for the Nusselt (Nu) number distribution downstream the grid shows reasonable agreement for the standard grid design. However, there is discrepancy in the decay trend of Nu number between the prediction and measurement for the split-vane pair gird. This would be improved by adopting the finer mesh (y⁺ < 1) simulation and Low-Reynolds form turbulence model, which is our future research work.     

Evaluation of computational fluid dynamic methods for reactor safety analysis (ECORA)

The objective of the ECORA project is the evaluation of computational fluid dynamics (CFD) software for reactor safety applications, resulting in best practice guidelines (BPG) for an efficient use of CFD for reactor safety problems. The project schedule is as follows: (i) establishment of BPGs for use of CFD codes, for judgement of CFD calculations and for assessment of experimental data; (ii) assessment of CFD simulations for three-dimensional flows in LWR primary systems and containments; (iii) quality-controlled CFD simulations for selected UPTF and SETH PANDA test cases; and (iv) demonstration of CFD code customisation for PTS analysis by implementation and validation of improved turbulence and two-phase flow models.The project started in October 2001 and is for a period of 36 months. The project consortium consists of 12 partners combining thermal-hydraulic experts, code developers, safety experts and engineers from nuclear industry and research organizations. At mid-term, the following results were achieved: (i) BPGs are available for simulations of reactor safety relevant flows. These BPGs have found interest in the European projects FLOMIX-R, ASTAR and ITEM; (ii) important flow phenomena for PTS and containment flows have been identified; (iii) experimental data featuring these phenomena have been selected and described in a standardised manner suitable for simulation with CFD methods; (iii) surveys of existing CFD calculations and experimental data for containment and primary loop flows have been performed and documented; (iv) first results for simulations of PTS-relevant single-phase and two-phase flow cases are available.Documentation is available via the internet at http://domino.grs.de/ecora/ecora.nsf. The models developed within the project are implemented in industrial and commercial CFD software packages and are therefore accessible by industry and research institutions.     

Erratum to “Experimental study on the air/water counter-current flow limitation in a model of the hot leg of a pressurized water reactor”

This paper replaces the paper published in the journal by Deendarlianto et al. (2008). Because of an error in the implementation of the air flow meter some of the data given by Deendarlianto et al. (2008) are wrong. They are corrected within the present paper. The general results and conclusions remain unchanged.An experimental investigation on the air/water counter-current two-phase flow in a horizontal rectangular channel connected to an inclined riser has been conducted. This test-section representing a model of the hot leg of a pressurized water reactor is mounted between two separators in a pressurized experimental vessel. The cross-section and length of the horizontal part of the test-section are (0.25 m × 0.05 m) and 2.59 m, respectively, whereas the inclination angle of the riser is 50°. The flow was captured by a high speed camera in the bended region of the hot leg, delivering a detailed view of the stratified interface as well as of dispersed structures like bubbles and droplets. Counter-current flow limitation (CCFL), or the onset of flooding, was found by analyzing the water levels measured in the separators. The counter-current flow limitation is defined as the maximum air mass flow rate at which the discharged water mass flow rate is equal to the inlet water mass flow rate.From the high-speed observations it was found that the initiation of flooding coincides with the formation of slug flow. Furthermore, a slight hysteresis was noticed between flooding and deflooding. The CCFL data was compared with similar experiments and empirical correlations available in the literature. Therefore, the Wallis-parameter was calculated for the rectangular cross-sections by using the channel height as length, instead of the diameter. The agreement of the CCFL curve is good, but the zero liquid penetration was found at lower values of the Wallis parameter than in most of the previous work. This deviation can be attributed to the special rectangular geometry of the hot leg model of HZDR, since the other investigations were done for pipes.     

CHF experiment and CFD analysis in a 2 × 3 rod bundle with mixing vane

In this study, the CHF enhancement using various mixing vanes is evaluated and the flow characteristics are investigated through the CHF experiments and CFD analysis.CHF tests were performed using 2 × 2 and 2 × 3 rod bundles and with R-134a as the working fluid. The test section geometry was identical to that of commercial PWR fuel assembly not including the heated length (1.125 m) and number of fuel rods. From the CHF tests, it was found that the CHF enhancement using mixing vanes under higher mass flux (1400 kg/m² s) and lower pressure (15 bar) conditions is larger than the CHF enhancements under other conditions. Among the mixing vanes used in this study, the swirl vane showed the best performance under relatively low pressure (15 bar) and mass flux (300-1000 kg/m² s) conditions and the hybrid vane performed best near the PWR operating conditions.The detailed flow characteristics were also investigated by CFD analysis using the same conditions as the CHF tests. To calculate the subcooled boiling flow, the wall partitioning model was applied to the wall boundary and various two-phase parameters were also considered. The reliability of the CFD analysis in the boiling analysis was confirmed by comparing the average void fractions of the analysis and the experiments: the results agreed well. From the CFD analysis, the void fraction flattening as a result of the lateral velocity induced by the mixing vane was observed. By the lateral motion of the liquid, the void fraction in the near wall was decreased and that of the core region was increased resulting in the void fraction flattening. The decrease of the void fraction in the near wall region promoted liquid supply to the wall and consequently the CHF increased. For the quantification of the void flatness, an index was developed and the applicability of the index in the CHF assessment was confirmed.     

One-dimensional three-field model of condensation in horizontal countercurrent flow with supercritical liquid velocity

A one-dimensional three-field model was developed to predict the flow of liquid and vapor that results from countercurrent flow of water injected into the hot leg of a PWR and the oncoming steam flowing from the upper plenum. The model solves the conservation equations for mass, momentum, and energy in a continuous-vapor field, a continuous-liquid field, and a dispersed-liquid (entrained-droplet) field. Single-effect experiments performed in the upper plenum test facility (UPTF) of the former SIEMENS KWU (now AREVA) at Mannheim, Germany, were used to validate the countercurrent flow limitation (CCFL) model in case of emergency core cooling water injection into the hot legs. Subcooled water and saturated steam flowed countercurrent in a horizontal pipe with an inside diameter of 0.75 m. The flow of injected water was varied from 150 kg/s to 400 kg/s, and the flow of steam varied from 13 kg/s to 178 kg/s. The subcooling of the liquid ranged from 0 K to 104 K. The velocity of the water at the injection point was supercritical (greater than the celerity of a gravity wave) for all the experiments. The three-field model was successfully used to predict the experimental data, and the results from the model provide insight into the mechanisms that influence the flows of liquid and vapor during countercurrent flow in a hot leg. When the injected water was saturated and the flow of steam was small, all or most of the injected water flowed to the upper plenum. Because the velocity of the liquid remained supercritical, entrainment of droplets was suppressed. When the injected water was saturated and the flow of steam was large, the interfacial shear stress on the continuous liquid caused the velocity in the liquid to become subcritical, resulting in a hydraulic jump. Entrainment ensued, and the flow of liquid to the end of the hot leg was greatly reduced.The influence of condensation on the transition from supercritical to subcritical flow as observed in the experimental data is also predicted with the three-field model. When the injected water was subcooled, condensation on the flow of continuous liquid caused a reduction in the flow of vapor and, consequently, a reduction in the interfacial shear stress. Therefore, the flow of liquid remained supercritical to the end of the hot leg at the upper plenum. The entire flow of injected water flowed to the end of the hot leg at higher flows of steam when the injected water was subcooled than when it was saturated. When the flow of vapor was large enough to cause a hydraulic jump in the subcooled liquid, the rate of entrained droplets was greatly increased. The interfacial surface area of the droplets was several orders of magnitude greater than for the continuous-liquid field, and condensation rate on the droplet field was also several orders of magnitude greater. When the flow of vapor from the upper plenum was at its greatest, most of the flow in the continuous liquid was entrained before reaching the upper plenum. The large flow of subcooled droplets caused three-quarters of the steam to condense.     

A thermal-fluid assessment of a cooled-vessel concept for a VHTR

Flow distribution and thermal analyses of a conceptual design of a cooled vessel for a very high temperature reactor (VHTR), which has a forced vessel cooling with an internal coolant path through a permanent side reflector, have been performed. A computational fluid dynamics (CFD) code was employed to investigate flow distributions at inlet and upper plenums of the proposed cooled-vessel concept. Thermal-fluid analyses of the cooled vessel during a normal operation were carried out by using the CFD code with the boundary conditions provided by the GAMMA system analysis code. The transient analyses during postulated accidents were conducted by the GAMMA code itself. According to the results, the flow deviation at the riser holes due to a change of the inlet flow path to the core inlet is about ±20% which results in about a 3-7% core flow deviation from the average value depending on the upper plenum height. The pressure drops in the inlet and upper plenums are estimated to be from 13 to 25 kPa with a change of the upper plenum height. A cooling flow of more than 4 kg/s is sufficient to maintain the RPV temperature within the required limit during a normal operation. Transient analysis reveals that the reactor vessel is exposed to a temperature above its limit of 371 °C but this duration is shorter than the allowable time for a creep region with a sufficient safety margin. The results suggest that the cooled-vessel concept considered in this paper has the potential to be used for a VHTR but further and more detailed studies are required to realize the proposed concept.