Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluating Compliance Inspection Data to Understand Emissions from Marginal Conventional Oil & Gas Wells in New York State

It is estimated that over half (~66%) of active oil and gas wells are marginally productive, producing less than 15 barrels of oil equivalent per day. This presentation showcases an evaluation of New York State compliance inspection reports and discusses the current capabilities of machine learning to identify high-risk MCWs that could be targeted during methane emission reduction efforts and future model deployments. Findings are valuable to the New York Department of Environmental Conservation and other state oil and gas regulatory agencies seeking to prioritize the spending of Methane Emission Reduction Program funds to plug marginal conventional wells with emissions.

compliance inspection data↗
DOE OSTI2025

Hourly Electricity Demand Profiles for Each County in the Contiguous United States

This dataset provides estimated hourly electricity demand for each county in the contiguous United States from 2016-2023. The demand profiles represent the sum of two components: (1) Weighted averages of reported hourly demand profiles for North American Electric Reliability Corporation balancing authority (BA) regions and subregions, scaled to match annual estimates of county-level retail sales and direct use of electricity and weighted by the estimated percentage of county load served by each BA region or subregion. (2) Weighted averages of modeled hourly, county- and sector-level distributed photovoltaic (DPV) capacity factor profiles, scaled to match annual estimates of on-site consumption of DPV-generated electricity for each county and weighted by the percentage of consumption attributable to each sector Annual county-level retail sales are estimated by aggregating utility-reported sales to the state level and allocating the results to counties according to each county's share of state population. Annual county-level direct use is calculated by aggregating power plant-reported direct use values. Annual county-level on-site consumption of DPV-generated electricity is estimated by aggregating utility-reported net metering data to determine the amount of DPV-generated electricity sold back to the grid for each state, subtracting those values from modeled state-level DPV generation estimates, and allocating the results to counties according to each county's share of statewide modeled DPV generation. The open-source Python code used to develop this dataset is available at "Historical Load Data Repository" link below.

14 SOLAR ENERGY↗
DOE OSTI2020

A Privacy Preserving Model-Free Optimization and Control Framework for Demand Response from Residential Thermal Loads

We consider the problem of optimizing the cost of procuring electricity for a large collection of homes managed by a load serving entity, by pre-cooling or pre-heating the thermal inertial loads in the homes to avoid procuring power during periods of peak electricity pricing. We would like to accomplish this objective in a completely privacy-preserving and model-free manner, that is, without direct access to the state variables (temperatures or power consumption) or the dynamical models (thermal characteristics) of individual homes, while guaranteeing personal comfort constraints of the consumers. We propose a two-stage optimization and control framework to address this problem. In the first stage, we use a long short-term memory (LSTM) network to predict hourly electricity prices, based on historical pricing data and weather forecasts. Given the hourly price forecast and thermal models of the homes, the problem of designing an optimal power consumption trajectory that minimizes the total electricity procurement cost for the collection of thermal loads can be formulated as a large-scale integer program (with millions of variables) due to the on-off cyclical dynamics of such loads. We provide a simple heuristic relaxation to make this large-scale optimization problem model-free and computationally tractable. In the second stage, we translate the results of this optimization problem into distributed open-loop control laws that can be implemented at individual homes without measuring or estimating their state variables, while simultaneously ensuring consumer comfort constraints. We demonstrate the performance of this approach on a large-scale test case comprising of 500 homes in the Houston area and benchmark its performance against a direct model-based optimization and control solution.

Sivaranjani, S.↗
DOE OSTI2020

Optimization-Based Fast-Frequency Estimation and Control of Low-Inertia Microgrids

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.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Variational Quantum Fidelity Estimation

Computing quantum state fidelity will be important to verify and characterize states prepared on a quantum computer. In this work, we propose novel lower and upper bounds for the fidelity F ( ρ , σ ) based on the ``truncated fidelity'' F ( ρ m , σ ) , which is evaluated for a state ρ m obtained by projecting ρ onto its m -largest eigenvalues. Our bounds can be refined, i.e., they tighten monotonically with m . To compute our bounds, we introduce a hybrid quantum-classical algorithm, called Variational Quantum Fidelity Estimation, that involves three steps: (1) variationally diagonalize ρ , (2) compute matrix elements of σ in the eigenbasis of ρ , and (3) combine these matrix elements to compute our bounds. Our algorithm is aimed at the case where σ is arbitrary and ρ is low rank, which we call low-rank fidelity estimation, and we prove that no classical algorithm can efficiently solve this problem under reasonable assumptions. Finally, we demonstrate that our bounds can detect quantum phase transitions and are often tighter than previously known computable bounds for realistic situations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Memristive Behavior Enabled by Amorphous–Crystalline 2D Oxide Heterostructure

The emergence of memristive behavior in amorphous–crystalline 2D oxide heterostructures, which are synthesized by atomic layer deposition (ALD) of a few-nanometer amorphous Al 2 O 3 layers onto atomically thin single-crystalline ZnO nanosheets, is demonstrated. The conduction mechanism is identified based on classic oxygen vacancy conductive channels. ZnO nanosheets provide a 2D host for oxygen vacancies, while the amorphous Al 2 O 3 facilitates the generation and stabilization of the oxygen vacancies. The conduction mechanism in the high-resistance state follows Poole–Frenkel emission, and in the the low-resistance state is fitted by the Mott–Gurney law. From the slope of the fitting curve, the mobility in the low-resistance state is estimated to be ≈2400 cm 2 V -1 s -1 , which is the highest value reported in semiconductor oxides. When annealed at high temperature to eliminate oxygen vacancies, Al is doped into the ZnO nanosheet, and the memristive behavior disappears, further confirming the oxygen vacancies as being responsible for the memristive behavior. The 2D heterointerface offers opportunities for new design of high-performance memristor devices.

2D heterostructures↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2022

Kernel learning backward SDE filter for data assimilation

In this paper, we develop a kernel learning backward SDE filter method to estimate the state of a stochastic dynamical system based on its partial noisy observations. A system of forward backward stochastic differential equations is used to propagate the state of the target dynamical model, and Bayesian inference is applied to incorporate the observational information. Further, to characterize the dynamical model in the entire state space, we introduce a kernel learning method to learn a continuous global approximation for the conditional probability density function of the target state by using discrete approximated density values as training data. Numerical experiments demonstrate that the kernel learning backward SDE is highly effective.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Plasma-parameter dependence of ro-vibrational temperatures for H 2 in LHD divertor

We analyzed a thousand visible spectra of Fulcher-α band measured for divertor plasmas in Large Helical Device. With a coronal model and Bayesian inference, the population distribution of hydrogen molecule in the electronic ground state was estimated. The non-thermal population distribution was recovered with a two-temperature model, which has two sets of rotational and vibrational temperatures, as well as their mixture coefficient. The lower rotational temperature significantly changes according to the plasma parameters. Its nearly linear dependence on the electron density was found, which is consistent with previous works. The lower vibrational temperature also shows a small density dependence, as reported by a previous work. On the other hand, the higher rotational and vibrational temperatures as well as the mixture coefficient only show slight changes over the broad range of plasma parameters. These population parameters show a significant correlation; with higher electron density, all the temperatures and the fraction of the higher-temperature component increase simultaneously. In conclusion, this suggests that the electron impact plays an important role to determine the population distribution.

42 ENGINEERING↗
DOE OSTI2023

California's harvested wood products: A time-dependent assessment of life cycle greenhouse gas emissions

Following life-cycle assessment (LCA) methodology, this study presents a state-level estimation of embodied carbon of wood products harvested in 2019 from California and subsequently processed, manufactured, transported, used, and disposed at the end-of-life (EoL). In a conventional static approach to LCA, all GHG emissions were aggregated and considered to occur at year 0 of the given time horizon (500 years in this study) and used a static characterization factor (CF). In dynamic LCA, GHG emissions occurring in different years were considered, and their global warming impact (GWI) was determined using a time-dependent CF over the selected time horizon of 500 years. Four scenarios were developed to examine the impact of EoL choices on GWI. It was found that dynamic GWI for all scenarios ranged from 0.27 to 0.93 million tonne CO₂e, which were 45–73 % lower than those estimated with static LCA approach, indicating that the static LCA approach could lead to an underestimation of the benefits of substituting wood for non-wood products, compared to those based on dynamic LCA approach. This analysis also demonstrated that the choice of EoL treatment option is a key factor affecting the estimated GWI as it directly determines the annual emission of GHGs released into atmosphere and subsequently their warming effect depending on the time harvested wood products (HWPs) spend in the horizon of assessment. Altogether, the dynamic LCA performed in this study enabled more robust interpretations of embodied carbon by including temporal boundaries associated with the HWPs life cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Physical Properties of Candidate X-ray Detector Material Rb 4 Ag 2 BiBr 9

Recently, metal halide perovskites have emerged as promising semiconductor candidates for sensitive X-ray photon detection due to their suitable band gap energies, excellent charge transport properties, and low material cost afforded by their low-temperature solution-processing preparation. Furthermore, we report an improved methodology for single crystal growth and thermal and electrical properties of a two-dimensional (2D) layered halide material Rb 4 Ag 2 BiBr 9 , which has been identified as a potential candidate for X-ray radiation detection applications. The measured heat capacity for Rb 4 Ag 2 BiBr 9 implies that there are no structural phase transitions upon cooling. Temperature dependence of thermal transport measurements further suggests remarkably low thermal conductivities of Rb 4 Ag 2 BiBr 9 that are comparable to the lowest reported in literature. The bulk crystal resistivity is determined to be 2.59 × 10 9 Ω·cm from the current–voltage (I–V) curve. Density of trap states is estimated to be ~10 10 cm –3 using the space-charge-limited-current measurements. The fabricated Rb 4 Ag 2 BiBr 9 -based X-ray detector shows good operational stability with no apparent current drift, which may be ascribed to the 2D crystal structure of Rb 4 Ag 2 BiBr 9 . Finally, by varying the X-ray tube current to change the corresponding dose rate, the Rb4Ag2BiBr9 X-ray detector sensitivity is determined to be 222.03 μC Gy –1 cm –2 (at an electric field of E = 24 V/mm).

36 MATERIALS SCIENCE↗
DOE OSTI2025

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermoelastic Properties of Fe 3+ ‐Rich Jeffbenite and Application to Superdeep Diamond Barometry

Abstract Jeffbenite (Mg 3 Al 2 Si 3 O 12 ) is a tetragonal phase found in so far only in superdeep diamonds, and its thermoelastic parameters are a prerequisite for determining entrapment pressures as it is regarded as a potential indicator for superdeep diamonds. In this study, the thermoelastic properties of synthetic Fe 3+ ‐jeffbenite were measured up to 33.7 GPa and 750 K. High‐temperature static compression data were fitted, giving ( ∂K T 0 / ∂ T ) P = −0.0107 (4) GPa/K and α T = 3.50 (3) × 10 −5 K −1 . The thermoelastic properties and phase stability are applied to modeling isomekes, or P‐T paths intersecting possible conditions of entrapment in diamond. We calculate that under ideal exhumation, jeffbenite entrapped at mantle transition zone conditions will exhibit a high remnant pressure at 300 K ( P inc ) of ∼5.0 GPa. Elastic geobarometry on future finds of jeffbenite inclusions can use the new equation of state to estimate entrapment pressures for this phase with still highly uncertain stability field in the mantle.

Geology↗
DOE OSTI2021

Ultraviolet spectroscopy of pressurized and supercritical carbon dioxide

Abstract Carbon dioxide (CO 2 ) is prevalent in planetary atmospheres and sees use in a variety of industrial applications. Despite its ubiquitous nature, its photochemistry remains poorly understood. In this work we explore the density dependence of pressurized and supercritical CO 2 electronic absorption spectra by vacuum ultraviolet spectroscopy over the wavelength range 1455-2000 Å. We show that the lowest absorption band transition energy is unaffected by a density increase up to and beyond the thermodynamic critical point (137 bar, 308 K). However, the diffuse vibrational structure inherent to the spectrum gradually decreases in magnitude. This effect cannot be explained solely by collisional broadening and/or dimerization. We suggest that at high densities close proximity of neighboring CO 2 molecules with a variety of orientations perturbs the multiple monomer electronic state potential energy surfaces, facilitating coupling between binding and dissociative states. We estimate a critical radius of ~4.1 Å necessary to cause such perturbations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗
DOE OSTI2023

Shock-induced transformation of nitinol shape memory alloy: Effect of stress state on transformation

Due to its numerous practical applications and intriguing phase transformation behavior, shape memory alloys (SMAs) have garnered significant research and development interests. In the past, most studies on the mechanical behavior of SMAs have been conducted under uniaxial stress loadings. Limited research on SMAs under shock loading has not provided conclusive results regarding their transformation behavior and transformation stress under such loading. Additionally, there is a lack of comprehensive understanding regarding the effects of different stress states on transformation behavior. The main objectives of this study are to address these issues. To achieve these objectives, a series of shock wave experiments were designed and conducted. Additionally, quasi-static and dynamic uniaxial stress experiments were carried out to establish a baseline for comparison. The results revealed that the transformation stress under dynamic uniaxial strain shock loading was approximately 1.92 GPa in contrast to 0.5 GPa (quasi-static) to 0.8 GPa (dynamic) observed in uniaxial stress loading. The transformation behavior exhibited noticeable rate sensitivity for both types of loading. There appeared to be a critical strain rate above which the austenite phase was driven to a metastable state. This estimated critical axial strain rate along the loading direction was approximately 2 × 10 3 /s–4 × 10 3 /s for uniaxial stress loading and approximately 2 × 10 6 /s for uniaxial strain loading. The apparent high transformation stress for uniaxial strain loading can likely be attributed to a combination of high-pressure confinement and high strain rate. Furthermore, determining their relative contributions remains an open issue.

36 MATERIALS SCIENCE↗