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.
Photoelastic measurements of residual strains are used extensively in the QC and inspection of transparent materials. A new method of measurements, based on Spectral Contents Analysis, is described in this paper. The method uses a personal computer for photoelastic data acquisition, eliminating personal skill and subjectivity. the new tool should make the measurements of residual strains for QC simpler and more reliable.
A method of computing the power spectral densities of the lateral response of airplanes to random atmospheric turbulence was adapted to an electronic digital computer. By use of this program, the power spectral densities of the lateral roll, yaw, and sideslip angular displacement of several conventional and STOL airplanes were computed. The results show that for the conventional airplanes, the roll response is more prominent than that for yaw or sideslip response. For the STOL airplanes, on the other hand, the yaw and sideslip responses were larger than the roll response. The response frequency of the STOL airplanes generally is higher than that for the conventional airplanes. This combination of greater sensitivity of the STOL airplanes in yaw and sideslip and the frequency at which they occur could be a factor causing the poor riding qualities of this class of airplanes.
A method for selecting the set of spectral characteristics that provides the smallest increase in prediction error is of interest to those using hyperspectral imaging (HSI) to monitor water quality. The spectral characteristics of interest to these applications are spectral bandwidth and location. Three water quality constituents of interest that are detectable via remote sensing are chlorophyll (CHL), total suspended solids (TSS), and colored dissolved organic matter (CDOM). Hyperspectral data provides a rich source of information regarding the content and composition of these materials, but often provides more data than an analyst can manage. This study addresses the spectral characteristics need for water quality monitoring for two reasons. First, determination of the greatest contribution of these spectral characteristics would greatly improve computational ease and efficiency. Second, understanding the spectral capabilities of different spectral resolutions and specific regions is an essential part of future system development and characterization. As new systems are developed and tested, water quality managers will be asked to determine sensor specifications that provide the most accurate and efficient water quality measurements. We address these issues using data from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and a set of models to predict constituent concentrations.
A method to remotely sense marine stratus clouds, that are underneath the boundary layer inversion produced by strong subsidence, is developed with the help of the measurements made by the Infrared Interferometer Spectrometer (IRIS) on board the Nimbus 4 satellite. These measurements, which had a spectral resolution of 2.8/cm and a field view of about 100 km in the region 870-980/cm, allow one to compare the weak absorption due to water vapor lines with that of the weak nu3-nu1 band of CO2. From such a comparison it is possible to sense the presence of these boundary layer stratus clouds. Maps of these clouds over subsidence regions of the tropical and subtropical oceans deduced from IRIS data show the usefulness of the technique to climatological studies.
On top of machine learning (ML) models, uncertainty quantification (UQ) functions as an essential layer of safety assurance that could lead to more principled decision making by enabling sound risk assessment and management. The safety and reliability improvement of ML models empowered by UQ has the potential to significantly facilitate the broad adoption of ML solutions in high-stakes decision settings, such as healthcare, manufacturing, and aviation, to name a few. In this tutorial, we aim to provide a holistic lens on emerging UQ methods for ML models with a particular focus on neural networks and the applications of these UQ methods in tackling engineering design as well as prognostics and health management problems. Towards this goal, we start with a comprehensive classification of uncertainty types, sources, and causes pertaining to UQ of ML models. Next, we provide a tutorial-style description of several state-of-the-art UQ methods: Gaussian process regression, Bayesian neural network, neural network ensemble, and deterministic UQ methods focusing on spectral-normalized neural Gaussian process. Established upon the mathematical formulations, we subsequently examine the soundness of these UQ methods quantitatively and qualitatively (by a toy regression example) to examine their strengths and shortcomings from different dimensions. Then, we review quantitative metrics commonly used to assess the quality of predictive uncertainty in classification and regression problems. Afterward, we discuss the increasingly important role of UQ of ML models in solving challenging problems in engineering design and health prognostics. In conclusion, two case studies with source codes available on GitHub are used to demonstrate these UQ methods and compare their performance in the life prediction of lithium-ion batteries at the early stage (case study 1) and the remaining useful life prediction of turbofan engines (case study 2).
Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.
The method of Maximum Likelihood (ML) is used to estimate the spectral parameters of an assumed broken power law energy spectrum from simulated detector responses. This methodology, which requires the complete specificity of all cosmic-ray detector design parameters, is shown to provide approximately unbiased, minimum variance, and normally distributed spectra information for events detected by an instrument having a wide range of commonly used detector response functions. The ML procedure, coupled with the simulated performance of a proposed space-based detector and its planned life cycle, has proved to be of significant value in the design phase of a new science instrument. The procedure helped make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope. This ML methodology is then generalized to estimate broken power law spectral parameters from real cosmic-ray data sets.
This study proposes a novel method for signal detection and feature extraction based on the spectral correlation function, enabling improved characterization of grid-signal distortions. Our approach differs from existing treatments of signal distortion in its analysis of the varied spectral content of signals observed in real-world scenarios. The method we propose has state-of-the-art discriminative power that provides meaningful and understandable characterizations of various grid events and anomalies. In conclusion, to validate the approach, we use real world data from the Grid Event Signature Library, which is maintained jointly by Oak Ridge National Laboratory and Lawrence Livermore National Laboratory.
A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.
Few field methods are available for characterizing source zones contaminated with aqueous film forming foam (AFFF). Non-invasive geophysical characterization of AFFF source zone contamination in situ could assist with the delineation and characterization of these sites, allowing for more informed sampling regimes aimed at quantifying subsurface poly- and perfluoroalkyl substances (PFAS) contamination. In this study, we present initial results from the investigation of the sensitivity of two existing surface and borehole-deployable geophysical technologies, spectral induced polarization (SIP), and low field nuclear magnetic resonance (NMR), to soils impacted with AFFF. To investigate the sensitivity of these methods to AFFF-impacted soil, bench-scale column experiments were conducted on samples consisting of natural and synthetic soils and groundwater. While our findings do not show strong evidence of NMR sensitivity to soil PFAS contamination, we do find evidence that SIP has sufficient sensitivity to detect PFAS sorption to soils. This finding is based on evidence that AFFF constituents associated with the pore surface produce a measurable polarization response in both freshly contaminated synthetic soils and in soils historically impacted with AFFF.
Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.
The Angstrom Absorption Exponent (AAE) is a parameter commonly used to characterize the wavelength-dependence of aerosol absorption optical depth (AAOD). It is closely related to aerosol composition. Black carbon (BC) containing aerosols yield AAE values near unity whereas Organic carbon (OC) aerosol particles are associated with values larger than 2. Even larger AAE values have been reported for desert dust aerosol particles. Knowledge of spectral AAOD is necessary for the calculation of direct radiative forcing effect of aerosols and for inferring aerosol composition. We have developed a satellitebased method of determining the spectral AAOD of absorbing aerosols. The technique uses multi-spectral measurements of upwelling radiation from scenes where absorbing aerosols lie above clouds as indicated by the UV Aerosol Index. For those conditions, the satellite measurement can be explained, using an approximations of Beer's Law (BL), as the upwelling reflectance at the cloud top attenuated by the absorption effects of the overlying aerosol layer. The upwelling reflectance at the cloud-top in an aerosol-free atmospheric column is mainly a function of cloud optical depth (COD). In the proposed method of AAE derivation, the first step is determining COD which is retrieved using a previously developed color-ratio based approach. In the second step, corrections for molecular scattering effects are applied to both the observed ad the calculated cloud reflectance terms, and the spectral AAOD is then derived by an inversion of the BL approximation. The proposed technique will be discussed in detail and application results making use of OMI multi-spectral measurements in the UV-Vis. will be presented.
Galaxy evolution has been studied by interpreting the spectral energy distribution of galaxies using spectral synthesis codes. This method has been crucial in discovering different pillars of modern galaxy evolution theories. However, this analysis was mostly carried out using spectral synthesis codes that are purely stellar, that is, they assume that the nebular contribution to the total continuum is negligible. The code FADO is the first publicly available population spectral synthesis tool that treats the contribution from ionised gas to the observed emission self-consistently. This is expected to have a particularly strong effect in star-forming (SF) galaxies. We study the impact of the nebular contribution on the determination of the star formation rate (SFR), stellar mass, and consequent effect on the star-forming main sequence (SFMS) at low redshift. We applied FADO to the spectral database of the SDSS to derive the physical properties of galaxies. As a comparison, we used the data in the MPA-JHU catalogue, which contains the properties of SDSS galaxies derived without the nebular contribution. We selected a sample of SF galaxies with Hα and Hβ flux measurements, and we corrected the fluxes for the nebular extinction through the Balmer decrement. We then calculated the Hα luminosity to estimate the SFR. Then, by combining the stellar mass and SFR estimates from FADO and MPA-JHU, the SFMS was obtained. The Hα flux estimates are similar between FADO and MPA-JHU. Because the Hα flux was used as tracer of the SFR, FADO and MPA-JHU agree in their SFR. The stellar mass estimates are slightly higher for FADO than for MPA-JHU on average. However, considering the uncertainties, the differences are negligible. With similar SFR and stellar mass estimates, the derived SFMS is also similar between FADO and MPA-JHU. Our results show that for SDSS normal SF galaxies, the additional modelling of the nebular contribution does not affect the retrieved fluxes and consequentially also does not influence SFR estimators based on the extinction-corrected Hα luminosity. For the stellar masses, the results point to the same conclusion. These results are a consequence of the fact that the vast majority of normal SF galaxies in the SDSS have a low nebular contribution. However, the obtained agreement might only hold for local SF galaxies, but higher-redshift galaxies might show different physical properties when FADO is used. This would then be an effect of the expected increased nebular contribution.
Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.
A method is presented for biasing spectral estimates to enhance detection of tonal signals against a background of broadband noise. In this method, a nonlinear average of an ensemble of individual spectral estimates is made where broadband noise energy is biased downward, pure tone energy is unbiased, and a mixture of the two is biased by an amount that depends on the ratio of tonal energy to broadband energy. The method is analyzed to provide estimates of the extent of tonal signal detection enhancement.
The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.
Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.
The infrared reflexion-spectra of minerals and rocks are used for remote sensing of targets. The reflexion-spectra of silicate rocks vary quite significantly from mineral to mineral in the wave length region from 8 to 12 micrometers. The rock forming minerals like quartz, feldspar, mica and the clay minerals show very different spectral shapes and positions of their maximum of the spectral reflexion. The presence of a good atmospherical window in that spectral region makes the method of differential-reflexion measurement feasible, for remote sensing application. A tunable CO2-laser was used as transmitter for infrared radiation. Laboratory tests showed the feasibility of the method under different simulated environmental conditions. Because of the very narrow bandwidth of the laser-emission lines, reflexion-spectra with extremely high spectral resolution were obtained.