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 1,009 records · Page 56

Secondary battery management

A method of managing a battery system, the battery system including at least one battery cell, at least one sensor configured to measure at least one characteristic of the battery cell, and a battery management system including a microprocessor and a memory, the method comprising receiving by the battery management system, from the at least one sensor at least one measured characteristic of the battery cell at a first time and at least one measured characteristic of the battery cell at a second time. The battery management system estimating, at least one state of the battery cell by applying a physics-based battery model, the physics based battery model being based on differential algebraic equations; and regulating by the battery management system, at least one of charging or discharging of the battery cell based on the at least one estimated state.

25 ENERGY STORAGE↗

Lunar orbit determination in the presence of unmodeled accelerations

A technique for estimating the state of an artificial satellite in the presence of unmodeled accelerations is presented. The unmodeled acceleration is approximated by a first-order Gauss-Markov sequence which can be separated into a timewise-correlated component and a purely random component. Using this approximation, a sequential procedure for estimating the position, velocity, and the unmodeled acceleration is developed. The method is evaluated by reducing range-rate observations obtained by tracking the Apollo 10 and 11 spacecraft during the lunar-orbit phase of the mission. Numerical results are presented which show that the observation residual pattern lies within the observation noise standard deviation. The values of the estimated components of the unmodeled acceleration are repeatable from orbit to orbit within a given mission and from mission to mission when the same ground track is covered. Finally, the variation in the radial component of the unmodeled acceleration shows a high correlation with the reported location of the lunar surface mascons.

Ingram, D. S.↗

Minimum Detectable Quantity Calculation for Radiation Portal Monitors

The minimum detectable quantity of a radiation portal monitor is the smallest amount of radioactive material that can be detected passing through the monitor (with a specified detection probability). The minimum detectable quantity is a function of many factors, including isotope, velocity, distance between pillars, height of portal and vehicle, source distribution and position, background radiation, background suppression, detector volume and positions, detector efficiency, decision metrics and algorithms, the presence of NORM (naturally occurring radioactive material), and alarm thresholds. Experimentally testing a radiation portal monitor to failure for all combinations of these factors is extremely time-consuming and often cost prohibitive. This document outlines a straightforward method for quickly estimating the minimum detectable quantity for moving sources over a large variety of conditions with a minimal number of static measurements. The document also describes a proof of concept software application that the International Atomic Energy Agency has designed based on the described method that can be used by Member States to estimate a monitor’s minimum detectable quantity.

Blessinger, Christopher S.↗

An Assessment of Technical Hydropower Potential at Non-Powered Dams in the United States

Historically, dams have been constructed for a variety of purposes, such as providing a more secure and reliable water supply, mitigating impacts from variations in river flow, allowing continuous navigability, and harnessing mechanical power. A relatively small portion of dams have been designed to store or regulate flows for the purpose of generating electricity (roughly 3% of nationally inventoried dams in the US and 17% of the dams included in the World Register of Dams). The remaining population of existing non-powered dams (NPDs) presents both an opportunity to generate renewable energy and a need to modernize aging infrastructure. This report describes an assessment of more than 2,600 NPDs in the US that have a collective potential of nearly 4 GW in new power capacity. Previous national-scale assessments were aimed at evaluating the theoretical maximum power potential at existing dams in the United States. These estimates were based on the best available information at the time for water availability, hydraulic head, and representative regional capacity factors. This study revisits a subset of 3,299 dams identified in the most recent theoretical resource assessment and uses more detailed and updated hydrologic data to produce estimates of technical potential. These improvements in data enable estimates that more realistically reflect what is physically possible given simple assumptions about the existing structure and constraints on flow and head (Figure 1).

13 HYDRO ENERGY↗

Lanczos algorithm for lattice QCD matrix elements

Recent work [M. L. Wagman, Lanczos, the transfer matrix, and the signal-to-noise problem, .] found that an analysis formalism based on the Lanczos algorithm allows energy levels to be extracted from Euclidean correlation functions with faster ground-state convergence than effective masses, convergent estimators for multiple states from a single correlator, and two-sided error bounds. After filtering out spurious eigenvalues and using outlier-robust estimators within a nested bootstrap framework, Lanczos estimators behave more like multistate fit results than effective masses—but without involving statistical fitting. We extend this formalism to the determination of matrix elements from three-point correlation functions and provide a physical picture of “spurious-state filtering” involving restriction to a Hermitian subspace. We demonstrate similar advantages for matrix elements as for spectroscopy through example applications to noiseless mock-data and (bare) forward matrix elements of the strange scalar current between both ground and excited states with the quantum numbers of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time Domain Tool Validation Using ARES I-X Flight Data

The ARES I-X vehicle was launched from NASA's Kennedy Space Center (KSC) on October 28, 2009 at approximately 11:30 EDT. ARES I-X was the first test flight for NASA s ARES I launch vehicle, and it was the first non-Shuttle launch vehicle designed and flown by NASA since Saturn. The ARES I-X had a 4-segment solid rocket booster (SRB) first stage and a dummy upper stage (US) to emulate the properties of the ARES I US. During ARES I-X pre-flight modeling and analysis, six (6) independent time domain simulation tools were developed and cross validated. Each tool represents an independent implementation of a common set of models and parameters in a different simulation framework and architecture. Post flight data and reconstructed models provide the means to validate a subset of the simulations against actual flight data and to assess the accuracy of pre-flight dispersion analysis. Post flight data consists of telemetered Operational Flight Instrumentation (OFI) data primarily focused on flight computer outputs and sensor measurements as well as Best Estimated Trajectory (BET) data that estimates vehicle state information from all available measurement sources. While pre-flight models were found to provide a reasonable prediction of the vehicle flight, reconstructed models were generated to better represent and simulate the ARES I-X flight. Post flight reconstructed models include: SRB propulsion model, thrust vector bias models, mass properties, base aerodynamics, and Meteorological Estimated Trajectory (wind and atmospheric data). The result of the effort is a set of independently developed, high fidelity, time-domain simulation tools that have been cross validated and validated against flight data. This paper presents the process and results of high fidelity aerospace modeling, simulation, analysis and tool validation in the time domain.

Hough, Steven↗

A Comparison of Filter-based Approaches for Model-based Prognostics

Model-based prognostics approaches use domain knowledge about a system and its failure modes through the use of physics-based models. Model-based prognosis is generally divided into two sequential problems: a joint state-parameter estimation problem, in which, using the model, the health of a system or component is determined based on the observations; and a prediction problem, in which, using the model, the stateparameter distribution is simulated forward in time to compute end of life and remaining useful life. The first problem is typically solved through the use of a state observer, or filter. The choice of filter depends on the assumptions that may be made about the system, and on the desired algorithm performance. In this paper, we review three separate filters for the solution to the first problem: the Daum filter, an exact nonlinear filter; the unscented Kalman filter, which approximates nonlinearities through the use of a deterministic sampling method known as the unscented transform; and the particle filter, which approximates the state distribution using a finite set of discrete, weighted samples, called particles. Using a centrifugal pump as a case study, we conduct a number of simulation-based experiments investigating the performance of the different algorithms as applied to prognostics.

prognostics↗

Estimation of unmodeled forces on a low-thrust space vehicle

The application of a sequential estimation algorithm, which compensates for random errors in the dynamic model, to the problem of estimating the state of a continuously thrusting solar electric propulsion space vehicle is investigated. The dynamic model errors, due to random anomalies in the propulsion system, are approximated successfully by both first order and second order Gauss-Markov processes to obtain a more accurate and stable orbit determination algorithm. The importance of correct dynamic and measurement modeling in achieving accurate estimates is demonstrated.

Tapley, B. D.↗

Studying QED 3 with radial quantization on the lattice: Free limit

To investigate the three-dimensional quantum electrodynamics in the radial quantization on the lattice, the lattice action is constructed and the free limit is studied on 𝑆 2 × ℝ. With the overlap fermion, it is numerically verified that the important symmetries of the theory can be realized on the lattice. The analytic correlators are derived and compared to the lattice results, which agree including the overall normalization. The 𝑂⁡(𝑎 2 ) scaling is confirmed toward the analytic value in the continuum limit, and the number of reproduced excited states is estimated heuristically for the first few refinement levels. Our study helps us identify the features of the theory that we can study on the icosahedral lattice without fine-tuning.

Boyle, Peter A. [Brookhaven National Laboratory (B↗

A comparison of the effects of small nonlinearities on several estimation schemes.

A variety of techniques are available for estimating the states of nonlinear dynamic systems from noisy data. These procedures are generally equivalent when applied to linear systems. This paper investigates the difference between several of these procedures in the presence of small dynamic and observational nonlinearities. In particular, it examines one least square batch processing algorithm, and three recursive algorithms similar to the Kalman filter. To first order, all the estimators have the same covariance. Expressions for the means, however, show that each estimator has a different bias. The examples presented show that the biases can be a strong function of such parameters as initial covariances and number of data points being considered.

Conrad, B.↗

Whose Gas is it anyway? Differentiating the Source of a Large Soil Vapor Plume beneath Two Adjacent Waste Sites - 20487

DOE contractor CH2M Hill Plateau Remediation Company is currently responsible for conducting groundwater contamination monitoring at several RCRA treatment, storage, and disposal units located on the Hanford Site in Richland, Washington State. The Nonradioactive Dangerous Waste Landfill treatment, storage, and disposal unit presents a distinct groundwater monitoring problem because of a large multi-contaminant soil vapor plume beneath it that is a likely source of low-level volatile organic compound groundwater contamination. Adjacent to Nonradioactive Dangerous Waste Landfill is the Solid Waste Landfill. Volatile organic compounds are inventory components of both the Nonradioactive Dangerous Waste Landfill and the Solid Waste Landfill. Therefore, it is possible that both sites could be contributing to the soil vapor plume. For regulatory purposes, it is important to differentiate which site is the primary contributor of volatile organic compounds to the plume. An approach was developed to identify the primary volatile organic compound source of the soil vapor plume beneath Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. The site conceptual model hypothesis of vapor-phase volatile organic compound transport to the dissolved phase in groundwater was tested by a simple mathematical model of vapor/liquid equilibrium concentrations at the groundwater/air interface. Once it was shown that vapor-phase volatile organic compound transport to groundwater was a valid conceptual model for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill, spatial and statistical methods were used to determine the primary site contributing to the majority of volatile organic compounds to the soil vapor plume. Average groundwater chloroform, tetrachloroethene, and trichloroethene concentrations from Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring network wells were plotted on maps of the facilities and immediate vicinities and compared to soil vapor sampling probe locations. Principal component analysis and mixing ratios were used to identify source contributions of each treatment, storage, and disposal unit to the plume. Results of the vapor/liquid equilibrium concentrations mathematical model showed that transport phenomena outweigh steady-state equilibria. Estimated vapor/liquid equilibrium concentrations were considerably lower than soil vapor measurements. The results indicate that dynamic vadose zone and groundwater factors such as decreased vapor concentrations with depth, vapor dilution from dispersion in the vadose zone, and advective and diffusional volatile organic compound dilution in groundwater result in groundwater volatile organic compound concentrations much less than would be measured under steady-state equilibrium conditions. Site source contribution differentiation by principal component analysis and mixing ratios was inconclusive using actual soil gas data because of the similarity in concentration values in both datasets for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. Similar data populations suggest mixing of the vapor contributions from both sites by dispersion through the soil matrix pore spaces. However, when groundwater volatile organic compound data were compared between the Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring networks, Solid Waste Landfill mean concentrations were higher, suggesting more vapor-phase volatile organic compound transport to groundwater at those locations. Simulated volatile organic compound soil vapor and groundwater datasets created to test the methods developed for this study show that the method can be successful in source differentiation when significantly different datasets are compared. This paper will describe a method of testing a conceptual model for vapor-phase contaminant transport to groundwater and for differentiating site sources of contaminants comprising a mixed-constituent soil vapor plume. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Patterns, drivers, and a predictive model of dam removal cost in the United States

Given the burgeoning dam removal movement and the large number of dams approaching obsolescence in the United States, cost estimating data and tools are needed for dam removal prioritization, planning, and execution. We used the list of removed dams compiled by American Rivers to search for publicly available reported costs for dam removal projects. Total cost information could include component costs related to project planning, dam deconstruction, monitoring, and several categories of mitigation activities. We compiled reported costs from 455 unique sources for 668 dams removed in the United States from 1965 to 2020. The dam removals occurred within 571 unique projects involving 1–18 dams. When adjusted for inflation into 2020 USD, cost of these projects totaled $\$1.522$ billion, with per-dam costs ranging from $\$1$ thousand (k) to $\$268.8$ million (M). The median cost for dam removals was $\$157$k, $\$823$k, and $\$6.2$M for dams that were< 5 m, between 5–10 m, and > 10 m in height, respectively. Geographic differences in total costs showed that northern states in general, and the Pacific Northwest in particular, spent the most on dam removal. The Midwest and the Northeast spent proportionally more on removal of dams less than 5 m in height, whereas the Northwest and Southwest spent the most on larger dam removals > 10 m tall. We used stochastic gradient boosting with quantile regression to model dam removal cost against potential predictor variables including dam characteristics (dam height and material), hydrography (average annual discharge and drainage area), project complexity (inferred from construction and sediment management, mitigation, and post-removal cost drivers), and geographic region. Dam height, annual average discharge at the dam site, and project complexity were the predominant drivers of removal cost. The final model had an R 2 of 57% and when applied to a test dataset model predictions had a root mean squared error of $\$5.09$M and a mean absolute error of $\$1.45$M, indicating its potential utility to predict estimated costs of dam removal. We developed a R shiny application for estimating dam removal costs using customized model inputs for exploratory analyses and potential dam removal planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Multi-area parameter error identification for large power systems

Power grid model parameters may contain errors due to various reasons. Detecting and correcting parameter errors typically requires significant computational effort due to the size and complexity of the parameter database. While the normalized Lagrange multiplier (NLM) method can effectively detect, identify and correct parameter errors, its computational burden could rapidly grow with increasing system size. This paper addresses this issue by proposing a multi-area parameter error identification method. Each area has its own outlier detection tool for detecting the incorrect parameters and measurements within the area. On the other hand, due to the reduced redundancy at area boundaries, parameter errors on branches incident to boundary buses may not be detected. Such errors are subsequently detected by a coordination level estimator completing the system-wide parameter detection procedure. In conclusion, performance of the developed method is demonstrated using the IEEE 118-bus and 2000-bus Texas synthetic systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Distributed quantum sensing of multiple phases with fewer photons

Abstract Distributed quantum metrology has drawn intense interest as it outperforms the optimal classical counterparts in estimating multiple distributed parameters. However, most schemes so far have required entangled resources consisting of photon numbers equal to or more than the parameter numbers, which is a fairly demanding requirement as the number of nodes increases. Here, we present a distributed quantum sensing scenario in which quantum-enhanced sensitivity can be achieved with fewer photons than the number of parameters. As an experimental demonstration, using a two-photon entangled state, we estimate four phases distributed 3 km away from the central node, resulting in a 2.2 dB sensitivity enhancement from the standard quantum limit. Our results show that the Heisenberg scaling can be achieved even when using fewer photons than the number of parameters. We believe our scheme will open a pathway to perform large-scale distributed quantum sensing with currently available entangled sources.

Science & Technology - Other Topics↗

Nearly optimal state preparation for quantum simulations of lattice gauge theories

Here, we present several improvements to the recently developed ground-state preparation algorithm based on the quantum eigenvalue transformation for unitary matrices (QETU), apply this algorithm to a lattice formulation of U(1) gauge theory in (2+1) dimensions, as well as propose an alternative application of QETU, a highly efficient preparation of Gaussian distributions. The QETU technique was originally proposed as an algorithm for nearly optimal ground-state preparation and ground-state energy estimation on early fault-tolerant devices. It uses the time-evolution input model, which can potentially overcome the large overall prefactor in the asymptotic gate cost arising in similar algorithms based on the Hamiltonian input model. We present modifications to the original QETU algorithm that significantly reduce the cost for the cases of both exact and Trotterized implementation of the time evolution circuit. We use QETU to prepare the ground state of a U(1) lattice gauge theory in two spatial dimensions, explore the dependence of computational resources on the desired precision and system parameters, and discuss the applicability of our results to general lattice gauge theories. We also demonstrate how the QETU technique can be utilized for preparing Gaussian distributions and wave packets in a way which outperforms existing algorithms for as little as n q ≳ 2–5 qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Stabilizer Framework for the Contextual Subspace Variational Quantum Eigensolver and the Noncontextual Projection Ansatz

Quantum chemistry is a promising application for noisy intermediate-scale quantum (NISQ) devices. However, quantum computers have thus far not succeeded in providing solutions to problems of real scientific significance, with algorithmic advances being necessary to fully utilize even the modest NISQ machines available today. We discuss a method of ground state energy estimation predicated on a partitioning of the molecular Hamiltonian into two parts: one that is noncontextual and can be solved classically, supplemented by a contextual component that yields quantum corrections obtained via a Variational Quantum Eigensolver (VQE) routine. This approach has been termed Contextual Subspace VQE (CS-VQE); however, there are obstacles to overcome before it can be deployed on NISQ devices. The problem we address here is that of the ansatz, a parametrized quantum state over which we optimize during VQE; it is not initially clear how a splitting of the Hamiltonian should be reflected in the CS-VQE ansätze. We propose a “noncontextual projection” approach that is illuminated by a reformulation of CS-VQE in the stabilizer formalism. This defines an ansatz restriction from the full electronic structure problem to the contextual subspace and facilitates an implementation of CS-VQE that may be deployed on NISQ devices. We validate the noncontextual projection ansatz using a quantum simulator and demonstrate chemically precise ground state energy calculations for a suite of small molecules at a significant reduction in the required qubit count and circuit depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of Disturbance History on Forest Carbon Stocks and Fluxes: Merging Satellite Disturbance Mapping with Forest Inventory Data in a Carbon Cycle Model Framework

Forest carbon stocks and fluxes are highly dynamic following stand-clearing disturbances from severe fire and harvest and this presents a significant challenge for continental carbon budget assessments. In this work we use forest inventory data to parameterize a carbon cycle model to represent post-disturbance carbon trajectories of carbon pools and fluxes for specific forest types growing in high and low site productivity class settings. We then apply these trajectories to landscapes and regions based on forest age distributions derived from either the FIA data or from Landsat time series stacks (1985-2006) for 54 representative scenes throughout most of the conterminous United States. We estimate the net carbon uptake in forests caused by post-disturbance growth and decomposition ("regrowth sink") for forested regions across the country. At the landscape scale, the prevailing condition of positive net ecosystem productivity (NEP) is in stark contrast to local patches with large sources, particularly in the west where fires and clear cuts create contiguous disturbed patches. At the continental scale, regional differences in disturbance rates reflect management patterns of high disturbance rates in the Southeastern and South Central states, and lower disturbance rates in the Northeast and Northern Lakes States. Despite low contemporary disturbance rates in the Northeast and Northern Lakes States (0.61 and 0.74% y(exp −1)), the regrowth sink there remains of moderate to large strength (88 and 57 g C m(exp −2) y(exp −1)) owing to the continued legacy from historical clearing. Large regrowth sinks are also found in the Southeast, South Central, and Pacific Southwest regions (85, 86, and 95 g C m(exp −2) y(exp −1)) where disturbance rates also tend to be higher (1.59, 1.38, and 0.93% y(exp −1)). Overall, the Landsat-derived disturbance rates are elevated relative to FIA-derived rates (1.19 versus 0.93% y(exp −1)) particularly for western regions. The differences only modestly adjust regional- and continental-scale carbon budgets, reducing NEP from forest regrowth by about 8%.

carbon↗

Evaluating the photophysical and photochemical characteristics of green-emitting cerium(III) mono-cyclooctatetraenide complexes

The photophysical and photochemical reactivity of the organometallic mono(cyclooctatetraenide)-Ce(III) complexes: [(C 8 H 8 )Ce(μ-X)(THF) 2 ] 2 X = O 3 SCF 3 - (1), Cl - (2), were studied. Here, in the course of these studies, a new polymorph of [(C 8 H 8 )Ce(μ-O 3 SCF 3 )(THF) 2 ] 2 (1) was described. Photoluminescence (PL) studies of 1 and 2 showed characteristic excitation and emission bands in the visible region, with bright green light emission under light irradiation from states corresponding to the 5d → 4f transition. Complexes 1 and 2 exhibit relatively long lifetime of 205.4 ± 0.2 and 145.8 ± 0.4 ns respectively. And the quantum yields (Φ) for 1 and 2 are 0.180 and 0.068 respectively. Electrochemical studies were performed on complexes 1 and 2 with a reversible Ce(III)/Ce(IV) redox couple recorded at $E_{1/2}$ = -1.53 V versus the Fc/Fc + for 2. Complex 1 shows an irreversible Ce(III/IV) oxidation wave. Complexes 1 and 2 revealed strongly-reducing, estimated, excited state reduction potentials ($E^*_{1/2}$) of -3.31 and -4.02 V versus the Fc/Fc + respectively, with small Stokes shifts of 0.12 eV. With their associated relatively long lifetimes, small Stokes shifts and large, negative $E^*_{1/2}$ values, both complexes were evaluated as potential photosensitisers for halogen atom transfer (XAT) using a test reaction of the dehalogenation of benzyl chloride.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗