Pulsed power and plasma physics are topics of great study at both Sandia National Laboratories (SNL or Sandia) and the University of New Mexico (UNM). The goal of this research is to further knowledge and understanding of these fields using the resources of both SNL and UNM in three ways. The first way is through the comprehension, application, and testing of theory. Reading and analytically deriving theoretical solutions of problems both real-world and simplified will allow for a fresh perspective and the furthering of the theory. One such theory is Ottinger's generalized theory for voltage measurement in magnetically insulated transmission lines (MITLs). By working through the math, a deeper understanding of the theory is gained from which one may add more physically accurate and/or more detailed physics into the theory. Additionally, understanding the theory lays a good foundation from which one can analyze, test, and compare results to the theory in the following two ways that will advance the fields of pulsed power and plasma physics. The second way is through the modeling and simulation of real-world and simplified problems that utilize and test the afore mentioned theories. Theory can be applied to a simulation domain by using the unstructured time-domain electromagnetic (UTDEM) codes EMPHASIS and EMPIRE as well as the physical modeling software CUBIT, all of which were developed at SNL. Problems such as the modeling and design of the extended MITL on HERMES III, the understanding of space-charge-limited emission from vacuum cathodes, and the interaction between a relativistic electron beam and an ideal gas can all be modeled, simulated, and analyzed with this set of codes. Here the advantage is three-fold. Firstly, theory that describes our understanding of these problems can be put to the test and advanced through iterative simulation and analysis. Secondly, the understanding of these problems will have a positive impact on national security through the advancement of the technological capability of the United States of America. Thirdly, and not unrelated to the prior advantage, is the validation and verification of EMPIRE and EMPHASIS. This segues into the third way, which is through experiment and the comparison of experiment to simulated and theoretical results. Performing experimental comparisons completes the scientific method and grounds all of the work in reality. Being able to physically test theory and simulation is necessary for any real conclusions to be drawn. Another advantage for carrying out experimental work is to advance the physical testing capabilities of SNL. Several systems will be developed and tested through the course of this work that positively impact technological advancement of Sandia National Labs. Lastly, all of the above work will converge to yield a well-rounded perspective that ties the three categories of research together.
The azimuth of an incoming acoustic wave cannot be determined using microbarometers on a free floating balloon. A single observation of infrasound-induced acceleration on a large zero pressure balloon suggested that a motion sensing "aeroseismometer" could fill this gap. Here, a flight test of prototype balloon-borne aeroseismometers is presented. Two balloons, each carrying accelerometers and IMUs, recorded three sets of chemical explosions. The resulting balloon motion time series allows the explosive source to be geolocated. The future of this technology is discussed, along with a planned publication. Finally, recommendations and lessons learned from the campaign are discussed.
A phenomenological model of cavitation is presented, based on observations that both large relative negative pressures and large negative time derivatives of pressure are required for cavitation onset. We simulated two cavitation experiments to generate cavitation scaling parameters for relative pressure drop and rate of pressure drop. Our results show the model, while simple, is effective at reproducing results from laboratory experiments of cavitation. The parameters were then used in conjunction with a human surrogate computational model to predict, at any position within the head, the probability of intracranial cavitation caused by exposure to a blast event. The results suggest that the magnitude of blast overpressure observed in field data is sufficient to cause intracranial cavitation. Our analysis indicates that the helmeted head, when compared to the unhelmeted head configuration, results in a decrease but not elimination of cavitation exposure. When density functions of cavitation probability versus cumulative brain volume are combined with an injury severity model, the results show helmet efficacy at low and moderate risk levels. However, the convergence of unhelmeted and helmeted probability density functions at high-to-excessive risk thresholds indicates the helmet offers diminishing protection at elevated exposure levels, relative to the unhelmeted baseline. Future investigation and collaboration with neuroscience subject matter experts are needed to contextualize the current computational results. While the present work contributes specific and quantified predictions of intracranial cavitation location and severity, more research is required to apply our results to clinical settings with population-based brain injury subjects and controls. The relationship between our intracranial cavitation predictions with their anticipated clinical sequelae remains a topic in need of exploration.
In October of 2020, dust samples were collected from the surface of spent nuclear fuel (SNF) dry storage canisters during an inspection at an inland Independent Spent Fuel Storage Installation, the second inland site at which surface deposits have been sampled. The purpose of the sampling was to assess the composition and abundance of the soluble salts present on the canister surface, information which provides a metric for potential corrosion risks. The samples were delivered to Sandia National laboratories for analysis. At Sandia, the soluble salts were leached from the dust and quantified by ion chromatography. In addition, subsamples of the dust were taken for scanning electron microscopy to determine the texture and mineralogy of the dust and salts. The results of those analyses are presented in this report.
Pyrazole, a member of the structural class of azoles, exhibits molecular properties of interest in pharmaceuticals and materials chemistry, owing to the two adjacent nitrogen atoms in the five-membered ring system. The weakly basic nitrogen atoms of deprotonated pyrazoles have been applied in coordination chemistry, particularly to access coordination polymers and metal-organic frameworks, and homocoupling reactions can in principle provide facile access to bipyrazole ligands. In this context, we summarize recent advances in homocoupling reactions of pyrazoles and other types of azoles (imidazoles, triazoles and tetrazoles) to highlight the utility of homocoupling reactions in synthesizing symmetric bi-heteroaryl systems compared with traditional synthesis. Metal-free reactions and transition-metal catalyzed homocoupling reactions are discussed with reaction mechanisms in detail.
This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.
A stripmap Synthetic Aperture Radar (SAR) image is a long SAR image along some centerline, and formed from multiple synthetic apertures. At issue is that the centerline in the image actually corresponds to an arc on a round earth, and multiple strategies exist for fitting the image centerline to the round earth. Some of those strategies involve Rhumb lines, great circle paths, and great ellipse paths. Some are better than others in polar regions. Notions of parallel flight paths for the radar during data collection also require careful consideration of the geometry of a round earth.
This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 14959and 18019. Revision 14959 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 18019. Along with the newly updated MELCOR Users Guide and Reference Manual, users are aware and able to assess the new capabilities for their modeling and analysis applications.
The motivation for this report is to discuss and present some realistic tree models employed in computational electromagnetics (EM) simulations to study foliage penetration (FOPEN) at Ku-band. The detail obtained in these trees is unprecedented in FOPEN modeling since many studies in this area focus on lower frequencies where precise tree parameters are not required due to the associated large wavelengths relative to the tree dimensions. The focus of this study is in the Ku-band range where the wavelength is notably smaller and the details of the trees have more of an influence on EM waves (i.e. scattering, attenuating, reflecting, diffracting etc.). Therefore, explicit tree parameters are modeled. Also, moderate foliage is of most interest because with less dense foliage t here is a higher percentage of Ku-band transmission. The EM wave and foliage interaction s are simulated with the computational electromagnetics (CEM) Altair FEKO software. The realistic tree model s implemented for simulations are created in the computer-aided design (CAD) software Arbaro and the module CADFEKO that is offered in FEKO. Details of these tree models are provided, and EM simulation results will be discussed in a follow-on report
This report describes the design, development, and testing of a prototype 100 kWt particle-to-supercritical CO2 (sCO2) heat exchanger. An analytic hierarchy process was implemented to compare and evaluate alternative heat-exchanger designs (fluidized bed, shell-and-plate moving packed bed, and shell-and-tube moving packed bed) that could meet the high pressure (≥ 20 MPa) and high temperature (≥ 700 °C) operational requirements associated with sCO2 power cycles. Cost, heat-transfer coefficient, structural reliability, manufacturability, parasitics and heat losses, scalability, compatibility, erosion and corrosion, transient operation, and inspection ease were considered in the evaluation. A 100 kWt shell-and-plate design was selected for construction and integration with Sandia’s falling particle receiver system that heats the particles using concentrated sunlight. Sandia worked with industry to design and construct the moving packed-bed shell-and-plate heat exchanger. Tests were performed to evaluate its performance using both electrical heating and concentrated sunlight to heat the particles. Overall heat transfer coefficients at off-design conditions (reduced operating temperatures and only three stainless steel banks in the counter-crossflow heat exchanger) were measured to be approximately ~25 - 70 W/m2-K, significantly lower than simulated values of >100 W/m2-K. Tests using the falling particle receiver to heat the particles with concentrated sunlight yielded overall heat transfer coefficients of ~35 – 80 W/m2-K with four banks (including a nickel-alloy bank above the three stainless steel banks). The overall heat transfer coefficient was observed to decrease with increasing particle inlet temperatures, which contrasted the results of simulations that showed an increase in heat transfer coefficient with temperature due to increased effective particle-bed thermal conductivity from radiation. The likely cause of the discrepancy was particle-flow maldistributions and funnel flow within the heat exchanger caused by internal ledges and cross-bracing, which could have been exacerbated by increased particle-wall friction at higher temperatures. Additional heat loss at higher temperatures may also contribute to a lower overall heat-transfer coefficient. Design challenges including pressure drop, particle and sCO2 flow maldistribution, and reduced heat transfer coefficient are discussed with approaches for mitigation in future designs. Lessons learned regarding instrumentation, performance characterization, and operation of particle components and sCO2 flow loops are also discussed. Finally, a 200 MWt commercial-scale shell-and-plate heat-exchanger design based on the concepts investigated in this report is proposed.
Dynamic grain growth is demonstrated to be much faster than static grain growth in a body-centered-cubic, interstitial-free steel sheet material at 850∘C. Dynamic grain growth occurs during concurrent plastic deformation at elevated temperature, whereas static grain growth occurs during static annealing. Grain growth during steady-state plastic flow in tension at 850∘C to a true strain of 0.2 at a true-strain rate of 10 - 4 s - 1 doubled grain size, while static annealing for the same time produced no increase in grain size. This is described as dynamic normal grain growth (DNGG) because no abnormally large grains were observed. The recrystallized microstructure of the steel demonstrated a log-normal distribution of grain sizes. DNGG produced bimodal grain size distributions that deviate from the theoretical expectation of a simple shift to larger sizes during normal growth. The bimodal distributions contained a remnant of small grains that were not consumed during grain growth. DNGG produced a crystallographic texture that is unique from both the recrystallized material and that produced by lattice rotation alone. DNGG strengthened the { 111 } ⟨ 110 ⟩ and { 111 } ⟨ 112 ⟩ components of the strong γ-fiber component in the original recrystallization texture. Lattice rotation from tensile deformation, by contrast, strengthened the α-fiber components that intersect the original γ-fiber.
The fatigue strength of additively manufactured metallic parts in their as-built surface condition is mainly dominated by the surface roughness. Post-processing is often inevitable to reduce surface roughness effects even though post-processing diminishes the main advantage of additive manufacturing, which is net-shaped direct-to-service production. This study investigates the underlying mechanisms responsible for fatigue failure of additively manufactured 304L stainless steel parts in as-built and machined/polished surface conditions. Both strain- and force-controlled, fully reversed fatigue tests were conducted to gain a comprehensive understanding of surface roughness effects on fatigue behavior. The sensitivity to surface roughness is shown to be dependent on the control mode, with stress-based fatigue tests showing greater sensitivity than strain-based fatigue tests. Moreover, the fatigue life estimation for as-built specimens was performed based on surface roughness parameters as well as the fatigue properties of the specimens in machined/polished surface condition of the material without using any fatigue data of specimens in as-built surface condition. Accordingly, the effect of surface roughness on the fatigue behavior could be estimated reasonably well in both strain-life and stress-life approaches.
We present in this paper the results from a recent study on the stability of the quantum Hall skyrmions state at a Landau level filling factor (ν) close to ν = 1 in a narrow GaAs quantum well. Consistent with previous work, a resonant behavior is observed in the resistively detected NMR measurements. In the subsequent current-voltage (I-V) measurements to examine its breakdown behavior under radio frequency radiations, we observe that the critical current assumes the largest value right at the 75As nuclear resonant frequency. We discuss possible origin for this unexpectedly enhanced stability.
Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. A multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure build-ups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.
We have derived a gas ingress and egress equation from the first principles of ideal gases. This work is intended to benefit the hermetic microelectronics packaging community, but it may be applied to other fields that require a deep understanding of gas ingress and egress dynamics. The equation outlined herein encompasses package material properties, package characteristics, hermetic testing conditions, and service conditions. It serves as a practical utility for calculating package pressure changes due to gas ingress and egress and, therefore, a power tool for component and system service life predictions.
Orthorhombic RMnO3 (R = rare-earth cation) compounds are type-II multiferroics induced by inversion-symmetry-breaking of spin order. They hold promise for magneto-electric devices. However, no spontaneous room-temperature ferroic property has been observed to date in orthorhombic RMnO3. Here, using 3D straining in nanocomposite films of (SmMnO3)0.5((Bi,Sm)2O3)0.5, we demonstrate room temperature ferroelectricity and ferromagnetism with TC,FM ~ 90 K, matching exactly with theoretical predictions for the induced strain levels. Large in-plane compressive and out-of-plane tensile strains (−3.6% and +4.9%, respectively) were induced by the stiff (Bi,Sm)2O3 nanopillars embedded. The room temperature electric polarization is comparable to other spin-driven ferroelectric RMnO3 films. Also, while bulk SmMnO3 is antiferromagnetic, ferromagnetism was induced in the composite films. The Mn-O bond angles and lengths determined from density functional theory explain the origin of the ferroelectricity, i.e. modification of the exchange coupling. Our structural tuning method gives a route to designing multiferroics.
A 250kW hydrogen electrolysis facility was recently installed at the Natural Energy Laboratory of Hawaii Authority's (NELHA's) campus. This facility that will begin operation in 2020 to produce hydrogen for fuel cell buses on the island to demonstrate of the application of hydrogen to decarbonize transportation. Given the size of the electrolysis station, it has the potential to significantly increase electricity costs for the campus, which is subject to energy and peak demand charges from the local utility. In this paper, we analyze the cost of hydrogen production at NELHA given the rate structure options available from the utility. Production costs are estimated using optimal versus constant scheduling of the facility to meet the buses’ demand. A model of the electrolysis station is used to capture changes in production efficiency over the power range in the optimization routine. The effects of combining the station and campus load versus standalone operation and increasing solar generation are also explored. The analyses surrounding this scenario show the importance of multiple factors on the potential profitability of hydrogen production in behind-the-meter applications and show trends that could have implications for other similar installations.
Understanding the mechanistic relationship between the environment, microstructure, and local kinetics of atmospheric corrosion damage remains a central challenge. To address this challenge, this study used laboratory-based X-ray tomography to directly observe attack in-operando over an extended period, enabling insights into the evolving growth kinetics and morphology of individual pits over months of exposure. Damage progression associated with nine pits in a 99.9% pure aluminum wire exposed to chloride salts in humid air was characterized. Most pits grew at a nominally linear rate up until pit death, which occurred within 12–24 h of nucleation. Exceptions to this were observed, with three pits exhibiting bimodal growth kinetics and growing for 40 or more hours. This was explained by secondary droplets that formed near the pits, increasing the cathode area. A corrosion-driven drying mechanism likely contributed to pit death in both cases. Pits first grew into the material followed by lateral expansion.
We present evidence of inverse Hall-Petch behavior for a single-phase high entropy alloy (CoCrFeMnNi) in ultra-high vacuum and show that it is associated with low friction coefficients (~0.3). Grain size measurements by STEM validate a recently proposed dynamic amorphization model that accurately predicts grain size-dependent shear strength in the inverse Hall-Petch regime. Wear rates in the initially soft (coarse grained) material were shown to be remarkably low (~10–6 mm3/N-m), the lowest for any HEA tested in an inert environment where oxidation and the formation of mixed metal-oxide films is mitigated. The combined high wear resistance and low friction are linked to the formation of an ultra-nanocrystalline near-surface layer. The dynamic amorphization model was also used to predict an average high angle grain boundary energy (0.87 J/m2). This value was used to explain cavitation-induced nanoporosity found in the highly deformed surface layer, a phenomenon that has been linked to superplasticity.
This paper presents a multifidelity uncertainty quantification framework called MFNets. We seek to address three existing challenges that arise when experimental and simulation data from different sources are used to enhance statistical estimation and prediction with quantified uncertainty. Specifically, we demonstrate that MFNets can (1) fuse heterogeneous data sources arising from simulations with different parameterizations, e.g simulation models with different uncertain parameters or data sets collected under different environmental conditions; (2) encode known relationships among data sources to reduce data requirements; and (3) improve the robustness of existing multi-fidelity approaches to corrupted data. MFNets construct a network of latent variables (LVs) to facilitate the fusion of data from an ensemble of sources of varying credibility and cost. These LVs are posited as explanatory variables that provide the source of correlation in the observed data. Furthermore, MFNets provide a way to encode prior physical knowledge to enable efficient estimation of statistics and/or construction of surrogates via conditional independence relations on the LVs. We highlight the utility of our framework with a number of theoretical results which assess the quality of the posterior mean as a frequentist estimator and compare it to standard sampling approaches that use single fidelity, multilevel, and control variate Monte Carlo estimators. We also use the proposed framework to derive the Monte Carlo-based control variate estimator entirely from the use of Bayes rule and linear-Gaussian models -- to our knowledge the first such derivation. Finally, we demonstrate the ability to work with different uncertain parameters across different models.
Journal of Physics B: Atomic, Molecular and Optical Physics
Chauhan, Dhaval; Limbachiya, Chetan G.; Tokesi, Karoly; Champion, Christophe; Weck, Philippe F.; Tribedi, Lokesh C.
We report the measurement of the absolute double differential cross sections (DDCS) of secondary electrons emitted due to the ionization of N2 molecule in collisions with fast electrons having energies between 3 and 5 keV. The emitted electrons with energies from 1-500 eV have been measured for different forward and backward emission angles. The measured DDCS have been compared with the state-of-the-art first Born approximation with correct boundary condition (CB1) model calculations as well as with the classical trajectory Monte Carlo (CTMC) method. From the measured DDCS, the single differential cross sections (SDCS) as a function of the emission energies have been computed and eventually the total ionization cross sections (TCS) have been derived. The TCS values are also compared with a semi-empirical calculation, namely, the CSP-ic (complex scattering potential-ionization contribution) model.
Dynamic instability of microtubules is characterized by stochastically alternating phases of growth and shrinkage and is hypothesized to be controlled by the conformation and nucleotide state of tubulin dimers within the microtubule lattice. Specifically, conformation changes (compression) in the tubulin dimer following the hydrolysis of GTP have been suggested to generate stress and drive depolymerization. In the present study, molecular dynamics simulations were used in tandem with in vitro experiments to investigate changes in depolymerization based on the presence of islands of uncompressed (GMPCPP) dimers in the microtubule lattice. Both methods revealed an exponential decay in the kinetic rate of depolymerization corresponding to the relative level of uncompressed (GMPCPP) dimers, beginning at approximately 20% incorporation. This slowdown was accompanied by a distinct morphological change from unpeeling “ram’s horns” to blunt-ended dissociation at the microtubule end. Collectively these data demonstrated that islands of uncompressed dimers can alter the mechanism and kinetics of depolymerization in a manner consistent with promoting rescue events.
Parallel implementations of linear iterative solvers generally alternate between phases of data exchange and phases of local computation. Increasingly large problem sizes and more heterogeneous compute architectures make load balancing and the design of low latency network interconnects that are able to satisfy the communication requirements of linear solvers very challenging tasks. In particular, global communication patterns such as inner products become increasingly limiting at scale. We explore the use of asynchronous communication based on one-sided Message Passing Interface primitives in the context of domain decomposition solvers. In particular, a scalable asynchronous two-level Schwarz method is presented. We discuss practical issues encountered in the development of a scalable solver and show experimental results obtained on a state-of-the-art supercomputer system that illustrate the benefits of asynchronous solvers in load balanced as well as load imbalanced scenarios. Using the novel method, we can observe speedups of up to four times over its classical synchronous equivalent.
The failure forces and fracture strengths of polysilicon microelectromechanical system (MEMS) components in the form of stepped tensile bars with shoulder fillets were measured using a sequential failure chain methodology. Approximately 150 specimens for each of four fillet geometries with different stress concentration factors were tested. The resulting failure force and strength distributions of the four geometries were related by a common sidewall flaw population existing within different effective stressed lengths. The failure forces, strengths, and flaw population were well described by a weakest-link based analytical framework. Finite element analysis was used to verify body-force based expressions for the stress concentration factors and to provide insight into the variation of specimen effective length with fillet geometry. Monte Carlo simulations of flaw size and location, based on the strength measurements, were also used to provide insight into fillet shape and size effects. The successful description of the shoulder fillet specimen strengths provides further empirical support for application of the strength and flaw framework in MEMS fabrication and design optimization.
Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive, and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid freeform fabrication (SFF). Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror, and thereafter, an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is now becoming ubiquitous. This article is intended to highlight the most recent advances in the computation benefit specific to the area of precision SFF as these systems are traversing through the journey of digitalization and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.
One of the most severe obstacles to increasing the longevity of tungsten-based plasma facing components, such as divertor tiles, is the surface deterioration driven by sub-surface helium bubble formation and rupture. Supported by experimental observations at PISCES, this work uses molecular dynamics simulations to identify the microscopic mechanisms underlying suppression of helium bubble formation by the introduction of plasma-borne beryllium. Simulations of the initial surface material (crystalline W), early-time Be exposure (amorphous W-Be) and final WBe2 intermetallic surfaces were used to highlight the effect of Be. Significant differences in He retention, depth distribution and cluster size were observed in the cases with beryllium present. Helium resided much closer to the surface in the Be cases with nearly 80% of the total helium inventory located within the first 2 nm. Moreover, coarsening of the He depth profile due to bubble formation is suppressed due to a one-hundred fold decrease in He mobility in WBe2, relative to crystalline W. This is further evidenced by the drastic reduction in He cluster sizes even when it was observed that both the amorphous W-Be and WBe2 intermetallic phases retain nearly twice as much He during cumulative implantation studies.
Cathodic kinetics in magnesium chloride (MgCl2) solutions were investigated on platinum (Pt) and stainless steel 304 L (SS304 L). Density, viscosity, and dissolved oxygen concentration for MgCl2 solutions were also measured. A 2-electron transfer for oxygen reduction reaction (ORR) on Pt was determined using a rotating disk electrode. SS304 L displayed non-Levich behavior for ORR and, due to ORR suppression and buffering of near surface pH by Mg-species precipitation, the primary cathodic reaction was the hydrogen evolution reaction (HER) in saturated MgCl2. Furthermore, non-carbonate precipitates were found to be kinetically favored. Implications of HER are discussed through atmospheric corrosion and stress corrosion cracking.
Optical polarizers encompass a class of anisotropic materials that pass-through discrete orientations of light and are found in wide-ranging technologies, from windows and glasses to cameras, digital displays and photonic devices. The wire-grids, ordered surfaces, and aligned nanomaterials used to make polarized films cannot be easily reconfigured once aligned, limiting their use to stationary cross-polarizers in, for example, liquid crystal displays. Here we describe a supramolecular material set and patterning approach where the polarization angle in stand-alone films can be precisely defined at the single pixel level and reconfigured following initial alignment. This capability enables new routes for non-binary information storage, retrieval, and intrinsic encryption, and it suggests future technologies such as photonic chips that can be reconfigured using non-contact patterning.
Lozanovski, Bill; Downing, David; Tino, Rance; Du Plessis, Anton; Tran, Phuong; Jakeman, John D.; Shidid, Darpan; Emmelmann, Claus; Qian, Ma; Choong, Peter; Brandt, Milan; Leary, Martin
Additive Manufacturing (AM), commonly referred to as 3D printing, offers the ability to not only fabricate geometrically complex lattice structures but parts in which lattice topologies in-fill volumes bounded by complex surface geometries. However, current AM processes produce defects on the strut and node elements which make up the lattice structure. This creates an inherent difference between the as-designed and as-fabricated geometries, which negatively affects predictions (via numerical simulation) of the lattice's mechanical performance. Although experimental and numerical analysis of an AM lattice's bulk structure, unit cell and struts have been performed, there exists almost no research data on the mechanical response of the individual as-manufactured lattice node elements. This research proposes a methodology that, for the first time, allows non-destructive quantification of the mechanical response of node elements within an as-manufactured lattice structure. A custom-developed tool is used to extract and classify each individual node geometry from micro-computed tomography scans of an AM fabricated lattice. Voxel-based finite element meshes are generated for numerical simulation and the mechanical response distribution is compared to that of the idealised computer-aided design model. The method demonstrates compatibility with Uncertainty Quantification methods that provide opportunities for efficient prediction of a population of nodal responses from sampled data. Overall, the non-destructive and automated nature of the node extraction and response evaluation is promising for its application in qualification and certification of additively manufactured lattice structures.
Driscoll, Darren M.; Dandu, Naveen K.; Hahn, Nathan H.; Seguin, Trevor J.; Persson, Kristin A.; Zavadil, Kevin R.; Curtiss, Larry A.; Balasubramanian, Mahalingam
Ca-ion electrochemical systems have been pushed to the forefront of recent multivalent energy storage advances due to their use of earth-abundant redox materials and their high theoretical specific densities in relation to monovalent or even other more widely explored multivalent-charge carriers. However, significant pitfalls in metal plating and stripping arise from electrolyte decomposition and can be related to the coordination environment around Ca2+ with both the negatively charged anion and the organic-aprotic solvent. In this study, we apply multiple spectroscopic techniques in conjunction with density functional theory to evaluate the coordination environment of Ca2+ across a class of ethereal solvents. Through the combination of X-ray absorption fine structure and time-dependent density functional theory, descriptive measures of the local geometry, coordination, and electronic structure of Ca-ethereal complexes provide distinct structural trends depending on the extent of the Ca2+-solvent interaction. Finally, we correlate these findings with electrochemical measurements of calcium tetrakis(hexafluoroisopropoxy)borate (CaBHFIP2) salts dissolved within this class of solvents to provide insight into the preferred structural configuration of Ca2+ electrolytic solutions for optimized electrochemical plating and stripping.
Sandia National Laboratories’ (Sandia’s) Electrical Sciences Group (1350) provides a spectrum of solutions for electromagnetic (EM) environment effects on electrical systems to advance physical understanding including experiments, analytical and numerical model development leveraging first-principles physics, and high-performance computational modeling and simulation (M&S). The six departments in the Electrical Sciences Group provide analysis, design guidance, and experiments in support of nuclear weapons qualification in normal, abnormal, and hostile environments, as well as research and development for advanced electrical systems that can operate through these environments. This document is intended to provide a quick summary of Sandia’s EM experimental facilities, most of which provide unique national capabilities.
By 2030 about half of all spent nuclear fuel (SNF) arising from the current fleet of commercial power plants will be in dual-purpose canisters (DPCs), which are designed for storage and transportation but not for disposal. As an alternative to complete repackaging of the fuel for disposal, considerable cost savings and lower worker dose could be realized by directly disposing of this SNF in DPCs. The principal technical consideration is criticality control in a geologic repository, because the DPCs are large and depend on neutron absorbing basket components for criticality control. Neutron absorbing materials are generally aluminum-based, and under disposal conditions can degrade after a few hundred years contact with ground water. Simple modifications to the SNF assemblies or the DPC baskets could help to achieve direct disposal, and this is one of the approaches being studied to address the possibility of disposal criticality (SNL 2020a). Five fuel/basket modification concepts have been proposed (SNL 2020b) and a virtual workshop was conducted to solicit review and feedback on these concepts. The proposed solutions are: 1) zone loading of DPCs to limit reactivity, 2) replacing absorber plates with advanced neutron absorbing (ANA) material, 3) adding disposal control rods to pressurized water reactor (PWR) assemblies, 4) rechanneling boiling water reactor (BWR) assemblies with ANA material, and 5) basket insert plates (chevron inserts) made from ANA material. The presentations from the workshop are provided in this report, and the workshop discussions are summarized. This information includes prioritization of the proposed fuel/basket modification solutions, and prioritization of the associated model development, validation testing, and quality assurance activities. Information documented in this report will help to steer research and development efforts at Sandia National Laboratories, Oak Ridge National Laboratory, and Idaho National Laboratory that support the U.S. Department of Energy, Office of Nuclear Energy, Spent Fuel and Waste Science and Technology program
Alemazkoor, Negin; Rachunok, Benjamin; Chavas, Daniel R.; Staid, Andrea S.; Louhghalam, Arghavan; Nateghi, Roshanak; Tootkaboni, Mazdak
Nine in ten major outages in the US have been caused by hurricanes. Long-term outage risk is a function of climate change-triggered shifts in hurricane frequency and intensity; yet projections of both remain highly uncertain. However, outage risk models do not account for the epistemic uncertainties in physics-based hurricane projections under climate change, largely due to the extreme computational complexity. Instead they use simple probabilistic assumptions to model such uncertainties. Here, we propose a transparent and efficient framework to, for the first time, bridge the physics-based hurricane projections and intricate outage risk models. We find that uncertainty in projections of the frequency of weaker storms explains over 95% of the uncertainty in outage projections; thus, reducing this uncertainty will greatly improve outage risk management. We also show that the expected annual fraction of affected customers exhibits large variances, warranting the adoption of robust resilience investment strategies and climate-informed regulatory frameworks.
This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.
Hot-electron generation has been a topic of intense research for decades for numerous applications ranging from photodetection and photochemistry to biosensing. Recently, the technique of hot-electron generation using non-radiative decay of surface plasmons excited by metallic nanoantennas, or meta-atoms, in a metasurface has attracted attention. These metasurfaces can be designed with thicknesses on the order of the hot-electron diffusion length. The plasmonic resonances of these ultrathin metasurfaces can be tailored by changing the shape and size of the meta-atoms. One of the fundamental mechanisms leading to generation of hot-electrons in such systems is optical absorption, therefore, optimization of absorption is a key step in enhancing the performance of any metasurface based hot-electron device. Here we utilized an artificial intelligence-based approach, the genetic algorithm, to optimize absorption spectra of plasmonic metasurfaces. Using genetic algorithm optimization strategies, we designed a polarization insensitive plasmonic metasurface with 90% absorption at 1550 nm that does not require an optically thick ground plane. We fabricated and optically characterized the metasurface and our experimental results agree with simulations. Finally, we present a convolutional neural network that can predict the absorption spectra of metasurfaces never seen by the network, thereby eliminating the need for computationally expensive simulations. Our results suggest a new direction for optimizing hot-electron based photodetectors and sensors.
Suri, Pranav K.; Nathaniel, James E.; Li, Nan; Baldwin, Jon K.; Wang, Yongqiang; Hattar, Khalid M.; Taheri, Mitra L.
Gold is a noble metal typically stable as a solid in a face-centered cubic (FCC) structure under ambient conditions; however, under particular circumstances aberrant allotropes have been synthesized. In this work, we document the phase transformation of 25 nm thick nanocrystalline (NC) free-standing gold thin-film via in situ ion irradiation studied using atomic-resolution transmission electron microscopy (TEM). Utilizing precession electron diffraction (PED) techniques, crystallographic orientation and the radiation-induced relative strains were measured and furthermore used to determine that a combination of surface and radiation-induced strains lead to an FCC to hexagonal close packed (HCP) crystallographic phase transformation upon a 10 dpa radiation dose of Au4+ ions. Contrary to previous studies, HCP phase in nanostructures of gold was stabilized and did not transform back to FCC due to a combination of size effects and defects imparted by damage cascades.
Tools are now available that enable measurement electromagnetic radiation (EMR) from active electronics in an item. This radiation may be intended WIFI or cellular network links, for example or unintended such as the switching noise generated by DC-to-DC converters. It would be extremely valuable to have the capability to discriminate between the low-voltage DC-to-DC converters or other digital noise prevalent in most modern electronics, versus the high-voltage DC-to-DC converters used in utility firesets. Previous work performed under a Sandia Laboratory Directed Research and Development (LDRD) project on Charge State Detection using a deep neural network has been continued in this effort. A state-of-the-art supervised machine learning algorithm has not only been extended to discriminate between low and high voltage converters but has been validated in determining a converters make and model.
The vertical borehole array at Farnsworth Unit, TX is used to monitor microseismic activity in the subsurface around the Carbon Capture and Sequestration (CCS) reservoir. The array consists of 16 3-component seismometers spaced vertically in a single borehole. Tube or borehole waves traveling up or down the borehole can corrupt signals of interest, such as microseismic events. A denoising convolutional neural network (DCNN) was trained to remove borehole waves from seismic waveforms of microseismic events for the purpose of reducing unwanted signal detections and better characterizing events of interest. This R&D leverages the work of Sandia colleague Rigo Tibi, who used a DCNN developed by Greg Beroza at Stanford University to improve the signal-to-noise ratio (SNR) of teleseismic events detected by the International Monitoring System.
This report summarizes the key contributions and lessons learned from SNL experience in technical reviews of Controls awardees in the DOE SPA program from 2013 - 2020. The purpose of this report is to provide observations and technical suggestions that are likely to be beneficial to the WEC industry as a whole. Over the course of the SPA FOA program, SNL has engaged in technical review for a total of 5 different Controls awardees. The awardees represent a diversity of WEC devices and the application of different control design approaches. The report begins with a summary of key performance metrics results reported by the 5 Controls awardees. This is followed by a summary of observations and lessons learned distilled from the technical reviews of the awardees . The report concludes with a list of general technical suggestions for future WEC controls projects.
Natural and induced fractures are potential preferential pathways for migration of radioactive gases to earths surface from underground nuclear explosions (UNEs). This report documents X-ray computed tomography (XRCT) imaging on 26 samples of rock core that was collected to support the Underground Nuclear Explosion Signatures Experiment (UNESE) program. The XRCT datasets are intended to help fill a data gap on the three-dimensional (3D) characteristics of natural and/or induced fractures at the centimeter and smaller scale, which may strongly influence multiphase fluid flow and transport properties of preferential flow paths and interaction with the matrix of the surrounding host rock. Pre- and post-UNE rock samples were carefully chosen to enable comparison of fractures as a function of lithologic and petrophysical properties, as well as distance to the past UNEs. This report serves as documentation for the data, including an introduction with the research motivation, a methods and materials section, descriptions of the XRCT datasets without post-processing, and recommendations for 3D quantification via image analysis and digital rock physics.
Vaccination and the alternative behavior, vaccine refusal, are a classic example of manifesting behaviors driven by social norms and norm violations. Establishing how norms emerge, and under what circumstances people choose to violate them are key issues to understand in modeling epidemics. Interactions between individuals can lead to large-scale patterning of behavior (emergent phenomena). As norm violations are revealed through human behavior, drawing on psychological theory and principles to predict those violations is a viable approach for more human-constrained epidemiological models. As an example of the implications at scale, vaccine refusal is correlated with the spread of mis/disinformation about vaccine side-effects. Considering the complexities of network dynamics, the downstream effects means that if even a small group within a population are persuaded against vaccination, there is a reservoir from which disease and disease outbreaks can propagate. This work will attempt to identify those psychological indicators, to define circumstances that predict health behaviors, and identify potentially modifiable antecedents of health behavior, and factors that influence changes toward health protective behaviors.
This report summarizes design optimization and performance evaluation studies for a new prototype CONFIDANTE (CONfirmation using a Fast-neutron Imaging Detector with Anti-image Null-positive Time Encoding) warhead confirmation system. It was found that a spherical mask geometry and a 2” diameter cylindrical central detector is expected to best discriminate between the parametrically varied source distributions that were evaluated. The optimized design as fabricated and its performance was evaluated in a series of laboratory measurements. The performance was in good agreement with the design studies, with demonstrated discrimination between objects with differences in scale on the order of 5 centimeters or better at 1 meter stand-off
In this paper, two modifications are introduced for improving the accuracy, versatility, and robustness of a class of hybrid methods for radiation transport. In general, such methods are constructed by splitting the radiative flux into collided and uncollided components to which low- and high-resolution angular approximations are applied, respectively. In this work we focus on discrete ordinates discretizations of high and low order. The first modification we introduce changes the way in which the collided component is mapped into the uncollided component at the end of each time step in a simulation. The new mapping is a Nyström-type reconstruction that is applicable to arbitrary discrete ordinates quadratures, is guaranteed to preserve positivity of the solution provided that all ordinate weights are positive, is significantly more accurate than previous methods, and can be readily extended to other discretizations such as moment methods, finite element methods, and diffusion approximations. The second modification leverages integral deferred correction (IDC) to iteratively correct for the splitting error introduced by the inconsistency in angular discretization between the collided and uncollided components, in addition to improving the accuracy of the low-order temporal error that is treated by traditional IDC methods. Numerical tests in one- and two-dimensional geometries are used to demonstrate the increased accuracy and efficiency of the proposed modifications. It is found that the two techniques combined yield methods with solution accuracy and memory requirements comparable to that of monolithic discrete ordinates methods while reducing runtime by as much as a factor of between two and ten, depending on the problem.
Dispensers are the top cause of maintenance events and down-time at hydrogen fueling stations. In an effort to help characterize and enable improvements in dispenser reliability, an extensive accelerated lifetime testing set-up was designed and built at NREL involving components typically part of dispensing operations at fueling stations. Device Under Test (DUTs) included different components such as normally open valves, normally closed valves, fueling nozzles, breakaways devices and filters. Conditions of testing included pressures, and flow rates similar to light duty fuel cell electric vehicles fueling at -40°C, and -20°C for thousands of cycles in hydrogen. Tested components (failed and non-failed) were disassembled at SNL and polymeric O-rings were carefully retrieved and cataloged for chemical and physical characterization. Data collected was compared to similar O-rings from unexposed or non-tested components for hydrogen effects, and failure modes. Degradation analyses, based on select polymer chemistries common across all component types, their location within components, visual assessment of damage coupled with strong hydrogen effects from chemical characterization, was completed and presented to NREL and DOE. Overall, the failure rate amongst the components was not as high as expected for the test conditions. Among the component types tested, breakaways were the most susceptible to damage under these test conditions, with fueling nozzles a close second. The proper combination of selection of the right polymer and optimum component design was found to make a strong difference in component reliability under severe dispenser operating conditions. Physical degradation of polymers, rather than chemical changes due to low temperature hydrogen exposure, is more prevalent as failure mode for these test conditions. The nature and the extent of the degradation was much less at -20°C as compared to -40°C. The damage and failure rates were higher at lower temperatures than at higher test temperatures. As expected, increasing the number of cycles at the lowest test temperature (-40°C) increased damage. This indicates that cycling at the low temperature of -40°C required by SAE J2601 can reduce component life in fuel dispensing operations
Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Our research question is, "Microscopically, how does the hydration/dehydration process occur?"
Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Microscopically, how does the hydration/dehydration process occur?
This paper presents a topology optimization formulation for frequency-domain dynamics to reduce solution dependence upon initial guess and considered loading conditions. Due to resonance phenomena in undamped steady-state dynamics, objectives measuring dynamic response possess many local minima that may represent poor solutions to a design problem, an issue exacerbated for design with respect to multiple frequencies. We propose an extension of the modified error-in-constitutive-equations (MECE) method, used previously in material identification inverse problems, as a new approach for frequency-domain dynamics topology optimization to mitigate these issues. The main idea of the proposed framework is to incorporate an additional penalty-like term in the objective function that measures the discrepancy in the constitutive relations between stresses and strains and between inertial forces and displacements. Then, the design problem is cast within a PDE-constrained optimization formulation in which we seek displacements, stresses, inertial forces, and a density-field solution that minimize our new objective subject to conservation of linear momentum plus some additional constraints. We show that this approach yields superior designs to conventional gradient-based optimization approaches that solely use a functional of displacements as the objective, while strictly enforcing the constitutive equations. The MECE strategy integrates into a density-based topology optimization scheme for void–solid or two-phase material structural design. We highlight the merits of our approach in a variety of scenarios for direct frequency response design, considering multiple frequency load cases and structural objectives.
Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.
The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO2 has been captured in a tamped Richtmyer–Meshkov instability (RMI) environment at a pressure regime of 4–10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker–Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.
In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.
Grid-level energy storage systems are needed to enable intermittent renewables. Li-ion, Pb-acid battery systems have been implemented but pose safety and environmental risks. Successful grid storage must be safe, reliable, and low-cost.
Grid-level energy storage systems are needed to enable intermittent renewables. Li-ion, Pb-acid battery systems have been implemented but pose safety and environmental risks. Successful grid storage must be safe, reliable, and low-cost.
The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.
Cryogenic environments make superconducting computing possible by reducing thermal noise, electrical resistance and heat dissipation. Heat generated by the electronics and thermal conductivity of electrical transmission lines to the outside world constitute two main sources of thermal load in such systems. As a result, higher data rates require additional transmission lines which come at an increasingly higher cooling power cost. Hybrid or monolithic integration of silicon photonics with the electronics can be the key to higher data rates and lower power costs in these systems. We present a 4-channel wavelength division multiplexing photonic integrated circuit (PIC) built from modulators in the AIM Photonics process development kit (PDK) that operate at 25 Gbps at room temperature and 10 Gbps at 40 K. We further demonstrate 2-channel operation for 20 Gbps aggregate data rate at 40 K using two different modulators/wavelengths, with the potential for higher aggregate bit rates by utilizing additional channels.
This report provides an analysis of the clad barrier function associated with the direct disposal of dual purpose canisters (DPCs) under hypothetical conditions in a shale repository and in an alluvial repository, including the effect of a postulated criticality event inside a disposed DPC. Should a postulated criticality event occur in a hypothetical shale repository, cladding will primarily degrade by general corrosion. Stress corrosion cracking, hydride cracking, creep failure, pitting and crevice corrosion, rod pressurization, and clad unzipping are calculated to have little impact on cladding persistence. At the higher temperature expected during a postulated criticality event in a saturated shale repository, general corrosion of cladding would be rapid - on the order of 0.034 microns/yr. A few hundred years after onset of a postulated criticality event in a shale repository complete general corrosion of fuel assembly grid spacer walls and guide tubes will likely result in settling of fuel rods upon each other. This rod consolidation should displace the water moderator and possibly terminate a postulated criticality. The primary potential degradation pathway for cladding in a hypothetical alluvial repository is localized corrosion by fluoride, which cannot occur in a shale repository. Fluoride-enhanced corrosion of cladding would be accelerated under the slightly higher (< 100°C) temperatures associated with a postulated criticality event. The impact of criticality in both cases (shale and alluvial) would be to increase the amount of failed cladding. But it would require very specialized transport pathways.
This report provides a high-level test plan for deploying three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUHOMS Advanced Horizontal Storage Modules (AHSM) from Orano (formerly Transnuclear Inc.). The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or adapt in response to conflicting requirements. The goal of the testing is to collect highly defensible and detailed surface deposition measurements from the surface of dry storage systems in a marine coastal environment to guide chloride-induced stress corrosion crack (CISCC) research. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the thermal-hydraulic environment. Instrumentation throughout the canister, storage module, and environment will provide an extensive amount of information for the use of model validation. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer a high-fidelity quantification of the conditions experienced in a harsh yet realistic environment.
Sandia National Laboratories Radiation Protection department did not successfully address the causes leading to the issue NTS-SS-SNL-NMSITE-2017-0004, “Contamination Monitor Left in Service After Inadequate Calibration”, nor identify through the verification and validation process that corrective actions were not fully implemented and effective at addressing the issue(s). This report is the third causal in four years associated with the iPCM12 instrument and the Radiation Protection Instrumentation (RPI) organization within Radiation Protection (00628).
Tranchida, Julien G.; Dos Santos, Gonzalo; Aparicio, Romina; Linares, D.; Miranda, E.N.; Pastor, Gustavo M.; Bringa, Eduardo M.
In this paper, the magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics simulations coupled to molecular dynamics, using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M (T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M (T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering. Analytical approximations to the magnetization as a function of temperature and size are proposed.
Krisman, Alex; Meagher, Patrick; Zhao, Xinyu; Park, Ji-Woong; Lu, Tianfeng; Chen, Jacqueline H.
The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.
On June 5th, 2020 Sandia National Laboratories signed and recommended for closure the evidence package for the Noncompliance Tracking System NTS—SS-SNL-NMSITE-2017-0004, Contamination Monitor Left in Service After Inadequate Calibration. On July 16th, 2020 the Sandia Field Office sent a follow-up question list to the Radiation Protection Instrumentation Program Lead to clarify some points while reviewing the evidence package. The follow-up questions readily identified deficiencies with all three calibrations of the iPCM12 systems. On July 30th, 2020 Radiation Protection hosted a meeting to better understand the exact nature of the concerns, determine if the iPCM12 machines were out of calibration and therefore would have to be taken out of service. One result of this meeting was to determine if SNL could withdraw the evidence package. SNL transmitted a letter on August 26th, 2020 recalling the evidence package and committed to performing a thorough causal analysis and invoke the required issues management processes to ensure that corrective actions are sustained.
Tsyrenova, Ayuna; Farooq, Muhammad Q.; Anthony, Stephen M.; Mollaeian, Keyvan; Li, Yifan; Liu, Fei; Miller, Kyle; Ren, Juan; Anderson, Jared L.; Jiang, Shan
This study reveals the unique role on Janus particles of the solid-solid interface at the boundary in determining particle interactions and assembly. In an aqueous ionic liquid (IL) solution, Janus spheres adopt intriguing orientations with their boundaries pinned on the glass substrate. It was further discovered that the orientation was affected by the particle amphiphilicity as well as the chemical structure and concentration of the IL. Further characterization suggests that the adsorption on the hydrophilic side is due to both an electrostatic interaction and hydrogen bonding, while adsorption on the hydrophobic side is due to hydrophobic attraction. Through the concerted interplay of all these interactions, the amphiphilic boundary may attract an excessive amount of IL cations, which guide the unique orientations of the Janus spheres. The results highlight the importance of the Janus boundary that has not been recognized previously. Adsorption at the solid-solid interfaces may inspire new applications in areas such as separation and catalysis.
A new scatter calculation algorithm has been implemented in the Gamma Detector Response and Analysis Software (GADRAS) package that accounts for spectral effects of scattering materials not in the line of sight of the detector and the source. Previously, GADRAS would only apply scattering effects due to materials that fall between the source and detector. This new routine will allow better modeling of various scenarios including gamma imagers, collimated detectors, or traditional gamma detectors where scattering materials are present.
Producing and distributing COVID-19 vaccine during the pandemic is a major logistical challenge requiring careful planning and efficient execution. This report presents information on logistical, policy and technical issues relevant to rapidly fielding a COVID-19 vaccination program. For this study we (a) conducted literature review and subject matter expert elicitation to understand current vaccine manufacturing and distribution capabilities and vaccine allocation strategies, (b) designed a baseline vaccine distribution strategy and modeling strategy to provide insight into the potential for targeted distribution of limited initial vaccine supplies, and (c) developed parametric interfaces to enable vaccine distribution scenarios to be analyzed in depth with Sandias Adaptive Recovery Model that will allow us evaluate the additional sub- populations and alternative distribution scenarios from a public health benefit and associated economic disruption Principal issues, challenges, and complexities that complicate COVID-19 vaccine delivery identified in our literature and subject matter expert investigation include these items: The United States has not mounted an urgent nationwide vaccination campaign in recent history. The existing global manufacturing and distribution infrastructure are not able to produce enough vaccine for the population immediately. Vaccines, once available will be scarce resources. Prioritization for vaccine allocation will be built on existing distribution networks. Vaccine distribution may not have a universal impact on disease transmission and morbidity because of scarcity, priority population demographics, and underlying disease transmission rates. Considerations for designing a vaccine distribution strategy are discussed. A baseline distribution strategy is designed and tested using the Adaptive Recovery Model, which couples a deterministic compartmental epidemiological model and a stochastic network model. We show the impact of this vaccine distribution strategy on hospitalizations, mortality, and contact tracing requirements. This model can be used to quantitatively evaluate alternative distribution scenarios, guiding policy decisions as vaccine candidates are narrowed down.
I am a staff scientist at Sandia National Laboratories (SNL). I work on multiple fundamental-science projects and lead modeling/data-analysis for the stellar opacity experiments. I also work with theorists to refine plasma material-property calculations to make NNSA simulations more predictable. These are challenging and important problems for national security. In fact, our experiments raised questions about “opacity”, the property of matter that controls energy transport inside stars. This work was published by Nature due to its serious implication over broad applications, and our team received NNSA Defense Program Award of Excellence, both in 2015. Work at SNL is collegial and full of learning through interactions with the world’s finest experimentalists and theorists. When I started my Ph.D. program, working at a national lab was a dream position. Today, I collaborate with top scientists to solve mission-critical problems. The Stewardship Science Academic Programs (SSAP) and National Laser Users’ Facility (NLUF) played a vital role in the efficient learning and seamless transition to SNL.
Antz, Hartwig; Boman, Erik G.; Gates, Mark; Kruger, Scott; Li, Sherry; Loe, Jennifer A.; Osei-Kuffuor, Daniel; Tomov, Stan; Tsai, Yaohung M.; Meier Yang, Ulrike
The use of multiple types of precision in mathematical software has the potential to increase its performance on new heterogeneous architectures. The xSDK project focuses both on the investigation and development of multiprecision algorithms as well as their inclusion into xSDK member libraries. This report summarizes current efforts on including and/or using mixed precision capabilities in the math libraries Ginkgo, heFFTe, hypre, MAGMA, PETSc/TAO, SLATE, SuperLU, and Trilinos, including KokkosKernels. It contains both numerical results from libraries that already provide mixed precision capabilities, as well as descriptions of the strategies to incorporate multiprecision into established libraries.
Neuromorphic computing is a critical future technology for the computing industry, but it has yet to achieve its promise and has struggled to establish a cohesive research community. A large part of the challenge is that full realization of the potential of brain inspiration requires advances in both device hardware, computing architectures, and algorithms. This simultaneous development across technology scales is unprecedented in the computing field. This article presents a strategy, framed by market and policy pressures, for moving past these current technological and cultural hurdles to realize its full impact across technology. Achieving the full potential of brain-derived algorithms as well as post-complementary metal-oxide-semiconductor (CMOS) scaling neuromorphic hardware requires appropriately balancing the near-term opportunities of deep learning applications with the long-term potential of less understood opportunities in neural computing.
The equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50000 K, for densities between -ü=3 and 15g/cm3, with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.
The thirty-year non-stationary historical trends in the wave energy climate for United States coastal waters between 1980 and 2009 are investigated using spectral partitioned wave data generated from a WaveWatch III® (version 5.05) hindcast. In addition to historical trends in the omni-directional wave power, frequency and directionally resolved wave power, frequency and directional spreading, and seasonal variability, are examined for the first time, including their geographical distribution. These historical wave energy climate trends are linked to changes to the dominant wave systems and commensurate trends in the historical wind climate. Total wave power trends are consistent with other studies, but the present study identifies regions, and specific frequency and direction bands, where significant wave energy climate changes have. For some regions, changes to wave climate parameters exceeded one-percent annually, more than thirty-percent over the study period. Non-stationary trends of this magnitude have significant implications for ocean and coastal engineering projects designed assuming wave climates are stationary and warrant their consideration in design practices.
DeJong, Stephanie D.; Van Benthem, Mark H.; Keller, Timothy J.; Gillispie, Gregory D.
This work presents a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.
High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO), hydroperoxyalkyl peroxy (OOQOOH), and ketohydroperoxide (HOOPO, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. With the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO, OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO, OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O2 concentration, and temperature on the competition among radical decomposition and second O2 addition pathways.
Steady-state photocapacitance (SSPC) was conducted on nonpolar m-plane GaN n-type Schottky diodes to evaluate the defects induced by inductively coupled plasma (ICP) dry etching in etched-and-regrown unipolar structures. An ∼10× increase in the near-midgap Ec - 1.9 eV level compared to an as-grown material was observed. Defect levels associated with regrowth without an etch were also investigated. The defects in the regrown structure (without an etch) are highly spatially localized to the regrowth interface. Subsequently, by depth profiling an etched-and-regrown sample, we show that the intensities of the defect-related SSPC features associated with dry etching depend strongly on the depth away from the regrowth interface, which is also reported previously [Nedy et al., Semicond. Sci. Technol. 30, 085019 (2015); Fang et al., Jpn. J. Appl. Phys. 42, 4207-4212 (2003); and Cao et al., IEEE Trans. Electron Devices 47, 1320-1324 (2000)]. A photoelectrochemical etching (PEC) method and a wet AZ400K treatment are also introduced to reduce the etch-induced deep levels. A significant reduction in the density of deep levels is observed in the sample that was treated with PEC etching after dry etching and prior to regrowth. An ∼2× reduction in the density of Ec - 1.9 eV level compared to a reference etched-and-regrown structure was observed upon the application of PEC etching treatment prior to the regrowth. The PEC etching method is promising for reducing defects in selective-area doping for vertical power switching structures with complex geometries [Meyers et al., J. Electron. Mater. 49, 3481-3489 (2020)].
Machine learning (ML) techniques are being used to detect increasing amounts of malware and variants. Despite successful applications of ML, we hypothesize that the full potential of ML is not realized in malware analysis (MA) due to a semantic gap between the ML and MA communities-as demonstrated in the data that is used. Due in part to the available data, ML has primarily focused on detection whereas MA is also interested in identifying behaviors. We review existing open-source malware datasets used in ML and find a lack of behavioral information that could facilitate stronger impact by ML in MA. As a first step in bridging this gap, we label existing data with behavioral information using open-source MA reports-1) altering the analysis from identifying malware to identifying behaviors, 2)~aligning ML better with MA, and 3)~allowing ML models to generalize to novel malware in a zero/few-shot learning manner. We classify the behavior of a malware family not seen during training using transfer learning from a state-of-the-art model for malware family classification and achieve 57%-84% accuracy on behavioral identification but fail to outperform the baseline set by a majority class predictor. This highlights opportunities for improvement on this task related to the data representation, the need for malware specific ML techniques, and a larger training set of malware samples labeled with behaviors.
The swelling of clay at high temperature and pressure is important for applications including nuclear waste storage but is not well understood. A molecular dynamics study of the swelling of Na montmorillonite in water at several temperatures (T = 298, 400, and 500 K) and water environment pressures (Pe = 5 and 100 MPa) is reported here. Adopting a rarely used setup that enables swelling pressure to be resolved with an accuracy of ~1 MPa, the swelling pressure was computed at basal spacings 1.6–2.6 nm (or 2–5 water layers between neighboring clay sheets), which has not been widely studied before. At T = 298 K and Pe = 5 MPa, swelling pressure Ps oscillates at d-spacing d smaller than 2.2 nm and decays monotonically as d increases. Increasing T to 500 K but keeping Pe at 5 MPa, Ps remains oscillatory at small d, but its repulsive peak at d = 2.2 nm shifts to ~2.0 nm and Ps at different d-spacings can grow more attractive or repulsive. At d > 2.0 nm, Ps is weakened greatly. Keeping T at 500 K and increasing Pe to 100 MPa, Ps recovers toward that at T = 298 K and Pe = 5 MPa, however, the repulsive peak at d = 2.0 nm remains the same. The opposite effects of increasing temperature and pressure on the density and dielectric screening of water, which control ion correlations and thus double layer repulsion, are essential for understanding the observed swelling pressure at elevated temperatures and its response to environment pressures.