scholarly journals A simulation of a pebble bed reactor core by the MCNP-4C computer code

2009 ◽  
Vol 24 (3) ◽  
pp. 177-182 ◽  
Author(s):  
Moshkbar Bakhshayesh ◽  
Naser Vosoughi

Lack of energy is a major crisis of our century; the irregular increase of fossil fuel costs has forced us to search for novel, cheaper, and safer sources of energy. Pebble bed reactors - an advanced new generation of reactors with specific advantages in safety and cost - might turn out to be the desired candidate for the role. The calculation of the critical height of a pebble bed reactor at room temperature, while using the MCNP-4C computer code, is the main goal of this paper. In order to reduce the MCNP computing time compared to the previously proposed schemes, we have devised a new simulation scheme. Different arrangements of kernels in fuel pebble simulations were investigated and the best arrangement to decrease the MCNP execution time (while keeping the accuracy of the results), chosen. The neutron flux distribution and control rods worth, as well as their shadowing effects, have also been considered in this paper. All calculations done for the HTR-10 reactor core are in good agreement with experimental results.

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.


2017 ◽  
Vol 2017 ◽  
pp. 1-6
Author(s):  
Jingyu Zhang ◽  
Fu Li ◽  
Yuliang Sun

The pebble-bed reactor HTR-PM is being built in China and is planned to be critical in one or two years. At present, one emphasis of engineering design is to determine the fuel management scheme of the initial core and running-in phase. There are many possible schemes, and many factors need to be considered in the process of scheme evaluation and analysis. Based on the experience from the constructed or designed pebble-bed reactors, the fuel enrichment and the ratio of fuel spheres to graphite spheres are important. In this paper, some relevant physical considerations of the initial core and running-in phase of HTR-PM are given. Then a typical scheme of the initial core and running-in phase is proposed and simulated with VSOP code, and some key physical parameters, such as the maximum power per fuel sphere, the maximum fuel temperature, the refueling rate, and the discharge burnup, are calculated. Results of the physical parameters all satisfy the relevant design requirements, which means the proposed scheme is safe and reliable and can provide support for the fuel management of HTR-PM in the future.


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.


Author(s):  
Geoffrey J. Peter

High Temperature Gas Cooled Reactor (HTGR) development and operation is expanding in the United Kingdom, Russia, USA (Generation IV Reactors), and France (Pebble Bed Modular Reactor, PBMR). A prototype pebble bed reactor producing 10 MW thermal, High Temperature Reactor (HTR-10) is in operation in China. However, the general public remains skeptical of the safety and the perceived dangers of possible accidents. Of particular concern are blockages caused by local variations in flow and heat transfer that lead to hot spots within the bed. This paper models the accident scenario resulting from blockages due to the retention of dust in the coolant gas or from the rupture of one or more fuel particles used in the High Temperature Gas Cooled (Pebble Bed) Nuclear Reactors using the commercially available computer code COMSOL. Numerical modeling of flow and heat transfer in a packed bed produces an Elliptical Non-Linear Partial Differential equation that requires custom made computer codes. Previously published results obtained from the use of a custom-made verified computer code are limited to one accident scenario and involve considerable modification to study different accident scenarios. Thus the use of a commercially available computer code that can simulate many different accident scenarios is of considerable advantage. Further, this paper compares numerical solutions obtained from custom-made computer code with COMSOL simulation and discusses the advantages and limitations of both codes.


Author(s):  
Hery Adrial ◽  
Amir Hamzah ◽  
Entin Hartini

GAMMA DOSE RATE ANALYSIS IN BIOLOGICAL SHIELDING OF HTGR-10 MWth PEBBLE BED REACTOR. HTGR-10 MWth is a high-temperature gas-cooled reactor. The fuel and moderator are pebble shaped with a radius of 3 cm. One fuel pebble consists of thousands of UO2 kernels with a density of 10.4 gram/cc and the enrichment rate of 17%. The core of HTGR-10 MWth is the center of origin of neutrons and gamma radiation resulting from the interaction of neutrons with pebble fuel, moderator and biological shield. The various types of radiations generated from such nuclear reactions should be monitored to ensure the safety of radiation workers. This research was conducted using MCNP-6 Program package with the aim to calculate and analyze gamma radiation dose in biological shield of HTGR-10 MWth. In this study, the biological shield is divided into 10 equal segments. The first step of the research is to benchmark the created program against the critical height of HTR-10. The results of the benchmarking show an error rate of ± 1.1327%, while the critical core height of HTGR 10 MWth for the ratio of pebble fuel and pebble moderator (F:M) of 52: 48 occurs at a height of 134 cm. The rate of gamma dose at the core is 3.0052E + 05 mSv/hr. On the biological shield made of regular concrete with a density of 2.3 grams/cc, the rate of gamma dose decreases according to an equation y = 0.0042 e-0.03x. Referring to Perka Bapeten no 4 of 2013, the safe limits for workers and radiation protection officers will be achieved if the minimum thickness of biological shield is 115 cm with gamma dose rate of 0 mSv/hour.Keywords: Gamma dose rate, HTGR 10 MWth, biological shield, pebble


Author(s):  
Ramatsemela Mphahlele ◽  
Abderrafi M. Ougouag ◽  
Kostadin N. Ivanov ◽  
Hans D. Gougar

A practical methodology is developed for the determination of spectral zones in Pebble Bed Reactors (PBR). The methodology involves the use of spectral indices based on few-group diffusion theory whole core calculations. In this work a spectral zone is defined as made up of a number of nodes whose characteristics are collectively similar and that are assigned the same few-group diffusion constants. Therefore, spectral indices that reflect the physical behaviors of interest can be used to characterize said behaviors within each zone and thus to identify and distinguish the spectral zones. Several plausible spectral indices have been investigated in this work. Special emphasis and focus was placed on the trend or behavior of the spectral index at various positions along the radial and axial dimensions in the core. The ratio of group-wise surface currents to total surface fluxes, has been used to successfully identify spectral zone boundaries. A plot of the absolute value of this ratio versus position in the reactor exhibits a series of minima and maxima points. These extrema correlate with regions of significant spectral changes, and therefore are identified as plausible spectral zone boundaries.


Sign in / Sign up

Export Citation Format

Share Document