Electric Vehicles (EV) present a unique challenge to electric power system (EPS) operations because of the potential magnitude and timing of load increases due to EV charging. Time-of-Use (TOU) electricity pricing is an established way to reduce peak system loads. It is effective at shifting the timing of some customer-activated residential loads – such as dishwashers, washing machines, or HVAC systems – to off-peak periods. EV charging, though, can be larger than typical residential loads (up to 19.2 kW) and may have on-board controls that automatically begin charging according to a pre-set schedule, such as when off-peak periods begin. To understand and quantify the potential impact of EV charging's response to TOU pricing, this paper simulates 10 distribution feeders with predicted 2030 EV adoption levels. The simulation results show that distribution EPS experience an increase in peak demand as high as 20% when a majority of the charging begins immediately after on-peak times end, as might occur if EV charging is automatically scheduled. However, if charging start times are randomized within the off-peak period, EV charging is spread out and the simulations showed a decrease in the peak load to be 5% lower than results from simulations that did not implement TOU rates.
The value of long-term wave hindcasts for investigating wave climates, wave energy resources, and extreme wave conditions has motivated research developing, calibrating and validating wave hindcast models. Past hindcast model validation studies examined the accuracy in modeling bulk wave parameters of overall sea states without considering the dependency of the model's skill within different sea states. In the present study, a framework for wave hindcast model validation is developed by examining the model accuracy for the most frequently occurring sea states, sea states contributing the most energy to total wave power, sea states associated with hurricane events, and those with the largest model error. Validations using bulk wave parameters and frequency-directional spectra at these key sea states and extreme wave conditions based on univariate and bivariate-contour methods provide insights to improve model accuracy, identifying the model's strong and weak points, and pathways for improvement, e.g., modeling wave-current interactions and adjusting wind data. This study adds to a growing body of research demonstrating that a carefully calibrated and verified spectral wave hindcast model can be used to estimate key wave energy parameters over a wide range of wave energy climates, as well as their spatial, temporal, frequency, directional, and probabilistic distributions.
While the use of machine learning (ML) classifiers is widespread, their output is often not part of any follow-on decision-making process. To illustrate, consider the scenario where we have developed and trained an ML classifier to find malicious URL links. In this scenario, network administrators must decide whether to allow a computer user to visit a particular website, or to instead block access because the site is deemed malicious. It would be very beneficial if decisions such as these could be made automatically using a trained ML classifier. Unfortunately, due to a variety of reasons discussed herein, the output from these classifiers can be uncertain, rendering downstream decisions difficult. Herein, we provide a framework for: (1) quantifying and propagating uncertainty in ML classifiers; (2) formally linking ML outputs with the decision-making process; and (3) making optimal decisions for classification under uncertainty with single or multiple objectives.
Numerical Methods for Partial Differential Equations
D'Elia, Marta D.; Aulisa, Eugenio; Capodaglio, Giacomo; Chierici, Andrea
In this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.
Single particle aerosol mass spectrometry (SPAMS), an analytical technique for measuring the size and composition of individual micron-scale particles, is capable of analyzing atmospheric pollutants and bioaerosols much more efficiently and with more detail than conventional methods which require the collection of particles onto filters for analysis in the laboratory. Despite SPAMS’ demonstrated capabilities, the primary mechanisms of ionization are not fully understood, which creates challenges in optimizing and interpreting SPAMS signals. In this paper, we present a well-stirred reactor model for the reactions involved with the laser-induced vaporization and ionization of an individual particle. The SPAMS conditions modeled in this paper include a 248 nm laser which is pulsed for 8 ns to vaporize and ionize each particle in vacuum. The ionization of 1 μm, spherical Al particles was studied by approximating them with a 0-dimensional plasma chemistry model. The primary mechanism of absorption of the 248 nm photons was pressure-broadened direct photoexcitation to Al(y2D). Atoms in this highly excited state then undergo superelastic collisions with electrons, heating the electrons and populating the lower energy excited states. We found that the primary ionization mechanism is electron impact ionization of various excited state Al atoms, especially Al(y2D). Because the gas expands rapidly into vacuum, its temperature decreases rapidly. The rate of three-body recombination (e− + e− + Al+ → Al + e−) increases at low temperature, and most of the electrons and ions produced recombine within several μs of the laser pulse. The importance of the direct photoexcitation indicates that the relative peak heights of different elements in SPAMS mass spectra may be sensitive to the available photoexcitation transitions. The effects of laser intensity, particle diameter, and expansion dynamics are also discussed.
Sangoleye, Fisayo; Johnson, Jay; Chavez, Adrian R.; Tsiropoulou, Eirini E.; Marton, Nicholas L.; Hentz, Charles R.; Yannarelli, Albert
Microgrids require reliable communication systems for equipment control, power delivery optimization, and operational visibility. To maintain secure communications, Microgrid Operational Technology (OT) networks must be defensible and cyber-resilient. The communication network must be carefully architected with appropriate cyber-hardening technologies to provide security defenders the data, analytics, and response capabilities to quickly mitigate malicious and accidental cyberattacks. In this work, we outline several best practices and technologies that can support microgrid operations (e.g., intrusion detection and monitoring systems, response tools, etc.). Then we apply these recommendations to the New Jersey TRANSITGRID use case to demonstrate how they would be deployed in practice.
Book, Cameron; Hoffman, Matthew J.; Kachuck, Samuel B.; Hillebrand, Trevor R.; Price, Stephen F.; Perego, Mauro P.; Bassis, Jeremy N.
Viscoelastic rebound of the solid Earth upon the removal of ice loads has the potential to inhibit marine ice sheet instability, thereby forestalling ice-sheet retreat and global mean sea-level rise. The timescale over which the solid Earth - ice sheet system responds to changes in ice thickness and bedrock topography places a strong control on the spatiotemporal influence of this negative feedback mechanism. In this study, we assess the impact of solid-earth rheological structure on model projections of the retreat of Thwaites Glacier, West Antarctica, and the concomitant sea-level rise by coupling the dynamic ice sheet model MALI to a regional glacial isostatic adjustment (GIA) model. We test the sensitivity of model projections of ice-sheet retreat and associated sea-level rise across a range of four different solid-earth rheologies, forced by standard ISMIP6 ocean and atmospheric datasets for the RCP8.5 climate scenario. These model parameters are applied to 500-year, coupled ice-sheet - GIA simulations. For the mantle viscosity best supported by observations, the negative GIA feedback leads to a reduction in mass loss that remains above 20% after about a hundred years. Mass-loss reduction peaks at 50% around 2300, which is when a control simulation without GIA experiences its maximum rate of retreat. For a weaker solid-earth rheology that is unlikely but compatible with observational uncertainty, mass loss reduction remains above 50% after 2150. At 2100, mass loss reduction is 10% for the best-fit rheology and 25% for the weakest rheology. At the same time, we estimate that water expulsion from the rebounding solid Earth beneath the ocean near Thwaites Glacier may increase sea-level rise by up to 20% at five centuries. Additionally, the reduction in ice-sheet retreat caused by GIA is substantially reduced under stronger climate forcings, suggesting that the stabilizing feedback of GIA will also be an indirect function of emissions scenario. We hypothesize that feedbacks between the solid Earth - ice sheet system are controlled by a competition between the spatial extent and timescale of bedrock uplift relative to the rate of grounded ice retreat away from the region of most rapid unloading. Although uncertainty in solid-earth rheology leads to large uncertainty in future sea-level rise contribution from Thwaites Glacier, under all plausible parameters the GIA effects are too large to be ignored for future projections of Thwaites Glacier of more than a century.
A mathematical framework is provided for a substructuring-based domain decomposition (DD) approach for nonlocal problems that features interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation (PDE) problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal DD methods.
The Sandia National Laboratories site sustainability plan and its associated DOE Sustainability Dashboard data entries encompass Sandia National Laboratories contributions toward meeting the DOE sustainability goals. This site sustainability plan fulfills the contractual requirement for National Technology & Engineering Solutions of Sandia, LLC, the management and operating contractor for Sandia National Laboratories, to deliver an annual sustainability plan to the DOE National Nuclear Security Administration Sandia Field Office.
In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model's input parameters. We develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model's response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.
This report examines the localization of high frequency electromagnetic fields in general three-dimensional cavities along periodic paths between opposing sides of the cavity. The focus is on the case where the mirrors at the ends of the orbit are concave and have two different radii of curvature. The cases where these orbits lead to unstable localized modes are known as scars. The ellipsoidal coordinate system is utilized in the construction of the scarred modes. The field at the interior foci is examined as well as trigonometric projections along the periodic scarred ray path.
Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. In this work we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. We also introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.
One detonation test was monitored for impulse noise at Thunder Range on November 6, 2020. The TNT equivalency for this shot was 100 lbs. Ground zero for this test was located at an open area of Range 7. The Armag, where MOW were remotely located, was just south of the east end of the TR shock tube The purpose of this sampling event was to characterize the noise attenuation provided by the Armag which will be used as the primary firing location in a future test series. To determine attenuation provided by the Armag, one sound level monitor was placed inside and a second monitor was placed outside the hardened structure at the same distance from ground zero. The Armag was located 1300 feet from ground zero. Essential personnel performing these tests were remotely located inside the Armag and wore hearing protection with a minimum NRR of 23. Members of the Workforce (MOW) who are exposed to noise levels above 140 dBC, regardless of hearing protection worn, are required to be enrolled into the SNL Hearing Conservation Program which includes audiometric testing, online training (HCP100) and wearing hearing protection.
Niobium doped lead-tin-zirconate-titanate ceramics near the PZT 95/5 orthorhombic AFE – rhombohedral FE morphotropic phase boundary Pb1-0.5y(Zr0.865-xTixSn0.135)1-yNbyO3 were prepared according to a 22+1 factorial design with x = 0.05, 0.07 and y = 0.0155, 0.0195. The ceramics were prepared by a traditional solid-state synthesis route and sintered to near full density at 1250°C for 6 h. All compositions were ∼98% dense with no detectable secondary phases by X-ray diffraction. The ceramics exhibited equiaxed grains with intergranular porosity, and grain size was ∼5 µm, decreasing with niobium substitution. Compositions exhibited remnant polarization values of ∼32 µC/cm2, increasing with Ti substitution. Depolarization by the hydrostatic pressure induced FE-AFE phase transition was drastically affected by variation of the Ti and Nb substitution, increasing at a rate of 113 MPa /1% Ti and 21 MPa/1% Nb. Total depolarization output was insensitive to the change in Ti and Nb substitution, ∼32.8 µC/cm2 for the PSZT ceramics. The R3c-R3m and R3m-Pm3m phase transition temperatures on heating ranged from 90 to 105°C and 183 to 191°C, respectively. Ti substitution stabilized the R3c and R3m phases to higher temperatures, while Nb substitution stabilized the Pm3m phase to lower temperatures. Thermal hysteresis of the phase transitions was also observed in the ceramics, with transition temperature on cooling being as much as 10°C lower.
This manual describes the use of the Xyce™ Parallel Electronic Simulator. Xyce™ has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. (2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. (3) Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using modern coding practices. Xyce™ is a parallel code in the most general sense of the phrase—a message passing parallel implementation—which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel eficiency is achieved as the number of processors grows.
Cyber security has been difficult to quantify from the perspective of defenders. The effort to develop a cyber-attack with some ability, function, or consequence has not been rigorously investigated in Operational Technologies. This specification defines a testing structure that allows conformal and repeatable cyber testing on equipment. The purpose of the ETE is to provide data necessary to analyze and reconstruct cyber-attack timelines, effects, and observables for training and development of Cyber Security Operation Centers. Standardizing the manner in which cyber security on equipment is investigated will allow a greater understanding of the progression of cyber attacks and potential mitigation and detection strategies in a scientifically rigorous fashion.
Members of the Nuclear Criticality Safety (NCS) Program at Sandia National Laboratories (SNL) have updated the suite of benchmark problems developed to validate MCNP6 Version 2.0 for use in NCS applications. The updated NCS benchmark suite adds approximately 600 new benchmarks and includes peer review of all input files by two different NCS engineers (or one NCS engineer and one candidate NCS engineer). As with the originally released benchmark suite, the updated suite covers a broad range of fissile material types, material forms, moderators, reflectors, and neutron energy spectra. The benchmark suite provides a basis to establish a bias and bias uncertainty for use in NCS analyses at SNL.
A mathematical framework is provided for a substructuring-based domain decomposition (DD) approach for nonlocal problems that features interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation (PDE) problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal DD methods.
Elastomeric rubbers serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose-liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. Cyclic (de)pressurization is known to degrade these materials through the process of cavitation. This readily visible failure mode occurs as a fracture or rupture of the material and is due to the oversaturated gas localizing to form gas bubbles. Computational modeling in the Hydrogen Materials Compatibility Program (H-Mat), co-led by Sandia National Laboratories and Pacific Northwest National Laboratory, employs multi-scale simulation efforts to build a predictive understanding of hydrogen-induced damage in materials. Modeling efforts within the project aim to provide insight into how to formulate materials that are less sensitive to high-pressure hydrogen-induced failure. In this document, we summarize results from atomistic molecular dynamics simulations, which make predictive assessments of the effects of compositional variations in the commonly used elastomer, ethylene propylene diene monomer (EPDM).
The recent discovery of bright, room-temperature, single photon emitters in GaN leads to an appealing alternative to diamond best single photon emitters given the widespread use and technological maturity of III-nitrides for optoelectronics (e.g. blue LEDs, lasers) and high-speed, high-power electronics. This discovery opens the door to on-chip and on-demand single photon sources integrated with detectors and electronics. Currently, little is known about the underlying defect structure nor is there a sense of how such an emitter might be controllably created. A detailed understanding of the origin of the SPEs in GaN and a path to deterministically introduce them is required. In this project, we develop new experimental capabilities to then investigate single photon emission from GaN nanowires and both GAN and AlN wafers. We ion implant our wafers with the ion implanted with our focused ion beam nanoimplantation capabilities at Sandia, to go beyond typical broad beam implantation and create single photon emitting defects with nanometer precision. We've created light emitting sources using Li+ and He+, but single photon emission has yet to be demonstrated. In parallel, we calculate the energy levels of defects and transition metal substitutions in GaN to gain a better understanding of the sources of single photon emission in GaN and AlN. The combined experimental and theoretical capabilities developed throughout this project will enable further investigation into the origins of single photon emission from defects in GaN, AlN, and other wide bandgap semiconductors.
The wide variety of inverter control settings for solar photovoltaics (PV) causes the accurate knowledge of these settings to be difficult to obtain in practice. This paper addresses the problem of determining inverter reactive power control settings from net load advanced metering infrastructure (AMI) data. The estimation is first cast as fitting parameterized control curves. We argue for an intuitive and practical approach to preprocess the AMI data, which exposes the setting to be extracted. We then develop a more general approach with a data-driven reactive power disaggregation algorithm, reframing the problem as a maximum likelihood estimation for the native load reactive power. These methods form the first approach for reconstructing reactive power control settings of solar PV inverters from net load data. The constrained curve fitting algorithm is tested on 701 loads with behind-the-meter (BTM) PV systems with identical control settings. The settings are accurately reconstructed with mean absolute percentage errors between 0.425% and 2.870%. The disaggregation-based approach is then tested on 451 loads with variable BTM PV control settings. Different configurations of this algorithm reconstruct the PV inverter reactive power timeseries with root mean squared errors between 0.173 and 0.198 kVAR.
Geologic Disposal Safety Assessment Framework is a state-of-the-art simulation software toolkit for probabilistic post-closure performance assessment of systems for deep geologic disposal of nuclear waste developed by the United States Department of Energy. This paper presents a generic reference case and shows how it is being used to develop and demonstrate performance assessment methods within the Geologic Disposal Safety Assessment Framework that mitigate some of the challenges posed by high uncertainty and limited computational resources. Variance-based global sensitivity analysis is applied to assess the effects of spatial heterogeneity using graph-based summary measures for scalar and time-varying quantities of interest. Behavior of the system with respect to spatial heterogeneity is further investigated using ratios of water fluxes. This analysis shows that spatial heterogeneity is a dominant uncertainty in predictions of repository performance which can be identified in global sensitivity analysis using proxy variables derived from graph descriptions of discrete fracture networks. New quantities of interest defined using water fluxes proved useful for better understanding overall system behavior.
The use of simple models for response prediction of building structures is preferred in earthquake engineering for risk evaluations at regional scales, as they make computational studies more feasible. The primary impediment in their gainful use presently is the lack of viable methods for quantifying (and reducing upon) the modeling errors/uncertainties they bear. This study presents a Bayesian calibration method wherein the modeling error is embedded into the parameters of the model. The method is specifically described for coupled shear-flexural beam models here, but it can be applied to any parametric surrogate model. The major benefit the method offers is the ability to consider the modeling uncertainty in the forward prediction of any degree-of-freedom or composite response regardless of the data used in calibration. The method is extensively verified using two synthetic examples. In the first example, the beam model is calibrated to represent a similar beam model but with enforced modeling errors. In the second example, the beam model is used to represent the detailed finite element model of a 52-story building. Both examples show the capability of the proposed solution to provide realistic uncertainty estimation around the mean prediction.
Photovoltaic (PV) performance is affected by reversible and irreversible losses. These can typically be mitigated through responsive and proactive operations and maintenance (O&M) activities. However, to generate profit, the cost of O&M must be lower than the value of the recovered electricity. This value depends both on the amount of recovered energy and on the electricity prices, which can vary significantly over time in spot markets. The present work investigates the impact of the electricity price variability on the PV profitability and on the related O&M activities in Italy, Portugal, and Spain. It is found that the PV revenues varied by 1.6 × to 1.8 × within the investigated countries in the last 5 years. Moreover, forecasts predict higher average prices in the current decade compared to the previous one. These will increase the future PV revenues by up to 60% by 2030 compared to their 2015–2020 mean values. These higher revenues will make more funds available for better maintenance and for higher quality components, potentially leading to even higher energy yield and profits. Linearly growing or constant price assumptions cannot fully reproduce these expected price trends. Furthermore, significant price fluctuations can lead to unexpected scenarios and alter the predictions.
On October 1, 2022, sound level measurements were taken at various locations throughout Kirtland Air Force Base (KAFB) and Southeastern Albuquerque. The purpose was to support sound propagation modeling predictions and sound regulations for public exposure during the detonation of an approximately 300-pound energetic experiment. Ground Zero was located on Range 7 of Sandia Thunder Range (06647). A total of 8 measurement locations were identified (e.g., 5 on KAFB and 3 in the Southeastern Albuquerque neighborhoods).
Members of the Workforce (MOW) who are exposed to noise levels above 140 dBC, regardless of hearing protection worn, are required to be enrolled into the SNL Hearing Conservation Program which includes audiometric testing, online training (HCP100) and wearing hearing protection. Based on the area impact noise sample results, the attenuation provided by the MFCP was protective for mitigating noise to levels below the ACGIH TLV of 140 dBC. The results also validated the scaled distance equation in an open-air environment as the results at K635 (864 feet) were below 140 dBC.
The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, natural gas, and autogas systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying leak sizes and rates, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 5.0. The most significant change for HyRAM+ version 5.0 from HyRAM+ version 4.1 is the ability to model blends of different fuels. HyRAM+ was previously only suitable for use with hydrogen, methane, or propane, with users having the ability to use methane as a proxy for natural gas and propane as a proxy for autogas/liquefied petroleum gas. In version 5.0, real natural gas or autogas compositions can be modeled as the fuel, or even blends of natural gas with hydrogen. These blends can be used in the standalone physics models, but not yet in the quantitative risk assessment mode of HyRAM+.
Self-determination has been an on-going effort for Native American people and gained much traction with the passing of The Energy Policy Act of 2005, which included the Indian Tribal Energy Development and Self-Determination Act. Congress passed this act to assist Native American tribes and Alaska Native villages with planning, development, and assistance to achieve their energy goals. The Ute Mountain Ute Tribe (UMUT) has relied on oil and natural gas for economic support the last 70 years. Burning fossil fuels, along with oil and gas development, decreases the quality of air and leads to increased greenhouse gas emissions. Subsequently, the burning of fossil fuels to produce energy is now more costly than many renewable energy sources, including solar photovoltaic (PV) systems. Environmental stewardship, along with the need to maintain revenue generation, has led UMUT’s efforts to achieve energy self-determinism employing PV and exploring other technology. In the past, the tribe completed a 1 megawatt PV project near Towaoc, Colorado, which serves as a case study on the tribe’s energy goals: a future where renewables will dominate their energy landscape. This paper explores UMUT’s past and on-going efforts toward energy independence and how it relates to the broader landscape of Native American energy sovereignty.
We present a proof-of-concept demonstration of a narrow linewidth $^{87}$Rb magneto-optical trap (MOT) operating on the narrow linewidth $5S_{1/2}$ → $6P_{3/2}$ transition at 420 nm. We stabilized the absolute frequency of the 420 nm laser to an atomic transition in $^{87}$Rb and demonstrate a MOT using 420 nm light driving the $5S_{1/2}$, $F = 2$ → $6P_{3/2}, F' = 3$ transition. We then use tome-of-flight measurements to characterize the 420 nm MOT temperature, observing a minimum temperature of about $T^{(420)}_{horizontal}$ = 150μK and $T^{(420)}_{vertical}$ = 250μK before the opportunity to perform significant characterization and optimization. Although this temperature is significantly higher then the expected 420 nm Doppler cooling limit ($T_D^{(420)}$ ≈ 34 μK), these are already approaching the Doppler limit of a standard 780 nm MOT ($T_D^{(780)}$ ≈ 146 μK). We believe that with further optimization the Doppler cooling limit of ≈ 34 μK can be achieved. This initial result answers our key research question and demonstrates the viability of applying narrow linewidth laser cooling as a robust technique for future fieldable quantum sensors.
This document is a reference guide to the Xyce™ Parallel Electronic Simulator, and is a companion document to the Xyce™ Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce™. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce™ Users' Guide.
Mishra, Umakant; Yeo, Kyongmin; Adhikari, Kabindra; Riley, William J.; Hoffman, Forrest M.; Hudson, Corey
Accurate representation of environmental controllers of soil organic carbon (SOC) stocks in Earth System Model (ESM) land models could reduce uncertainties in future carbon–climate feedback projections. Using empirical relationships between environmental factors and SOC stocks to evaluate land models can help modelers understand prediction biases beyond what can be achieved with the observed SOC stocks alone. In this study, we used 31 observed environmental factors, field SOC observations (n = 6,213) from the continental United States, and two machine learning approaches (random forest [RF] and generalized additive modeling [GAM]) to (a) select important environmental predictors of SOC stocks, (b) derive empirical relationships between environmental factors and SOC stocks, and (c) use the derived relationships to predict SOC stocks and compare the prediction accuracy of simpler model developed with the machine learning predictions. Out of the 31 environmental factors we investigated, 12 were identified as important predictors of SOC stocks by the RF approach. In contrast, the GAM approach identified six (of those 12) environmental factors as important controllers of SOC stocks: potential evapotranspiration, normalized difference vegetation index, soil drainage condition, precipitation, elevation, and net primary productivity. The GAM approach showed minimal SOC predictive importance of the remaining six environmental factors identified by the RF approach. Our derived empirical relations produced comparable prediction accuracy to the GAM and RF approach using only a subset of environmental factors. The empirical relationships we derived using the GAM approach can serve as important benchmarks to evaluate environmental control representations of SOC stocks in ESMs, which could reduce uncertainty in predicting future carbon–climate feedbacks.
Barium titanate (BTO) is a ferroelectric material used in capacitors because of its high bulk dielectric constant. However, the impact of the size of BTO on its dielectric constant is not yet fully understood and is highly contested. Here, we present an investigation into the dielectric constant of BTO nanoparticles with diameters ranging between 50 and 500 nm. BTO nanoparticles were incorporated into acrylonitrile butadiene styrene and injection molded into parallel plate capacitors, which were used to determine nanocomposite dielectric constants. The dielectric constants of BTO nanoparticles were obtained by combining experimental measurements with computational results from COMSOL simulations of ABS-matrix nanocomposites containing BTO. The dielectric constant of BTO was observed to be relatively constant at nanoparticle diameters as small as 200 nm but sharply declined at smaller nanoparticle sizes. These results will be useful in the development of improved energy storage and power conditioning systems utilizing BTO nanoparticles.
Battery storage systems are increasingly being installed at photovoltaic (PV) sites to address supply-demand balancing needs. Although there is some understanding of costs associated with PV operations and maintenance (O&M), costs associated with emerging technologies such as PV plus storage lack details about the specific systems and/or activities that contribute to the cost values. This study aims to address this gap by exploring the specific factors and drivers contributing to utility-scale PV plus storage systems (UPVS) O&M activities costs, including how technology selection, data collection, and related and ongoing challenges. Specifically, we used semi-structured interviews and questionnaires to collect information and insights from utility-scale owners and operators. Data was collected from 14 semi-structured interviews and questionnaires representing 51.1 MW with 64.1 MWh of installed battery storage capacity within the United States (U.S.). Differences in degradation rate, expected life cycle, and capital costs are observed across different storage technologies. Most O&M activities at UPVS related to correcting under-performance. Fires and venting issues are leading safety concerns, and owner operators have installed additional systems to mitigate these issues. There are ongoing O&M challenges due the lack of storage-specific performance metrics as well as poor vendor reliability and parts availability. Insights from this work will improve our understanding of O&M consideration at PV plus storage sites.