• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3073-3084
• /
• 2022

The propagation of radiation source uncertainties in spent nuclear fuel (SNF) cask shielding calculations is presented in this paper. The uncertainty propagation employs the depletion and source term outputs of the deterministic code STREAM as input to the transport simulation of the Monte Carlo (MC) codes MCS and MCNP6. The uncertainties of dose rate coming from two sources: nuclear data and modeling parameters, are quantified. The nuclear data uncertainties are obtained from the stochastic sampling of the cross-section covariance and perturbed fission product yields. Uncertainties induced by perturbed modeling parameters consider the design parameters and operating conditions. Uncertainties coming from the two sources result in perturbed depleted nuclide inventories and radiation source terms which are then propagated to the dose rate on the cask surface. The uncertainty analysis results show that the neutron and secondary photon dose have uncertainties which are dominated by the cross section and modeling parameters, while the fission yields have relatively insignificant effect. Besides, the primary photon dose is mostly influenced by the fission yield and modeling parameters, while the cross-section data have a relatively negligible effect. Moreover, the neutron, secondary photon, and primary photon dose can have uncertainties up to about 13%, 14%, and 6%, respectively.

• Nuclear Engineering and Technology
• /
• v.54 no.10
• /
• pp.3803-3810
• /
• 2022

OpenMC is a community-driven open-source Monte Carlo neutron and photon transport simulation code. The Weight Window Mesh (WWM) function and an automatic Global Variance Reduction (GVR) method was recently developed and implemented in a developmental branch of OpenMC. This WWM function and GVR method broaden OpenMC's usage in general purposes deep penetration shielding calculations. However, the Local Variance Reduction (LVR) method, which suits the source-detector problem, is still missing in OpenMC. In this work, the Weight Window Generator (WWG) function has been developed and benchmarked for the same branch. This WWG function allows OpenMC to generate the WWM for the source-detector problem on its own. Single-material cases with varying shielding and sources were used to benchmark the WWG function and investigate how to set up the particle histories utilized in WWG-run and WWM-run. Results show that there is a maximum improvement of WWM generated by WWG. Based on the above results, instructions on determining the particle histories utilized in WWG-run and WWM-run for optimal computation efficiency are given and tested with a few multi-material cases. These benchmarks demonstrate the ability of the OpenMC WWG function and the above instructions for the source-detector problem. This developmental branch will be released and merged into the main distribution in the future.

• Radiation Oncology Journal
• /
• v.19 no.1
• /
• pp.66-73
• /
• 2001

Purpose : For the research of Boron Neutron Capture Therapy (BNCT), fast neutrons generated from the MC-50 cyclotron with maximum energy of 34.4 MeV in Korea Cancer Center Hospital were moderated by 70 cm paraffin and then the dose characteristics were investigated. Using these results, we hope to establish the protocol about dose measurement of epi-thermal neutron, to make a basis of dose characteristic of epi-thermal neutron emitted from nuclear reactor, and to find feasibility about accelerator-based BNCT. Method and Materials : For measuring the absorbed dose and dose distribution of fast neutron beams, we used Unidos 10005 (PTW, Germany) electrometer and IC-17 (Far West, USA), IC-18, ElC-1 ion chambers manufactured by A-150 plastic and used IC-l7M ion chamber manufactured by magnesium for gamma dose. There chambers were flushed with tissue equivalent gas and argon gas and then the flow rate was S co per minute. Using Monte Carlo N-Particle (MCNP) code, transport program in mixed field with neutron, photon, electron, two dimensional dose and energy fluence distribution was calculated and there results were compared with measured results. Results : The absorbed dose of fast neutron beams was $6.47\times10^{-3}$ cGy per 1 MU at the 4 cm depth of the water phantom, which is assumed to be effective depth for BNCT. The magnitude of gamma contamination intermingled with fast neutron beams was $65.2{\pm}0.9\%$ at the same depth. In the dose distribution according to the depth of water, the neutron dose decreased linearly and the gamma dose decreased exponentially as the depth was deepened. The factor expressed energy level, $D_{20}/D_{10}$, of the total dose was 0.718. Conclusion : Through the direct measurement using the two ion chambers, which is made different wall materials, and computer calculation of isodose distribution using MCNP simulation method, we have found the dose characteristics of low fluence fast neutron beams. If the power supply and the target material, which generate high voltage and current, will be developed and gamma contamination was reduced by lead or bismuth, we think, it may be possible to accelerator-based BNCT.

• Nuclear Engineering and Technology
• /
• v.50 no.7
• /
• pp.1051-1059
• /
• 2018

Fast periodic pulsed reactor is a type of reactor in which the fission bursts are formed entirely with external reactivity modulation with a specified time periodicity. This type of reactors could generate much larger intensity of neutron beams for experimental use, compared with the steady state reactors. In the design of fast periodic pulsed reactors, the time dependent simulation of the power pulse is majorly based on a point kinetic model, which is known to have limitations. A more accurate calculation method is desired for the design analyses of fast periodic pulsed reactors. Monte Carlo computer code MCNP6 is used for this task due to its three dimensional transport capability with a continuous energy library. Some new routines were added to simulate the rotation of the movable reflector parts in the time dependent calculation. Fast periodic pulsed reactor IBR-2M was utilized to validate the new routines. This reactor is periodically in prompt supercritical state, which lasts for ${\sim}400{\mu}s$, during the equilibrium state. This generates long neutron fission chains, which requires tremendously large amount of computation time during Monte Carlo simulations. Russian Roulette was applied for these very long neutron chains in MCNP6 calculation, combined with other approaches to improve the efficiency of the simulations. In the power pulse of the IBR-2M at equilibrium state, there is some discrepancy between the experimental measurements and the calculated results using the point kinetics model. MCNP6 results matches better the experimental measurements, which shows the merit of using MCNP6 calculation relative to the point kinetics model.

• Nuclear Engineering and Technology
• /
• v.54 no.2
• /
• pp.764-769
• /
• 2022

The United Arab Emirates is currently building and operating four units of the APR-1400 developed by a South Korean vendor, Korea Electric Power Corporation (KEPCO). This paper attempts to perform APR-1400 reactor core analysis by using the well-known two-step method. The two-step method was applied to the APR-1400 first cycle using the open-source nodal diffusion code, KOMODO. In this study, the group constants were generated using CASMO-4 fuel transport lattice code. The simulation was performed in Hot Zero Power (HZP) at steady-state and transient conditions. Some typical parameters necessary for the Nuclear Design Report (NDR) were evaluated in this paper, such as effective neutron multiplication factor, control rod worth, and critical boron concentration for steady-state analysis. Other parameters such as reactivity insertion, power, and fuel temperature changes during the Reactivity Insertion Accident (RIA) simulation were evaluated as well. The results from KOMODO were verified using PARCS and SIMULATE-3 nodal core simulators. It was found that KOMODO gives an excellent agreement.

• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4731-4742
• /
• 2022

RBMK fuel assemblies differ from other LWR FA due to a specific arrangement of the fuel rods, the low enrichment, and the used burnable absorber - erbium. Therefore, there is a challenge to adapt modeling tools, developed for other LWR types, to solve RBMK problems. A set of 10 different depletion simulation schemes were tested to estimate the impact on reactivity and spent fuel composition of possible SCALE code options for the neutron transport modelling and the use of different nuclear data libraries. The simulations were performed using cross-section libraries based on both, VII.0 and VII.1, versions of ENDF/B nuclear data, and assuming continuous energy and multigroup simulation modes, standard and user-defined Dancoff factor values, and employing deterministic and Monte Carlo methods. The criticality analysis with burn-up credit was performed for the SFP loaded with RBMK-1500 FA. Spent fuel compositions were taken from each of 10 performed depletion simulations. The criticality of SFP is found to be overestimated by up to 0.08% in simulation cases using user-defined Dancoff factors comparing the results obtained using the continuous energy library (VII.1 version of ENDF/B nuclear data). It was shown that such discrepancy is determined by the higher U-235 and Pu-239 isotopes concentrations calculated.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1326-1338
• /
• 2017

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of Consortium for Advanced Simulation of Light Water Reactors (CASL) tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated good agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle hot full power equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics/subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of $0.99994{\pm}6.8E-6$ (95% confidence interval). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.688-694
• /
• 2021

The aim of this work is the analysis, design and hardware implementation of the fractional-order point kinetics (FNPK) model along with its closed-loop controller. The stability and closed-loop control of FNPK models are critical issues. The closed-loop stability of the controller-plant structure is established. Further, the designed PI/PD controllers are implemented in real-time on a DSP processor. The simulation and real-time hardware studies confirm that the designed PI/PD controllers result in a damped stable closed-loop response.

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.414-422
• /
• 2021

The aim of this work considers a second order point reactor kinetics model based on the P1 approximation of transport theory, called in this work as P1 point reactor model. The P1 point reactor model implicitly considers the time derivative of the neutron source which has not been thus considered previously. The inverse method to calculate the reactivity in nuclear reactors -chosen because its high accuracy- Matrix Formulation. The numerical results shown that the Matrix Formulation for the reactivity estimation constitutes a method with insignificant calculation errors.

• Nuclear Engineering and Technology
• /
• v.13 no.4
• /
• pp.264-276
• /
• 1981

A simulation procedure which can represent time-dependent nuclear characteristics of TRIGA Mark-III reactor is developed. CITATION, a multi-group diffusion-depletion program, has been utilized as calculational tool. The group structure employed in this study consists of 7 groups: -3-fast and 4-thermal-which is conventionally utilized in TRIGA type reactor analysis. Three-dimensional nuclear characteristics are synthesized by combining results from two-dimensional plane calculation and two-dimensional cylinder calculation, since direct three-dimensional approach is not yet possible. An effort ia made to develope a method which can extract effective zone and group dependent bucklings by neutron diffusion theory rather than conventional zone and/or group independent Ducklings by neutron transport theory, since neutron leakage is quite high for small core such as research reactors. It is turned out that the method developed in this study gives satisfactory results. The calculation is performed under assumptions that all control rods are fully withdrawn, that no samples are inserted in the irradiation holes and that the core is located in the center of the reactor pool. Burnup-dependent variation of core excess reactivity, time dependent change of Xe-135 poisoning and reactivity worth of rotary specimen rack are calculated and compared with operation records. Neutron flux and power distribution as well as neutron spectrum in each irradiation .facility are presented.