scholarly journals Hydrodynamic behavior of a bare rod bunle

1977 ◽  
Author(s):  
Ιωάννης Μπαρτζής

Temperature distributions within the rod bundle of a nuclear reactor is of major importance in nuclear reactor design. However temperature information presupposes knowledge of the hydrodynamic behavior of the coolant which is the most difficult part of the problem due to complexity of the turbulence phenomena. In the present work a 2-equation turbulence model - a strong candidate for analyzing actual three dimensional turbulent flows - has been used to predict fully developed flow of infinite bare rod bundle of various aspect ratios (P/D). The model has been modified to take into account anisotropic effects of eddy viscosity. Secondary flow calculations have been also performed although the mo^el seems to be too rough to predict the secondary flow correctly. Heat Transfer calculations have been performed to confirm the importance of anisotropic viscosity in temperature predictions. All numerical calculations for flow and heat have been performed by two computer codes developed in the present work which were based on the TEACH code [71]· Also experimental measurements of the distribution of axial velocity, turbulent axial velocity, turbulent kinetic energy and radial Reynoldsstresses were performed in the developing and fully developed regions. A 2-channel Laser Doppler Anemometer working on the Reference mode with forward scattering was used to perform the measurements in a simulated interior subchannel of a triangular rod array with P/D=1.124. Comparisons between the analytical results and the results of this experiment as well as other experimental data in rod bundle array available in literature are presented. The predictions are in good agreement with the results for the high Reynolds numbers.

1979 ◽  
Vol 101 (4) ◽  
pp. 628-634 ◽  
Author(s):  
J. G. Bartzis ◽  
N. E. Todreas

Temperature distribution within the rod bundle of a nuclear reactor is of major importance in nuclear reactor design. However temperature information presupposes knowledge of the hydrodynamic behavior of the coolant which is the most difficult part of the problem due to the complexity of the turbulence phenomena. In the present work a two equation turbulence model—a strong candidate for analyzing actual three dimensional turbulent flows—has been used to predict fully developed flow of infinite bare rod bundle of various aspect ratios (P/D). The model has been modified to take into account anisotropic effects of eddy viscosity. Secondary flow calculations have been also performed although the model seems to be too rough to predict the secondary flow correctly. Heat transfer calculations have been performed to confirm the importance of anisotropic viscosity in temperature predictions. Experimental measurements of the distribution of axial velocity, turbulent axial velocity, turbulent kinetic energy and radial Reynolds stresses were performed in the developing and fully developed regions. A two channel Laser Doppler Anemometer working in the reference mode with forward scattering was used to perform the measurements in a simulated interior subchannel of a triangular rod array with P/D = 1.124. Comparisons between the analytical results and the results of this experiment as well as other experimental data in rod bundle arrays available in the literature are presented. The predictions are in good agreement with the results for high Reynolds numbers.


2021 ◽  
Vol 925 ◽  
Author(s):  
Wei He ◽  
Sebastian Timme

This article uses triglobal stability analysis to address the question of shock-buffet unsteadiness, and associated modal dominance, on infinite wings at high Reynolds number, expanding upon recent biglobal work, aspiring to elucidate the flow phenomenon's origin and characteristics. Infinite wings are modelled by extruding an aerofoil to finite aspect ratios and imposing a periodic boundary condition without assumptions on spanwise homogeneity. Two distinct steady base flows, spanwise uniform and non-uniform, are analysed herein on straight and swept wings. Stability analysis of straight-wing uniform flow identifies both the oscillatory aerofoil mode, linked to the chordwise shock motion synchronised with a pulsation of its downstream shear layer, and several monotone (non-oscillatory), spatially periodic shock-distortion modes. Those monotone modes become outboard travelling on the swept wing with their respective frequencies and phase speeds correlated with the sweep angle. In the limiting case of very small wavenumbers approaching zero, the effect of sweep creates branches of outboard and inboard travelling modes. Overall, triglobal results for such quasi-three-dimensional base flows agree with previous biglobal studies. On the contrary, cellular patterns form in proper three-dimensional base flow on straight wings, and we present the first triglobal study of such an equilibrium solution to the governing equations. Spanwise-irregular modes are found to be sensitive to the chosen aspect ratio. Nonlinear time-marching simulations reveal the flow evolution and distinct events to confirm the insights gained through dominant modes from routine triglobal stability analysis.


1976 ◽  
Vol 98 (2) ◽  
pp. 262-268 ◽  
Author(s):  
P. Carajilescov ◽  
N. E. Todreas

Reactor fuel elements generally consist of rod bundles with the coolant flowing axially through the bundles in the space between the rods. Heat transfer calculations form an important part in the design of such elements, which can only be carried out if information of the velocity field is available. A one-equation statistical model of turbulence is applied to compute the detailed description of velocity field (axial and secondary flows) and the wall shear stress distribution of steady, fully developed turbulent flows with incompressible, temperature-independent fluid, flowing through triangular arrays of rods with different aspect ratios (P/D). Also experimental measurements of the distributions of the axial velocity, turbulence kinetic energy, and Reynolds stresses were performed using a laser Doppler anemometer (LDA), operating in a “fringe” mode with forward scattering, in a simulated interior subchannel of a triangular rod array with P/D = 1.123 and L/DH = 77. From the experimental results, a new mixing length distribution is proposed. Comparisons between the analytical results and the results of this experiment as well as other experimental data available in the literature are presented. The results are in good agreement.


1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


Author(s):  
Vladimir Kriventsev ◽  
Hiroyuki Ohshima ◽  
Akira Yamaguchi ◽  
Hisashi Ninokata

A new model of turbulence is proposed for the estimation of Reynolds stresses in turbulent fully-developed flow in a wall-bounded straight channel of an arbitrary shape. Ensemble-averaged Navier-Stokes, or Reynolds, equations are considered to be sufficient and practical enough to describe the turbulent flow in complex geometry of rod bundle array. We suggest the turbulence is a process of developing of external perturbations due to wall roughness, inlet conditions and other factors. We also assume that real flows are always affected by perturbations of any possible scale lower than the size of the channel. Thus, turbulence can be modeled in the form of internal or “turbulent” viscosity. The main idea of a Multi-Scale Viscosity (MSV) model can be expressed in the following phenomenological rule: A local deformation of axial velocity can generate the turbulence with the intensity that keeps the value of the local turbulent Reynolds number below some critical one. Therefore, in MSV, the only empirical parameter is the critical Reynolds number. From analysis of dimensions, some physical explanations of Reynolds number are possible. We can define the local turbulent Reynolds number in two ways: i) simply as Re = ul/v, where u is a local velocity deformation within the local scale l and v is total accumulated molecular and turbulent viscosity of all scales lower then 1. ii) Re = K/W, where K is kinetic energy and W is work of friction/dissipation forces. Both definitions above have been implemented in the calculation of samples of basic fully-developed turbulent flows in straight channels such as a circular tube and annular channel. MSV has been also applied to prediction of turbulence-driven secondary flow in elementary cell of the infinitive hexagonal rod array. It is known that the nature of these turbulence-driven motions is originated in anisotropy of turbulence structure. Due to the lack of experimental data up to date, numerical analysis seems to be the only way to estimate intensity of the secondary flows in hexagonal fuel assemblies of fast breeder reactors (FBR). Since MSV can naturally predict turbulent viscosity anisotropy in directions normal and parallel to the wall, it is capable to calculate secondary flows in the cross-section of the rod bundle. Calculations have shown that maximal intensity of secondary flow is about 1% of the mean axial velocity for the low-Re flows (Re = 8170), while for higher Reynolds number (Re = 160,100) the intensity of secondary flow is as negligible as 0.2%.


2005 ◽  
Vol 9 (3) ◽  
pp. 58-63
Author(s):  
L.X. L.X. ◽  
B. Wang ◽  
F.H. She

A three-dimensional model is developed to evaluate the effect of secondary flow generated from strongly bent duct profiles and turbulent flow of high Reynolds number on fibre movement in a bent channel. The fibre configuration is more complex in a three-dimensional model with the introduction of secondary flow. The strategies of mesh generation for threedimensional problems are discussed. The flow characteristic in the transfer-channel of a rotor spinning machine is predicted.


Kerntechnik ◽  
2021 ◽  
Vol 86 (3) ◽  
pp. 210-216
Author(s):  
M. P. Sharma ◽  
A. Moharana

Abstract Subchannel analysis codes are widely used for the thermal-hydraulic design of nuclear reactor rod bundle. The effectiveness of subchannel analysis codes depends on turbulent mixing between these subchannels. Turbulent mixing has no direct contribution to the axial mass flow rate through subchannel but it will cause exchange of momentum and energy between the neighboring subchannels. Thus, it is important to evaluate the turbulent mixing coefficient for reactor rod bundle as it is a significant factor in the lateral energy and momentum equation for subchannel analysis codes like COBRA IIIC, COBRA-IV and MATRA LMR-FB. With the rapid developments in computational fluid dynamics and computer performance, three-dimensional analyses of turbulent flows occurring in the nuclear rod bundle have become more prominent. Several numerical analyses have already been attempted to investigate the flow behavior in rod bundles of different reactors. Much of these are dedicated to find out the structure of turbulence in rod bundle but a few analyses has been done to evaluate the magnitude of the turbulent mixing coefficient. In view of this, CFD analyses were carried out to determine the turbulent mixing coefficient in the simulated sub-channels of the reactor rod bundle. Previous studies on the structure of turbulence reveals that it is highly anisotropic. Hence, the Reynolds Stress Model (RSM), finer mesh and near wall distance ( y + ≤ 2) is required to capture turbulent mixing phenomena. The validation of results is done by comparing with subchannel mixing experiments.


2014 ◽  
Vol 16 (6) ◽  
pp. 1250-1264 ◽  
Author(s):  
Yuchuan Bai ◽  
Xiaolong Song ◽  
Shuxian Gao

Turbulent flow in meandering open channels is one of the most complicated and unpredictable turbulent flows as the interaction of various forces, such as pressure gradient, centrifugal force, and wall shear stresses severely affect the flow pattern. In order to improve significance in engineering application, understanding the overall flow characteristic is the focus. This paper presents the results of numerical and experimental investigations of flow in a 180° mild bend, which is close to criticality with curvature ratio R/B = 3. Considering the characteristic of various models, three-dimensional (3D) re-normalization group (RNG) k–ε model was adopted to simulate the flow efficiently. Governing equations of the flow were solved with a finite-volume method. The pressure-based coupled algorithm was used to compute the pressure. The flow velocities were measured experimentally with Micro acoustic Doppler velocimeter. Good agreement between the numerical results and measurements indicated that RNG k–ε model can successfully predict this flow phenomenon. The flow pattern in this bend is influenced widely by the secondary flow. The variations of velocity components, streamlines, secondary flow, and wall shear stresses are analysed in the study. Some newly discovered phenomenon in this special state are worth noting.


2003 ◽  
Vol 478 ◽  
pp. 227-235 ◽  
Author(s):  
J. D. GIBBON ◽  
Charles R. DOERING

Dissipation-range intermittency was first observed by Batchelor & Townsend (1949) in high Reynolds number turbulent flows. It typically manifests itself in spatio-temporal binary behaviour which is characterized by long, quiescent periods in the signal which are interrupted by short, active ‘events’ during which there are large excursions away from the average. It is shown that Leray's weak solutions of the three-dimensional incompressible Navier–Stokes equations can have this binary character in time. An estimate is given for the widths of the short, active time intervals, which decreases with the Reynolds number. In these ‘bad’ intervals singularities are still possible. However, the average width of a ‘good’ interval, where no singularities are possible, increases with the Reynolds number relative to the average width of a bad interval.


1991 ◽  
Vol 113 (2) ◽  
pp. 400-406 ◽  
Author(s):  
Y. T. Lin ◽  
M. Choi ◽  
R. Greif

A study has been made of the heat transfer, flow, and particle deposition relative to the modified chemical vapor deposition (MCVD) process. The effects of variable properties, buoyancy, and tube rotation have been included in the study. The resulting three-dimensional temperature and velocity fields have been obtained for a range of conditions. The effects of buoyancy result in asymmetric temperature and axial velocity profiles with respect to the tube axis. Variable properties cause significant variations in the axial velocity along the tube and in the secondary flow in the region near the torch. Particle trajectories are shown to be strongly dependent on the tube rotation and are helices for large rotational speeds. The component of secondary flow in the radial direction is compared to the thermophoretic velocity, which is the primary cause of particle deposition in the MCVD process. Over the central portion of the tube the radial component of the secondary flow is most important in determining the motion of the particles.


Sign in / Sign up

Export Citation Format

Share Document