Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear material accountability”

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 343 records · Page 19

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Residual Insoluble Solids Expected from Saltcake Dissolution

A review of previous SRS saltcake dissolution and characterization studies was performed to inform the analysis of the amount of residual insoluble solids that may result from saltcake dissolution. A realistic value covering most cases for the maximum volume of a settled layer of residual insoluble solids that would result from dissolution of a volume of saltcake was calculated as 0.08 mL/mL (8 vol%). Applying factors to account for dissolution test uncertainty and differences between in-tank settling and laboratory testing, the recommended value for a settled layer of residual insoluble solids that would result from dissolution of a volume of saltcake is 0.16 mL/mL (16 vol%). This should apply to dissolution of saltcake with typical minimal amounts of entrained sludge but may not apply to saltcake with significant sludge. As it is composed of primarily non-radioactive isotopes, the precipitated aluminum hydroxide material does not contribute to the heat load.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of PBR Spent Fuel Criticality and Dose Rate Compliance for Storage and Transportation

Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗

Non-nuclear Component Signatures for Warhead Dismantlement Confirmation

The verification of warhead dismantlement is expected to be an important component in future arms reduction treaties. Historic approaches developed with future arms control treaty verification in mind often involve intrusive measurements, process monitoring, and/or inspector presence to provide confidence that an authentic warhead has been dismantled. This work explores the possibility of reducing the negative impacts of these invasive approaches while also delivering a method that is more likely to provide non-sensitive data that can be shared with not only other nuclear weapons states but also non-nuclear weapons states partners. This work explores a novel approach for verifying dispositioned non-nuclear weapon components, providing confidence post-dismantlement that a treaty accountable item that was dismantled was in fact a treaty-relevant nuclear weapon system as declared. This method provides an alternative to intrusive inspection processes in nuclear weapons production environments, which would require significant changes to the host’s operational behaviors. It achieves this by identifying intrinsic neutron-induced signatures of non-nuclear components to determine their authenticity and estimate the duration they were exposed within a nuclear weapons system using technologies that are already in use for other national security applications. Intrinsic radiation effects studies are already a part of the stockpile aging and surveillance evaluations. However, none of these technologies and approaches have been previously considered for verification applications of non-nuclear component disposition. In this report, we introduce modeling studies that have been used to identify the most promising candidate parts and materials with signatures that are measurable and actionable. These models have been validated with laboratory measurements of signatures induced by the exposure of candidate materials to neutrons over a range of times. Predictive modeling then demonstrates the methodology for estimating exposure times and/or limits. Laboratory measurements of authentic non-nuclear parts from a dismantled warhead demonstrate the feasibility of employing these signature measurements. And finally, a concept of operations (CONOPS) for the potential use of this methodology is presented.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Criticality Control Overpack Fire Testing (Phase III)

This report will describe the one test conducted during phase III of the Pipe Overpack Container (POC) test campaign, present preliminary results from these tests, and discuss implications for the Criticality Control Overpack (CCO). The goal of this test was to see if aerosol surrogate material inside the Criticality Control Container (CCC) gets released when the drum lid of the CCO comes off during a thirty-minute long, fully-engulfing, fire test. As expected from POC tests conducted in Phase I and II of this test campaign, the CCO drum lid is ejected about one minute after the drum is exposed to fully-engulfing flames. The remaining pressure inside the drum is high enough to eject the top plywood dunnage a considerable distance from the drum. Subsequently, most of the bottom plywood dunnage supporting the CCC burns off during and after the fire. High pressure buildup inside the CCC and inside two primary containers holding the surrogate powder also results in damage to the filter media of the CCC and the filter-house, thread attachment of the primary canisters. No discernable release of surrogate powder material was detected from the two primary containers when pre- and post-test average mass were compared. However, when the average masses are corrected to account for possible uncertainties in mass measurements, error overlap does not preclude the possibility that some surrogate powder mass may have been lost from these primary canisters. Still, post-test conditions of the secondary canisters enclosing these two primary canisters suggest it is very unlikely this mass loss would have escaped into the CCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗

Passive and active neutron signatures of 233 U for nondestructive assay

The thorium fuel cycle is emerging as an attractive alternative to conventional nuclear fuel cycles, as it does not require the enrichment of uranium for long-term sustainability. The operating principle of this fuel cycle is the irradiation of 232 Th to produce 233 U, which is fissile and sustains the fission chain reaction. 233 U poses unique challenges for nuclear safeguards, as it is associated with a uniquely extreme γ-ray environment from 232 U contamination, which limits the feasibility of the γ-ray-based assay, as well as more conservative accountability requirements than for 235 U set by the International Atomic Energy Agency. Consequently, instrumentation used for safeguarding 235 U in traditional fuel cycles may be inapplicable. It is essential that the nondestructive signatures of 233 U be characterized so that nuclear safeguards can be applied to thorium fuel-cycle facilities as they come online. In this work, a set of 233 U 3 O 8 plates, containing 984 g 233 U, was measured at the National Criticality Experiments Research Center. A high-pressure 4 He gaseous scintillation detector, which is insensitive to γ-rays, was used to perform a passive fast neutron spectral signature measurement of 233 U 3 O 8 , and was used in conjunction with a pulsed deuterium-tritium neutron generator to demonstrate the differential die-away signature of this material. Furthermore, an array of 3 He detectors was used in conjunction with the same neutron generator to measure the delayed neutron time profile of 233 U, which is unique to this nuclide. These measurements provide a benchmark for future nondestructive assay instrumentation development, and demonstrate a set of key neutron signatures to be leveraged for nuclear safeguards in the thorium fuel cycle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear quantum effects on the quasiparticle properties of the chloride anion aqueous solution within the GW approximation

Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral abinitio molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared with the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the Cl – anion and water. Furthermore, our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lifetime estimation of a BWR core shroud in terms of IGSCC

The continued operation of aging Boiling Water Reactors (BWRs) worldwide requires gradually increasing cost of inspection, maintenance, and repair. Intergranular Stress Corrosion Cracking (IGSCC) in sensitized austenitic stainless steel piping first became a major issue for BWRs in the 1980s, resulting in recognition of the susceptibility of reactor internals to IGSCC. Shroud cracking identified in 1993–1994 confirmed that IGSCC of internals is a significant issue for BWRs. IGSCC is a time-dependent, material degradation process, which is caused and accelerated by the presence of residual stresses, material sensitization, irradiation, cold work, elevated temperature, and corrosive environments. This paper emphasizes the importance of accounting for corrosive environments, or more exactly, electrochemical phenomena in modelling IGSCC and predicting the service life of BWR in-vessel components, and stresses the necessity of performing such modeling not just for a single state point under full power conditions, but for the whole operating history of the reactor, including startups and shutdowns. Overall, this paper demonstrates that ignoring electrochemical considerations may result in underestimating component lifetime and lead to unnecessary expenses for inspection and repair.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Rapid Sequential Separation Method for Determination of Actinides and TENORM in Fracking Wastes - 20223

With recent advances in unconventional drilling technology in the US and around the world, the risks of environmental contamination and exposure due to technologically enhanced naturally occurring radioactive materials (TENORM) from fracking wastes have also surfaced. Although the Permian Basin (southeastern New Mexico, western Texas) has long been a modest producer of oil, the advent of fracking technology a decade ago turned it into hot property for producers. A 2018 assessment of undiscovered, technically recoverable continuous oil and gas resources by the U.S. Geological Survey estimates an average 46.3 billion barrels of oil and 281 trillion cubic feet of gas. This is likely to increase the use of unconventional drilling and TENORM generation. TENORM concentrations in oil and gas exploration and production waste can be several hundred to several thousand times higher than background TENORM concentrations. Both the thorium and uranium decay chains contribute to airborne radionuclides arising from the radon gas escaping the ground and subsequently decaying as airborne particulates. Additionally, radon decay products {sup 210}Pb and {sup 210}Po can build up in scale on the internal surfaces of oil and gas handling pipes and in sludge in refineries, becoming potential inhalation and ingestion hazards for workers. Furthermore, southeastern NM is also home to world's only licensed and operating transuranic nuclear waste repository, the DoE's Waste Isolation Pilot Plant (WIPP). Plutonium isotopes ({sup 239+240}Pu) and {sup 241}Am, are expected to account for more than 99% of the total radioactivity scheduled for disposal in the WIPP repository. Thus, an improved understanding of the environmental fate and transport of TENORM, liberated by unconventional drilling, is essential to assess how best to protect individuals and the environment. In this context, accurate measurement of TENORM and actinides in environmental and biological samples is essential. In this presentation, a new sequential method for the separation and pre-concentration of actinides (Pu, Am, Np) and TENORM (Po, U and Th) derived from oil and gas exploration is proposed. The TEVA method involves a rapid co-precipitation step to remove matrix interferences from the samples, followed by plutonium oxidation state adjustment to Pu (IV) and an incubation period of ∼ 1 hour at 50-60 deg. C to allow the resultant Po (II) to oxidize into Po(IV). The polonium, neptunium, thorium and plutonium are then separated on a TEVA column, while americium and uranium are separated on a TRU column. After separation, the alpha counting source was prepared by micro-precipitation with copper sulfide (CuS) for polonium and neodymium fluoride (NdF{sub 3}) for actinides. The efficiency and reliability of the procedures were tested by analyzing filter, drinking water and frack sand samples. The method is simple, robust and can be performed quickly with excellent removal of interferences, high chemical recovery and very good alpha peak resolution. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Thermodynamic non-ideality and disorder heterogeneity in actinide silicate solid solutions

Non-ideal thermodynamics of solid solutions can greatly impact materials degradation behavior. We have investigated an actinide silicate solid solution system (USiO 4 –ThSiO 4 ), demonstrating that thermodynamic non-ideality follows a distinctive, atomic-scale disordering process, which is usually considered as a random distribution. Neutron total scattering implemented by pair distribution function analysis confirmed a random distribution model for U and Th in first three coordination shells; however, a machine-learning algorithm suggested heterogeneous U and Th clusters at nanoscale (~2 nm). The local disorder and nanosized heterogeneous is an example of the non-ideality of mixing that has an electronic origin. Partial covalency from the U/Th 5f–O 2p hybridization promotes electron transfer during mixing and leads to local polyhedral distortions. The electronic origin accounts for the strong non-ideality in thermodynamic parameters that extends the stability field of the actinide silicates in nature and under typical nuclear waste repository conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear propagating modes beyond the phonons in fluorite-structured crystals

The vibrational energy of crystals is known to propagate in a fixed number of phonon branches allowed by symmetry. In the realm of nonlinear dynamics, however, additional nonlinear propagating modes are possible. Nonlinear propagating modes have unique properties that are important in many disciplines including optical communications, conducting polymers, biology, magnetism, and nuclear physics. Yet, despite the crucial importance of crystal lattice vibrations in fundamental and applied science, such additional propagating modes have not been observed in ordinary crystals. Here, we show that propagating modes exist beyond the phonons in fluorite-structured thoria, urania, and natural calcium fluoride using neutron scattering and first-principles calculations. These modes are observed at temperatures ranging from 5 K up to 1200 K, extend to frequencies 30–40% higher than the maximum phonon frequency, and travel at velocities comparable to or higher than the fastest phonon. The nonlinear origin of the modes is explained in part via perturbation theory, which approximately accounts for nonlinearity. Given that these modes are still clearly observed at 5 K, they are likely an inherent feature of the quantum ground state. The existence of these waves in three-dimensional crystals may have ramifications for material properties.

36 MATERIALS SCIENCE↗

Molten Salt Reactor Experiment Benchmark Evaluation

As one of the six advanced reactor concepts selected by the 2002 Generation IV roadmap as a technology meriting future research, Molten Salt Reactors (MSRs) attracted broad attention and multiple private and public entities are working towards its commercialization. Facing stringent regulatory requirements, validations of computational codes used to calculate and prove the safety characteristics of MSRs play a key role. This project developed the world-first, MSR-related reactor physics benchmark based on a series of zero-power experiments conducted at the Molten Salt Reactor Experiment (MSRE). The benchmark was successfully reviewed by the International Reactor Physics Experiment Evaluation Project (IRPhEP) committee and is available in the IRPhEP handbook starting form the 2019 edition. To benchmark experimental data from the MSRE, a three-dimensional high-fidelity model was developed using the Monte Carlo neutron transport code Serpent 2 and ad-hoc methods have been implemented to account for the unique feature of fuel salt motion. The calculated effective multiplication factor, k eff , for the first criticality experiment, when 235 U was progressively added to the fuel salt in order to achieve criticality with stationary salt and isothermal conditions, was 1.02132 (±3 pcm). The total uncertainty for experimental k eff was estimated to be 420 pcm. The calculated k eff is 2.154% larger than the experimental and benchmark model value, which is approximately five times the benchmark model uncertainty. It is to be noted that, for systems containing large volume of graphite (or other carbonaceous materials), Monte Carlo codes tend to overestimate the keff of the benchmark model by 1% to 2%. The bias is possibly due to uncertainties in the impurity content of the graphite blocks, uncertainty in the neutron capture cross section of carbon or uncertainty in the properties of nuclear-grade graphite. The calculated reactivity coefficient of 235 U concentration (0.2228 ± 0.0014, represented as the change of reactivity over the relative change of 235 U mass in loop) matches well with the experimental value (0.223 ± 0.006), strengthening the confidence on the accurate representation of the fuel salt composition in the MSRE benchmark. Additional reactivity effect calculations, including the control rod bank worth, reactivity effects of fuel circulation and isothermal and fuel temperature coefficients show good agreement with the experimental values (within 1σ) as well.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗