Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

We show how to sample uniformly within the three-sided region bounded by a circle, a radial ray, and a tangent, called a “chock.” By dividing a 2D planar rectangle into a background grid, and subtracting Poisson disks from grid squares, we are able to represent the available region for samples exactly using triangles and chocks. Uniform random samples are generated from chock areas precisely without rejection sampling. This provides the first implemented algorithm for precise maximal Poisson-disk sampling in deterministic linear time. We prove O(n · M(b) log b), where n is the number of samples, b is the bits of numerical precision and M is the cost of multiplication. Prior methods have higher time complexity, take expected time, are non-maximal, and/or are not Poisson-disk distributions in the most precise mathematical sense. We fill this theoretical lacuna.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Process parameter selection in laser powder bed fusion (LPBF) controls the as-printed dimensional tolerances, pore formation, surface quality and microstructure of printed metallic structures. Measuring the stochastic mechanical performance for a wide range of process parameters is cumbersome both in time and cost. In this study, we overcome these hurdles by using high-throughput tensile (HTT) testing of over 250 dogbone samples to examine process-driven performance of strut-like small features, ~1 mm2 in austenitic stainless steel (316 L). The output mechanical properties, porosity, surface roughness and dimensional accuracy were mapped across the printable range of laser powers and scan speeds using a continuous wave laser LPBF machine. Tradeoffs between ductility and strength are shown across the process space and their implications are discussed. While volumetric energy density deposited onto a substrate to create a melt-pool can be a useful metric for determining bulk properties, it was not found to directly correlate with output small feature performance.

Abstract not provided.

Abstract not provided.

Fault detection and classification in photovoltaic (PV) systems through real-time monitoring is a fundamental task that ensures quality of operation and significantly improves the performance and reliability of operating systems. Different statistical and comparative approaches have already been proposed in the literature for fault detection; however, accurate classification of fault and loss incidents based on PV performance time series remains a key challenge. Failure diagnosis and trend-based performance loss routines were developed in this work for detecting PV underperformance and accurately identifying the different fault types and loss mechanisms. The proposed routines focus mainly on the differentiation of failures (e.g., inverter faults) from irreversible (e.g., degradation) and reversible (e.g., snow and soiling) performance loss factors based on statistical analysis. The proposed routines were benchmarked using historical inverter data obtained from a 1.8 MWp PV power plant. The results demonstrated the effectiveness of the routines for detecting failures and loss mechanisms and the capability of the pipeline for distinguishing underperformance issues using anomaly detection and change-point (CP) models. Finally, a CP model was used to extract significant changes in time series data, to detect soiling and cleaning events and to estimate both the performance loss and degradation rates of fielded PV systems.

Abstract not provided.

The Advanced Reactor Safeguards (ARS) program was established in 2020 as part of appropriations for the Advanced Reactor Demonstration Program (ARDP) through the Office of Nuclear Energy in the Department of Energy. The goal of this program is to help address near term challenges that advanced nuclear reactor vendors face in meeting domestic Material Control and Accountancy (MC&A) and Physical Protection System (PPS) requirements for U.S. construction. The technical work in the program is meant to (1) support nuclear reactor vendors with advanced MC&A and PPS designs for next generation reactors, (2) provide technical bases for the regulator, and (3) promote the integration of Safeguards and Security by Design early in the design process. Existing domestic regulations for safeguards and security, as outlined in the Code of Federal Regulations, were written for large light water reactors, and rule-making efforts are underway to develop regulations more suited to different reactor designs. The ARS program seeks to remove roadblocks in the deployment of new and advanced reactors by solving regulatory challenges, reducing safeguards and security costs, and utilizing the latest technologies and approaches for robust plant monitoring and protection. This roadmap discusses the goals of the ARS program, current research, and program plan for the next five years.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Crossed-field diodes (CFDs) are used in multiple high-power applications and are characterized by an applied magnetic field orthogonal to the electric field, induced by the applied voltage across the anode-cathode gap. In vacuum, the Hull cutoff magnetic field (HCMF) represents the maximum applied magnetic field for which an electron from the cathode can reach the anode. This study investigates the effects of non-vacuum conditions on electron trajectories by introducing electron mobility, which represents particle collisions. We used numerical solutions of the electron force law and particle-in-cell simulations (XPDP1) to assess electron motion for various electron mobilities. For magnetic fields above the HCMF in vacuum, reducing the electron mobility increases the time for an electron emitted from the cathode to reach the anode. Reducing mobility below 22 C s/kg eliminates the HCMF for any magnetic field, meaning that an emitted electron will always cross the gap. We derived the magnetic field, mobility, and electron transit time corresponding to this condition by solving for the condition when the electron velocity in the direction across the anode-cathode gap going to zero at the anode. A parametric study of these conditions using theory and XPDP1 is performed under different gap distances, voltages, and magnetic fields.

Abstract not provided.

We propose developing numerical methods informed by novel experimental diagnostics that transition from solid-to-fluid, while accurately predicting the stress and deformation regardless of phase.

This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-22SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), numerical model comparison (DECOVALEX) and an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process. Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations. The work summarized in this annual update has occurred during the COVID-19 pandemic, and little international or domestic travel has occurred. Most of the collaborations have been conducted via email or as virtual meetings, but a slow return to travel and in-person meetings has begun.

This report documents the updated seismic shake table test plan. The report describes the shake table inputs (ground motions), test hardware, shake table facility, friction experiment, and proposed instrumentation.

CMOS-based computing systems that employ the von Neumann architecture are relatively limited when it comes to parallel data storage and processing. In contrast, the human brain is a living computational signal processing unit that operates with extreme parallelism and energy efficiency. Although numerous neuromorphic electronic devices have emerged in the last decade, most of them are rigid or contain materials that are toxic to biological systems. In this work, we report on biocompatible bilayer graphene-based artificial synaptic transistors (BLAST) capable of mimicking synaptic behavior. The BLAST devices leverage a dry ion-selective membrane, enabling long-term potentiation, with ~50 aJ/µm² switching energy efficiency, at least an order of magnitude lower than previous reports on two-dimensional material-based artificial synapses. The devices show unique metaplasticity, a useful feature for generalizable deep neural networks, and we demonstrate that metaplastic BLASTs outperform ideal linear synapses in classic image classification tasks. With switching energy well below the 1 fJ energy estimated per biological synapse, the proposed devices are powerful candidates for bio-interfaced online learning, bridging the gap between artificial and biological neural networks.

This SAND Report Guide offers support to authors, technical writers, principal investigators, and others involved in the process of creating, formatting, or refining a SAND Report. It details what you need to know before you begin compiling a SAND Report, directs you to the SAND Report templates, outlines the order of elements in a SAND Report, and explains what to do when your report is completed and ready for Review and Approval and subsequent distribution. Supporting information is provided in the appendix, such as where to get technical assistance, trademarks, Microsoft Word, and equations.

It has been demonstrated that grid cells in the brain are encoding physical locations using hexagonally spaced, periodic phase-space representations. We explore how such a representation may be computationally advantageous for related engineering applications. Theories of how the brain decodes from a phase-space representation have been developed based on neuroscience data. However, theories of how sensory information is encoded into this phase space are less certain. Here we show a method for how a navigation-relevant input space such as elevation trajectories may be mapped into a phase-space coordinate system that can be decoded using previously developed theories. We also consider how such an algorithm may then also be mapped onto neuromrophic systems. Just as animals can tell where they are in a local region based on where they have been, our encoding algorithm enables the localization to a position in space by integrating measurements from a trajectory over a map. In this paper, we walk through our approach with simulations using a digital elevation model.

It has been demonstrated that grid cells in the brain are encoding physical locations using hexagonally spaced, periodic phase-space representations. We explore how such a representation may be computationally advantageous for related engineering applications. Theories of how the brain decodes from a phase-space representation have been developed based on neuroscience data. However, theories of how sensory information is encoded into this phase space are less certain. Here we show a method for how a navigation-relevant input space such as elevation trajectories may be mapped into a phase-space coordinate system that can be decoded using previously developed theories. We also consider how such an algorithm may then also be mapped onto neuromrophic systems. Just as animals can tell where they are in a local region based on where they have been, our encoding algorithm enables the localization to a position in space by integrating measurements from a trajectory over a map. In this paper, we walk through our approach with simulations using a digital elevation model.

We are currently witnessing the dawn of hydrogen (H₂) economy, where H₂ will soon become a primary fuel for heating, transportation, and longdistance and long-term energy storage. Among diverse possibilities, H₂ can be stored as a pressurized gas, a cryogenic liquid, or a solid fuel via adsorption onto porous materials. Metal–organic frameworks (MOFs) have emerged as adsorbent materials with the highest theoretical H₂ storage densities on both a volumetric and gravimetric basis. However, a critical bottleneck for the use of H₂ as a transportation fuel has been the lack of densification methods capable of shaping MOFs into practical formulations while maintaining their adsorptive performance. Here, we report a high-throughput screening and deep analysis of a database of MOFs to find optimal materials, followed by the synthesis, characterization, and performance evaluation of an optimal monolithic MOF (monoMOF) for H₂ storage. After densification, this _monoMOF stores 46 g L^–1 H₂ at 50 bar and 77 K and delivers 41 and 42 g L^–1 H₂ at operating pressures of 25 and 50 bar, respectively, when deployed in a combined temperature– pressure (25–50 bar/77 K → 5 bar/160 K) swing gas delivery system. This performance represents up to an 80% reduction in the operating pressure requirements for delivering H₂ gas when compared with benchmark materials and an 83% reduction compared to compressed H₂ gas. Our findings represent a substantial step forward in the application of high-density materials for volumetric H₂ storage applications.

Dirac semimetals have attracted a great deal of current interests due to their potential applications in topological quantum computing, low-energy electronic devices, and single photon detection in the microwave frequency range. Herein are results from analyzing the low magnetic (B) field weak-antilocalization behaviors in a Dirac semimetal Cd₃As₂ thin flake device. At high temperatures, the phase coherence length lΦ first increases with decreasing temperature (T) and follows a power law dependence of l_Φ ∝ T^–0.4. Below ~3 K, l_Φ tends to saturate to a value of ~180 nm. Another fitting parameter α, which is associated with independent transport channels, displays a logarithmic temperature dependence for T > 3 K, but also tends to saturate below ~3 K. The saturation value, ~1.45, is very close to 1.5, indicating three independent electron transport channels, which we interpret as due to decoupling of both the top and bottom surfaces as well as the bulk. This result, to our knowledge, provides first evidence that the surfaces and bulk states can become decoupled in electronic transport in Dirac semimetal Cd₃As₂.

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager-Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the "ground-truth" molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ38from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

The selective amorphization of SiGe in Si/SiGe nanostructures via a 1 MeV Si+ implant was investigated, resulting in single-crystal Si nanowires (NWs) and quantum dots (QDs) encapsulated in amorphous SiGe fins and pillars, respectively. The Si NWs and QDs are formed during high-temperature dry oxidation of single-crystal Si/SiGe heterostructure fins and pillars, during which Ge diffuses along the nanostructure sidewalls and encapsulates the Si layers. The fins and pillars were then subjected to a 3 × 1015 ions/cm2 1 MeV Si+ implant, resulting in the amorphization of SiGe, while leaving the encapsulated Si crystalline for larger, 65-nm wide NWs and QDs. Interestingly, the 26-nm diameter Si QDs amorphize, while the 28-nm wide NWs remain crystalline during the same high energy ion implant. This result suggests that the Si/SiGe pillars have a lower threshold for Si-induced amorphization compared to their Si/SiGe fin counterparts. However, Monte Carlo simulations of ion implantation into the Si/SiGe nanostructures reveal similar predicted levels of displacements per cm3. Molecular dynamics simulations suggest that the total stress magnitude in Si QDs encapsulated in crystalline SiGe is higher than the total stress magnitude in Si NWs, which may lead to greater crystalline instability in the QDs during ion implant. The potential lower amorphization threshold of QDs compared to NWs is of special importance to applications that require robust QD devices in a variety of radiation environments.

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1-6 nm and pH ∼2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump-probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Nanothermite NiO-Al is a promising material system for low gas emission heat sources; yet, its reactive properties are highly dependent on material processing conditions. In the current study, sputter deposition is used to fabricate highly controlled nanolaminates comprised of alternating NiO and Al layers. Films having an overall stoichiometry of 2Al to 3NiO were produced with different bilayer thicknesses to investigate how ignition and self-sustained, high temperature reactions vary with changes to nanometer-scale periodicity and preheat conditions. Ignition studies were carried out with both hot plate and laser irradiation and compared to slow heating studies in hot-stage x-ray diffraction. Ignition behavior has bilayer thickness and heating rate dependencies. The 2Al/3NiO with λ ≤ 300 nm ignited via solid/solid diffusion mixing (activation energy, Ea = 49 ± 3 kJ/mole). Multilayers having λ≥ 500 nm required a more favorable mixing kinetics of solid/liquid dissolution into molten Al (Ea = 30 ± 4 kJ/mole). This solid/liquid dissolution Ea is a factor of 5 lower than that of the previously reported powder compacts due to the elimination of a passivating Al oxide layer present on the powder. The reactant mixing mechanism between 300 and 500 nm bilayer thicknesses was dependent on the ignition source's heating rate. The self-propagating reaction velocities of 2Al/3NiO multilayers varied from 0.4 to 2.5 m/s. Pre-heating nanolaminates to temperatures below the onset reaction temperatures associated with forming intermediate nickel aluminides at multilayer interfaces led to increased propagation velocities, whereas pre-heating samples above the onset temperatures inhibited subsequent attempts at laser ignition.

A study was conducted to investigate the outgassing characteristics of four thermal filler materials. The purpose of this screening was to identify any outgassing products that might be considered reactive, specifically compounds that could result in corrosion in the systems where these materials are used. A range of compounds was observed in the sample headspaces, though most do not stand out as being known reactive species of concern. However, several halogenated compounds and sulfurous compounds- classes compounds known to facilitate corrosion reactions under certain conditions- were observed in low concentrations. The TFLEX 760 exhibited the highest total outgassing, while the GR130 had the lowest. Therm-a-Gap75 and the Si thermal grease exhibited very similar outgassing profiles. It is difficult to predict the extent to which any given compound observed in an analysis of this type might pose a risk in an actual system; factors such as temperature, system geometry, concentration, and gas conductance all play a role in the kinetics governing chemical reactions. It is recommended that the results of these analyses are shared with pertinent materials SMEs familiar with the system(s) in question to evaluate potential risks.

This report deals with the development and evaluation of a numerical model to examine applied coating to a metal substrate subjected to a ballistic impact. The numerical model will be used to examine the benefit of the coating in resisting penetration due to the impact. For a detailed examination the Retch-Ipson curve is used as a metric. The numerical data is plotted and then fit to the Retch-Ipson curve and error calculations are used to compare the difference between the numerical output and the experimental data. This initial study is an examination of a few shortcomings of the standard material models used, and demonstrate the future work that is needed to understand the ballistic behavior of materials.

Next generation ion traps will likely need to support tens if not hundreds of ions in order to achieve several logical qubits. As we scale to those sizes, the same problems we face now – rf dissipation, control I/O, and optical access – will only grow and become more complicated. While many of these challenges can potentially be solved with technology integration, independently researching the feasibility of that integration and other solutions may help reduce the time and risk in scaling up to larger traps, by testing on smaller less complex devices. We should also consider other fabrication techniques that may help scale to larger devices, such as: through-substrate-vias (TSVs), different metal coatings, exotic rf routing, on chip laser sources, or even a secondary macroscopic trap to reload ions from. To have these technologies ready for full scale integration when we need them, ion traps with some of these capabilities need to be produced now. Developing the rigorous fabrication methods for producing reliable traps takes time and experimentation. We propose developing larger ion traps and reliable integrated technology in conjunction to make both available faster.

This summary report, prepared by Condor Country Consulting Inc. (CCCI), covers all construction monitoring of the Swainson’s hawk (Buteo swainsoni, SWHA) nest that occurred from April 2022 to July 2022 at the Sandia National Laboratories, California (SNL/CA). SNL/CA personnel first observed the SWHA pair near the pine tree along East Avenue that contained last year’s nest on April 19th, 2022. SWHA are a threatened species in the state of California and require a quarter-mile non-disturbance buffer. SNL/CA received approval from the California Department of Fish and Wildlife (CDFW) to continue work within the buffer area with a biological monitor on site to observe the SWHA for construction related disturbances. Monitoring began on April 20th, 2022 and concluded on July 1st, 2022.

Large-scale CO2 sequestration into geological formations has been suggested to reduce CO2 emissions from industrial activities. However, much like enhanced geothermal stimulation and wastewater injection, CO2 sequestration has a potential to induce earthquake along weak faults, which can be considered a negative impact on safety and public opinion. This study shows the physical mechanisms of potential seismic hazards along basement faults driven by CO2 sequestration under variation in geological and operational constraints. Specifically we compare the poroelastic behaviors between multiphase flow and single-phase flow cases, highlighting specific needs of evaluating induced seismicity associated with CO2 sequestration. In contrast to single-phase injection scenario, slower migration of the CO2 plume than pressure pulse may delay accumulation of pressure and stress along basement faults that may not be mitigated immediately by shut-in of injection. The impact of multiphase flow system, therefore, needs to be considered for proper monitoring and mitigation strategies.

The mesostructure of porous electrodes used in lithium-ion batteries strongly influences cell performance. Accurate imaging of the distribution of phases in these electrodes would allow this relationship to be better understood through simulation. However, imaging the nanoscale features in these components is challenging. While scanning electron microscopy is able to achieve the required resolution, it has well established difficulties imaging porous media. This is because the flat imaging planes prepared using focused ion beam milling will intersect with the pores, which makes the images hard to interpret as the inside walls of the pores are observed. It is common to infiltrate porous media with resin prior to imaging to help resolve this issue, but both the nanoscale porosity and the chemical similarity of the resins to the battery materials undermine the utility of this approach for most electrodes. In this study, a technique is demonstrated which uses in situ infiltration of platinum to fill the pores and thus enhance their contrast during imaging. Reminiscent of the Japanese art of repairing cracked ceramics with precious metals, this technique is referred to as the kintsugi method. The images resulting from applying this technique to a conventional porous cathode are presented and then segmented using a multi-channel convolutional method. We show that while some cracks in active material particles were empty, others appear to be filled (perhaps with the carbon binder phase), which will have implications for the rate performance of the cell. Energy dispersive X-ray spectroscopy was used to validate the distribution of phases resulting from image analysis, which also suggested a graded distribution of the binder relative to the carbon additive. The equipment required to use the kintsugi method is commonly available in major research facilities and so we hope that this method will be rapidly adopted to improve the imaging of electrode materials and porous media in general.

Erosion and flooding impacts to Arctic coastal environments are intensifying with nearshore oceanographic conditions acting as a key environmental driver. Robust and comprehensive assessment of the nearshore oceanographic conditions require knowledge of the following boundary conditions: incident wave energy, water level, incident wind energy, ocean temperature and salinity, bathymetry, and shoreline orientation. The number of offshore oceanographic boundary conditions can be large, requiring a significant computational investment to reproduce nearshore conditions. This present study develops location-independent typologies to reduce the number of boundary conditions needed to assess nearshore oceanographic environments in both a Historical (2007–2019) and Future (2020–2040) timespan along the Alaskan North Slope. We used WAVEWATCH III® and Delft3D Flexible Mesh model output from six oceanographic sites located along a constant ∼50 m bathymetric line spanning the Chukchi to Beaufort Seas. K-means clustering was applied to the energy-weighted joint-probability distribution of significant wave height (Hs) and peak period (Tp). Distributions of wave and wind direction, wind speed, and water level associated with location-independent centroids were assigned single values to describe a reduced order, typological rendition of offshore oceanographic conditions. Reanalysis data (e.g., ASRv2, ERA5, and GOFS) grounded the historical simulations while projected conditions were obtained from downscaled GFDL-CM3 forced under RCP8.5 conditions. Location-dependence for each site is established through the occurrence joint-probability distribution in the form of unique scaling factors representing the fraction of time that the typology would occupy over a representative year. As anticipated, these typologies show increasingly energetic ocean conditions in the future. They also enable computationally efficient simulation of the nearshore oceanographic environment along the North Slope of Alaska for better characterization of coastal processes (e.g., erosion, flooding, or sediment transport).

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

Multi-phase problems have so many more unknowns, we’d like to have a tool to constrain some open questions related to microstructure and twin & dislocation behavior. We want an atomistic scale perspective on aspects of strength. Some multi-scale questions accessible to atomistic study: What lattice-specific behavior influences dislocation production/mobility and/or twinning? Do the phase transformations wipe-out, modify or preserve grain size and orientation? Does plastic strain reset at phase transition? If so under what conditions? Tin is the material chosen for the effort because it is non-hazardous and has multiple accessible solid phases at relatively low pressures.

Accurate diagnosis of areal density ($\textit{ρR}$) is critical for the inference of performance metrics in inertial confinement fusion implosions. One potential source of error in this diagnosis is the existence of low mode perturbations in the imploding target, which lead to asymmetries in the inference of the $\textit{ρR}$ from different lines of sight. Here, the error accrued as a result of limited coverage of the sphere due to a finite number of detectors is quantified, and the development of a forward scatter measurement from the OMEGA neutron time-of-flight detectors is motivated. In this work, a method by which the 1D-equivalent 4 π-averaged $\langleρR\rangle$ can be reconstructed, if accurate mode information can be diagnosed by other means, is validated.

Infrasound is generated by a variety of natural and anthropogenic sources. Infrasonic waves travel through the dynamic atmosphere, which can change on the order of minutes to hours. Infrasound propagation largely depends on the wind and temperature structure of the atmosphere. Numerical weather prediction models are available to provide atmospheric specifications, but uncertainties in these models exist and they are computationally expensive to run. Machine learning has proven useful in predicting tropospheric weather using Long Short-Term Memory (LSTM) networks. An LSTM network is utilized to make atmospheric specification predictions up to ~30 km for three different training and testing scenarios: (a) the model is trained and tested using only radiosonde data from the Albuquerque, NM, USA station, (b) the model is trained on radiosonde stations across the contiguous US, excluding the Albuquerque, NM, USA station, which was reserved for testing, and (c) the model is trained and tested on radiosonde stations across the contiguous US. Long Short-Term Memory predictions are compared to a state-of-the-art reanalysis model and show cases where the LSTM outperforms, performs equally as well, or underperforms in comparison to the state-of-the-art. Regional and temporal trends in model performance across the US are also discussed. Results suggest that the LSTM model is a viable tool for predicting atmospheric specifications for infrasound propagation modeling.

We report peridynamics is a relatively new field in continuum mechanics that has developed over the past 20 years. This paper studies the evolution of collaborations in the field of peridynamics since its inception using social network analysis. For this purpose, we construct a network for each year from 2000 to 2019 describing co-authorship between scientists in peridynamics. In these networks, each node represents a scientist and each link connects two co-authoring scientists with a link weight representing the frequency and strength of their collaboration; each network as a whole can be thought of as a graph representation of the peridynamics community for the given year. By constructing a network for each year, we are able to analyze the evolution of the network in time and discuss the implications of this evolution for the peridynamics community. Our study demonstrates that the peridynamics community has been growing exponentially in size in recent years. Centrality metrics are also used to identify the most collaborative scientists in the community. Moreover, we compute link recommendations based on both elevating a scientist’s position in the network with respect to certain centrality metrics or closing structural holes in the network identified with persistent homology. We further extend the analysis to higher-order networks whose nodes represent groups of scientists in the community and whose links connect collaborating groups. In some sense, our work studies the past, present, and future of the peridynamics community.

Some past and future supercomputer nodes incorporate High- Bandwidth Memory (HBM). Compared to standard DRAM, HBM has similar latency, higher bandwidth and lower capacity. In this paper, we evaluate algorithms for managing High- Bandwidth Memory automatically. Previous work suggests that, in the worst case, performance is extremely sensitive to the policy for managing the channel to DRAM. Prior theory shows that a priority-based scheme (where there is a static strict priority-order among p threads for channel access) is O(1)-competitive, but FIFO is not, and in the worst case is ?(p) competitive. Following this theoretical guidance would be a disruptive change for vendors, who currently use FIFO variants in their DRAMcontroller hardware. Our goal is to determine theoretically and empirically whether we can justify recommending investment in priority-based DRAM controller hardware. In order to experiment with DRAM channel protocols, we chose a theoretical model, validated it against real hardware, and implemented a basic simulator. We corroborated the previous theoretical results for the model, conducted a parameter sweep while running our simulator on address traces from memory bandwidth-bound codes (GNU sort and TACO sparse matrix-vector product), and designed better channel-access algorithms.

Symbolic regression (SR) with a multi-gene genetic program has been used to elucidate new empirical equations describing diffusion in Lennard-Jones (LJ) fluids. Examples include equations to predict self-diffusion in pure LJ fluids and equations describing the finite-size correction for self-diffusion in binary LJ fluids. The performance of the SR-obtained equations was compared to that of both the existing empirical equations in the literature and to the results from artificial neural net (ANN) models recently reported. It is found that the SR equations have improved predictive performance in comparison to the existing empirical equations, even though employing a smaller number of adjustable parameters, but show an overall reduced performance in comparison to more extensive ANNs.