Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 487 records · Page 27

Distributed Sensing and Computer Vision Methods for Advanced Air Mobility Approach and Landing

Advanced Air Mobility (AAM) aircraft require precision approach and landing systems (PALS) in several types of environments such as urban, suburban, and rural. It is difficult to implement current state-of-the-art methods approved for automated approach and landing for AAM operations. However, existing technology and systems that use vision, IR, radar, and GPS methods provide baseline perception and sensing requirements for AAM aircraft approach and landing. This paper focuses on vision-based PAL and computer vision feature correspondence methods to demonstrate a baseline navigation system while adhering to the Federal Aviation Administration requirements and regulations. The coplanar algorithm determines pose estimation, which feeds into an Extended Kalman filter that combines IMU with vision to create a sensor fusion navigation solution for GPS-denied environments. The state estimate leads to glideslope and localizer error computations, which will be pertinent for designing and deriving guidance laws and control laws for AAM PALS. The IMU and vision navigation solution provides promising simulation results for AAM PALS. This paper builds on previous work by incorporating high fidelity simulations with computer graphics rendering to demonstrate a distributed sensor network to track an AAM aircraft during approach and landing to compare with the aircraft's onboard navigation solution.

Evan Kawamura↗

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Online Monitoring Applications Enabled by Phasor Measurement Units: Technical Assistance to the Power Sectors of Southeast Asia

In this report, an overview of several online applications enabled by PMU measurements is provided. Besides a brief technical background for each application, the report also discusses control room displays and alarming methodologies used by North American organizations, and applicable standards set by the North American Electric Reliability Corporation (NERC). The five applications discussed herein are inertia monitoring, linear state estimation, voltage stability monitoring, small-signal stability monitoring, and forced oscillation monitoring.

24 POWER TRANSMISSION AND DISTRIBUTION↗

CalderaCast User Manual Version 1.0

CalderaCast is a user-friendly web-based tool for electrical-load forecast, providing stakeholders with a fully customizable decision-support framework that estimates the likely power draw from a possible future electric-vehicle (EV) charging station at a given location on a given day along an alternative fuel corridor (AFC). These EV charging profiles are accurately modeled in CalderaCast using the Caldera software framework developed by Idaho National Laboratories (INL), reflecting the realistic charging levels observed in actual charge events. This tool was developed as part of the National Electric Vehicle Infrastructure (NEVI) program, which is quickly generating substantial interest from would-be charging station operators (CSO), large and small electric utilities, and state transportation planners, some of whom had not seriously considered EV charging previously. All these entities—with or without background in EV infrastructure—must estimate the electricity load that a proposed charging station will generate. This load forecast is critically important for a utility to properly assess the capacity of their distribution network to support the proposed station or properly size grid upgrades for potential load growth due to future EV adoption, vehicle technology improvements, or station growth. This document describes each aspect of the CalderaCast tool and provides guidance to users who are interested in utilizing the tool for their work.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Tensor decompositions for count data that leverage stochastic and deterministic optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the maximum likelihood estimator (MLE) of the Poisson CPD model. Here, this work presents two new algorithms that extend state-of-the-art local methods for Poisson CPD. Hybrid GCP-CPAPR combines Generalized Canonical Decomposition (GCP) with stochastic optimization and CP Alternating Poisson Regression (CPAPR), a deterministic algorithm, to increase the probability of converging to the MLE over either method used alone. Restarted CPAPR with SVDrop uses a heuristic based on the singular values of the CPD model unfoldings to identify convergence toward optimizers that are not the MLE and restarts within the feasible domain of the optimization problem, thus reducing overall computational cost when using a multi-start strategy. We provide empirical evidence that indicates our approaches outperform existing methods with respect to converging to the Poisson CPD MLE.

CPAPR↗

Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) (Final Scientific/Technical Report)

This report summarizes the activities conducted under the DOE VTO funded project DE- EE0008451, where ABB Inc. (ABB), in collaboration with Idaho National Laboratory (INL), APS Global (APS), and XOS Trucks (XOS) pursued the development of a cyber-resilient extreme fast charging (XFC) management system. This project entitled Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) focuses on a resilient architecture for smart charging EV Supply Equipment (EVSE) device control and Coordinated Anomaly Detection System (CADS) features that can be added at the charging site depot level to increase cybersecurity. The project was split into two budget periods focused first on developing the threat model and resilient control concepts and second on testing, improving, and validating those developed resilient control algorithms and features with a focus on key vulnerabilities identified during the threat assessment portion of the project. During the first budget period of the CyberX project, the ABB led team focused on activities to identify, model, and quantitatively prioritize high-impact attack scenarios with potential cyber-physical effects while also modeling and developing concepts for a resilient control system that could securely address integration of DERs and other resources with EV charging. Development of the security focused XFC management system (XMS) was accomplished first by offline simulation using a developed XFC station or depot with 480V input level and simulating measurement inputs to monitoring and control systems in concept development. A representative distribution grid model was developed supporting an EV charging site model with BESS and 6 general EV charging models. These EV charging models allowed multiple configurations of charging level, multiple connected protection and measurement devices, and simulation function to show general compromise of EV, BESS, and protection features based on parallel threat analysis. During the second budget period, the EV site and supporting systems model was developed in more detail and converted from offline model to real-time to real-time with EV charging hardware in the loop (HIL). The resilient control architecture developed as concept in the first part of the project was further tested and validated for integration of local energy resources and XFC charging station site equipment while maintaining cybersecure operating principles. The proposed resilient architecture for smart charging and cybersecurity features consists of two main concepts developed and tested within the project. The first concept is an XFC management system (XMS) consisting of a hardware gateway, software platform, and Supervisory Control and Data Acquisition (SCADA) or Distribution Management System integration components. The second concept is a Coordinated Anomaly Detection System (CADS) which forms a primarily software-related subsystem of the total CyberX solution focused on monitoring system measurements, estimation of measurement states, and predicting current at the utility point of interaction based on machine learning for anomaly detection.

33 ADVANCED PROPULSION SYSTEMS↗

ARPA-E PNDIODES "PN Junctions by Ion Implantation" Project (Final Report)

Electricity accounts for nearly 40% of the energy generated in the United States, of which 65% is lost to electrical system inefficiency before reaching end-use sectors. Power electronics will play a growing role in this distribution and end-use, with up to 80% of electricity estimated to pass through power electronics from generation to consumption by 2030. This gives technical improvements in power electronics enormous potential to drive efficiency gains throughout the U.S. economy. Current Si-based power architectures including MOSFETs and IGBTs are limited by high losses, low switching frequency, and poor high-temperature performance. Given the material properties of GaN, devices can exhibit a given turn-on resistance and breakdown voltage with a minute fraction of the area of a Si device—all with higher switching speeds. While n-type doping of III-nitrides has reached a high level of maturity and integration into commercialized products, p-type doping presents distinct challenges that requires continued research. The success of the next generation of GaN-based power devices depends in large part on the ability to form selective-area p-type regions, which can be accomplished by ion implantation of the prevailing acceptor dopant, Mg. This approach is relatively inexpensive and is a fab-ready processing step for creation of current-generation power electronic devices including BJTs, IGBTs, and diode rectifiers. However, ion implantation induces lattice damage and creates point defects within the material, which negatively impacts dopant activation annealing. These forms of damage may be removed by high-temperature annealing, however short annealing time and higher pressures are required to prevent decomposition of the material. We approach this challenge through an innovative combination of approaches toward obtaining high-conductivity p-type GaN. By the novel employment of a gyrotron, a microwave source that can be shaped into a highly localized beam to dissipate very high power in a semiconductor target, a wafer of GaN can be annealed in short (<2 s) ‘pulses’ at temperatures of 1450 °C (2640 °F) to achieve p-type activation without appreciable surface degradation. The p-type activation of Mg-implanted GaN can be further aided by co-implantation with N, believed to enhance the removal of point defects that compensate the acceptors dopant. Elevation to such high temperatures for even a short time induces degradation of GaN-to overcome, we have implemented additional measures to retain high-quality crystalline GaN while removing implant damage and activating Mg dopant. A bilayer cap of AlN is deposited by a combination of chemical and physical vapor deposition to protect the GaN surface. To prevent decomposition of GaN loss of N during annealing, a custom vessel—to contain pressurized N while allowing injection of microwave energy—is employed. By detailed study of the interaction of various measures that are implemented to protect the sample surface (GaN), activate the implanted dopants, and mitigate damage of the lattice heating damage during annealing, we are able to achieve p-type GaN and observe recovery of the crystal and removal of substantial density of defects created during the implantation process.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Power system monitoring and source control of the Space Station Freedom DC power system testbed

Unlike a terrestrial electric utility which can purchase power from a neighboring utility, the Space Station Freedom (SSF) has strictly limited energy resources; as a result, source control, system monitoring, system protection, and load management are essential to the safe and efficient operation of the SSF Electric Power System (EPS). These functions are being evaluated in the DC Power Management and Distribution (PMAD) Testbed which NASA LeRC has developed at the Power System Facility (PSF) located in Cleveland, Ohio. The testbed is an ideal platform to develop, integrate, and verify power system monitoring and control algorithms. State Estimation (SE) is a monitoring tool used extensively in terrestrial electric utilities to ensure safe power system operation. It uses redundant system information to calculate the actual state of the EPS, to isolate faulty sensors, to determine source operating points, to verify faults detected by subsidiary controllers, and to identify high impedance faults. Source control and monitoring safeguard the power generation and storage subsystems and ensure that the power system operates within safe limits while satisfying user demands with minimal interruptions. System monitoring functions, in coordination with hardware implemented schemes, provide for a complete fault protection system. The objective of this paper is to overview the development and integration of the state estimator and the source control algorithms.

Kimnach, Greg L.↗

Power system monitoring and source control of the Space Station Freedom dc-power system testbed

Unlike a terrestrial electric utility which can purchase power from a neighboring utility, the Space Station Freedom (SSF) has strictly limited energy resources; as a result, source control, system monitoring, system protection, and load management are essential to the safe and efficient operation of the SSF Electric Power System (EPS). These functions are being evaluated in the dc Power Management and Distribution (PMAD) Testbed which NASA LeRC has developed at the Power System Facility (PSF) located in Cleveland, Ohio. The testbed is an ideal platform to develop, integrate, and verify power system monitoring and control algorithms. State Estimation (SE) is a monitoring tool used extensively in terrestrial electric utilities to ensure safe power system operation. It uses redundant system information to calculate the actual state of the EPS, to isolate faulty sensors, to determine source operating points, to verify faults detected by subsidiary controllers, and to identify high impedance faults. Source control and monitoring safeguard the power generation and storage subsystems and ensure that the power system operates within safe limits while satisfying user demands with minimal interruptions. System monitoring functions, in coordination with hardware implemented schemes, provide for a complete fault protection system. The objective of this paper is to overview the development and integration of the state estimator and the source control algorithms.

Kimnach, Greg L.↗

Stationary Plasma Thruster Plume Emissions

The emission spectrum from a xenon plasma produced by a Stationary Plasma Thruster provided by the Ballistic Missile Defense Organization (BMDO) was measured. Approximately 270 individual Xe I, Xe II, and XE III transitions were identified. A total of 250 mW of radiated optical emission was estimated from measurements taken at the thruster exit plane. There was no evidence of erosion products in the emission signature. Ingestion and ionization of background gas at elevated background pressure was detected. The distribution of excited states could be described by temperatures ranging from fractions of 1 eV to 4 eV with a high degree of uncertainty due to the nonequilibrium nature of this plasma. The plasma was over 95 percent ionized at the thruster exit plane. Between 10 and 20 percent of the ions were doubly charged. Two modes of operation were identified. The intensity of plasma emission increased by a factor of two during operation in an oscillatory mode. The transfer between the two modes of operation was likely related to unidentified phenomena occurring on a time scale of minutes.

Manzella, David H.↗

Hot, cold, or just right? An infrared biometric sensor to improve occupant comfort and reduce overcooling in buildings via closed-loop control

To improve occupant comfort and save energy in buildings, we have developed a closed-loop air conditioning (AC) sensor-controller that predicts occupant thermal sensation from the thermographic measurement of skin temperature distribution, then uses this information to reduce overcooling (cooling-energy overuse that discomforts occupants) by regulating AC output. Taking measures to protect privacy, it combines thermal-infrared (TIR) and color (visible spectrum) cameras with machine vision to measure the skin-surface temperature profile. Since the human thermoregulation system uses skin blood flow to maintain thermoneutrality, the distribution of skin temperature can be used to predict warm, neutral, and cool thermal states. We conducted a series of human-subject thermal-sensation trials in cold-to-hot environments, measuring skin temperatures and recording thermal sensation votes. We then trained random-forest classification machine-learning models (classifiers) to estimate thermal sensation from skin temperatures or skin-temperature differences. The estimated thermal sensation was input to a proportional integral (PI) control algorithm for the AC, targeting a sensation level between neutral and warm. Our sensor-controller includes a sensor assembly, server software, and client software. The server software orients the cameras and transmits images to the client software, which in turn assesses occupant skin temperature distribution, estimates occupant thermal sensation, and controls AC operation. A demonstration conducted in a conference room in an office building near Houston, TX showed that our system reduced overcooling, decreasing AC load by 42% when the room was occupied while improving occupant comfort (fraction of “comfortable” votes) by 15 percentage points.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Bayesian Structural Time Series for Behind-the-Meter Photovoltaic Disaggregation

Distributed photovoltaic (PV) generation often occurs ``behind the meter": a grid operator can only observe the net load, which is the sum of the gross load and distributed PV generation. This lack of observability poses a challenge to system operation at both bulk level and distribution level. The lack of real-time or near-future disaggregated estimates of gross load and PV generation will lead to over scheduling of energy production and regulation reserves, reliability constraints violations, wear and tear of controller devices, and potentially cascading failures of a system. In this paper we propose the use of a Bayesian Structural Time Series (BSTS) model with local solar irradiance measurements to disaggregate the summed PV generation and gross load signals at a downstream measurement site. BSTSs are a highly expressive model class that blends classic time series models with the powerful Bayesian state space estimation framework. Disaggregation is done probabilistically, which automatically quantifies the uncertainties of the estimated PV generation and gross load consumption. Depending on the data availability in real-time, it can be used to disaggragate PV and gross load at customer site, or can be used at the feeder level. In this paper, we focus on solving the problem at feeder level. We compare the performance of a BSTS model as well as a handful of state-of-the-art methods on a Pecan Street AMI dataset, using the National Solar Radiation Database (NSRDB) to estimate local irradiance.

bayesian structural time series↗

An estimation of aerodynamic forces and moments on an airplane model under steady state spin conditions

The problem of estimating the aerodynamic characteristics of a wing body configuration of a spinning airplane at extreme angles of attack and spin rates is approached by the application of strip theory. Semiempirical methods are used to estimate the aerodynamic force distributions over each component using static, nonrotational wind-tunnel test data. To these predictions, corrections are applied to account for the primary rotational flow effects which are significant at large spin rates. The results of this modified strip theory are shown to be in good agreement with spin tunnel rotary balance test data.

Pamadi, B. N.↗

A Nonstationary and Non-Gaussian Moving Average Model for Solar Irradiance

Historically, power has flowed from large power plants to customers. Increasing penetration of distributed energy resources such as solar power from rooftop photovoltaic has made the distribution network a two-way-street with power being generated at the customer level. The incorporation of renewables introduces additional uncertainty and variability into the power grid. Distribution network operation studies are being adapted to include renewables; however, such studies require high quality solar irradiance data that adequately reflect realistic meteorological variability. Data from satellite-based products are spatially complete, but temporally coarse, whereas solar irradiances exhibit high frequency variation at very fine timescales. We propose a new stochastic method for temporally downscaling global horizontal irradiance (GHI) to 1 min resolution, but we do not consider the spatial aspect due to limited availability of the in situ irradiance measurements. Solar irradiance's first and second-order structures vary diurnally and seasonally, and our model adapts to such nonstationarity. Empirical irradiance data exhibits highly non-Gaussian behavior; we develop a nonstationary and non-Gaussian moving average model that is shown to capture realistic solar variability at multiple timescales. We also propose a new estimation scheme based on Cholesky factors of empirical autocovariance matrices, bypassing difficult and inaccessible likelihood-based approaches. The model is demonstrated for a case study of three locations that are located in diverse climates through the United States. The model is compared against competitors from the literature and is shown to provide better uncertainty and variability quantification on testing data.

Cholesky factor↗

Rapid attribution analysis of the extraordinary heat wave on the Pacific coast of the US and Canada in June 2021

Towards the end of June 2021, temperature records were broken by several degrees Celsius in several cities in the Pacific Northwest areas of the US and Canada, leading to spikes in sudden deaths and sharp increases in emergency calls and hospital visits for heat-related illnesses. Here we present a multi-model, multi-method attribution analysis to investigate the extent to which human-induced climate change has influenced the probability and intensity of extreme heat waves in this region. Based on observations, modelling and a classical statistical approach, the occurrence of a heat wave defined as the maximum daily temperature (TXx) observed in the area 45–52 °N, 119–123 °W, was found to be virtually impossible without human-caused climate change. The observed temperatures were so extreme that they lay far outside the range of historical temperature observations. This makes it hard to state with confidence how rare the event was. Using a statistical analysis that assumes that the heat wave is part of the same distribution as previous heat waves in this region led to a first-order estimation of the event frequency of the order of once in 1000 years under current climate conditions. Using this assumption and combining the results from the analysis of climate models and weather observations, we found that such a heat wave event would be at least 150 times less common without human-induced climate change. Also, this heat wave was about 2 °C hotter than a 1-in-1000-year heat wave would have been in 1850–1900, when global mean temperatures were 1.2 °C cooler than today. Looking into the future, in a world with 2 °C of global warming (0.8 °C warmer than today), a 1000-year event would be another degree hotter. Our results provide a strong warning: our rapidly warming climate is bringing us into uncharted territory with significant consequences for health, well-being and livelihoods. Adaptation and mitigation are urgently needed to prepare societies for a very different future.

54 ENVIRONMENTAL SCIENCES↗

A study of fracture mechanisms in ATD roller bearing

The purpose was to investigate how microstructures, especially anisotropy, affects internal stresses and the overall mechanical response of bearings. Samples with the stress axis along the aligned carbide direction possessed high modulus values compared to those with their axis perpendicular to the carbide stringers. The difference in the modulus was found to be more than a factor of two. A series of experiments was conducted on rolled samples to further investigate this effect; the two sets of results were consistent with each other. The degree of anisotropy of the microstructure in terms of the carbide and matrix orientations was determined using x-ray diffraction. The stress state determination was conducted using neutron diffraction. It was found that there was little variation in the distribution of the internal stresses amongst different samples, indicating that small changes in the processing and geometrical conditions did not result in significant variations in the internal stress. A nominal tensile hoop stress of 39 ksi was obtained for the inner raceway. Furthermore, during the course of the investigation, it was apparent that there was a need to determine the stress state induced by the shrink fitting process. Therefore, a series of experiments was conducted using strain gages to identify the stress distribution in a shrink fitting process in three different types of geometries. Correlations were obtained to estimate the highest stress values in the outer and inner groove geometry. A finite element program based on the ANSYS system was developed to compute the stress distribution in the inner raceway geometry. This analysis indicates that the highest tensile stress in the system occurs at the ID of the ring with a stress value of over 5 times that of the applied radial stress. Results from all these facets were correlated with one another. It appears that the material does not fail as a result of any one single factor, but results from a combination of the various factors investigated.

Zee, Ralph H.↗

A Method to Retrieve Rainfall Rate Over Land from TRMM Microwave Imager Observations

Over tropical land regions, rain rate maxima in mesoscale convective systems revealed by the Precipitation Radar (PR) flown on the Tropical Rainfall Measuring Mission (TRMM) satellite are found to correspond to thunderstorms, i.e., Cbs. These Cbs are reflected as minima in the 85 GHz brightness temperature, T85, observed by the TRMM Microwave Imager (TMI) radiometer. Because the magnitude of TMI observations do not discriminate satisfactorily convective and stratiform rain, we developed here a different TMI discrimination method. In this method, two types of Cbs, strong and weak, are inferred from the Laplacian of T85 at minima. Then, to retrieve rain rate, where T85 is less than 270 K, a weak (background) rain rate is deduced using T85 observations. Furthermore, over a circular area of 10 km radius centered at the location of each T85 minimum, an additional Cb component of rain rate is added to the background rain rate. This Cb component of rain rate is estimated with the help of (T19-T37) and T85 observations. Initially, our algorithm is calibrated with the PR rain rate measurements from 20 MCS rain events. After calibration, this method is applied to TMI data taken from several tropical land regions. With the help of the PR observations, we show that the spatial distribution and intensity of rain rate over land estimated from our algorithm are better than those given by the current TMI-Version-5 Algorithm. For this reason, our algorithm may be used to improve the current state of rain retrievals on land.

Prabhakara, C.↗