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.
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.
Cooper, Samuel J.; Roberts, Scott A.; Liu, Zhao; Winiarski, Bartlomiej
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.
Eymold, William K.; Flanary, Christopher; Erikson, Li; Nederhoff, Kees; Chartrand, Chris C.; Jones, Craig; Kasper, Jeremy; Bull, Diana L.
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.
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.
A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.
Security engineering approaches can often focus on a particular domain—physical security, cyber security, or personnel security, for example. Yet, security systems engineering consistently faces challenges requiring socio-technical solutions to address evolving and dynamic complexity. While some drivers of this complexity stem from complex risk environments, innovative adversaries, and disruptive technologies, other drivers are endogenous and emerge from the interactions across security engineering approaches. In response, INCOSE's Systems Security Working Group identified the need to better coordinate “disparate security solutions [that] operate independently” as one of eleven key concepts in their IS21 FuSE Security Roadmap. From this perspective, this need for “security orchestration” aligns with the perspective that security is a property that emerges from interactions within complex systems. Current efforts at Sandia National Laboratories are developing a systems security engineering approach that describes high consequence facility (HCF) security as a multidomain set of interacting layers. The result is a multilayered network (MLN)-based approach that captures the interactions between infrastructure, physical components, digital components, and humans in nuclear security systems. This article will summarize the MLN-based approach to HCF security and describe two preliminary results demonstrating potential benefits from incorporating interactions across disparate security solutions. Here, leveraging the logical structure of networks, this MLN model-based approach provides an example of how security orchestration provides enhanced systems security engineering solutions.
Spark-ignition (SI) engine efficiency can be increased by operating lean and with increased compression ratio (CR), but both of these measures tend to increase the propensity for undesirable acoustic knock generation. It is well known that increased CR makes the engine more prone to knock due to increased combustion pressures and temperatures, but it may be less well understood why lean operation would exacerbate knock generation. For typical gasoline-range fuels, the laminar flame speed becomes very low (roughly only 20% compared to stoichiometric conditions) for an air-excess ratio (λ) of 2. Indirectly, this exacerbates the knock challenge in two ways; a) it may necessitate operation with a combustion phasing near Top Dead Center (TDC) to complete the combustion before expansion cooling occurs, b) it increases cycle-to-cycle variations, making it more challenging to operate near the knock limits. In addition, the high intake pressure required for lean operation (nearly a factor of two higher for λ = 2 compared to λ = 1) increases the oxygen concentration which promotes end-gas autoignition and knock generation. Towards overcoming these challenges of lean combustion, this study aims to provide a better understanding of fuel autoignition under various conditions. First, to reveal the octane appetite under lean conditions, this experimental work utilized fuels of varying Research Octane Number (RON) and octane sensitivity (S). It was found that lean operation favored fuels that have high RON and high S since those were less knock limited. However, two compositionally different fuels with similarly high RON and S exhibited notable difference in knock limits under lean operation, indicating that RON and S may fail to accurately rank order fuels' knock propensity. Second, the experiments show that under boosted conditions end-gas autoignition becomes sensitive to the level of trapped residual nitric oxide (NO), which in turn is very sensitive to variations of both actual λ and combustion phasing, among other factors. The results suggest that strong knock-suppression benefits could be realized if single-ppm NO mole fraction can be consistently maintained in the reactants. Finally, it is noted that maintaining knock-free operation is particularly important for lean operation because the lower peak combustion temperatures lower the speed of sound, which in turn shifts the frequency content of the in-cylinder knock to a lower frequency range. Lower knock frequencies can increase the transmission efficiency from the combustion chamber to the outer surfaces of the engine, potentially increasing engine noise levels if knock occurs.
Wall impingement and fuel film deposition in gasoline direct injection engines under cold start conditions are major concerns for emissions reduction. However, it is challenging to study the dynamics of film deposition under realistic conditions because of the difficulty of measuring the thicknesses of these microscale films. Low-coherence interferometry provides a quantitative optical film thickness measurement technique that can be applied to study this problem. This work presents the first high-speed spectral low-coherence interferometry measurements of impinging gasoline direct injection sprays. The feasibility and practical concerns associated with high-speed low-coherence interferometry systems are explored. Two approaches to spectral low-coherence interferometry: Michelson interferometry and Fizeau interferometry, were implemented and are compared. The results show that Fizeau interferometry is the better option for measurements of impinging sprays in closed spray vessels. The high-speed low-coherence interferometry system was applied in the Fizeau configuration to measure time-resolved film thickness of impinging sprays under engine-relevant conditions to demonstrate its capabilities.
The purpose of this report is to document the probable origin and cause of the 2022 Lurance Fire. As weather was mostly clear skies and lightning was not a factor, the fire had to be the result of human activity.
Brillouin based distributed fiber sensors present a unique set of characteristics amongst fiber sensing architectures. They are able to measure absolute strain and temperature over long distances, with high spatial resolution, and very large dynamic range in off-the-shelf fiber. However, Brillouin sensors traditionally provide only modest sensitivity due to the weak dependence of the Brillouin frequency on strain and the high signal to noise ratio required to identify the resonance’s peak frequency to within a small fraction of its linewidth. Recently, we introduced a technique which substantially improves the precision of Brillouin fiber sensors by exciting a series of lasing modes in a fiber loop cavity that experience Brillouin amplification at discrete locations in the fiber. The narrow-linewidth and high intensity of the lasing modes enabled ultra-low noise Brillouin sensors with large dynamic range. However, our initial demonstration was only modestly distributed: measuring strain at 40, non-contiguous positions along a 400 m fiber. In this work, we greatly extend this methodology to enable fully distributed sensing at 1000 contiguous locations along 3.5 km of fiber—an order of magnitude increase in sensor count and range. This highly-multiplexed Brillouin fiber laser sensor provides a strain noise as low as 34 nε/√Hz and we analyze the limiting factors in this approach.