Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed state 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Modeling and Evaluation of the Global Sea-Salt Aerosol Distribution: Sensitivity to Emission Schemes and Resolution Effects at Coastal/Orographic Sites

One of the major sources of uncertainty in model estimates of the global sea-salt aerosol distribution is the emission parameterization. We evaluate a new sea-salt aerosol life cycle module coupled to the online multi-scale chemical transport model NMMB/BSC-CTM. We compare 5 year global simulations using five state-of-the-art sea-salt open-ocean emission schemes with monthly averaged coarse aerosol optical depth (AOD) from selected AERONET sun photometers, surface concentration measurements from the University of Miami's Ocean Aerosol Network, and measurements from two NOAA/PMEL cruises (AEROINDOEX and ACE1). Model results are highly sensitive to the introduction of sea-surface-temperature (SST)-dependent emissions and to the accounting of spume particles production. Emission ranges from 3888 teragrams per year to 8114 teragrams per year, lifetime varies between 7.3 hours and 11.3 hours, and the average column mass load is between 5.0 teragrams and 7.2 teragrams. Coarse AOD is reproduced with an overall correlation of around 0.5 and with normalized biases ranging from +8.8 percent to +38.8 percent. Surface concentration is simulated with normalized biases ranging from minus 9.5 percent to plus 28 percent and the overall correlation is around 0.5. Our results indicate that SST-dependent emission schemes improve the overall model performance in reproducing surface concentrations. On the other hand, they lead to an overestimation of the coarse AOD at tropical latitudes, although it may be affected by uncertainties in the comparison due to the use of all-sky model AOD, the treatment of water uptake, deposition and optical properties in the model and/or an inaccurate size distribution at emission.

Bias↗
NASA NTRS2016

The Pilatus Unmanned Aircraft System for Lower Atmospheric Research

This paper presents details of the University of Colorado (CU) "Pilatus" unmanned research aircraft, assembled to provide measurements of aerosols, radiation and thermodynamics in the lower troposphere. This aircraft has a wingspan of 3.2 m and a maximum take-off weight of 25 kg, and it is powered by an electric motor to reduce engine exhaust and concerns about carburetor icing. It carries instrumentation to make measurements of broadband up- and downwelling shortwave and longwave radiation, aerosol particle size distribution, atmospheric temperature, relative humidity and pressure and to collect video of flights for subsequent analysis of atmospheric conditions during flight. In order to make the shortwave radiation measurements, care was taken to carefully position a high-quality compact inertial measurement unit (IMU) and characterize the attitude of the aircraft and its orientation to the upward-looking radiation sensor. Using measurements from both of these sensors, a correction is applied to the raw radiometer measurements to correct for aircraft attitude and sensor tilt relative to the sun. The data acquisition system was designed from scratch based on a set of key driving requirements to accommodate the variety of sensors deployed. Initial test flights completed in Colorado provide promising results with measurements from the radiation sensors agreeing with those from a nearby surface site. Additionally, estimates of surface albedo from onboard sensors were consistent with local surface conditions, including melting snow and bright runway surface. Aerosol size distributions collected are internally consistent and have previously been shown to agree well with larger, surface-based instrumentation. Finally the atmospheric state measurements evolve as expected, with the near-surface atmosphere warming over time as the day goes on, and the atmospheric relative humidity decreasing with increased temperature. No directional bias on measured temperature, as might be expected due to uneven heating of the sensor housing over the course of a racetrack pattern, was detected. The results from these flights indicate that the CU Pilatus platform is capable of performing research-grade lower tropospheric measurement missions.

de Boer, Gijs↗
DOE OSTI2024

Forest Carbon Storage in the Western United States: Distribution, Drivers, and Trends

Abstract Forests are a large carbon sink and could serve as natural climate solutions that help moderate future warming. Thus, establishing forest carbon baselines is essential for tracking climate‐mitigation targets. Western US forests are natural climate solution hotspots but are profoundly threatened by drought and altered disturbance regimes. How these factors shape spatial patterns of carbon storage and carbon change over time is poorly resolved. Here, we estimate live and dead forest carbon density in 19 forested western US ecoregions with national inventory data (2005–2019) to determine: (a) current carbon distributions, (b) underpinning drivers, and (c) recent trends. Potential drivers of current carbon included harvest, wildfire, insect and disease, topography, and climate. Using random forests, we evaluated driver importance and relationships with current live and dead carbon within ecoregions. We assessed trends using linear models. Pacific Northwest (PNW) and Southwest (SW) ecoregions were most and least carbon dense, respectively. Climate was an important carbon driver in the SW and Lower Rockies. Fire reduced live and increased dead carbon, and was most important in the Upper Rockies and California. No ecoregion was unaffected by fire. Harvest and private ownership reduced carbon, particularly in the PNW. Since 2005, live carbon declined across much of the western US, likely from drought and fire. Carbon has increased in PNW ecoregions, likely recovering from past harvest, but recent record fire years may alter trajectories. Our results provide insight into western US forest carbon function and future vulnerabilities, which is vital for effective climate change mitigation strategies.

Environmental Sciences & Ecology↗
NASA NTRS2015

Modelling Groundwater Depletion at Regional and Global Scales: Present State and Future Prospects.

Except for frozen water in ice and glaciers, groundwater is the world's largest distributed store of freshwater and has strategic importance to global food and water security. In this paper, the most recent advances quantifying groundwater depletion (GWD) are comprehensively reviewed. This paper critically evaluates the recently advanced modeling approaches estimating GWD at regional and global scales, and the evidence of feedbacks to the Earth system including sea-level rise associated with GWD. Finally, critical challenges and opportunities in the use of groundwater are identified for the adaption to growing food demand and uncertain climate.

Groundwater Depletion (GWD)↗
NASA NTRS1975

The state of meteoritic material on the moon

Velocity distributions obtained from earth-based meteor observations and corrected to constant meteoroid mass are transformed to meteoroid impact velocity distributions at the moon. Average velocities at the moon thus obtained range from 13 to 18 km/sec. This information is then combined with impact heating estimates to determine the amounts of resident lunar material and incoming meteoritic material that is melted or vaporized. It is found that about 5% of the resulting melt and about 25% of the impact generated vapor should be of meteoritic content. The uncertainties in these estimates are rather large, however, because of uncertainties in the various experimental results required to establish the parameters used in this analysis. The impact generated melts and impact generated vapor deposits, if the latter can be found, should provide additional means for chemically separating out the meteoritic component in the lunar soil.

Zook, H. A.↗
DOE OSTI2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗
DOE OSTI2021

6D Simulations of PIP-II Booster Injection

The PIP-II superconducting linac will deliver 2 mA average H⁻ beam current at 800 MeV to the existing Booster synchrotron over a period of 0.55 ms (285 turns). As a result, the injected beam power will quadruple to 17 kW. Safe operation at the increased beam power implies careful attention to the origin, magnitude, and distribution of both controlled and uncontrolled losses. Uncontrolled losses are due to neutral ions in excited states stripped in downstream magnets and large angle scattered protons from parasitic foil hits. The relative magnitudes of these loss mechanisms is used to determine the optimal foil thickness. A transverse painting scheme involving closed orbit motion will be used to mitigate space charge effects and minimize parasitic foil hits. Using a detailed full 6D simulation of the injection process, we compute large angle scattering losses and compare results to back of the envelope estimates. We investigate possible impact of space charge on the emittance and beam distribution both during and at the conclusion of the injection period.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Spatial Heterogeneity in Cloud Condensation Nuclei (SHINE) Field Campaign Report

Multiple land-atmosphere coupling processes influence the spatial distribution of boundary-layer energy, moisture, aerosol, or cloud condensation nuclei (CCN) properties. Accurately quantifying the spatial and temporal distributions of CCN concentrations is essential for investigating aerosol-cloud interactions and reducing uncertainties in aerosol radiative forcing estimates of Earth's system.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Using probabilistic analysis to assess the reliability of predicted SRB aft-skirt stresses

Probabilistic failure analysis is a tool to predict the reliability of a part or system. Probabalistic techniques were used to predict critical stresses which occur in the solid rocket booster aft-skirt during main engine buildup, immediately prior to lift-off. More than any other hold down post (HDP) load component, the Z loads are sensitive to variations in strains and calibration constants. Also, predicted aft-skirt stresses are strongly affected by HDP load variations. Therefore, the instrumented HDP are not effective load transducers for Z loads, and, when used with aft skirt stress indicator equations, yield estimates with large uncertainty. Monte Carlo simulation proved to be a straight forward way of studying the overlapping effects of multiple parameters on predicted equipment performance. An advantage of probabilistic analysis is the degree of uncertainty of each parameter as stated explicitly by its probability distribution. It was noted, however, that the choice of parameter distribution had a large effect on the simulation results. Many times these distributions must be assumed. The engineer who is designing the part should be responsible for the choice of parameter distribution.

Richardson, James A.↗
NASA NTRS2018

Improving the SMAP Level-4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level 4 Soil Moisture (L4 SM) product. The L4 SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. The L4 SM data are published with a mean latency of approx. 2.5 days from the time of observation and provide global, three-hourly, 9 km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. In 2018, the product was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high quality space-borne remote sensing observations. Land cover inputs were updated to the GlobCover2009 product, which is based on satellite observations from the Medium Resolution Imaging Spectrometer. Topographic statistics now rely on observations from the Shuttle Radar Topography Mission. Finally, vegetation height inputs are derived from space-borne Lidar measurements. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of soil moisture from the models root-zone excess reservoir into the surface excess reservoir. Specifically, the replenishment of soil moisture near the surface from below under non-equilibrium conditions was substantially reduced, which brings the models surface soil moisture more in line with the SMAP Level 2 and in situ soil moisture. Finally, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4. Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not, and on balance, the overall improvement is modest at best. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation minus forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4 SM product and the lessons learned from the effort.

Reichle, Rolf↗
NASA NTRS2007

Quaternion Averaging

Many applications require an algorithm that averages quaternions in an optimal manner. For example, when combining the quaternion outputs of multiple star trackers having this output capability, it is desirable to properly average the quaternions without recomputing the attitude from the the raw star tracker data. Other applications requiring some sort of optimal quaternion averaging include particle filtering and multiple-model adaptive estimation, where weighted quaternions are used to determine the quaternion estimate. For spacecraft attitude estimation applications, derives an optimal averaging scheme to compute the average of a set of weighted attitude matrices using the singular value decomposition method. Focusing on a 4-dimensional quaternion Gaussian distribution on the unit hypersphere, provides an approach to computing the average quaternion by minimizing a quaternion cost function that is equivalent to the attitude matrix cost function Motivated by and extending its results, this Note derives an algorithm that deterniines an optimal average quaternion from a set of scalar- or matrix-weighted quaternions. Rirthermore, a sufficient condition for the uniqueness of the average quaternion, and the equivalence of the mininiization problem, stated herein, to maximum likelihood estimation, are shown.

Markley, F. Landis↗
DOE OSTI2022

Adversarial autoencoder ensemble for fast and probabilistic reconstructions of few-shot photon correlation functions for solid-state quantum emitters

Second-order photon correlation measurements [g (2) (τ) functions] are widely used to classify single-photon emission purity in quantum emitters or to measure the multiexciton quantum yield of emitters that can simultaneously host multiple excitations – such as quantum dots – by evaluating the value of g (2) (τ = 0). Accumulating enough photons to accurately calculate this value is time consuming and could be accelerated by fitting of few-shot photon correlations. Here, we develop an uncertainty-aware, deep adversarial autoencoder ensemble (AAE) that reconstructs noise-free g (2) (τ) functions from noise-dominated, few-shot inputs. The model is trained with simulated g (2) (τ) functions that are facilely generated by Poisson sampling time bins. The AAE reconstructions are performed orders-of-magnitude faster, with reconstruction errors and estimates of g (2) (τ = 0) that are lower in variance and similar in accuracy compared to Maximum likelihood estimation and Levenberg-Marquardt least-squares fitting approaches, for simulated and experimentally measured few-shot g (2) (τ) functions (~100 two-photon events) of InP/ZnS/ZnSe and CdS/CdSe/CdS quantum dots. The deep-ensemble model comprises eight individual autoencoders, allowing for probabilistic reconstructions of noise-free g (2) (τ) functions, and we show that the predicted variance scales inversely with number of shots, with comparable uncertainties to computationally intensive Markov chain Monte Carlo sampling. Furthermore, this work demonstrates the advantage of machine learning models to perform uncertainty-aware, fast, and accurate reconstructions of simple Poisson-distributed photon correlation functions, allowing for on-the-fly reconstructions and accelerated materials characterization of solid-state quantum emitters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Causal CCF Parameter Estimations 2020

This report documents the quantitative results of the causal common-cause failure (CCF) parameter estimations for the failure cause groups “component,” “design,” “environment,” “human,” and “other,” based on CCF data through 2020 in the U.S. Nuclear Regulatory Commission (NRC) CCF database: https://rads.inl.gov/Pages/CCF.aspx. This report utilizes the same data period (2006–2020) and CCF templates as INL/EXT-21-62940, Revision 1, CCF Parameter Estimations, 2020 Update. The 2015 causal CCF prior distributions for the specific failure cause groups (instead of the 2015 generic CCF prior distributions) were used in this report to estimate the associated causal CCF parameters. All the 2015 causal CCF prior distributions and generic CCF prior distributions were developed in INL/EXT-21-43723, Developing Generic Prior Distributions for Common Cause Failure Alpha Factors and Causal Alpha Factors, using CCF data from 1997 to 2015. These quantitative results were developed to support the causal alpha factor model and should be used as appropriate in probabilistic risk assessment (PRA) studies such as the NRC Significance Determination Process for commercial nuclear power plants in the United States.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗
DOE OSTI2023

Constraining the $\overline{K}$N coupled channel dynamics using femtoscopic correlations at the LHC

The interaction of K – with protons is characterised by the presence of several coupled channels, systems like $\overline{K}$ 0 n and πΣ with a similar mass and the same quantum numbers as the K – p state. The strengths of these couplings to the K – p system are of crucial importance for the understanding of the nature of the Λ(1405) resonance and of the attractive K – p strong interaction. In this article, we present measurements of the K – p correlation functions in relative momentum space obtained in pp collisions at $\sqrt{s}$=13 Te, in p–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te, and (semi)peripheral Pb–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te. The emitting source size, composed of a core radius anchored to the K + p correlation and of a resonance halo specific to each particle pair, varies between 1 and 2 fm in these collision systems. The strength and the effects of the $\overline{K}$ 0 n and πΣ inelastic channels on the measured K – p correlation function are investigated in the different colliding systems by comparing the data with state-of-the-art models of chiral potentials. A novel approach to determine the conversion weights ω, necessary to quantify the amount of produced inelastic channels in the correlation function, is presented. In this method, particle yields are estimated from thermal model predictions, and their kinematic distribution from blast-wave fits to measured data. The comparison of chiral potentials to the measured K – p interaction indicates that, while the πΣ –K – p dynamics is well reproduced by the model, the coupling to the $\overline{K}$ 0 n channel in the model is currently underestimated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Development of a Comparison Framework for Evaluating Environmental Contours of Extreme Sea States

Environmental contours of extreme sea states are often utilized for the purposes of reliability-based offshore design. Many methods have been proposed to estimate environmental contours of extreme sea states, including, but not limited to, the traditional inverse first-order reliability method (I-FORM) and subsequent modifications, copula methods, and Monte Carlo methods. These methods differ in terms of both the methodology selected for defining the joint distribution of sea state parameters and in the method used to construct the environmental contour from the joint distribution. It is often difficult to compare the results of proposed methods to determine which method should be used for a particular application or geographical region. The comparison of the predictions from various contour methods at a single site and across many sites is important to making environmental contours of extreme sea states useful in practice. The goal of this paper is to develop a comparison framework for evaluating methods for developing environmental contours of extreme sea states. This paper develops generalized metrics for comparing the performance of contour methods to one another across a collection of study sites, and applies these metrics and methods to develop conclusions about trends in the wave resource across geographic locations, as demonstrated for a pilot dataset. These proposed metrics and methods are intended to judge the environmental contours themselves relative to other contour methods, and are thus agnostic to a specific device, structure, or field of application. The metrics developed and applied in this paper include measures of predictive accuracy, physical validity, and aggregated temporal performance that can be used to both assess contour methods and provide recommendations for the use of certain methods in various geographical regions. The application and aggregation of the metrics proposed in this paper outline a comparison framework for environmental contour methods that can be applied to support design analysis workflows for offshore structures. This comparison framework could be extended in future work to include additional metrics of interest, potentially including those to address issues pertinent to a specific application area or analysis discipline, such as metrics related to structural response across contour methods or additional physics-based metrics based on wave dynamics.

54 ENVIRONMENTAL SCIENCES↗