Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistics of extremes”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Cloud and Radiation Studies during SAFARI 2000

Though the emphasis of the Southern Africa Regional Science Initiative 2000 (SAFARI-2000) dry season campaign was largely on emission sources and transport, the assemblage of aircraft (including the high altitude NASA ER-2 remote sensing platform and the University of Washington CV-580, UK MRF C130, and South African Weather Bureau JRA in situ aircrafts) provided a unique opportunity for cloud studies. Therefore, as part of the SAFARI initiative, investigations were undertaken to assess regional aerosol-cloud interactions and cloud remote sensing algorithms. In particular, the latter part of the experiment concentrated on marine boundary layer stratocumulus clouds off the southwest coast of Africa. Associated with cold water upwelling along the Benguela current, the Namibian stratocumulus regime has received limited attention but appears to be unique for several reasons. During the dry season, outflow of continental fires and industrial pollution over this area can be extreme. From below, upwelling provides a rich nutrient source for phytoplankton (a source of atmospheric sulphur through DMS production as well as from decay processes). The impact of these natural and anthropogenic sources on the microphysical and optical properties of the stratocumulus is unknown. Continental and Indian Ocean cloud systems of opportunity were also studied during the campaign. Aircraft flights were coordinated with NASA Terra Satellite overpasses for synergy with the Moderate Resolution Imaging Spectroradiometer (MODIS) and other Terra instruments. An operational MODIS algorithm for the retrieval of cloud optical and physical properties (including optical thickness, effective particle radius, and water path) has been developed. Pixel-level MODIS retrievals (11 km spatial resolution at nadir) and gridded statistics of clouds in th SAFARI region will be presented. In addition, the MODIS Airborne Simulator flown on the ER-2 provided high spatial resolution retrievals (50 m at nadir). These retrievals will be discussed and compared with in situ observations.

Platnick, Steven↗
DOE OSTI2024

Topological Characterization and Uncertainty Visualization of Atmospheric Rivers

Atmospheric rivers (ARs) are long, narrow regions of water vapor in the Earth's atmosphere that transport heat and moisture from the tropics to the mid-latitudes. ARs are often associated with extreme weather events in North America and contribute significantly to water supply and flood risk. However, characterizing ARs has been a major challenge due to the lack of a universal definition and their structural variations. Existing AR detection tools (ARDTs) produce distinct AR boundaries for the same event, making the risk assessment of ARs a difficult task. Understanding these uncertainties is crucial to improving the predictability of AR impacts, including their landfall areas and associated precipitation, which could cause catastrophic flooding and landslides over the coastal regions. In this work, we develop an uncertainty visualization framework that captures boundary and interior uncertainties, i.e., structural variations, of an ensemble of ARs that arise from a set of ARDTs. We first provide a statistical overview of the AR boundaries using the contour boxplots of Whitaker et al. that highlight the structural variations of AR boundaries based on their nesting relationships. We then introduce the topological skeletons of ARs based on Morse complexes that characterize the interior variation of an ensemble of ARs. We propose an uncertainty visualization of these topological skeletons, inspired by MetroSets of Jacobson et al. that emphasizes the agreements and disagreements across the ensemble members. Through case studies and expert feedback, we demonstrate that the two approaches complement each other, and together they could facilitate an effective comparative analysis process and provide a more confident outlook on an AR's shape, area, and onshore impact.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Guiding the choice of informatics software and tools for lipidomics research applications

Progress in mass spectrometry lipidomics has led to a rapid proliferation of studies across biology and biomedicine. These generate extremely large raw datasets requiring sophisticated solutions to support automated data processing. To address this, numerous software tools have been developed and tailored for specific tasks. However, for researchers, deciding which approach best suits their application relies on ad hoc testing, which is inefficient and time consuming. Here we first review the data processing pipeline, summarizing the scope of available tools. Next, to support researchers, LIPID MAPS provides an interactive online portal listing open-access tools with a graphical user interface. This guides users towards appropriate solutions within major areas in data processing, including (1) lipid-oriented databases, (2) mass spectrometry data repositories, (3) analysis of targeted lipidomics datasets, (4) lipid identification and (5) quantification from untargeted lipidomics datasets, (6) statistical analysis and visualization, and (7) data integration solutions. Detailed descriptions of functions and requirements are provided to guide customized data analysis workflows.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2010

MG Isotopic Measurement of FIB-Isolated Presolar Silicate Grains

The majority of presolar oxide and silicate grains are ascribed to origins in low-mass red giant and asymptotic giant branch (AGB) stars based on their O isotopic ratios. However, a minor population of these grains (< 10%) has O isotopic ratios incompatible with these sources. Two principle alternative sources are higher-than-solar metallicity (Z) stars or, more likely, supernovae (SN) [1-3]. These rare (Group 4) grains [3] are characterized by enrichments in O-18, and typically also enrichments in O-17. An even rarer subset of grains with extremely large enrichments in O-17 and smaller depletions in O-18 were suggested to come from binary star systems [2]. To establish the origins of these isotopically unusual grains, it is necessary to examine isotopic systems in addition to O. Presolar silicates offer several elements diagnostic of their stellar sources and nuclear processes, including O, Si, Mg, Fe and Ca. However, the database for minor element isotopic compositions in silicates is seriously lacking. To date only two silicate grains have been analyzed for Mg [4] or Fe [5]. One major complicating factor is their small size (average ~230 nm), which greatly limits the number of measurements that can be performed on any one grain and makes it more difficult to obtain statistically relevant data. This problem is compounded because the grains are identified among isotopically solar silicates, which contribute a diluting signal in isotopic measurements [1]. Thus, relatively small isotopic anomalies are missed due to this dilution effect. By applying focused ion beam (FIB) milling, we obtain undiluted Mg isotopic ratios of isolated rare presolar silicate grains to investigate their sources.

Messenger, Scott R.↗
NASA NTRS2012

Uncertainty Quantification for Polynomial Systems via Bernstein Expansions

This paper presents a unifying framework to uncertainty quantification for systems having polynomial response metrics that depend on both aleatory and epistemic uncertainties. The approach proposed, which is based on the Bernstein expansions of polynomials, enables bounding the range of moments and failure probabilities of response metrics as well as finding supersets of the extreme epistemic realizations where the limits of such ranges occur. These bounds and supersets, whose analytical structure renders them free of approximation error, can be made arbitrarily tight with additional computational effort. Furthermore, this framework enables determining the importance of particular uncertain parameters according to the extent to which they affect the first two moments of response metrics and failure probabilities. This analysis enables determining the parameters that should be considered uncertain as well as those that can be assumed to be constants without incurring significant error. The analytical nature of the approach eliminates the numerical error that characterizes the sampling-based techniques commonly used to propagate aleatory uncertainties as well as the possibility of under predicting the range of the statistic of interest that may result from searching for the best- and worstcase epistemic values via nonlinear optimization or sampling.

Crespo, Luis G.↗
DOE OSTI2021

Proactive Operations and Investment Planning via Stochastic Optimization to Enhance Power Systems’ Extreme Weather Resilience

We present scalable stochastic optimization approaches for improving power systems’ resilience to extreme weather events. We consider both proactive redispatch and transmission line hardening as alternatives for mitigating expected load shed due to extreme weather, resulting in large-scale stochastic linear programs (LPs) and mixed-integer linear programs (MILPs). We solve these stochastic optimization problems with progressive hedging (PH), a parallel, scenario-based decomposition algorithm. Our computational experiments indicate that our proposed method for enhancing power system resilience can provide high-quality solutions efficiently. With up to 128 scenarios on a 2,000-bus network, the operations (redispatch) and investment (hardening) resilience problems can be solved in approximately 6 min and 2 h of wall-clock time, respectively. Additionally, we solve the investment problems with up to 512 scenarios, demonstrating that the approach scales very well with the number of scenarios. Moreover, the method produces high quality solutions that result in statistically significant reductions in expected load shed. Our proposed approach can be augmented to incorporate a variety of other operational and investment resilience strategies, or a combination of such strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2003

Software Reliability 2002

In FY01 we learned that hardware reliability models need substantial changes to account for differences in software, thus making software reliability measurements more effective, accurate, and easier to apply. These reliability models are generally based on familiar distributions or parametric methods. An obvious question is 'What new statistical and probability models can be developed using non-parametric and distribution-free methods instead of the traditional parametric method?" Two approaches to software reliability engineering appear somewhat promising. The first study, begin in FY01, is based in hardware reliability, a very well established science that has many aspects that can be applied to software. This research effort has investigated mathematical aspects of hardware reliability and has identified those applicable to software. Currently the research effort is applying and testing these approaches to software reliability measurement, These parametric models require much project data that may be difficult to apply and interpret. Projects at GSFC are often complex in both technology and schedules. Assessing and estimating reliability of the final system is extremely difficult when various subsystems are tested and completed long before others. Parametric and distribution free techniques may offer a new and accurate way of modeling failure time and other project data to provide earlier and more accurate estimates of system reliability.

Wallace, Dolores R.↗
DOE OSTI2020

Power and Wind Shear Implications of Large Wind Turbine Scenarios in the US Central Plains

Continued growth of wind turbine physical dimensions is examined in terms of the implications for wind speed, power and shear across the rotor plane. High-resolution simulations with the Weather Research and Forecasting model are used to generate statistics of wind speed profiles for scenarios of current and future wind turbines. The nine-month simulations, focused on the eastern Central Plains, show that the power scales broadly as expected with the increase in rotor diameter (D) and wind speeds at hub-height (H). Increasing wind turbine dimensions from current values (approximately H = 100 m, D = 100 m) to those of the new International Energy Agency reference wind turbine (H = 150 m, D = 240 m), the power across the rotor plane increases 7.1 times. The mean domain-wide wind shear exponent (α) decreases from 0.21 (H = 100 m, D = 100 m) to 0.19 for the largest wind turbine scenario considered (H = 168 m, D = 248 m) and the frequency of extreme positive shear (α > 0.2) declines from 48% to 38% of 10-min periods. Thus, deployment of larger wind turbines potentially yields considerable net benefits for both the wind resource and reductions in fatigue loading related to vertical shear.

17 WIND ENERGY↗
DOE OSTI2025

A Unified Model of Cosmic-Ray Propagation and Radio Extreme Scattering Events from Intermittent Interstellar Structures

Intermittent magnetic structures are a plausible candidate for explaining cosmic-ray (CR) diffusion rates derived from observed CR energy spectra. Independently, studies of extreme scattering events (ESEs) of radio quasars and pulsar scintillation have hinted that very straight, large aspect ratio magnetic current sheets may be responsible for the localized large scattering of radio waves. The required shortest axis of the typical structures producing ESEs is of the same scale (∼au) as the gyroradii of ∼GeV CRs. In this Letter, we propose that the same magnetic/density sheets can produce large scattering of both CRs and radio waves. We demonstrate that the geometry and volume-filling factor of the sheets derived from quasar ESEs can explain the observed mean free path of GeV CRs without introducing free parameters. The model places constraints on the sheet geometry, such as straightness and large aspect ratio, and assumes that the statistics of the sheets are similar throughout the Galactic volume. We therefore discuss observational tests of the sheet model, which includes observations of echoes in pulsars and fast radio bursts, gravitationally lensed quasars, the distribution of ESE durations, and spatial correlations between ESE events and rotation measure fluctuations. Such tests will be enabled by upcoming wide-field radio instruments, including the Canadian Hydrogen Observatory and Radio-transient Detector and Deep Synoptic Array 2000 Antennas.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗
NASA NTRS1994

The complex X-ray spectra of two high redshift quasars observed with ASCA

We report preliminary results of Advanced Satellite for Cosmology and Astrophysics (ASCA) observations of two high-redshift quasars (zeta approximately 3) PKS 0438-436 and PKS 2126-158. We find spectral flattening towards low energies in both sources, first found with Roentgen Satellite (ROSAT) and interpreted as being due to excess X-ray absorption. However, the ROSAT data lacked the bandpass and sensitivity to unambiguously support this interpretation. We find that for PKS 2126-158, absorption in intervening low-redshift material along the line-of-sight provides a statistically better description of the data than absorption in material intrinsic to the quasar. The redshift of the putative absorber is required to be less than a few tenths. Given the redshift constraints and the fact that absorption is not common in lower redshift quasars, this result is not easily explained by a simple model. Either the absorber may be complex or else the intrinsic X-ray spectrum may be complex, or both. Compton-reflection dominated models are ruled out but a broken power-law model provides good fits to the data without requiring excess absorption. The break energy is at approximately 6-8 keV in the quasar frame and the energy index of the soft component is required to be extremely flat. This may mean that the X-ray emission in these objects is dominated by synchrotron losses as in BL Lac objects.

Serlemitsos, Peter↗
NASA NTRS2011

Fermi and Swift Gamma-Ray Burst Afterglow Population Studies

The new and extreme population of GRBs detected by Fermi -LAT shows several new features in high energy gamma-rays that are providing interesting and unexpected clues into GRB prompt and afterglow emission mechanisms. Over the last 6 years, it has been Swift that has provided the robust dataset of UV/optical and X-ray afterglow observations that opened many windows into components of GRB emission structure. The relationship between the LAT detected GRBs and the well studied, fainter, less energetic GRBs detected by Swift -BAT is only beginning to be explored by multi-wavelength studies. We explore the large sample of GRBs detected by BAT only, BAT and Fermi -GBM, and GBM and LAT, focusing on these samples separately in order to search for statistically significant differences between the populations, using only those GRBs with measured redshifts in order to physically characterize these objects. We disentangle which differences are instrumental selection effects versus intrinsic properties, in order to better understand the nature of the special characteristics of the LAT bursts.

Racusin, Judith L.↗
NASA NTRS2011

Fermi and Swift Gamma-Ray Burst Afterglow Population Studies

The new and extreme population of GRBs detected by Fermi-LAT shows several new features in high energy gamma-rays that are providing interesting and unexpected clues into GRB prompt and afterglow emission mechanisms. Over the last 6 years, it has been Swift that has provided the robust data set of UV/optical and X-ray afterglow observations that opened many windows into components of GRB emission structure. The relationship between the LAT GRBs and the well studied, fainter, less energetic GRBs detected by Swift-BAT is only beginning to be explored by multiwavelength studies. We explore the large sample of GRBs detected by BAT only, BAT and Fermi-GBM, and GBM and LAT, focusing on these samples separately in order to search for statistically significant differences between the populations, using only those GRBs with measured redshifts in order to physically characterize these objects. We disentangle which differences are instrumental selection effects versus intrinsic properties, in order to better understand the nature of the special characteristics of the LAT bursts.

Racusin, J. L.↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Calomplification — the power of generative calorimeter models

Motivated by the high computational costs of classical simulations, machine-learned generative models can be extremely useful in particle physics and elsewhere. They become especially attractive when surrogate models can efficiently learn the underlying distribution, such that a generated sample outperforms a training sample of limited size. This kind of GANplification has been observed for simple Gaussian models. We show the same effect for a physics simulation, specifically photon showers in an electromagnetic calorimeter.

47 OTHER INSTRUMENTATION↗
NASA NTRS1991

Development of a calibrated software reliability model for flight and supporting ground software for avionic systems

The object of this project was to develop and calibrate quantitative models for predicting the quality of software. Reliable flight and supporting ground software is a highly important factor in the successful operation of the space shuttle program. The models used in the present study consisted of SMERFS (Statistical Modeling and Estimation of Reliability Functions for Software). There are ten models in SMERFS. For a first run, the results obtained in modeling the cumulative number of failures versus execution time showed fairly good results for our data. Plots of cumulative software failures versus calendar weeks were made and the model results were compared with the historical data on the same graph. If the model agrees with actual historical behavior for a set of data then there is confidence in future predictions for this data. Considering the quality of the data, the models have given some significant results, even at this early stage. With better care in data collection, data analysis, recording of the fixing of failures and CPU execution times, the models should prove extremely helpful in making predictions regarding the future pattern of failures, including an estimate of the number of errors remaining in the software and the additional testing time required for the software quality to reach acceptable levels. It appears that there is no one 'best' model for all cases. It is for this reason that the aim of this project was to test several models. One of the recommendations resulting from this study is that great care must be taken in the collection of data. When using a model, the data should satisfy the model assumptions.

Lawrence, Stella↗
NASA NTRS2012

Coordinated STEM/FIB/NanoSIMS Analyses of Presolar Silicates in Comet Dust and Primitive Meteorites

Silicate grains were among the most abundant mineralogical building blocks of our Solar System. These grains were the detritus from earlier generations of stars that have been recycled in the early solar nebula. Rare sub-micrometer survivors of this processing have been identified in meteorites, micrometeorites and interplanetary dust particles (IDPs). These silicate grains are recognized as presolar in origin because of their extremely anomalous isotopic compositions that reflect nucleosynthetic processes in their stellar sources (evolved stars, novae and supernovae). We perform coordinated chemical, mineralogical and isotopic studies of these grains to determine their origins and histories. We examine the complex mineralogy and petrography of presolar silicates using imaging, diffraction and chemical data obtained from thin sections with the JSC JEOL 2500 field-emission STEM equipped with a Noran thin window energy dispersive x-ray (EDX) spectrometer and a Gatan Tridiem GIF. Quantitative element x-ray maps (spectrum images) are acquired by rastering a 4 nm incident probe whose dwell time is minimized to avoid beam damage and element diffusion during mapping. Successive image layers are acquired and combined in order to achieve approx 1% counting statistics for major elements. The IDP samples are prepared by ultramicrotomy of particles embedded in epoxy or elemental sulfur. After EDX mapping, the sections are subjected to C, N, and O isotopic imaging with the JSC NanoSIMS 50L ion microprobe. We prepare sections of some meteorite grains using the JSC FEI Quanta 3D focused ion beam (FIB) instrument. The specimen surface is protected from the FIB milling process by layers of electron beam-deposited C and Pt followed by an ion-deposited Pt layer. We also use the FIB to preferentially remove surrounding grains to reduce the background in subsequent NanoSIMS measurements. For mineralogical studies, we again employ the FIB instrument to deposit a protective cap over the grain of interest and then extract the grain and thin it to electron transparency for TEM analysis.

Keller, Lindsay↗
NASA NTRS2008

On Certain New Methodology for Reducing Sensor and Readout Electronics Circuitry Noise in Digital Domain

NASA Hubble Space Telescope (HST) and upcoming cosmology science missions carry instruments with multiple focal planes populated with many large sensor detector arrays. These sensors are passively cooled to low temperatures for low-level light (L3) and near-infrared (NIR) signal detection, and the sensor readout electronics circuitry must perform at extremely low noise levels to enable new required science measurements. Because we are at the technological edge of enhanced performance for sensors and readout electronics circuitry, as determined by thermal noise level at given temperature in analog domain, we must find new ways of further compensating for the noise in the signal digital domain. To facilitate this new approach, state-of-the-art sensors are augmented at their array hardware boundaries by non-illuminated reference pixels, which can be used to reduce noise attributed to sensors. There are a few proposed methodologies of processing in the digital domain the information carried by reference pixels, as employed by the Hubble Space Telescope and the James Webb Space Telescope Projects. These methods involve using spatial and temporal statistical parameters derived from boundary reference pixel information to enhance the active (non-reference) pixel signals. To make a step beyond this heritage methodology, we apply the NASA-developed technology known as the Hilbert- Huang Transform Data Processing System (HHT-DPS) for reference pixel information processing and its utilization in reconfigurable hardware on-board a spaceflight instrument or post-processing on the ground. The methodology examines signal processing for a 2-D domain, in which high-variance components of the thermal noise are carried by both active and reference pixels, similar to that in processing of low-voltage differential signals and subtraction of a single analog reference pixel from all active pixels on the sensor. Heritage methods using the aforementioned statistical parameters in the digital domain (such as statistical averaging of the reference pixels themselves) zeroes out the high-variance components, and the counterpart components in the active pixels remain uncorrected. This paper describes how the new methodology was demonstrated through analysis of fast-varying noise components using the Hilbert-Huang Transform Data Processing System tool (HHT-DPS) developed at NASA and the high-level programming language MATLAB (Trademark of MathWorks Inc.), as well as alternative methods for correcting for the high-variance noise component, using an HgCdTe sensor data. The NASA Hubble Space Telescope data post-processing, as well as future deep-space cosmology projects on-board instrument data processing from all the sensor channels, would benefit from this effort.

Kizhner, Semion↗
NASA NTRS2011

Fermi and Swift Gamma-Ray Burst Afterglow Population Studies

The new and extreme population of GRBs detected by Fermi-LAT shows several new features in high energy gamma-rays that are providing interesting and unexpected clues into GRB prompt and afterglow emission mechanisms. Over the last 6 years, it has been Swift that has provided the robust data set of UV/optical and X-ray afterglow observations that opened many windows into components of GRB emission structure. The relationship between the LAT GRBs and the well studied, fainter, less energetic GRBs detected by Swift-BAT is only beginning to be explored by multi-wavelength studies. We explore the large sample of GRBs detected by BAT only, BAT and Fermi-GBM, and GBM and LAT, focusing on these samples separately in order to search for statistically significant differences between the populations, using only those GRBs with measured redshifts in order to physically characterize these objects. We disentangle which differences are instrumental selection effects versus intrinsic properties, in order to better understand the nature of the special characteristics of the LAT bursts.

Racusin, Judith I.↗