A vital part of the licensing process for advanced (non-LWR) nuclear reactor developers in the United States is the assessment of the reactor’s source term, i.e., the potential release of radionuclides from the reactor system to the environment during normal operations and accident sequences. In comparison to source term assessments which follow a bounding approach with conservative assumptions, a mechanistic approach to modeling radionuclide transport, which realistically accounts for transport and retention phenomena, is expected to be used for advanced reactor systems. As the designs of advanced reactors increase in maturity and progress towards licensing, there is a need to advance modeling and simulation capabilities in analyzing the mechanistic source term (MST) of a prospective reactor concept. In the present work, a survey is provided of existing computational capabilities for the modeling of advanced reactors MSTs. The following reactors are considered: high temperature gas reactors (HTGR); molten salt reactors (MSR) which include salt-fueled reactors and fluoride salt-cooled high temperature reactors (FHR); and sodium- and lead-cooled fast reactors (SFR, LFR). A review of relevant codes which may be useful in providing information to MST analyses is also completed, including codes that have been used for source term analyses of LWRs, as well as those being developed for other aspects of advanced reactor system modeling such as reactor physics, thermal hydraulics, and chemistry. A discussion of MST modeling capabilities for each reactor type is provided with additional focus on important phenomena and functional requirements. Additionally, a comprehensive survey is provided of tools for consequence modeling such as atmospheric transport and dispersion (ATD).
Graphene plasmons provide a compelling avenue toward chip-scale dynamic tuning of infrared light. Dynamic tunability emerges through controlled alterations in the optical properties of the system defining graphene's plasmonic dispersion. Typically, electrostatic induced alterations of the carrier concentration in graphene working in conjunction with mobility have been considered the primary factors dictating plasmonic tunability. We find here that the surrounding dielectric environment also plays a primary role, dictating not just the energy of the graphene plasmon but so too the magnitude of its tuning and spectral width. To arrive at this conclusion, poles in the imaginary component of the reflection coefficient are used to efficiently survey the effect of the surrounding dielectric on the tuning of the graphene plasmon. By investigating several common polar materials, optical phonons (i.e., the Reststrahlen band) of the dielectric substrate are shown to appreciably affect not only the plasmon's spectral location but its tunability, and its resonance shape as well. In particular, tunability is maximized when the resonances are spectrally distant from the Reststrahlen band, whereas sharp resonances (i.e., high-Q) are achievable at the band's edge. These observations both underscore the necessity of viewing the dielectric environment in aggregate when considering the plasmonic response derived from two-dimensional materials and provide heuristics to design dynamically tunable graphene-based infrared devices.
The final review for the FY20 Advanced Simulation and Computing (ASC) Integrated Codes (IC) L2 Milestone #7181 was conducted on August 31, 2020 at Sandia National Laboratories in Albuquerque, New Mexico. The review panel unanimously agreed that the milestone has been successfully completed. Roshan Quadros (1543) led the milestone team and various members from the team presented the results. The review panel was comprised of staff from Sandia National Laboratories Albuquerque and California that are involved with computational engineering modeling and analysis. The panel consisted of experts in the fields of solid modeling, discretization, meshing, simulation workflows, and computational analysis including personnel Brett Clark (1543, Chair); Jay Foulk (8363); Jackie Moore (1553); Ron Kensek (1341); Ed Hoffman (8753); Dan Ibanez (1443). The presentation documented the technical approach of the team and summarized the results with sufficient detail to demonstrate both the value and the completion of the milestone. A separate SAND report was also generated with more detail to supplement the presentation. The purpose of the milestone was to advance capabilities for automatically finding, displaying, and resolving geometric overlaps in CAD models.
Residual stress measurements using neutron diffraction and the contour method were performed on a valve housing made from 316 L stainless steel powder with intricate three-dimensional internal features using laser powder-bed fusion additive manufacturing. The measurements captured the evolution of the residual stress fields from a state where the valve housing was attached to the base plate to a state where the housing was cut free from the base plate. Making use of this cut, thus making it a non-destructive measurement in this application, the contour method mapped the residual stress component normal to the cut plane (this stress field is completely relieved by cutting) over the whole cut plane, as well as the change in all stresses in the entire housing due to the cut. The non-destructive nature of the neutron diffraction measurements enabled measurements of residual stress at various points in the build prior to cutting and again after cutting. Good agreement was observed between the two measurement techniques, which showed large, tensile build-direction residual stresses in the outer regions of the housing. The contour results showed large changes in multiple stress components upon removal of the build from the base plate in two distinct regions: near the plane where the build was cut free from the base plate and near the internal features that act as stress concentrators. These observations should be useful in understanding the driving mechanisms for builds cracking near the base plate and to identify regions of concern for structural integrity. Neutron diffraction measurements were also used to show the shear stresses near the base plate were significantly lower than normal stresses, an important assumption for the contour method because of the asymmetric cut.
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials' design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.
The term “in situ processing” has evolved over the last decade to mean both a specific strategy for visualizing and analyzing data and an umbrella term for a processing paradigm. The resulting confusion makes it difficult for visualization and analysis scientists to communicate with each other and with their stakeholders. To address this problem, a group of over 50 experts convened with the goal of standardizing terminology. This paper summarizes their findings and proposes a new terminology for describing in situ systems. An important finding from this group was that in situ systems are best described via multiple, distinct axes: integration type, proximity, access, division of execution, operation controls, and output type. Here, they discuss these axes, evaluate existing systems within the axes, and explore how currently used terms relate to the axes.
Determining a process–structure–property relationship is the holy grail of materials science, where both computational prediction in the forward direction and materials design in the inverse direction are essential. Problems in materials design are often considered in the context of process–property linkage by bypassing the materials structure, or in the context of structure–property linkage as in microstructure-sensitive design problems. However, there is a lack of research effort in studying materials design problems in the context of process–structure linkage, which has a great implication in reverse engineering. In this work, given a target microstructure, we propose an active learning high-throughput microstructure calibration framework to derive a set of processing parameters, which can produce an optimal microstructure that is statistically equivalent to the target microstructure. The proposed framework is formulated as a noisy multi-objective optimization problem, where each objective function measures a deterministic or statistical difference of the same microstructure descriptor between a candidate microstructure and a target microstructure. Furthermore, to significantly reduce the physical waiting wall-time, we enable the high-throughput feature of the microstructure calibration framework by adopting an asynchronously parallel Bayesian optimization by exploiting high-performance computing resources. Case studies in additive manufacturing and grain growth are used to demonstrate the applicability of the proposed framework, where kinetic Monte Carlo (kMC) simulation is used as a forward predictive model, such that for a given target microstructure, the target processing parameters that produced this microstructure are successfully recovered.
We report the atomic- and nanosecond-scale quantification of kinetics of a shock-driven phase transition in Zr metal. We uniquely make use of a multiple shock-and-release loading pathway to shock Zr into the β phase and to create a quasisteady pressure and temperature state shortly after. Coupling shock loading with in situ time-resolved synchrotron x-ray diffraction, we probe the structural transformation of Zr in the steady state. Our results provide a quantified expression of kinetics of formation of β-Zr phase under shock loading: transition incubation time, completion time, and crystallization rate.
