Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:data analysis”

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 1,099 records · Page 61

Classifying and analyzing small-angle scattering data using weighted k nearest neighbors machine learning techniques

A consistent challenge for both new and expert practitioners of small-angle scattering (SAS) lies in determining how to analyze the data, given the limited information content of said data and the large number of models that can be employed. Machine learning (ML) methods are powerful tools for classifying data that have found diverse applications in many fields of science. Here, ML methods are applied to the problem of classifying SAS data for the most appropriate model to use for data analysis. The approach employed is built around the method of weighted k nearest neighbors (wKNN), and utilizes a subset of the models implemented in the SasView package (https://www.sasview.org/) for generating a well defined set of training and testing data. The prediction rate of the wKNN method implemented here using a subset of SasView models is reasonably good for many of the models, but has difficulty with others, notably those based on spherical structures. A novel expansion of the wKNN method was also developed, which uses Gaussian processes to produce local surrogate models for the classification, and this significantly improves the classification accuracy. Further, by integrating a stochastic gradient descent method during post-processing, it is possible to leverage the local surrogate model both to classify the SAS data with high accuracy and to predict the structural parameters that best describe the data. The linking of data classification and model fitting has the potential to facilitate the translation of measured data into results for both novice and expert practitioners of SAS.

97 MATHEMATICS AND COMPUTING↗

Influence of Averaging Method on the Evaluation of a Coastal Ocean Color Event on the U.S. Northeast Coast

Application of appropriate spatial averaging techniques is crucial to correct evaluation of ocean color radiometric data, due to the common log-normal or mixed log-normal distribution of these data. Averaging method is particularly crucial for data acquired in coastal regions. The effect of averaging method was markedly demonstrated for a precipitation-driven event on the U.S. Northeast coast in October-November 2005, which resulted in export of high concentrations of riverine colored dissolved organic matter (CDOM) to New York and New Jersey coastal waters over a period of several days. Use of the arithmetic mean averaging method created an inaccurate representation of the magnitude of this event in SeaWiFS global mapped chl a data, causing it to be visualized as a very large chl a anomaly. The apparent chl a anomaly was enhanced by the known incomplete discrimination of CDOM and phytoplankton chlorophyll in SeaWiFS data; other data sources enable an improved characterization. Analysis using the geometric mean averaging method did not indicate this event to be statistically anomalous. Our results predicate the necessity of providing the geometric mean averaging method for ocean color radiometric data in the Goddard Earth Sciences DISC Interactive Online Visualization ANd aNalysis Infrastructure (Giovanni).

Acker, James G.↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Evidence against naive truncations of the OPE from e + e − → hadrons below charm

The operator product expansion (OPE), truncated in dimension, is employed in many contexts. An example is the extraction of the strong coupling, $\alpha_s$, from hadronic $\tau$-decay data, using a variety of analysis methods based on finite-energy sum rules. Here, we reconsider a long-used method, which parametrizes non-perturbative contributions to the $I=1$ vector and axial vacuum polarizations with the OPE, setting several higher-dimension coefficients to zero in order to implement the method in practice. The assumption that doing this has a negligible effect on the value of $\alpha_s$ is tantamount to the assumption that the low-dimension part of the OPE converges rapidly with increasing dimension near the $\tau$ mass. Were this assumption valid, it would certainly have to be valid at energies above the $\tau$ mass as well. It follows that the method can be tested using data obtained from $e^+e^-\to\mbox{hadrons}$, as they are not limited by the kinematic constraints of $\tau$ decays. We carry out such an investigation using a recent high-precision compilation for the $R$ ratio, arguing that it provides insights into the validity of the strategy, even if it probes a different, though related channel. We find that $e^+e^-$-based tests call into question the implied assumption of rapid convergence of the low-dimension part of the OPE around the $\tau$ mass, and thus underscore the need to restrict finite-energy sum-rule analyses to observables which receive only contributions from lower-order terms in the OPE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation, Analysis, and Application of Internal Strain-Gage Balance Data

Experimental processes, analytical methods, and numerical algorithms are described that may be used to predict the forces and moments of an internal strain-gage balance during a wind tunnel test. First, the control volume model of a strain-gage balance and the concepts of load state, load space, and output space are introduced. These important abstractions provide a better understanding of fundamental characteristics of different balance load prediction approaches. Then, the description of strain-gage balance data and the definition of the primary gage sensitivity is discussed. Afterwards, basic elements of the calibration of a typical six-component balance are reviewed. Two fundamentally different balance load prediction methods, the processing of check loads, and related topics are also discussed. Three real-world balance data examples are reviewed in great detail to illustrate typical analysis results for a variety of strain-gage balance designs. Finally, important observations are summarized and recommendations are provided. Additional information and detailed mathematical derivations can be found in the appendices of the document. They include the following topics: balance terminology, definitions of important statistical metrics, balance axis system conventions, balance load transformations, the combined load diagram, electrical output format options, bi-directional gage output characteristics, determination of the natural zeros, derivation of two balance load prediction methods, description of two tare load iteration algorithms, modeling of balance temperature effects, basics of three-component moment balances, definition of the percent contribution, detection of linear and near-linear dependencies in balance calibration data, a regression model term selection algorithm, and other related topics.

wind tunnel test↗

Two-dimensional wave-number spectral analysis techniques for phase contrast imaging turbulence imaging data on large helical device

An analysis method for unfolding the spatially resolved wave-number spectrum and phase velocity from the 2D CO{sub 2} laser phase contrast imaging system on the large helical device is described. This is based on the magnetic shear technique which identifies propagation direction from 2D spatial Fourier analysis of images detected by a 6 × 8 detector array. Because the strongest modes have wave-number at the lower end of the instrumental k range, high resolution spectral techniques are necessary to clearly resolve the propagation direction and hence the spatial distribution of fluctuations along the probing laser beam. Multiple-spatial point cross-correlation averaging is applied before calculating the spatial power spectrum. Different methods are compared, and it is found that the maximum entropy method (MEM) gives best results. The possible generation of artifacts from the over-narrowing of spectra are investigated and found not to be a significant problem. The spatial resolution Δρ (normalized radius) around the peak wave-number, for conventional Fourier analysis, is ∼0.5, making physical interpretation difficult, while for MEM, Δρ ∼ 0.1.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of the sound speed in dense fluid deuterium along the cryogenic liquid Hugoniot

Many experimental studies, spanning several decades of research and using various dynamic compression schemes, have been conducted to investigate cryogenic liquid deuterium under strong shock compression. The consensus emerging from these studies is that of a progressive dissociation of the D 2 molecules into an electrically conducting, atomic plasma, when subjected to shock pressures exceeding ~50 GPa. While state-of-the-art numerical simulations based on density-functional-theory or quantum Monte-Carlo techniques capture this behavior quite well, subtle differences subsist between these simulations and the available experimental data regarding the pressure-density compressibility. Here, leveraging a recently developed analysis method for high-resolution Doppler interferometric velocity data, we present Eulerian sound speed measurements in compressed deuterium to shock pressures between 50 and 200 GPa. These results, extracted from laser-driven shockwave experiments, are found to agree with several of the most accurate equation of state models for deuterium at those conditions up to ~150 GPa. However, the data indicate that these models fail to reproduce the experimentally observed sound speed at higher pressures, approaching 200 GPa. In particular, we unveil a discrepancy between the experimental results and the equation of state model that is most commonly used in inertial confinement fusion at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of three-component aeromagnetic data

The numerical method of obtaining three field components from total field measurements, using double Fourier series expansion, is presented. The expressions for moments of the anomalous field components over a finite area are given. The magnitude and direction of the magnetization vector indicate that the vertical component of the magnetic field calculated from total field observations is more accurate at higher geomagnetic latitudes than at lower latitudes. The opposite is true for the horizontal components. The error in determining the magnetization vector directions are significantly large over most of the range of variation of declination and inclination of the vector, demonstrating the practical limitations of computing field components from total field data even under the best of conditions.

Bhattacharyya, B. K.↗

Recent advances in thermal-structural analysis and design

A continuing effort to develop a comprehensive capability for thermal-structural analysis and automated design (sizing) is reviewed. A principal role in the activity is played by the finite element program SPAR which contains both an efficient structural and thermal analysis capability. The benefit of having thermal and structural analyses in the same finite element program is illustrated by the application of SPAR to design calculations for the National Transonic Facility - a cryogenic wind tunnel under construction at Langley. Some experience with large-scale thermal structural analysis problems, particularly the space shuttle orbiter, has led to the identification of some analysis needs. Those needs include automated model generation and data output for lumped parameter thermal analysis, faster solution methods for nonlinear transient heat transfer, automated interpolation of temperature data from a thermal finite element model to a dissimilar structural finite model, and automated techniques to identify the times at which the critical combinations of transient heating and loads occur on a structure.

Adelman, H. M.↗

Scientific study of data analysis

We present a comparison between two numerical methods for the extrapolation of nonlinear force-free magnetic fields, the Iterative Method (IM) and the Progressive Extension Method (PEM). The advantages and disadvantages of these two methods are summarized and the accuracy and numerical instability are discussed. On the basis of this investigation, we claim that the two methods do resemble each other qualitatively.

Wu, S. T.↗

Cleanliness evaluation of rough surfaces with diffuse IR reflectance

Contamination on bonding surfaces has been determined to be a primary cause for degraded bond strength in certain solid rocket motor bondlines. Hydrocarbon and silicone based organic contaminants that are airborne or directly introduced to a surface are a significant source of contamination. Diffuse infrared (IR) reflectance has historically been used as an effective technique for detection of organic contaminants, however, common laboratory methods involving the use of a Fourier transform IR spectrometer (FTIR) are impractical for inspecting the large bonding surface areas found on solid rocket motors. Optical methods involving the use of acousto-optic tunable filters and fixed bandpass optical filters are recommended for increased data acquisition speed. Testing and signal analysis methods are presented which provide for simultaneous measurement of contamination concentration and roughness level on rough metal surfaces contaminated with hydrocarbons.

Pearson, L. H.↗

Method for Visually Integrating Multiple Data Acquisition Technologies for Real Time and Retrospective Analysis

A system for display on a single video display terminal of multiple physiological measurements is provided. A subject is monitored by a plurality of instruments which feed data to a computer programmed to receive data, calculate data products such as index of engagement and heart rate, and display the data in a graphical format simultaneously on a single video display terminal. In addition live video representing the view of the subject and the experimental setup may also be integrated into the single data display. The display may be recorded on a standard video tape recorder for retrospective analysis.

Bogart, Edward H.↗

Analysis Methods and Results for Weak Gamma-Ray Bursts in the BATSE Data

We report initial results on the statistical properties of the dimmest gamma-ray bursts (GRBs) observed with the Burst and Transient Source Experiment (BATSE), using new ground-based methods to obtain a sample of GRBs from 502 days of BATSE data. Using the most sensitive ground-based detection of GRBs, the sample extends to GRBs much fainter than those detected by the on-board trigger, but because of the temporal resolution of the data, the sample is limited to GRBs of duration of at least 2(approx.)s. For each detected event, Bayesian probabilities are calculated for the event to belong to each of seven classes of differing physical origins. The sample of GRB candidates is defined by the requirement that the Bayesian probability for belonging to the GRB class is higher than 0.5. The intensity distribution of the GRB sample is corrected using a Monte Carlo simulation of the post-flight detection efficiency. The dimmest BATSE bursts of the sample continue the hardness-intensity trend seen in brighter GREs and are consistent with isotropy.

Mitrofanov, I. G.↗

Stochastic Richardson Extrapolation Based Numerical Error Estimation for Kinetic Plasma Simulations

We present a numerical error estimation technique specifically formulated to deal with stochastic code output with multiple discretization parameters. This method is based on multiple fits to an error model with arbitrary convergence rates and cross-coupling terms, performed using nonlinear optimization. The fitting approach varies by the type of residual norm which influences the importance of outliers, and weights which influences the relative importance of data points in the coarse and refined regions of discretization parameter space. To account for the influence of stochastic noise, these fits are performed on multiple bootstrap values based on the underlying response data set. Using an automated discretization domain selection scheme, the fits are performed on a series of reduced sets of discretization levels in order to find an optimal fully-converged result estimate in the minimum variance sense; this automated approach enables straightforward analysis of multiple quantities of interest and/or time and spatially-dependent response data. The overall numerical error analysis method is useful for verification and validation problems for stochastic simulation methods and forms a key component in the overall uncertainty quantification process. The method was demonstrated for steady and unsteady electron diode problems simulated using a particle-in-cell kinetic plasma code, demonstrating excellent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

How reliable is distribution of relaxation times (DRT) analysis? A dual regression-classification perspective on DRT estimation, interpretation, and accuracy

The distribution of relaxation times (DRT) has gained increasing attention and adoption in recent years as a versatile method for analyzing electrochemical impedance spectroscopy (EIS) data obtained from complex devices like fuel cells, electrolyzers, and batteries. The DRT deconvolutes the impedance without a priori specification of a generative model, which is especially useful for interpretation and model selection when the governing principles of the system under study are not fully understood. However, DRT estimation is an ill-posed inversion problem that must be addressed with a subjective choice of regularization and tuning, which leaves substantial risk of misleading interpretations of EIS data. In this work, we suggest a new classification view of the DRT inversion to clarify DRT estimation and interpretation. We introduce a dual regression-classification framework that unifies the classification and regression views of the DRT inversion with wide-reaching implications for DRT analysis. The dual framework is employed to demonstrate a new kind of DRT inversion algorithm and develop novel evaluation metrics that capture previously ignored aspects of DRT accuracy. These approaches are applied to both synthetic data and experimental spectra collected from a protonic ceramic fuel cell and a lithium-ion battery to illustrate their broad utility. The dual inversion algorithm shows promising performance for accurate DRT estimation and autonomous model identification, while the dual evaluation approach produces metrics that meaningfully assess the strengths and risks of DRT algorithms. Here this work provides valuable insight for both practical application of the DRT to experimental data and further development of EIS analysis methods.

36 MATERIALS SCIENCE↗

Cultivating an Emergent Earth Observation Analytics Ecosystem in the Cloud

A diverse set of data analytics systems for Earth Observations are sprouting up in the Earth Science community, with a wealth of processing algorithms and analysis methods. There is a similar wealth of data resources available via myriad data providers and clearinghouses, including large institutional systems like the Earth Observing System Data and Information System, Comprehensive Large Scale Array-data Stewardship System, and Federated Earth Observation Missions gateway. With Earth system science driving a need to work with more datasets together, and the community developing more analysis tools (some of them dataset-specific), how can we develop analysis workflows that incorporate far-flung datasets and leverage analysis resources from multiple organizations? Cloud computing points the way toward a solution in two different respects. Firstly, the access to and abstraction of virtually unlimited storage and computing power provides an environment that enables more straightforward means of pulling datasets and analysis resources together. Just as importantly, however, cloud computing serves as an example of an "ecosystem" of interoperating services, since the essence of cloud computing is the presentation of all resources as a service, from hardware to infrastructure to platform to software. This enables the combination of off-the-shelf, diverse services to construct entire systems that emerge out of an equally diverse community of architects and developers. This approach can be similarly applied to the data and analysis resources in the Earth Observation community. By exposing these resources via well understood services, and consuming resources in the same way, different organizations can construct bespoke analysis workflows and systems for their own purposes. The key leap the community needs to make is to develop analysis systems in components that interact with other components via services. The result would be a rich ecosystem of analytics components that can be combined to analyze datasets at scale and in conjunction with other datasets from other sources.

chaos↗

Fluorescence Correlation Spectroscopy with Photobleaching Correction in Slowly Diffusing Systems

Fluorescence correlation spectroscopy (FCS) is a powerful tool to quantitatively study the diffusion of fluorescently labeled molecules. It allows in principle important questions of macromolecular transport and supramolecular aggregation in living cells to be addressed. However, the crowded environment inside the cells slows diffusion and limits the reservoir of labeled molecules, causing artifacts that arise especially from photobleaching and limit the utility of FCS in these applications. We present a method to compute the time correlation function from weighted photon arrival times, which compensates computationally during the data analysis for the effect of photobleaching. We demonstrate the performance of this method using numerical simulations and experimental data from model solutions. Using this technique, we obtain correlation functions in which the effect of photobleaching has been removed and in turn recover quantitatively accurate mean-square displacements of the fluorophores, especially when deviations from an ideal Gaussian excitation volume are accounted for by using a reference calibration correlation function. This allows quantitative FCS studies of transport processes in challenging environments with substantial photobleaching like in living cells in the future.

74 ATOMIC AND MOLECULAR PHYSICS↗