• 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.

• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1097-1104
• /
• 2023

A preliminary criticality analysis for novel pyrochemical apparatuses for the reprocessing of mixed uranium-plutonium nitride spent nuclear fuel from the BREST-OD-300 reactor was performed. High-temperature processing apparatuses, "metallization" electrolyzer, refinery remelting apparatus, refining electrolyzer, and "soft" chlorination apparatus are considered in this work. Computational models of apparatuses for two neutron radiation transport codes (MCU-FR and MCNP) were developed and calculations for criticality were completed using the Monte Carlo method. The criticality analysis was performed for different loads of fissile material into the apparatuses including overloading conditions. Various emergency situations were considered, in particular, those associated with water ingress into the chamber of the refinery remelting apparatus. It was revealed that for all the considered computational models nuclear safety rules are satisfied.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.19 no.1
• /
• pp.27-35
• /
• 2023

Pressure tubes are the main components of PHWR core and serve as the pressure boundary of the primary heat transport system. However, because pressure tubes have changed their geometrical dimensions under the severe operating conditions of high temperature, high pressure and neutron irradiation according to the increase of operation time, all dimensional changes should be predicted to ensure that dimensions remain within the allowable design ranges during the operation. Among the deformations, the diameter expansion due to creep leads to the increase of bypass flow which may not contribute to the fuel cooling, the decrease of critical channel power and finally the deration of the power to maintain the operational safety margin. This study is focused on the modeling of the expansion of the pressure tube diameter based on the operating conditions and measured diameter data. The pressure tube diameter expansion was modeled using the neutron flux and temperature distributions of each fuel channel and each fuel bundle as well as the measured diameter data. Although the basic concept of the current modeling approach is simple, the diameter prediction results using the developed methodology showed very good agreement with the real data, compared to the existing methodology.

• Nuclear Engineering and Technology
• /
• v.50 no.6
• /
• pp.820-828
• /
• 2018

The growing nuclear threat has amplified the need for developing diverse and accurate nuclear forensics analysis techniques to strengthen nuclear security measures. The work presented here is part of a research effort focused on developing a methodology for reactor-type discrimination of weapons-grade plutonium. To verify the developed methodology, natural $UO_2$ fuel samples were irradiated in a thermal neutron spectrum at the University of Missouri Research Reactor (MURR) and produced approximately $20{\mu}g$ of weapons-grade plutonium test material. Radiation transport simulations of common thermal reactor types that can produce weapons-grade plutonium were performed, and the results are presented here. These simulations were needed to verify whether the plutonium produced in the natural $UO_2$ fuel samples during the experimental irradiation at MURR was a suitable representative to plutonium produced in common thermal reactor types. Also presented are comparisons of fission product and plutonium concentrations obtained from computational simulations of the experimental irradiation at MURR to the nondestructive and destructive measurements of the irradiated natural $UO_2$ fuel samples. Gamma spectroscopy measurements of radioactive fission products were mostly within 10%, mass spectroscopy measurements of the total plutonium mass were within 4%, and mass spectroscopy measurements of stable fission products were mostly within 5%.

• 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.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.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.230-237
• /
• 2020

In this part, an implicit time dependent solution is presented for the Boltzmann transport equation discretized by the analytic coarse mesh finite difference method (ACMFD) over the spatial domain as well as the simplified P₃ (SP₃) for the angular variable. In the first part of this work we proposed a SP₃-ACMFD approach to solve the static eigenvalue equations which provide the initial conditions for temp dependent equations. Having solved the 3D multi-group SP₃-ACMFD static equations, an implicit approach is resorted to ensure stability of time steps. An exponential behavior is assumed in transverse integrated equations to establish a relationship between flux moments and currents. Also, analytic integration is benefited for the time-dependent solution of precursor concentration equations. Finally, a multi-channel one-phase thermal hydraulic model is coupled to the proposed methodology. Transient equations are then solved at each step using the GMRES technique. To show the sufficiency of proposed transient SP₃-ACMFD approximation for a full core analysis, a comparison is made using transport peers as the reference. To further demonstrate superiority, results are compared with a 3D multi-group transient diffusion solver developed as a byproduct of this work. Outcomes confirm that the idea can be considered as an economic interim approach which is superior to the diffusion approximation, and comparable with transport in results.

• Nuclear Engineering and Technology
• /
• v.53 no.4
• /
• pp.1079-1087
• /
• 2021

NTP-ERSN is a package developed for solving the multigroup form of the discrete ordinates, characteristics and collision probability of the Boltzmann transport equation in one-dimensional cartesian geometry, by combining pin cells. In this work, C5G7 MOX benchmark is used to verify the accuracy and efficiency of NTP-ERSN package, by treating reactor core problems without spatial homogenization. This benchmark requires solutions in the form of normalized pin powers as well as the vectors and the eigenvalue. All NTP-ERSN simulations are carried out with appropriate spatial and angular approximations. A good agreement between NTP-ERSN results with those obtained with OpenMC calculation code for seven energy groups. In addition, our studies about angular and mesh refinements are carried out to produce better quality solution. Moreover, NTP-ERSN GUI has also been updated and adapted to python 3 programming language.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.469-476
• /
• 2020

Research reactors in-core experimental facilities are designed to provide the highest steady state flux for user's irradiation requirements. However, fuel conversion from highly enriched uranium (HEU) to low enriched uranium (LEU) driven by the ongoing effort to diminish proliferation risk, will impact reactor physics parameters. Preserving the reactor capability to produce the needed flux to perform its intended research functions, determines the conversion feasibility. This study investigates the neutron flux in the central experimental facility of two material test reactors (MTR), the IAEA generic10 MW benchmark reactor and the 22 MW s Egyptian Test and Research Reactor (ETRR-2). A 3D full core model with three uranium enrichment of 93%, 45%, and 20% was constructed utilizing the OpenMC particle transport Monte Carlo code. Neutronics calculations were performed for fresh fuel, the beginning of life cycle (BOL) and end of life cycle (EOL) for each of the three enrichments for both the IAEA 10 MW generic reactor and core 1/98 of the ETRR-2 reactor. Criticality calculations of the effective multiplication factor (K_eff) were executed for each of the twelve cases; results show a reasonable agreement with published benchmark values for both reactors. The thermal, epithermal and fast neutron fluxes were tallied across the core, utilizing the mesh tally capability of the code and are presented here. The axial flux in the central experimental facility was tallied at 1 cm intervals, for each of the cases; results for IAEA 10 MW show a maximum reduction of 14.32% in the thermal flux of LEU to that of the HEU, at EOL. The reduction of the thermal flux for fresh fuel was between 5.81% and 9.62%, with an average drop of 8.1%. At the BOL the thermal flux showed a larger reduction range of 6.92%-13.58% with an average drop of 10.73%. Furthermore, the fission reaction rate was calculated, results showed an increase in the peak fission rate of the LEU case compared to the HEU case. Results for the ETRR-2 reactor show an average increase of 62.31% in the thermal flux of LEU to that of the HEU due to the effect of spectrum hardening. The fission rate density increased with enrichment, resulting in 34% maximum increase in the HEU case compared to the LEU case at the assemblies surrounding the flux trap.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.18 no.2
• /
• pp.43-49
• /
• 2022

The purpose of this study is to evaluate the radiation effect on damage when the external shield of the spent nuclear fuel transport cask is damaged due to impact as the cause of an unexpected accident. The neutron and gamma-ray intensities and spectra are calculated using the ORIGEN-Arp module in the SCALE 6.2.4 code package⁽¹⁾ and then using MCNP6.2⁽²⁾ code calculate the dose rate. In order to evaluate the radiation dose according to the size of damage caused by external impact, various sized holes of 0.3~13.7% are assumed in the outer shield of the cask to evaluate the sensitivity to the dose. In the case of radiation source leakage, damage to the nuclear fuel assembly is assumed to be up to 6% based on overseas test cases. When only the outer shield is damaged, the maximum surface dose is calculated as 3.12E+03 mSv/hr. However, if the radiation source is leaked due to damage to the nuclear fuel assembly, it becomes 7.00E+05 mSv/hr which is about 200 times greater than the former case.