• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1860-1865
• /
• 2019

The confinement degradation of the energetic particles during RMP would be a key issue in success of realizing the successful energy production using fusion plasma, because a 3.5 MeV energetic alpha particle should be able to sustain the burning plasma after the ignition. As KSTAR recent results indicate the generation of high-performance plasma(${\beta}_p{\sim}3$), the confinement of the energetic particles is also an important key aspect in neutral beam driven plasma. In general, the measured absolute value of the neutron intensity is generally used for to estimating the confinement time of energetic particles by comparing it with the theoretical value based on transport calculations. However, the availability of, but for its calculation process, many accurate diagnostic data of plasma parameters such as thermal and incident fast ion density, are essential to the calculation process. In this paper, the time evolution of the neutron signal from an He3 counter during the beam blank has permitted to facilitate the estimation of the slowing down time of energetic particles and the method is applied to investigate the fast ion effect on ELM suppressed KSTAR plasma which is heated by high energy deuterium neutral beams.

• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1769-1780
• /
• 2022

A new method of dose calculation algorithm, called GPU-accelerated Monte Carlo and collapsed cone Convolution (GMCC) was developed to improve the calculation speed of BNCT treatment planning system. The GPU-accelerated Monte Carlo routine in GMCC is used to simulate the neutron transport over whole energy range and the Collapsed Cone Convolution method is to calculate the gamma dose. Other dose components due to alpha particles and protons, are calculated using the calculated neutron flux and reaction data. The mathematical principle and the algorithm architecture are introduced. The accuracy and performance of the GMCC were verified by comparing with the FLUKA results. A water phantom and a head CT voxel model were simulated. The neutron flux and the absorbed dose obtained by the GMCC were consistent well with the FLUKA results. In the case of head CT voxel model, the mean absolute percentage error for the neutron flux and the absorbed dose were 3.98% and 3.91%, respectively. The calculation speed of the absorbed dose by the GMCC was 56 times faster than the FLUKA code. It was verified that the GMCC could be a good candidate tool instead of the Monte Carlo method in the BNCT dose calculations.

• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1071-1083
• /
• 2023

The on-the-fly energy release per fission (OTFK) model is implemented in STREAM to continuously update the Kappa values during the depletion calculation. The explicit neutron and photon energy distribution, which has not been considered in previous STREAM versions, is incorporated into the existing on-the-fly model. The impacts of the modified OTFK model with explicit neutron and photon heating in STREAM on the power distribution, fuel temperature, and other core parameters during depletion with feedback calculations are studied using several problems from the VERA benchmark suit. Overall, the explicit heating calculation provides a better power map for the feedback calculations particularly when strong gamma emitters are present. Generally, the fuel temperature decreases when neutron and photon heating is employed because fission neutrons and gamma rays are transported away from their points of generation. This energy release model in STREAM indicates that gamma energy accounts for approximately 9.5%-10% of the total energy released, and approximately 2.4%-2.6% of the total energy released will be deposited in the coolant for the VERA 5, NuScale, and Yonggwang Unit 3 2D cores.

• Journal of Radiation Protection and Research
• /
• v.38 no.1
• /
• pp.37-43
• /
• 2013

In case of neutron dose measurement using TLDs (thermo-luminescence dosimeters), because the neutron energy dependence of the TLD is very high, the calibration of the energy response according to the characteristics of the neutron spectrum of workplace is required. In the present study, the ambient dose equivalent rates inside and around the Long-Counter (neutron detector) with narrow and complex inside in the neutron field of $^{252}Cf$ were evaluated. The calibration factors to account for the neutron energy dependence of TLDs were established for both the bare and $D_2O$ modulated $^{252}Cf$ neutron beams, respectively. The values of the TLD's measurement were compared with the computational results of the MCNPX (Monte Carlo N-Particles transport code). When using the two calibration factors of the TLD than a single calibration factor, the measured and the calculated values at the point of verification outside and inside the Long-Counter were in more good agreement. This results show that TLD should be calibrated in the reference neutron field similar to workplace situation.

• Nuclear Engineering and Technology
• /
• v.55 no.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.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.182-183
• /
• 2004

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.2018-2025
• /
• 2019

This study systematically compared two hybrid deterministic/Monte Carlo transport codes, ADVANTG/MCNP and MAVRIC, in solving a difficult shielding problem for a real-world spent fuel storage cask. Both hybrid codes were developed based on the consistent adjoint driven importance sampling (CADIS) methodology but with different implementations. The dose rate distributions on the cask surface were of primary interest and their predicted results were compared with each other and with a straightforward MCNP calculation as a baseline case. Forward-Weighted CADIS was applied for optimization toward uniform statistical uncertainties for all tallies on the cask surface. Both ADVANTG/MCNP and MAVRIC achieved substantial improvements in overall computational efficiencies, especially for gamma-ray transport. Compared with the continuous-energy ADVANTG/MCNP calculations, the coarse-group MAVRIC calculations underestimated the neutron dose rates on the cask's side surface by an approximate factor of two and slightly overestimated the dose rates on the cask's top and side surfaces for fuel gamma and hardware gamma sources because of the impact of multigroup approximation. The fine-group MAVRIC calculations improved to a certain extent and the addition of continuous-energy treatment to the Monte Carlo code in the latest MAVRIC sequence greatly reduced these discrepancies. For the two continuous-energy calculations of ADVANTG/MCNP and MAVRIC, a remaining difference of approximately 30% between the neutron dose rates on the cask's side surface resulted from inconsistent use of thermal scattering treatment of hydrogen in concrete.

• Journal of Radiation Protection and Research
• /
• v.26 no.3
• /
• pp.225-229
• /
• 2001

In Korea, an advanced reactor system of 330MWt power called SMART (System integrated Modular Advanced ReacTor) is being developed by KAERI to supply energy for seawater desalination as well as electricity generation. A shielding design of the SMART core with 49 CEDM is established by a two-dimensional discrete ordinates radiation transport analyses. The DORT two-dimensional discrete ordinates transport code is used to evaluate the SMART shielding designs. Three axial regions represent the SMART reactor assembly, each of which is modeled in the R-Z geometry. The BUGLE-96 library is used in the analyses, which consists of 47 neutron and 20 gamma energy groups. The results indicate that the maximum neutron fluence at the bottom of reactor vessel is $5.89 {\times} 10^{17}\;n/cm^2$ and that on the radial surface of reactor vessel is $4.49 {\times} 10^[16}\;n/cm^2$. These results meet the requirement, $1.0 {\times} 10^{20}\;n/cm^2$, in 10 CFR 50.61 and the integrity of SMART reactor vessel during the lifetime of the reactor is confirmed.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.319-328
• /
• 1992

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1240-1249
• /
• 2017

An essential component of the neutron transport solver is the resonance self-shielding calculation used to determine equivalence cross sections. The neutron transport code, MPACT, is currently using the subgroup self-shielding method, in which the method of characteristics (MOC) is used to solve purely absorbing fixed-source problems. Recent efforts incorporating multigroup kernels to the MOC solvers in MPACT have reduced runtime by roughly $2{\times}$. Applying the same concepts for self-shielding and developing a novel lumped parameter approach to MOC, substantial improvements have also been made to the self-shielding computational efficiency without sacrificing any accuracy. These new multigroup and lumped parameter capabilities have been demonstrated on two test cases: (1) a single lattice with quarter symmetry known as VERA (Virtual Environment for Reactor Applications) Progression Problem 2a and (2) a two-dimensional quarter-core slice known as Problem 5a-2D. From these cases, self-shielding computational time was reduced by roughly $3-4{\times}$, with a corresponding 15-20% increase in overall memory burden. An azimuthal angle sensitivity study also shows that only half as many angles are needed, yielding an additional speedup of $2{\times}$. In total, the improvements yield roughly a $7-8{\times}$ speedup. Given these performance benefits, these approaches have been adopted as the default in MPACT.