ICONE19-43215 NUMERICAL RESEARCH ON THE THERMAL HYDRAULICS OF THE COOLANT IN A PEBBLE BED REACTOR CORE BY CFD

2011 ◽  
Vol 2011.19 (0) ◽  
pp. _ICONE1943-_ICONE1943
Author(s):  
Hua Li ◽  
Suizheng Qiu ◽  
Guanghui Su ◽  
Wenxi Tian ◽  
Youjia Zhang
Author(s):  
Xiang Zhao ◽  
Trent Montgomery ◽  
Sijun Zhang

This paper presents combined computational fluid dynamics (CFD) and discrete element method (DEM) simulations of fluid flow and relevant heat transfer in the pebble bed reactor core. In the pebble bed reactor core, the coolant passes highly complicated flow channels, which are formed by thousands of pebbles in a random way. The random packing structure of pebbles is crucial to CFD simulations results. The realistic packing structure in an entire pebble bed reactor (PBR) is generated by discrete element method (DEM). While in CFD calculations, selection of the turbulence models have great importance in accuracy and capturing the details of the flow features, in our numerical simulations both large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) models are employed to investigate the effects of different turbulence models on gas flow field and relevant heat transfer. The calculations indicate the complex flow structure within the voids between the pebbles.


2021 ◽  
Vol 927 (1) ◽  
pp. 012018
Author(s):  
Nicholas Sidharta ◽  
Almanzo Arjuna

Abstract Pebble bed reactor with a once-through-then-out fuelling scheme has the advantage of simplifying the refueling system. However, the core upper-level power density is relatively higher than the bottom, producing an asymmetric core axial power distribution. Several burnable poison (BP) configurations are used to flatten the peak power density and improve power distribution while suppressing the excess core reactivity at the beginning of the burnup cycle. This study uses HTR-PM, China’s pebble bed reactor core, to simulate several burnable poison (BP) configurations. Serpent 2 coupled with Octave and a discrete element method simulation is used to model and simulate the pebble bed reactor core. It is found that erbium needs a large volumetric fraction in either QUADRISO or distributed BP to perform well. On the other hand, gadolinium and boron need a smaller volumetric fraction but perform worse in radial power distribution criteria in the fuel sphere. This study aims to verify the effect of BP added fuel pebbles on an OTTO refueling scheme HTR-PM core axial power distribution and excess reactivity.


2014 ◽  
Vol 2014 ◽  
pp. 1-16 ◽  
Author(s):  
Shixiong Song ◽  
Xiangzhou Cai ◽  
Yafen Liu ◽  
Quan Wei ◽  
Wei Guo

The present paper systematically investigated pore scale thermal hydraulics characteristics of molten salt cooled high temperature pebble bed reactor. By using computational fluid dynamics (CFD) methods and employing simplified body center cubic (BCC) and face center cubic (FCC) model, pressure drop and local mean Nusselt number are calculated. The simulation result shows that the high Prandtl number molten salt in packed bed has unique fluid-dynamics and thermodynamic properties. There are divergences between CFD results and empirical correlations’ predictions of pressure drop and local Nusselt numbers. Local pebble surface temperature distributions in several default conditions are investigated. Thermal removal capacities of molten salt are confirmed in the case of nominal condition; the pebble surface temperature under the condition of local power distortion shows the tolerance of pebble in extreme neutron dose exposure. The numerical experiments of local pebble insufficient cooling indicate that in the molten salt cooled pebble bed reactor, the pebble surface temperature is not very sensitive to loss of partial coolant. The methods and results of this paper would be useful for optimum designs and safety analysis of molten salt cooled pebble bed reactors.


Nukleonika ◽  
2019 ◽  
Vol 64 (4) ◽  
pp. 131-138
Author(s):  
◽  
Topan Setiadipura ◽  
Jim C. Kuijper ◽  

Abstract As a crucial core physics parameter, the control rod reactivity has to be predicted for the control and safety of the reactor. This paper studies the control rod reactivity calculation of the pebble-bed reactor with three scenarios of UO2, (Th,U)O2, and PuO2 fuel type without any modifications in the configuration of the reactor core. The reactor geometry of HTR-10 was selected for the reactor model. The entire calculation of control rod reactivity was done using the MCNP6 code with ENDF/B-VII library. The calculation results show that the total reactivity worth of control rods in UO2-, (U,Th)O2-, and PuO2-fueled cores is 15.87, 15.25, and 14.33%Δk/k, respectively. These results prove that the effectiveness of total control rod in thorium and uranium cores is almost similar to but higher than that in plutonium cores. The highest reactivity worth of individual control rod in uranium, thorium and plutonium cores is 1.64, 1.44, and 1.53%Δk/k corresponding to CR8, CR1, and CR5, respectively. The other results demonstrate that the reactor can be safely shutdown with the control rods combination of CR3+CR5+CR8+CR10, CR2+CR3+CR7+CR8, and CR1+CR3+CR6+CR8 in UO2-, (U,Th)O2-, and PuO2-fueled cores, respectively. It can be concluded that, even though the calculation results are not so much different, however, the selection of control rods should be considered in the pebble-bed core design with different scenarios of fuel type.


2014 ◽  
Vol 2014 ◽  
pp. 1-10 ◽  
Author(s):  
Xinlong Jia ◽  
Xingtuan Yang ◽  
Nan Gui ◽  
Yu Li ◽  
Jiyuan Tu ◽  
...  

The experimental method (side area method) and DEM simulation have been carried out to analyse the stagnant zone in the quasi-two-dimensional silos. The side area method is a phenomenological method by means of investigating the interface features of different areas composed of different coloured pebbles. Two methods have been discussed to define the stagnant zone. In particular, the area of the stagnant zone has been calculated with the mean-streamline method, and the tracking time of different marking pebbles has been investigated with the stagnant time method to explore the kinematics characteristics of the pebbles. The stagnant zone is crucial for the safety of the pebble-bed reactor, and the practical reactor core must avoid the existence of the stagnant zone. Furthermore, this paper also analyses the effects of bed configuration (the bed height, the base angle, and the friction coefficient) on stagnant zone with the two methods mentioned above. In detail, the bed height shows little impact on the stagnant zones when the bed height exceeds a certain limit, while the base angle has negative prominent correlation with the stagnant zone. The friction coefficient effect seems complicated and presents the great nonlinearity, which deserves to be deeply investigated.


Author(s):  
Gang Zhao ◽  
Ping Ye ◽  
Toru Obala

Spherical fuel elements are distributed randomly in the pebble bed reactor core and helium flow through the pebble bed to remove nuclear reaction heat. Pebble bed reactor core is usually treated as a uniform porous media flow in thermal hydraulic research. However the porosity distribution is nonuniform and the porosity near the wall increase sharply. A new random model is developed in this paper to investigate thermal hydraulic characteristics of pebble bed reactor core. Porosity assumption is based on porosity measurement of other research. Porosity simulation is divided into three parts according to the distance from wall. In the center of core, porosity is assumed to obey normal distribution, where average porosity is from the experimental relation based on statistical results. The mean and standard deviation of porosity distribution near the wall will increase because of the wall effect, where the distance from the wall is about three times of fuel ball’s diameter. The third part is zone from three times to five times of ball’s diameter departed from the wall. The wall effect of this zone is between center and the wall zone. Based on above assumption, a random porosity simulation is completed to apply in this research. COMSOL Multiphysics 3.5a software is used in this research. COMSOL Multiphysics are a calculation platform using proven Finite Elements Methods (FEM). In this research, Brinkman equation for porous media flow is applied in the simulation. Non-thermal Balance model is used in local heat transfer research between gas and pebble bed. A geometry model is built to simulate HTR-10. Temperature profile of variant porosity is gained from stationary analysis and comparison with uniform porosity is also discussed in the paper. For transient analysis, four cases simulation is carried out in the research. Case 1 and 2 simulate heat transfer phenomena with forced cooling system and with passive cooling system after reactor shut down. Way-Wigner-curve is used in Case 1 and Case 2 to simulate decay heat in the calculation. Case 3 and Case 4 simulate ATWS phenomena with natural convection and without natural convection system when blower is trip off in normal operation. Simulation results also are compared with some ATWS experiments and some discussion is done in the paper. From the results, it can be seen that random porosity will affect temperature distribution near the wall and make outlet temperature non-uniform greatly. The maximum temperature of variant porosity is much greater than the maximum temperature of uniform porosity at the same condition. Transient analyses of variant model show that passive cooling system can remove residual heat even in accident conditions when the blower trip off whether reactor shut down or not and the analyses results correspond substantially with experimental results. In general, variant porosity should be considered in the thermal hydraulic research of pebble bed reactor core. Variant porosity model can provide good prediction of heat transfer phenomena than uniform porosity model. Especially it can explain some transient analysis results.


Sign in / Sign up

Export Citation Format

Share Document