Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

Selective electrified polyethylene upcycling by pore-modulated pyrolysis

Plastic waste is a increasing problem, accumulating in landfills and the environment. Pyrolysis is a promising and industrially relevant approach for transforming plastic waste into value-added chemicals. However, the selectivity and yield of traditional plastic pyrolysis are poor, with products featuring broad molar mass distributions. Here we report a highly selective, energy-efficient and catalyst-free pyrolysis method that can upcycle plastic into value-added chemicals via pore-modulated pyrolysis. Using a Joule-heated carbon column, we demonstrate the pivotal role of the reactor’s graded porous structure in decreasing the polydispersity of the reaction intermediates, enabling high product selectivity and yield. The decreasing pore size of the reactor modulates the mass transport in an apparent gating effect—preventing high-molar-mass species from exiting the reactor before sufficient pyrolysis has occurred. Using polyethylene as a model reactant, we demonstrate a high yield of 65.9 ± 5.2% and up to 80.8% selectivity toward value-added aviation fuel precursor (C8–C18 hydrocarbons) without the use of any catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated Process Testing of MSR Salt Spill Accidents

Part of the licensing process for new nuclear reactors requires vendors to assess the potential consequences of identified accident scenarios using accident progression modeling. The accident scenario that will likely be evaluated by all molten salt reactor (MSR) developers is a spill of radionuclide-bearing fuel salt onto the reactor containment floor (i.e., a salt spill accident). The development of accident progression models requires experimental data to inform which processes to incorporate, to enable the calculation of parameters to model these processes, and to validate the model predictions. The data should quantify the sensitivities of key processes (e.g., molten salt spreading, heat transfer, containment structure corrosion, radionuclide vaporization, and aerosol generation) towards the initial conditions of the spill, the ambient environment, and the features of the containment. In addition, results from integrated process tests that quantify coupled processes are required to validate systems-level models. This report documents results from integrated process tests conducted on simulated molten salt spill accidents. The generated experimental data simultaneously quantify the heat transfer behavior of the spilled salt, compositional changes to the bulk salt, and the release of surrogate fission products from the spilled salt as aerosol particles. All tests that were conducted used FLiNaK doped with surrogate fission products, and the variables that were evaluated included the initial salt temperature and the concentration of surrogate fission products present in the salt. The major accomplishments of this work include identifying surrogate fuel salt compositions that provide insight into the dispersal behavior of radionuclides of potential significance to the source term, employing previously developed methods and measurement techniques to simultaneously measure key processes, generating data on the coupled processes of molten salt heat transfer and surrogate fission product release as aerosol particles, demonstrating new test methods for real-time monitoring of the flow rate of the spill and aerosol size quantification in an argon atmosphere, and developing a mass transfer model for cesium and iodine release from molten FLiNaK to provide insight into aerosol formation by vapor condensation. The same methodology applied herein can be employed to study different salt compositions of interest to MSR developers, different environmental conditions, and other variables that are relevant to postulated accident scenarios. The insights gained from these integrated process tests conducted at a laboratory scale will be incorporated into future integral effects tests conducted at an engineering scale.

20 FOSSIL-FUELED POWER PLANTS↗

Response Time Analysis and Test of Protection System Instrument Channels for APR1400 and OPR1000

Safety limits are required to maintain the integrity of physical barriers designed to prevent the uncontrolled release of radioactive materials in nuclear power plants. The safety analysis establishes two critical constraints that include an analytical limit in terms of a measured or calculated variable, and a specific time after the analytical limit is reached to begin protective action. Keeping with the nuclear regulations and industry standards, satisfying these two requirements will ensure that the safety limit will not be exceeded during the design basis event, either an anticipated operational occurrence or a postulated accident. Various studies on the setpoint determination methodology for the safety-related instrumentation have been actively performed to ensure that the requirement of the analytical limit is satisfied. In particular, the protection setpoint methodology for the advanced power reactor 1400 (APP1400) and the optimized power reactor 1000 (OPR1000) has been recently developed to cover both the design basis event and the beyond design basis event. The developed setpoint methodology has also been quantitatively validated using specific computer programs and setpoint calculations. However, the safety of nuclear power plants cannot be fully guaranteed by satisfying the requirement of the analytical limit. In spite of the response time verification requirements of nuclear regulations and industry standards, it is hard to find the studies on the systematically integrated methodology regarding the response time evaluation. In cases of APR1400 and OPR1000, the response time analysis for the plant protection system is partially included in the setpoint calculation and the response time test is separately performed via the specific plant procedure. The test technique has a drawback which is the difficulty to demonstrate completeness of timing test. The analysis technique has also a demerit of resulting in extreme times that not actually possible. Thus, the establishment of the systematic response time evaluation methodology is needed to justify the conformance to the response time requirement used in the safety analysis. This paper proposes the response time evaluation methodology for APR1400 and OPR1000 using the combined analysis and test technique to confirm that the plant protection system can meet the analytical response time assumed in the safety analysis. In addition, the results of the quantitative evaluation performed for APR1400 and OPR1000 are presented in this paper. The proposed response time analysis technique consists of defining the response time requirement, determining the critical signal path for the trip parameter, allocating individual response time to each component on the signal path, and analyzing the total response time for the trip parameter, and demonstrates that the total analyzed response time does not exceed the response time requirement. The proposed response time test technique is composed of defining the response time requirement, determining the critical signal path for the trip parameter, determining the test method for each component on the signal path, performing the response time test, and demonstrates that the total test result does not exceed the response time requirement. The total response time should be tested in a single test that covers from the sensor to the final actuation device on the instrument channel. When the total channel is not tested in a single test, separate tests on groups of components or single components including the total instrument channel shall be combined to verify the total channel response. For APR1400 and OPR1000, the ramp test technique is used for the pressure and differential pressure transmitters and the step function testing technique is applied to the signal processing equipment and final actuation device. As a result, it can be demonstrated that the response time requirement is satisfied by the combined analysis and test technique. Therefore, the proposed methodology in this paper plays a crucial role in guaranteeing the safety of the nuclear power plants systematically satisfying one of two critical requirements from the safety analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of ORIGEN for VVER-440 Spent Fuel with Application to Fork Detector Safeguards Measurements

The US Department of Energy National Nuclear Security Administration and the European Atomic Energy Community (Euratom) are collaborating with the Radiation and Nuclear Safety Authority in Finland to assess spent fuel verification methods for potential application to measurements at the planned spent fuel encapsulation plant and geological repository prior to final disposal in Finland. The fork detector (FDET) used to measure the neutron and gamma radiation from a given spent fuel assembly is an existing technology widely used by the International Atomic Energy Agency (IAEA) and Euratom for the safeguards verification of spent fuel operator declaration. Recently, an FDET data analysis software, referred to as the ORIGEN module, was developed and incorporated into the IAEA/Euratom Integrated Review and Analysis Program. This module uses the ORIGEN burnup code to calculate the nuclide concentrations and the neutron and gamma source terms in the spent fuel based on operator declarations. Then, to predict the expected FDET signals for a given fuel assembly, the module combines the source terms calculated by ORIGEN with FDET response functions that were pre-generated using MCNP models. The inspector can compare the calculated signals to the measured values in real time to identify possible anomalies in fuel assembly’s operator declaration, integrity, or measurement. In this work, the capability of the ORIGEN module was extended from application to typical light water reactor fuels to use for VVER-440 assembly types. The accuracy of ORIGEN for this type of fuel assembly calculations was assessed by comparing calculated nuclide concentrations against destructive assay measurements for VVER-440 spent fuel. The performance of the ORIGEN module for FDET safeguards verification in routine inspections was assessed using FDET measurement data for 13 spent VVER-440 assemblies that were measured at the Loviisa Nuclear Power Plant in Finland, by comparing the calculated signals to the measured quantities. The results show that the calculated FDET neutron and gamma detector signals are generally within 12% of the measurements except for one assembly. Thus, these results are applicable to future safeguards verification in the planned Finland encapsulation plant, such as for passive neutron albedo reactivity safeguards measurements.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

The DD Neutron Generator as an Alternative to Am(Li) Isotopic Neutron Source in the Uranium Neutron Coincidence Collar

The 235 U mass assay of bulk uranium items, such as oxide canisters, fuel pellets, and fuel assemblies, is not achievable by traditional gamma-ray assay techniques because of the limited penetration of the item by the characteristic 235 U gamma rays. Instead, fast neutron interrogation methods such as active neutron coincidence counting must be used. For international safeguards applications, the most commonly used active neutron systems, the Active Well Coincidence Counter (AWCC), Uranium Neutron Coincidence Collar (UNCL), and 252 Cf Shuffler, rely on fast neutron interrogation using an isotopic neutron source [i.e., 252 Cf or Am(Li)] to achieve better measurement accuracies than are possible using gamma-ray techniques for high-mass, high-density items. However, the Am(Li) sources required for the AWCC and UNCL systems are no longer manufactured, and newly produced systems rely on limited supplies of sources salvaged from disused instruments. The 252 Cf shuffler systems rely on the use of high-output 252 Cf sources, which while still available have become extremely costly for use in routine operations and require replacement every 5–7 years. Lack of a suitable alternative neutron interrogation source would leave a potentially significant gap in the safeguarding of uranium processing facilities. Previously we examined the performance of Oak Ridge National Laboratory’s Large Volume Active Well Coincidence Counter modified to accept a commercially available deuterium-deuterium (DD) neutron generator to examine the potential of the DD neutron generator as an alternative to the isotopic sources [1]. That work indicated that the DD generator was a viable replacement for the AWCC applications, but because of the large dimensions of that generator, the DD/AWCC would have to be run as a delayed neutron counting system to provide equivalent measurement precision. In this work we examine the performance of a standard UNCL modified to accept the DD neutron generator for the assay of 235 U based on the results of Monte Carlo N-Particle (MCNP) simulations and measurements of depleted uranium and highly enriched uranium items. The projected performance of the modified UNCL for fresh pressurized water reactor fuel assemblies operated in both traditional coincidence counting and delayed neutron counting modes is presented.

07 ISOTOPE AND RADIATION SOURCES↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Light collection and pulse-shape discrimination in elongated scintillator cells for the PROSPECT reactor antineutrino experiment

A meter-long, 23-liter EJ-309 liquid scintillator detector has been constructed to study the light collection and pulse-shape discrimination performance of elongated scintillator cells for the PROSPECT reactor antineutrino experiment. The magnitude and uniformity of light collection and neutron-gamma discrimination power in the energy range of antineutrino inverse beta decay products have been studied using gamma and spontaneous fission calibration sources deployed along the cell axis. We also study neutron-gamma discrimination and light collection abilities for differing PMT and reflector configurations. As a result, key design features for optimizing MeV-scale response and background rejection capabilities are identified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Prediction of Adsorption Energies for Chemical Species on Metal Catalyst Surfaces Using Machine Learning

Computational catalyst screening has the potential to significantly accelerate heterogeneous catalyst discovery. Typically, this involves developing microkinetic reactor models that are based on parameters obtained from density functional theory and transition-state theory. To reduce the large computational cost involved in computing various adsorption and transition-state energies of all possible surface states on a large number of catalyst models, linear scaling relations for surface intermediates and transition states have been developed that only depend on a few, typically one or two descriptors, such as the carbon atom adsorption energy. As a result, only the descriptor values have to be computed for various active site models to generate volcano curves in activity or selectivity. Unfortunately, for more complex chemistries the predictability of linear scaling relations is unknown. Also, the selection of descriptors is essentially a trial and error process. Here, using a database of adsorption energies of the surface species involved in the decarboxylation and decarbonylation of propionic acid over eight monometalic transition-metal catalyst surfaces (Ni, Pt, Pd, Ru, Rh, Re, Cu, Ag), we tested if nonlinear machine learning (ML) models can outperform the linear scaling relations in prediction accuracy when predicting the adsorption energy for various species on a metal surface based on data from the rest of the metal surfaces. We found linear scaling relations to hold well for predictions across metals with a mean-absolute error of 0.12 eV, and ML methods being unable to outperform linear scaling relations when the training dataset contains a complete set of energies for all of the species on various metal surfaces. Only when the training dataset is incomplete, namely, contains a random subset of species’ energies for each metal, a currently unlikely scenario for catalyst screening, do kernel-based ML models significantly outperform linear scaling relations. We also found that simple coordinate-free species descriptors, such as bond counts, achieve as good results as sophisticated coordinate-based descriptors. Finally, we propose an approach for automatic discovery of appropriate metal descriptors using principal component analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Neutronics of Advanced-Fuel Fusion Experiments with DT Triggering

Nuclear fusion studies are currently devoted to the Deuterium-Tritium (DT) fuel cycle, the easiest way to reach ignition. The recent stress on safety by the world community has stimulated the research on other fuel cycles, based on advanced reactions, such as Deuterium-Helium-3 (DHe). With DHe, it is not necessary to breed and fuel tritium. The DHe cycle has a very low presence of energetic fusion neutrons, due to DD and DT side reactions, and it has the possibility to obtain electrical power by direct energy conversion of the protons produced by the fusion reactions. DHe fusion has its own set of problems, such as the availability of {sup 3}He and the attainment of the higher plasma parameters. To explore the possibilities of DHe plasmas, a burning plasma experiment at high magnetic field and high plasma densities is particularly attractive. Ignitor is a proposed compact high-magnetic field tokamak, aimed at studying plasma-burning conditions in Deuterium-Tritium plasmas. The Ignitor experiment is designed to reach ignition under controlled DT burning conditions. The machine parameters have been established, on the basis of existing knowledge of the confinement properties of high density plasmas. The scientific goal of the Ignitor experiment is to approach, for the first time, the ignition conditions of a magnetically confined D-T plasma. However, the plasma density limit in Ignitor being well above the optimal density for DT ignition, is potentially suitable to the higher densities required for DHe burning. In fact, Ignitor has been designed also to explore conditions where 14.7-MeV protons and 3.6-MeV alpha particles produced by the DHe reactions can supply a significant amount of thermal energy to a well-confined plasma. In particular, Ignitor can sustain plasma currents exceeding those required to confine proton orbits at birth, and has more than sufficiently high densities so that the slowing-down time of both the protons and the DT alpha particles is shorter than the electron energy replacement time of the thermal plasma in which they are produced. Preliminary analyses show that a fusion power PF around 2 MW, may be reached. As a start, Ignitor can perform initial studies at the level of approximately 1 MW of power in charged particles from the DHe reaction in a mostly DT plasma. CANDOR, a design evolution of Ignitor in the direction of a reactor using a DHe fuel cycle, has been proposed; it is a feasibility study of a high-field DHe experiment of larger dimensions and higher fusion power than Ignitor, however still based on the core Ignitor technologies. The main characteristics of the CANDOR machine are the following: the major radius Ro is about double than Ignitor, plasma currents up to 25 MA with toroidal magnetic fields B{sub T} = 13 T can be produced. Unlike Ignitor, CANDOR would operate with higher values of poloidal beta around unity and the central part of the plasma column in the Second Stability region. The characteristic times over which the plasma discharge can be sustained are longer by more than a factor of 4 than those of Ignitor. To produce such a strong magnetic field in a larger machine, the toroidal field coils have been divided into two sets of coils and the central solenoid (air core transformer) is placed between them in the inboard part. In CANDOR, the DHe burning regime can be reached, by a combination of ICRF heating and alpha particle heating due to DT fusion reactions, which take the role of a sort of trigger. Thanks to this fact, and unlike other proposed DHe fusion experiments, CANDOR is capable of reaching DHe ignition using the existing technologies, and present knowledge of plasma confinement phenomena. With the initial use of DT, where tritium comprises some 50% of an initial lower density DT plasma, but is not added thereafter, the need for an intense auxiliary heating, which is one of the main technological drawbacks of DHe ignition, would be considerably alleviated, becoming feasible with the present technology. However, this method has the disadvantage of using tritium and of presenting a higher neutron flux (due to DT reactions) than 'pure' DHe plasmas, generating a neutron flux transient when passing from the initial DT trigger reaction to the final DHe burning plasma. However, for an experimental machine one-of-a-kind of relatively small dimensions, this is not an issue. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Material Assessment for ITER's Collective Thomson Scattering first mirror

The International Thermonuclear Energy Reactor (ITER) Collective Thomson Scattering (CTS) system is a diagnostic instrument that measures plasma density and velocity through Thomson scattering of microwave radiation. Some of the key components of the CTS are quasi-optical mirrors that are used to produce astigmatic beam patterns, which have impact on the strength and spatial resolution of the diagnostic signal. The mirrors are exposed to neutron radiation, which may alter the quality of the signal received. In this work, three different materials (molybdenum (Mo), stainless steel 316 (SS-316) and tungsten (W)) are considered for the first mirror of the CTS. The objective is to access which of the material studied are best suited for this mirror, considering different neutron radiation loads simulated scenarios defined by ITER, based on the resultant stresses and temperature distributions. For it, the neutron irradiation, and subsequent heat-load on the mirrors are simulated using the Monte Carlo N-Particle (MCNP) code. Based on the MCNP heat-load results, transient thermal-structural Finite Element Analysis (FEA) of the mirror over a 400 s discharge, with and without cooling on the rear side, are conducted using in commercial FEA software ANSYS. Results show that of the tested materials Mo and W are the most suitable material for this application. Even though, this study does not yet consider the variation of the material properties with temperature, it presents a quick initial satisfactory assessment that may be considered in subsequent and more complex analysis. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ-ZrHx and ϵ-ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of θ 13 in Double Chooz using neutron captures on hydrogen with novel background rejection techniques

We observed a measurement of the Double Chooz collaboration and the neutrino mixing angle θ 13 using reactor $\bar{v}$ e via the inverse beta decay reaction in which the neutron is captured on hydrogen. Our measurement is based on 462.72 live days data, approximately twice as much data as in the previous such analysis, collected with a detector positioned at an average distance of 1050 m from two reactor cores. Several novel techniques have been developed to achieve significant reductions of the backgrounds and systematic uncertainties. Accidental coincidences, the dominant background in this analysis, are suppressed by more than an order of magnitude with respect to our previous publication by a multi-variate analysis. Furthermore, these improvements demonstrate the capability of precise measurement of reactor $\bar{v}$ e without gadolinium loading. Spectral distortions from the $\bar{v}$ e reactor flux predictions previously reported with the neutron capture on gadolinium events are confirmed in the independent data sample presented here. A value of sin 2 2θ 13 = 0.095 0.039 +0.038 (stat+syst) is obtained from a fit to the observed event rate as a function of the reactor power, a method insensitive to the energy spectrum shape. A simultaneous fit of the hydrogen capture events and of the gadolinium capture events yields a measurement of sin 2 2θ 13 = 0.088 ± 0.033(stat+syst).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Automated Resonance Fitting for Nuclear Data Evaluation

Global and national efforts to deliver high-quality nuclear data to users have a wide-ranging impact, affecting applications in national security, reactor operations, basic science, medicine, and more. Cross section evaluation is a major part of this effort, combining theory and experimentation to produce recommended values and uncertainties for reaction probabilities. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant manual effort and months of time from expert scientists. In this article, non-convex non-linear optimization methods are combined with concepts of inferential statistics to infer a resonance model from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and the potential for bias from prior assumptions, while enhancing reproducibility and documentation, thereby addressing well-known challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Correlated Production and Analog Transport of Fission Neutrons and Photons Using Fission Models FREYA, FIFRELIN, and the Monte Carlo Code TRIPOLI-4

Fission modeling in general-purpose Monte Carlo transport codes often relies on average nuclear data provided by the international evaluation libraries. As such, only average fission multiplicities are available and correlations between fission neutrons and photons are missing. Whereas uncorrelated fission physics is usually sufficient for standard reactor core and radiation shielding calculations, correlated fission secondaries are required for specialized nuclear instrumentation and detector modeling. For coincidence counting detector optimization for instance, the precise simulation of fission neutrons and photons that remain correlated in time from birth to detection is essential. New developments were recently integrated into the Monte Carlo transport code TRIPOLI-4 to model fission physics more precisely, the purpose being to access event-by-event fission events from two different fission models: Fission Reaction Event Yield Algorithm (FREYA) and Fission Fragment Evaporation Leading to an Investigation of Nuclear Data (FIFRELIN). TRIPOLI-4 simulations can now be performed, either by connecting via an application programming interface to the Lawrence Livermore National Laboratory fission library including FREYA, or by reading the external fission event data files produced by FIFRELIN beforehand. These new capabilities enable us to easily compare results from the Monte Carlo transport calculations using the two fission models in a nuclear instrumentation application. In the first part of this paper, the broad underlying principles of the two fission models are recalled. We then present the experimental measurements of neutron angular correlations for 252 Cf(sf) and 240 Pu(sf). The correlations were measured for several neutron kinetic energy thresholds. In the latter part of this paper, simulation results are compared with the experimental data. Spontaneous fissions in 252 Cf and 240 Pu are modeled by FREYA or FIFRELIN. Emitted neutrons and photons are subsequently transported to an array of scintillators by TRIPOLI-4 in the analog mode to preserve their correlations. Angular correlations between fission neutrons obtained independently from these TRIPOLI-4 simulations, using either FREYA or FIFRELIN, are compared with the experimental results. In conclusion, for 240 Pu(sf), the measured correlations were used to tune one of the FREYA model parameters. Similarly, different sets of parameters have been tested in FIFRELIN to try to improve the agreement with the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Evaluation of Changes to the MTR Fuel Dissolution Flowsheet to Increase the Number of L-Bundles Charged to a Dissolver

Material Test Reactor (MTR) fuels are dissolved in the H-Canyon facility using a Hg-catalyzed, HNO₃ dissolution flowsheet. The flowsheet is based on experimental work in which laboratory-scale experiments were used to provide bounding values for both offgas and H₂ generation rates. The generation of H₂ gas, primarily from the dissolution of Al, is controlled by limiting the surface area exposed to the solution. The number of L-Bundles charged and their submergence level defines the exposed surface area which is directly proportional to the H₂ generation rate. For the MTR fuel dissolution flowsheet, the maximum number of L-Bundles which can be charged to a dissolver was specified as a function of the Al concentration in the solution and the submergence level of the bundle. The MTR fuels are currently dissolved in the 6.1D dissolver. To reach the desired terminal Al concentration, a total of 12 L-Bundles of fuel are dissolved per batch. Depending on the submergence level being used in the dissolver, the 12 bundles are dissolved in three charges of four bundles per charge or in a sequence of three, four, and five bundles. However, if six L-Bundles could be dissolved in the initial charge, a dissolver batch could be completed in two dissolutions instead of three. To assess the feasibility of dissolving six bundles in the first charge, the Savannah River National Laboratory evaluated changes to the MTR flowsheet which increase the number of bundles which can be charged to a dissolver. Since the development of the existing MTR fuel dissolution flowsheet, improved offgas measurement methods have been developed to more accurately calculate the offgas composition and generation rate during laboratory-scale dissolutions. In the first part of this study, improvements in the calibration model for the Raman spectrometer (used for offgas characterization) and a more effective data smoothing routine were applied to existing data. The uncertainty in the offgas composition data used to define the MTR fuel dissolution flowsheet were subsequently re-evaluated which resulted in an approximate 50% reduction in the uncertainty of the H₂ gas concentration measurements. The reduced uncertainty was used in new calculations to define the maximum number of L-Bundles which can be charged to either the 6.1D or 6.4D dissolver as a function of the Al concentration. With the reduced uncertainty in the H₂ concentration measurements, a flowsheet was defined which supported initially charging 6 bundles to a dissolver. This flowsheet is based on the use of 0.002 M Hg, an L-Bundle immersion depth of 54 in., an air sparge/purge flowrate of 40 scfm, an iodine reactor operating temperature of 200 °C, and not exceeding 60% of the H₂ LFL in the dissolver offgas. In the second phase of this study, calculations were performed to assess the impact of changing various operating parameters to increase the number of L-Bundles charged to an H-Canyon dissolver. Applying the improved offgas measurement techniques including the reduced uncertainty in the H₂ concentration measurements to the existing data, the number of bundles of MTR fuel which can be charged to either the 6.1D or 6.4D dissolver were calculated based on varying the bundle immersion depth (44-60 in.), air sparge/purge flowrate (40-60 scfm), iodine reactor operating temperature (150-200 °C), and the percentage of the H₂ LFL which must not be exceeded (60-76%). Qualitatively, the percentage of the LFL had the largest influence on the maximum number of bundles which can be charged to a dissolver followed by the air sparge/purge flowrate, L-Bundle immersion depth, and the iodine reactor operating temperature. During the evaluation of each variable, the other variables were held at the baseline flowsheet conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Core-Shell Oxidative Aromatization Catalysts for Single Step Liquefaction of Distributed Shale Gas (Final Technical Report)

The objective of this project was to design and demonstrate a core-shell structured multifunctional catalyst to convert the light (dry) components of shale gas into liquid aromatic compounds (primarily benzene and toluene) in a single step. Operated in a modular oxidative aromatization system (OAS) under a cyclic redox scheme, the novel catalyst and process can significantly improve the value and transportability of distributed shale gas. Since the project started, each quarter addressed a different set of tasks related to the completion of the milestone detailed in the project award. The yearly summaries of these tasks are summarized below: Q1-Q4: • Conducted project planning and literature search. • Investigated a number of SHC redox catalysts using thermogravimetric analysis and fixed-bed reactor experiments. • Initiated process modeling towards generating two process models for the methane DHA base case and OAS process. • Developed DHA catalysts capable of producing >500 g/kg-cat-hr aromatics at 80% or greater aromatics selectivity at 700°C. Q5-Q8: • Developed alternative approaches with sequential bed configurations to enhance the aromatic yields based on OCM+DHA • Improved the zeolite synthesis efficiency by using the microwave-assisted technique and investigated the synthesis conditions on the zeolite yield, crystalline structure and morphology • Constructed a set of Aspen Plus process models with significant energy savings for OAS as compared to the base case non-oxidative DHA. • Adapted conventional hydrothermal method to be applicable to the microwave synthesizer unit for more efficient catalyst synthesis. • Studied the structure of the OCM catalyst and the dispersion of the carbonate in the redox reactions and in methane flow with Raman Spectroscopy. Q9-Q12: • Scaled up the catalyst synthesis with the microwave synthesis method. Based on its performance, procedural characterizations and catalytic performance testing were further conducted for the new microwave synthesized catalysts with the newly-developed product analysis procedure. • Developed the reaction system setup for the C2-DHA or OCM+DHA reaction product and achieved a better product collection-analysis method for the aromatic products with an improved carbon balance. The product from the OCM reaction exhibited complicated effects on the DHA catalyst. • Conducted additional OCM catalyst characterization using Near Ambient Pressure X-ray Photoelectron Spectroscopy and in situ Raman characterization • Validated the significant energy savings for OAS as compared to the base case non-oxidative DHA. Successfully set up the simulation model for the OCM+DHA+SHC reaction system based on the updated experimental results from NCSU. Q13-End of project: • Synthesized new zeolite catalysts by the microwave method, conducted characterizations (XRD, SEM, and TEM) and catalytic behavior testing. • Explored the “wet” C 2 H 6 and C 2 H 4 DHA reactions with using steam co-feed. A subsequent reduction as the regeneration step can regenerate the DHA catalyst and recover 99% activity of the fresh performance. • Achieved a 15.3% single-pass aromatic yield from methane by rationally combining the OCM and DHA at different temperatures. • Conducted a 105-hour stability test with an improved regeneration procedure, with an average aromatic yield of 13.8%. • Developed new catalyst and achieved a record-high 23.2% yield.

03 NATURAL GAS↗