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 253 records · Page 14

Temporal Wind Pairs for Space Launch Vehicle Capability Assessment and Risk Mitigation

Space launch vehicles incorporate upper-level wind assessments to determine wind effects on the vehicle and for a commit to launch decision. These assessments make use of wind profiles measured hours prior to launch and may not represent the actual wind the vehicle will fly through. Uncertainty in the winds over the time period between the assessment and launch introduces uncertainty in assessment of vehicle controllability and structural integrity that must be accounted for to ensure launch safety. Temporal wind pairs are used in engineering development of allowances to mitigate uncertainty. Five sets of temporal wind pairs at various times (0.75, 1.5, 2, 3 and 4-hrs) at the United States Air Force Eastern Range and Western Range, as well as the National Aeronautics and Space Administration's Wallops Flight Facility are developed for use in upper-level wind assessments on vehicle performance. Historical databases are compiled from balloon-based and vertically pointing Doppler radar wind profiler systems. Various automated and manual quality control procedures are used to remove unacceptable profiles. Statistical analyses on the resultant wind pairs from each site are performed to determine if the observed extreme wind changes in the sample pairs are representative of extreme temporal wind change. Wind change samples in the Eastern Range and Western Range databases characterize extreme wind change. However, the small sample sizes in the Wallops Flight Facility databases yield low confidence that the sample population characterizes extreme wind change that could occur.

Decker, Ryan K.↗

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Passive Inter-Photon Imaging

Digital camera pixels measure image intensities by converting incident light energy into an analog electrical current, and then digitizing it into a fixed-width binary representation. This direct measurement method, while conceptually simple, suffers from limited dynamic range and poor performance under extreme illumination --- electronic noise dominates under low illumination, and pixel full-well capacity results in saturation under bright illumination. We propose a novel intensity cue based on measuring inter-photon timing, defined as the time delay between detection of successive photons. Based on the statistics of inter-photon times measured by a time-resolved single-photon sensor, we develop theory and algorithms for a scene brightness estimator which works over extreme dynamic range; we experimentally demonstrate imaging scenes with a dynamic range of over ten million to one. The proposed techniques, aided by the emergence of single-photon sensors such as single-photon avalanche diodes (SPADs) with picosecond timing resolution, will have implications for a wide range of imaging applications: robotics, consumer photography, astronomy, microscopy and biomedical imaging.

97 MATHEMATICS AND COMPUTING↗

The hidden structure of hydrodynamic transport in random fracture networks

We study the large-scale dynamics and prediction of hydrodynamic transport in random fracture networks. The flow and transport behaviour is characterized by first passage times and displacement statistics, which show heavy tails and anomalous dispersion with a strong dependence on the injection condition. The origin of these behaviours is investigated in terms of Lagrangian velocities sampled equidistantly along particle trajectories, unlike classical sampling strategies at a constant rate. The velocity series are analysed by their copula density, the joint distribution of the velocity unit scores, which reveals a simple, albeit hidden, correlation structure that can be described by a Gaussian copula. Based on this insight, we derive a Langevin equation for the evolution of equidistant particle speeds. In this framework, particle motion is quantified by a stochastic time-domain random walk, the joint density of particle position, and speed satisfies a Klein–Kramers equation. The upscaled theory quantifies particle motion in terms of the characteristic fracture length scale and the distribution of Eulerian flow velocities. That is, it is predictive in the sense that it does not require the a priori knowledge of transport attributes. The upscaled model captures non-Fickian transport features, and their dependence on the injection conditions in terms of the velocity point statistics and average fracture length. It shows that the first passage times and displacement moments are dominated by extremes occurring at the first step. The presented approach integrates the interaction of flow and structure into a predictive model for large-scale transport in random fracture networks.

42 ENGINEERING↗

Winds aloft statistical analysis in support of day of launch Shuttle systems evaluation

In connection with the development of the Meteorological Interactive Data Display System (MIDDS) for utilization by the Launch Systems Evaluation Advisory Team (LSEAT), requirements have been established to expand the pre-launch analysis of winds aloft for the Space Shuttle. Statistical analyses developed for the system include: comparison of pre-launch wind component profiles to wind component extremes at each altitude calculated from launch site historical data; conditional probability ellipses for wind vectors at a future time given the wind vector at an initial time; comparison of observed extreme wind shear and associated wind speed with launch site historical data utilizing the bivariate extreme value (Gumbel) distribution; estimation of extremes of wind speed or wind shear at a future time given the extremes of either variable at an initial time, utilizing the conditional extreme value distribution; power spectrum analysis for tracking wind perturbation energy in sequential pre-launch Jimsphere wind profiles.

Adelfang, S. I.↗

What the Long-Term Sunspot Record Tells Us About Space Climate

Direct observations of sunspots span the nearly 400 years since the time of Galileo. Dedicated observing programs at several observatories over the last 150 years have provided detailed information not only on the number of sunspots but on their sizes and positions as well. The data acquired by those original observers, and by those who have more recently brought those observations to light, provide important clues about the nature of the solar cycle and its contribution to space climate. The period of the cycle, the equator-ward drift of the active latitudes, the asymmetry between the rise to maximum and the fill to minimum, shifting asymmetries between northern and southern hemisphere activity, the tilt of active regions, and the increasing amplitude of the cycles since the Maunder Minimum are all well established. Other, less well established characteristics such as multi-cycle and short-term periodicities, often depend upon the method of data analysis. The strong correlation between sunspot statistics and other measures of solar activity, coupled with the length of the sunspot record, make these observations extremely valuable for characterizing and understanding space climate.

Hathaway, D. H.↗

A method for predicting failure statistics for steady state elevated temperature structural components

This paper presents the initial development of a high temperature life prediction method that accounts for the variability in the material properties of Grade 91 steel. The method accounts for material variability by fitting a variable 3-parameter Weibull distribution to experimental rupture data and accounts for the variability of creep deformation on the steady-state stresses via a Monte Carlo approach. To ensure reasonable computational times, the model represents the material as an extremely viscous Stokes fluid with a non-Newtonian viscosity, therefore solving the stress relaxation problem with a steady, static, instead of transient, analysis. Furthermore, the complete statistical analysis combines this model for creep deformation with a probabilistic model for creep rupture to evaluate the probability of premature failure for a set of sample problems, comparing the predicted failure statistics to the design life predicted by the ASME Boiler and Pressure Vessel Code rules.

42 ENGINEERING↗

Characteristics of Extreme Auroral Charging Events

The highest level spacecraft charging observed in low Earth orbit (LEO) occurs when spacecraft are exposed to energetic auroral electrons. Since auroral charging has been identified as a mechanism responsible for on-orbit anomalies and even possible satellite failures it is important to consider extreme auroral charging events as design and test environments for spacecraft to be used in high inclination LEO orbits. This paper will report on studies of extreme auroral charging events using data from the SSJ/4 and SSJ/5 precipitating electron and ion sensors on the Defense Meteorology Satellite Program (DMSP) satellites. Early studies of DMSP charging to negative potentials ≥100 V focused on statistics of the electron environment responsible for charging. Later statistical studies of auroral charging have generally focused on solar cycle dependence of charging behavior and magnitude of the maximum potential and duration of the charging events. We extend these studies to focus on more detailed investigations of extreme charging event characteristics that are required to evaluate potential threats to spacecraft systems. A collection of example auroral charging events is assembled from the DMSP data set using the criteria that "extreme auroral charging" is defined as periods with spacecraft negative potentials ≥400 V. Specific characteristics to be treated include (but are not limited to) maximum and mean potentials, time history of spacecraft potentials through the events, total charging duration and the time potentials exceed voltage thresholds, frame charging/discharging rates, and information on geographic and geomagnetic latitudes at which the events are observed. Finally, we will comment on the implications of these studies for potential auroral charging risks to the International Space Station.

Minow, Joseph I.↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

Big Data Analytics for Long-Term Meteorological Observations at Hanford Site

A growing number of physical objects with embedded sensors with typically high volume and frequently updated data sets has accentuated the need to develop methodologies to extract useful information from big data for supporting decision making. This study applies a suite of data analytics and core principles of data science to characterize near real-time meteorological data with a focus on extreme weather events. To highlight the applicability of this work and make it more accessible from a risk management perspective, a foundation for a software platform with an intuitive Graphical User Interface (GUI) was developed to access and analyze data from a decommissioned nuclear production complex operated by the U.S. Department of Energy (DOE, Richland, USA). Exploratory data analysis (EDA), involving classical non-parametric statistics, and machine learning (ML) techniques, were used to develop statistical summaries and learn characteristic features of key weather patterns and signatures. The new approach and GUI provide key insights into using big data and ML to assist site operation related to safety management strategies for extreme weather events. Specifically, this work offers a practical guide to analyzing long-term meteorological data and highlights the integration of ML and classical statistics to applied risk and decision science.

54 ENVIRONMENTAL SCIENCES↗

Statistical analysis of multivariate atmospheric variables

Topics covered include: (1) estimation in discrete multivariate distributions; (2) a procedure to predict cloud cover frequencies in the bivariate case; (3) a program to compute conditional bivariate normal parameters; (4) the transformation of nonnormal multivariate to near-normal; (5) test of fit for the extreme value distribution based upon the generalized minimum chi-square; (6) test of fit for continuous distributions based upon the generalized minimum chi-square; (7) effect of correlated observations on confidence sets based upon chi-square statistics; and (8) generation of random variates from specified distributions.

Tubbs, J. D.↗

Estimating maximum instantaneous distortion from inlet total pressure rms and PSD measurements

An inexpensive method is developed to determine the extreme values of instantaneous inlet distortion. This method also provides insight into the basic mechanics of unsteady inlet flow and the associated engine reaction. The analysis is based on fundamental fluid dynamics and statistical methods to provide an understanding of the turbulent inlet flow and quantitatively relate the rms level and power spectral density (PSD) function of the measured time variant total pressure fluctuations to the strength and size of the low pressure regions. The most probable extreme value of the instantaneous distortion is then synthesized from this information in conjunction with the steady state distortion. Results of the analysis show the extreme values to be dependent upon the steady state distortion, the measured turbulence rms level and PSD function, the time on point, and the engine response characteristics. Analytical projections of instantaneous distortion are presented and compared with data obtained by a conventional, highly time correlated, 40 probe instantaneous pressure measurement system.

Melick, H. C., Jr.↗

Determination of Non-Thermal Velocity Distributions from SERTS Linewidth Observations

Non-thermal velocities obtained from the measurement of coronal Extreme Ultraviolet (EUV) linewidths have been consistently observed in solar EUV spectral observations and have been theorized to result fro m many plausible scenarios including wave motions, turbulence, or magnetic reconnection. Constraining these velocities can provide a physical limit for the available energy resulting from unresolved motions in the corona. We statistically determine a series of non-thermal velocity distributions from linewidth measurements of 390 emission lines from a wide array of elements and ionization states observed during the Solar Extreme Ultraviolet Research Telescope and Spectrograph 1991-1997 flights covering the spectral range 174-418 A and a temperature range from 80,000 K to 12.6 MK. This sample includes 248 lines from active regions, 101 lines from quiet-Sun regions, and 41 lines were observed from plasma off the solar limb. We find a strongly peaked distribution corresponding to a non-thermal velocity of 19-22 km/s in all three of the quiet-Sun, active region, and off-limb distributions. For the possibility of Alfven wave resonance heating, we find thai velocities in the core of these distributions do not provide sufficient energy, given typical densities and magnetic field strengths for the coronal plasma, to overcome the estimated coronal energy losses required to maintain the corona at the typical temperatures working as the sole mechanism. We find that at perfect efficiency 50%-60% of the needed energy flux can be produced from the non-thermal velocities measured.

Coyner, Aaron J.↗

Decadal predictability of North Atlantic blocking and the NAO

Can multi-annual variations in the frequency of North Atlantic atmospheric blocking and mid-latitude circulation regimes be skilfully predicted? Recent advances in seasonal forecasting have shown that mid-latitude climate variability does exhibit significant predictability. However, atmospheric predictability has generally been found to be quite limited on multi-annual timescales. New decadal prediction experiments from NCAR are found to exhibit remarkable skill in reproducing the observed multi-annual variations of wintertime blocking frequency over the North Atlantic and of the North Atlantic Oscillation (NAO) itself. This is partly due to the large ensemble size that allows the predictable component of the atmospheric variability to emerge from the background chaotic component. The predictable atmospheric anomalies represent a forced response to oceanic low-frequency variability that strongly resembles the Atlantic Multi-decadal Variability (AMV), correctly reproduced in the decadal hindcasts thanks to realistic ocean initialization and ocean dynamics. The occurrence of blocking in certain areas of the Euro-Atlantic domain determines the concurrent circulation regime and the phase of known teleconnections, such as the NAO, consequently affecting the stormtrack and the frequency and intensity of extreme weather events. Therefore, skilfully predicting the decadal fluctuations of blocking frequency and the NAO may be used in statistical predictions of near-term climate anomalies, and it provides a strong indication that impactful climate anomalies may also be predictable with improved dynamical models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semi-automatic crop inventory from sequential ERTS-1 imagery

The detection of a newly introduced crop into the Imperial (California) Valley by sequential ERTS-1 imagery is proving that individual crop types can be identified by remote sensing techniques. Initial results have provided an extremely useful product for water agencies. A system for the identification of field conditions enables the production of a statistical summary within two to three days of receipt of the ERTS-1 imagery. The summary indicates the total acreage of producing crops and irrigated planted crops currently demanding water and further indicates freshly plowed fields that will be demanding water in the near future. Relating the field conditions to the crop calendar of the region by means of computer techniques will provide specific crop identification for the 8000 plus fields.

Johnson, C. W.↗

NASA Tech Briefs, June 2014

Topics include: Real-Time Minimization of Tracking Error for Aircraft Systems; Detecting an Extreme Minority Class in Hyperspectral Data Using Machine Learning; KSC Spaceport Weather Data Archive; Visualizing Acquisition, Processing, and Network Statistics Through Database Queries; Simulating Data Flow via Multiple Secure Connections; Systems and Services for Near-Real-Time Web Access to NPP Data; CCSDS Telemetry Decoder VHDL Core; Thermal Response of a High-Power Switch to Short Pulses; Solar Panel and System Design to Reduce Heating and Optimize Corridors for Lower-Risk Planetary Aerobraking; Low-Cost, Very Large Diamond-Turned Metal Mirror; Very-High-Load-Capacity Air Bearing Spindle for Large Diamond Turning Machines; Elevated-Temperature, Highly Emissive Coating for Energy Dissipation of Large Surfaces; Catalyst for Treatment and Control of Post-Combustion Emissions; Thermally Activated Crack Healing Mechanism for Metallic Materials; Subsurface Imaging of Nanocomposites; Self-Healing Glass Sealants for Solid Oxide Fuel Cells and Electrolyzer Cells; Micromachined Thermopile Arrays with Novel Thermo - electric Materials; Low-Cost, High-Performance MMOD Shielding; Head-Mounted Display Latency Measurement Rig; Workspace-Safe Operation of a Force- or Impedance-Controlled Robot; Cryogenic Mixing Pump with No Moving Parts; Seal Design Feature for Redundancy Verification; Dexterous Humanoid Robot; Tethered Vehicle Control and Tracking System; Lunar Organic Waste Reformer; Digital Laser Frequency Stabilization via Cavity Locking Employing Low-Frequency Direct Modulation; Deep UV Discharge Lamps in Capillary Quartz Tubes with Light Output Coupled to an Optical Fiber; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems, Version II; Advanced Sensor Technology for Algal Biotechnology; High-Speed Spectral Mapper; "Ascent - Commemorating Shuttle" - A NASA Film and Multimedia Project DVD; High-Pressure, Reduced-Kinetics Mechanism for N-Hexadecane Oxidation; Method of Error Floor Mitigation in Low-Density Parity-Check Codes; X-Ray Flaw Size Parameter for POD Studies; Large Eddy Simulation Composition Equations for Two-Phase Fully Multicomponent Turbulent Flows; Scheduling Targeted and Mapping Observations with State, Resource, and Timing Constraints;

Source record↗

First detection of VHE gamma-ray emission from TXS 1515–273, study of its X-ray variability and spectral energy distribution

ABSTRACT We report here on the first multiwavelength (MWL) campaign on the blazar TXS 1515–273, undertaken in 2019 and extending from radio to very-high-energy gamma-rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotron-peaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6σ. The average integral VHE flux of the source is 6 ± 1 per cent of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMM–Newton, and NuSTAR. The long continuous X-ray observations were separated by ∼9 h, both showing clear hour scale flares. In the XMM–Newton data, both the rise and decay time-scales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability time-scales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme synchrotron-peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band.

Acciari, V. A.↗

Investigating the potential of a global precipitation forecast to inform landslide prediction

Extreme rainfall events within landslide-prone areas can be catastrophic, resulting in loss of property, infrastructure, and life. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model provides routine near-real time estimates of landslide hazard using Integrated Multi-Satellite Precipitation Retrievals for the Global Precipitation Mission (IMERG). However, it does not provide information on potential landslide hazard in the future. Forecasting potential landslide events at a global scale presents an area of open research. This study compares a global precipitation forecast provided by NASA's Goddard Earth Observing System (GEOS) with near-real time satellite precipitation estimates. The Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) reference is used to assess the performance of both satellite and model-based precipitation products over the contiguous United States (CONUS). The forecast lead time of 24hrs is considered, with a focus on extreme precipitation events. The performance of IMERG and GEOS-Forecast products is assessed in terms of the probability of detection, success ratio, critical success index and hit bias as well as continuous statistics. The results show that seasonality influences the performance of both satellite and model-based precipitation products. Comparison of IMERG and GEOS-Forecast globally as well as in several event case studies (Colombia, southeast Asia, and Tajikistan) reveals that GEOS-Forecast detects extreme rainfall more frequently relative to IMERG for these specific analyses. For recent landslide points across the globe, the 24hr accumulated precipitation forecast >100 mm corresponds well with near-real time daily accumulated IMERG precipitation estimates. GEOS-Forecast and IMERG precipitation match more closely for tropical cyclones than for other types of storms. The main intention of this study is to assess the viability of using a global forecast for landslide predictions and understand the extent of the variability between these products to inform where we would expect the landslide modeling results to most prominently diverge. Results of this study will be used to inform how forecasted precipitation estimates can be incorporated into the LHASA model to provide the first global predictive view of landslide hazards.

S. Khan↗