Neural Network Based Deep Learning Method for Multi-Dimensional Neutron Diffusion Problems with Novel Treatment to Boundary

In this paper, the artificial neural networks (ANN) based deep learning (DL) techniques were developed to solve the neutron diffusion problems for the continuous neutron flux distribution without domain discretization in advance. Due to its mesh-free property, the DL solution can easily be extended to complicated geometries. Two specific realizations of DL methods with different boundary treatments are developed and compared for accuracy and efficiency, including the boundary independent method (BIM) and boundary dependent method (BDM). The performance comparison on analytic benchmark indicates BDM being the preferred DL method. Novel constructions of trial function are proposed to generalize the application of BDM. For a more in-depth understanding of the BDM on diffusion problems, the influence of important hyper-parameters is further investigated. Numerical results indicate that the accuracy of BDM can reach hundreds of times higher than that of BIM on diffusion problems. This work can provide a new perspective for applying the DL method to nuclear reactor calculations.

Simulation of VVER-1000 Guillotine Large Break Loss of Coolant Accident Using RELAP5/SCDAPSIM/MOD3.5

The safety performance of nuclear power plants (NPPs) is a very important factor in evaluating nuclear energy sustainability. Safety analysis of passive and active safety systems have a positive influence on reactor transient mitigation. One of the common transients is primary coolant leg rupture. This study focused on guillotine large break loss of coolant (LB-LOCA) in one of the reactor vessels, in which cold leg rupture occurred, after establishment of a steady-state condition for the VVER-1000. The reactor responses and performance of emergence core cooling systems (ECCSs) were investigated. The main safety margin considered during this simulation was to check the maximum value of the clad surface temperature, and it was then compared with the design licensing limit of 1474 K. The calculations of event progression used the engineering-level RELAP5/SCDAPSIM/MOD3.5 thermal-hydraulic program, which also provide a more detailed treatment of coolant system thermal hydraulics and core behavior. The obtained results show that actuation of ECCSs at their actuation set points provided core cooling by injecting water into the reactor pressure vessel, as expected. The peak cladding temperature did not overpass the licensing limit during this LB-LOCA transient. The primary pressure above the core decreased rapidly from 15.7 MPa to 1 MPa in less than 10 s, then stabilizes up to the end of transient. The fuel temperature decreased from 847 K to 378 K during the first 30 s of the transient time. The coolant leakage reduced from 9945 kg/s to approximately 461 kg/s during the first 190 s in the transient. Overall, the study shows that, within the design of the VVER-1000, safety systems of the have inherent robustness of containing guillotine LB-LOCA.

A Rate Theory Model of Radiation-Induced Swelling in an Austenitic Stainless Steel

Many rate theory models of cavity (void) swelling have been published over the past 50 years, all having the same, or similar, structures. A rigorous validation of the models has not been possible because of the dearth of information concerning the microstructures that correspond with the swelling data. Whereas the lack of microstructure information is still an issue for historical swelling data, in the past 10–20 years data have been published on the evolution of the microstructure (point defect yields from collision cascades, cavity number densities, and dislocation densities/yield strengths) allowing certain gaps in information to be filled when considering historic swelling data. With reasonable estimates of key microstructure parameters, a standard rate theory model can be applied, and the model parameter space explored, in connection with historical swelling data. By using published data on: (i) yield strength as a function of dose and temperature (to establish an empirical expression for dislocation density evolution); (ii) cavity number densities as a function of temperature; and (iii) freely migrating defect (FMD) production as a function of primary knock-on atom (PKA) spectrum, the necessary parameter and microstructure inputs that were previously unknown can be used in model development. This paper describes a rate-theory model for void swelling of 316 stainless steel irradiated in the EBR-2 reactor as a function of irradiation temperature and neutron dose.

Perspectives on a Severe Accident Consequences—10 Years after the Fukushima Accident

Scientific issues that draw international attention from the public and experts during the last 10 years after the Fukushima accident are discussed. An assessment of current severe accident analysis methodology, impact on the views of nuclear reactor safety, dispute on the safety of fishery products, discharge of radioactive water to the ocean, status of decommissioning, and needs for long-term monitoring of the environment are discussed.

Nuclear Data Sensitivity Study for the EBR-II Fast Reactor Benchmark Using SCALE with ENDF/B-VII.1 and ENDF/B-VIII.0

The EBR-II benchmark, which was recently included in the International Handbook of Evaluated Reactor Physics Benchmark Experiments, served as a basis for assessing the performance of the SCALE code system for fast reactor analyses. A reference SCALE model was developed based on the benchmark specifications. Great agreement was observed between the eigenvalue calculated with this SCALE model and the benchmark eigenvalue. To identify potential gaps and uncertainties of nuclear data for the simulation of various quantities of interest in fast spectrum systems, sensitivity and uncertainty analyses were performed for the eigenvalue, reactivity effects, and the radial power profile of EBR-II using the two most recent ENDF/B nuclear data library releases. While the nominal results are consistent between the calculations with the different libraries, the uncertainties due to nuclear data vary significantly. The major driver of observed uncertainties is the uncertainty of the 235U (n,γ) reaction. Since the uncertainty of this reaction is significantly reduced in the ENDF/B-VIII.0 library compared to ENDF/B-VII.1, the obtained output uncertainties tend to be smaller in ENDF/B-VIII.0 calculations, although the decrease is partially compensated by increased uncertainties in 235U fission and ν¯.

Proximal Policy Optimization for Radiation Source Search

Rapid search and localization for nuclear sources can be an important aspect in preventing human harm from illicit material in dirty bombs or from contamination. In the case of a single mobile radiation detector, there are numerous challenges to overcome such as weak source intensity, multiple sources, background radiation, and the presence of obstructions, i.e., a non-convex environment. In this work, we investigate the sequential decision making capability of deep reinforcement learning in the nuclear source search context. A novel neural network architecture (RAD-A2C) based on the advantage actor critic (A2C) framework and a particle filter gated recurrent unit for localization is proposed. Performance is studied in a randomized 20×20 m convex and non-convex simulation environment across a range of signal-to-noise ratio (SNR)s for a single detector and single source. RAD-A2C performance is compared to both an information-driven controller that uses a bootstrap particle filter and to a gradient search (GS) algorithm. We find that the RAD-A2C has comparable performance to the information-driven controller across SNR in a convex environment. The RAD-A2C far outperforms the GS algorithm in the non-convex environment with greater than 95% median completion rate for up to seven obstructions.

The Stability of Linear Diffusion Acceleration Relative to CMFD

Coarse Mesh Finite Difference (CMFD) is a widely-used iterative acceleration method for neutron transport problems in which nonlinear terms are introduced in the derivation of the low-order CMFD diffusion equation. These terms, including the homogenized diffusion coefficient, the current coupling coefficients, and the multiplicative prolongation constant, are subject to numerical instability when a scalar flux estimate becomes sufficiently small or negative. In this paper, we use a suite of contrived problems to demonstrate the susceptibility of CMFD to failure for each of the vulnerable quantities of interest. Our results show that if a scalar flux estimate becomes negative in any portion of phase space, for any iterate, numerical instability can occur. Specifically, the number of outer iterations required for convergence of the CMFD-accelerated transport problem can increase dramatically, or worse, the iteration scheme can diverge. An alternative Linear Diffusion Acceleration (LDA) scheme addresses these issues by explicitly avoiding local nonlinearities. Our numerical results show that the rapid convergence of LDA is unaffected by the very small or negative scalar flux estimates that can adversely affect the performance of CMFD. Therefore, our results demonstrate that LDA is a robust alternative to CMFD for certain sensitive problems in which CMFD can exhibit reduced effectiveness or failure.

Neutronic Characteristics of ENDF/B-VIII.0 Compared to ENDF/B-VII.1 for Light-Water Reactor Analysis

The Evaluated Nuclear Data File (ENDF)/B-VIII.0 data library was released in 2018. To assess the new data during development and shortly after release, many validation calculations were performed with static, room-temperature benchmarks. Recently, when performing validation of ENDF/B-VIII.0 for pressurized water reactor depletion calculations, a regression in performance compared to ENDF/B-VII.1 was observed. This paper documents extensive benchmark calculations for light-water reactors performed using continuous energy Monte Carlo code with ENDF/B-VII.1 and -VIII.0 and neutronic characteristics of ENDF/B-VIII.0 are discussed and compared to those of ENDF/B-VII.1. It is our recommendation that ENDF/B data library assessment should include reactor-specific benchmark assessments, including depletion cases, such that these types of regressions may be caught earlier in the data development cycle.

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) Using CTF

With the resurgence of interest in molten salt reactors, there is a need for new experiments and modeling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory to be placed in the Advanced Test Reactor (ATR). One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this has on the fuel-salt properties, source term generation, and corrosion control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system, and noble metals (Mo, Tc, Te, etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in CTF (Coolant-Boiling in Rod Arrays—Two Fluids) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design. A noble metal–helium bubble mass transfer model is coupled with the thermal-hydraulic results from CTF to determine the effectiveness of this insoluble fission product (IFP) extraction method for VESIL. Amounts of IFP species extracted to the off-gas system and species distributions in VESIL after a 60-day ATR cycle are reported.

Fourth-Order Adjoint Sensitivity and Uncertainty Analysis of an OECD/NEA Reactor Physics Benchmark: I. Computed Sensitivities

This work extends the investigation of higher-order sensitivity and uncertainty analysis from 3rd-order to 4th-order for a polyethylene-reflected plutonium (PERP) OECD/NEA reactor physics benchmark. Specifically, by applying the 4th-order comprehensive adjoint sensitivity analysis methodology (4th-CASAM) to the PERP benchmark, this work presents the numerical results of the most important 4th-order sensitivities of the benchmark’s total leakage response with respect to the benchmark’s 180 microscopic total cross sections, which includes 180 4th-order unmixed sensitivities and 360 4th-order mixed sensitivities corresponding to the largest 3rd-order ones. The numerical results obtained in this work reveal that the number of 4th-order relative sensitivities that have large values (e.g., greater than 1.0) is far greater than the number of important 1st-, 2nd- and 3rd-order sensitivities. The majority of those large sensitivities involve isotopes 1H and 239Pu contained in the PERP benchmark. Furthermore, it is found that for most groups of isotopes 1H and 239Pu of the PERP benchmark, the values of the 4th-order relative sensitivities are significantly larger than the corresponding 1st-, 2nd- and 3rd-order sensitivities. The overall largest 4th-order relative sensitivity S(4)σt,6g=30,σt,6g=30,σt,6g=30,σt,6g=30=2.720×106 is around 291,000 times, 6350 times and 90 times larger than the corresponding largest 1st-order, 2nd-order and 3rd-order sensitivities, respectively, and the overall largest mixed 4th-order relative sensitivity S(4)σt,630,σt,630,σt,630,σt,530=2.279×105 is also much larger than the largest 2nd-order and 3rd-order mixed sensitivities. The results of the 4th-order sensitivities presented in this work have been independently verified with the results obtained using the well-known finite difference method, as well as with the values of the corresponding symmetric 4th-order sensitivities. The 4th-order sensitivity results obtained in this work will be subsequently used on the 4th-order uncertainty analysis to evaluate their impact on the uncertainties they induce in the PERP leakage response.