Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed 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,243 records · Page 69

Synergistic Retrievals of Ice in High Clouds From Elastic Backscatter Lidar, Ku-band Radar and Submillimeter Wave Radiometer Observations

In this study, we investigate the synergy of elastic backscatter lidar, Ku-band radar, and sub-millimeter-wave radiometer measurements in the retrieval of ice from satellite observations. The synergy is analyzed through the generation of a large dataset of IceWater Content (IWC) profiles and simulated lidar, radar and radiometer observations. The characteristics of the instruments e.g. frequencies, sensitivities, etc. are set based on the expected characteristics of instruments of the Atmosphere Observing System (AOS) mission. A hold-out validation methodology is used to assess the accuracy of the IWC profiles retrieved from various combinations of observations from the three instruments. Specifically, the IWC and associated observations are randomly divided into two datasets, one for training and the other for evaluation. The training dataset is used to train the retrieval algorithm, while the evaluation dataset is used to assess the retrieval performance. The dataset of IWC profiles is derived from CloudSat reflectivity and CALIOP lidar observations. The retrieval of the ice water content IWC profiles from the computed observations is achieved in two steps. In the first step, a class, out of 18 potential classes characterized by different vertical distribution of IWC, is estimated from the observations. The 18 classes are predetermined based on the k-Means clustering algorithm. In the second step, the IWC profile is estimated using an Ensemble Kalman Smoother (EKS) algorithm that uses the estimated class as a priori information. The results of the study show that the synergy of lidar, radar, and radiometer observations is significant in the retrieval of the IWC profiles. Nevertheless, it should be mentioned that this synergy was found under idealized conditions, and additional work might be required to materialize it in practice. The inclusion of the lidar backscatter observations in the retrieval process has a larger impact on the retrieval performance than the inclusion of the radar observations. As ice clouds have a significant impact on atmospheric radiative processes, this work is relevant to ongoing efforts to reduce uncertainties in climate analyses and projections.

Mircea Grecu↗

Recommended Practice for Use of ExB Probes in Electric Propulsion Testing

An ExB probe is an electric propulsion plume diagnostic that has at its core a Wien filter. Ions entering the instrument experience perpendicular electric and magnetic fields and are filtered based on their velocity. While the ExB probe can be used to estimate the ion velocity distribution function, in the electric propulsion community, it is commonly used to measure species fractions. There are numerous probe designs, implementation and operation procedures, and data analysis approaches described and in use across the community. This paper provides recommendations and descriptions of best practices for design, implementation, and data analysis for the ExB probe, with a particular emphasis on Hall and ion thrusters. This work contributes to the broader community goal to standardize the use of diagnostics in electric propulsion testing.

ExB probe↗

Constraints on neutrino oscillation parameters with the NOvA experiment

NOvA is a long-baseline neutrino oscillation experiment, which consists of two finely-segmented liquid-scintillator detectorsoperating 14 mrad off-axis from the NuMI muon neutrino beam. With an 810 km baseline, the combined measurements of muon neutrino disappearance and electron neutrino appearance allow the determination of neutrino oscillation unknowns, namely the mass hierarchy, the octant of the largest neutrino mixing angle, and the CP violating phase.In this dissertation, I present the joint analysis of $\nu_\mu\to\nu_\mu$ and $\nu_mu\to\nu_e$ oscillation data with an exposure of$8.85\times10^{20}$ protons on target on the 14 kton detector. It includes the estimation of neutrino energy distributions in the fardetector using near detector data constraints, details of the implementation of the simultaneous fit of \numu and \nue samples, andthe effect of systematic uncertainties on the measurement of oscillation parameters. I also discuss NOvA's projected sensitivity to determine the mass hierarchy and discover CP violation in future analyses with increased exposure and the addition ofantineutrino datasets.With 66 electron neutrino and 126 muon neutrino candidates, the bestfit to the data corresponds to the normal mass hierarchy, $\Delta m^2_{32}=2.44\times 10^{-3} eV^2/c^4$, $\sin^2\theta_{23}=0.56$,and $\delta_{CP} = 1.21\pi$, with the allowed regions $\Delta m^2_{32} \in[2.37,2.52]\times 10^{-3} eV^2/c^4$, $\sin^2\theta_{23} \in [0.43,0.51] \cup[0.52,0.60]$ and $\delta_{CP} \in [0,0.12\pi] \cup [0.91\pi,2\pi]$ at the68.3\% C.L. Our data disfavor maximal mixing by 0.8$\sigma$, $\delta_{CP}= \pi/2$ in the inverted hierarchy at higher than 3$\sigma$, and theentire inverted mass hierarchy at the 95\% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BETO 2021 Peer Review - Feedstock Evaluation and Biofuels Production Potential

The goal of this project is to provide foundational data, strategic analyses, and outreach related to waste-to-energy (WTE) resources to support further development of the WTE industry. It builds on previous project outcomes, e.g., estimates of quantity, geographic distribution, and prices of wet WTE resources (food waste, sludge, manure, waste fats and oils). The project provides better understanding of the WTE resource potential and economic viability to enable development of new technologies and support strategic decisions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Rain Fallspeeds and Rates Derived from Airborne Nadir-Pointing Doppler Radar Measurements

The use of vertical-incidence Doppler velocity in addition to radar reflectivity may yield information on drop size distribution and therefore result in better rainrate estimates. Doppler velocity can provide useful information on the raindrop size distribution. Doppler velocities from a zenith-pointing radar represent the sum of the mean reflectivity-weighted hydrometeor fallspeed and the vertical air motion. Dual-parameter rain estimation methods using the Doppler velocity, require that the latter can be removed, or is negligible. Atlas et al. (1972) derived relations between Doppler velocity, reflectivity, and rain rate assuming an exponential size distribution for rain. Ulbrich (1994) expanded on this work by deriving the relation between the Doppler velocity and the reflectivity assuming a Gamma size distribution. This distribution provides a more realistic representation of the small rain drops. To get accurate information on raindrop size distributions with the above method, the air motions must be removed from the observed Doppler velocities

Heymsfield, Gerald M.↗

Modeling Building Use and Population Distribution Opportunity Using Open Geosocial Data in Urban Areas

Knowing where people are and knowing how the land in urban areas is being used is critical for assessing human mobility patterns, monitoring population health, planning for transportation services, and responding to natural disasters and emergency situations. While geophysical phenomena like land cover are more easily observed through remote sensing, social phenomena like how a city block or building is being used by the population is generally unobservable through standard remote sensing practices. Current methods and land use data by their nature are also difficult to scale to national or global levels. While land use data are often available in large cities, there are abundant data gaps globally. Furthermore, land use data can often appear relatively coarse in its spatial and semantic resolution (i.e. 1km grid cells rather than building level, and "commercial" rather than "coffee shop"). In particular, there is a need to harvest scalable data and develop methods to create global land use maps. To this end, in this work we combine building footprints and Points of Interest (POI) to estimate land use at the building level, including its mixed-use intensity. Next, we propose a method for defining semantic temporal zones referring to Morning, Day, Evening, and Night times unique to a region, based on the cumulative hourly change in the POIs' hours of operation. We then use these POIs' hours of operation to estimate non-residential population distribution opportunity and land use change over a 24-hour period. We validate our results through manual checks of randomly sampled individual buildings using multiple sources. Finally, we discuss the impact of our results on identifying population distribution opportunity and the reliability of estimating building use.

Sparks, Kevin↗

Dynamically stratified Monte Carlo forecasting

A new method for performing Monte Carlo forecasts is introduced. The method, called dynamic stratification, selects initial perturbations based on a stratification of the error distribution. A simple implementation is presented in which the error distribution used for the stratification is estimated from a linear model derived from a large ensemble of 12-h forecasts with the full dynamic model. The stratification thus obtained is used to choose a small subsample of initial states with which to perform the dynamical Monte Carlo forecasts. Several test cases are studied using a simple two-level general circulation model with uncertain initial conditions. It is found that the method provides substantial reductions in the sampling error of the forecast mean and variance when compared to the more traditional approach of choosing the initial perturbations at random. The degree of improvement, however, is sensitive to the nature of the initial error distribution and to the base state. In practice the method may be viable only if the computational burden involved in obtaining an adequate estimate of the error distribution is shared with the data-assimilation procedure.

Schubert, Siegfried↗

Uncertainties for two-dimensional models of solar rotation from helioseismic eigenfrequency splitting

Observed solar p-mode frequency splittings can be used to estimate angular velocity as a function of position in the solar interior. Formal uncertainties of such estimates depend on the method of estimation (e.g., least-squares), the distribution of errors in the observations, and the parameterization imposed on the angular velocity. We obtain lower bounds on the uncertainties that do not depend on the method of estimation; the bounds depend on an assumed parameterization, but the fact that they are lower bounds for the 'true' uncertainty does not. Ninety-five percent confidence intervals for estimates of the angular velocity from 1986 Big Bear Solar Observatory (BBSO) data, based on a 3659 element tensor-product cubic-spline parameterization, are everywhere wider than 120 nHz, and exceed 60,000 nHz near the core. When compared with estimates of the solar rotation, these bounds reveal that useful inferences based on pointwise estimates of the angular velocity using 1986 BBSO splitting data are not feasible over most of the Sun's volume. The discouraging size of the uncertainties is due principally to the fact that helioseismic measurements are insensitive to changes in the angular velocity at individual points, so estimates of point values based on splittings are extremely uncertain. Functionals that measure distributed 'smooth' properties are, in general, better constrained than estimates of the rotation at a point. For example, the uncertainties in estimated differences of average rotation between adjacent blocks of about 0.001 solar volumes across the base of the convective zone are much smaller, and one of several estimated differences we compute appears significant at the 95% level.

Genovese, Christopher R.↗

National soil organic carbon estimates can improve global estimates

Global inventories of spatial and vertical distribution of soil organic carbon (SOC) stocks are being used in national and global initiatives targeted to mitigate climate change and land degradation impacts. Yet, national level high-resolution estimates of SOC stocks can be useful for improving the accuracy of global SOC inventories. We estimated spatially resolved SOC stocks of surface 0–30 cm and subsurface 30–100 cm layers at a spatial resolution of 30m in tropical Island, Sri Lanka using a legacy harmonized soil database of 122 soil profiles. The national estimates were compared with two global estimates derived from WISE30sec and SoilGrids250m. The tropical Island (land area=64,610 km2) occupying 0.03% of global land area showed a considerable heterogeneity in SOC stocks ranging from 2.0-342.5 Mg ha -1 and 2.7-391.7 Mg ha -1 in the surface and subsurface soil layers, respectively. We found, elevation, precipitation and slope angle as main environmental controllers of the spatial distribution of SOC stocks under tropical climate. Incorporating the pedogenic information (derived from soil series level legacy map, soil orders and suborders) with environmental controllers resulted in better regression models of predicting surface (R 2 =0.61) and subsurface (R 2 =0.81) SOC stocks. Geographically weighted regression kriging derived maps of SOC stocks revealed that 0–100 cm soil layer of the tropical Island stored 500 Tg C contributing for 0.04% of the global SOC stocks. The validation results of our estimates showed low Mean Estimation Error (MEE: surface -1.6 and subsurface -1.6 Mg ha -1 ) and Root Mean Square Error (RMSE: surface 29.5 and subsurface 24.9 Mg ha -1 ) indicating a low bias and satisfactory predictions. The relative improvement of the prediction accuracy of the SOC stocks of our geospatial estimates in the 0–30 cm layer in comparison to SoilGrids250m and WISE30sec data derived SOC stocks were 51.7% and 35.2%, respectively. The SOC stocks predictions of the 30-100 cm soil layer showed even better relative improvement compared to SoilGrids250m (78.4%) and WISE30sec (57.4%) SOC estimates. Compared to estimates of total SOC stocks resulted in this study, WISE30sec data derived SOC stock maps showed 30% over estimation of the C stock in surface 0-30 cm (332 Tg C) and 41% overestimation in 30-100 cm layer (343 Tg C). The over estimation of total SOC stocks by the SoilGrids250 SOC stocks map for the surface 0-30 cm layer was 122% (567 Tg C) and for the 30-100 cm layer it was 209% (750 Tg C). We conclude that the fusion of legacy soil information of SOC stocks with appropriate environmental covariates and pedogenic information derived from legacy area-class soil maps at national level can produce more accurate inventories of spatial and vertical distribution of SOC stocks. As a result, these national inventories have a great potential of upgrading global inventories of SOC stocks.

54 ENVIRONMENTAL SCIENCES↗

Regularization via f -Divergence: An Application to Multi-Oxide Spectroscopic Analysis

In this paper, we explore the application of convolutional neural networks (CNNs) for predicting the chemical composition of complex geologic samples in a simulated Martian atmospheric environment. Specifically, we aim to characterize oxide weight percentages (wt.%) of rock samples analyzed by remote Laser-Induced Breakdown Spectroscopy (LIBS), framing the problem as a multi-target regression task . Neural networks trained on LIBS spectra are prone to overfitting due to high spectral complexity, limited labeled data, and measurement noise. While regularization is critical for improving generalization, common methods (e.g., ℓ 2 regularization) impose constraints not directly tied to data distribution properties. We propose a novel regularization method based on a specific ƒ-divergence induced by a graph-based estimator, designed to constrain the distributional discrepancy between predictions and targets. This regularizer serves a dual purpose: (a) mitigating overfitting by enforcing a constraint on the distributional difference between predictions and noisy targets, and (b) acting as an auxiliary loss that penalizes large divergences. To enable backpropagation, we develop a differentiable approximation of this particular ƒ-divergence, making the method feasible for neural networks. Experiments on ChemCam and SuperCam LIBS calibration spectra show that mathematical equation-divergence regularization outperforms or matches standard regularization methods (ℓ 1 , ℓ 2 , dropout) and the classical baseline, partial least squares (PLS). Combining ƒ-divergence regularization with standard regularization yields further performance gains, indicating that distributional regularization is useful in this context giving a promising direction for robust model training in planetary science applications. Source code is publicly available at Klein and Li (2025), https://doi.org/10.11578/dc.20250530.7.

58 GEOSCIENCES↗

Stationary Distribution and Perturbation Bounds for a Stochastic Inventory Model

In this paper, we use the theory of generalized inverses to compute the stationary distribution of the (s, S) inventory model directly from the transition matrix of the underlying Markov chain and to provide perturbation bounds. Indeed, approximations are employed to build a tractable model, and statistical methods are used to estimate the unknown parameters and distributions. Hence, the system is subject to perturbations that may cause deviation in the characteristics. The proposed perturbation bound provides a means to estimate the impact of the perturbations on the performance measures of the considered inventory system.

97 MATHEMATICS AND COMPUTING↗

A semiempirical model for heavy element depletion in the interstellar medium

An empirical relationship is derived which reproduces the observed depletion factors in interstellar gas by relating them to some preshock distribution factor and the binding energy of the elements in appropriate compounds. By fitting this equation to the data, estimates are obtained for both the distribution factor and the critical binding energy, which is related to the shock strength at which grains are completely destroyed. The critical binding energy value therefore provides an estimate of the energy involved in the most recent shock event, while the preshock distribution factor is related to the preshock elemental distribution between gas and dust. It is noted that the minerals most likely to survive are the refractory ones with the largest binding energies.

Trivedi, B. M. P.↗

Experimental determination of the yield of “scission” neutrons from the spontaneous fission of {sup 252}Cf

Angular and energy distributions of prompt neutrons from the spontaneous fission of {sup 252}Cf have been measured. The yield of “scission” neutrons has been estimated by comparing the measured distributions with calculations within the model of emission of neutrons from completely accelerated fragments. It has been shown that the best description of measured angular and energy distributions of fission neutrons is achieved under the assumption that neutrons with a higher (about 6–8%) probability are emitted along the fission axis in the center-of-mass system of fission fragments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Graph-Learning-Assisted State and Event Tracking for Solar-Penetrated Power Grids with Heterogeneous Data Sources

Unlike transmission systems, distribution systems do not typically contain sufficient metering to enable real-time state estimation. The lack of sufficient real-time measurements prohibits accurate and timely monitoring of the state of distribution systems. As a result, control and optimal operation of distribution systems, especially those containing large numbers of renewable generation units are not possible without proper data and information about the current state of the system. The main motivation of this project is to address this shortcoming by developing an approach which provides “predicted” real-time measurements so that they can be used to execute a distribution system state estimator. Thus, the objective of the project is to make the distribution systems fully observable, such that the hosting capacity for solar generation can be accurately estimated, and unnecessary solar curtailments can be avoided. In order to accomplish this goal, the project investigated the use of a grid-model-informed machine learning (ML) tool which integrates heterogeneous data streams obtained from AMI meters, SCADA as well as PMU measurements and created synchronous measurement snapshots for the state estimator (SE); and developed a hybrid robust SE which provides not only accurate state estimates but also real-time feedback for the ML model refinement.

14 SOLAR ENERGY↗

Role of the Reference Bus in Three-Phase State Estimation

As more and more renewable energy sources are being connected to the power grid at lower voltage levels, closer monitoring of distribution networks is becoming a necessity. This is accomplished by state estimation and deployment of Distribution Management Systems (DMS). Distribution feeders are typically untransposed with unbalance loads. Thus, their modeling and analysis show significant differences compared to those for transmission networks. This paper focuses on the state estimation problem and its formulation for unbalanced distribution systems. Specifically, it is concerned about the modeling and solution of three-phase networks without a fully balanced bus to be selected as the angle reference. It is shown that best results will be obtained via the use of a virtual symmetrical bus. Validation of the proposed formulation is carried out by simulations.

Langner, A. L.↗

Observing System Simulations for Small Satellite Formations Estimating Bidirectional Reflectance

The bidirectional reflectance distribution function (BRDF) gives the reflectance of a target as a function of illumination geometry and viewing geometry, hence carries information about the anisotropy of the surface. BRDF is needed in remote sensing for the correction of view and illumination angle effects (for example in image standardization and mosaicing), for deriving albedo, for land cover classification, for cloud detection, for atmospheric correction, and other applications. However, current spaceborne instruments provide sparse angular sampling of BRDF and airborne instruments are limited in the spatial and temporal coverage. To fill the gaps in angular coverage within spatial, spectral and temporal requirements, we propose a new measurement technique: Use of small satellites in formation flight, each satellite with a VNIR (visible and near infrared) imaging spectrometer, to make multi-spectral, near-simultaneous measurements of every ground spot in the swath at multiple angles. This paper describes an observing system simulation experiment (OSSE) to evaluate the proposed concept and select the optimal formation architecture that minimizes BRDF uncertainties. The variables of the OSSE are identified; number of satellites, measurement spread in the view zenith and relative azimuth with respect to solar plane, solar zenith angle, BRDF models and wavelength of reflection. Analyzing the sensitivity of BRDF estimation errors to the variables allow simplification of the OSSE, to enable its use to rapidly evaluate formation architectures. A 6-satellite formation is shown to produce lower BRDF estimation errors, purely in terms of angular sampling as evaluated by the OSSE, than a single spacecraft with 9 forward-aft sensors. We demonstrate the ability to use OSSEs to design small satellite formations as complements to flagship mission data. The formations can fill angular sampling gaps and enable better BRDF products than currently possible.

small satellite↗

Determination of the area and mass distribution of orbital debris fragments

A technique is described to estimate the area-to-mass ratio of debris fragments using orbital fragments obtained by radar. The area-to-mass ratio of about 2600 fragments arising from the breakup of 24 artificial satellites was determined; an analysis of the data on about 200 objects with known mass, size, and shape has been made, and a calibration of the observed radar cross-section (RCS) to the effective area of these objects has provided a method to estimate the effective area of debris fragments. From the knowledge of the effective area and the estimated area-to-mass ratio, the mass and area distribution of each of the known breakup has been obtained. As a function of time, the orbital elements can be used to invert any propagation algorithm to yield the area-to-mass ratio of an orbiting object.

Badhwar, Gautam D.↗