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 937 records · Page 52

An Example Crossover Experiment for Testing New Vicarious Calibration Techniques for Satellite Ocean Color Radiometry

Vicarious calibration of ocean color satellites involves the use of accurate surface measurements of water-leaving radiance to update and improve the system calibration of ocean color satellite sensors. An experiment was performed to compare a free-fall technique with the established MOBY measurement. We found in the laboratory that the radiance and irradiance instruments compared well within their estimated uncertainties for various spectral sources. The spectrally averaged differences between the NIST values for the sources and the instruments were less than 2.5% for the radiance sensors and less than 1.5% for the irradiance sensors. In the field, the sensors measuring the above-surface downwelling irradiance performed nearly as well as they had in the laboratory, with an average difference of less than 2%. While the water-leaving radiance, L(sub w) calculated from each instrument agreed in almost all cases within the combined instrument uncertainties (approximately 7%), there was a relative bias between the two instrument classes/techniques that varied spectrally. The spectrally averaged (400 nm to 600 nm) difference between the two instrument classes/techniques was 3.1 %. However the spectral variation resulted in the free fall instruments being 0.2% lower at 450 nm and 5.9% higher at 550 nm. Based on the analysis of one matchup, the bias in the L(sub w), was similar to that observed for L(sub u)(1 m) with both systems, indicating the difference did not come from propagating L(sub u)(1 m) to L(sub w).

Voss, Kenneth J.↗

Recent developments of the NESSUS probabilistic structural analysis computer program

The NESSUS probabilistic structural analysis computer program combines state-of-the-art probabilistic algorithms with general purpose structural analysis methods to compute the probabilistic response and the reliability of engineering structures. Uncertainty in loading, material properties, geometry, boundary conditions and initial conditions can be simulated. The structural analysis methods include nonlinear finite element and boundary element methods. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. The scope of the code has recently been expanded to include probabilistic life and fatigue prediction of structures in terms of component and system reliability and risk analysis of structures considering cost of failure. The code is currently being extended to structural reliability considering progressive crack propagation. Several examples are presented to demonstrate the new capabilities.

Millwater, H.↗

Computational Modeling of Meteor-Generated Ground Pressure Signatures

We present a thorough validation of a computational approach to predict infrasonic signatures of centimeter-sized meteoroids. We assume that the energy deposition along the meteor trail is dominated by atmospheric drag and simulate the steady, inviscid flow of air in thermochemical equilibrium to compute the meteoroid's near-body pressure signature. This signature is then propagated through a stratified and windy atmosphere to the ground using a methodology adapted from aircraft sonic-boom analysis. An assessment of the numerical accuracy of the near field and the far field solver is presented. The results show that when the source of the signature is the cylindrical Mach-cone, the simulations closely match the observations. The prediction of the shock rise-time, the zero-peak amplitude of the waveform, and the duration of the positive pressure phase are consistently within 10% of the measurements. Uncertainty in the shape of the meteoroid results in a poorer prediction of the trailing part of the waveform. Overall, our results independently verify energy deposition estimates deduced from optical observations.

Computational↗

Intercomparison of Middle Atmospheric Meteorological Analyses for the Northern Hemisphere Winter 2009-2010

Detailed meteorological analyses based on observations extending through the middle atmosphere (~15 to 100 km altitude) can provide key information to whole atmosphere modelling systems regarding the physical mechanisms linking day-to-day changes in ionospheric electron density to meteorological variability near the Earth’s surface. However, the extent to which independent middle atmosphere analyses differ in their representation of wave-induced coupling to the ionosphere is unclear. To begin to address this issue, we present the first intercomparison among four such analyses, JAGUAR-DAS, MERRA-2, NAVGEM-HA, and WACCMX+DART, focusing on the Northern Hemisphere (NH) 2009-2010 winter, which includes a major sudden stratospheric warming (SSW). This intercomparison examines the altitude, latitude, and time dependences of zonal mean zonal winds and temperatures among these four analyses over the 1 December 2009 – 31 March 2010 period, as well as latitude and altitude dependences of monthly mean amplitudes of the diurnal and semidiurnal migrating solar tides, the eastward propagating diurnal zonal wave number 3 nonmigrating tide, and traveling planetary waves associated with the quasi-5-day and quasi-2-day Rossby modes. Our results show generally good agreement among the four analyses up to the stratopause (~50 km altitude). Large discrepancies begin to emerge in the mesosphere and lower thermosphere owing to (1) differences in the types of satellite data assimilated by each system and (2) differences in the details of the global atmospheric models used by each analysis system. The results of this intercomparison provide initial estimates of uncertainty in analyses commonly used to constrain middle atmospheric meteorological variability in whole atmosphere model simulations.

John P Mccormack↗

Probabilistic Prediction of Geomagnetic Storms and the K p index

Geomagnetic activity is often described using summary indices to summarize the likelihood of space weather impacts, as well as when parameterizing space weather models. The geomagnetic index K p in particular, is widely used for these purposes. Current state-of-the-art forecast models provide deterministic K p predictions using a variety of methods – including empirically-derived functions, physics-based models, and neural networks – but do not provide uncertainty estimates associated with the forecast. This paper provides a sample methodology to generate a 3-hour-ahead K p prediction with uncertainty bounds and from this provide a probabilistic geomagnetic storm forecast. Specifically, we have used a two-layered architecture to separately predict storm (K p ≥ 5 – ) and non-storm cases. As solar wind-driven models are limited in their ability to predict the onset of transient-driven activity we also introduce a model variant using solar X-ray flux to assess whether simple models including proxies for solar activity can improve the predictions of geomagnetic storm activity with lead times longer than the L1-to-Earth propagation time. By comparing the performance of these models we show that including operationally-available information about solar irradiance enhances the ability of predictive models to capture the onset of geomagnetic storms and that this can be achieved while also enabling probabilistic forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗

Preliminary analyses of space radiation protection for lunar base surface systems

Radiation shielding analyses are performed for candidate lunar base habitation modules. The study primarily addresses potential hazards due to contributions from the galactic cosmic rays. The NASA Langley Research Center's high energy nucleon and heavy ion transport codes are used to compute propagation of radiation through conventional and regolith shield materials. Computed values of linear energy transfer are converted to biological dose-equivalent using quality factors established by the International Commision of Radiological Protection. Special fluxes of heavy charged particles and corresponding dosimetric quantities are computed for a series of thicknesses in various shield media and are used as an input data base for algorithms pertaining to specific shielded geometries. Dosimetric results are presented as isodose contour maps of shielded configuration interiors. The dose predictions indicate that shielding requirements are substantial, and an abbreviated uncertainty analysis shows that better definition of the space radiation environment as well as improvement in nuclear interaction cross-section data can greatly increase the accuracy of shield requirement predictions.

Nealy, John E.↗

Optical PAyload for Lasercomm Science (OPALS) link validation

Recently several day and nighttime links under diverse atmospheric conditions were completed using the Optical Payload for Lasercomm Science (OPALS) flight system on-board the International Space Station (ISS). In this paper we compare measured optical power and its variance at either end of the link with predictions that include atmospheric propagation models. For the 976 nm laser beacon mean power transmitted from the ground to the ISS the predicted mean irradiance of 10's of microwatts per square meter close to zenith and its decrease with range and increased air mass shows good agreement with predictions. The irradiance fluctuations sampled at 100 Hz also follow the expected increase in scintillation with air mass representative of atmospheric coherence lengths at zenith at 500 nm in the 3-8 cm range. The downlink predicted power of 100's of nanowatts was also reconciled within the uncertainty of the atmospheric losses. Expected link performance with uncoded bit-error rates less than 1E-4 required for the Reed-Solomon code to correct errors for video, text and file transmission was verified. The results of predicted and measured powers and fluctuations suggest the need for further study and refinement.

ISS↗

Isomer-Dependent Reaction Mechanisms of Cyclic Ether Intermediates: cis -23-Dimethyloxirane and trans -23-Dimethyloxirane.

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-˙QOOH) or from reactions of HÒO with alkenes. The cis- and trans-isomers of 2,3-dimethyloxirane are intermediates of n-butane oxidation, and while rate coefficients for β-˙QOOH → 2,3-dimethyloxirane + OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3-dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomerdependence of 2,3-dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl-initiated oxidation of both cis-2,3-dimethyloxirane and trans-2,3-dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ˙R + O 2 surfaces were computed at the ccCA-PS3 level of theory. In total, 28 barrier heights were computed on the 2,3-dimethyloxiranylperoxy surfaces. Two notable aspects are low-lying pathways that form resonance-stabilized ketohydroperoxide-type radicals caused by ˙QOOH ring-opening when the unpaired electron is localized adjacent to the ether group, and cis-trans isomerization of ˙R and ˙QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3-dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ˙R + O 2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide-type radicals. The present work provides the first analysis of 2,3-dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relativity time-delay experiments utilizing 'Mariner' spacecraft

Relativity predicts that the transit time of a signal propagated from the earth to a spacecraft and retransmitted back to earth ought to exhibit an additional, variable time delay. The present work describes some of the analytical techniques employed in experiments using Mariner spacecraft designed to test the accuracy of this prediction. Two types of data are analyzed in these relativity experiments; these include phase-coherent, two-way Doppler shift and round-trip, transit-time measurements. Results of Mariner 6 and 7 relativistic time-delay experiments are in agreement with Einstein's theory of general relativity with an uncertainty of 3%.

Esposito, P. B.↗

Searches for neutrinos from cosmic-ray interactions in the Sun using seven years of IceCube data

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SAνs) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SAνs also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. In this work, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SAν flux is observed. After inclusion of systematic uncertainties, we set a 90% upper limit of 1.02 +0.20 -0.18 ·10 -13 GeV-1cm -2 s -1 at 1 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaporation and propagation of liquid drop streams at vacuum pressures: Experiments and modeling

Evaporation of streams of liquid droplets in environments at vacuum pressures below the vapor pressure has not been widely studied. Here, experiments and simulations are reported that quantify the change in droplet diameter when a steady stream of ≈100 μm diameter drops is injected into a chamber initially evacuated to <10 –8 bar. In experiments, droplets fall through the center of a 0.8 m long liquid nitrogen cooled shroud, simulating infinity radiation and vapor mass flux boundary conditions. Furthermore, experimentally measured changes in drop diameters vary from ≈0 to 6 μm when the initial vapor pressure is increased from 10 –6 to 10 –3 bar by heating the liquid. Measured diameter changes are predicted by a model based on the Hertz-Knudsen equation. One uncertainty in the calculation is the “sticking coefficient” β. Assuming a constant β for all conditions studied here, predicted diameter changes best match measurements with β≈0.3. This value falls within the range of β reported in the literature for organic liquids. Finally, at the higher vapor pressure conditions considered here, the drop stream disperses transverse to the main flow direction. This spread is attributed to forces imparted by an absolute pressure gradient produced by the evaporating stream.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Are cosmic rays effective for ionization of protoplanetary disks?

The principal uncertainty in studying the magnetic properties of protoplanetary disks concerns the ionization levels of the disk's gas. The low gas temperature precludes thermal ionization, leaving cosmic rays as the dominant source of ionization. It has been shown that the resulting electrical conductivity is just high enough for a MHD dynamo to produce contemporaneously a magnetic field in most of the extended parts of a turbulent protoplanetary disk. Here we argue that the effectiveness of cosmic rays to ionize the bulk of the gas is impaired by the influence of the generated magnetic field on the propagation of cosmic rays within a disk. Cosmic rays scatter on magnetic inhomogeneities, and, as a result, their penetration depth decreases to only a fraction of the disk half-thickness, resulting in a severe depletion of free charge from the midplane regions of a disk. That, in turn, undercuts the dynamo mechanism, so extended parts of the disk are free from a dynamically significant magnetic field. We also point out that any additional, even small, in situ source of ionization, such as radioactive Al-26, may again make a dynamo a viable regeneration process capable of producing a dynamically important magnetic field.

Dolginov, Arkady Z.↗

Toward estimation of H(sub 0) from VLBI observations of the gravitational lens system 0957+561

We have analyzed two epochs of 6 cm very long base interferometry (VLBI) observations of the two brightest images, A and B, of the gravitationally lensed quasar 0957+561 in a search for correlated structural evolution from which to estimate delta tau(sub BA), the difference in propagation time for light traveling between the quasar and Earth along the paths corresponding to these two images. However, we detected only marginal increases in the separation between the core and the innermost jet component in each image over the 2 year interval 1987-1989. Because monitoring image brightness variations has since led to apparently reliable estimates of delta tau(sub BA), we intend to use spatial variations in the relative magnification field across the extent of the images to constrain further the mass distribution of the lens, which is now the principal contributor to the uncertainty in estimates of H(sub 0) from this lens system.

Campbell, R. M.↗

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base (FY2020 L2 Milestone MRT 7129) [Slides]

Advanced reactive burn models enable phenomena to be modeled that aren’t possible with programmed burn. Programmed burn propagates a burn front at a prescribed speed from a prescribed initiation point. Reactive burn uses the hydro variables (pressure, temperature, etc.) to determine when, how rapidly, and to what extent a zone detonates. This enables a completely new set of problems to be simulated. Our goal in this milestone has been to improve the usability of reactive burn modes. Thus, the closure criteria are: 1. Perform code and solution verification of the reactive burn models using appropriate solutions and simplified HE setups. 2. Perform validation assessments using small-scale experiments. 3. Evaluate sensitivities and quantify uncertainties due to model form variations and mesh dependencies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Retrieval Assimilation and Modeling of Atmospheric Water Vapor from Ground- and Space-Based GPS Networks: Investigation of the Global and Regional Hydrological Cycles

Uncertainty over the response of the atmospheric hydrological cycle (particularly the distribution of water vapor and cloudiness) to anthropogenic forcing is a primary source of doubt in current estimates of global climate sensitivity, which raises severe difficulties in evaluating its likely societal impact. Fortunately, a variety of advanced techniques and sensors are beginning to shed new light on the atmospheric hydrological cycle. One of the most promising makes use of the sensitivity of the Global Positioning System (GPS) to the thermodynamic state, and in particular the water vapor content, of the atmosphere through which the radio signals propagate. Our strategy to derive the maximum benefit for hydrological studies from the rapidly increasing GPS data stream will proceed in three stages: (1) systematically analyze and archive quality-controlled retrievals using state-of-the-art techniques; (2) employ both currently available and innovative assimilation procedures to incorporate these determinations into advanced regional and global atmospheric models and assess their effects; and (3) apply the results to investigate selected scientific issues of relevance to regional and global hydrological studies. An archive of GPS-based estimation of total zenith delay (TZD) data and water vapor where applicable has been established with expanded automated quality control. The accuracy of the GPS estimates is being monitored; the investigation of systematic errors is ongoing using comparisons with water vapor radiometers. Meteorological packages have been implemented. The accuracy and utilization of the TZD estimates has been improved by implementing a troposphere gradient model. GPS-based gradients have been validated as real atmospheric moisture gradients, establishing a link between the estimated gradients and the passage of weather fronts. We have developed a generalized ray tracing inversion scheme that can be used to analyze occultation data acquired from space- or land-based receivers. The National Center for Atmospheric Research mesoscale model (version MM5) has been adapted for Southern California, and assimilation studies are underway. Additional information is contained in the original.

Dickey, Jean O.↗

Statistical Orbit Determination using the Particle Filter for Incorporating Non-Gaussian Uncertainties

The tracking of space objects requires frequent and accurate monitoring for collision avoidance. As even collision events with very low probability are important, accurate prediction of collisions require the representation of the full probability density function (PDF) of the random orbit state. Through representing the full PDF of the orbit state for orbit maintenance and collision avoidance, we can take advantage of the statistical information present in the heavy tailed distributions, more accurately representing the orbit states with low probability. The classical methods of orbit determination (i.e. Kalman Filter and its derivatives) provide state estimates based on only the second moments of the state and measurement errors that are captured by assuming a Gaussian distribution. Although the measurement errors can be accurately assumed to have a Gaussian distribution, errors with a non-Gaussian distribution could arise during propagation between observations. Moreover, unmodeled dynamics in the orbit model could introduce non-Gaussian errors into the process noise. A Particle Filter (PF) is proposed as a nonlinear filtering technique that is capable of propagating and estimating a more complete representation of the state distribution as an accurate approximation of a full PDF. The PF uses Monte Carlo runs to generate particles that approximate the full PDF representation. The PF is applied in the estimation and propagation of a highly eccentric orbit and the results are compared to the Extended Kalman Filter and Splitting Gaussian Mixture algorithms to demonstrate its proficiency.

Mashiku, Alinda↗

Contribution of Asteroid Generated Tsunami to the Impact Hazard

The long-standing uncertainty about the importance of asteroid-generated tsunami was addressed at a workshop in August 2016, co-sponsored by NASA and NOAA. Experts from NASA, NOAA, the DoE tri-labs (LLNL, SNL, and LANL), DHS, FEMA, and academia addressed the hazard of tsunami created by asteroid impacts, focusing primarily on NEAs with diameter less than 250m. Participants jointly identified key issues and shared information for nearly a year to coordinate their results for discussion at the workshop. They used modern computational tools to examine 1) Near-field wave generation by the impact; 2) Long-distance wave propagation; 3) Damage from coastal run-up and inundation, and associated hazard. The workshop resulted in broad consensus that the asteroid impact tsunami threat is not as great as previously thought.

Asteroid↗