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,423 records · Page 79

The rate of binary black hole mergers inferred from advanced LIGO observations surrounding GW150914

A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on 2015 September 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 days around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false-alarm rate (FAR) of <4.9×10{sup −6} yr{sup −1}, yielding a p-value for GW150914 of <2×10{sup −7}. Parameter estimation follow-up on this trigger identifies its source as a binary black hole (BBH) merger with component masses (m{sub 1},m{sub 2})=(36{sub −4}{sup +5},29{sub −4}{sup +4}) M{sub ⊙} at redshift z=0.09{sub −0.04}{sup +0.03} (median and 90% credible range). Here, we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between 2--53 Gpc{sup −3} yr{sup −1} (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from 13--600 Gpc{sup −3} yr{sup −1} depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range 2--600 Gpc{sup −3} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Geometric Tail Approximation for Reliability and Survivability

A common problem in developing high-reliability systems is estimating the reliability for a population of components that cannot be 100% tested. The radiation survivability of a population of components is often estimated by testing a very small sample to some multiple of the required specification level, known as an overtest. Given a successful test with a sufficient overtest margin, the population of components is assumed to have the required survivability or radiation reliability. However, no mathematical justification for such claims has been crafted without making aggressive assumptions regarding the statistics of the unknown distribution. Here we illustrate a new approach that leverages geometric bounding arguments founded on relatively modest distribution assumptions to produce conservative estimates of component reliability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Comparative Study on Computation of Cumulative Distribution Function in Predicting Time of Failure of Engineering Systems

Estimating accurate Time-of-Failure (ToF) of a system is key in making the decisions that impact operational safety and optimize cost. In this context, it is interesting to note that different approaches have been explored to tackle the problem of estimating ToF. The difference is in part characterized by different definitions of the hazard zones. The conventional definition for the cumulative distribution function (CDF) calculation is assumed to have well-defined hazard zones, that is, hazard zones defined as a function of the system state trajectory. An alternate method suggests the use of hazard zones defined as a function of the system state at time , instead of hazard zones defined as a function of system state up to and including time k (Acuna and Orchard 2018, 2017). This paper explores these differences and their impact on ToF estimation. Results for the conventional CDF definition indicated that, (i) the cumulative distribution function is always an increasing function of time, even when realizations of the degradation process are not monotonic, (ii) the sum of all probabilities is always 1 and does not need to be normalized, and (iii) all probabilities are positive and less than or equal to 1. Similar results are not observed for CDF calculation with hazard zones defined as a function only of the system state at time . Results for ToF estimation using Acuna's definition differ, suggesting that there is an underlying assumption of independence in the hazard zone definition. Therefore, we present an alternate definition of hazard zone which guarantees the properties of a well-defined CDF with a more straightforward ToF definition.

Hazard Zone↗

On the Path to SunShot. Emerging Issues and Challenges in Integrating Solar with the Distribution System

This report analyzes distribution-integration challenges, solutions, and research needs in the context of distributed generation from PV (DGPV) deployment to date and the much higher levels of deployment expected with achievement of the U.S. Department of Energy's SunShot targets. Recent analyses have improved estimates of the DGPV hosting capacities of distribution systems. This report uses these results to statistically estimate the minimum DGPV hosting capacity for the contiguous United States using traditional inverters of approximately 170 GW without distribution system modifications. This hosting capacity roughly doubles if advanced inverters are used to manage local voltage and additional minor, low-cost changes could further increase these levels substantially. Key to achieving these deployment levels at minimum cost is siting DGPV based on local hosting capacities, suggesting opportunities for regulatory, incentive, and interconnection innovation. Already, pre-computed hosting capacity is beginning to expedite DGPV interconnection requests and installations in select regions; however, realizing SunShot-scale deployment will require further improvements to DGPV interconnection processes, standards and codes, and compensation mechanisms so they embrace the contributions of DGPV to system-wide operations. SunShot-scale DGPV deployment will also require unprecedented coordination of the distribution and transmission systems. This includes harnessing DGPV's ability to relieve congestion and reduce system losses by generating closer to loads; minimizing system operating costs and reserve deployments through improved DGPV visibility; developing communication and control architectures that incorporate DGPV into system operations; providing frequency response, transient stability, and synthesized inertia with DGPV in the event of large-scale system disturbances; and potentially managing reactive power requirements due to large-scale deployment of advanced inverter functions. Finally, additional local and system-level value could be provided by integrating DGPV with energy storage and 'virtual storage,' which exploits improved management of electric vehicle charging, building energy systems, and other large loads. Together, continued innovation across this rich distribution landscape can enable the very-high deployment levels envisioned by SunShot.

14 SOLAR ENERGY↗

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions

Abstract This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three‐dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R ‐squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically‐informed estimates. The results of this work are used to inform a follow‐on study that tests additional data‐driven approaches in multiple basins within the contiguous United States.

54 ENVIRONMENTAL SCIENCES↗

The diffuse galactic gamma radiation - The Compton contribution and component separation by energy interval and galactic coordinates

The diffuse high-energy galactic gamma radiation to be expected from cosmic ray interactions with matter and photons is considered with particular emphasis on the contribution of Compton radiation from cosmic ray electrons. The intensity, spectrum and spatial distribution of the expected galactic gamma radiation are estimated based on models of the matter, cosmic ray and photon distributions to take into account the contributions of bremsstrahlung, high-energy cosmic-ray nucleon and interstellar matter interactions as well as Compton interactions between cosmic ray electrons and background photons. Results suggest that the Compton gamma ray contribution from cosmic ray electron interactions with galactic visible and infrared photons is substantially larger than previously believed. Analysis of the energy spectra and latitude dependence of the various sources reveals that the Compton radiation, bremsstrahlung and nuclear cosmic ray-matter interaction radiation should be separable, with Compton radiation dominating at energies from 10 to 100 MeV at galactic latitudes greater than several degrees. Results demonstrate the potential of gamma ray observations in studies of galactic structure, cosmic ray electrons and galactic photon density.

Kniffen, D. A.↗

A Hybrid-Learning Algorithm for Online Dynamic State Estimation in Multimachine Power Systems

With the increasing penetration of distributed generators in the smart grids, having knowledge of rapid real-time electromechanical dynamic states has become crucial to system stability control. Conventional Supervisory Control and Data Acquisition (SCADA)-based dynamic state estimation (DSE) techniques are limited by the slow sampling rates, while the emerging phasor measurement units (PMUs) technology enables rapid real-time measurements at network nodes. Using generator bus terminal voltages, we propose a hybrid-learning DSE (HL-DSE) algorithm to estimate the synchronous machine rotor angle and speed in real time. The HL-DSE takes the power system model into account and trains neuroestimators with real-time data in an online manner. Compared with traditional DSE methods, the HL-DSE overcomes limitations by using a data-driven approach in conjunction with the physical power system model. The time efficiency, accuracy, convergence, and robustness of the proposed algorithm are tested under noises and fault conditions in both small- and large-scale test systems. Simulation results show that the proposed HL-DSE is much more computationally efficient than widely used Kalman filter (KF)-based methods while maintaining comparable accuracy and robustness. In particular, HL-DSE is over 100 times faster than square-root unscented KF (SR-UKF) and 80 times faster than extended KF (EKF). The advantages and challenges of the HL-DSE are also discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Energy Release upon Bombardment of a Tungsten Surface by Deuterium Atoms

The depth distribution of the energy release (linear energy transfer) occurring upon bombardment of a tungsten target by deuterium atoms with energies ranging from 100 eV to 10 MeV is calculated. It is demonstrated that, contrary to common belief, the maximum energy release at energies below 100 keV is localized near the surface of the solid. At energies above 100 keV, the shape of the distribution changes, and a Bragg maximum emerges near the particle stopping point. The depth distribution of the energy release in tungsten under conditions typical of the ITER tokamak reactor is determined. This distribution provides an opportunity to estimate the wall heating by incident plasma atoms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Midinfrared spectral investigations of carbonates: Analysis of remotely sensed data

Recent airborne thermal infrared observations of Mars from the Kuiper Airborne Observatory (KAO) have provided evidence for the presence of carbonates, sulfates, and hydrates. Using the optical properties of calcite and anhydrite, it was estimated that CO3's and SO4's constituted about 1 to 3 and 10 to 15 wt. percent, repectively of the materials composing the atmospheric dust. Using the derived value as an estimate of total CO3 abundance, and making an assumption that the CO3's were uniformly distributed within the Martian regolith, it was estimated that such a CO3 reservoir could contain roughly 2 to 5 bars of CO2. While the results indicate that several volatile-bearing materials are present on Mars, the observations from the KAO are inherently limited in their ability to determine the spatial distributions of these materials. However, previous spacecraft observations of Mars provide both the spectral coverage necessary to identify these materials, as well as the potential for investigating their spatial variability. This has prompted us to pursue a reinvestigation of the Mariner 6 and 7 infrared spectrometer and Mariner 9 infrared interferometer spectrometer observations. The former data have been recently made available in digital format and calibration of wavelengths and intensities are almost complete. Additionally, we are pursuing the derivation of optical constants of more appropriate carbonates and sulfates.

Roush, T.↗

Effect of graphite morphology on the tensile strength and thermal conductivity of cast iron

Highlights: • The graphite morphology is characterized from 3D structure using X–ray tomography. • The size, distribution, shape and interconnectivity of graphite are estimated. • The influences of graphite morphology on tensile properties and thermal conductivity are evaluated. • The mechanism of tensile strength and the model of thermal conductivity are discussed. - Abstract: Graphite morphology is a very important parameter affecting the properties of cast iron. In this study, the graphite morphology was characterized using both 3D and 2D structures. The experimental results revealed that when the graphite was less spheroidal, the graphite had greater interconnectivity in the 3D structure of the cast iron. The size, distribution, volume fraction, and shape of the graphite were estimated. Moreover, the influences of the graphite's morphology, size, and distribution on the tensile property and thermal conductivity at different temperatures were estimated. The experimental results revealed that with an increase in temperature, the thermal conductivity initially increased, and then decreased, whereas the tensile strength decreased slightly at first, and then decreased rapidly. In addition, with an increase in the percentage of compacted graphite, the tensile strength decreased and the thermal conductivity increased. The fracture surfaces of tensile specimens were also estimated. Finally, the mechanism of tensile strength and the model of thermal conductivity were investigated.

36 MATERIALS SCIENCE↗

SU-F-T-159: Monte Carlo Simulation Studies of Three-Dimensional Dose Distribution for Polymer Gel Dosimeter and Radiochromic Gel Dosimeter in a Proton Beam

Purpose: The purpose of this simulation study is to evaluate the proton detectability of gel dosimeters, and estimate the three-dimensional dose distribution of protons in the radiochromic gel and polymer gel dosimeter compared with the dose distribution in water. Methods: The commercial composition ratios of normoxic polymer gel and LCV micelle radiochromic gel were included in this simulation study. The densities of polymer and radiochromic gel were 1.024 and 1.005 g/cm3, respectively. The 50, 80 and 140 MeV proton beam energies were selected. The dose distributions of protons in the polymer and radiochromic gel were simulated using Monte Carlo radiation transport code (MCNPX 2.7.0, Los Alamos Laboratory). The water equivalent depth profiles and the dose distributions of two gel dosimeters were compared for the water. Results: In case of irradiating 50, 80 and 140 MeV proton beam to water phantom, the reference Bragg-peak depths are represented at 2.22, 5.18 and 13.98 cm, respectively. The difference in the water equivalent depth is represented to about 0.17 and 0.37 cm in the radiochromic gel and polymer gel dosimeter, respectively. The proton absorbed doses in the radiochromic gel dosimeter are calculated to 2.41, 3.92 and 6.90 Gy with increment of incident proton energies. In the polymer gel dosimeter, the absorbed doses are calculated to 2.37, 3.85 and 6.78 Gy with increment of incident proton energies. The relative absorbed dose in radiochromic gel (about 0.47 %) is similar to that of water than the relative absorbed dose of polymer gel (about 2.26 %). In evaluating the proton dose distribution, we found that the dose distribution of both gel dosimeters matched that of water in most cases. Conclusion: As the dosimetry device, the radiochromic gel dosimeter has the potential particle detectability and is feasible to use for quality assurance of proton beam therapy beam.

07 ISOTOPE AND RADIATION SOURCES↗

Estimating crop proportions from remotely sensed data

The classification/pixel-count method for estimating the proportion of wheat in each segment is theoretically biased even if all distributional assumptions are met. Alternative ways to estimate crop proportions are examined and their performance testing is considered. Topics covered include general linear functional estimates, the method of moments, and maximum likelihood estimators.

Feiveson, A. H.↗

Static shape control for flexible structures

An integrated methodology is described for defining static shape control laws for large flexible structures. The techniques include modeling, identifying and estimating the control laws of distributed systems characterized in terms of infinite dimensional state and parameter spaces. The models are expressed as interconnected elliptic partial differential equations governing a range of static loads, with the capability of analyzing electromagnetic fields around antenna systems. A second-order analysis is carried out for statistical errors, and model parameters are determined by maximizing an appropriate defined likelihood functional which adjusts the model to observational data. The parameter estimates are derived from the conditional mean of the observational data, resulting in a least squares superposition of shape functions obtained from the structural model.

Rodriguez, G.↗

Seasonal Changes in Titan's Surface Temperatures

Cassini's extended mission has provided the opportunity to search for seasonal variations on Titan. In particular, surface temperatures are expected to have shifted significantly in latitude during the completed portion of the mission. Spectra recorded by the Composite Infrared Spectrometer (CIRS) during the nominal mission (2004-08) and the Equinox mission. (2008-10) have already shown changes in temperature. CIRS has detected a seasonal shift in the latitudinal distribution of surface brightness temperatures by comparing zonal averages from two time segments, one period in late northern winter centered on L(sub s) approximately 335 deg and a second period centered on the equinox (L(sub s) approximately 0 deg.). The earlier period had a meridional distribution similar to that previously reported: 93.5 K at the equator, 91.7 K at 85 S and 899 K at 85 N. The newly measured distribution near equinox shows a cooling in the south and a warming in the north, both by about 0.5 K. We estimate that. the centroid of the distribution moved from approximately 16 S to 7 S between the two periods. This gives a seasonal lag behind insolation of delta L(sub s) approximately 13 deg. The CIRS equinox results are consistent with those of Voyager IRIS, which encountered Titan in November 1980, just following the previous northern equinox (L(sub s) = 10 deg.). When compared with predictions from general circulation models, seasonal variations of surface temperature can help constrain the identification of surface materials. Our measurements most closely match the case of a porous ice regolith treated by Tokano, but with some apparent differences between the northern and southern hemispheres. CIRS will extend its study of seasonal variations in surface temperature on Titan as Cassini continues through northern spring.

Jennings, Donald E.↗

The Two-Particle Correlation Function for Systems with Long-Range Interactions

In this paper, we study the truncated two-particle correlation function in particle systems with long range interactions. For Coulombian and soft potentials, we define and give well-posedness results for the equilibrium correlations. In the Coulombian case, we prove the onset of the Debye screening length in the equilibrium correlations, for suitable velocity distributions. Additionally, we give precise estimates on the effective range of interaction between particles. In the case of soft potential interaction the equilibrium correlations and their fluxes in the space of velocities are shown to be linearly stable.

97 MATHEMATICS AND COMPUTING↗

Coseismic fault slip associated with the 1992 M(sub w) 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence

Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Bennett, Richard A.↗

Constraining Black Carbon Aerosol over Asia using OMI Aerosol Absorption Optical Depth and the Adjoint of GEOS-Chem

Accurate estimates of the emissions and distribution of black carbon (BC) in the region referred to here as Southeastern Asia (70degE-l50degE, 11degS-55degN) are critical to studies of the atmospheric environment and climate change. Analysis of modeled BC concentrations compared to in situ observations indicates levels are underestimated over most of Southeast Asia when using any of four different emission inventories. We thus attempt to reduce uncertainties in BC emissions and improve BC model simulations by developing top-down, spatially resolved, estimates of BC emissions through assimilation of OMI observations of aerosol absorption optical depth (AAOD) with the GEOS-Chem model and its adjoint for April and October of 2006. Overwhelming enhancements, up to 500%, in anthropogenic BC emissions are shown after optimization over broad areas of Southeast Asia in April. In October, the optimization of anthropogenic emissions yields a slight reduction (1-5%) over India and parts of southern China, while emissions increase by 10-50% over eastern China. Observational data from in situ measurements and AERONET observations are used to evaluate the BC inversions and assess the bias between OMI and AERONET AAOD. Low biases in BC concentrations are improved or corrected in most eastern and central sites over China after optimization, while the constrained model still underestimates concentrations in Indian sites in both April and October, possibly as a. consequence of low prior emissions. Model resolution errors may contribute up to a factor of 2.5 to the underestimate of surface BC concentrations over northern India. We also compare the optimized results using different anthropogenic emission inventories and discuss the sensitivity of top-down constraints on anthropogenic emissions with respect to biomass burning emissions. In addition, the impacts of brown carbon, the formulation of the observation operator, and different a priori constraints on the optimization are investigated. Overall, despite these limitations and uncertainties, using OMI AAOD to constrain BC sources improves model representation of BC distributions, particularly over China.

Aerosols↗

Evolution of deep-bed filtration of engine exhaust particulates with trapped mass

Micro-scale filtration experiments were performed on cordierite filter samples using particulate matter (PM) generated by a spark-ignition direct-injection (SIDI) engine fueled with tier II EEE certification gasoline. Size-resolved mass and number concentrations were obtained from several engine operating conditions. The resultant mass-mobility relationships showed weak dependence on the operating condition. An integrated particle size distribution (IPSD) method was used estimate the PM mass concentration in the exhaust stream from the SIDI engine and a heavy duty diesel (HDD) engine. The average estimated mass concentration between all conditions was ~77****** % of the gravimetric measurements performed on Teflon filters. Despite the relatively low elemental carbon fraction (~0.4 to 0.7), the IPSD mass for stoichiometric SIDI exhaust was ~83±38 % of the gravimetric measurement. Identical cordierite filter samples with properties representative of diesel particulate filters were sequentially loaded with PM from the different SIDI engine operating conditions, in order of increasing PM mass concentration. Simultaneous particle size distribution measurements upstream and downstream of the filter sample were used to evaluate filter performance evolution and the instantaneous trapped mass within the filter for two different filter face velocities. The evolution of filtration performance for the different samples was sensitive only to trapped mass, despite using PM from a wide range of operating conditions. Higher filtration velocity resulted in a more rapid shift in the most penetrating particle size towards smaller mobility diameters.

33 ADVANCED PROPULSION SYSTEMS↗