Compared to the classical Lanczos algorithm, the s-step Lanczos variant has the potential to improve performance by asymptotically decreasing the synchronization cost per iteration. However, this comes at a price; despite being mathematically equivalent, the s-step variant may behave quite differently in finite precision, potentially exhibiting greater loss of accuracy and slower convergence relative to the classical algorithm. It has previously been shown that the errors in the s-step version follow the same structure as the errors in the classical algorithm, but are amplified by a factor depending on the square of the condition number of the O(s)-dimensional Krylov bases computed in each outer loop. As the condition number of these s-step bases grows (in some cases very quickly) with s, this limits the s values that can be chosen and thus can limit the attainable performance. In this work, we show that if a select few computations in s-step Lanczos are performed in double the working precision, the error terms then depend only linearly on the conditioning of the s-step bases. This has the potential for drastically improving the numerical behavior of the algorithm with little impact on per-iteration performance. Our numerical experiments demonstrate the improved numerical behavior possible with the mixed precision approach, and also show that this improved behavior extends to mixed precision s-step CG. Here, we present preliminary performance results on NVIDIA V100 GPUs that show that the overhead of extra precision is minimal if one uses precisions implemented in hardware.
The objective of this report is to accept or reject the hypothesis that the experiments conducted under TP 08-02 Revision 0 (Ismail et al., 2008) were affected by CO2(g) intrusion and sample contamination. The test of the hypothesis is accomplished by comparing the experimental data collected under the protocols of TP 08-02 Revision O and TP 20-01 Revision O (Kirkes and Zhang, 2020). The protocols of TP 20-01 Revision 0 minimize the possibilities of CO2(g) intrusion and sample contamination. The experimental data sets obtained under both TPs will be assessed statistically to see if they are identical or not.
This report analyzes experimental data from Test Plans TP 08-02, TP 12-02, and TP 20-01 to add new log K values and Pitzer interaction parameters for Fe, Pb, Mg, Nd and B reactions to the WIPP geochemical thermodynamic database, data0.fm 1.
The Dynamic Networks Experiment 2018 (DNE18) was a collaborative effort between Los Alamos National Laboratory (LANL), Sandia National Laboratories (SNL), Lawrence Livermore National Laboratory (LLNL) and Pacific Northwest National Laboratory (PNNL) designed to evaluate methodologies for multi-modal data ingestion and processing. One component of this virtual experiment was a quantitative assessment of current capabilities for infrasound data processing, beginning with the establishment of a baseline for infrasound signal detection. To produce such baselines, SNL and LANL exploited a common dataset of infrasound data recorded across a regional network in Utah from December 2010 through February 2011. We utilize two automated signal detectors, the Adaptive F-Detector (AFD) and the Multivariate Adaptive Learning Detector (MALD) to produce automated signal detection catalogs and an analyst-produced catalog. Comparisons indicate that automatic detectors may be able to identify small amplitude, low SNR events that cannot be identified by analyst review. We document detector performance in terms of precision and recall, demonstrating that the AFD is more precise, but the MALD has higher recall. We use a synthetic dataset of signals embedded in pink noise in order to highlight shortcomings in assessing detection algorithms for low signal to noise ratio signals which are commonly of interest to the nuclear monitoring community. For comparisons utilizing the synthetic dataset, the AFD has higher recall while precision is equal for both detectors. These results indicate that both detectors perform well across a variety of background noise environments; however, both detectors fail to identify repetitive, short duration signals arriving from similar backazimuths. These failures represent specific scenarios that could be targeted for further detector development.
Characteristic modes on infinite periodic structures are studied using spectral dyadic Green’s functions. This formulation demonstrates that, in contrast to the modal analysis of finite structures, the number of radiating characteristic modes is limited by unit cell size and incident wave vector (i.e., scan angle or phase shift per unit cell). Here, the reflection tensor is decomposed into modal contributions from radiating modes, indicating that characteristic modes are a predictably sparse basis in which to study reflection phenomena.
The interplay of stress, disorder, and Coulomb screening dictating the mobility of doped cadmium oxide (CdO) is examined using Raman spectroscopy to identify the mechanisms driving dopant incorporation and scattering within this emerging infrared optical material. Specifically, multi-wavelength Raman and UV-vis spectroscopies are combined with electrical Hall measurements on a series of yttrium (X = Y) and indium (X = In) doped X:CdO thin-films. Hall measurements confirm n-type doping and establish carrier concentrations and mobilities. Spectral fitting along the low-frequency Raman combination bands, especially the TA+TO(X) mode, reveals that the evolution of strain and disorder within the lattice as a function of dopant concentration is strongly correlated with mobility. Coupling between the electronic and lattice environments was examined through analysis of first- and second-order longitudinal-optical phonon-plasmon coupled modes that monotonically decrease in energy and asymmetrically broaden with increasing dopant concentration. By fitting these trends to an impurity-induced Fröhlich model for the Raman scattering intensity, exciton-phonon and exciton-impurity coupling factors are quantified. These coupling factors indicate a continual decrease in the amount of ionized impurity scattering with increasing dopant concentration and are not as well correlated with mobility. This shows that lattice strain and disorder are the primary determining factors for mobility in donor-doped CdO. In aggregate, the study confirms previously postulated defect equilibrium arguments for dopant incorporation in CdO while at the same time identifying paths for its further refinement.
High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.
Borehole cement is used across the range of energy technologies to stabilize casing, to serve as a barrier to behind-casing fluid movement. Cement debonding and other flaws, both at cement interfaces and within the cement itself, can create leakage pathways that can threaten safety to personnel, and wellbore performance, with economic and regulatory consequences. A primary method to assess cement health and wellbore integrity is via acoustic methods. This project was designed with three aims: demonstrate a significant improvement in the interpretation of cement acoustic behavior, both during curing, and in interpreting effects of flaws and evolving interfaces; develop sensor technologies to improve signal-noise ratios and cement acoustic responses; and lastly, provide a borehole demonstration of at least one of these technologies. We have accomplished the first two objectives, and the third, delayed by pandemic health concerns, is proceeding as of this writing via a technology partner with the University of Texas Advanced Energy Consortium.
Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na+ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.
This project is intended to support the development of new traction drive systems that meet the targets of 100 kW/L for power electronics and 50 kW/L for electric machines with reliable operation to 300,000 miles. To meet these goals, new designs must be identified that make use of state-of-the-art and next-generation electronic materials and design methods. Designs must exploit synergies between components, for example converters designed for high-frequency switching using wide band gap devices and ceramic capacitors. This project includes: (1) a survey of available technologies; (2) the development of design tools that consider the converter volume and performance; (3) exercising the design software to evaluate performance gaps and predict the impact of certain technologies and design approaches, i.e. GaN semiconductors, ceramic capacitors, and select topologies; and (4) building and testing hardware prototypes to validate models and concepts. Early instantiations of the design tools enable co-optimization of the power module and passive elements and provide some design guidance; later instantiations will enable the co-optimization of inverter and machine. Prototype testing begins with evaluation of simpler conversion topologies (i.e. the half-bridge boost converter) and progresses with fabrication of prototype inverter drives.
This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. This study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.
The global energy system is undergoing significant changes, including a shift in energy generating technologies to more renewable energy sources. However, the dependence of renewable energy sources on local environmental conditions could also increase disruptions in service through exposures to compound, extreme weather events. By fusing three diverse datasets (operations and maintenance tickets, weather data, and production data), this analysis presents a novel methodology to identify and evaluate performance impacts arising from extreme weather events across diverse geographical regions. Text analysis of maintenance tickets identified snow, hurricanes, and storms as the leading extreme weather events affecting photovoltaic plants in the United States. Statistical techniques and machine learning were then implemented to identify the magnitude and variability of these extreme weather impacts on site performance. Impacts varied between event and non-event days, with snow events causing the greatest reductions in performance (54.5%), followed by hurricanes (12.6%) and storms (1.1%). Machine learning analysis identified key features in determining if a day is categorized as low performing, such as low irradiance, geographic location, weather features, and site size. This analysis improves our understanding of compound, extreme weather event impacts on photovoltaic systems. These insights can inform planning activities, especially as renewable energy continues to expand into new geographic and climatic regions around the world.
Novel materials based on the aluminum oxyhydroxide boehmite phase were prepared using a glycothermal reaction in 1,4-butanediol. Under the synthesis conditions, the atomic structure of the boehmite phase is altered by the glycol solvent in place of the interlayer hydroxyl groups, creating glycoboehmite. The structure of glycoboehmite was examined in detail to determine that glycol molecules are intercalated in a bilayer structure, which would suggest that there is twice the expansion identified previously in the literature. This precursor phase enables synthesis of two new phases that incorporate either polyvinylpyrrolidone or hydroxylpropyl cellulose nonionic polymers. These new materials exhibit changes in morphology, thermal properties, and surface chemistry. All the intercalated phases were investigated using PXRD, HRSTEM, SEM, FT-IR, TGA/DSC, zeta potential titrations, and specific surface area measurement. These intercalation polymers are non-ionic and interact through wetting interactions and hydrogen bonding, rather than by chemisorption or chelation with the aluminum ions in the structure.
International Conference for High Performance Computing, Networking, Storage and Analysis, SC
Mullowney, Paul; Li, Ruipeng; Thomas, Stephen; Ananthan, Shreyas; Sharma, Ashesh; Rood, Jon S.; Williams, Alan B.; Sprague, Michael A.
The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices the result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.
This work focuses on the space-time reduced-order modeling (ROM) method for solving large-scale uncertainty quantification (UQ) problems with multiple random coefficients. In contrast with the traditional space ROM approach, which performs dimension reduction in the spatial dimension, the space-time ROM approach performs dimension reduction on both the spatial and temporal domains, and thus enables accurate approximate solutions at a low cost. We incorporate the space-time ROM strategy with various classical stochastic UQ propagation methods such as stochastic Galerkin and Monte Carlo. Numerical results demonstrate that our methodology has significant computational advantages compared to state-of-the-art ROM approaches. By testing the approximation errors, we show that there is no obvious loss of simulation accuracy for space-time ROM given its high computational efficiency.
The performance of solid-state electrochemical systems is intimately tied to the potential and lithium distributions across electrolyte-electrode junctions that give rise to interface impedance. Here, we combine two operando methods, Kelvin probe force microscopy (KPFM) and neutron depth profiling (NDP), to identify the rate-limiting interface in operating Si-LiPON-LiCoO2 solid-state batteries by mapping the contact potential difference (CPD) and the corresponding Li distributions. The contributions from ions, electrons, and interfaces are deconvolved by correlating the CPD profiles with Li-concentration profiles and by comparisons with first-principles-informed modeling. We find that the largest potential drop and variation in the Li concentration occur at the anode-electrolyte interface, with a smaller drop at the cathode-electrolyte interface and a shallow gradient within the bulk electrolyte. Correlating these results with electrochemical impedance spectroscopy following battery cycling at low and high rates confirms a long-standing conjecture linking large potential drops with a rate-limiting interfacial process.
The Galactic/Extra Ultra long Duration Balloon Spectroscopic-Stratospheric Terahertz Observatory (GUSTO), is a NASA balloon-borne project and is scheduled for launch in late 2022. The balloon will carry a spectroscopic telescope that will detect three brightest emission lines from interstellar medium. GUSTO measurements will shed light on the life-cycle of the gases in the Milky Way and Large Magellanic Cloud (LMC). In this study, we will discuss the details of a quantum cascade laser used in the local oscillator for detecting the oxygen line at 4.74 THz.
Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain Monte Carlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density ne and ion temperature Ti, since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the ne-Ti ambiguity. The measured line shapes indicate densities that are 80-100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. The experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.
Batteries are an enabling technology for addressing sustainability through the electrification of various forms of transportation (1) and grid storage. (2) Batteries are truly multi-scale, multi-physics devices, and accordingly various theoretical descriptions exist to understand their behavior (3-5) ranging from atomistic details to techno-economic trends. As we explore advanced battery chemistries (6,7) or previously inaccessible aspects of existing ones, (8-10) new theories are required to drive decisions. (11-13) The decisions are influenced by the limitations of the underlying theory. Advanced theories used to understand battery phenomena are complicated and require substantial effort to reproduce. However, such constraints should not limit the insights from these theories. We can strive to make the theoretical research verifiable such that any battery stakeholder can assess the veracity of new theories, sophisticated simulations or elaborate analyses. We distinguish verifiability, which amounts to “Can I trust the results, conclusions and insights and identify the context where they are relevant?”, from reproducibility, which ensures “Would I get the same results if I followed the same steps?” With this motivation, we propose a checklist to guide future reports of theoretical battery research in Table 1. We hereafter discuss our thoughts leading to this and how it helps to consistently document necessary details while allowing complete freedom for creativity of individual researchers. Given the differences between experimental and theoretical studies, the proposed checklist differs from its experimental counterparts. (14,15) This checklist covers all flavors of theoretical battery research, ranging from atomic/molecular calculations (16-19) to mesoscale (20,21) and continuum-scale interactions, (9,22) and techno-economic analysis. (23,24) Finally, as more and more experimental studies analyze raw data, (25) we feel this checklist would be broadly relevant.
An intriguing new class of two-dimensional (2D) materials based on metal-organic frameworks (MOFs) has recently been developed that displays electrical conductivity, a rarity among these nanoporous materials. The emergence of conducting MOFs raises questions about their fundamental electronic properties, but few studies exist in this regard. Here, we present an integrated theory and experimental investigation to probe the effects of metal substitution on the charge transport properties of M-HITP, where M = Ni or Pt and HITP = 2,3,6,7,10,11-hexaiminotriphenylene. The results show that the identity of the M-HITP majority charge carrier can be changed without intentional introduction of electronically active dopants. We observe that the selection of the metal ion substantially affects charge transport. Using the known structure, Ni-HITP, we synthesized a new amorphous material, a-Pt-HITP, which although amorphous is nevertheless found to be porous upon desolvation. Importantly, this new material exhibits p-type charge transport behavior, unlike Ni-HITP, which displays n-type charge transport. These results demonstrate that both p- and n-type materials can be achieved within the same MOF topology through appropriate choice of the metal ion.
A new scaling is derived that yields a Reynolds-number-independent profile for all components of the Reynolds stress in the near-wall region of wall-bounded flows, including channel, pipe and boundary layer flows. The scaling demonstrates the important role played by the wall shear stress fluctuations and how the large eddies determine the Reynolds number dependence of the near-wall turbulence behaviour.
A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.
In this study, several uncertainty quantification and sensitivity analysis methods are used to determine the most sensitive geometric and material input parameters of a cantilevered pipeline conveying fluid when uncertainty is introduced to the system at the onset of instability. The full nonlinear equations of motion are modeled using the extended Hamilton’s principle and then discretized using Galerkin’s method. A parametric study is first performed, and the Morris elementary effects are calculated to obtain a preliminary understanding of how the onset speed changes when each parameter is introduced to a ± 5% uncertainty. Then, four different input uncertainty distributions, mainly, uniform and Gaussian distribution, are chosen to investigate how input distributions affect uncertainty in the output. A convergence analysis is used to determine the number of samples needed to maintain simulation accuracy while saving the most computational time. Then, Monte Carlo simulations are run, and the output distributions for each input distribution at ± 1%, ± 3% and ± 5% input uncertainty range are found and discussed. Additionally, the Pearson correlation coefficients are evaluated for different uncertainty ranges. A final Monte Carlo study is performed in which single parameters are held constant while all others still have uncertainty. Overall, the flow speed at the onset of instability is the most sensitive to changes in the outer diameter of the pipe.
For strongly charged polyelectrolytes in salt-free solutions, we use molecular dynamics simulations of a coarse-grained bead-spring model to calculate overlap concentrations c∗ and chain structure for polymers containing N = 10 to 1600 monomers. Over much of this range, we find that the end-to-end distance R∗ at c∗ increases faster than linearly with increasing N, as chains at the overlap concentration approach strongly extended conformations. This trend results in the overlap concentration c∗ decreasing as a stronger function of N than the classical prediction c∗ ∼N-2. This stronger dependence can be fit either by a logarithmic correction to scaling or by an apparent scaling c∗ ∼N-m, with m > 2.
Dalbey, Keith R.; Eldred, Michael S.; Geraci, Gianluca; Jakeman, John D.; Maupin, Kathryn A.; Monschke, Jason A.; Seidl, Daniel T.; Tran, Anh; Menhorn, Friedrich; Zeng, Xiaoshu
The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.
This report describes research and development (R&D) activities conducted during fiscal year 2021 (FY21) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). The R&D activities focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Generic Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc.
There are several different calculation approaches and tools that can be used to evaluate the risk of hydrogen energy applications. A comparative study of Air Liquide’s ALDEA (Air Liquide Dispersion and Explosion Assessment) tools suite and Sandia’s HyRAM (Hydrogen Risk Assessment Models) toolkit has been conducted. The purpose of this study was to understand and evaluate the differences between the two calculation approaches, and identify areas for model improvements. There were several scenarios examined in this effort regarding hydrogen release dynamics. These scenarios include free jet release cases at varying pressures, vessel blowdown, and hydrogen build-up scenarios with and without ventilation. For each scenario, the input and output of the HyRAM calculations are documented, along with a comparison to the ALDEA results. Generally, the results from the two different tools were reasonably aligned. However, there were fundamental differences in evaluation methodology and functional limitations in HyRAM that caused discrepancies in some calculations.
Ninad, Nayeem; Apablaza-Arancibia, Estefan; Bui, Michel; Johnson, Jay
As more countries seek solutions to their de-carbonization targets using renewable energy (RE) technologies, interconnection standards and national grid codes for distributed energy resources (DER) are being updated to support higher penetrations of RE and improve grid stability. Common grid-code revisions mandate DER devices, such as solar inverters and energy storage systems, ride-through (RT) voltage and frequency disturbances. This is necessary because as the percentage of generation from DER increases, there is a greater risk power system faults will cause many or all DER to trip, triggering a substantial load-generation imbalance and possible cascading blackout. This paper demonstrates for the first time a methodology to verify commercial DER devices are compliant to new voltage, frequency, and rate of change of frequency (ROCOF) RT requirements established in IEEE Std. 1547-2018. The methodology incorporates a software automation tool, called the SunSpec System Validation Platform (SVP), in combination with a hardware-in-the-loop (HIL) system to execute the IEEE Std. 1547.1-2020 RT test protocols. In this paper, the approach is validated with two commercial photovoltaic inverters, the test results are analyzed for compliance, and improvements to the test procedure are suggested.
We developed a simplistic physics-based model of an all-optical neural network that mimics the encoder part of an autoencoder neural network for image compression. Our approach relies on the generation of a MATLAB-based model for both data compression and decompression and utilizes MATLAB's built-in autoencoder networks in combination with simple propagation of optical fields between layers constituting phase elements via Fourier transform. We optimize the phase elements using the particle swarm optimization technique and using our model, we demonstrate a compression ratio of 25% for 2828-pixel input images containing numeric digits from 0 to 9.
The tearing parameter criterion and failure propagation method currently used in the multilinear elastic-plastic constitutive model was added as an option to modular failure capabilities. Currently, this implementation is only available to the J2 plasticity model due to the formulation of the failure propagation approach. The implementation was verified against analytical solutions for both a uniaxial tension and a pure shear boundary-value problem. Possible improvements to, and necessary generalizations of, the failure method to extend it as a modular option for all plasticity models are highlighted.
The Sandia Solution Mining Code (SANSMIC) has been used for many years to examine the development of salt cavern geometry, both in a confirmatory manner with comparisons made to real-world sonar data and in a predictive manner when updated sonar data are not available. SANSMIC models require some modeling choices in order to incorporate real-world data. Key modeling choices include the vertical resolution of cavern geometry to implement, as well as how to incorporate daily raw water injection data into the SANSMIC model. This report documents five studies that address the impact of the modeling choices on the predicted cavern geometries and calculated leaching efficiencies. In most cases, hypothetical cylindrical initial cavern geometries are used to provide a common baseline against which to test the systematic variation of input variables including cavern radius, oil-brine-interface (OBI) depth, vertical cell size, raw water injection rate, raw water injection duration, workover time, and number of leaching stages. The use of smaller cell sizes is recommended moving forward to provide a better one-to-one relationship between sonar data and the modeled cavern. A new methodology for incorporating raw water injection data is also recommended, in order to more closely model real-world injection and workover times. Overall, the systematic studies performed here have increased our confidence in previous SANSMIC model results, as well future use of the code for predicting leaching effects on cavern geometries. Some minor changes to modeling choices are recommended, which can easily be applied with the version of SANSMIC currently under development.
This project will test the coupling of light emitted from silicon vacancy and nitrogen vacancy defects in diamond into additively manufactured photonic wire bonds toward integration into an "on-chip quantum photonics platform". These defects offer a room-temperature solid state solution for quantum information technologies but suffer from issues such as low activation rate and variable local environments. Photonic wire bonding will allow entanglement of pre-selected solid-state defects alleviating some of these issues and enable simplified integration with other photonic devices. These developments could prove to be key technologies to realize quantum secured networks for national security applications.
Tonopah Test Range (TTR), in support of its testing mission and modernization effort acquired a fleet of new gimballed tracking mounts (GTMs) manufactured by BAE Systems. The new GTMs can be operated remotely during flight tests and provide near real-time target tracking data. Furthermore, test vehicle Time-Space-Position-Information (TSPI) is evaluated using post-test synchronized imagery and pointing angle measurements acquired from each tracking mount. To comply with the Nuclear Enterprise Assurance Program (NEAP), all measurements devices must be certified. In keeping with the NEAP program, qualification of the new GTMs have been assessed to confirm that their pointing angle measurements produce acceptable TSPI results. This study only evaluated the four GTMs as a stand-alone solution and found that the GTMs meet their performance requirement of 0.006 degrees RMS error (or less) for post-processed pointing angles and produced TSPI solution with error volumes on the order of one meter or less. The new GTMs will be utilized in combination with existing optical tracking mounts, which will only improve the accuracy of the resulting TSPI data product. Details regarding the approach, analysis, summary results, and conclusions are presented.
A new liquid sample adapter design for the Explosive Destruction Systems has been developed. The design features a semi-transparent fluoropolymer tube coupled to the vessel high pressure sample valve with a closing quick connect fitting. The sample tubes are the pressure-limiting component. The tubes were hydrostatically tested to establish failure characteristics and pressure limits at ambient and operational temperatures. A group of tubes from two manufacturing lots were tested to determine the consistency of the commercial part. An upper pressure limit was determined for typical operations.
The TChem open-source software is a toolkit for computing thermodynamic properties, source term, and source term’s Jacobian matrix for chemical kinetic models that involve gas and surface reactions.
Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.
D'Elia, Marta; Howard, Amanda; Kirby, Michael R.; Kutz, Nathan; Tartakovsky, Alexandre M.; Viswanathan, Hari S.
A hallmark of the scientific process since the time of Newton has been the derivation of mathematical equations meant to capture relationships between observables. As the field of mathematical modeling evolved, practitioners specifically emphasized mathematical formulations that were predictive, generalizable, and interpretable. Machine learning’s ability to interrogate complex processes is particularly useful for the analysis of highly heterogeneous, anisotropic materials where idealized descriptions often fail. As we move into this new era, we anticipate the need to leverage machine learning to aid scientists in extracting meaningful, but yet sometimes elusive, relationships between observed quantities.
Heavy-Duty diesel engine manufacturers are continuously in pursuit of simple and low-cost technologies that can reduce emissions. Ducted fuel injection (DFI) and Cooled Spray (CS) technologies are two technologies that continue to show promise for significant particulate emissions reductions. These technologies represent a breakthrough in diesel engine combustion from the potential of nearly sootless diesel combustion. This can provide a significant decrease in harmful PM emissions and may enable further system optimization for reduced NOx emissions and increased efficiency. Combustion vessel experiments and engine demonstrations at Sandia, together with the large bore engine tests performed by Wabtec show that this technology may be applicable to heavy duty diesel engines across a wide range of engine sizes and speeds representing the majority of off-road diesel engines. However, very little is known about the ideal geometry, scaling properties or effectiveness of these technologies over the engine operating map. This project will address those uncertainties through a series of experiments performed in an optical and a metal single-cylinder engine.
Regulatory drivers and market demands for lower pollutant emissions, lower carbon dioxide emissions, and lower fuel consumption motivate the development of cleaner and more fuel-efficient engine operating strategies. Most current production heavy-duty diesel engines use a combination of both in-cylinder and exhaust emissions-control strategies to achieve these goals. The emissions and efficiency performance of in-cylinder strategies depend strongly on flow and mixing processes that can be influenced by using multiple fuel injections. Past work performed under this project showed that adding a second injection can reduce soot to levels below what would have been produced by an unchanged first injection, thereby increasing load while decreasing soot and potentially reducing brake specific fuel consumption. Information characterizing the important in-cylinder processes with multiple injections has been gleaned from ensemble-averaged planar laser-induced incandescence (PLII) imaging visualizing the soot cloud and planar induced fluorescence (PLIF) of OH characterizing the soot oxidation regions. PLII showed a consistent disruption of the first injection soot cloud by the second injection. In conjunction with OH-PLIF, differences in soot oxidation patterns for multiple injections compared to single injections were observed. This understanding was further enhanced in FY20, when high-speed imaging resolving the above-mentioned effects in a single cycle were combined with direct numerical simulations investigating the multiple-injection ignition process on the microscopic level of turbulence and chemistry interaction. In FY21, these findings in conjunction with findings from other researchers published in the scientific literature were composed into a preliminary multiple-injection conceptual model of fuel-mixing, injection and ignition processes. Remaining key research questions were also highlighted. In addition, wall heat flux was investigated experimentally and with numerical simulations to understand the potential of multiple injections to reduce the engine heat losses and further enhance the efficiency.
Graph algorithms enable myriad large-scale applications including cybersecurity, social network analysis, resource allocation, and routing. The scalability of current graph algorithm implementations on conventional computing architectures are hampered by the demise of Moore’s law. We present a theoretical framework for designing and assessing the performance of graph algorithms executing in networks of spiking artificial neurons. Although spiking neural networks (SNNs) are capable of general-purpose computation, few algorithmic results with rigorous asymptotic performance analysis are known. SNNs are exceptionally well-motivated practically, as neuromorphic computing systems with 100 million spiking neurons are available, and systems with a billion neurons are anticipated in the next few years. Beyond massive parallelism and scalability, neuromorphic computing systems offer energy consumption orders of magnitude lower than conventional high-performance computing systems. We employ our framework to design and analyze new spiking algorithms for shortest path and dynamic programming problems. Our neuromorphic algorithms are message-passing algorithms relying critically on data movement for computation. For fair and rigorous comparison with conventional algorithms and architectures, which is challenging but paramount, we develop new models of data-movement in conventional computing architectures. This allows us to prove polynomial-factor advantages, even when we assume a SNN consisting of a simple grid-like network of neurons. To the best of our knowledge, this is one of the first examples of a rigorous asymptotic computational advantage for neuromorphic computing.
In this project we endeavored to improve the state-of-the-art in UV lasers diodes. We made important advancements in several fronts from modeling, to epitaxial growth, to fabrication, and testing. Throughout the project it became clear that polarization doping would be able to help advance the state of laser diode design in terms of electrical performance, but the optical design would need to be investigated to ensure that a 2D guided mode would be supported. New capability in optical modeling using commercial software demonstrated that the new polarization doped structures would be viable. New capability in pulsed testing was established to reach the current and voltage required. Our fabricated devices had some parasitic electrical paths which hindered performance that we were ultimately unable to overcome in the project timeframe. We do believe that future projects will be able to leverage the advancements made under this project.
The need to reduce the carbon footprint from medium- and heavy-duty diesel engines is clear; low-carbon biofuels are a powerful means to achieve this. Liquid fuels are rapidly deployed because existing infrastructure can be utilized for their production, transport, and distribution. Their impact is unique as they can decrease the greenhouse gas (GHG) emissions of existing vehicles and in applications resistant to electrification. However, introducing new diesel-like bio-blends into the market is very challenging. At a minimum, it requires a comprehensive understanding of the life-cycle GHG emissions of the fuels, the implications for refinery optimization and economics, the fuel’s impact on the infrastructure, the effect on the combustion performance of current and future vehicle fleets, and finally the implications for exhaust aftertreatment systems and compliance with emissions regulations. Such understanding is sought within the Co-Optima project.
At the Z Facility at Sandia National Laboratories, the magnetized liner inertial fusion (MagLIF) program aims to study the inertial confinement fusion in deuterium-filled gas cells by implementing a three-step process on the fuel: premagnetization, laser preheat, and Z-pinch compression. In the laser preheat stage, the Z-Beamlet laser focuses through a thin polyimide window to enter the gas cell and heat the fusion fuel. However, it is known that the presence of the few μm thick window reduces the amount of laser energy that enters the gas and causes window material to mix into the fuel. These effects are detrimental to achieving fusion; therefore, a windowless target is desired. The Lasergate concept is designed to accomplish this by "cutting"the window and allowing the interior gas pressure to push the window material out of the beam path just before the heating laser arrives. In this work, we present the proof-of-principle experiments to evaluate a laser-cutting approach to Lasergate and explore the subsequent window and gas dynamics. Further, an experimental comparison of gas preheat with and without Lasergate gives clear indications of an energy deposition advantage using the Lasergate concept, as well as other observed and hypothesized benefits. While Lasergate was conceived with MagLIF in mind, the method is applicable to any laser or diagnostic application requiring direct line of sight to the interior of gas cell targets.
Adams, Brian M.; Bohnhoff, William J.; Dalbey, Keith R.; Ebeida, Mohamed S.; Eddy, John P.; Eldred, Michael S.; Hooper, Russell W.; Hough, Patricia D.; Hu, Kenneth T.; Jakeman, John D.; Khalil, Mohammad; Maupin, Kathryn A.; Monschke, Jason A.; Ridgway, Elliott M.; Rushdi, Ahmad A.; Seidl, Daniel T.; Stephens, John A.; Winokur, Justin G.
The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers.
