High-resolution neutronics model for ^(238)Pu production in high-flux reactors

We proposed and compared three methods(filter burnup,single energy burnup,and burnup extremum analysis)to build a high-resolution neutronics model for 238Pu production in high-flux reactors.The filter burnup and single energy burnup methods have no theoretical approximation and can achieve a spectrum resolution of up to~1 eV,thereby constructing the importance curve and yield curve of the full energy range.The burnup extreme analysis method combines the importance and yield curves to consider the influence of irradiation time on production efficiency,thereby constructing extreme curves.The three curves,which quantify the transmutation rate of the nuclei in each energy region,are of physical significance because they have similar distributions.A high-resolution neutronics model for ^(238)Pu production was established based on these three curves,and its universality and feasibility were proven.The neutronics model can guide the neutron spectrum optimization and improve the yield of ^(238)Pu by up to 18.81%.The neutronics model revealed the law of nuclei transmutation in all energy regions with high spectrum resolution,thus providing theoretical support for high-flux reactor design and irradiation production of ^(238)Pu.

Neutronics analysis of a subcritical blanket system driven by a gas dynamic trap-based fusion neutron source for ^(99)Mo production

Gamma-emitting radionuclide ^(99m)Tc is globally used for the diagnosis of various pathological conditions owing to its ideal single-photon emission computed tomography (SPECT) characteristics.However,the short half-life of ^(99m)Tc (T_(1/2)=6 h)makes it difficult to store or transport.Thus,the production of ^(99m)Tc is tied to its parent radionuclide ^(99)Mo (T_(1/2)=66 h).The major production paths are based on accelerators and research reactors.The reactor process presents the potential for nuclear proliferation owing to its use of highly enriched uranium (HEU).Accelerator-based methods tend to use deuterium–tritium(D–T) neutron sources but are hindered by the high cost of tritium and its challenging operation.In this study,a new ^(99)Mo production design was developed based on a deuterium–deuterium (D–D) gas dynamic trap fusion neutron source (GDT-FNS) and a subcritical blanket system (SBS) assembly with a low-enriched uranium (LEU) solution.GDT-FNS can provide a relatively high-neutron intensity,which is one of the advantages of ^(99)Mo production.We provide a Monte Carlo-based neutronics analysis covering the calculation of the subcritical multiplication factor (k_(s)) of the SBS,optimization design for the reflector,shielding layer,and ^(99)Mo production capacity.Other calculations,including the neutron flux and nuclear heating distributions,are also provided for an overall evaluation of the production system.The results demonstrated that the SBS meets the nuclear critical safety design requirement (k_(s)<0.97) and maintained a high ^(99)Mo production capacity.The proposed system can generate approximately 157 Ci ^(99)Mo for a stable 24 h operation with a neutron intensity of 1×10^(14) n/s,which can meet 50%of China’s demand in 2025.

Neutronics analysis of a stacked structure for a subcritical system with LEU solution driven by a D-T neutron source for~(99)Mo production

The utilization of neutrons markedly affects the medical isotope yield of a subcritical system driven by an external D-T neutron source.The general methods to improve the utilization of neutrons include moderating multiplying,and reflecting neutrons,which ignores the use of neutrons that backscatter to the source direction.In this study,a stacked structure was formed by assembling the multiplier and the low-enriched uranium solution to enable the full use of neutrons that backscatter to the source direction and further improve the utilization of neutrons.A model based on SuperMC was used to evaluate the neutronics and safety behavior of the subcritical system,such as the neutron effective multiplication factor,neutron energy spectrum,medical isotope yield,and heat deposition.Based on the calculation results,when the intensity of the neutron source was 59×10^(13)n/s,the optimized design with a stacked structure could increase the yield of ^(99)Mo to182 Ci/day,which is approximately 16% higher than that obtained with a single-layer structure.The inlet H_(2)O coolant velocity of 1.0 m/s and initial temperature of 20℃ were also found to be sufficient to prevent boiling of the fuel solution.

Updated neutronics analyses of a water cooled ceramic breeder blanket for the CFETR

The water cooled ceramic breeder（WCCB） blanket employing pressurized water as a coolant is one of the breeding blanket candidates for the China Fusion Engineering Test Reactor（CFETR）.Some updating of neutronics analyses was needed, because there were changes in the neutronics performance of the blanket as several significant modifications and improvements have been adopted for the WCCB blanket, including the optimization of radial build-up and customized structure for each blanket module. A 22.5 degree toroidal symmetrical torus sector 3 D neutronics model containing the updated design of the WCCB blanket modules was developed for the neutronics analyses. The tritium breeding capability, nuclear heating power, radiation damage,and decay heat were calculated by the MCNP and FISPACT code. The results show that the packing factor and^6 Li enrichment of the breeder should both be no less than 0.8 to ensure tritium self-sufficiency. The nuclear heating power of the blanket under 200 MW fusion power reaches201.23 MW. The displacement per atom per full power year（FPY） of the plasma-facing component and first wall reach 0.90 and 2.60, respectively. The peak H production rate reaches150.79 appm/FPY and the peak He production reaches 29.09 appm/FPY in blanket module #3.The total decay heat of the blanket modules is 2.64 MW at 1 s after shutdown and the average decay heat density can reach 11.09 kW m^（-3） at that time. The decay heat density of the blanket modules slowly decreases to lower than 10 W m^（-3） in more than ten years.

Implementation of high-fidelity neutronics and thermal–hydraulic coupling calculations in HNET

To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.

Neutronics analysis for MSR cell with different fuel salt channel geometries

The neutronic properties of molten salt reactors(MSRs)differ from those of traditional solid fuel reactors owing to their nuclear fuel particularity.Based on the Monte-Carlo N particle transport code,the effects of the size and shape of the fuel salt channel on the neutron physics of an MSR cell are investigated systematically in this study.The results show that the infinite multiplication factor(k?)first increases and then decreases with the change in the graphite cell size under certain fuel volume fraction(FVF)conditions.For the same FVF and average chord length,when the average chord length is relatively small,the k?values for different fuel salt channel shapes agree well.When the average chord length is relatively large,the k?values for different fuel salt channel shapes differ significantly.In addition,some examples of practical applications of this study are presented,including cell selection for the core and thermal expansion displacement analysis of the cell.

Neutronics Optimization of LiPb-He Dual-Cooled Fuel Breeding Blanket for the Fusion-Driven sub-critical System

The concept of the liquid Li17Pb83 and Helium gas dual-cooled Fuel Breeding Blanket (FBB) for the Fusion-Driven sub-critical System (FDS) is presented and analyzed. Taking self-sustaining tritium (TBR >1.05) and annual output of 100 kg or more fissile 239Pu (FBR > 0.238) as objective parameters, and based on the three-dimensional Monte Carlo neutron-photon transport code MCNP/4A, a neutronics-optimizated calculation of different cases was carried out and the concept is proved feasible. In addition, the total breeding ratio ( BR = TBR + FBR ) is listed corresponding to different cases.

Neutronics Optimization of Tritium Breeding Blan-ket for the FDS

Neutronics optimization calculations have been performed for the tritium breed-ing blankets with solid ceramic breeder Li2O and 1iquid eutectic breeder Lil7Pb83, respectively,based on a 2-D geometrical configuration using the Monte Carlo neutron-photon transport codeMCNP/4B. The effects of beryllium, 6Li enrichment and various structural materials on TritiumBreeding Ratio have been systematically analyzed.

Analysis of Neutronics and Thermal-Hydraulic Behavior in a Fuel Pin of Pressurized Water Reactor (PWR)

This paper presents a comparative analysis of different parameters such as enthalpy, moderator temperature, moderator density, flow velocity, pressure, and fuel temperature profile at the fuel pin cell level of PWR. Moreover, in this paper pitches to fuel pin radius ratio are varied from 2.3 to 4. The methods and implementation strategy are such that the coupled neutronic and thermal-hydraulic analysis is executed in a fully one dimensional (1D) manner. The thermal hydraulic is based on moderator/coolant mass and enthalpy equation together with one group diffusion equation for fuel pin. Modelling of fuel pin cell and subchannel is executed in two steps. First, the governing equations are derived assuming that all the parameters appearing in the equations are temperature independent. Fuel pin centerline temperature and radially averaged temperature equations are derived from Fourier laws of thermal conductivity. Finally, diffusion coefficient, fission cross-section and absorbing cross-section are evaluated with respect to the fuel pin temperature. The outcome will be helpful for further neutronics and thermal analysis of PWR. Thermal hydraulics parameter varies the maximum 30 percentage from the lowermost value.

Neutronics Comparison Analysis of the Water Cooled Ceramics Breeding Blanket for CFETR

China Fusion Engineering Test Reactor（CFETR） is an ITER-like fusion engineering test reactor that is intended to fill the scientific and technical gaps between ITER and DEMO.One of the main missions of CFETR is to achieve a tritium breeding ratio that is no less than 1.2to ensure tritium self-sufficiency.A concept design for a water cooled ceramics breeding blanket（WCCB） is presented based on a scheme with the breeder and the multiplier located in separate panels for CFETR.Based on this concept,a one-dimensional（1D） radial built breeding blanket was first designed,and then several three-dimensional models were developed with various neutron source definitions and breeding blanket module arrangements based on the 1D radial build.A set of nuclear analyses have been carried out to compare the differences in neutronics characteristics given by different calculation models,addressing neutron wall loading（NWL）,tritium breeding ratio（TBR）,fast neutron flux on inboard side and nuclear heating deposition on main in-vessel components.The impact of differences in modeling on the nuclear performance has been analyzed and summarized regarding the WCCB concept design.

Preliminary performance analysis and optimization based on 1D neutronics model for Indian DEMO HCCB blanket

India,under its breeding blanket R&D program for DEMO,is focusing on the development of two tritium breeding blanket concepts;namely the lead-lithium-cooled ceramic breeder and the helium-cooled ceramic breeder(HCCB).The study presented in this paper focuses on the neutronic design analysis and optimization from the tritium breeding perspective of the HCCB blanket.The Indian concept has an edge-on configuration and is one of the variants of the helium-cooled solid breeder blanket concepts proposed by several partner countries in ITER.The Indian HCCB blanket having lithium titanate(Li2TiO3)as the tritium breeder and beryllium(Be)as the neutron multiplier with reduced-activation ferritic/martensitic steel structure aims at utilizing the low-energy neutrons at the rear part of the blanket.The aim of the optimization study is to minimize the radial blanket thickness while ensuring tritium self-sufficiency and provide data for further neutronic design and thermal-hydraulic layout of the HCCB blanket.It is found that inboard and outboard blanket thicknesses of 40 cm and 60 cm,respectively,can give a tritium breeding ratio(TBR)>1.3,with 60%6Li enrichment,which is assumed to be sufficient to cover potential tritium losses and associated uncertainties.The results also demonstrated that the Be packing fraction(PF)has a more profound impact on the TBR as compared to 6Li enrichment and the PF of Li2TiO3.

Application of local Monte Carlo method in neutronics calculation of EAST radial neutron camera

The Local Monte Carlo(LMC)method is used to solve the problems of deep penetration and long time in the neutronics calculation of the radial neutron camera(RNC)diagnostic system on the experimental advanced superconducting tokamak(EAST),and the radiation distribution of the RNC and the neutron flux at the detector positions of each channel are obtained.Compared with the results calculated by the global variance reduction method,it is shown that the LMC calculation is reliable within the reasonable error range.The calculation process of LMC is analyzed in detail,and the transport process of radiation particles is simulated in two steps.In the first step,an integrated neutronics model considering the complex window environment and a neutron source model based on EAST plasma shape are used to support the calculation.The particle information on the equivalent surface is analyzed to evaluate the rationality of settings of equivalent surface source and boundary.Based on the characteristic that only a local geometric model is needed in the second step,it is shown that the LMC is an advantageous calculation method for the nuclear shielding design of tokamak diagnostic systems.

An improved analysis method for assessing the nuclear-heating impact on the stability of toroidal field magnets in fusion reactors

The superconducting magnet system of a fusion reactor plays a vital role in plasma confinement,a process that can be dis-rupted by various operational factors.A critical parameter for evaluating the temperature margin of superconducting magnets during normal operation is the nuclear heating caused by D-T neutrons.This study investigates the impact of nuclear heat-ing on a superconducting magnet system by employing an improved analysis method that combines neutronics and thermal hydraulics.In the magnet system,toroidal field(TF)magnets are positioned closest to the plasma and bear the highest nuclear-heat load,making them prime candidates for evaluating the influence of nuclear heating on stability.To enhance the modeling accuracy and facilitate design modifications,a parametric TF model that incorporates heterogeneity is established to expedite the optimization design process and enhance the accuracy of the computations.A comparative analysis with a homogeneous TF model reveals that the heterogeneous model improves accuracy by over 12%.Considering factors such as heat load,magnetic-field strength,and cooling conditions,the cooling circuit facing the most severe conditions is selected to calculate the temperature of the superconductor.This selection streamlines the workload associated with thermal-hydraulic analysis.This approach enables a more efficient and precise evaluation of the temperature margin of TF magnets.Moreover,it offers insights that can guide the optimization of both the structure and cooling strategy of superconducting magnet systems.

Preliminary safety analysis for heavy water-moderated molten salt reactor

The heavy water-moderated molten salt reactor(HWMSR)is a newly proposed reactor concept,in which heavy water is adopted as the moderator and molten salt dissolved with fissile and fertile elements is used as the fuel.Issues arising from graphite in traditional molten salt reactors,including the positive temperature coefficient and management of highly radio-active spent graphite waste,can be addressed using the HWMSR.Until now,research on the HWMSR has been centered on the core design and nuclear fuel cycle to explore the viability of the HWMSR and its advantages in fuel utilization.However,the core safety of the HWMSR has not been extensively studied.Therefore,we evaluate typical accidents in a small modular HWMSR,including fuel salt inlet temperature overcooling and overheating accidents,fuel salt inlet flow rate decrease,heavy water inlet temperature overcooling accidents,and heavy water inlet mass flow rate decrease accidents,based on a neutronics and thermal-hydraulics coupled code.The results demonstrated that the core maintained safety during the investigated accidents.

Steady state investigation on neutronics of a molten salt reactor considering the flow effect of fuel salt

The Molten Salt Reactor （MSR）, one of the ‘Generation Ⅳ＇ concepts, is a liquid-fuel reactor, which is different from the conventional reactors using solid fissile materials due to the flow effect of fuel salt. The study on its neutronics considering the fuel salt flow, which is the base of the thermal-hydraulic calculation and safety analysis, must be done. In this paper, the theoretical model on neutronics under steady condition for a single-liquid-fueled MSR is conducted and calculated by numerical method. The neutronics model consists of two group neutron diffusion equations for fast and thermal neutron fluxes, and balance equations for six-group delayed neutron precursors considering the flow effect of fuel salt. The spatial discretization of the above models is based on the finite volume method, and the discretization equations are computed by the source iteration method. The distributions of neutron fluxes and the distributions of the delayed neutron precursors in the core are obtained. The numerical calculated results show that, the fuel salt flow has little effect on the distribution of fast and thermal neutron fluxes and the effective multiplication factor; however, it affects the distribution of the delayed neutron precursors significantly, especially the long-lived one. In addition, it could be found that the delayed neutron precursors influence the neutronics slightly under the steady condition.

Performance optimization of the neutron-sensitive image intensifier used in neutron imaging

As a non-destructive testing technology,neutron imaging plays an important role in various fields,including material science,nuclear engineering,and fundamental science.An imaging detector with a neutron-sensitive image intensifier has been developed and demonstrated to achieve good spatial resolution and timing resolution.However,the influence of the working voltage on the performance of the neutron-sensitive imaging intensifier has not been studied.To optimize the performance of the neutron-sensitive image intensifier at different voltages,experiments have been performed at the China Spallation Neutron Source(CSNS)neutron beamline.The change in the light yield and imaging quality with different voltages has been acquired.It is shown that the image quality benefits from the high gain of the microchannel plate(MCP)and the high accelerating electric field between the MCP and the screen.Increasing the accelerating electric field is more effective than increasing the gain of MCPs for the improvement of the imaging quality.Increasing the total gain of the MCP stack can be realized more effectively by improving the gain of the standard MCP than that of the n MCP.These results offer a development direction for image intensifiers in the future.

Design and implementation of the monochromator shielding for the cold neutron spectrometers XINGZHI and BOYA

An innovative monochromator shielding is designed and implemented for the cold neutron spectrometers XINGZHI and BOYA operated by Renmin University of China at China Advanced Research Reactor.Via Monte Carlo simulations and careful mechanical designs,a shielding configuration has been successfully developed to satisfy safety requirements of below 3μSv/h dose rate at its exterior,meanwhile fulfilling space,floor load and nonmagnetic requirements.Composite materials are utilized to form the sandwich-type shielding walls:the inner layer of boron carbide rubber,the middle layer of steel-encased lead and the outer layer of borated polyethylene.Special-shaped liftable shielding blocks are incorporated to facilitate a continuous adjustment of the neutron energy while preventing radiation leakage.Our work has demonstrated that by utilizing composite shielding materials,along with the sandwich structure and liftable shielding blocks,a compact and lightweight shielding solution can be achieved.This enables the realization of advanced neutron scattering instruments that provide expanded space of measurement,larger energy and momentum coverage,and higher flux on the sample.This shielding represents the first of its kind in neutron scattering instruments in China.Following its successful operation,it has been subsequently employed by other neutron instruments across the country.

Benchmark experiment on slab^(238)U with D-T neutrons for validation of evaluated nuclear data

A benchmark experiment on^(238)U slab samples was conducted using a deuterium-tritium neutron source at the China Institute of Atomic Energy.The leakage neutron spectra within energy levels of 0.8-16 MeV at 60°and 120°were measured using the time-of-flight method.The samples were prepared as rectangular slabs with a 30 cm square base and thicknesses of 3,6,and 9 cm.The leakage neutron spectra were also calculated using the MCNP-4C program based on the latest evaluated files of^(238)U evaluated neutron data from CENDL-3.2,ENDF/B-Ⅷ.0,JENDL-5.0,and JEFF-3.3.Based on the comparison,the deficiencies and improvements in^(238)U evaluated nuclear data were analyzed.The results showed the following.(1)The calculated results for CENDL-3.2 significantly overestimated the measurements in the energy interval of elastic scattering at 60°and 120°.(2)The calculated results of CENDL-3.2 overestimated the measurements in the energy interval of inelastic scattering at 120°.(3)The calculated results for CENDL-3.2 significantly overestimated the measurements in the 3-8.5 MeV energy interval at 60°and 120°.(4)The calculated results with JENDL-5.0 were generally consistent with the measurement results.

Nondestructive technique for identifying nuclides using neutron resonance transmission analysis at CSNS Back-n

Nondestructive and noninvasive neutron assays are essential applications of neutron techniques.Neutron resonance transmission analysis(NRTA)is a powerful nondestructive method for investigating the elemental composition of an object.The back-streaming neutron line(Back-n)is a newly built time-of-flight facility at the China Spallation Neutron Source(CSNS)that provides neutrons in the eV to 300 MeV range.A feasibility study of the NRTA method for nuclide identification was conducted at the CSNS Back-n via two test experiments.The results demonstrate that it is feasible to identify different elements and isotopes in samples using the NRTA method at Back-n.This study reveals its potential future applications.

An efficient parallel algorithm of variational nodal method for heterogeneous neutron transport problems

The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.