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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An automated procedure built on MTEX for reconstructing deformation twin hierarchies from electron backscattered diffraction datasets of heavily twinned microstructures

Here we present a set of algorithms built on the MTEX and MATLAB graph toolboxes for automatic reconstruction of deformation twin hierarchies from Electron Backscatter Diffraction (EBSD) datasets with a focus on developing methods for heavily twinned microstructures (twin fractions >0.5). The algorithms address key issues arising at large strains, mainly: missing twin relationships, grouping of heavily deformed grain fragments into families of similar orientation originating from a single initial grain, identification of parent fragments for large twin volume fractions, and classification of families having twin relationships with multiple families. To facilitate the development of these algorithms, large-grained ultra-high purity α-Ti deformed in compression along two directions is investigated. Graphs are utilized to handle non-local geometric merging and to represent relationships throughout the reconstruction process. When determining if a grain fragment is from the undeformed microstructure, the combined metrics of the fragment's orientation volume fraction in the initial texture and the directed graph centrality measure of out-closeness (the number of nodes reached in a graph from a given node) are essential. To address automation in reconstructing the sequence of twinning and relating fragments originating from a single grain in the initial microstructure, the twin family tree is formulated as a minimum spanning tree emanating from the initial grain family. A scheme constructing the distances associated with twin relationship comprising the spanning tree is developed, and a novel quasi-directional Prim spanning tree algorithm is used to determine the twin family tree. The procedure is demonstrated to significantly improve the level of automation in reconstructing twin hierarchies in heavily twinned microstructure compared to other methodologies in literature. The procedure can readily be applied to analyses of twinning in metals, as well as provide an approach for routinely extracting twin statistics at larger deformation levels than previously possible. Significantly, the procedure is demonstrated to be capable of identifying third generation twinning in α-Ti microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tracking atomic structure evolution during directed electron beam induced Si-atom motion in graphene via deep machine learning

Using electron beam manipulation, we enable deterministic motion of individual Si atoms in graphene along predefined trajectories. Structural evolution during the dopant motion was explored, providing information on changes of the Si atom neighborhood during atomic motion and providing statistical information of possible defect configurations. The combination of a Gaussian mixture model and principal component analysis applied to the deep learning-processed experimental data allowed disentangling of the atomic distortions for two different graphene sublattices. This approach demonstrates the potential of e-beam manipulation to create defect libraries of multiple realizations of the same defect and explore the potential of symmetry breaking physics. The rapid image analytics enabled via a deep learning network further empowers instrumentation for e-beam controlled atom-by-atom fabrication. Here, the analysis described in the paper can be reproduced via an interactive Jupyter notebook at https://git.io/JJ3Bx.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Field-Dependent Ionic Conductivities from Generalized Fluctuation-Dissipation Relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian, and subsequently, the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian, and the conductivity increases with applied field. Furthermore, this nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigation into the instantaneous centre of rotation for enhanced design of floating offshore wind turbines

The dynamic behaviour of floating offshore wind turbines (FOWTs) involves complex interactions of multivariate loads from wind, waves, and currents, which result in complex motion characteristics. Although methods for analysing global motion responses are well-established, the time- and location-dependent kinematics remain underexplored. This paper investigates the instantaneous centre of rotation (ICR), a point of zero velocity at a time instance of general plane motion. Understanding and strategically positioning the ICR can reduce the dynamic motion in critical structural locations, enhancing the performance and structural robustness of FOWTs. The paper presents a method for computing the ICR using time-domain simulation results and proposes a statistical analysis approach suitable for design studies. Building on prior research, it examines the sensitivity of the ICR to external loading and design features, providing insights into how these factors influence motion response and how the motion response influences the statistics of the ICR, structural loads, and other performance metrics of interest. The study explores two FOWT configurations, a spar and a semisubmersible, identifying design variables that most effectively control the ICR statistics and identifying the ICR statistics most correlated with the responses of interest. Finally, through two case studies, we demonstrate how to apply these new insights in a practical design scenario. By adjusting the design variables most correlated with the ICR (fairlead vertical position and centre of mass for the spar and mooring line length and offset column diameter for the semisubmersible), we successfully modified the designs of the floating support structures to reduce the loads in the mooring lines, tower base, and blade roots, improving the ultimate strength and fatigue characteristics compared to the original designs.

17 WIND ENERGY↗
NASA NTRS1977

A reduction in ag/residential signature conflict using principal components analysis of LANDSAT temporal data

Methods to accurately delineate the types of land cover in the urban-rural transition zone of metropolitan areas were considered. The application of principal components analysis to multidate LANDSAT imagery was investigated as a means of reducing the overlap between residential and agricultural spectral signatures. The statistical concepts of principal components analysis were discussed, as well as the results of this analysis when applied to multidate LANDSAT imagery of the Washington, D.C. metropolitan area.

Williams, D. L.↗
NASA NTRS1977

A reduction in ag./residential signature conflict using principal components analysis of Landsat temporal data

One important objective of a cooperative project between the U.S. Bureau of Census and NASA is to develop the ability to accurately delineate the types of land cover in the urban-rural transition zone of metropolitan areas. The application of principal components analysis to multidate Landsat imagery is being investigated as a method of reducing the overlap between residential and agricultural spectral signatures. The statistical concepts of principal components analysis are discussed, as well as the results of this analysis when applied to multidate Landsat imagery of the Washington, D.C. metropolitan area.

Williams, D. L.↗
NASA NTRS1983

FUNSTAT and statistical image representations

General ideas of functional statistical inference analysis of one sample and two samples, univariate and bivariate are outlined. ONESAM program is applied to analyze the univariate probability distributions of multi-spectral image data.

Parzen, E.↗
NASA NTRS1987

Satellite-derived attributes of cloud vortex systems and their application to climate studies

Defense Meteorological Satellite Program (DMSP) visible and infrared mosaics are analyzed in conjunction with synoptic meteorological observations of sea level pressure (SLP) and upper-air height to derive composite patterns of cyclonic cloud vortices for the Northern Hemisphere. The patterns reveal variations in the structure and implied dynamics of cyclonic systems at different stages of development that include: (1) increasing vertical symmetry of the lower-level and upper-air circulations and (2) decreasing lower-tropospheric thicknesses and temperature advection, associated with increasing age of the vortex. Cloud vortices are more intense in winter than in summer and typically reach maximum intensity in the short-lived prespiral signature stage. There are major structural differences among frontal wave, polar air, and 'instant occlusion' cyclogenesis types. Cyclones in the dissipation stage may reintensify (deepen), as denoted by the appearance in the imagery of an asymmetric cloud band or a tightened spiral vortex. The satellite-derived statistics on cloud vortex intensity, which are seasonal- and latitude- as well as type-dependent, are applied to a preliminary examination of the synoptic manifestations of seasonal climate variability. An apparently close relationship is found, for two winter and spring seasons, between Northern Hemisphere cyclonic activity and variations in cryosphere variables, particularly the extent of Arctic sea ice. The results may indicate that increased snow and ice extent accompany a southward displacement of cyclonic activity and/or a predominance of deeper systems. However, there is also a strong regional dependence to the ice-synoptics feedback. This study demonstrates the utility of high resolution meteorological satellite imagery for studies of climate variations (climate dynamics).

Carleton, Andrew M.↗
NASA NTRS1988

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗
NASA NTRS1993

Study of multi-dimensional radiative energy transfer in molecular gases

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical arrow band model with an exponential-tailed inverse intensity distribution. Consideration of spectral correlation results in some distinguishing features of the Monte Carlo formulations. Validation of the Monte Carlo formulations has been conducted by comparing results of this method with other solutions. Extension of a one-dimensional problem to a multi-dimensional problem requires some special treatments in the Monte Carlo analysis. Use of different assumptions results in different sets of Monte Carlo formulations. The nongray narrow band formulations provide the most accurate results.

Liu, Jiwen↗
NASA NTRS1992

Correlation, Cost Risk, and Geometry

The geometric viewpoint identifies the choice of a correlation matrix for the simulation of cost risk with the pairwise choice of data vectors corresponding to the parameters used to obtain cost risk. The correlation coefficient is the cosine of the angle between the data vectors after translation to an origin at the mean and normalization for magnitude. Thus correlation is equivalent to expressing the data in terms of a non orthogonal basis. To understand the many resulting phenomena requires the use of the tensor concept of raising the index to transform the measured and observed covariant components into contravariant components before vector addition can be applied. The geometric viewpoint also demonstrates that correlation and covariance are geometric properties, as opposed to purely statistical properties, of the variates. Thus, variates from different distributions may be correlated, as desired, after selection from independent distributions. By determining the principal components of the correlation matrix, variates with the desired mean, magnitude, and correlation can be generated through linear transforms which include the eigenvalues and the eigenvectors of the correlation matrix. The conversion of the data to a non orthogonal basis uses a compound linear transformation which distorts or stretches the data space. Hence, the correlated data does not have the same properties as the uncorrelated data used to generate it. This phenomena is responsible for seemingly strange observations such as the fact that the marginal distributions of the correlated data can be quite different from the distributions used to generate the data. The joint effect of statistical distributions and correlation remains a fertile area for further research. In terms of application to cost estimating, the geometric approach demonstrates that the estimator must have data and must understand that data in order to properly choose the correlation matrix appropriate for a given estimate. There is a general feeling by employers and managers that the field of cost requires little technical or mathematical background. Contrary to that opinion, this paper demonstrates that a background in mathematics equivalent to that needed for typical engineering and scientific disciplines at the masters or doctorate level is appropriate within the field of cost risk.

Dean, Edwin B.↗
NASA NTRS2010

A Comparison of Tropical Storm (TS) and Non-TS Gust Factors for Assessing Peak Wind Probabilities at the Eastern Range

Peak wind speed is an important forecast element to ensure the safety of personnel and flight hardware at Kennedy Space Center (KSC) and the Cape Canaveral Air Force Station (CCAFS) in East-Central Florida. The 45th Weather Squadron (45 WS), the organization that issues forecasts for the KSC/CCAFS area, finds that peak winds are more difficult to forecast than mean winds. This difficulty motivated the 45 WS to request two independent studies. The first (Merceret 2009) was the development of a reliable model for gust factors (GF) relating the peak to the mean wind speed in tropical storms (TS). The second (Lambert et al. 2008) was a climatological study of non-TS cool season (October-April) mean and peak wind speeds by the Applied Meteorology Unit (AMU; Bauman et al. 2004) without the use of GF. Both studies presented their statistics as functions of mean wind speed and height. Most of the few comparisons of TS and non-TS GF in the literature suggest that non-TS GF at a given height and mean wind speed are smaller than the corresponding TS GF. The investigation reported here converted the non-TS peak wind statistics calculated by the AMU to the equivalent GF statistics and compared them with the previous TS GF results. The advantage of this effort over all previously reported studies of its kind is that the TS and non-TS data were taken from the same towers in the same locations. This eliminates differing surface attributes, including roughness length and thermal properties, as a major source of variance in the comparison. The goal of this study is two-fold: to determine the relationship between the non-TS and TS GF and their standard deviations (GFSD) and to determine if models similar to those developed for TS data in Merceret (2009) could be developed for the non-TS environment. The results are consistent with the literature, but include much more detailed, quantitative information on the nature of the relationship between TS and non-TS GF and GFSD as a function of height and mean wind speed.

Merceret, Francis J.↗
NASA NTRS2011

Quantitative Examination of a Large Sample of Supra-Arcade Downflows in Eruptive Solar Flares

Sunward-flowing voids above post-coronal mass ejection flare arcades were first discovered using the soft X-ray telescope aboard Yohkoh and have since been observed with TRACE (extreme ultraviolet (EUV)), SOHO/LASCO (white light), SOHO/SUMER (EUV spectra), and Hinode/XRT (soft X-rays). Supra-arcade downflow (SAD) observations suggest that they are the cross-sections of thin flux tubes retracting from a reconnection site high in the corona. Supra-arcade downflowing loops (SADLs) have also been observed under similar circumstances and are theorized to be SADs viewed from a perpendicular angle. Although previous studies have focused on dark flows because they are easier to detect and complementary spectral data analysis reveals their magnetic nature, the signal intensity of the flows actually ranges from dark to bright. This implies that newly reconnected coronal loops can contain a range of hot plasma density. Previous studies have presented detailed SAD observations for a small number of flares. In this paper, we present a substantial SADs and SADLs flare catalog. We have applied semiautomatic detection software to several of these events to detect and track individual downflows thereby providing statistically significant samples of parameters such as velocity, acceleration, area, magnetic flux, shrinkage energy, and reconnection rate. We discuss these measurements (particularly the unexpected result of the speeds being an order of magnitude slower than the assumed Alfven speed), how they were obtained, and potential impact on reconnection models.

Savage, Sabrina L.↗
NASA NTRS2012

Representation of Probability Density Functions from Orbit Determination using the Particle Filter

Statistical orbit determination enables us to obtain estimates of the state and the statistical information of its region of uncertainty. In order to obtain an accurate representation of the probability density function (PDF) that incorporates higher order statistical information, we propose the use of nonlinear estimation methods such as the Particle Filter. The Particle Filter (PF) is capable of providing a PDF representation of the state estimates whose accuracy is dependent on the number of particles or samples used. For this method to be applicable to real case scenarios, we need a way of accurately representing the PDF in a compressed manner with little information loss. Hence we propose using the Independent Component Analysis (ICA) as a non-Gaussian dimensional reduction method that is capable of maintaining higher order statistical information obtained using the PF. Methods such as the Principal Component Analysis (PCA) are based on utilizing up to second order statistics, hence will not suffice in maintaining maximum information content. Both the PCA and the ICA are applied to two scenarios that involve a highly eccentric orbit with a lower apriori uncertainty covariance and a less eccentric orbit with a higher a priori uncertainty covariance, to illustrate the capability of the ICA in relation to the PCA.

Mashiku, Alinda K.↗
NASA NTRS2015

Characterization of Machine Variability and Progressive Heat Treatment in Selective Laser Melting of Inconel 718

The absence of an economy of scale in spaceflight hardware makes additive manufacturing an immensely attractive option for propulsion components. As additive manufacturing techniques are increasingly adopted by government and industry to produce propulsion hardware in human-rated systems, significant development efforts are needed to establish these methods as reliable alternatives to conventional subtractive manufacturing. One of the critical challenges facing powder bed fusion techniques in this application is variability between machines used to perform builds. Even with implementation of robust process controls, it is possible for two machines operating at identical parameters with equivalent base materials to produce specimens with slightly different material properties. The machine variability study presented here evaluates 60 specimens of identical geometry built using the same parameters. 30 samples were produced on machine 1 (M1) and the other 30 samples were built on machine 2 (M2). Each of the 30-sample sets were further subdivided into three subsets (with 10 specimens in each subset) to assess the effect of progressive heat treatment on machine variability. The three categories for post-processing were: stress relief, stress relief followed by hot isostatic press (HIP), and stress relief followed by HIP followed by heat treatment per AMS 5664. Each specimen (a round, smooth tensile) was mechanically tested per ASTM E8. Two formal statistical techniques, hypothesis testing for equivalency of means and one-way analysis of variance (ANOVA), were applied to characterize the impact of machine variability and heat treatment on six material properties: tensile stress, yield stress, modulus of elasticity, fracture elongation, and reduction of area. This work represents the type of development effort that is critical as NASA, academia, and the industrial base work collaboratively to establish a path to certification for additively manufactured parts. For future flight programs, NASA and its commercial partners will procure parts from vendors who will use a diverse range of machines to produce parts and, as such, it is essential that the AM community develop a sound understanding of the degree to which machine variability impacts material properties.

Prater, T.↗