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,711 records · Page 95

A global model for rain attenuation prediction

A model was developed for the prediction of rain attenuation statistics for 8 rain rate climate regions spanning all land regions from the equator to the poles. The model uses empirically determined rain rate distributions for each climate region and an empirical model to estimate path averaged rain rate for given values of point rain rate. Attenuation distribution observations were compared with model predictions. Good agreement was obtained for earth-space paths at 0.1 percent of the year. At the extremes, 0.01 and 1 percent of the year the observations were generally within the model estimate bounds but the agreement is not as good as at 0.1 percent of the year. The discrepancies are attributed to the limited duration of the attenuation observations and to the limited number of rain rate distribution observations.

Crane, R. K.↗

Rare Events via Cross-Entropy Population Monte Carlo

Rare events are events that happen with very low frequency. Estimating rare event probabilities using Monte Carlo techniques is computationally expensive, often to the point of intractability, and special methods are required. Importance sampling (IS) is a well known technique that uses a proposal distribution in place of a target distribution to lower the variance of the estimator. Key to the success of IS methods is the choice of a proposal distribution, or the parameters governing the distribution. Adaptive importance sampling improves the parameters of a family or population of proposal distributions iteratively through trials. We present a novel cross-entropy population Monte Carlo algorithm, which adapts the parameters of proposals through the cross-entropy method. The proposed method stands apart from previous work in that we are not optimizing a mixture distribution. Instead, we leverage deterministic mixture weights and optimize the distributions individually through a reinterpretation of the typical derivation of the cross-entropy method. Demonstrations on rare event examples show that the algorithm can outperform existing resampling based population Monte Carlo methods, especially for higher-dimensional problems. Finally, we also demonstrate efficacy on a conjunction analysis problem.

97 MATHEMATICS AND COMPUTING↗

Characterizing the System Impulse Response Function from Photon-Counting LiDAR Data

NASA's Multiple Altimeter Beam Experimental LiDAR (MABEL) is an aircraft-based photon-counting laser altimeter designed as a simulator to test measurement techniques and algorithms for Advanced Topographic Laser Altimeter System (ATLAS), the sole instrument on NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission. By measuring the time of flight, pointing angle, and absolute position for individual photons, ICESat-2 provides detailed elevation measurements of earth's surface. Calculating accurate and precise elevations requires an understanding of how photons interact with surfaces, and characterization of the photon distribution after returning from surfaces. Neither MABEL nor ATLAS records the transmitted laser pulse shape, relying instead on aggregating several pulses worth of photons, often using histograms, to characterize the pulse shape. In this paper, we assess the limitations of using histograms and propose a more robust method to describe MABEL's system impulse-response function using an exponentially modified Gaussian distribution. We also provide standard error estimates for the arithmetic mean and standard deviation calculations, and for exponentially modified Gaussian parameters using a Monte Carlo sensitivity analysis. We apply this method to photon returns from a sea ice lead and from a dry salt lake bed as case studies for estimating the standard error associated with sample size for the arithmetic mean and standard deviation, and for the exponentially modified Gaussian parameters. We use these standard errors to calculate the minimum number of photons required to find both Gaussian and exponentially modified Gaussian distribution parameters within 3 cm of their parent population values.

photoncounting↗

Temporal Fluctuations in a 100-A LaB6 Hollow Cathode

The temporal fluctuations in the near plume of a 100-Amperes LaB6 (Lanthanum hexaboride) hollow cathode are experimentally investigated. At high currents, turbulent oscillations may contribute to two of the anomalous processes of hollow cathode operation - anomalous resistivity and the production of energetic ions. A detailed study of the properties of the oscillations in a high current cathode is necessary to determine the impact of the fluctuation spectrum on these two effects. In this investigation, a probe array is employed to measure the amplitude and dispersion of axial modes in the plume while a retarding potential analyzer yields estimates of the radial ion energy distribution. The onset of the ion acoustic turbulence (IAT) is observed at high current values, and the character of the turbulent spectrum is shown to agree with weak turbulent theory: the amplitude of the spectrum decreases with flow rate but increases with discharge current. Estimates of the anomalous collision frequency based on experimental observations indicate that the IAT collision frequency can exceed the classical collision frequency at sufficiently high discharge current densities. Additionally, the onset of the IAT is shown to be correlated with the appearance of a high energy ion tail and that the energy in this tail is comparable to the energy in the experimentally-observed IAT.

production of ions↗

Characteristics of the vertical profiles of dual-frequency, dual-polarization radar data in stratiform rain

Airborne dual-wavelength and dual-polarization radar data are analyzed for measurements taken in stratiform rain in the western Pacific during September 1990. The focus of the paper is on the vertical profiles of the linear depolarization ratio, LDR (10 GHz); the reflectivity factor, dBZ (10 GHz); and the dual-frequency ratio, DFR (10, 34.45 GHz). Statistical characterizations of the maxima of these quantities and the relative locations at which they occur suggest that the eccentricity of the melting particles is fairly large and that the shape and size of the particles are correlated. To try to explain these features, two types of simulation are presented. In the first, a set of measured drop size distributions is used in the context of a standard model of the melting layer. Variations in snow density, as well as shape, size, and orientation distributions are used to study the relationship between these parameters and the radar measurements. To reduce the amount of ambiguity in the estimation, a second type of simulation is described in which the size distribution of the snow is estimated. Comparisons between the simulated and measured profiles indicate that radar measurements can be used to derive certain characteristics of the particle size and shape distributions in the melting layer.

Meneghini, R.↗

The Chandra cosmos legacy survey: clustering of X-ray-selected AGNs at 2.9 ≤ z ≤ 5.5 using photometric redshift probability distribution functions

We present the measurement of the projected and redshift-space two-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy active galactic nucleus (AGN) at 2.9 ≤ z ≤ 5.5 (〈L{sub bol}〉 ∼ 10{sup 46} erg s{sup −1}) using the generalized clustering estimator based on phot-z probability distribution functions in addition to any available spec-z. We model the projected 2pcf, estimated using π {sub max} = 200 h{sup −1} Mpc with the two-halo term and we derive a bias at z ∼ 3.4 equal to b = 6.6{sub −0.55}{sup +0.60}, which corresponds to a typical mass of the hosting halos of log M {sub h} = 12.83{sub −0.11}{sup +0.12} h{sup −1} M {sub ⊙}. A similar bias is derived using the redshift-space 2pcf, modeled including the typical phot-z error σ {sub z} = 0.052 of our sample at z ≥ 2.9. Once we integrate the projected 2pcf up to π {sub max} = 200 h{sup −1} Mpc, the bias of XMM and Chandra COSMOS at z = 2.8 used in Allevato et al. is consistent with our results at higher redshifts. The results suggest only a slight increase of the bias factor of COSMOS AGNs at z ≳ 3 with the typical hosting halo mass of moderate-luminosity AGNs almost constant with redshift and equal to log M {sub h} = 12.92{sub −0.18}{sup +0.13} at z = 2.8 and log M {sub h} = 12.83{sub −0.11}{sup +0.12} at z ∼ 3.4, respectively. The observed redshift evolution of the bias of COSMOS AGNs implies that moderate-luminosity AGNs still inhabit group-sized halos at z ≳ 3, but slightly less massive than observed in different independent studies using X-ray AGNs at z ≤ 2.

79 ASTRONOMY AND ASTROPHYSICS↗

Satellite-Derived Distributions, Inventories and Fluxes of Dissolved and Particulate Organic Matter Along the Northeastern U.S. Continental Margin

Estuaries and the coastal ocean experience a high degree of variability in the composition and concentration of particulate and dissolved organic matter (DOM) as a consequence of riverine and estuarine fluxes of terrigenous DOM, sediments, detritus and nutrients into coastal waters and associated phytoplankton blooms. Our approach integrates biogeochemical measurements, optical properties and remote sensing to examine the distributions and inventories of organic carbon in the U.S. Middle Atlantic Bight and Gulf of Maine. Algorithms developed to retrieve colored DOM (CDOM), Dissolved (DOC) and Particulate Organic Carbon (POC) from NASA's MODIS-Aqua and SeaWiFS satellite sensors are applied to quantify the distributions and inventories of DOC and POC. Horizontal fluxes of DOC and POC from the continental margin to the open ocean are estimated from SeaWiFS and MODIS-Aqua distributions of DOC and POC and horizontal divergence fluxes obtained from the Northeastern North Atlantic ROMS model. SeaWiFS and MODIS imagery reveal the importance of estuarine outflow to the export of CDOM and DOC to the coastal ocean and a net community production of DOC on the shelf.

Mannino, A.↗

Clusterdv: a simple density-based clustering method that is robust, general and automatic

Abstract Motivation How to partition a dataset into a set of distinct clusters is a ubiquitous and challenging problem. The fact that data vary widely in features such as cluster shape, cluster number, density distribution, background noise, outliers and degree of overlap, makes it difficult to find a single algorithm that can be broadly applied. One recent method, clusterdp, based on search of density peaks, can be applied successfully to cluster many kinds of data, but it is not fully automatic, and fails on some simple data distributions. Results We propose an alternative approach, clusterdv, which estimates density dips between points, and allows robust determination of cluster number and distribution across a wide range of data, without any manual parameter adjustment. We show that this method is able to solve a range of synthetic and experimental datasets, where the underlying structure is known, and identifies consistent and meaningful clusters in new behavioral data. Availability and implementation The clusterdv is implemented in Matlab. Its source code, together with example datasets are available on: https://github.com/jcbmarques/clusterdv. Supplementary information Supplementary data are available at Bioinformatics online.

Marques, João C.↗

Observations of Ocean Primary Productivity Using MODIS

Measuring the magnitude and variability of oceanic net primary productivity (NPP) represents a key advancement toward our understanding of the dynamics of marine ecosystems and the role of the ocean in the global carbon cycle. MODIS observations make two new contributions in addition to continuing the bio-optical time series begun with Orbview-2's SeaWiFS sensor. First, MODIS provides weekly estimates of global ocean net primary productivity on weekly and annual time periods, and annual empirical estimates of carbon export production. Second, MODIS provides additional insight into the spatial and temporal variations in photosynthetic efficiency through the direct measurements of solar-stimulated chlorophyll fluorescence. The two different weekly productivity indexes (first developed by Behrenfeld & Falkowski and by Yoder, Ryan and Howard) are used to derive daily productivity as a function of chlorophyll biomass, incident daily surface irradiance, temperature, euphotic depth, and mixed layer depth. Comparisons between these two estimates using both SeaWiFS and MODIS data show significant model differences in spatial distribution after allowance for the different integration depths. Both estimates are strongly dependence on the accuracy of the chlorophyll determination. In addition, an empirical approach is taken on annual scales to estimate global NPP and export production. Estimates of solar stimulated fluorescence efficiency from chlorophyll have been shown to be inversely related to photosynthetic efficiency by Abbott and co-workers. MODIS provides the first global estimates of oceanic chlorophyll fluorescence, providing an important proof of concept. MODIS observations are revealing spatial patterns of fluorescence efficiency which show expected variations with phytoplankton photo-physiological parameters as measured during in-situ surveys. This has opened the way for research into utilizing this information to improve our understanding of oceanic NPP variability. Deriving the ocean bio-optical properties places severe demands on instrument performance (especially band to band precision) and atmospheric correction. Improvements in MODIS instrument characterization and calibration over the first 16 mission months have greatly improved the accuracy of the chlorophyll input fields and FLH, and therefore the estimates of NPP and fluorescence efficiency. Annual estimates now show the oceanic NPP accounts for 40-50% of the global total NPP, with significant interannual variations related to large scale ocean processes. Spatial variations in ocean NPP, and exported production, have significant effects on exchange of CO2 between the ocean and atmosphere. Further work is underway to improve both the primary productivity model functions, and to refine our understanding of the relationships between fluorescence efficiency and NPP estimates. We expect that the MODIS instruments will prove extremely useful in assessing the time dependencies of oceanic carbon uptake and effects of iron enrichment, within the global carbon cycle.

Esaias, Wayne E.↗

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗

Systematic comparison of local approaches for isotopically nonstationary metabolic flux analysis

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determine cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA) that relies on the patterns of isotope labeling of metabolites in the network. For metabolic systems, typical for plants, where all potentially labeled atoms effectively have only one source atom pool, only isotopically nonstationary MFA can provide information about intracellular fluxes. There are several global approaches that implement MFA for an entire metabolic network and estimate, at once, a steady-state flux distribution for all reactions with identifiable fluxes in the network. In contrast, local approaches deal with estimation of fluxes for a subset of reactions, with smaller data demand for flux estimation. Here we present a systematic comparative review and benchmarking of the existing local approaches for isotopically nonstationary MFA. The comparison is conducted with respect to the required data and underlying computational problems solved on a synthetic network example. Furthermore, we benchmark the performance of these approaches in estimating fluxes for a subset of reactions using data obtained from the simulation of nitrogen fluxes in the Arabidopsis thaliana core metabolism. The findings pinpoint practical aspects that need to be considered when applying local approaches for flux estimation in large-scale plant metabolic networks.

59 BASIC BIOLOGICAL SCIENCES↗

Radar-based Snow Property Profile Retrievals for NSA

The product contains estimated snow microphysical properties retrieved primarily from vertically-profiling radar observations using a priori constraints that are dependent on atmospheric state. The product provides data at 1-minute time resolution, including snowfall rate at the surface and vertical profiles at 300-meter resolution of snow particle size distribution parameters and snow water contents. Estimated uncertainties for the snowfall rates and snow water contents are also provided. The product is obtained using an optimal estimation retrieval in which radar forward model and a priori state characteristics are based on expected snow particle microphysical and radar scattering properties evaluated from field experiment data.

54 ENVIRONMENTAL SCIENCES↗

Neck Rupture and Scission Neutrons in Nuclear Fission

Just before a nucleus undergoes fission, a neck is formed between the emerging fission fragments. It is widely accepted that this neck undergoes a rather violent rupture, despite the absence of unambiguous experimental evidence. The main difficulty in addressing the neck rupture and saddle-to-scission stages of fission is that both are highly nonequilibrium processes. Here, we present the first fully microscopic characterization of the scission mechanism, along with the spectrum and the spatial distribution of scission neutrons (SNs), and some upper limit estimates for the emission of charged particles. The spectrum of SNs has a distinct angular distribution, with neutrons emitted in roughly equal numbers in the equatorial plane and along the fission axis. They carry an average energy around 3±0.5 MeV for the fission of 236 U, 240 Pu, and 252 Cf, and a maximum of 16–18 MeV. We estimate a conservative lower bound of 9%–14% of the total emitted neutrons are produced at scission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutrino Tridents at Near Detectors [Slides]

A full 4PS computation is required to obtain correct estimates for the number of events. Let us stress that the EPA gives a reasonable result for the dimuon channel due to a serendipitous behavior of the Form Factor. EPA can artificially suppress the coherent scattering contribution and increase the diffractive one giving rise to an incorrect rate and distributions of observable quantities. We have estimated the background for each trident channel via a Monte Carlo simulation using GENIE, and identified the dominant contributions arising primarily from particle misidentification. Reduction of ~6 orders of magnitude x CC in the background is necessary to observe trident events at DUNE ND. A more careful analysis is needed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-approach assessment of the spatial distribution of the specific yield: application to the Crau plain aquifer, France

Spatially distributed values of the specific yield, a fundamental parameter for transient groundwater mass balance calculations, were obtained by means of three independent methods for the Crau plain, France. In contrast to its traditional use to assess recharge based on a given specific yield, the water-table fluctuation (WTF) method, applied using major recharging events, gave a first set of reference values. Then, large infiltration processes recorded by monitored boreholes and caused by major precipitation events were interpreted in terms of specific yield by means of a one-dimensional vertical numerical model solving Richards’ equations within the unsaturated zone. Finally, two gravity field campaigns, at low and high piezometric levels, were carried out to assess the groundwater mass variation and thus alternative specific yield values. The range obtained by the WTF method for this aquifer made of alluvial detrital material was 2.9– 26%, in line with the scarce data available so far. The average spatial value of specific yield by the WTF method (9.1%) is consistent with the aquifer scale value from the hydro-gravimetric approach. In this investigation, an estimate of the hitherto unknown spatial distribution of the specific yield over the Crau plain was obtained using the most reliable method (the WTF method). A groundwater mass balance calculation over the domain using this distribution yielded similar results to an independent quantification based on a stable isotope-mixing model. This agreement reinforces the relevance of such estimates, which can be used to build a more accurate transient hydrogeological model.

58 GEOSCIENCES↗

A General State Estimation Formulation for Three-Phase Unbalanced Power Systems

Almost all of the three-phase state estimation algorithms assume existence of a reference bus whose phase angles are perfectly balanced. This assumption is quite realistic for transmission systems, and also for most distribution systems that are connected to a strong transmission system where transmission side can modeled by a balanced reference bus. However, for distribution systems having high penetration of renewable sources or for microgrids operating in islanded mode, the assumption of a balanced reference bus will not be realistic. While there are recent publications focusing on this problem, formulation of the three-phase unbalanced state estimation problem with proper treatment of the reference bus remains unaddressed. In this paper, a new formulation will be described where an accurate state estimation solution can be obtained for any unbalanced threephase system irrespective of its operating conditions (balanced or highly unbalanced), configuration (isolated microgrid, connected to transmission system, etc.) and whether or not it contains any synchronous generators. Validation of the proposed formulation will be carried out via simulations.

State Estimation, Distribution System, Unbalanced ↗

Microscopic theory of angular momentum distributions across the full range of fission fragments

Modern nuclear theory provides qualitative insights into the fundamental mechanisms of nuclear fission and is increasingly capable of making reliable quantitative predictions. Most quantities of interest pertain to the primary fission fragments, whose subsequent decay is typically modeled using statistical reaction models. Consequently, a key objective of fission theory is to inform these models by predicting the initial conditions of the primary fragments. Here, in this work, we employ a framework that combines joint angular momentum and particle number projection with time-dependent configuration mixing to calculate the angular momentum distributions of primary fragments. Focusing on the benchmark cases of neutron-induced fission of 235 U and 239 Pu , we predict—for the first time—microscopic angular momentum distributions for all fragments observed in experiments. Our results reveal a pronounced sawtooth pattern in the average angular momentum as a function of fragment mass, consistent with recent measurements. Additionally, we observe substantial variations in angular momentum distributions along isobaric chains, indicating that commonly used empirical formulas lack sufficient accuracy. We also quantify a strong correlation between the angular momentum and the deformation of the fragments at scission, and a weak correlation in the magnitude of the angular momentum between fragment partners. The generated data will enable estimation of the impact of microscopic distributions on fission spectra, paving the way toward fission modeling based on microscopic inputs.

Physics - Nuclear physics and radiation physics↗

On the implementation of modal filters for control of structures

The most common technique for the control of structures is modal control. In modal control, the differential equations in terms of actual coordinates are replaced by a set of ordinary differential equations in terms of the modal coordinates known as modal equations. In designing feedback controls in conjunction with the modal equations, one must know the modal states for the modes targeted for control. The sensors measure actual states, however. The modal states can be estimated by means of a Luenberger observer or modal filters. The modal filters produce estimates of the modal states from distributed measurements of the states. If distributed measurements are not available, then they can be reconstructed from measurements at discrete points via interpolation. This paper examines various questions associated with the implementation of modal filters, such as the effect of choice of interpolation functions and sensors locations, as well as of measurement errors, on the state estimation process. The method is demonstrated by means of two numerical examples.

Meirovitch, L.↗