• Nuclear Engineering and Technology
• /
• 제41권7호
• /
• pp.893-906
• /
• 2009

The present work presents the development of a Large Eddy Simulation (LES) methodology viable for complex geometries and suitable for the simulation of rod-bundles. The use of LES and Direct Numerical Simulation (DNS) allows for a deeper analysis of the flow field and the use of stochastical tools in order to obtain additional insight into rod-bundle hydrodynamics. Moreover, traditional steady-state CFD simulations fail to accurately predict distributions of velocity and temperature in rod-bundles when the pitch (P) to diameter (D) ratio P/D is smaller than 1.1 for triangular lattices of cylindrical pins. This deficiency is considered to be due to the failure to predict large-scale coherent structures in the region of the gap. The main features of the code include multi-block capability and the use of the fractional step algorithm. As a Sub-Grid-Scale (SGS) model, a Dynamic Smagorinsky model has been used. The code has been tested on plane channel flow and the flow in annular ducts. The results are in excellent agreement with experiments and previous calculations.

• Nuclear Engineering and Technology
• /
• 제55권7호
• /
• pp.2534-2546
• /
• 2023

Gas-cooled space reactor, which adopts He-Xe gas mixture as working fluid, is a better choice for megawatt power generation. In this paper, thermal-hydraulic characteristics of He-Xe gas mixture in 2×2 rod bundle wrapped with helical wires is numerically investigated. The velocity, pressure and temperature distribution of the coolant are obtained and analyzed. The results show that the existence of helical wires forms the vortexes and changes the velocity and temperature distribution. Hot spots are found at the contact corners between helical wires and fuel rods. The highest temperature of the hot spots reach 1600K, while the mainstream temperature is less than 400K. The helical wire structure increases the friction pressure drop by 20%-50%. The effect extent varies with the pitch and the number of helical wires. The helical wire structure leads to the reduction of Nusselt number. Comparing thermal-hydraulic performance ratios (THPR) of different structures, the THPR values are all less than 1. It means that gas-cooled space reactor adopting helical wires could not strengthen the core heat removal performance. This work provides the thermal-hydraulic design basis for He-Xe gas cooled space nuclear reactor.

• Nuclear Engineering and Technology
• /
• 제55권10호
• /
• pp.3775-3786
• /
• 2023

Mixing Vane Grid (MVG) is one of the most important structures in fuel assembly due to its high performance in mixing the coolant and ultimately increasing Critical Heat Flux (CHF), which avoids the temperature rising suddenly of fuel rods. To evaluate the mixing performance of the MVG, a Total Diffusion Coefficient (TDC) mixing coefficient is defined in the subchannel analysis code. Conventionally, the TDC of the spacer grid is obtained from the combination of experiments and subchannel analysis. However, the processing of obtaining and determine a reasonable TDC is much challenging, it is affected by boundary conditions and MVG geometries. In is difficult to perform all the large and costing rod bundle tests. In this paper, the CFD method was applied in TDC analysis. A typical 5 × 5 MVG was simulated and validated to estimate the mixing performance of the MVG. The subchannel code was used to calculate the TDC. Firstly, the CFD method was validated from the aspect of pressure drop and lateral temperature distribution in the subchannels. Then the effect of boundary conditions including the inlet temperature, inlet velocities, heat flux ratio between hot and cold rods and the arrangement of hot and cold rods on MVG mixing and TDC were studied. The geometric effects on mixing are also carried out in this paper. The effect of vane pattern on mixing was investigated to determine which one is the best to represent the grid's mixing performance.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2728-2735
• /
• 2022

The research on the flow and heat transfer characteristics of lead bismuth(LBE) is significant for the thermal-hydraulic calculation, safety analysis and practical application of lead-based fast reactors(LFR). In this paper, a new CFD model is proposed to solve the thermal-hydraulic analysis of LBE. The model includes two parts: turbulent model and turbulent Prandtl, which are the important factors for LBE. In order to find the best model, the experiment data and design of 19-pin hexagonal rod bundle with spacer grid, undertaken at the Karlsruhe Liquid Metal Laboratory (KALLA) are used for CFD calculation. Furthermore, the turbulent model includes SST k - 𝜔 and k - 𝜀; the turbulent Prandtl includes Cheng-Tak and constant (Pr_t =1.5,2.0,2.5,3.0). Among them, the combination between SST k - 𝜔 and Cheng-Tak is more suitable for the experiment. But in the low Pe region, the deviation between the experiment data and CFD result is too much. The reason may be the inlet-effect and when Pe is in a low level, the number of molecular thermal diffusion occupies an absolute advantage, and the buoyancy will enhance. In order to test and verify versatility of the model, the NCCL performed by the Nuclear Thermal-hydraulic Laboratory (Nuthel) of Xi'an Jiao tong University is used for CFD to calculate. This paper provides two verification examples for the new universal model.

• Nuclear Engineering and Technology
• /
• 제41권5호
• /
• pp.677-690
• /
• 2009

A previous study demonstrated that the two-row fuel assembly has much more favorable neutron-physical and thermal-hydraulic behavior than the conventional one-row fuel assemblies. Based on the newly developed two-row fuel assembly, an SCWR core is proposed and analyzed. The performance of the proposed core is investigated with 3-D coupled neutron-physical and thermal-hydraulic calculations. During the coupling procedure, the thermal-hydraulic behavior is analyzed using a sub-channel analysis code and the neutron-physical performance is computed with a 3-D diffusion code. This paper presents the main results achieved thus far related to the distribution of some neutronic and thermal-hydraulic parameters. It shows that with adjustment of the coolant and moderator mass flow in different assemblies, promising neutron-physical and thermal-hydraulic behavior of the SCWR core is achieved. A sensitivity study of the heat transfer correlation is also performed. Since the pin power in fuel assemblies can be non-uniform, a sub-channel analysis is necessary in order to investigate the detailed distribution of thermal-hydraulic parameters in the hottest fuel assembly. The sub-channel analysis is performed based on the bundle averaged parameters obtained with the core analysis. With the sub-channel analysis approach, more precise evaluation of the hot channel factor and maximum cladding surface temperature can be achieved. The difference in the results obtained with both the sub-channel analysis and the fuel assembly homogenized method confirms the importance of the sub-channel analysis.

• Nuclear Engineering and Technology
• /
• 제49권1호
• /
• pp.267-276
• /
• 2017

In any reactor physics analysis, the instantaneous power distribution in the core can be calculated when the actual bundle-wise burnup distribution is known. Considering the fact that CANDU (Canada Deuterium Uranium) utilizes on-power refueling to compensate for the reduction of reactivity due to fuel burnup, in the CANDU fuel management analysis, snapshots of power and burnup distributions can be obtained by simulating and tracking the reactor operation over an extended period using various tools such as the $^*SIMULATE$ module of the Reactor Fueling Simulation Program (RFSP) code. However, for some studies, such as an evaluation of a conceptual design of a next-generation CANDU reactor, the preferred approach to obtain a snapshot of the power distribution in the core is based on the patterned-channel-age model implemented in the $^*INSTANTAN$ module of the RFSP code. The objective of this approach is to obtain a representative snapshot of core conditions quickly. At present, such patterns could be generated by using a program called RANDIS, which is implemented within the $^*INSTANTAN$ module. In this work, we present an alternative approach to derive the patterned-channel-age model where a simulated-annealing-based algorithm is used to find such patterns, which produce reasonable power distributions.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1532-1539
• /
• 2019

This paper analyzes a new method for reducing boiling water reactor fuel temperature during a Loss of Coolant Accident (LOCA). The method uses a device called Reverse Flow Restriction Device (RFRD) at the inlet of fuel bundles in the core to prevent coolant loss from the bundle inlet due to the reverse flow after a large break in the recirculation loop. The device allows for flow in the forward direction which occurs during normal operation, while after the break, the RFRD device changes its status to prevent reverse flow. In this paper, a detailed simulation of LOCA has been carried out using the U.S. NRC's TRACE code to investigate the effect of RFRD on the flow rate as well as peak clad temperature of BWR fuel bundles during three different LOCA scenarios: small break LOCA (25% LOCA), large break LOCA (100% LOCA), and double-ended guillotine break (200% LOCA). The results demonstrated that the device could substantially block flow reversal in fuel bundles during LOCA, allowing for coolant to remain in the core during the coolant blowdown phase. The device can retain additional cooling water after activating the emergency systems, which maintains the peak clad temperature at lower levels. Moreover, the RFRD achieved the reflood phase (when the saturation temperature of the clad is restored) earlier than without the RFRD.

• Nuclear Engineering and Technology
• /
• 제51권1호
• /
• pp.310-316
• /
• 2019

Fluidelastic instability is a key issue in steam generator tube bundles subjected in cross-flow. With a low flow velocity, a large amplitude vibration of the tube observed by many researchers. However, the mechanism of this vibration is seldom analyzed. In this paper, the mechanism of cross-flow effects on fluidelastic instability of tube bundles was investigated. Analysis reveals that when the system reaches the critical state, there would be two forms, with two critical velocities, and thus two expressions for the critical velocities were obtained. Fluidelastic instability experiment and numerical analysis were conducted to obtain the critical velocity. And, if system damping is small, with increases of the flow velocity, the stability behavior of tube array changes. At a certain flow velocity, the stability of tube array reaches the first critical state, a dynamic bifurcation occurs. The tube array returns to a stable state with continues to increase the flow velocity. At another certain flow velocity, the stability of tube array reaches the second critical state, another dynamic bifurcation occurs. However, if system damping is big, there is only one critical state with increases the flow velocity. Compared the results of experiments to numerical analysis, it shows a good agreement.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3166-3175
• /
• 2022

The most numerical investigations of the thermal-hydraulic phenomena following the loss of the residual heat removal capability during the mid-loop operation of the pressurized water reactor were performed according to simplifications and are not sufficiently accurate. To perform more accurate and more reliable predictions of thermal-hydraulic accidents in a nuclear power plant using computational fluid dynamics codes, a more detailed methodology is needed. Modelling results identified that thermal stratification and natural convection are observed. Temperatures of lower monitoring points remain low, while temperatures of upper monitoring points increase over time. The water in the heated region, in the upper unheated region and the pipe region was well mixed due to natural convection, meanwhile, there is no natural convection in the lower unheated region. Water temperature in the pipe region increased after a certain time delay due to circulation of flow induced by natural convection in the heated and upper unheated regions. The modelling results correspond to the experimental data. The developed computational fluid dynamics methodology could be applied for modelling of two-component single/two-phase natural convection and thermal stratification phenomena during the mid-loop operation of the pressurized water reactor or other nuclear and non-nuclear installations at similar conditions.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2714-2719
• /
• 2022

The present paper proposes some new computational methods and results in the framework of flow computation through congested domains seen as porous media, as it can be found in the core of a Pressurized Water Reactor (PWR). The flow is thus mostly governed by the distribution of pressure losses, both through the porous structures, such as fuel assemblies, and in the thin fluid layers between them. The purpose of the present paper is to consider the question of the interaction of a flow and a rod bundle from an analytical point of view gathering all the contributions through a set of equations as simple and representative as possible. It aims at demonstrating a sound understanding of the relevant phenomena governing the flow establishment in the geometry of interest instead of relying mainly on a posteriori observations obtained both experimentally and numerically. Comparison with two set of experimental results showed good agreement. The model proposed being analytical it appears easily implementable for studies needing an expression of fluid forces in a rod array as for fuel assembly bowing issue. It would be interesting to test the reliability of the model on other geometry with different P/R ratios.