SEARCH · Engineering Papers
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.
A brief overview of statistical tools needed to perform post flight/test reconstruction of state variables is given. Linear regression, recursive linear regression, and the exact connection between the Kalman filter and linear regression are discussed. The regression connection is expected to serve as an aid in the application of a recently developed analytical method of flight reconstruction to single engine test firing data.
Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble averaged results for a widely used 1-D turbulence approximation.
Delicate microgravity science is unlikely to succeed on the International Space Station if vibratory and transient disturbers corrupt the environment. An analytical approach to compute the on-orbit acceleration environment at science experiment locations within a standard payload rack resulting from these disturbers is presented. This approach has been grounded by correlation and comparison to test verified transfer functions. The method combines the results of finite element and statistical energy analysis using tested damping and modal characteristics to provide a reasonable approximation of the total root-mean-square (RMS) acceleration spectra at the interface to microgravity science experiment hardware.
Large and growing data resources on the spatial and temporal diversity and distribution of the more than 400 carbon-bearing mineral species reveal patterns of mineral evolution and ecology. Recent advances in analytical and visualization techniques leverage these data and are propelling mineralogy from a largely descriptive field into one of prediction within complex, integrated, multidimensional systems. These discoveries include: (1) systematic changes in the character of carbon minerals and their networks of coexisting species through deep time; (2) improved statistical predictions of the number and types of carbon minerals that occur on Earth but are yet to be discovered and described; and (3) a range of proposed and ongoing studies related to the quantification of network structures and trends, relation of mineral “natural kinds” to their genetic environments, prediction of the location of mineral species across the globe, examination of the tectonic drivers of mineralization through deep time, quantification of preservational and sampling bias in the mineralogical record, and characterization of feedback relationships between minerals and geochemical environments with microbial populations. These aspects of Earth’s carbon mineralogy underscore the complex coevolution of the geosphere and biosphere and highlight the possibility for scientific discovery in Earth and planetary systems.
The solar investment community is in need of up-to-date and data-derived system degradation rates to use in financial models that calculate the risk and expectation of energy revenue. To date, financial models choose degradation rates from past studies that are limited and may not be representative of solar projects in development today. As the solar industry expands in all sectors (residential, commercial, and utility), there is a growing amount of energy generation data available to analyze. With the help of the open source degradation analysis code (RdTools) it is now possible to derive degradation on an ongoing basis and continuously provide up-to-date durability statistics to the investment community. Ongoing and accurate statistics of fielded photovoltaic systems allows financial stakeholders to constrain energy revenue projections, lower financial risk, and increase the bankability of solar. This presentation will provide a mid-term progress update on a 2 year SETO-funded degradation project that involves applying RdTools analysis on kWh Analytics’ industry database. With results from the RdTools analysis in hand, machine learning methods are employed to quantify how much variance of system performance degradation can be explained by predictor variables such as environmental factors, equipment bill of materials characteristics, and system design information. The goal of this presentation is to promote the usage of RdTools on an ongoing to data owners and to solicit feedback from researchers and stakeholders on how to maximize impact of the SETO-funded kWh Analytics degradation project.
Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.
We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. Here, in this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case.
One of the challenging tasks in modern aquatic remote sensing is the retrieval of near-surface concentrations of Total Suspended Solids (TSS). This study aims to present a Statistical, inherent Optical property (IOP) -based, and muLti-conditional Inversion proceDure (SOLID) for enhanced retrievals of satellite-derived TSS under a wide range of in-water bio-optical conditions in rivers, lakes, estuaries, and coastal waters. In this study, using a large in situ database (N > 3500), the SOLID model is devised using a three-step procedure: (a) water-type classification of the input remote sensing reflectance (R(sub rs)), (b) retrieval of particulate backscattering (b(sub bp)) in the red or near-infrared (NIR) regions using semi-analytical, machine-learning, and empirical models, and (c) estimation of TSS from b(sub bp) via water-type-specific empirical models. Using an independent subset of our in situ data (N = 2729) with TSS ranging from 0.1 to 2626.8 [g/m (exp 3)], the SOLID model is thoroughly examined and compared against several state-of-the-art algorithms (Miller and McKee, 2004; Nechad et al., 2010; Novoa et al., 2017; Ondrusek et al., 2012; Petus et al., 2010). We show that SOLID outperforms all the other models to varying degrees, i.e., from 10 to > 100%, depending on the statistical attributes (e.g., global versus water-type-specific metrics). For demonstration purposes, the model is implemented for images acquired by the MultiSpectral Imager aboard Sentinel-2A/B over the Chesapeake Bay, San-Francisco-Bay-Delta Estuary, Lake Okeechobee, and Lake Taihu. To enable generating consistent, multimission TSS products, its performance is further extended to, and evaluated for, other missions, such as the Ocean and Land Color Instrument (OLCI), Moderate Resolution Imaging Spectroradiometer (MODIS), Visible Infrared Imaging Radiometer Suite (VIIRS), and Operational Land Imager (OLI). Sensitivity analyses on uncertainties induced by the atmospheric correction indicate that 10% uncertainty in Rrs leads to < 20% uncertainty in TSS retrievals from SOLID. While this study suggests that SOLID has a potential for producing TSS products in global coastal and inland waters, our statistical analysis certainly verifies that there is still a need for improving retrievals across a wide spectrum of particle loads.
The dominant paradigm for assessing ecological responses to climate change assumes that future states of individuals and populations can be predicted by current, species-wide performance variation across spatial climatic gradients. However, if the fates of ecological systems are better predicted by past responses to in situ climatic variation through time, this current analytical paradigm may be severely misleading. Empirically testing whether spatial or temporal climate responses better predict how species respond to climate change has been elusive, largely due to restrictive data requirements. Here, we leverage a newly collected network of ponderosa pine tree-ring time series to test whether statistically inferred responses to spatial versus temporal climatic variation better predict how trees have responded to recent climate change. When compared to observed tree growth responses to climate change since 1980, predictions derived from spatial climatic variation were wrong in both magnitude and direction. This was not the case for predictions derived from climatic variation through time, which were able to replicate observed responses well. Future climate scenarios through the end of the 21st century exacerbated these disparities. These results suggest that the currently dominant paradigm of forecasting the ecological impacts of climate change based on spatial climatic variation may be severely misleading over decadal to centennial timescales.
The purpose of this project was to utilize computational tools to understand the relationships between processing, microstructure, and properties for additively manufactured (AM) aluminum alloys for automotive applications, and to provide an engineering solution for helping to optimize process conditions. The project leverages ORNL developments in computational modeling, including AM process modeling, phase-field based microstructure evolution predictions, and data analytics techniques for mapping process conditions to material outcomes. The project utilized an Al-Cu-Mn-Zr alloy as a model material for studying formation of defects and microstructural features in response to variations in process conditions. Based on both pre-existing experimental data and simulation results, statistical process maps were constructed to identify regions of process space with minimal defect formation and advantageous microstructures and properties. The software tools used for this purpose were successful disseminated to GM, who were able to successful compile the relevant HPC codes within their own computing ecosystem and perform initial calculations to reproduce ORNL results.
We present three-dimensional Thomas-Fermi-Dirac calculations of lattice mixing energies of hydrogen with carbon and oxygen atoms, respectively. The results are used to derive effective interatomic potentials for use in liquid-state mixture calculations. We then use the potentials to derive analytic expressions for binary mixture-free energies and to map out the phase diagrams of mixtures of hydrogen with, respectively, helium, carbon, and oxygen, over a pressure range of about 5 to about 10 to the 3rd Mbar. Within this pressure range, all three of the latter elements are found to have unlimited solubility in metallic hydrogen over a temperature range which lies above their pure-element melting temperatures, and which includes likely interior temperatures in the Jovian planets.
Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second-order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single-path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble-averaged results for a widely used one-dimensional turbulence approximation.
Photon statistics and phase properties of two-mode squeezed number states are studied. It is shown that photon number distribution and Pegg-Barnett phase distribution for such states have similar (N + 1)-peak structure for nonzero value of the difference in the number of photons between modes. Exact analytical formulas for phase distributions based on different phase approaches are derived. The Pegg-Barnett phase distribution and the phase quasiprobability distribution associated with the Wigner function are close to each other, while the phase quasiprobability distribution associated with the Q function carries less phase information.
We investigate the general case of the photon distribution of a two-mode squeezed vacuum and show that the distribution of photons among the two modes depends on four parameters: two squeezing parameters, the relative phase between the two oscillators and their spatial orientation. The distribution of the total number of photons depends only on the two squeezing parameters. We derive analytical expressions and present pictures for both distributions.
A matrix of seven low-activity waste (LAW) glasses was fabricated and characterized to be used by the Hanford Waste Treatment and Immobilization Plant (WTP) as analytical instrument calibration and matrix interference standards when measuring LAW glass compositions. The matrix of glasses was designed to represent the range of LAW glass compositions that will be generated by the WTP. Samples from each of the seven glasses were measured with electron probe microanalysis (EPMA), ion chromatography (IC), and bulk inductively coupled plasma mass spectroscopy (ICP-MS) to determine the variability and accuracy across batches and within individually poured glass bar samples. High relative percent differences and standard deviations were calculated for components with low concentrations due to instrument detection limits. IC and ICP-MS percent differences and standard deviations were higher at low component concentrations. Otherwise, variability in measured compositions from IC and ICP-MS was generally lower than measurements from EPMA, with the exclusion of Cl, F, Si, and B. Overall, the two methods proved comparable and complementary. From EPMA, IC, and ICP-MS data, volatile elements were underrepresented compared to batched compositions, which suggested loss of these elements during melting. EPMA data was statistically analyzed to evaluate homogeneity within single bars and whole compositions for each of the seven glasses. Most of the variability within glass compositions occurred at the bar-to-bar level; however, no clear systematic trends were observed. Based on these analyses, when using these glasses as reference material, it is recommended that the overall mean and variance of each composition be used for performing instrument calibrations and analytical corrections.
We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.
Nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) spectroscopy of bulk quantum materials have provided insight into phenomena, such as quantum phase criticality, magnetism, and superconductivity. With the emergence of nanoscale 2D materials with magnetic phenomena, inductively detected NMR and NQR spectroscopy are not sensitive enough to detect the smaller number of spins in nanomaterials. The nitrogen-vacancy (NV) center in diamond has shown promise in bringing the analytic power of NMR and NQR spectroscopy to the nanoscale. However, due to depth-dependent formation efficiency of the defect centers, noise from surface spins, band bending effects, and the depth dependence of the nuclear magnetic field, there is ambiguity regarding the ideal NV depth for surface NMR of statistically polarized spins. In this work, we prepared a range of shallow NV ensemble layer depths and determined the ideal NV depth by performing NMR spectroscopy on statistically polarized 19 F in Fomblin oil on the diamond surface. We found that the measurement time needed to achieve a signal-to-noise ratio of 3 using XY8-N noise spectroscopy has a minimum at an NV ensemble depth of 5.5 ± 1.5 nm for ensembles activated from 100 ppm nitrogen concentration. To demonstrate the sensing capabilities of NV ensembles, we perform NQR spectroscopy on the 11 B of hexagonal boron nitride flakes. We compare our best diamond to previous work with a single NV and find that this ensemble provides a shorter measurement time with excitation diameters as small as 4 μm. Furthermore, this analysis provides ideal conditions for further experiments involving NMR/NQR spectroscopy of 2D materials with magnetic properties.
We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.