• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.482-487
• /
• 2001

A CFD study was conducted to evaluate the nuclear fuel assembly coolant mixing that is promoted by the flow-mixing vanes on the grid spacer. Four mixing vanes (split vane, swirl vane, twisted vane, hybrid vane) were chosen in this study. A single subchannel of one grid span is modeled using the flow symmetry. The three mixing vanes other than swirl vane generate a large crossflow between the subchannels and a skewed elliptic swirling flow in the subchannel near the grid spacer. The swirl vane induces a circular swirling flow in the subchannel and a negligible crossflow. The split vane and the twisted vane were predicted to result in relatively larger pressure drop across the grid spacer. Since the average turbulent kinetic energy in the subchannel rapidly decreases to a fully developed level downstream of the spacer, turbulent mixing caused by the mixing vanes appears to be not as effective as swirling flow mixing in the subchannel. In summary, the CFD analysis represented the overall characteristics of coolant mixing well in a nuclear fuel assembly with the flow mixing vanes on the grid spacer. The CFD study is therefore quite useful for the development of an advanced flow-mixing vane.

• Nuclear Engineering and Technology
• /
• v.48 no.2
• /
• pp.376-385
• /
• 2016

A hexagonally arrayed 37-pin wire-wrapped rod bundle has been chosen to provide the experimental data of the pressure loss and flow rate in subchannels for validating subchannel analysis codes for the sodium-cooled fast reactor core thermal/hydraulic design. The iso-kinetic sampling method has been adopted to measure the flow rate at subchannels, and newly designed sampling probes which preserve the flow area of subchannels have been devised. Experimental tests have been performed at 20-115% of the nominal flow rate and $60^{\circ}C$ (equivalent to Re ~ 37,100) at the inlet of the test rig. The pressure loss data in three measured subchannels were almost identical regardless of the subchannel locations. The flow rate at each type of subchannel was identified and the flow split factors were evaluated from the measured data. The predicted correlations and the computational fluid dynamics results agreed reasonably with the experimental data.

• Nuclear Engineering and Technology
• /
• v.41 no.3
• /
• pp.295-306
• /
• 2009

The subchannel grade void distributions in the NUPEC (Nuclear Power Engineering Corporation) BFBT (BWR Full-Size Fine-Mesh Bundle Tests) facility were evaluated with the subchannel analysis code MATRA and the system code MARS. Fifteen test series from five different test bundles were selected for an analysis of the steady-state subchannel void distributions. Two transient cases, a turbine trip without a bypass as a typical power transient and a re-circulation pump trip as a flow transient, were also chosen for this analysis. It was found that the steady-state void distributions calculated by both the MATRA and MARS codes coincided well with the measured data in the range of thermodynamic qualities from 5% to 25%. The results of the transient calculations were also similar and were highly feasible. However, the computational aspects of the two codes were clearly different.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.79-87
• /
• 2002

The present work analyzed the effect of mixing vane shape on the flow structure and heat transfer downstream of mixing vane in a subchannel of fuel assembly, by obtaining velocity and pressure fields, turbulent intensity, flow-mixing factors, heat transfer coefficient and friction factor using three-dimensional RANS analysis. NJl5, NJ25, NJ35, NJ45, which were designed by the authors, were tested to evaluate the performances in enhancing the heat transfer. Standard $\kappa-\epsilon$ model is used as a turbulence closure model, and, periodic and symmetry conditions are set as boundary conditions. The flow blockage ratio is kept constant, but the twist angle of mixing vane is changed. The results with three turbulence models( $\kappa-\epsilon$, $\kappa-\omega$, RSM) were compared with experimental data.

• Nuclear Engineering and Technology
• /
• v.27 no.4
• /
• pp.491-502
• /
• 1995

In accurate prediction of cross-flow based on detailed knowledge of the velocity field in subchannels of a nuclear fuel assembly is of importance in nuclear fuel performance analysis. In this study, the low-Reynolds number k-$\varepsilon$ turbulence model has been adopted in too adjacent subchannels with cross-flow. The secondary flow is accurately estimated by the anisotropic algebraic Reynolds stress model. This model was numerically calculated by the finite element method and has been verified successfully through comparison with existing experimental data. Finally, with the numerical analysis of the velocity Held in such subchannel domain, an analytical correlation of the lateral loss coefficient is obtained to predict the cross-flow rate in subchannel analysis codes. The correlation is expressed as a function of the ratio of the lateral How velocity to the donor subchannel axial velocity, recipient channel Reynolds number and pitch-to-diameter.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2630-2636
• /
• 1995

A theoretically based prediction for the determination of the subchannel friction factor at low pitch to the rod diameter ratio (P/D < 1.2) in the bare rod bundle flow has been developed. The present model assumes the validity of the Law of Wall over the entire flow area. The algebraic form of the Law of the Wall is integrated over the entire flow area and the local friction velocity variation along the rod periphery is considered in this study. The present method is applied to the rod bundles with P/D < 1.2, and the prediction results show good agreement with the available experimental data.

• Nuclear Engineering and Technology
• /
• v.31 no.3
• /
• pp.314-327
• /
• 1999

A subchannel analysis code MATRA applicable to PWRs and ALWRs has been developed to be run on an IBM PC or HP WS based on the existing CDC CYBER mainframe version of COBRA-Rf-1. This MATRA code is a thermal-hydraulic analysis code based on the subchannel approach for calculating the enthalpy and How distribution in fuel assemblies and reactor cores for both steady-state and transient conditions. HATRA has been provided with an improved structure, various functions, and models to give more convenient user environment and to enhance the code accuracy. Among them, the pressure drop model has been improved to be applied to non-square-lattice rod arrays, and the models for the lateral transport between adjacent subchannels have been improved to enhance the accuracy in predicting two-phase flow phenomena. The predictions of MATRA were compared with the experimental data on the flow and enthalpy distribution in some sample rod-bundle cases to evaluate the performance of MATRA. All the results revealed that the predictions of MATRA were better than those of COBRA-IV-I.

• Nuclear Engineering and Technology
• /
• v.51 no.3
• /
• pp.665-675
• /
• 2019

The wire-wrapped fuel bundle is an assembly design in a sodium-cooled fast reactor. A wire spacer is used to maintain a constant gap between rods and to enhance the mixing of coolants. The wire makes the flow complicated by creating a sweeping flow and vortex flow. The vortex affects the flow field and heat transfer inside the subchannels. However, studies on vortices in this geometry are limited. The purpose of this research is to investigate the vortex flow created in the wire-wrapped fuel bundle. For analysis, a RANS-based numerical analysis was conducted for a 37-pin geometry. The sensitivity study shows that simulation with the shear stress transport model is appropriate. For the case of Re of 37,100, the mechanisms of onset, periodicity, and rotational direction were analyzed. The vortex structures were reconstructed in a three-dimensional space. Vortices were periodically created in the interior subchannel three times for one wire rotation. In the edge subchannel, the largest vortex occurred. This large vortex structure blocked the swirl flow in the peripheral region. The small vortex formed in the corner subchannel was negligible. The results can help in understanding the flow field inside subchannels with sweeping flow and vortex structures.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.8-15
• /
• 2003

Three nonlinear κ-ε models with the wall function method are applied to the fully developed turbulent flow in a triangular subchannel of a bare rod bundle. Typical predicted quantities such as axial and secondary velocities, turbulent kinetic energy and wall shear stress are compared in details both qualitatively and quantitatively with both each other and experimental data. The nonlinear κ-ε models by Speziale[1] and Myong and Kasagi[2] are found to be capable of predicting accurately noncircular duct flows involving turbulence-driven secondary motion. The nonlinear κ-ε model by Shih et aL.[3] adopted in a commercial code is found to be unable to predict accurately noncircular flows with the prediction level of secondary flows one order less than that of the experiment.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.71-77
• /
• 2016

The CUPID code is a transient, three-dimensional, two-fluid, thermal-hydraulic code designed for a component-scale analysis of nuclear reactor components. The primary objective of this study is to assess the applicability of CUPID to single-phase turbulent flow analyses of $2{\times}2$ rod bundle subchannel. The bulk velocity at the inlet varies from 1.0 m/s up to 2.0 m/s which is equivalent to the fully turbulent flow with the range of Re=12,500 to 25,000. Adiabatic single-phase flow is assumed. The velocity profile at the exit region is quantitatively compared with both experimental measurement and commercial CFD tool. Three different boundary conditions are simulated and quantitatively compared each other. The calculation results of CUPID code shows a good agreement with the experimental data. It is concluded that the CUPID code has capability to reproduce the turbulent flow behavior for the $2{\times}2$ rod bundle geometry.