Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty propagation”

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 847 records · Page 47

Integrating deep neural networks with full-waveform inversion: Reparameterization, regularization, and uncertainty quantification

Full-waveform inversion (FWI) is an accurate imaging approach for modeling the velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong nonlinearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We have adopted a neural-network-based full-waveform inversion (NNFWI) method that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. Furthermore, the velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model determine that NNFWI can mitigate local minima, especially for imaging high-contrast features such as salt bodies, and it significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.

58 GEOSCIENCES↗

Effect of Discontinuities and Uncertainties on the Response and Failure of Composite Structures

The overall goal of this research was to assess the effect of discontinuities and uncertainties on the nonlinear response and failure of composite structures subjected to combined mechanical and thermal loads. The four key elements of the study were: (1) development of simple and efficient procedures for the accurate determination of transverse shear and transverse normal stresses in structural sandwiches as well as in unstiffened and stiffened composite panels and shells; (2) study the effects of transverse stresses on the response, damage initiation and propagation in composite and sandwich structures; (3) use of hierarchical sensitivity coefficients to identify the major parameters that affect the response and damage in each of the different levels in the hierarchy (micro-mechanical, layer, panel, subcomponent and component levels); and (4) application of fuzzy set techniques to identify the range and variation of possible responses. The computational models developed were used in conjunction with experiments, to understand the physical phenomena associated with the nonlinear response and failure of composite and sandwich structures. A toolkit was developed for use in conjunction with deterministic analysis programs to help the designer in assessing the effect of uncertainties in the different computational model parameters on the variability of the response quantities.

Noor, Ahmed K.↗

Challenges in Modeling the Sun-Earth System

The transfer of mass, energy and momentum through the coupled Sun-Earth system spans a wide range of scales in time and space. While profound advances have been made in modeling isolated regions of the Sun-Earth system, minimal progress has been achieved in modeling the end-to-end system. Currently, end-to-end modeling of the Sun-Earth system is a major goal of the National Space Weather and NASA Living With a Star (LWS) programs. The uncertainty in the underlying physics responsible for coupling contiguous regions of the Sun-Earth system is recognized as a significant barrier to progress. Our limited understanding of the underlying coupling physics is illustrated by the following example questions: how does the propagation of a typical CME/solar flare influence the measured properties of the solar wind at 1 AU? How does the solar wind compel the dynamic response of the Earth's magnetosphere? How is variability in the ionosphere-thermosphere system coupled to magnetospheric variations? Why do these and related important questions remain unanswered? What are the primary problems that need to be resolved to enable significant progress in comprehensive modeling of the Sun-Earth system? Which model/technique improvements are required and what new data coverage is required to enable full model advances? This poster opens the discussion for how these and other important questions can be addressed. A workshop scheduled for October 8-22, 2004 in Huntsville, Alabama, will be a forum for identifying ana exploring promising new directions and approaches for characterizing and understanding the system. To focus the discussion, the workshop will emphasize the genesis, evolution, propagation and interaction of high-speed solar wind streamers or CME/flares with geospace and the subsequent response of geospace from its outer reaches in the magnetosphere to the lower edge of the ionosphere-mesosphere-thermosphere. Particular emphasis will be placed on modeling the coupling aspects of these phenomena across boundaries between regions and on data analysis that guides and constrains model results. Specific topics to be addressed are: Corotating interaction regions, Coronal mass ejections, Energetic particles, System preconditioning, Extreme events and super storms, End-to-End modeling efforts.

Spann, James↗

Observations of the Li, Be, and B Isotopes and Constraints on Cosmic-ray Propagation

The abundance of Li, Be, and B isotopes in galactic cosmic rays (GCR) between E=50-200 MeV/nucleon has been observed by the Cosmic Ray Isotope Spectrometer (CRIS) on NASA's ACE mission since 1997 with high statistical accuracy. Precise observations of Li, Be, B can be used to constrain GCR propagation models. We find that a diffusive reacceleration model with parameters that best match CRIS results (e.g. B/C, Li/C, etc) are also consistent with other GCR observations. A approx. 15-20% overproduction of Li and Be in the model predictions is attributed to uncertainties in the production cross-section data. The latter becomes a significant limitation to the study of rare GCR species that are generated predominantly via spallation.

deNolfo, G. A.↗

The Influence of β-decay Rates on r-process Observables

Abstract The rapid neutron capture process ( r -process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via the r -process. Given that the r -process occurs far from nuclear stability, unmeasured beta-decay rates play an essential role in setting the timescale for the r -process. In an effort to better understand the sensitivity of kilonova modeling to different theoretical global beta-decay descriptions, we incorporate these into nucleosynthesis calculations. We compare the results of these calculations and highlight differences in kilonova nuclear energy generation and light-curve predictions, as well as final abundances and their implications for nuclear cosmochronometry. We investigate scenarios where differences in beta-decay rates are responsible for increased nuclear heating on timescales of days that propagates into a significantly increased average bolometric luminosity between 1 and 10 days post-merger. We identify key nuclei, both measured and unmeasured, whose decay rates directly impact nuclear heating generation on timescales responsible for light-curve evolution. We also find that uncertainties in beta-decay rates significantly impact age estimates from cosmochronometry.

79 ASTRONOMY AND ASTROPHYSICS↗

Satellite test of the isotropy of the one-way spe ed of light using ExTRAS

A test of the second postulate of special relativity, the universality of the speed of light, using the ExTRAS (Experiment on Timing Ranging and Atmospheric Sounding) payload to be flown on board a Russian Meteor-3M satellite (launch date January 1997) is proposed. The propagation time of a light signal transmitted from one point to another without reflection would be measured directly by comparing the phases of two hydrogen maser clocks, one on board and one on the ground, using laser or microwave time transfer systems. An estimated uncertainty budget of the proposed measurements is given, resulting in an expected sensitivity of the experiment of delta c/c is less than 8xl0(exp -10) which would be an improvement by a factor of approximately 430 over previous direct measurements and by a factor of approximately 4 over the best indirect measurement. The proposed test would require no equipment additional to what is already planned and so is of inherently low-cost. It could be carried out by anyone having access to a laser or microwave ground station and a hydrogen maser.

Wolf, Peter↗

Longitudinal Variations in Jupiter's Winds

Long-term studies of Jupiter's zonal wind field revealed temporal variations on the order of 20 to 40 m/s at many latitudes, greater than the typical data uncertainties of 1 to 10 m/s. No definitive periodicities were evident, however, though some latitudinally-confined signals did appear at periods relevant to the Quasi- Quadrennial Oscillation (Simon-Miller & Gierasch, Icarus, in press). As the QQO appears, from vertical temperature profiles, to propagate downward, it is unclear why a signal is not more obvious, unless other processes dominate over possibly weaker forcing from the QQO. An additional complication is that zonal wind profiles represent an average over some particular set of longitudes for an image pair and most data sets do not offer global wind coverage. Lien avoiding known features, such as the large anticyclonic vortices especially prevalent in the south, there can be distinct variations in longitude. We present results on the full wind field from Voyager and Cassini data, showing apparent longitudinal variations of up to 60 m/s or more. These are particularly obvious near disruptions such as the South Equatorial Disturbance, even when the feature itself is not clearly visible. These two dates represent very different states of the planet for comparison: Voyagers 1 & 2 flew by Jupiter shortly after a global upheaval, while many regions were in a disturbed state, while the Cassini view is typical of a more quiescent period present during much of the 1990s and early 2000s.

Simon-Miller, Amy A.↗

NASA’s Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA’s Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible-band photometry and near-infrared spectroscopy of exoplanet host stars. Here, we present an end-to-end simulation study quantifying Pandora’s ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of joint visible photometry (0.4–0.7 μm) and near-infrared spectroscopy (0.9–1.6 μm, R ≈ 120) recover photospheric temperatures with typical uncertainties of ≈30 K, with no significant bias. Models with two spectral components (i.e., a quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ≈0.3%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of 10 2 –10 3 ppm are reduced to ≲10 ppm—well below Pandora’s expected transmission spectroscopy precision (30–100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the 10 3 ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar–planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.

Astronomy and AstroPhysics↗

Jovian upper atmospheric temperature measurement by the Voyager 1 UV spectrometer

A measurement of the neutral temperature in the exosphere of Jupiter is reported. The observations were carried out by the Voyager 1 Ultraviolet Spectrometer during the solar occultation phase beginning 4.5 hours after the Jupiter encounter. The range of the exospheric temperature is found to be 1450 (+300, -250) K. The uncertainty in the determination of the temperature appears to be largely a consequence of the lack of knowledge of the brightness distribution on the solar disk at the time of the observations. The high exospheric temperature appears to result from a combination of upward-propagating inertia gravity waves, magnetospheric soft electrons, and auroral electrons, including those from the Io-plasma torus.

Atreya, S. K.↗

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗

Experimental test of general relativity using time-delay data from Mariner 6 and Mariner 7

Range and Doppler data from Mariner 6 and Mariner 7 have been analyzed for purposes of measuring the effect of the sun's gravity field on S-band radio transmissions to the spacecraft. The prediction of general relativity, that the round-trip time delays between the station and the spacecraft will be increased by about 200 microsec near superior conjunction, has been verified with an uncertainty of 3 per cent or less. The dominant error source on the experiment is the stochastic nongravitational forces which act on the spacecraft and which limit the accuracy of the determination of the Mariner orbits. Effects of free coronal electrons on the round-trip propagation also make the measurement difficult, but the coronal contribution to the error is only about 1 per cent. Separate analyses of data from Mariner 6 and Mariner 7 yield two values near unity for the parameter gamma* which differ by about 0.3 per cent from each other. This and other considerations suggest that the error in the experiment is perhaps less than 3 per cent.

Anderson, J. D.↗

Long-distance nuclear matrix elements for neutrinoless double-beta decay from lattice QCD

Neutrinoless double-beta ( 0 ν β β ) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles. Interpretations of the stringent experimental constraints on 0 ν β β -decay half-lives require calculations of nuclear matrix elements. This work presents the first lattice quantum chromodynamics (LQCD) calculation of the matrix element for 0 ν β β decay in a multinucleon system, specifically the n n → p p e e transition, mediated by a light left-handed Majorana neutrino propagating over nuclear-scale distances. This calculation is performed with quark masses corresponding to a pion mass of m π = 806 MeV at a single lattice spacing and volume. The statistically cleaner Σ − → Σ + e e transition is also computed in order to investigate various systematic uncertainties. The prospects for matching the results of LQCD calculations onto a nuclear effective field theory to determine a leading-order low-energy constant relevant for 0 ν β β decay with a light Majorana neutrino are investigated. This work, therefore, sets the stage for future calculations at physical values of the quark masses that, combined with effective field theory and nuclear many-body studies, will provide controlled theoretical inputs to experimental searches of 0 ν β β decay. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fully Bayesian Analysis With Model Inadequacy Correction For Nuclear Graphite Property Models With Hierarchical Variance Structure

Nuclear-grade graphites are extensively utilized in the core designs of various advanced nuclear reactors. Within the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, such as nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity, forecasting component performance over the reactor's lifespan, and developing design standards necessitate robust tools for predicting fracture initiation and propagation in graphite structural components within nuclear reactors. This code enables the Bayesian calibration of properties for nuclear-grade graphites. Using a hierarchical Bayesian approach, multiple experimental data sources are combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due inadequacies in the model and the inherent spread in the experimental data are quantified.

Dhulipala, Som Lakshmi NarasimhaLakshmi Narasimha ↗

Implications of Projected Climate Change for Groundwater Recharge in the Western United States

Existing studies on the impacts of climate change on groundwater recharge are either global or basin/ location-specific. The global studies lack the specificity to inform decision making, while the local studies do little to clarify potential changes over large regions (major river basins, states, or groups of states), a scale often important in the development of water policy. An analysis of the potential impact of climate change on groundwater recharge across the western United States (west of 100 degrees longitude) is presented synthesizing existing studies and applying current knowledge of recharge processes and amounts. Eight representative aquifers located across the region were evaluated. For each aquifer published recharge budget components were converted into four standard recharge mechanisms: diffuse, focused, irrigation, and mountain-systems recharge. Future changes in individual recharge mechanisms and total recharge were then estimated for each aquifer. Model-based studies of projected climate-change effects on recharge were available and utilized for half of the aquifers. For the remainder, forecasted changes in temperature and precipitation were logically propagated through each recharge mechanism producing qualitative estimates of direction of changes in recharge only (not magnitude). Several key patterns emerge from the analysis. First, the available estimates indicate average declines of 10-20% in total recharge across the southern aquifers, but with a wide range of uncertainty that includes no change. Second, the northern set of aquifers will likely incur little change to slight increases in total recharge. Third, mountain system recharge is expected to decline across much of the region due to decreased snowpack, with that impact lessening with higher elevation and latitude. Factors contributing the greatest uncertainty in the estimates include: (1) limited studies quantitatively coupling climate projections to recharge estimation methods using detailed, process-based numerical models; (2) a generally poor understanding of hydrologic flowpaths and processes in mountain systems; (3) difficulty predicting the response of focused recharge to potential changes in the frequency and intensity of extreme precipitation events; and (4) unconstrained feedbacks between climate, irrigation practices, and recharge in highly developed aquifer systems.

longitude↗

Tonal Noise Prediction of a Distributed Propulsion Unmanned Aerial Vehicle

Noise is expected to be a major barrier of unmanned aerial vehicles (UAVs) to public acceptance. A noise prediction scheme is introduced in this paper and applied to a specific vehicle configuration, namely, the Greased Lightning-10. Results herein will be used to demonstrate the feasibility of incorporating a noise constraint within the flight control system of a distributed electric propulsion vehicle by modifying commands for low-noise operation over sensitive areas, e.g., communities, schools, etc. Steady loading and thickness noise signatures of single propellers are computed using the Propeller Analysis System of the NASA Aircraft NOise Prediction Program. The individual signatures are then superposed at far field observers after applying corrections for spherical spreading and phase delays based on straight ray propagation. Two-propeller source fields are verified using analytical directivity patterns of monopoles. Notional effects of rotation rate, rotation direction, and relative propeller phase are given. Under ideal circumstances and equivalent RPM, random phasing, which occurs in most small UAVs, can produce up to 20 decibels uncertainty in the tonal sound pressure level at a given ground observer. Additionally, directivity modification via relative propeller phase control is shown to have great potential as a noise reduction technique. This paper will focus on the forward flight mode but will also briefly discuss the vertical flight mode.

Pascioni, Kyle A.↗

Locating Operational Events of the Cooling Tower of a Nuclear Reactor with a Very Local Seismic Network

Geolocation of emergent seismic signals is challenging at close distances. Here, we used three-component data from a seismic network and a targeted experiment at a research nuclear reactor to locate seismic sources. Utilizing known events collected during the targeted experiment, we were able to infer source locations with seismic amplitudes and polarization characteristics of the data. Although the resolution of the source location is not perfect, the seismic amplitudes and polarization analysis offer useful constraints. For the known events, the source region inferred with our analysis includes the true source locations. Synthetic tests indicate the resolution is largely due to limited data coverage and measurement uncertainties because the synthetic tests show similar results compared with the field data. We identified the source of the unknown event through spectrum cross correlation between the signals from the known events and an unknown event. Our findings were confirmed by operational staff at the facility. When the propagation medium properties (i.e., seismic velocity and quality factor for attenuation) are known, our analysis can be applied to continuous data from a seismic array to infer both source amplitude and location. If the medium properties are not known, a targeted experiment can be conducted to estimate them.

58 GEOSCIENCES↗

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗

An investigation of techniques for the measurement and interpretation of cosmic ray isotopic abundances

An instrument, the Caltech High Energy Isotope Spectrometer Telescope was developed to measure isotopic abundances of cosmic ray nuclei by employing an energy loss - residual energy technique. A detailed analysis was made of the mass resolution capabilities of this instrument. A formalism, based on the leaky box model of cosmic ray propagation, was developed for obtaining isotopic abundance ratios at the cosmic ray sources from abundances measured in local interstellar space for elements having three or more stable isotopes, one of which is believed to be absent at the cosmic ray sources. It was shown that the dominant sources of uncertainty in the derived source ratios are uncorrelated errors in the fragmentation cross sections and statistical uncertainties in measuring local interstellar abundances. These results were applied to estimate the extent to which uncertainties must be reduced in order to distinguish between cosmic ray production in a solar-like environment and in various environments with greater neutron enrichments.

Wiedenbeck, M. E.↗