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

• Journal of Radiation Protection and Research
• /
• v.32 no.4
• /
• pp.158-167
• /
• 2007

It is required that the neutron dosimetry be present to monitor the reactor vessel throughout its plant life. The Ex-vessel Neutron Dosimetry Systems which consist of sensor sets, radiometric monitors, gradient chains, and support hardware have been installed for 3-Loop plants after a complete withdrawal of all six in-vessel surveillance capsules. The systems have been installed in the reactor cavity annulus in order to characterize the neutron energy spectrum over the beltline region of the reactor vessel. The installed dosimetry were withdrawn and evaluated after a irradiation during one cycle and then compared to the cycle specific neutron transport calculations. The reaction rates from the measurement and calculation were compared and the results show good agreements each other.

• Journal of Radiation Protection and Research
• /
• v.41 no.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
• /
• v.48 no.6
• /
• pp.1291-1302
• /
• 2016

This paper presents a new and efficient scheme to determine the optimal neutron source position in a model near-equilibrium pressurized water reactor, which is based on the OPR1000 Hanul Unit 3 Cycle 7 configuration. The proposed scheme particularly assigns importance of source positions according to the local adjoint flux distribution. In this research, detailed pin-by-pin reactor adjoint fluxes are determined by using the Monte Carlo KENO-VI code from solutions of the reactor homogeneous critical adjoint transport equations. The adjoint fluxes at each allowable source position are subsequently ranked to yield four candidate positions with the four highest adjoint fluxes. The study next simulates ex-core detector responses using the Monte Carlo MAVRIC code by assuming a neutron source is installed in one of the four candidate positions. The calculation is repeated for all positions. These detector responses are later converted into an inverse count rate ratio curve for each candidate source position. The study confirms that the optimal source position is the one with very high adjoint fluxes and detector responses, which is interestingly the original source position in the OPR1000 core, as it yields an inverse count rate ratio curve closest to the traditional 1/M line. The current work also clearly demonstrates that the proposed adjoint flux-based approach can be used to efficiently determine the optimal geometry for a neutron source and a detector in a modern pressurized water reactor core.

• Corrosion Science and Technology
• /
• v.20 no.5
• /
• pp.308-314
• /
• 2021

Neutron dose level at bottom head of a reactor pressure vessel (RPV) was calculated using reactor vessel neutron transport for a Korean nuclear power plant A. At 34 EFPY with a 40-year (2042) design life after plating repair, irradiation fast neutron effect was 6.6x10¹⁵ n/cm². As helium(He) gas can be generated by Ni only at 1/10⁶ level of 5 × 10²¹ n/cm², He generation possibility in the Ni plating layer is very little during 40 years of operation (2042, 34 EFPY). Thermal neutrons can significantly affect the generation of He from Ni metal. At 10 years after a repair, He can be generated at a level of about 0.06 appm, a level that can add general welding repair without any consideration. After 40 years of repair, 9.8 appm of He may be generated. Although this is a rather high value, it is within the range of 0.1 to 10 appm when welding repair can be applied. Clad repair by Ni electroplating technology is expected to greatly improve the operation efficiency by improving the safety and shortening the maintenance period of the nuclear power plant.

• Journal of the Korean Society of Radiology
• /
• v.16 no.5
• /
• pp.527-536
• /
• 2022

The radiation shielding characteristic of neutron shielding material has been studied as the preliminary study in order to design cosmic-ray shielding material. Specially, Soft Magnetic Material, known to be effective in EMP and radiation shielding, has been investigated to check if the material would be applicable to cosmic-ray shielding. In this work, thermal neutron shielding experiment was conducted and the Monte Carlo N-Particle(MCNP) was applied to employ skymap.dat, which is cosmic-ray data embedded in MCNP. As a result, polyethylene, borated polyethylene, and carbon nano tube, containing carbon or hydrogen, have been found to be effective in reduction of neutron flux below 20 MeV (including thermal, epithermal, evaporation). In contrast, the materials composed of iron such as SS316 and Soft Magnetic Material show a good shielding performance in the cascade energy range (above 20 MeV). Since Soft Magnetic Material is consisting of 13% of boron, it can also decrease thermal neutron flux, so it is expected that it would show a significant reduction on the entire range of neutron energy if the Soft Magnetic Material is used with hydrogen and carbon, so called low Z material.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.485-498
• /
• 2020

A diffusion synthetic acceleration (DSA) technique for the S_N transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.223-229
• /
• 2020

The present work proposes a solution to the static Boltzmann transport equation approximated by the simplified P₃ (SP₃) on angular, and the analytic coarse mesh finite difference (ACMFD) for spatial variables. Multi-group SP₃-ACMFD equations in 3D rectangular geometry are solved using the GMRES solution technique. As the core time dependent analysis necessitates the solution of an eigenvalue problem for an initial condition, this work is hence devoted to development and verification of the proposed static SP₃-ACMFD solver. A 3D multi-group static diffusion solver is also developed as a byproduct of this work to assess the improvement achieved using the SP₃ technique. Static results are then compared against transport benchmarks to assess the proximity of SP₃-ACMFD solutions to their full transport peers. Results prove that the approach can be considered as an acceptable interim approximation with outputs superior to the diffusion method, close to the transport results, and with the computational costs less than the full transport approach. The work would be further generalized to time dependent solutions in Part II.

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

• Journal of Radiation Protection and Research
• /
• v.10 no.1
• /
• pp.41-49
• /
• 1985

Neutron Streaming analysis in 1300 MWe pressurized water reactor cavity was performed. In this calculation, the discrete ordinates transport codes, ANISN and DOT 3.5, and the Monte Carlo code, TRIPOLI-02 were used with the coupling code, DOTTRI. In this study IBM 3033 type computer was used. The calculated neutron fluxes and dose rates were compared with the measured data in a 900MWe pressurized water reactor cavity to show a good agreement, although some deviations in the results for each energy group were noticed. These results will be applied in the radiation shielding design of high capacity nuclear power reactors and, to the means of radiation protection in case of the reactor maintenance and the access of the reactor cavity.