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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗
NASA NTRS2011

Local Estimators for Spacecraft Formation Flying

A formation estimation architecture for formation flying builds upon the local information exchange among multiple local estimators. Spacecraft formation flying involves the coordination of states among multiple spacecraft through relative sensing, inter-spacecraft communication, and control. Most existing formation flying estimation algorithms can only be supported via highly centralized, all-to-all, static relative sensing. New algorithms are needed that are scalable, modular, and robust to variations in the topology and link characteristics of the formation exchange network. These distributed algorithms should rely on a local information-exchange network, relaxing the assumptions on existing algorithms. In this research, it was shown that only local observability is required to design a formation estimator and control law. The approach relies on breaking up the overall information-exchange network into sequence of local subnetworks, and invoking an agreement-type filter to reach consensus among local estimators within each local network. State estimates were obtained by a set of local measurements that were passed through a set of communicating Kalman filters to reach an overall state estimation for the formation. An optimization approach was also presented by means of which diffused estimates over the network can be incorporated in the local estimates obtained by each estimator via local measurements. This approach compares favorably with that obtained by a centralized Kalman filter, which requires complete knowledge of the raw measurement available to each estimator.

Fathpour, Nanaz↗
DOE OSTI2025

Preliminary Assessment of Renewable Energy Options at the Sravasti Abbey [Slides]

Through the U.S. Department of Energy's Energy to Communities (E2C) program, NREL, other national laboratory experts, and select organizations provide Expert Match - free, short-term technical assistance to address near-term clean energy challenges and questions. Expert Match is for community stakeholders who have decision-making power or influence in their community but need access to additional clean energy expertise to inform key upcoming decisions. This Expert Match request supported the Sravasti Abbey in Washington a preliminary assessment of distributed energy technology financial feasibility for their site. The outcome consists of low estimated savings from these technologies due to their already low energy consumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Scalable quantum computational science: A perspective from block-encodings and polynomial transformations

Significant developments made in quantum hardware and error correction recently have been driving quantum computing toward practical utility. However, gaps remain between abstract quantum algorithmic development and practical applications in computational sciences. In this perspective article, we propose several properties that scalable quantum computational science methods should possess. We further discuss how block-encodings and polynomial transformations can potentially serve as a unified framework with the desired properties. Recent advancements on these topics are presented, including the construction and assembly of block-encodings, and various generalizations of quantum signal processing (QSP) algorithms to perform polynomial transformations. The scalability of QSP methods on parallel and distributed quantum architectures is also highlighted. Promising applications in simulation and observable estimation in chemistry, physics, and optimization problems are presented. We hope this perspective serves as a gentle introduction to state-of-the-art quantum algorithms for the computational science community and inspires future development of scalable quantum computational science methodologies that bridge theory and practice.

Bayesian inference↗
DOE OSTI2022

Accelerating astronomical and cosmological inference with preconditioned Monte Carlo

ABSTRACT We introduce preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilizes a Normalizing Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison, respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster, respectively, than nested sampling (NS). We found that in higher dimensional applications, the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Integrated Demand Management: The FCA Balancing Algorithm

The Collaborative Trajectory Options Program (CTOP) is a new NextGen function that resides within the FAA’s Traffic Flow Management System (TFMS) software. CTOP was designed to be used in the United States’ national airspace to control the distribution of air traffic demand, both temporally and geographically, through a set of Flow Constrained Areas (FCAs). Each FCA is assigned capacity values that represent an estimated upper limit of manageable demand, and CTOP controls to that capacity by using Expect Departure Clearance Time (EDCT) air traffic control clearances, and/or rerouting flights to an alternate FCA or out of the CTOP completely. CTOP’s decisions about which control options best meet system and user needs are based on the Trajectory Options Set (TOS) submitted by airline operators for each flight. The FCA Balancing Algorithm (FBA) described in this document is a proposed decision support capability that generates CTOP capacity entries for FCAs controlling traffic to a common downstream constraint. The FBA provides an equitable distribution of impact, and addresses the added complexity of assigning FCA capacity values for this specific type of use case, greatly simplifying the FAA air traffic manager’s task.

IDM↗
DOE OSTI2023

Viability of Reclaiming Municipal Wastewater for Potential Microalgae-Based Biofuel Production in the U.S.

Reclaimed municipal wastewater is a crucial component in biofuel production, especially in regions experiencing increasing freshwater scarcity. However, accurately estimating the potential for fuel production is challenging because of the uneven distribution of biofuel feedstock regions and wastewater treatment plants (WWTPs). This study assesses the viability of using reclaimed municipal water for algal biomass production in pond systems co-located with WWTPs under scenarios driven by biomass production and based on water transport logistics. We performed state- and county-level analysis of reclaimed water resources throughout the United States based on WWTP facility data. We overlaid these data onto estimated algae facility sites and examined the temporal resource availability to address seasonal variations in cultivation demand. Our findings reveal that 2694 billion liters per year of reclaimed water could potentially be used to produce 42.2 million metric tons (ash-free dry weight) of algal biomass, equivalent to 29.2 billion liters of renewable diesel equivalent (RDe). The use of reclaimed water would double current national water reuse and expand such reuse significantly in 455 counties across the United States. However, when we limit the construction of algae facilities to counties that can fully meet their water demand in order to minimize water transport burdens, the available supply decreases by 80%, to 512 billion liters, resulting in annual production of 12.2 billion liters of RDe, which still doubles current biodiesel production. Our analysis highlights the degree to which the location and flow of WWTPs and water transport affect the deployment of algae biofuel facilities and tradeoffs. These findings underscore the importance of improving the current WWTP infrastructure for reclaimed water reuse, especially in southern states.

09 BIOMASS FUELS↗
DOE OSTI2025

Uncertainties in the effects of organic aerosol coatings on polycyclic aromatic hydrocarbon concentrations and their estimated health effects

We used the CAM5 model to examine how different particle-bound polycyclic aromatic hydrocarbon (PAH) degradation approaches affect the spatial distribution of benzo(a)pyrene (BaP). Three approaches were evaluated: NOA (no effect of OA coatings state on BaP), shielded (viscous OA coatings shield BaP from oxidation under cool and dry conditions) and ROI-T (viscous OA coatings slow BaP oxidation in response to temperature and humidity). Results show that BaP concentrations vary seasonally, influenced by emissions, deposition, transport and degradation approach, all of which are influenced by meteorological conditions. All simulations predict higher population-weighted global average (PWGA) fresh BaP concentrations during December–January–February (DJF) compared to June–July–August (JJA), due to increased emissions from household activities and reduced removal processes during colder months. The shielded and ROI-T approaches, which account for OA coatings, result in 2–6 times higher BaP concentrations in DJF compared to NOA. The shielded simulation predicts the highest PWGA fresh BaP concentration (1.3 ng m −3 ), with 90 % of BaP protected from oxidation. In contrast, the ROI-T approach forecasts lower concentrations in middle to low latitudes, as it assumes less effective OA coatings under warmer, more humid conditions. Evaluations against observed BaP concentrations show the shielded approach performs best, with a normalized mean bias (NMB) within ± 20 %. The combined incremental lifetime cancer risk (ILCR) for both fresh and oxidized PAHs is similar across simulations, emphasizing the importance of considering both forms in health risk assessments. This study highlights the critical role of accurate degradation approaches in PAH modeling.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Data-Driven Mean-Corrected Recursive Estimation-Based Optimal DER Dispatch for Distribution System Voltage Control

Recent advances in smart inverters offer opportunities to mitigate adverse grid impacts caused by high penetrations of distributed photovoltaics (PV) in distribution grids, such as voltage violations. Here, this paper proposes a novel measurement-driven optimal power flow (OPF)-based distributed energy resource management system (DERMS) voltage regulation via recursive sensitivity estimation informed coordinated control of distributed PV inverters. The proposed approach leverages available grid and controllable DER measurements, eliminating reliance on system model information while being adaptive and robust to volatile operating conditions. A mean-corrected recursive ridge regression (MCRRR) algorithm is proposed for sensitivity estimation, continuously refining the sensitivity model through a closed-form solution. It effectively manages varying grid operating conditions, such as changes in power injections and topology reconfiguration, to facilitate a time-varying update of the Load Sensitivity Factors (LSF). The proposed approach is formulated as a linear programming (LP) problem and is thus scalable to larger-scale distribution systems. Its effectiveness and efficiency are demonstrated on a realistic distribution feeder with high PV penetrations in Southern California, USA.

14 SOLAR ENERGY↗
DOE OSTI2024

On the connection between least squares, regularization, and classical shadows

Classical shadows (CS) offer a resource-efficient means to estimate quantum observables, circumventing the need for exhaustive state tomography. Here, we clarify and explore the connection between CS techniques and least squares (LS) and regularized least squares (RLS) methods commonly used in machine learning and data analysis. By formal identification of LS and RLS ``shadows'' completely analogous to those in CS---namely, point estimators calculated from the empirical frequencies of single measurements---we show that both RLS and CS can be viewed as regularizers for the underdetermined regime, replacing the pseudoinverse with invertible alternatives. Through numerical simulations, we evaluate RLS and CS from three distinct angles: the tradeoff in bias and variance, mismatch between the expected and actual measurement distributions, and the interplay between the number of measurements and number of shots per measurement. Compared to CS, RLS attains lower variance at the expense of bias, is robust to distribution mismatch, and is more sensitive to the number of shots for a fixed number of state copies---differences that can be understood from the distinct approaches taken to regularization. Conceptually, our integration of LS, RLS, and CS under a unifying ``shadow'' umbrella aids in advancing the overall picture of CS techniques, while practically our results highlight the tradeoffs intrinsic to these measurement approaches, illuminating the circumstances under which either RLS or CS would be preferred, such as unverified randomness for the former or unbiased estimation for the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2011

Contributions of Precipitation and Soil Moisture Observations to the Skill of Soil Moisture Estimates in a Land Data Assimilation System

The contributions of precipitation and soil moisture observations to the skill of soil moisture estimates from a land data assimilation system are assessed. Relative to baseline estimates from the Modern Era Retrospective-analysis for Research and Applications (MERRA), the study investigates soil moisture skill derived from (i) model forcing corrections based on large-scale, gauge- and satellite-based precipitation observations and (ii) assimilation of surface soil moisture retrievals from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E). Soil moisture skill is measured against in situ observations in the continental United States at 44 single-profile sites within the Soil Climate Analysis Network (SCAN) for which skillful AMSR-E retrievals are available and at four CalVal watersheds with high-quality distributed sensor networks that measure soil moisture at the scale of land model and satellite estimates. The average skill (in terms of the anomaly time series correlation coefficient R) of AMSR-E retrievals is R=0.39 versus SCAN and R=0.53 versus CalVal measurements. The skill of MERRA surface and root-zone soil moisture is R=0.42 and R=0.46, respectively, versus SCAN measurements, and MERRA surface moisture skill is R=0.56 versus CalVal measurements. Adding information from either precipitation observations or soil moisture retrievals increases surface soil moisture skill levels by IDDeltaR=0.06-0.08, and root zone soil moisture skill levels by DeltaR=0.05-0.07. Adding information from both sources increases surface soil moisture skill levels by DeltaR=0.13, and root zone soil moisture skill by DeltaR=0.11, demonstrating that precipitation corrections and assimilation of satellite soil moisture retrievals contribute similar and largely independent amounts of information.

Reichle, Rolf H.↗
DOE OSTI2023

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗
DOE OSTI2020

Finding wombling boundaries in LHC data with Voronoi and Delaunay tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Toward Improved Regional Hydrological Model Performance Using State-Of-The-Science Data-Informed Soil Parameters

Accurate soil moisture and streamflow data are an aspirational need of many hydrologically relevant fields. Model simulated soil moisture and streamflow hold promise but models require validation prior to application. Calibration methods are commonly used to improve model fidelity but misrepresentation of the true dynamics remains a challenge. In this study, we leverage soil parameter estimates from the Soil Survey Geographic (SSURGO) database and the probability mapping of SSURGO (POLARIS) to improve the representation of hydrologic processes in the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) over a central California domain. Our results show WRF-Hydro soil moisture exhibits increased correlation coefficients ( r ), reduced biases, and increased Kling-Gupta Efficiencies (KGEs) across seven in situ soil moisture observing stations after updating the model's soil parameters according to POLARIS. Compared to four well-established soil moisture data sets including Soil Moisture Active Passive data and three Phase 2 North American Land Data Assimilation System land surface models, our POLARIS-adjusted WRF-Hydro simulations produce the highest mean KGE (0.69) across the seven stations. More importantly, WRF-Hydro streamflow fidelity also increases, especially in the case where the model domain is set up with SSURGO-informed total soil thickness. The magnitude and timing of peak flow events are better captured, r increases across nine United States Geological Survey stream gages, and the mean KGE across seven of the nine gages increases from 0.12 to 0.66. Our pre-calibration parameter estimate approach, which is transferable to other spatially distributed hydrological models, can substantially improve a model's performance, helping reduce calibration efforts and computational costs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Seasonal representation of extreme precipitation indices over the United States in CMIP6 present-day simulations

Realistically representing the present-day characteristics of extreme precipitation has been a challenge for global climate models, which is due in part to deficiencies in model resolution and physics, but is also due to a lack of consistency in gridded observations. In this study, we use three observation datasets, including gridded rain gauge and satellite data, to assess historical simulations from sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models. We separately evaluate summer and winter precipitation over the United States (US) with a comprehensive set of extreme precipitation indices, including an assessment of precipitation frequency, intensity and spatial structure. The observations exhibit significant differences in their estimates of area-average intensity distributions and spatial patterns of the mean and extremes of precipitation over the US. In general, the CMIP6 multi-model mean performs better than most individual models at capturing daily precipitation distributions and extreme precipitation indices, particularly in comparison to gauge-based data. Also, the representation of the extreme precipitation indices by the CMIP6 models is better in the summer than winter. Although the ‘standard’ horizontal-resolution can vary significantly across CMIP6 models, from ~0.7° to ~2.8°, we find that resolution is not a good indicator of model performance. Overall, our results highlight common biases in CMIP6 models and demonstrate that no single model is consistently the most reliable across all indices.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

A new method for predicting response in complex linear systems. II

The paper describes convergence and response studies in the low frequency range of complex systems, particularly with low values of damping of different distributions, and reports on the modification of the relaxation procedure required under these conditions. A new method is presented for response estimation in complex lumped parameter linear systems under random or deterministic steady state excitation. The essence of the method is the use of relaxation procedures with a suitable error function to find the estimated response; natural frequencies and normal modes are not computed. For a 45 degree of freedom system, and two relaxation procedures, convergence studies and frequency response estimates were performed. The low frequency studies are considered in the framework of earlier studies (Kayser and Bogdanoff, 1975) involving the mid to high frequency range.

Bogdanoff, J. L.↗
NASA NTRS2013

Distributed Prognostics and Health Management with a Wireless Network Architecture

A heterogeneous set of system components monitored by a varied suite of sensors and a particle-filtering (PF) framework, with the power and the flexibility to adapt to the different diagnostic and prognostic needs, has been developed. Both the diagnostic and prognostic tasks are formulated as a particle-filtering problem in order to explicitly represent and manage uncertainties in state estimation and remaining life estimation. Current state-of-the-art prognostic health management (PHM) systems are mostly centralized in nature, where all the processing is reliant on a single processor. This can lead to a loss in functionality in case of a crash of the central processor or monitor. Furthermore, with increases in the volume of sensor data as well as the complexity of algorithms, traditional centralized systems become for a number of reasons somewhat ungainly for successful deployment, and efficient distributed architectures can be more beneficial. The distributed health management architecture is comprised of a network of smart sensor devices. These devices monitor the health of various subsystems or modules. They perform diagnostics operations and trigger prognostics operations based on user-defined thresholds and rules. The sensor devices, called computing elements (CEs), consist of a sensor, or set of sensors, and a communication device (i.e., a wireless transceiver beside an embedded processing element). The CE runs in either a diagnostic or prognostic operating mode. The diagnostic mode is the default mode where a CE monitors a given subsystem or component through a low-weight diagnostic algorithm. If a CE detects a critical condition during monitoring, it raises a flag. Depending on availability of resources, a networked local cluster of CEs is formed that then carries out prognostics and fault mitigation by efficient distribution of the tasks. It should be noted that the CEs are expected not to suspend their previous tasks in the prognostic mode. When the prognostics task is over, and after appropriate actions have been taken, all CEs return to their original default configuration. Wireless technology-based implementation would ensure more flexibility in terms of sensor placement. It would also allow more sensors to be deployed because the overhead related to weights of wired systems is not present. Distributed architectures are furthermore generally robust with regard to recovery from node failures.

Goebel, Kai↗