Jang, Taejin; Mishra, Lubhani; Roberts, Scott A.; Planden, Brady; Subramaniam, Akshay; Uppaluri, Maitri; Linder, David; Gururajan, Mogadalai P.; Zhang, Ji G.; Subramanian, Venkat R.
Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. The article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.
This report is a companion document to a series of six white papers, prepared jointly by the Proliferation Resistance and Physical Protection Working Group (PRPPWG) and the six System Steering Committees (SSCs) and provisional System Steering Committees (pSSCs). This publication is an update to a similar series published in 2011 presenting crosscutting Proliferation Resistance & Physical Protection (PR&PP) characteristics for the six systems selected by the Generation IV International Forum (GIF) for further research and development, namely: the Lead-cooled Fast Reactor (LFR), the Sodium-cooled fast Reactor (SFR), the Very high temperature reactor (VHTR), the gas-cooled fast reactor (GFR), the Molten salt reactor (MSR) and the Supercritical water–cooled reactor (SCWR).
Thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio-posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens-is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300-430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.
An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. We show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child-Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.
Many engineering design problems can be formulated as decisions between two possible options. This is the case, for example, when a quantity of interest must be maintained below or above some threshold. The threshold thereby determines which input parameters lead to which option, and creates a boundary between the two options known as the decision boundary. This report details a machine learning approach for estimating decision boundaries, based on support vector machines (SVMs), that is amenable to large scale computational simulations. Because it is computationally expensive to evaluate each training sample, the approach iteratively estimates the decision boundary in a manner that requires relatively few training samples to glean useful estimates. The approach is then demonstrated on three example problems from structural mechanics and heat transport.
Villa, Daniel V.; Schostek, Tyler; Bianchi, Carlo; Macmillan, Madeline; Carvallo, Juan P.
The Multi-scenario extreme weather simulator (MEWS) is a stochastic weather generation tool. The MEWS algorithm uses 50 or more years of National Oceanic and Atmospheric Association (NOAA) daily summaries [1] for maximum and minimum temperature and NOAA climate norms [2] to calculate historical heat wave and cold snap statistics. The algorithm takes these statistics and shifts them according to multiplication factors provided in the Intergovernmental Panel on Climate Change (IPCC) physical basis technical summary [3] for heat waves.
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 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.
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 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.
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.
We are currently witnessing the dawn of hydrogen (H2) economy, where H2 will soon become a primary fuel for heating, transportation, and longdistance and long-term energy storage. Among diverse possibilities, H2 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 H2 storage densities on both a volumetric and gravimetric basis. However, a critical bottleneck for the use of H2 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 H2 storage. After densification, this monoMOF stores 46 g L–1 H2 at 50 bar and 77 K and delivers 41 and 42 g L–1 H2 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 H2 gas when compared with benchmark materials and an 83% reduction compared to compressed H2 gas. Our findings represent a substantial step forward in the application of high-density materials for volumetric H2 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 Cd3As2 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 Cd3As2.
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.
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.
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.
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 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.
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.
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.
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).
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.
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.
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.
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.
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.
As the U.S. electrifies the transportation sector, cyberattacks targeting vehicle charging could impact several critical infrastructure sectors including power systems, manufacturing, medical services, and agriculture. This is a growing area of concern as charging stations increase power delivery capabilities and must communicate to authorize charging, sequence the charging process, and manage load (grid operators, vehicles, OEM vendors, charging network operators, etc.). The research challenges are numerous and complicated because there are many end users, stakeholders, and software and equipment vendors interests involved. Poorly implemented electric vehicle supply equipment (EVSE), electric vehicle (EV), or grid operator communication systems could be a significant risk to EV adoption because the political, social, and financial impact of cyberattacks — or public perception of such — would ripple across the industry and produce lasting effects. Unfortunately, there is currently no comprehensive EVSE cybersecurity approach and limited best practices have been adopted by the EV/EVSE industry. There is an incomplete industry understanding of the attack surface, interconnected assets, and unsecured inter faces. Comprehensive cybersecurity recommendations founded on sound research are necessary to secure EV charging infrastructure. This project provided the power, security, and automotive industry with a strong technical basis for securing this infrastructure by developing threat models, determining technology gaps, and identifying or developing effective countermeasures. Specifically, the team created a cybersecurity threat model and performed a technical risk assessment of EVSE assets across multiple manufacturers and vendors, so that automotive, charging, and utility stakeholders could better protect customers, vehicles, and power systems in the face of new cyber threats.
Efforts to streamline and codify wave energy resource characterization and assessment for regional energy planning and wave energy converter (WEC) project development have motivated the recent development of resource classification systems. Given the unique interplay between WEC absorption and resource attributes, viz, available wave power frequency, directionality, and seasonality, various consensus resource classification metrics have been introduced. However, the main international standards body for the wave energy industry has not reached consensus on a wave energy resource classification system designed with clear goals to facilitate resource assessment, regional energy planning, project site selection, project feasibility studies, and selection of WEC concepts or archetypes that are most suitable for a given wave energy climate. A primary consideration of wave energy generation is the available energy that can be captured by WECs with different resonant frequency and directional bandwidths. Therefore, the proposed classification system considers combinations of three different wave power classifications: the total wave power, the frequency-constrained wave power, and the frequency-directionally constrained wave power. The dominant wave period bands containing the most wave power are sub-classification parameters that provide useful information for designing frequency and directionally constrained WECs. The bulk of the global wave energy resource is divided into just 22 resource classes representing distinct wave energy climates that could serve as a common language and reference framework for wave energy resource assessment if codified within international standards.
Incipient melting is a phenomenon that can occur in aluminum alloys where solute rich areas, such as grain boundaries, can melt before the rest of the material; incipient melting can degrade mechanical and corrosion properties and is irreversible, resulting in material scrapping. After detecting indications of incipient melting as the cause of failure in 7075 aluminum alloy parts (AA7075), a study was launched to determine threshold temperature for incipient melting. Samples of AA7075 were solution annealed using temperatures ranging from 870-1090F. A hardness profile was developed to demonstrate the loss of mechanical properties through the progression of incipient melting. Additionally, Zeiss software Zen Core Intellesis was utilized to more accurately quantify the changes in microstructural properties as AA7075 surpassed the onset of incipient melting. The results from this study were compared with previous AA7075 material that demonstrated incipient melting.
Khan, Rabia; Miller, Michelle; Valerio, Analise; Willhite, Isabella; Hernandez, Zachary
A need for a standardized Division 10000 onboarding program for virtual hires was identified by management to formalize the way employees and interns are onboarded and trained into Division 10000. This white paper provides effective short and long-term suggestions in the efforts of improving virtual onboarding. Data suggests that remote work is going to become the forefront of many industry practices, which indicates the need of a standardized virtual onboarding practices. With our research, gap assessments, benchmarking, and conducting interviews both internally and externally, we found that clarity, culture, and connection proved to be the strongest solutions in order to maintain Sandia’s competitive edge and sustain workers both remote and in-person.
Sandia National Laboratories has tested and evaluated an updated version of the MB3a infrasound sensor, designed by CEA and manufactured by SeismoWave. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, full scale, self-noise, dynamic range, response, passband, sensitivity variation due to changes in barometric pressure and temperature, and sensitivity to acceleration. The MB3a infrasound sensors are being evaluated for use in the International Monitoring System (IMS) of the Preparatory Commission to the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO).
Sandia National Laboratories has tested and evaluated a new version of the Chaparral 64S infrasound sensor, designed and manufactured by Chaparral Physics. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, full scale, self-noise, dynamic range, response, passband, sensitivity variation due to changes in barometric pressure and temperature, and sensitivity to acceleration. The Chaparral 64S infrasound sensors are being evaluated for use in the International Monitoring System (IMS) of the Preparatory Commission to the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO).
Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5 × 0.5 for 1991-2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.
Accurate diagnosis of areal density (ρ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 ρ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. A method by which the 1D-equivalent 4π-averaged ⟨ ρ R ⟩ can be reconstructed, if accurate mode information can be diagnosed by other means, is validated.
Sandia National Laboratories performed tests to address the potential vulnerability concerns of a coupled High-Altitude Electromagnetic Pulse (HEMP) inducing secondary coupling onto critical instrumentation and control cables in a nuclear power plant, with specific focus on early-time HEMP. Three types of receiving cables in nine configurations were tested to determine transfer functions between two electrically separated cables referenced to the common mode input current on the transmitting cable. One type of transfer function related the input short circuit current and resulting open circuit voltage on the receiving cable. The other transfer function related the input short circuit current and the resulting short circuit current on the receiving cable. A 500 A standard HEMP waveform was input into the transfer functions to calculate peak coupling values on the receiving cables. The highest level of coupling using the standard waveform occurred when cables were in direct contact, with a peak short circuit current of 85 A and open circuit voltage of 9.8 kV, while configurations with separated cables predicted coupling levels of less than 5 A or 500 V.
The objective of this milestone was to finish integrating GenTen tensor software with combustion application Pele using the Ascent in situ analysis software, partnering with the ALPINE and Pele teams. Also, to demonstrate the usage of the tensor analysis as part of a combustion simulation.
Although there are increasing numbers of distributed energy resources (DERs) and microgrids being deployed, current IEEE and utility standards generally strictly limit their interconnection inside secondary networks. Secondary networks are low-voltage meshed (non-radial) distribution systems that create redundancy in the path from the main grid source to each load. This redundancy provides a high level of immunity to disruptions in the distribution system, and thus extremely high reliability of electric power service. There are two main types of secondary networks, called grid and spot secondary networks, both of which are used worldwide. In the future, primary networks in distribution systems that might include looped or meshed distribution systems at the primary-voltage (mediumvoltage) level may also become common as a means for improving distribution reliability and resilience. The objective of this multiyear project is to increase the adoption of microgrids in secondary networks and meshed distribution systems by developing novel protection schemes that allow for safe reliable operation of DERs in secondary networks. We will address these challenges by working with the appropriate stakeholders of secondary network operators, protection vendors, and standards committee. The outcomes of this project include: a) development and/or demonstration of candidate methods for enabling protection of secondary networks containing high levels of DER; b) development of modeling and testing tools for protection systems designed for use with secondary networks including DERs; and c) development of new industrial partnerships to facilitate widespread results dissemination and eventual commercialization of results as appropriate.
In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous operating modes which induce new stresses on the power electronics components. In this work, an experimental and theoretical framework is introduced which couples laboratory-measured component stress with advanced inverter functionality and derives a reduction in useful lifetime based on an applicable reliability model. Multiple DER devices were instrumented to calculate the additional component stress under multiple reactive power setpoints to estimate associated DER lifetime reductions. A clear increase in switch loss was demonstrated as a function of irradiance level and power factor. This is replicated in the system-level efficiency measurements, although magnitudes were different—suggesting other loss mechanisms exist. Using an approximate Arrhenius thermal model for the switches, the experimental data indicate a lifetime reduction of 1.5% when operating the inverter at 0.85 PF—compared to unity PF—assuming the DER failure mechanism thermally driven within the H-bridge. If other failure mechanisms are discovered for a set of power electronics devices, this testing and calculation framework can easily be tailored to those failure mechanisms.
The information impulse function (IIF), running Variance, and local Hölder Exponent are three conceptually different time-series evaluation techniques. These techniques examine time-series for local changes in information content, statistical variation, and point-wise smoothness, respectively. Using simulated data emulating a randomly excited nonlinear dynamical system, this study interrogates the utility of each method to correctly differentiate a transient event from the background while simultaneously locating it in time. Computational experiments are designed and conducted to evaluate the efficacy of each technique by varying pulse size, time location, and noise level in time-series. Our findings reveal that, in most cases, the first instance of a transient event is more easily observed with the information-based approach of IIF than with the Variance and local Hölder Exponent methods. While our study highlights the unique strengths of each technique, the results suggest that very robust and reliable event detection for nonlinear systems producing noisy time-series data can be obtained by incorporating the IIF into the analysis.
Mcadams, Daniel R.; Johnson, Erik B.; Christian, James F.; Blakeley, Richard; Weststrate, Evan; He, Jane; Bader, Sven; Crotts, Bradley; Reinke, Charles M.
Elmslie, Timothy A.; Startt, Jacob K.; Soto-Medina, Sujeily; Feng, Keke; Zappala, Emma; Frandsen, Benjamin A.; Meisel, Mark W.; Dingreville, Remi P.; Hamlin, James J.
Magnetic, specific heat, and structural properties of the equiatomic Cantor alloy system are reported for temperatures between 5 and 300 K, and up to fields of 70 kOe. Magnetization measurements performed on as-cast, annealed, and cold-worked samples reveal a strong processing history dependence and that high-temperature annealing after cold working does not restore the alloy to a "pristine"state. Measurements on known precipitates show that the two transitions, detected at 43 and 85 K, are intrinsic to the Cantor alloy and not the result of an impurity phase. Experimental and ab initio density functional theory computational results suggest that these transitions are a weak ferrimagnetic transition and a spin-glass-like transition, respectively, and magnetic and specific heat measurements provide evidence of significant Stoner enhancement and electron-electron interactions within the material.
We consider a bilevel attacker–defender problem to find the worst-case attack on the relays that control transmission grid components. The attacker infiltrates some number of relays and renders all of the components connected to them inoperable with the goal of maximizing load shed. The defender responds by minimizing the resulting load shed, redispatching using a DC optimal power flow (DCOPF) problem on the remaining network. Though worst-case interdiction problems on the transmission grid have been studied for years, there remains a need for exact and scalable methods. Methods based on using duality on the inner problem rely on the bounds of the dual variables of the defender problem in order to reformulate the bilevel problem as a mixed integer linear problem. Valid dual bounds tend to be large, resulting in weak linear programming relaxations and, hence, making the problem more difficult to solve at scale. Often smaller heuristic bounds are used, resulting in a lower bound. In this work, we also consider a lower bound, but instead of bounding the dual variables, we drop the constraints corresponding to Ohm’s law, relaxing DCOPF to capacitated network flow. We present theoretical results showing that, for uncongested networks, approximating DCOPF with network flow yields the same set of injections and, thus, the same load shed, which suggests that this restriction likely gives a high-quality lower bound in the uncongested case. Furthermore, we show that, in the network flow relaxation of the defender problem, the duals are bounded by one, so we can solve our restriction exactly. Finally, because the big-M values in the linearization are equal to one and network flow has a well-known structure, we see empirically that this formulation scales well computationally with increased network size. Through empirical experiments on 16 networks with up to 6,468 buses, we find that this bound is almost always as tight as we can get from guessing the dual bounds even for congested networks in which the theoretical results do not hold. In addition, calculating the bound is approximately 150 times faster than achieving the same bound with the reformulation guessing the dual bounds.
Light-matter interaction optimization in complex nanophotonic structures is a critical step towards the tailored performance of photonic devices. The increasing complexity of such systems requires new optimization strategies beyond intuitive methods. For example, in disordered photonic structures, the spatial distribution of energy densities has large random fluctuations due to the interference of multiply scattered electromagnetic waves, even though the statistically averaged spatial profiles of the transmission eigenchannels are universal. Classification of these eigenchannels for a single configuration based on visualization of intensity distributions is difficult. However, successful classification could provide vital information about disordered nanophotonic structures. Emerging methods in machine learning have enabled new investigations into optimized photonic structures. In this work, we combine intensity distributions of the transmission eigenchannels and the transmitted speckle-like intensity patterns to classify the eigenchannels of a single configuration of disordered photonic structures using machine learning techniques. Specifically, we leverage supervised learning methods, such as decision trees and fully connected neural networks, to achieve classification of these transmission eigenchannels based on their intensity distributions with an accuracy greater than 99%, even with a dataset including photonic devices of various disorder strengths. Simultaneous classification of the transmission eigenchannels and the relative disorder strength of the nanophotonic structure is also possible. Our results open new directions for machine learning assisted speckle-based metrology and demonstrate a novel approach to classifying nanophotonic structures based on their electromagnetic field distributions. These insights can be of paramount importance for optimizing light-matter interactions at the nanoscale.
Finfrock, Christopher B.; Ellyson, Benjamin; Likith, Sri R.J.; Smith, Douglas R.; Rietema, Connor; Saville, Alec I.; Thrun, Melissa M.; Becker, C.G.; Araujo, Ana L.; Pavlina, Erik J.; Hu, Jun; Park, Jun-Sang; Clarke, Amy J.; Clarke, Kester D.
Understanding the deformation-induced martensitic transformation (DIMT) is critical for interpreting the structure-property relationships that govern the performance of transformation-induced plasticity (TRIP) assisted steels. However, modern TRIP-assisted steels often exhibit DIMT kinetics that are not easily captured by existing empirical models based on bulk tensile strain. We address this challenge by combined bulk uniaxial tensile tests and in-situ high energy synchrotron X-ray diffraction, which resolved the phase volume fractions, stress-strain response, and microstructure evolution of each constituent phase. A modification of the Olson-Cohen model is implemented, which describes the martensitic transformation kinetics as a function of the estimated partitioned strain in austenite, rather than the bulk tensile strain. This DIMT kinetic model is used as a framework to clarify the root cause of an insufficiently understood toughness trough reported for TRIP-assisted steels during deformation at elevated temperatures. Here, the importance of the temperature-dependent toughness is discussed, based on the opportunity to modify deformation processes to tailor the DIMT kinetics and mechanical properties during forming and in service.
Failure detection methods are of significant interest for photovoltaic (PV) site operators to help reduce gaps between expected and observed energy generation. Current approaches for field-based fault detection, however, rely on multiple data inputs and can suffer from interpretability issues. In contrast, this work offers an unsupervised statistical approach that leverages hidden Markov models (HMM) to identify failures occurring at PV sites. Using performance index data from 104 sites across the United States, individual PV-HMM models are trained and evaluated for failure detection and transition probabilities. This analysis indicates that the trained PV-HMM models have the highest probability of remaining in their current state (87.1% to 93.5%), whereas the transition probability from normal to failure (6.5%) is lower than the transition from failure to normal (12.9%) states. A comparison of these patterns using both threshold levels and operations and maintenance (O&M) tickets indicate high precision rates of PV-HMMs (median = 82.4%) across all of the sites. Although additional work is needed to assess sensitivities, the PV-HMM methodology demonstrates significant potential for real-time failure detection as well as extensions into predictive maintenance capabilities for PV.
Ducted fuel injection (DFI) is a novel combustion strategy that has been shown to significantly attenuate soot formation in diesel engines. While previous studies have used optical diagnostics and optical filter smoke number methods to show that DFI reduces in-cylinder soot formation and engine-out soot emissions, respectively, this is the first study to measure solid particle number (PN) emissions in addition to particle mass (PM). Furthermore, this study quantitatively evaluates the use of transient particle instruments for measuring particles from skip-fired operation in an optical single cylinder research engine (SCRE). Engine-out PN was measured using an engine exhaust particle sizer following a catalytic stripper, and PM was measured using a photoacoustic analyzer. The study improves on earlier preliminary emissions studies by clearly showing that DFI reduces overall PM by 76%–79% and PN for particles larger than 23 nm by 77% relative to conventional diesel combustion at a 1200-rpm, 13.3-bar gross indicated mean effective pressure operating condition. The degree of engine-out PM reduction with DFI was similar across both particulate measurement instruments used in the work. Through the use of bimodal distribution fitting, DFI was also shown to reduce the geometric mean diameter of accumulation mode particles by 26%, similar to the effects of increased injection pressure in conventional diesel combustion systems. This work clearly shows the significant solid particulate matter reductions enabled by DFI while also demonstrating that engine-out PN can be accurately measured from an optical SCRE operating in a skip-fired mode. Based on these results, it is believed that DFI has the potential to enable fuel savings when implemented in multi-cylinder engines, both by lowering the required frequency of active diesel particulate filter regeneration, and by reducing the backpressure imposed by exhaust filtration systems.
Cunningham, W.S.; Riano, J.S.; Wang, Wenbo; Hwang, Sooyeon; Hattar, Khalid M.; Hodge, Andrea M.; Trelewicz, Jason R.
Targeted doping of grain boundaries is widely pursued as a pathway for combating thermal instabilities in nanocrystalline metals. However, certain dopants predicted to produce grain-boundary-segregated nanocrystalline configurations instead form small nanoprecipitates at elevated temperatures that act to kinetically inhibit grain growth. Here, thermodynamic modeling is implemented to select the Mo–Au system for exploring the interplay between thermodynamic and kinetic contributions to nanostructure stability. Using nanoscale multilayers and in situ transmission electron microscopy thermal aging, evolving segregation states and the corresponding phase transitions are mapped with temperature. The microstructure is shown to evolve through a transformation at lower homologous temperatures (<600 °C) where solute atoms cluster and segregate to the grain boundaries, consistent with predictions from thermodynamic models. An increase in temperature to 800 °C is accompanied by coarsening of the grain structure via grain boundary migration but with multiple pinning events uncovered between migrating segments of the grain boundary and local solute clustering. Direct comparison between the thermodynamic predictions and experimental observations of microstructure evolution thus demonstrates a transition from thermodynamically preferred to kinetically inhibited nanocrystalline stability and provides a general framework for decoupling contributions to complex stability transitions while simultaneously targeting a dominant thermal stability regime.