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 method of synthesizing an arbitrarily shaped transient time pulse on vibration exciters is described. The transient waveform control technique is based on recent developments in digital time series analysis, the real-time FFT processor. A brief description of the theory, error estimates and hardware/software implementation to the JPL Dynamic Environmental Testing Laboratory is presented.
The development of vector processing computers of the streaming array architectures has made possible a dramatic decrease in the time required for the solution of problems, i.e., for those algorithms which readily lend themselves to sequential operations on long vectors. There has concurrently been rapid growth in the applications of the techniques generally known as mathematical digital filtering, also called signal analysis, digital signal processing, time series analysis, and digital image processing. The best known applications of these techniques is to seismic data and two-dimensional images, data types consisting of very large collections of numbers. A major limitation for these cases is the size and speed of the computer available. Therefore, a move to the streaming array architecture can result in marked improvement in the data analysis techniques which may be employed.
The implications of using different control variables for the analysis of moisture observations in a global atmospheric data assimilation system are investigated. A moisture analysis based on either mixing ratio or specific humidity is prone to large extrapolation errors, due to the high variability in space and time of these parameters and to the difficulties in modeling their error covariances. Using the logarithm of specific humidity does not alleviate these problems, and has the further disadvantage that very dry background estimates cannot be effectively corrected by observations. Relative humidity is a better choice from a statistical point of view, because this field is spatially and temporally more coherent and error statistics are therefore easier to obtain. If, however, the analysis is designed to preserve relative humidity in the absence of moisture observations, then the analyzed specific humidity field depends entirely on analyzed temperature changes. If the model has a cool bias in the stratosphere this will lead to an unstable accumulation of excess moisture there. A pseudo-relative humidity can be defined by scaling the mixing ratio by the background saturation mixing ratio. A univariate pseudo-relative humidity analysis will preserve the specific humidity field in the absence of moisture observations. A pseudorelative humidity analysis is shown to be equivalent to a mixing ratio analysis with flow-dependent covariances. In the presence of multivariate (temperature-moisture) observations it produces analyzed relative humidity values that are nearly identical to those produced by a relative humidity analysis. Based on a time series analysis of radiosonde observed-minus-background differences it appears to be more justifiable to neglect specific humidity-temperature correlations (in a univariate pseudo-relative humidity analysis) than to neglect relative humidity-temperature correlations (in a univariate relative humidity analysis). A pseudo-relative humidity analysis is easily implemented in an existing moisture analysis system, by simply scaling observed-minus background moisture residuals prior to solving the analysis equation, and rescaling the analyzed increments afterward.
This method is especially well suited for analyzing time-series data that represent nonstationary and nonlinear physical phenomena. The method is based principally on the concept of empirical mode decomposition (EMD), according to which any complicated signal (as represented by digital samples) can be decomposed into a finite number of functions, called "intrinsic mode functions" (IMFs), that admit well-behaved Hilbert transforms. The local energies and the instantaneous frequencies derived from the IMFs through Hilbert transforms can be used to construct an energy-frequency-time distribution, denoted a Hilbert spectrum.
Orthostatic Intolerance (OI) is the propensity to develop symptoms of fainting during upright standing. OI is associated with changes in heart rate, blood pressure and other measures of cardiac function. Problem: NASA astronauts have shown increased susceptibility to OI on return from space missions. Current methods for counteracting OI in astronauts include fluid loading and the use of compression garments. Multivariate trajectory spread is greater as OI increases. Pairwise comparisons at the same time within subjects allows incorporation of pass/fail outcomes. Path length, convex hull area, and covariance matrix determinant do well as statistics to summarize this spread Missing data problems Time series analysis need many more time points per OTT session treatment of trend? how incorporate survival information?
Explore the source record for details and available documents.
Earth Observation data are a vital resource for studying long term changes, but the large data volumes can be challenging to analyze. Time series analysis in particular is hampered by the typical thin-time-slice file organization. We examine several potential solutions inspired in large part by the data-parallel methods that have arisen with cloud computing. These solutions include various combinations of data re-organization, spatial indexing, distributed storage and pre-computation that we term "Analytics Optimized Data Stores" (AODS). We find that even simple solutions (such as a data cube) produce more than an order of magnitude improvement; the best provide two to three orders of magnitude improvement. The most performant solutions have tradeoffs in terms of generality or storage footprint, but may nonetheless be useful components in data analytics frameworks where performance is critical.
Results from Landsat satellite image times series analysis since 1983 of this study area showed gradual, statistically significant increases in the normalized difference vegetation index (NDVI) in more than 90% of the (predominantly second-growth) evergreen forest locations sampled.
The Framework for Optimization of Resources and Economics (FORCE) tool suite is the U.S. Department of Energy’s Nuclear Integrated Energy Systems (IES) Program flagship tool suite for technoeconomic IES analysis of IES. This tool suite is useful for analysis designed to evaluate and improve the technoeconomics of energy production systems, particularly for systems including nuclear technology. In this report, we document the development activity for the FORCE tool suite to extend its capabilities as performed during fiscal year 2024. In addition to reliability and accessibility, capability is one of the three standards guiding the development of the FORCE tool suite and the software codes that are its constituent parts. Extending the capabilities of the FORCE tool suite allows analysis both within the IES program as well as industry, university, and laboratory partners to perform analysis with more accuracy, insight, and impactful narrative. Four areas of capability development were the focus of activity this year: economic parameter uncertainty quantification, multiresolution analysis, components-to-optimization workflow automation, and statespace construction workflows for real-time optimal control. In economic parameter uncertainty quantification, the ability of HERON to capture risk due to scenarios (weather and energy demand uncertainty) was expanded to also include uncertainties in financial parameters such as capital cost or operation and maintenance costs. By including these sources of uncertainty, which are sometimes very large compared with scenario uncertainty, HERON is better able to capture the risk posed by investment in various IES technology. Because of this, analysts can also consider the reduction in risks that can be realized by choice of some technologies. In multiresolution analysis, development activity extended on work completed previously. In fiscal year 2023, methods for decomposing time series signals, such as demand, solar and wind availability, and price profiles, were analyzed and down-selected to those most effective at splitting signals into different resolutions. These resolutions allow considering the influence of different energy demand and supply behaviors across different time scales. For example, energy demand might be divided into seasonal, weekly, and hourly profiles. In fiscal year 2024, this preliminary work was extended and implemented within the Risk Analysis Virtual Environment (RAVEN) risk and uncertainty analysis platform, which is used throughout the FORCE framework. This development of the “multi-resolution time series analysis” (MR-TSA) module in RAVEN allows training synthetic history generators on complex time series. These synthetic history generators can then be used in HERON for generating scenarios that represent possible market and weather scenarios that can be analyzed on different time scales. We envision completing this work in the future, implementing multiresolution dispatch optimization strategies that can make the most beneficial use of these stratified time histories. In components-to-optimization workflow development, workflows for translating user inputs of components into algorithms for algebraic optimization were selected and implemented. Similar algorithms within the Holistic Energy Resource Optimization Network (HERON) were separated from the main code base of HERON and gathered with the components-to-optimization workflows in the new Dispatch Optimization Variable Engine (DOVE) software library. This modularization allows FORCE users to analyze dispatch optimization and energy system duty cycles independently of HERON, which previously was a burdensome task. Additionally, these dispatch optimization algorithms, set up in an independent library, can now be used across all software applications within FORCE, especially including the real-time optimal control software Optimization of Real-time Capacity Allocation (ORCA). Allowing FORCE software to share dispatch optimization algorithms within a single library allows for improved software maintenance and reliability. In statespace characterization workflow development, alternative workflows for optimizing dispatch with additional technical accuracy was the focus, particularly to improve the real-time optimization decision making in ORCA. Using algorithms and workflows initially developed for the Feasible Actuator Range Modifier (FARM), workflows for determining the statespace representation of IES were identified and demonstrated. The resulting dispatch optimization required a more robust optimization algorithm than that originally used in HERON (and moved to DOVE), which required adding an alternate workflow to DOVE that can more accurately match the behavior of physical systems using a partial differential equation representation. In conclusion, capability developments in the FORCE tool suite in fiscal year 2024 have improved the ability of the FORCE tool suite to perform
This paper presents a new automatic algorithm to derive optical snow grain size (SGS) at 1 km resolution using Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Differently from previous approaches, snow grains are not assumed to be spherical but a fractal approach is used to account for their irregular shape. The retrieval is conceptually based on an analytical asymptotic radiative transfer model which predicts spectral bidirectional snow reflectance as a function of the grain size and ice absorption. The analytical form of solution leads to an explicit and fast retrieval algorithm. The time series analysis of derived SGS shows a good sensitivity to snow metamorphism, including melting and snow precipitation events. Preprocessing is performed by a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which includes gridding MODIS data to 1 km resolution, water vapor retrieval, cloud masking and an atmospheric correction. MAIAC cloud mask (CM) is a new algorithm based on a time series of gridded MODIS measurements and an image-based rather than pixel-based processing. Extensive processing of MODIS TERRA data over Greenland shows a robust performance of CM algorithm in discrimination of clouds over bright snow and ice. As part of the validation analysis, SGS derived from MODIS over selected sites in 2004 was compared to the microwave brightness temperature measurements of SSM\I radiometer, which is sensitive to the amount of liquid water in the snowpack. The comparison showed a good qualitative agreement, with both datasets detecting two main periods of snowmelt. Additionally, MODIS SGS was compared with predictions of the snow model CROCUS driven by measurements of the automatic whether stations of the Greenland Climate Network. We found that CROCUS grain size is on average a factor of two larger than MODIS-derived SGS. Overall, the agreement between CROCUS and MODIS results was satisfactory, in particular before and during the first melting period in mid-June. Following detailed time series analysis of SGS for four permanent sites, the paper presents SGS maps over the Greenland ice sheet for the March-September period of 2004.
Bayesian Blocks are optimal piecewise linear representations (step function fits) of light-curves. The simple algorithm implementing this idea, using dynamic programming, has been extended to include more data modes and fitness metrics, multivariate analysis, and data on the circle (Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations, Scargle, Norris, Jackson and Chiang 2013, ApJ, 764, 167), as well as new results on background subtraction and refinement of the procedure for precise timing of transient events in sparse data. Example demonstrations will include exploratory analysis of the Kepler light curve archive in a search for "star-tickling" signals from extraterrestrial civilizations. (The Cepheid Galactic Internet, Learned, Kudritzki, Pakvasa1, and Zee, 2008, arXiv: 0809.0339; Walkowicz et al., in progress).
Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.
This is a tutorial on how to use the SYSAI (System Analysis using Statistical AI), a flexible statistical learning framework for the V&V and analysis of complex and high-dimensional Aerospace systems with DNN and AI components. SYSAI provides functionality for a variety of analyses and V&V tasks, including statistical data analysis, high dimensional safety-envelope and time-series analysis, property checking, as well as intelligent test-case generation. The tutorial will demonstrate SYSAI with our industrial partner’s Autonomous Centerline Tracking system, which uses a DNN to enable autonomous aircraft taxiing as an example. Video & Tutorial
The software for configuring a Laser Velocimeter (LV) counter processor system was developed using structured design. The LV system includes up to three counter processors and a rotary encoder. The software for configuring and testing the LV system was developed, tested, and included in an overall software package for data acquisition, analysis, and reduction. Error handling routines respond to both operator and instrument errors which often arise in the course of measuring complex, high-speed flows. The use of networking capabilities greatly facilitates the software development process by allowing software development and testing from a remote site. In addition, high-speed transfers allow graphics files or commands to provide viewing of the data from a remote site. Further advances in data analysis require corresponding advances in procedures for statistical and time series analysis of nonuniformly sampled data.
The software for configuring an LV counter processor system has been developed using structured design. The LV system includes up to three counter processors and a rotary encoder. The software for configuring and testing the LV system has been developed, tested, and included in an overall software package for data acquisition, analysis, and reduction. Error handling routines respond to both operator and instrument errors which often arise in the course of measuring complex, high-speed flows. The use of networking capabilities greatly facilitates the software development process by allowing software development and testing from a remote site. In addition, high-speed transfers allow graphics files or commands to provide viewing of the data from a remote site. Further advances in data analysis require corresponding advances in procedures for statistical and time series analysis of nonuniformly sampled data.
The IUE spectra of bright F-K dwarfs were used to investigate rotational modulation in chromospheric and transition region emission lines, with Stellingwerf's phase dispersion minimization (PDM) technique as the basis of time series analysis. Sampling-related aliases are found to dominate the PDM spectra in these sparsely sampled, irregularly spaced time series. The star chosen to illustrate the techniques is epsilon Eridani, a young, chromospherically active K2 dwarf. Strong indications of a rotation period close to 2.8 days in epsilon Eridani from 2 seasons of observations, with evidence for the persistence in longitude of the active area(s) over several hundred rotation cycles are found.
Orbital signals are being discovered in pre-Pleistocene sediments. Due to their hierarchical nature these cycle patterns are complex, and the imprecision of geochronology generally makes the assignment of stratigraphic cycles to specific orbital cycles uncertain, but in sequences such as the limnic Newark Group under study by Olsen and pelagic Cretaceous sequence worked on by our Italo-American group the relative frequencies yield a definitive match to the Milankovitch hierarchy. Due to the multiple ways in which climate impinges on depositional systems, the orbital signals are recorded in a multiplicity of parameters, and affect different sedimentary facies in different ways. In platform carbonates, for example, the chief effect is via sea-level variations (possibly tied to fluctuating ice volume), resulting in cycles of emergence and submergence. In limnic systems it finds its most dramatic expression in alternations of lake and playa conditions. Biogenic pelagic oozes such as chalks and the limestones derived from them display variations in the carbonate supplied by planktonic organisms such as coccolithophores and foraminifera, and also record variations in the aeration of bottom waters. Whereas early studies of stratigraphic cyclicity relied mainly on bedding variations visible in the field, present studies are supplementing these with instrumental scans of geochemical, paleontological, and geophysical parameters which yield quantitative curves amenable to time-series analysis; such analysis is, however, limited by problems of distorted time-scales. My own work has been largely concentrated on pelagic systems. In these, the sensitivity of pelagic organisms to climatic-oceanic changes, combined with the sensitivity of botton life to changes in oxygen availability (commonly much more restricted in the Past than now) has left cyclic patterns related to orbital forcing. These systems are further attractive because (1) they tend to offer depositional continuity, and (2) presence of abundant microfossils yields close ties to geochronology. A tantalizing possibility that stratigraphy may yield a record of orbital signals unrelated to climate has turned up in magnetic studies of our Cretaceous core. Magnetic secular variations here carry a strong 39 ka periodicity, corresponding to the theoretical obliquity period of that time - Does the obliquity cycle perhaps have some direct influence on the magnetic field?
Time-sequenced IUE spectra of a number of a bright F-K dwarfs have been used to investigate rotational modulation in several important chromospheric and transition emission lines. Stellingwerf's phase dispersion minimization (PDM) technique forms the basis of this time-series analysis of these sparsely sampled, irregularly spaced observational sequences, which consist of no more than 13 spectra in any one season. In such highly sample-limited time series, the resultant PDM spectrum is greatly complicated by overlapping alias orders which invariably produce strong interference effects, but the deliberately irregular sampling has proven of value in the development of techniques to identify the intrinsic stellar flux modulation period. Two stars are selected to illustrate these new techniques in detail.