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 505 records · Page 28

A study of explosive-induced fracture in polymethyl methacrylate (PMMA)

The fracture response of geologic materials is of interest for applications, including geothermal energy harnessing and containment of underground explosions. To better understand the explosively induced fracture response of geomaterials, polymethyl methacrylate (PMMA) was used as a transparent rock surrogate to allow imaging of internal shock propagation and fracture growth processes. Experiments were conducted using high-speed shadowgraphy and photon Doppler velocimetry (PDV), which were compared to numerical simulations. Experiments measured fractures produced in 304.8mm × 304.8mm × 304.8mm PMMA cubes with two simultaneously initiated detonators. The cubes were subjected to varying amounts and directions of externally applied uniaxial stresses, including no stress, 2 MPa stress, and 20 MPa stress. The fracture radius as a function of time was extracted from the high-speed videos. Post-test images of the PMMA cubes aided in the determination of three-dimensional effects not directly imaged by the cameras. The surface velocity history and the shock response captured in PDV and the high-speed videos were compared to the simulated explosive-induced shock response. The simulation results indicate that the shock drives the fracture for the first 20 μs corresponding to a fracture radius of approximately 15 mm in the experiments. The gas-driven fracture extent was estimated analytically using an equilibrium stress distribution calculated after the shock wave propagation through the sample. Reduction in the gas pressure due to the leakage of the explosive products through the crack as a function of time was accounted for. In conclusion, the estimated fracture lengths were in agreement with the experimentally observed fracture lengths.

15 GEOTHERMAL ENERGY↗

Understanding Uncertainties in Tropical Cyclone Rainfall Hazard Modeling Using Synthetic Storms

Tropical cyclone (TC) rainfall hazard assessment is subject to the bias in TC climatology estimation from climate simulations or synthetic downscaling. In this study, we investigate the uncertainty in TC rainfall hazard assessment induced by this bias using both rain gauge and radar observations and synthetic-storm-model-coupled TC rainfall simulations. We identify the storm’s maximum intensity, impact duration, and minimal distance to the site to be the three most important storm parameters for TC rainfall hazard, and the relationship between the important storm parameters and TC rainfall can be well captured by a physics-based TC rainfall model. The uncertainty in the synthetic rainfall hazard induced by the bias in TC climatology can be largely explained by the bias in the important storm parameters simulated by the synthetic storm model. Correcting the distribution of the most biased parameter may significantly improve rainfall hazard estimation. Bias correction based on the joint distribution of the important parameters may render more accurate rainfall hazard estimations; however, the general technical difficulties in resampling from high-dimensional joint probability distributions prevent more accurate estimations in some cases. The results of the study also support future investigation of the impact of climate change on TC rainfall hazards through the lens of future changes in the identified important storm parameters.

54 ENVIRONMENTAL SCIENCES↗

Distributed Energy Resource (DER) Reliability for Backup Electric Power Systems

Hospitals, emergency services, military bases, ports, airports, industries, commercial facilities, and others rely on backup power systems to provide electricity for their critical loads during grid outages. The purpose of this report it to provide accurate reliability information on commonly deployed distributed energy resources (DERs) to improve quantitative estimates for the reliability of these backup power systems during a grid outage. A backup power system consists of DERs, an electric distribution system with its associated switches and other devices, and mechanisms to control and manage the flow of electricity. Too often, facilities and campuses fail to properly quantify the reliability of their backup power systems. DERs are assumed to be 100% reliable, with the only concern being the availability of fuel. Such assumptions can lead to gross errors in the backup system's reliability estimates, particularly for long-duration outages. This report provides a set of estimates for reliability of emergency diesel generators (EDGs), natural gas prime generators and combined heat and power (CHP) prime movers, solar photovoltaics (PV), wind turbines, and Li-ion battery energy storage systems (BESS). The estimates are derived from empirical data when available and supplemented by modeling results when needed. These reliability estimates are for the DERs ability to provide power during a grid outage, ranging from an hour to 2 weeks.

14 SOLAR ENERGY↗

Exclusive J/$\psi$ detection and physics with ECCE

The EIC Comprehensive Chromodynamics Experiment (ECCE) detector has been recommended as a reference design for the proposed Electron-Ion Collider (EIC) program. Here, this paper presents simulation studies of exclusive J/$\psi$ detection and selected physics impact results in EIC using the projected ECCE detector concept. Exclusive quarkonium photoproduction is one of the most popular processes in EIC, which has a large cross section and a simple final state. Due to the gluonic nature of the exchange Pomeron, this process can be related to the gluon distributions in the nucleus. Preliminary results estimate the excellent statistics benefited from the large cross section of J/$\psi$ photoproduction and superior performance of ECCE detector concept. The precise measurement of exclusive J/$\psi$ photoproduction at EIC will help us to more deeply understand nuclear gluon distributions, near threshold production mechanism and nucleon mass structure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Velocity-dependent J-factors for Milky Way dwarf spheroidal analogues in cosmological simulations

We study the impact of the dark matter velocity distribution modelling on signals from velocity-dependent dark matter annihilation in Milky Way dwarf spheroidal galaxies. Using the high resolution APOSTLE simulations, we identify analogues corresponding to Milky Way dwarf spheroidal galaxies, and from these directly determine the dark matter pair-wise relative velocity distribution, and compare to best-fitting Maxwell-Boltzmann distribution models. For three velocity-dependent annihilation models, p-wave, d-wave, and the Sommerfeld model, we quantify the errors introduced when using the Maxwell-Boltzmann parameterization. We extract a simple power-law relation between the maximum circular velocity of the dwarf spheroidal analogue and the peak speed of the Maxwell-Boltzmann distribution. We show that this relation can be used to accurately calculate the dark matter relative velocity distribution, and find that it allows us to estimate the dark matter annihilation signal without the need to directly calculate the relative velocity distribution for each galaxy. The scatter in the J-factors calculated from the analogues dominates the uncertainty obtained when compared to the J-factor as determined from the observational data for each dwarf spheroidal, with the largest scatter from d-wave models and the smallest from Sommerfeld models.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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↗

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↗

Integrating Overhead Camera Imagery with Reinforcement Learning to Improve Fuel Economy Through Adaptive Traffic Control

Overhead cameras have been deployed for some time as a sensor for traffic control, with many advantages over conventional inductive loop sensors in terms of robustness, ease of installation, and simplicity with respect to maintenance. We explore the utility of such imagery (in particular, fisheye lens cameras) for intelligent traffic control in real-time using reinforcement learning, with a goal of improving fuel economy for the overall intersection and grid of intersections. We summarize our initial work in data collection and approach to ascertain the ability to visually sense vehicular fuel consumption using real-world data. We then describe our deep reinforcement learning topology to test adaptive control informed by the camera capabilities. Next, we tested the proof-of-concept through a variety of simulation scenarios using the Simulation of Urban Mobility (SUMO) traffic simulator, which included different estimates on vehicle classifications, different traffic distributions, and different traffic densities. Six different control policies were tested ranging from simple timing, to heuristic visual-based methods, to actual policies learned via the reinforcement learning approach. Our results demonstrate the feasibility of the approach for controlling traffic lights for better fuel efficiency based on visual vehicle estimates from commercial, fisheye lens cameras.

Karnowski, Thomas↗

Home range and resource selection of Virginia opossums in the rural southeastern United States

The Virginia opossum (Didelphis virginiana) has a rapidly expanding distribution in North America, but many aspects of its ecology remain relatively understudied, particularly in rural areas of its core range. We collected GPS telemetry data from 93 opossums in a rural, non-agricultural landscape in South Carolina, USA (2018–2019) to examine factors influencing space use and resource selection. Estimated male home ranges (99% utilization distributions) were on average 50% larger than those of females (mean home range 115.9 ± 103.7 ha vs 76.7 ± 75.0 ha). The home range size decreased on average by 20% with each 20% increase in deciduous land cover but was not affected by season or other landscape factors. Core area sizes (65% utilization distributions) were not influenced by sex (mean core area size 29.1 ± 23.7 ha and 22.4 ha ± 13.8 for males and females, respectively) or season, but the core area size decreased by 14% with each 400 m increase in distance from a permanent water source. Resource selection by opossums primarily occurred at the landscape level. Both males and females generally selected for wetlands while avoiding pine forests and developed/open areas, likely the result of differences in resource availability and predation risk between habitats. Opossums also tended to select for linear features such as unpaved roads and edge habitat, which may facilitate movement across the landscape. Finally, the home ranges we documented are among the largest recorded for opossums in the USA, likely the result of the relatively low resource abundance throughout our study area due to comparatively minimal anthropogenic influence.

59 BASIC BIOLOGICAL SCIENCES↗

Uranium particle age dating, aggregation, and model age best estimators

We present important aspects of uranium particle age dating by Large-Geometry Secondary Ion Mass Spectrometry (LG-SIMS) that can introduce bias and increase model age uncertainties, especially for small, young, and/or low-enriched particles. This metrology is important for applications related to International Nuclear Safeguards. We explore influential factors related to model age estimation, including the effects of evolving surface chemistry on inter-element measurements of particles (e.g., Th and U), detector background, and aggregation methods using simulated and actual particle samples. We introduce a new model age estimator, called “mid68”, that supplements 95% confidence intervals, providing a “best estimate” and uncertainty about the most likely age. The mid68 estimator can be calculated using the Feldman and Cousins method or Bayesian methods and provides a value with a symmetric uncertainty that can be used for calculations and approximate aggregation of processed model age values when the raw data and correction factors are not available. For particles yielding low 230 Th counts amidst nonzero detector background, their underlying model age probability distributions are asymmetric, so the mid68 estimator provides additional robust information regarding the underlying model age likelihood. This study provides a comprehensive and timely examination of critical aspects of uranium particle age dating as more laboratories establish particle chronometry capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GPS Spoofing Mitigation and Timing Risk Analysis in Networked Phasor Measurement Units via Stochastic Reachability

To address phasor measurement unit (PMU) vulnerability to spoofing, we propose the use of a set-valued state estimation technique known as stochastic reachability (SR)-based distributed Kalman filter (DKF) that computes secure global positioning system (GPS) timing across a network of receivers. Utilizing SR, we estimate not only GPS time but also its stochastic reachable set, which is parameterized by probabilistic zonotope (p-Zonotope). While requiring known measurement error bounds in only non-spoofed conditions, we designed a two-tiered approach. We first performed measurement-level spoofing mitigation via deviation of a measurement innovation from its expected p-Zonotope. We then performed state-level timing risk analysis via a determination of the intersection probability of the estimated p-Zonotope with an unsafe set that violates IEEE C37.118.1a-2014 standards. Finally, we validated our SR-DKF algorithm by subjecting it to a simulated receiver network to coordinate signal-level spoofing. We demonstrate improved timing accuracy and successful spoofing mitigation via the use of our SR-DKF algorithm. We also validated the robustness of the estimated timing risk as the number of receivers were varied.

47 OTHER INSTRUMENTATION↗

Fast Automatic Knot Placement Method for Accurate B-spline Curve Fitting

The choice of knot vector has immense influence on the resulting accuracy of a B-spline approximation of a curve. However, despite the significance of this problem and the various solutions that were proposed in the literature, optimizing the number and placement of knots remains a difficult task. This paper presents a novel method for the approximation of a curve by a B-spline of arbitrary order, which automatically determines a knot vector that achieves high approximation quality. Additionally, at the core of our approach is a feature function that characterizes the amount and spatial distribution of geometric details in the input curve by estimating its derivatives. Knots are then selected in such a way as to evenly distribute the feature contents across their intervals. A comparison to the state of the art for a wide variety of curves shows that our method is faster and achieves more accurate reconstruction results, while typically reducing the number of necessary knots.

97 MATHEMATICS AND COMPUTING↗