• Nuclear Engineering and Technology
• /
• v.50 no.8
• /
• pp.1339-1348
• /
• 2018

The aim of this paper is to determine neutronic performances of the light water reactor (LWR) spent fuel mixed with fertile thorium fuel in a FFHR. Time dependent three dimensional calculations for major technical data, such as blanket energy multiplication, tritium breeding ratio, cumulative fissile fuel enrichment and burnup have been performed by using Monte Carlo Neutron-Particle Transport code MCNP5 1.4, coupled with a novel interface code MCNPAS, which is developed by our research group. A self-sustaining tritium breeding ratio (TBR>1.05) has been kept throughout the calculations. The study has shown that the fissile fuel quality will be improved in the course of the transmutation of the LWR spent in the FFHR. The latter has gained the reusable fuel enrichment level conventional LWRs between one and two years. Furthermore, LWR spent fuel - thorium mixture provides higher burn-up values than in light water reactors.

• 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.55 no.1
• /
• pp.324-329
• /
• 2023

Different types of deterministic solution methods were used to solve neutron transport equations corresponding to half-space and slab albedo problems. In these types of solution methods, in addition to the error of the numerical solutions, the obtained results contain truncation and discretization errors. In the present work, a non-analog Monte Carlo method is provided to simulate the half-space and slab albedo problems with Inönü and Anlı-Güngör strongly anisotropic scattering functions. For each scattering function, the sampling method of the direction of the scattered neutrons is presented. The effects of different beams with different angular dependencies and the effects of different scattering parameters on the reflection probability are investigated using the developed Monte Carlo method. The validity of the Monte Carlo method is also confirmed through the comparison with the published data.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1210-1217
• /
• 2023

The assessment of radiation fields in the lithium loop pipes and dump tank during the operation were performed for International Fusion Materials Irradiation Facility - DEMO-Oriented NEutron Source (IFMIF-DONES) in order to obtain the radiation dose-rate maps in the component surroundings. Variance reduction techniques such as weight window mesh (produced with the ADVANTG code) were applied to bring the statistical uncertainty down to a reasonable level. The biological dose was given in the study, and potential shielding optimization is suggested and more thoroughly evaluated. The MCNP Monte Carlo was used to simulate a gamma particle transport for radiation shielding purposes for the current Li Systems' design. In addition, the shielding efficiency was identified for the Impurity Control System components and the dump tank. The analysis reported in this paper takes into account the radiation decay source from and activated corrosion products (ACPs), which is created by d-Li interaction. As a consequence, the radiation (resulting from ACPs and Be-7) shielding calculations have been carried out for safety considerations.

• Journal of Radiation Protection and Research
• /
• v.18 no.2
• /
• pp.47-59
• /
• 1993

The criticality of the shipping cask(KSC-7) for transportion of 7PWR spent fuel assemblies has been calculated and analysised on the basis of neutron transport theory. For criticality analysis, effects of the rod pitches, the fixed neutron absorbers(borated sus+boral) were considered. The effective multiplication factor has been calculated by KENO-Va, Mote Carlo method computer code, with the HANSEN-ROACH 16 group cross section set, which was made for personal computer system. The criticality for the KSC-7 cask was calculated in terms of the fresh fuel which was conservative for the aspects of nuclear critility. From the results of criticality analysis, the calculated Keff is proved to be lower than subcritical limit during normal transportation and under hypothetical accident condition. The maximum calculated criticalities of the KSC-7 were lower the safety criticality limit 1.0 recommended by US 10CFR71 both under normal and hypothetical accident condition. Also, to verify the KSC-7 criticality calculation results by using KENO-Va, it was carried out benchmark calculation with experimental data of B & W(Bobcock and Wilcox) company. From the 3s series of calculation of the KSC-7 cask and benchmark calculation, the cask was safely designed in nuclear criticality, respectively.

• Nuclear Engineering and Technology
• /
• v.54 no.4
• /
• pp.1253-1260
• /
• 2022

This paper reports the control method for the core power of the China initiative Accelerator Driven System (CiADS) facility. In the CiADS facility, an intense external neutron source provided by a proton accelerator coupled to a spallation target is used to drive a sub-critical reactor. Without any control rod inside the sub-critical reactor, the core power is controlled by adjusting the proton beam intensity. In order to continuously change the beam intensity, an adjustable aperture is considered to be used at the Low Energy Beam Transport (LEBT) line of the accelerator. The aperture size is adjusted based on the Proportional Integral Derivative (PID) controllers, by comparing either the setting beam intensity or the setting core power with the measured value. To evaluate the proposed control method, a CiADS core model is built based on the point reactor kinetics model with six delayed neutron groups. The simulations based on the CiADS core model have indicated that the core power can be controlled stably by adjusting the aperture size. The response time in the adjustment of the core power depends mainly on the adjustment time of the beam intensity.

• Nuclear Engineering and Technology
• /
• v.54 no.2
• /
• pp.523-531
• /
• 2022

In order to study the sensitivity and the uncertainty of the Moroccan research reactor TRIGA Mark II, a model of this reactor has been developed in our ERSN laboratory for use with the N-Particle MCNP Monte Carlo transport codes (version 6). In this article, the sensitivities of the effective multiplication factor of this reactor are evaluated using the ENDF/B-VII.0, ENDF/B-VII.1 and JENDL-4.0 libraries and in 44 energy groups, for the cross sections of the fuel (U-235 and U-238) and the moderator (H-1 and O-16). However, the quantification of the uncertainty of the nuclear data is performed using the nuclear code NJOY99 for the generation and processing of covariance matrices. On the one hand, the highest uncertainty deviations, calculated using the ENDFB-VII.1 and JENDL4.0 evaluations, are 2275, 386 and 330 pcm respectively for the reactions U₂₃₅(n, f), $ U_{235}(n\bar{\nu})$ and H₁(n, γ). On the other hand, these differences are very small for the neutron reactions of O-16 and U-238. Regarding the neutron spectra, in CT-mid plane, they are very close for the three evaluations (ENDF/B-VII.0, ENDF/B-VII.1 and JENDL-4.0). These spectra present two peaks (thermal and fission) around the energies 0.05 eV and 1 MeV.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2586-2591
• /
• 2022

In the present work, we have irradiated the DI-BSCCO superconducting tapes with the 100 keV deuterium ions to investigate the effect of ion irradiation on their critical current (Ic). The damage simulations are carried out using the binary collision approximation method to get the spatial distribution and depth profile of the damage events in the high temperature superconducting (HTS) tape. The point defects are formed near the surface of the HTS tape. These point defects change the vortex profile in the superconducting tape. Due to the long-range interaction of vortices with each other, the Ic of the tape degrades at the 77 K and self magnetic field. The radiation dose of 2.90 MGy degrades the 44% critical current of the tape. The results of the displacement per atom (dpa) and dose deposited by the deuterium ions are used to fit an empirical relation for predicting the degradation of the Ic of the tape. We include the dpa, dose and columnar defect terms produced by the incident particles in the empirical relation. The fitted empirical relation predicts that light ion irradiation degrades the Ic in the DI-BSCCO tape at the self field. This empirical relation can also be used in neutron irradiation to predict the lifetime of the DI-BSCCO tape. The change in the Ic of the DI-BSCCO tape due to deuterium irradiation is compared with the other second-generation HTS tape irradiated with energetic radiation.

• Journal of Radiation Protection and Research
• /
• v.45 no.3
• /
• pp.130-141
• /
• 2020

Background: Concrete activation in cyclotron vaults is a major concern associated with their decommissioning because a considerable amount of activated concrete is generated by secondary neutrons during the operation of cyclotrons. Reducing the amount of activated concrete is important because of the high cost associated with radioactive waste management. This study aims to investigate the capability of the neutron absorbing materials to reduce concrete activation. Materials and Methods: The Particle and Heavy Ion Transport code System (PHITS) code was used to simulate a cyclotron target and room. The dimensions of the room were 457 cm (length), 470 cm (width), and 320 cm (height). Gd₂O₃, B₄C, polyethylene (PE), and borated (5 wt% ^natB) PE with thicknesses of 5, 10, and 15 cm and their different combinations were selected as neutron absorbing materials. They were placed on the concrete walls to determine their effects on thermal neutrons. Thin B₄C and Gd₂O₃ were placed between the concrete wall and additional PE shield separately to decrease the required thickness of the additional shield, and the thermal neutron flux at certain depths inside the concrete was calculated for each condition. Subsequently, the optimum combination was determined with respect to radioactive waste reduction, price, and availability, and the total reduced radioactive concrete waste was estimated. Results and Discussion: In the specific conditions considered in this study, the front wall with respect to the proton beam contained radioactive waste with a depth of up to 64 cm without any additional shield. A single layer of additional shield was inefficient because a thick shield was required. Two-layer combinations comprising 0.1- or 0.4-cm-thick B₄C or Gd₂O₃ behind 10 cm-thick PE were studied to verify whether the appropriate thickness of the additional shield could be maintained. The number of transmitted thermal neutrons reduced to 30% in case of 0.1 cm-thick Gd₂O₃+10 cm-thick PE or 0.1 cm-thick B₄C+10 cm-thick PE. Thus, the thickness of the radioactive waste in the front wall was reduced from 64 to 48 cm. Conclusion: Based on price and availability, the combination of the 10 cm-thick PE+0.1 cmthick B₄C was reasonable and could effectively reduce the number of thermal neutrons. The amount of radioactive concrete waste was reduced by factor of two when considering whole concrete walls of the PET cyclotron vault.

• Nuclear Engineering and Technology
• /
• v.20 no.1
• /
• pp.65-70
• /
• 1988

A Shielding analysis of the shipping cask designed conceptually, of which shielding material are lead and resin, for containing 4 PWR spent fuel assemblies, has been made with the help of a computer code, ANISN. The shielding materials being used in the cask have been selected and arranged to minimize cask weight while maintaining an overall shielding effectiveness. Radiation source terms have been calculated by means of ORIGIN-2 code under the assumptions of 38,000 MWD/MTU burnup and 3-year cooling time. A calculation of gamma-ray and neutron dose rates on the cask surface and 1m from the surface has been done. It is revealed that the total dose rates under the normal transport and hypothetical accident conditions meet the standards specified.