Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

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 865 records · Page 48

Overview of the Digitization Workflow Post Image Acquisition of Apollo Lunar and Antarctic Meteorite Samples Using Agisoft Photoscan for the NASA 3D Astromaterials Virtual Samples Collection

The 3D Virtual Astromaterials Samples (3DVAS) collection is a multi-year funded project to create a digital database of sixty Apollo Lunar and Antarctic Meteorite samples following non-destructive documentation conservation protocols. After initial image processing, the photos are evaluated and processed using unique structure-from-motion photogrammetric techniques in a high performance modelling software designed to create a 3D model from 2D images: Agisoft Photoscan Pro. Agisoft Photoscan Pro uses image processing algorithms and techniques originating in computer vision to resolve 3D models for accurate and detailed visualization of a subject. The software provides a stepwise process that is tailored per model based on spatial and specular reflectance properties, for example. The process includes: photo alignment, creation of a dense point cloud, mesh, and finally texture. Photo alignment is dependent on model properties. The 3DVAS process requires a special rotation platform with calibrated photogrammetric targets, specific distance rotation protocols, and a contrasting background for alignment and scale accuracy. As a result of the photographic process, alignment will complete with two mirrored hemispheres that, in a sense, represent the 2D images overlapping to create a 3D model. Each dense point cloud is analyzed with provided statistical measures in a gradual selection process to eliminate outliers. The point cloud is reduced to include only data valuable to the final model. When a precise dense point cloud is achieved, a mesh and texture are applied. Each model is scaled with scale bar accuracies within 100 microns. Each sample has its own intimate process for modelling; there is no standard for the parameters required in the final creation of a high resolution model. By processing multiple samples, a skill is gained in practice to allow a close definition of the original sample and will result in the most detailed version of the sample shell. This process completes one-fifth of the 3DVAS protocol for providing accurate digital documentation. Each model shell is merged with X-ray Computed Tomography data to create a full volumetric sample. All 3DVAS data will be served on NASA's Astromaterials Acquisition and Curation website with an early subset of data available in 2019 and the 3D Virtual Astromaterials Samples Collection launch in 2020.

Thomas, Andi B.↗

Statistical analysis of trigonometric parallaxes

The present project was initiated with two specific medium-term goals: first, to develop a novel approach, based on global modeling and maximum likelihood, to the study of databases of stellar data, with specific attention to the results to be obtained by the Hipparcos mission; and second, to apply and test the methodology developed to existing ground-based data. The long-term goal was of course to integrate the methodology and the ground-based data into a global study of the results of the Hipparcos mission, which were expected to be made available in 1995 or 1996. The closing of the project is due to the recent change of home institution of the PI, who has moved from the University of Illinois to Johns Hopkins University. As a consequence of this move, the three scientists involved in this project will be in close proximity, thereby hopefully improving their ability to collaborate and increasing the productivity of the project. Continued funding for the project has been requested from Johns Hopkins University. This report describes briefly the results obtained so far both on the technical aspect of software development and on the scientific side of applications to existing ground-based data. Significant progress has been made on both counts, with several papers published in (or submitted to) refereed journals and in conference proceedings. Because it is hoped that the project can be successfully continued with NASA support, the report of the progress in each area includes also an assessment of how the current results fit in the expected continuation of the project. Our results to date include: code development (essentially completed); a detailed study of the kinematics and dynamics of stars escaping from the Hyades cluster, relevant to the question of membership in the Hyades; a study of the kinematics and luminosity calibration of nearby dwarfs; an assessment of the quality of the photometry included in the Hipparcos Input Catalog; and two studies of properties of nearby clusters, including a moving-cluster determination of the distance to Praesepe. The bibliography includes three papers submitted to refereed journals, two of which have already been published, and four contributions to conference proceedings. Finally, the work so far has also provided a very good introduction to stellar dynamics and astrometry for an undergraduate student, with educational benefits that had not been foreseen in the original proposal.

Casertano, Stefano↗

Dual adaptive control: Design principles and applications

The design of an actively adaptive dual controller based on an approximation of the stochastic dynamic programming equation for a multi-step horizon is presented. A dual controller that can enhance identification of the system while controlling it at the same time is derived for multi-dimensional problems. This dual controller uses sensitivity functions of the expected future cost with respect to the parameter uncertainties. A passively adaptive cautious controller and the actively adaptive dual controller are examined. In many instances, the cautious controller is seen to turn off while the latter avoids the turn-off of the control and the slow convergence of the parameter estimates, characteristic of the cautious controller. The algorithms have been applied to a multi-variable static model which represents a simplified linear version of the relationship between the vibration output and the higher harmonic control input for a helicopter. Monte Carlo comparisons based on parametric and nonparametric statistical analysis indicate the superiority of the dual controller over the baseline controller.

Mookerjee, Purusottam↗

A method of measuring fault latency in a digital flight control system

This paper describes the motivation, conduct, and analysis of some 2500 low-level hardware fault cases applied in automated testing at the NASA Ames Reconfigurable Digital Flight Control System Facility. Fault detection was correlated with hardware and software fault monitoring and, in limited cases, with sensitivity to flight program execution modes. The results are statistically assessed to ascertain system-level reliability implications based on a single-fault model. Extension to multiple-fault models is addressed. The overall methodology/facility itself is judged to be a promising enhancement to current practice.

Mcgough, John↗

Autoclass: An automatic classification system

The task of inferring a set of classes and class descriptions most likely to explain a given data set can be placed on a firm theoretical foundation using Bayesian statistics. Within this framework, and using various mathematical and algorithmic approximations, the AutoClass System searches for the most probable classifications, automatically choosing the number of classes and complexity of class descriptions. A simpler version of AutoClass has been applied to many large real data sets, has discovered new independently-verified phenomena, and has been released as a robust software package. Recent extensions allow attributes to be selectively correlated within particular classes, and allow classes to inherit, or share, model parameters through a class hierarchy. The mathematical foundations of AutoClass are summarized.

Stutz, John↗

Bayesian classification theory

The task of inferring a set of classes and class descriptions most likely to explain a given data set can be placed on a firm theoretical foundation using Bayesian statistics. Within this framework and using various mathematical and algorithmic approximations, the AutoClass system searches for the most probable classifications, automatically choosing the number of classes and complexity of class descriptions. A simpler version of AutoClass has been applied to many large real data sets, has discovered new independently-verified phenomena, and has been released as a robust software package. Recent extensions allow attributes to be selectively correlated within particular classes, and allow classes to inherit or share model parameters though a class hierarchy. We summarize the mathematical foundations of AutoClass.

Hanson, Robin↗

Recent Advances in Statistical Energy Analysis

Statistical Energy Analysis (SEA) has traditionally been developed using modal summation and averaging approach, and has led to the need for many restrictive SEA assumptions. The assumption of 'weak coupling' is particularly unacceptable when attempts are made to apply SEA to structural coupling. It is now believed that this assumption is more a function of the modal formulation rather than a necessary formulation of SEA. The present analysis ignores this restriction and describes a wave approach to the calculation of plate-plate coupling loss factors. Predictions based on this method are compared with results obtained from experiments using point excitation on one side of an irregular six-sided box structure. Conclusions show that the use and calculation of infinite transmission coefficients is the way forward for the development of a purely predictive SEA code.

K H Heron↗

The use of the AOA TTR-4P GPS receiver in operation at the BIPM for real-time restitution of GPS time

The Global Positioning System is an outstanding tool for the dissemination of time. Using mono-channel C/A-code GPS time receivers, the restitution of GPS time through the satellite constellation presents a peak-to-peak discrepancy of several tens of nanoseconds without SA but may be as high as several hundreds of nanoseconds with SA. As a consequence, civil users are more and more interested in implementing hardware and software methods for efficient restitution of GPS time, especially in the framework of the project of a real-time prediction of UTC (UTCp) which could be available in the form of time differences (UTCp - GPS time). Previous work, for improving the real-time restitution of GPS time with SA, to the level obtained without SA, focused on the implementation of a Kalman filter based on past data and updated at each new observation. An alternative solution relies upon the statistical features of the noise brought about by SA; it has already been shown that the SA noise is efficiently reduced by averaging data from numerous satellites observed simultaneously over a sufficiently long time. This method was successfully applied to data from a GPS time receiver, model AOA TTR-4P, connected to the cesium clock kept at the BIPM. This device, a multi-channel, dual frequency, P-code GPS time receiver, is one of the first TTR-4P units in operation in a civil laboratory. Preliminary comparative studies of this new equipment with conventional GPS time receivers are described in this paper. The results of an experimental restitution of GPS time, obtained in June 1993, are also detailed: 3 to 6 satellites were observed simultaneously with a sample interval of 15 s, an efficient smoothing of SA noise was realized by averaging data on all observed satellites over more than 1 hour. When the GPS system is complete in 1994, 8 satellites will be observable continuously from anywhere in the world and the same level of uncertainty will be obtained using a shorter averaging time.

Thomas, Claudine↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

[Compton GRO Cycle 7]

The goals of this project were to use the improved location accuracy of combined BATSE/Ulysses-IPN gamma-ray burst (GRB) data to search for small-angle clustering within the burst population and to search for possible correlations between bursts and catalogs of known astrophysical objects. The research was completed in several phases due to the length of time needed to process the basic BATSE/Ulysses data and to develop the appropriate analysis algorithms. In the first phase, BATSE/Ulysses data were gathered from their respective sources and processed to yield, final BATSE burst catalog locations and final BATSE/Ulysses IPN timing annuli. The final data are published. In the next phase, traditional statistical analysis techniques were used to search a subset of the data for small-angle clustering and. correlations, and to produce upper limits. Next, the unique analysis techniques and software described in the proposal were developed and,applied to the final set of data (corresponding to the 4B catalog). The final search did not find any significant self-clustering or correlations with supernovae, AGN, or near by galaxy clusters. The investigation ruled-out earlier findings that indicated (with marginal significance) such results, and. provided stringent limits on the allowable fraction of bursts that are self-clustered and associated with selected known objects. An added benefit of this research came in 1998 when it was suggested that some gamma-ray bursts may be associated with supernova (SN) explosions. We responded rapidly and applied the tools of this investigation to study this possibility. Using combined BATSE/Ulysses data, we were able to completely rule out any association between bright bursts and known supernovae, thereby questioning the validity of any GRB/SN association. If a GRB/SN association exists, it must be only with a small fraction of weak bursts.

Kippen, R. Marc↗

Modeling of polychromatic attenuation using computed tomography reconstructed images

This paper presents a procedure for estimating an accurate model of the CT imaging process including spectral effects. As raw projection data are typically unavailable to the end-user, we adopt a post-processing approach that utilizes the reconstructed images themselves. This approach includes errors from x-ray scatter and the nonidealities of the built-in soft tissue correction into the beam characteristics, which is crucial to beam hardening correction algorithms that are designed to be applied directly to CT reconstructed images. We formulate this approach as a quadratic programming problem and propose two different methods, dimension reduction and regularization, to overcome ill conditioning in the model. For the regularization method we use a statistical procedure, Cross Validation, to select the regularization parameter. We have constructed step-wedge phantoms to estimate the effective beam spectrum of a GE CT-I scanner. Using the derived spectrum, we computed the attenuation ratios for the wedge phantoms and found that the worst case modeling error is less than 3% of the corresponding attenuation ratio. We have also built two test (hybrid) phantoms to evaluate the effective spectrum. Based on these test phantoms, we have shown that the effective beam spectrum provides an accurate model for the CT imaging process. Last, we used a simple beam hardening correction experiment to demonstrate the effectiveness of the estimated beam profile for removing beam hardening artifacts. We hope that this estimation procedure will encourage more independent research on beam hardening corrections and will lead to the development of application-specific beam hardening correction algorithms.

NASA Discipline Musculoskeletal↗

Launch Vehicle Propulsion Parameter Design Multiple Selection Criteria

The optimization tool described herein addresses and emphasizes the use of computer tools to model a system and focuses on a concept development approach for a liquid hydrogen/liquid oxygen single-stage-to-orbit system, but more particularly the development of the optimized system using new techniques. This methodology uses new and innovative tools to run Monte Carlo simulations, genetic algorithm solvers, and statistical models in order to optimize a design concept. The concept launch vehicle and propulsion system were modeled and optimized to determine the best design for weight and cost by varying design and technology parameters. Uncertainty levels were applied using Monte Carlo Simulations and the model output was compared to the National Aeronautics and Space Administration Space Shuttle Main Engine. Several key conclusions are summarized here for the model results. First, the Gross Liftoff Weight and Dry Weight were 67% higher for the design case for minimization of Design, Development, Test and Evaluation cost when compared to the weights determined by the minimization of Gross Liftoff Weight case. In turn, the Design, Development, Test and Evaluation cost was 53% higher for optimized Gross Liftoff Weight case when compared to the cost determined by case for minimization of Design, Development, Test and Evaluation cost. Therefore, a 53% increase in Design, Development, Test and Evaluation cost results in a 67% reduction in Gross Liftoff Weight. Secondly, the tool outputs define the sensitivity of propulsion parameters, technology and cost factors and how these parameters differ when cost and weight are optimized separately. A key finding was that for a Space Shuttle Main Engine thrust level the oxidizer/fuel ratio of 6.6 resulted in the lowest Gross Liftoff Weight rather than at 5.2 for the maximum specific impulse, demonstrating the relationships between specific impulse, engine weight, tank volume and tank weight. Lastly, the optimum chamber pressure for Gross Liftoff Weight minimization was 2713 pounds per square inch as compared to 3162 for the Design, Development, Test and Evaluation cost optimization case. This chamber pressure range is close to 3000 pounds per square inch for the Space Shuttle Main Engine.

Shelton, Joey Dewayne↗

Thermal and Mechanical Property Characterization of the Advanced Disk Alloy LSHR

A low solvus, high refractory (LSHR) powder metallurgy disk alloy was recently designed using experimental screening and statistical modeling of composition and processing variables on sub-scale disks to have versatile processing-property capabilities for advanced disk applications. The objective of the present study was to produce a scaled-up disk and apply varied heat treat processes to enable full-scale demonstration of LSHR properties. Scaled-up disks were produced, heat treated, sectioned, and then machined into specimens for mechanical testing. Results indicate the LSHR alloy can be processed to produce fine and coarse grain microstructures with differing combinations of strength and time-dependent mechanical properties, for application at temperatures exceeding 1300 F.

Gabb, Timothy P.↗

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.↗

Improved Hierarchical Optimization-Based Classification of Hyperspectral Images Using Shape Analysis

A new spectral-spatial method for classification of hyperspectral images is proposed. The HSegClas method is based on the integration of probabilistic classification and shape analysis within the hierarchical step-wise optimization algorithm. First, probabilistic support vector machines classification is applied. Then, at each iteration two neighboring regions with the smallest Dissimilarity Criterion (DC) are merged, and classification probabilities are recomputed. The important contribution of this work consists in estimating a DC between regions as a function of statistical, classification and geometrical (area and rectangularity) features. Experimental results are presented on a 102-band ROSIS image of the Center of Pavia, Italy. The developed approach yields more accurate classification results when compared to previously proposed methods.

Classification↗

Remote Sensing, Uncertainty Quantification, and a Theory of Data Systems; Workshop Report

The purpose of the workshop was to invite statisticians, applied mathematicians, computer scientists, data system architects, experts in remote sensing technology, and Climate and Earth System scientists to review, discuss, and plan research on issues related to large-scale, efficient analysis of distributed data using spatial statistical methods. Our motivation in organizing this event was to catalyze interchange among experts on the fast-emerging problem of analysis of distributed data. As part of SAMSI's 2017-2018 Program on Mathematical and Statistical Methods for Climate and the Earth System, a Working Group on Remote Sensing was established to address statistical and mathematical research problems in the analysis of remote sensing data. The Working Group has five subgroups: 1) Spatial Retrieval Methodology (the so-called \Spatial-X" subgroup); 2) Spatial Analysis for Hyperspectral Data (the so-called \Spatial-Y" subgroup); 3) Emulators for Complex Forward Models; 4) Optimization for Remote Sensing Retrievals; and 5) Theory of Data Systems (ToDS). The ToDS subgroup spent the first half of this academic year formulating a framework in which to consider the joint problem of a) optimizing statistical methods for environments where data are distributed and too large to move to a central location, and b) the design of data system infrastructures within which to implement those statistical methods. To x ideas, the Workshop focused on spatial statistical methods. To date there are many new spatial statistical methods designed with massive data sets in mind, in the literature. However, very few have been implemented for remote sensing data, and none have been implemented in operational settings like those used by NASA and NOAA. A major impediment to their use in these cases is that the data are not only massive, but are stored in different physical locations. These data must be brought together in some way in order to estimate spatial covariance functions, but moving data to a central location for analysis is tedious at best and impossible at worst. Some remote data reduction is almost certainly necessary, but how much? What are the consequences for inference? The fundamental issue underlying these questions is how to navigate the trade-space between costs and uncertainty in the estimates or inferences that are ultimately produced.

Braverman, Amy↗

POLARRIS: A POLArimetric Radar Retrieval and Instrument Simulator

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Zh, Zdr, Kdp, ρhv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Zdr, Kdp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Toshi Matsui↗

Effects of Debulking on the Fiber Microstructure and Void Distribution in Carbon Fiber Reinforced Plastics

Carbon Fiber Reinforced Plastics (CFRPs) are widely used due to their high stiffness to weight ratios. A common process manufacturers use to increase the strength to weight ratio is debulking. Debulking is the process of compacting a dry fibrous reinforcement prior to resin infusion. This process is meant to decrease the average inter-fiber distance, effectively increasing the fiber volume fraction of the sample. While this process is widely understood macroscopically its effects on fibrous microstructures have not yet been well characterized. The aim of this work is to compare the microstructures of three CFRP laminates, varying only the debulking step in the manufacturing process. High resolution serial sections of all three laminates were taken for analysis. Using these scans, the fiber positions were reconstructed. Statistical descriptors such as local fiber and void volume fractions, fiber orientation, and void distribution and morphology were then generated for each sample. Fiber clusters present within the material were identified and analyzed for each level of debulking applied. Using these descriptors, the effects of debulking on the morphology and organization of the composite microstructure was evaluated.

carbon fiber↗