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 793 records · Page 44

Small impact craters in the lunar regolith - Their morphologies, relative ages, and rates of formation

The size frequency distributions, relative ages, and absolute age determination of lunar impact craters ranging in diameter from 1 to 100 m are discussed in relation to present-day and past impact rates. The size frequency distributions of craters produced by impacts on the lunar surface and of craters remaining after degradation with time and later impacts are examined, together with the steady state frequency curves for various crater morphologies. Relations between the relative ages of craters, their morphologies and their diameters are discussed, and the calibration of this kind of relationship with absolute crater ages by the location of individual craters with known absolute ages, diameters and frequency distribution statistics is considered. Estimates of present and past cratering rates based on crater size-frequency distributions are presented which suggest that the rate of impacting objects producing small lunar craters could have been constant over the last 100 million years.

Moore, H. J.↗

Quantifying Uncertainty in Particle Size Distribution Parameters Inferred from SAGE III/ISS Extinction Spectra

Stratospheric aerosols play key roles in the chemistry and radiation balance of the atmosphere and are a key input parameter for global chemistry and climate models. The degree to which aerosols impact chemistry and radiation balance depends primarily on their microphysical properties such as particle size distribution (PSD). The PSD is a mathematical description of the relative abundance of different sized particles within a sampling volume. If the PSD is accurately known then other key modeling parameters (e.g., surface area density) can be derived. Occultation observations from orbital instruments such as SAGE III/ISS have been used to infer these PSD parameters by inverting the extinction coefficient spectra. However, past efforts failed to address two key issues with this methodology: 1. The measurement uncertainty was not accounted for; 2. They assumed the PSDs to be single-mode only, while “real-world” PSDs are predominantly bi-modal. Accounting for both issues in the retrieval will yield an expanded solution space to the inferred PSD parameters; the question is “by how much?” To address this knowledge gap, we propose to carry out a series of simulations and, for every valid SAGE III/ISS data point, determine the range of PSD parameters that yield extinction spectra that are indistinguishable from the SAGE III/ISS data, within the limits of the reported uncertainty. Further, we will expand the solution space, for the first time, to include bimodal distributions. The results of this work will advance Earth system modeling/prediction capability through identifying the uncertainty of PSD parameter estimates using SAGE III/ISS data. The key benefits of this study over previous studies are twofold: 1. we will provide PSD estimates that include bimodal distributions in the solution space, 2. we will provide an uncertainty estimate for these parameters. The results of this study may be used directly in current and future climate and chemistry models.

SAGE III/ISS↗

Delensing the CMB with the cosmic infrared background: the impact of foregrounds

ABSTRACT The most promising avenue for detecting primordial gravitational waves from cosmic inflation is through measurements of degree-scale cosmic microwave background (CMB) B-mode polarization. This approach must face the challenge posed by gravitational lensing of the CMB, which obscures the signal of interest. Fortunately, the lensing effects can be partially removed by combining high-resolution E-mode measurements with an estimate of the projected matter distribution. For near-future experiments, the best estimate of the latter will arise from co-adding internal reconstructions (derived from the CMB itself) with external tracers such as the cosmic infrared background (CIB). In this work, we characterize how foregrounds impact the delensing procedure when CIB intensity, I, is used as the matter tracer. We find that higher point functions of the CIB and Galactic dust such as 〈BEI〉c and 〈EIEI〉c can, in principle, bias the power spectrum of delensed B-modes. To quantify these, we first estimate the dust residuals in currently available CIB maps and upcoming, foreground-cleaned Simons Observatory CMB data. Then, using non-Gaussian simulations of Galactic dust – extrapolated to the relevant frequencies, assuming the spectral index of polarized dust emission to be fixed at the value determined by Planck – we show that the bias to any primordial signal is small compared to statistical errors for ground-based experiments, but might be significant for space-based experiments probing very large angular scales. However, mitigation techniques based on multifrequency cleaning appear to be very effective. We also show, by means of an analytical model, that the bias arising from the higher point functions of the CIB itself ought to be negligible.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting Source Convergence Diagnostics in the KENO Monte Carlo Neutron Transport Codes [Abstract]

Monte Carlo criticality transport codes, which rely on the power iteration procedure, are a fundamental tool for nuclear criticality safety practitioners in assessing the neutron multiplication factor (k eff ) for problems involving fissile material. In these calculations, ensuring the convergence of both the fission source distributions and the k eff estimate for accurate results is crucial. However, a converged k eff estimate does not necessarily mean the fission source distribution is also converged because the fission source and flux distribution may continue to evolve even after k eff convergence. Therefore, most Monte Carlo transport criticality codes now offer various diagnostic tests to assess fission source convergence in addition to the k eff convergence by analyzing the trends of these quantities over multiple generations.

AZURE↗

Constraining the Milky Way Mass Profile with Phase-space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M={1.23}_{-0.18}^{+0.21}\times {10}^{12}{M}_{\odot }$ and a concentration of $c={9.4}_{-2.1}^{+2.8}$ with the prior based on the M–c relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of the kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. Finally, by varying the disruption rate of such satellites in this simulation, we estimate a ~5% scatter in the inferred MW halo mass among hydrodynamics-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Fusion techniques using distributed Kalman filtering for detecting changes in systems

A comparison is made of the performances of two detection strategies that are based on different data fusion techniques. The strategies detect changes in a linear system. One detection strategy involves combining the estimates and error covariance matrices of distributed Kalman filters, generating a residual from the used estimates, comparing this residual to a threshold, and making a decision. The other detection strategy involves a distributed decision process in which estimates from distributed Kalman filters are used to generate distributed residuals which are compared locally to a threshold. Local decisions are made and these decisions are then fused into a global decision. The performances of each of these detection schemes are compared, and it is concluded that better performance is achieved when local decisions are made and then fused into a global decision.

Belcastro, Celeste M.↗

LACIE performance predictor final operational capability program description, volume 3

The requirements and processing logic for the LACIE Error Model program (LEM) are described. This program is an integral part of the Large Area Crop Inventory Experiment (LACIE) system. LEM is that portion of the LPP (LACIE Performance Predictor) which simulates the sample segment classification, strata yield estimation, and production aggregation. LEM controls repetitive Monte Carlo trials based on input error distributions to obtain statistical estimates of the wheat area, yield, and production at different levels of aggregation. LEM interfaces with the rest of the LPP through a set of data files.

Source record↗

The terrestrial cratering record. II - The crater production rate

The terrestrial crater production rate is estimated from a compilation of the locations and ages of proven and probable terrestrial impact craters and compared with that calculated from astronomical observations and for other planets. The estimated size-frequency distribution of craters greater than 20 km in diameter formed during the Phanerozoic, a -2 power law, is similar to the distributions of lunar and Martian craters. The North American and East European cratons are estimated to have cratering rates of craters greater than 20 km of 0.36 x 10 to the -14/sq km yr and 0.33 x 10 to the -14/sq km yr, respectively, in good agreement with previous estimates and observations of Apollo bodies. Comparisons with the lunar rate indicate that the lunar and terrestrial rates overlap if the cratering flux had been constant over the last 3.4 billion years, however if a constant rate was not reached until 3.0 to 2.5 billion years ago, the lunar rate would have been half that of the terrestrial rate.

Grieve, R. A. F.↗

A validation methodology for fault-tolerant clock synchronization

A validation method for the synchronization subsystem of a fault-tolerant computer system is presented. The high reliability requirement of flight crucial systems precludes the use of most traditional validation methods. The method presented utilizes formal design proof to uncover design and coding errors and experimentation to validate the assumptions of the design proof. The experimental method is described and illustrated by validating an experimental implementation of the Software Implemented Fault Tolerance (SIFT) clock synchronization algorithm. The design proof of the algorithm defines the maximum skew between any two nonfaulty clocks in the system in terms of theoretical upper bounds on certain system parameters. The quantile to which each parameter must be estimated is determined by a combinatorial analysis of the system reliability. The parameters are measured by direct and indirect means, and upper bounds are estimated. A nonparametric method based on an asymptotic property of the tail of a distribution is used to estimate the upper bound of a critical system parameter. Although the proof process is very costly, it is extremely valuable when validating the crucial synchronization subsystem.

Johnson, S. C.↗

On adaptive weighted polynomial preconditioning for Hermitian positive definite matrices

The conjugate gradient algorithm for solving Hermitian positive definite linear systems is usually combined with preconditioning in order to speed up convergence. In recent years, there has been a revival of polynomial preconditioning, motivated by the attractive features of the method on modern architectures. Standard techniques for choosing the preconditioning polynomial are based only on bounds for the extreme eigenvalues. Here a different approach is proposed, which aims at adapting the preconditioner to the eigenvalue distribution of the coefficient matrix. The technique is based on the observation that good estimates for the eigenvalue distribution can be derived after only a few steps of the Lanczos process. This information is then used to construct a weight function for a suitable Chebyshev approximation problem. The solution of this problem yields the polynomial preconditioner. In particular, we investigate the use of Bernstein-Szego weights.

Fischer, Bernd↗

Simulation of radar reflectivity and surface measurements of rainfall

Raindrop size distributions (RSDs) are often estimated using surface raindrop sampling devices (e.g., disdrometers) or optical array (2D-PMS) probes. A number of authors have used these measured distributions to compute certain higher-order RSD moments that correspond to radar reflectivity, attenuation, optical extinction, etc. Scatter plots of these RSD moments versus disdrometer-measured rainrates are then used to deduce physical relationships between radar reflectivity, attenuation, etc., which are measured by independent instruments (e.g., radar), and rainrate. In this paper RSDs of the gamma form as well as radar reflectivity (via time series simulation) are simulated to study the correlation structure of radar estimates versus rainrate as opposed to RSD moment estimates versus rainrate. The parameters N0, D0 and m of a gamma distribution are varied over the range normally found in rainfall, as well as varying the device sampling volume. The simulations are used to explain some possible features related to discrepancies which can arise when radar rainfall measurements are compared with surface or aircraft-based sampling devices.

Chandrasekar, V.↗

Retrieval of Hydrometeor Drop Size Distributions from TRMM Field Campaign Profiler Doppler Velocity Spectra Observations

Consistent with the original proposal and work plan, this project focused on estimating the raindrop size distributions (DSDs) retrieved from vertically pointing Doppler radar profilers and analyzing the relationship of the retrieved DSDs with the dynamics of the precipitation processes. The first phase of this project focused on developing the model to retrieve the DSD from the observed Doppler velocity spectra. The second phase used this model to perform DSD retrievals from the profiler observations made during the TRMM Ground Validation Field Campaigns of TEFLUN-B, TRMM-LBA, and KWAJEX. The third phase of this project established collaborations with scientists involved with each field campaign in order to validate the profiler DSD estimates and to enable the profiler retrievals to be used in their research. Through these collaborations, the retrieved DSDs were placed into context with the dynamical processes of the observed precipitating cloud systems.

Williams, Christopher R.↗

Angular Distribution Models for Top-of-Atmosphere Radiative Flux Estimation from the Clouds and the Earth's Radiant Energy System Instrument on the Terra Satellite: Methodology - Part 1

The Clouds and Earth's Radiant Energy System (CERES) provides coincident global cloud and aerosol properties together with reflected solar, emitted terrestrial longwave and infrared window radiative fluxes. These data are needed to improve our understanding and modeling of the interaction between clouds, aerosols and radiation at the top of the atmosphere, surface, and within the atmosphere. This paper describes the approach used to estimate top-of-atmosphere (TOA) radiative fluxes from instantaneous CERES radiance measurements on the Terra satellite. A key component involves the development of empirical angular distribution models (ADMs) that account for the angular dependence of Earth's radiation field at the TOA. The CERES Terra ADMs are developed using 24 months of CERES radiances, coincident cloud and aerosol retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS), and meteorological parameters from the Global Modeling and Assimilation Office (GMA0) s Goddard Earth Observing System DAS (GEOS-DAS V4.0.3) product. Scene information for the ADMs is from MODIS retrievals and GEOS-DAS V4.0.3 properties over ocean, land, desert and snow, for both clear and cloudy conditions. Because the CERES Terra ADMs are global, and far more CERES data is available on Terra than was available from CERES on the Tropical Rainfall Measuring Mission (TRMM), the methodology used to define CERES Terra ADMs is different in many respects from that used to develop CERES TRMM ADMs, particularly over snow/sea-ice, under cloudy conditions, and for clear scenes over land and desert.

Loeb, N. G.↗

Predictive and Cooperative Voltage Control with Probabilistic Load and Solar Generation Forecasting

This paper proposes predictive cooperative voltage control method in a power system with high penetration of photovoltaic (PV) units. Cooperative distributed control of the reactive power output of PV inverters is coordinated with operation of voltage regulators (VRs) to maintain system voltages within an appropriate bandwidth. Probabilistic forecasting of the solar power generation and the loads is applied to estimate voltage changes which, in turn, are used to set the VR tap positions for preventing large voltage fluctuations with the lowest risk considering the voltage distribution estimation. The fine tuning of voltage adjustment is achieved by cooperative control of PV inverters to maintain a uniform voltage profile across the system. The proposed method is tested on the modified IEEE 123-node test feeder with high PV penetration using real insolation data and with constant loads replaced by several different load profiles. Simulation results demonstrate the effectiveness of the coordinated approach for voltage control with cooperative PV and predictive VR controls taking into account probabilistic load and solar power forecasts.

Cooperative Control↗

Simulation studies of wide and medium field of view earth radiation data analysis

A parameter estimation technique is presented to estimate the radiative flux distribution over the earth from radiometer measurements at satellite altitude. The technique analyzes measurements from a wide field of view (WFOV), horizon to horizon, nadir pointing sensor with a mathematical technique to derive the radiative flux estimates at the top of the atmosphere for resolution elements smaller than the sensor field of view. A computer simulation of the data analysis technique is presented for both earth-emitted and reflected radiation. Zonal resolutions are considered as well as the global integration of plane flux. An estimate of the equator-to-pole gradient is obtained from the zonal estimates. Sensitivity studies of the derived flux distribution to directional model errors are also presented. In addition to the WFOV results, medium field of view results are presented.

Green, R. N.↗

Nonlinear state estimation and feedback control of nonlinear and bilinear distributed parameter systems

This paper presents a theory of nonlinear state observers for nonlinear and bilinear distributed parameter systems. Convergence results are proved for these observers. Linear feedback control derived from such state observers is applied to the distributed parameter system and conditions are presented for closed-loop stability. The emphasis is on finite dimensional state observers and controllers (which can be implemented with on-line computers) and conditions for their successful operation with infinite dimensional distributed parameter systems.

Balas, M. J.↗

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Evaluation of Estimates Using Independent Data - Part 2

Rainfall rate estimates from space-borne k&ents are generally accepted as reliable by a majority of the atmospheric science commu&y. One-of the Tropical Rainfall Measuring Mission (TRh4M) facility rain rate algorithms is based upon passive microwave observations fiom the TRMM Microwave Imager (TMI). Part I of this study describes improvements in the TMI algorithm that are required to introduce cloud latent heating and drying as additional algorithm products. Here, estimates of surface rain rate, convective proportion, and latent heating are evaluated using independent ground-based estimates and satellite products. Instantaneous, OP5resolution estimates of surface rain rate over ocean fiom the improved TMI algorithm are well correlated with independent radar estimates (r approx. 0.88 over the Tropics), but bias reduction is the most significant improvement over forerunning algorithms. The bias reduction is attributed to the greater breadth of cloud-resolving model simulations that support the improved algorithm, and the more consistent and specific convective/stratiform rain separation method utilized. The bias of monthly, 2.5 deg. -resolution estimates is similarly reduced, with comparable correlations to radar estimates. Although the amount of independent latent heating data are limited, TMI estimated latent heating profiles compare favorably with instantaneous estimates based upon dual-Doppler radar observations, and time series of surface rain rate and heating profiles are generally consistent with those derived from rawinsonde analyses. Still, some biases in profile shape are evident, and these may be resolved with: (a) additional contextual information brought to the estimation problem, and/or; (b) physically-consistent and representative databases supporting the algorithm. A model of the random error in instantaneous, 0.5 deg-resolution rain rate estimates appears to be consistent with the levels of error determined from TMI comparisons to collocated radar. Error model modifications for non-raining situations will be required, however. Sampling error appears to represent only a fraction of the total error in monthly, 2S0-resolution TMI estimates; the remaining error is attributed to physical inconsistency or non-representativeness of cloud-resolving model simulated profiles supporting the algorithm.

Yang, Song↗

TLSLEAF: Automatic Leaf Angle Estimates From Single-Scan Terrestrial Laser Scanning

Leaf angle distribution (LAD) in forest canopies affects estimates of leaf area, light interception, and global-scale photosynthesis, but is often simplified to a single theoretical value. Here, we present TLSLeAF (Terrestrial Laser Scanning Leaf Angle Function), an automated open-source method of deriving LADs from terrestrial laser scanning. TLSLeAF produces canopy-scale leaf angle and LADs by relying on gridded laser scanning data. The approach increases processing speed, improves angle estimates, and requires minimal user input. Key features are automation, leaf–wood classification, beta parameter output, and implementation in R to increase accessibility for the ecology community. TLSLeAF precisely estimates leaf angle with minimal distance effects on angular estimates while rapidly producing LADs on a consumer-grade machine. We challenge the popular spherical LAD assumption, showing sensitivity to ecosystem type in plant area index and foliage profile estimates that translate to c. 25% and c. 11% increases in canopy net photosynthesis (c. 25%) and solar-induced chlorophyll fluorescence (c. 11%). TLSLeAF can now be applied to the vast catalog of laser scanning data already available from ecosystems around the globe. The ease of use will enable widespread adoption of the method outside of remote-sensing experts, allowing greater accessibility for addressing ecological hypotheses and large-scale ecosystem modeling efforts.

3D↗