• Nuclear Engineering and Technology
• /
• 제42권5호
• /
• pp.520-545
• /
• 2010

This paper presents the neutronics modeling capabilities of the fast reactor simulation system SHARP, which ANL is developing as part of the U.S. DOE's NEAMS program. We discuss the three transport solvers (PN2ND, SN2ND, and MOCFE) implemented in the UNIC code along with the multigroup cross section generation code $MC^2$-3. We describe the solution methods and modeling capabilities, and discuss the improvement needs for each solver, focusing on massively parallel computation. We present the performance test results against various benchmark problems and ZPR-6 and ZPPR critical experiments. We also discuss weak and strong scalability results for the SN2ND solver on the ZPR-6 critical assembly benchmarks.

• Nuclear Engineering and Technology
• /
• 제48권3호
• /
• pp.785-794
• /
• 2016

The Argonne National Laboratory of the United States and the Kharkov Institute of Physics and Technology of the Ukraine have been collaborating on the design, development and construction of a neutron source facility at Kharkov Institute of Physics and Technology utilizing an electron-accelerator-driven subcritical assembly. The electron beam power is 100 kW using 100-MeV electrons. The facility was designed to perform basic and applied nuclear research, produce medical isotopes, and train nuclear specialists. The biological shield of the accelerator building was designed to reduce the biological dose to less than 5.0e-03 mSv/h during operation. The main source of the biological dose for the accelerator building is the photons and neutrons generated from different interactions of leaked electrons from the electron gun and the accelerator sections with the surrounding components and materials. The Monte Carlo N-particle extended code (MCNPX) was used for the shielding calculations because of its capability to perform electron-, photon-, and neutron-coupled transport simulations. The photon dose was tallied using the MCNPX calculation, starting with the leaked electrons. However, it is difficult to accurately tally the neutron dose directly from the leaked electrons. The neutron yield per electron from the interactions with the surrounding components is very small, ~0.01 neutron for 100-MeV electron and even smaller for lower-energy electrons. This causes difficulties for the Monte Carlo analyses and consumes tremendous computation resources for tallying the neutron dose outside the shield boundary with an acceptable accuracy. To avoid these difficulties, the SOURCE and TALLYX user subroutines of MCNPX were utilized for this study. The generated neutrons were banked, together with all related parameters, for a subsequent MCNPX calculation to obtain the neutron dose. The weight windows variance reduction technique was also utilized for both neutron and photon dose calculations. Two shielding materials, heavy concrete and ordinary concrete, were considered for the shield design. The main goal is to maintain the total dose outside the shield boundary less than 5.0e-03 mSv/h during operation. The shield configuration and parameters of the accelerator building were determined and are presented in this paper.

• Journal of Radiation Protection and Research
• /
• 제41권4호
• /
• pp.344-349
• /
• 2016

Background: Carbon ion therapy has achieved satisfactory results. However, patients have a risk to get a secondary cancer. In order to estimate the risk, it is essential to understand particle transportation and nuclear reactions in the patient's body. The particle transport Monte Carlo simulation code is a useful tool to understand them. Since the code validation for heavy ion incident reactions is not enough, the experimental data of the elementary reaction processes are needed. Materials and Methods: We measured neutron production double-differential cross-sections (DDXs) on a carbon bombarded with 430 MeV/nucleon carbon beam at PH2 beam line of HIMAC facility in NIRS. Neutrons produced in the target were measured with NE213 liquid organic scintillators located at six angles of 15, 30, 45, 60, 75, and $90^{\circ}$. Results and Discussion: Neutron production double-differential cross-sections for carbon bombarded with 430 MeV/nucleon carbon ions were measured by the time-of-flight method with NE213 liquid organic scintillators at six angles of 15, 30, 45, 60, 75, and $90^{\circ}$. The cross sections were obtained from 1 MeV to several hundred MeV. The experimental data were compared with calculated results obtained by Monte Carlo simulation codes PHITS, Geant4, and FLUKA. Conclusion: PHITS was able to reproduce neutron production for elementary processes of carbon-carbon reaction precisely the best of three codes.

• Nuclear Engineering and Technology
• /
• 제55권4호
• /
• pp.1350-1364
• /
• 2023

2D/1D coupling method is an important neutron transport calculation method due to its high accuracy and relatively low computation cost. However, 2D/1D coupling method may diverge especially in small axial mesh size. To analyze the convergence behavior of 2D/1D coupling method, a Fourier analysis for k-eigenvalue neutron transport problems is implemented. The analysis results present the divergence problem of 2D/1D coupling method in small axial mesh size. Several common attempts are made to solve the divergence problem, which are to increase the number of inner iterations of the 2D or 1D calculation, and two times 1D calculations per outer iteration. However, these attempts only could improve the convergence rate but cannot deal with the divergence problem of 2D/1D coupling method thoroughly. Moreover, the choice of axial solvers, such as DGFEM SN and traditional SN, and its effect on the convergence behavior are also discussed. The results show that the choice of axial solver is a key point for the convergence of 2D/1D method. The DGFEM SN based 2D/1D method could converge within a wide range of optical thickness region, which is superior to that of traditional SN method.

• Nuclear Engineering and Technology
• /
• 제52권8호
• /
• pp.1817-1825
• /
• 2020

Background: Cosmic ray-induced particles can lead to failure of semiconductors packaged for export during air transport. This work performed MCNP 6.2 simulations to optimize shielding against neutrons and protons induced by cosmic radiation Methods and materials: The energy spectra of protons and neutrons by incident angle at the flight altitude were determined using atmospheric cuboid model. Various candidates for the shielding materials and the geometry of the Unit Load Device Container were evaluated to determine the conditions that allow optimal shielding at all sides of the container. Results: It was found that neutrons and protons, at the flight altitude, generally travel with a downward trajectory especially for the particles with high energy. This indicated that the largest number of particles struck the top of the container. Furthermore, the simulation results showed that, among the materials tested, borated polyethylene and stainless steel were the most optimal shielding materials. The optimal shielding structure was also determined with the weight limit of the container in consideration. Conclusions: Under the determined optimal shielding conditions, a significantly reduced number of neutrons and protons reach the contents inside the container, which ultimately reduces the possibility of semiconductor failure during air transport.

• Nuclear Engineering and Technology
• /
• 제55권10호
• /
• pp.3732-3753
• /
• 2023

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the k_eff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.

• 방사성폐기물학회지
• /
• 제16권3호
• /
• pp.359-367
• /
• 2018

Long-lived Small Modular Reactors are being promoted as an innovative way of catering to emerging markets and isolated regions. They can be operated continuously for decades without requiring additional fuel. A novel configuration of long-lived reactor core employs a mixed neutron spectrum, providing an improvement in nonproliferation metrics and in safety characteristics. Starting with a base sodium reactor design, moderating material is inserted in outer core assemblies to modify the fast spectrum. The assemblies are shuffled once during core lifetime to ensure that every fuel rod is exposed to the thermalized spectrum. The Mixed Spectrum Reactor is able to maintain a core lifetime over two decades while ensuring the plutonium it breeds is below the weapon-grade limit at the fuel discharge. The main drawbacks of the design are higher front-end fuel cycle costs and a 58% increase in core volume, although it is alleviated to some extent by a 48% higher power output.

• Nuclear Engineering and Technology
• /
• 제55권5호
• /
• pp.1593-1603
• /
• 2023

In deterministic and Monte Carlo transport codes, b-delayed emission is included using a group structure where all of the precursors are grouped together in 6 groups or families, but given the increase in computational power, nowadays there is no reason to keep this structure. Furthermore, there have been recent efforts to compile and evaluate all the available b-delayed neutron emission data and to measure new and improved data on individual precursors. In order to be able to perform a transient Monte Carlo simulation, data from individual precursors needs to be implemented in a transport code. This work is the first step towards the development of a tool to explore the effect of individual precursors in a fissile system. In concrete, individual precursor data is included by expanding the capabilities of the open source Monte Carlo code OpenMC. In the modified code - named Time Dependent OpenMC or OpenMC(TD)- time dependency related to β-delayed neutron emission was handled by using forced decay of precursors and combing of the particle population. The data for continuous energy neutron cross-sections was taken from JEFF-3.1.1 library. Regarding the data needed to include the individual precursors, cumulative yields were taken from JEFF-3.1.1 and delayed neutron emission probabilities and delayed neutron spectra were taken from ENDF-B/VIII.0. OpenMC(TD) was tested in a monoenergetic system, an energy dependent unmoderated system where the precursors were taken individually or in a group structure, and in a light-water moderated energy dependent system, using 6-groups, 50 and 40 individual precursors. Neutron flux as a function of time was obtained for each of the systems studied. These results show the potential of OpenMC(TD) as a tool to study the impact of individual precursor data on fissile systems, thus motivating further research to simulate more complex fissile systems.

• 한국환경과학회지
• /
• 제13권2호
• /
• pp.175-180
• /
• 2004

A prompt gamma-ray neutron activation (PGNA) system was simulated by the Monte Carlo N-Particle transport code (MCNP-4A) to estimate the level at which the scattered photon fluence rate, the absolute efficiency of the HPGe-detector, the volume of the concrete sample and the $^{35}$ /Cl(n, ${\gamma}$) reaction rate in this sample contribute to the count rate in the NaCl concentration measurement. The n- ${\gamma}$ fluence rates at the ST-2 beam tube exit of the HANARO reactor were used as input data, and the GAMMA-X type HPGe detector was modeled to tally 1.1649 MeV ${\gamma}$ -rays emitted from the $^{35}$ Cl(n, ${\gamma}$) reaction in the concrete sample. For three cylindrical concrete samples of 13.8, 46.8 and 157.1 ㎤ volumes, respectively, the relations between the NaCl weight fractions of 0.1, 1, 2 and 5 % in each of the concrete samples and the 1.1 649 MeV pulses created in the HPGe detector model were studied. As a result, it was found that the count rate at the same NaCl concentration nearly depends on the volume of the samples in a simulated condition of the same NaCl concentration samples, and that the linearities of the NaCl concentration calibration curves were reasonable in the narrow range of the NaCl weight fraction.

• Nuclear Engineering and Technology
• /
• 제44권2호
• /
• pp.161-176
• /
• 2012

McCARD is a Monte Carlo (MC) neutron-photon transport simulation code. It has been developed exclusively for the neutronics design of nuclear reactors and fuel systems. It is capable of performing the whole-core neutronics calculations, the reactor fuel burnup analysis, the few group diffusion theory constant generation, sensitivity and uncertainty (S/U) analysis, and uncertainty propagation analysis. It has some special features such as the anterior convergence diagnostics, real variance estimation, neutronics analysis with temperature feedback, $B_1$ theory-augmented few group constants generation, kinetics parameter generation and MC S/U analysis based on the use of adjoint flux. This paper describes the theoretical basis of these features and validation calculations for both neutronics benchmark problems and commercial PWR reactors in operation.