• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.106-108
• /
• 2002

Recent advances in radiation transport algorithms, computer hardware performance, and parallel computing make the clinical use of Monte Carlo based dose calculations possible. Monte Carlo treatment planning requires accurate beam information as input to generate accurate dose distributions. The procedures to obtain this accurate beam information are called "commissioning", which includes accelerator head modeling. In this study, we would like to investigate how much accurately Monte Carlo based dose calculations can predict the measured beam data in various conditions. The Siemens 6MV photon beam and the BEAM Monte Carlo code were used. The comparisons including the percentage depth doses and off-axis profiles of open fields and wedges, output factors will be presented.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.635-641
• /
• 2016

Because of prohibitive data storage requirements in large-scale simulations, the memory problem is an obstacle for Monte Carlo (MC) codes in accomplishing pin-wise three-dimensional (3D) full-core calculations, particularly for whole-core depletion analyses. Various kinds of data are evaluated and quantificational total memory requirements are analyzed based on the Reactor Monte Carlo (RMC) code, showing that tally data, material data, and isotope densities in depletion are three major parts of memory storage. The domain decomposition method is investigated as a means of saving memory, by dividing spatial geometry into domains that are simulated separately by parallel processors. For the validity of particle tracking during transport simulations, particles need to be communicated between domains. In consideration of efficiency, an asynchronous particle communication algorithm is designed and implemented. Furthermore, we couple the domain decomposition method with MC burnup process, under a strategy of utilizing consistent domain partition in both transport and depletion modules. A numerical test of 3D full-core burnup calculations is carried out, indicating that the RMC code, with the domain decomposition method, is capable of pin-wise full-core burnup calculations with millions of depletion regions.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.790-802
• /
• 2022

Large-scale reactor simulation often requires the use of Monte Carlo calculation techniques to estimate important reactor parameters. One drawback of these Monte Carlo calculation techniques is they inevitably result in some uncertainty in calculated quantities. The present study includes parametric uncertainty quantification (UQ) and sensitivity analysis (SA) on the Advanced Test Reactor Critical (ATRC) facility housed at Idaho National Laboratory (INL) and addresses some complications due to Monte Carlo uncertainty when performing these analyses. This approach for UQ/SA includes consideration of Monte Carlo code uncertainty in computed sensitivities, consideration of uncertainty from directly measured parameters and a comparison of results obtained from brute-force Monte Carlo UQ versus UQ obtained from a surrogate model. These methodologies are applied to the uncertainty and sensitivity of k_eff for two sets of uncertain parameters involving fuel plate geometry and fuel plate composition. Results indicate that the less computationally-expensive method for uncertainty quantification involving a linear surrogate model provides accurate estimations for k_eff uncertainty and the Monte Carlo uncertainty in calculated k_eff values can have a large effect on computed linear model parameters for parameters with low influence on k_eff.

• 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.53 no.9
• /
• pp.2830-2838
• /
• 2021

This work deals with the assessment of the burnup capabilities of the Serpent Monte Carlo code to predict spent nuclear fuel (SNF) isotopic concentrations for low-enriched uranium (LEU) fuel at different burnup levels up to 47 MWd/kgU. The irradiation of six UO₂ experimental samples in three different VVER-1000 reactor units has been simulated and the predicted concentrations of actinides up to ²⁴⁴Cm have been compared with the corresponding measured values. The results show a global good agreement between calculated and experimental concentrations, in several cases within the margins of the nuclear data uncertainties and in a few cases even within the reported experimental uncertainties. The differences in the performances of the JEFF3.1.1, ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries (NDLs) have also been assessed and the use of the newly released ENDF/B-VIII.0 library has shown an increased accuracy in the prediction of the C/E's for some of the actinides considered, particularly for the plutonium isotopes. This work represents a step forward towards the validation of advanced simulation tools against post irradiation experimental data and the obtained results provide an evidence of the capabilities of the Serpent Monte-Carlo code with the associated modern NDLs to accurately compute SNF nuclide inventory concentrations for VVER-1000 type reactors.

• Nuclear Engineering and Technology
• /
• v.46 no.3
• /
• pp.407-412
• /
• 2014

In this research, we investigated the burnup characteristics and the conversion of fertile $^{232}Th$ into fissile $^{233}U$ in the core of a Sodium-Cooled Fast Reactor (SFR). The SFR fuel assemblies were designed for burning $^{232}Th$ fuel (fuel pin 1) and $^{233}U$ fuel (fuel pin 2) and include mixed minor actinide compositions. Monte Carlo simulations were performed using Serpent Code1.1.19 to compare with CRAM (Chebyshev Rational Approximation Method) and TTA (Transmutation Trajectory Analysis) method in the burnup calculation mode. The total heating power generated in the system was assumed to be 2000 MWth. During the reactor operation period of 600 days, the effective multiplication factor (keff) was between 0.964 and 0.954 and peaking factor is 1.88867.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.477-484
• /
• 2020

The control rod drive mechanism needs to be debugged after reactor fresh fuel loading. It is of great importance to monitor the subcriticality of this process accurately. A modified method was applied to the subcriticality monitoring process, in which only a single control rod cluster was fully withdrawn from the core. In order to correct the error in the results obtained by Neutron Source Multiplication Method, which is based on one point reactor model, Monte Carlo neutron transport code was employed to calculate the fission neutron distribution, the iterated fission probability and the neutron flux in the neutron detector. This article analyzed the effect of a coarse mesh and a fine mesh to tally fission neutron distributions, the iterated fission probability distributions and to calculate correction factors. The subcriticality before and after modification is compared with the subcriticality calculated by MCNP code. The modified results turn out to be closer to calculation. It's feasible to implement the modified NSM method in large local reactivity addition process using Monte Carlo code based on 3D model.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.811-816
• /
• 2022

This work presents the methodology used to expand the capabilities of the Monte Carlo code OpenMC for the calculation of reactor kinetic parameters: effective delayed neutron fraction β_eff and neutron generation time L. The modified code, OpenMC(Time-Dependent) or OpenMC(TD), was then used to calculate the effective delayed neutron fraction by using the prompt method, while the neutron generation time was estimated using the pulsed method, fitting Λ to the decay of the neutron population. OpenMC(TD) is intended to serve as an alternative for the estimation of kinetic parameters when licensed codes are not available. The results obtained are compared to experimental data and MCNP calculated values for 18 benchmark configurations.

• Nuclear Engineering and Technology
• /
• v.53 no.9
• /
• pp.2788-2802
• /
• 2021

Multigroup cross section (MG XS) generation by the UNIST in-house Monte Carlo (MC) code MCS for fast reactor analysis using nodal diffusion codes is reported. The feasibility of the approach is quantified for two sodium fast reactors (SFRs) specified in the OECD/NEA SFR benchmark: a 1000 MW_th metal-fueled SFR (MET-1000) and a 3600 MW_th oxide-fueled SFR (MOX-3600). The accuracy of a few-group XSs generated by MCS is verified using another MC code, Serpent 2. The neutronic steady-state whole-core problem is analyzed using MCS/RAST-K with a 24-group XS set. Various core parameters of interest (core k_eff, power profiles, and reactivity feedback coefficients) are obtained using both MCS/RAST-K and MCS. A code-to-code comparison indicates excellent agreement between the nodal diffusion solution and stochastic solution; the error in the core k_eff is less than 110 pcm, the root-mean-square error of the power profiles is within 1.0%, and the error of the reactivity feedback coefficients is within three standard deviations. Furthermore, using the super-homogenization-corrected XSs improves the prediction accuracy of the control rod worth and power profiles with all rods in. Therefore, the results demonstrate that employing the MCS MG XSs for the nodal diffusion code is feasible for high-fidelity analyses of fast reactors.

• Structural Engineering and Mechanics
• /
• v.17 no.2
• /
• pp.267-279
• /
• 2004

This paper presents an improved Monte Carlo simulation for the probabilistic determination of initial cable forces of cable-stayed bridges under dead loads using the response surfaces method. A response surface (i.e. a quadratic response surface without cross-terms) is used to approximate structural response. The use of the response surface eliminates the need to perform a deterministic analysis in each simulation loop. In addition, use of the response surface requires fewer simulation loops than conventional Monte Carlo simulation. Thereby, the computation time is saved significantly. The statistics (e.g. mean value, standard deviation) of the structural response are calculated through conventional Monte Carlo simulation method. By using Monte Carlo simulation, it is possible to use the existing deterministic finite element code without modifying it. Probabilistic analysis of a truss demonstrates the proposed method' efficiency and accuracy; probabilistic determination of initial cable forces of a cable-stayed bridge under dead loads verifies the method's applicability.