Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time series 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Engine Icing Data - An Analytics Approach

Engine icing researchers at the NASA Glenn Research Center use the Escort data acquisition system in the Propulsion Systems Laboratory (PSL) to generate and collect a tremendous amount of data every day. Currently these researchers spend countless hours processing and formatting their data, selecting important variables, and plotting relationships between variables, all by hand, generally analyzing data in a spreadsheet-style program (such as Microsoft Excel). Though spreadsheet-style analysis is familiar and intuitive to many, processing data in spreadsheets is often unreproducible and small mistakes are easily overlooked. Spreadsheet-style analysis is also time inefficient. The same formatting, processing, and plotting procedure has to be repeated for every dataset, which leads to researchers performing the same tedious data munging process over and over instead of making discoveries within their data. This paper documents a data analysis tool written in Python hosted in a Jupyter notebook that vastly simplifies the analysis process. From the file path of any folder containing time series datasets, this tool batch loads every dataset in the folder, processes the datasets in parallel, and ingests them into a widget where users can search for and interactively plot subsets of columns in a number of ways with a click of a button, easily and intuitively comparing their data and discovering interesting dynamics. Furthermore, comparing variables across data sets and integrating video data (while extremely difficult with spreadsheet-style programs) is quite simplified in this tool. This tool has also gathered interest outside the engine icing branch, and will be used by researchers across NASA Glenn Research Center. This project exemplifies the enormous benefit of automating data processing, analysis, and visualization, and will help researchers move from raw data to insight in a much smaller time frame.

Engine Icing↗
DOE OSTI2023

Spatiotemporally Adaptive Compression for Scientific Dataset with Feature Preservation – A Case Study on Simulation Data with Extreme Climate Events Analysis

Scientific discoveries are increasingly constrained by limited storage space and I/O capacities. For time-series simulations and experiments, their data often need to be decimated over timesteps to accommodate storage and I/O limitations. In this paper, we propose a technique that addresses storage costs while improving post-analysis accuracy through spatiotemporal adaptive, error-controlled lossy compression. We investigate the trade-off between data precision and temporal output rates, revealing that reducing data precision and increasing timestep frequency lead to more accurate analysis outcomes. Additionally, we integrate spatiotemporal feature detection with data compression and demonstrate that performing adaptive error-bounded compression in higher dimensional space enables greater compression ratios, leveraging the error propagation theory of a transformation-based compressor. To evaluate our approach, we conduct experiments using the well-known E3SM climate simulation code and apply our method to compress variables used for cyclone tracking. Our results show a significant reduction in storage size while enhancing the quality of cyclone tracking analysis, both quantitatively and qualitatively, in comparison to the prevalent timestep decimation approach. Compared to three state-of-the-art lossy compressors lacking feature preservation capabilities, our adaptive compression framework improves perfectly matched cases in TC tracking by 26.4-51.3% at medium compression ratios and by 77.3-571.1% at large compression ratios, with a merely 5–11% computational overhead.

Gong, Qian↗
DOE OSTI2024

Learning continuous scattering length density profiles from neutron reflectivities using convolutional neural networks

Interpreting neutron reflectivity (NR) data using ad hoc multi-layer models and physics-based models provides information about spatially resolved neutron scattering length density (NSLD) profiles. Recent improvements in data acquisition systems have allowed acquiring thousands of NR curves in a couple of hours, which has led to a need for automated data analysis tools to interpret NR measurements in real-time. Here, we present a machine learning analysis workflow that uses a series of models, based on a convolutional neural network (CNN), to learn the relation between the NSLDs and the NRs, and subsequently produce continuous NSLD profiles directly from NRs. The usefulness of our CNN-based models is demonstrated by constructing NSLDs from NRs of several films containing homopolymer polyzwitterions and diblock copolymers mixed with different types of salts. Comparisons of the NSLDs with those constructed using ad hoc multi-layer models reveal a very good agreement, suggesting the potential of CNN-based models for real-time automated data analysis of NRs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

clif

SAND2022-12903 O clif contains code to perform climate fingerprinting. It includes capabilities for principal component analysis of climate data, Fourier analysis, a collection of pre-processing scripts, and time series statistical tools. It uses a templated object-oriented programming approach to integrate all classes with the operator in Scikit Learn. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Nichol, Jeffrey↗
DOE OSTI2022

When Do Different Scenarios of Future Electricity Demand Start to Meaningfully Diverge?

Climate and population change are known to influence electricity demand, but what is the impact of uncertainty in climate and population projections on electricity demand in the United States (U.S.) over the next 30 years? This question has important implications for investment decisions in the energy sector, which are typically made using a 15- to 30-year time horizon. Simply put, if future climate and population scenarios do not lead to meteorological conditions and subsequent energy demands that are distinctly different within the first 30 years then, for the purposes of investment decisions, it may not matter which potential pathway we are most likely on. The Integrated, Multisector, Multiscale Modeling (IM3) project has generated a wide, yet plausible range of 21st century high-resolution climate and population scenarios for the U.S. The IM3 projections span two climate scenarios (Representative Concentration Pathways 4.5 and 8.5) and two population scenarios (Shared Socioeconomic Pathways 3 and 5). For each of the climate scenarios, we reflect a range of climate model uncertainty by using boundary conditions from a set of climate models that are hotter and colder than the multi-model mean. In total there are eight (2 x 2 x 2) IM3 scenarios. This work explores a basic question: When do future climate and energy outcomes start to meaningfully diverge across the eight scenarios? The analysis is done on both the raw meteorology time-series from each scenario as well as a time-series of projected total electricity demand generated using IM3's Total ELectricity Loads (TELL) model. We also explore the spatial heterogeneity of the divergence signal to identify regions of the country that may be more or less path dependent.

Climate Change↗
NASA NTRS1988

Microwave remote sensing of sea ice

The long term objectives are: (1) to understand the physics of the multispectral microwave radiative characteristics of sea ice as it goes through different phases; (2) to improve characterization of sea ice cover using satellite microwave sensors; and (3) to study ice/ocean physical and biological processes associated with polynya formations and variability of the marginal sea ice region. Two field experiments were conducted to pursue these objectives. One involved measurements of radiative and physical characteristics of sea ice from a ship during a 3-month long cruise through the Weddell Sea ice pack during the Austral winter of 1986. The other involved similar measurements from two aircrafts and a submarine over the Central Arctic and Greenland Sea region. Preliminary results have already led to an enhanced understanding of the microwave signatures of pancake ice, nilas, first year ice, multiyear ice and effects of snow cover. Coastal and deep ocean polynyas and their role in bottom water formation and ocean circulation were studied using a time series of ice images from SMMR. An unsupervised cluster analysis of Arctic sea ice using SMMR and THIR emissivity and brightness temperature data was implemented. The analysis indicates the existence of several unique and persistent clusters in the Central Arctic region during winter and that the sum of the area of these clusters excluding those of first year ice is about 20 percent less than minimum ice cover area inferred from a previous summer data. This result is consistent with saline surface for some multiyear ice floes as observed during MIZEZ and suggests that a significant fraction of multiyear ice floes in the Arctic have first year ice signatures.

Comiso, J. C.↗
NASA NTRS1989

The role of refractive interstellar scintillation in the low frequency variability of extragalactic radio sources

A statistical analysis is presented of the 408 MHz flux density time series of 39 extragalactic radio sources, to study the role of refractive interstellar scintillation in the decimeter wavelength variability of compact extragalactic radio sources. It is found that the product of the scintillation index and the variability timescale is systematically larger for sources at lower galactic latitudes. This variation is obscured by much larger random variations. Possible explanations of these results are discussed.

Spangler, S.↗
NASA NTRS1993

Reconstruction of low-dimensional magnetospheric dynamics by singular spectrum analysis

Singular spectrum analysis was used to study global magnetospheric dynamics using time series auroral electrojet (AE) data. The calculated correlation dimension of approximately 2.5 confirms the low dimensionality of earlier results. This technique indicated that global magnetospheric dynamics can be described by three variables whose dynamical features were obtained from the AE data.

Sharma, A. S.↗
NASA NTRS1991

Sensor fusion techniques for predicting thunderstorm evolution using lightning and radar networks

A pattern recognition scheme that is used to generate initial seeds or 'first guess' fields for clustering the discrete lightning discharges into storm cells is presented. The clustering process is critically dependent on the prior accuracy of the lightning location estimates. The generation of subsequent storm life-cycle time series (for extrapolation forecasting) also relies on the cluster analysis procedure assigning the correct number of lightning discharges (objects) to the proper storms (groups). The advantages and limitations of different clustering strategies for storm identification and tracking are examined. Storm identification with lightning data alone is compared to storm identification with radar alone, and some synergies for sensor fusion are explored.

Goodman, Steven J.↗
NASA NTRS2000

Estimation of EOP From VLBI: Direct Approach

The currently adopted strategy of Earth Orientation Parameters (EOP) estimation from Very Long Baseline Interferometry (VLBI) is to estimate six parameters: Universal Time 1 (UT1), UT1 rate, pole positions, and nutation offsets for each 24-hour session independently. Then the resulting time series of raw EOP are filtered and a regression analysis is performed to obtain nutation coefficients, polhode of the pole, and other physical parameters. Thus, the latter parameters are obtained indirectly in two stages. An alternative approach of direct estimation of the final EOP is presented. Pole coordinates and UT1 are considered as a sum of three components: the low-period component that is modeled by a cubic spline, the harmonic component that includes forced nutation, precession and sub-daily variations of EOP, and the stochastic component that is modeled by a linear spline with segment length 1-2 hours. All parameters are obtained in a single LSQ solution using all available data.

POLAR MOTION↗
NASA NTRS2007

Dynamic and Regression Modeling of Ocean Variability in the Tide-Gauge Record at Seasonal and Longer Periods

Comparison of monthly mean tide-gauge time series to corresponding model time series based on a static inverted barometer (IB) for pressure-driven fluctuations and a ocean general circulation model (OM) reveals that the combined model successfully reproduces seasonal and interannual changes in relative sea level at many stations. Removal of the OM and IB from the tide-gauge record produces residual time series with a mean global variance reduction of 53%. The OM is mis-scaled for certain regions, and 68% of the residual time series contain a significant seasonal variability after removal of the OM and IB from the tide-gauge data. Including OM admittance parameters and seasonal coefficients in a regression model for each station, with IB also removed, produces residual time series with mean global variance reduction of 71%. Examination of the regional improvement in variance caused by scaling the OM, including seasonal terms, or both, indicates weakness in the model at predicting sea-level variation for constricted ocean regions. The model is particularly effective at reproducing sea-level variation for stations in North America, Europe, and Japan. The RMS residual for many stations in these areas is 25-35 mm. The production of "cleaner" tide-gauge time series, with oceanographic variability removed, is important for future analysis of nonsecular and regionally differing sea-level variations. Understanding the ocean model's strengths and weaknesses will allow for future improvements of the model.

Hill, Emma M.↗
DOE OSTI2026

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗
NASA NTRSDate not supplied

Impact of DORISx Post Seismic Deformation and Seasonal Corrections on the DPOD2020 Solution

As one of the tracking systems used for the altimeter missions (such as TOPEX/Poseidon, Envisat, Jason-1/-2/-3, CryoSat-2, Saral/AltiKa, Sentinel-3A/-3B, HY-2A/C/D, Jason-CS/Sentinel-6A, SWOT), the position of the DORIS tracking stations provides a fundamental reference for the estimation of the precise orbits and so, by extension is fundamental for the quality of the altimeter data analysis and derived products. Due to the time evolution of the DORIS ground network, stations included in ITRF2020 may have been decommissioned and stations that were added to the tracking network after 2021.0 are, by definition, not part of ITRF2020. Therefore, to satisfy operational requirements for precise orbit determination and routine delivery of geodetic products, the International DORIS Service maintains the DORIS terrestrial reference frame for Precise Orbit Determination (DPOD) solutions. Since 2016, the DPOD products are computed by the IDS Combination Center (CC) as a DORIS cumulative position and linear velocity solution aligned to the latest available ITRF and using the latest IDS weekly combined series. Mid-2023, the IDS CC started the computation of the second version of the DPOD2020 based on the IDS 20 weekly solutions from 1993.0 to 2023.0. With that new solution, the IDS CC will include Post Seismic Deformation (PSD) corrections for some DORIS sites such as Socorro Island. These DPOD PSD corrections are determined from the analysis of DORIS coordinate time series only. The second version will also include optional annual and semi-annual station corrections. The presentation will describe the analysis used for extracting the PSD and seasonal DORIS station corrections, compare such corrections to those derived for ITRF2020, and evaluate the impact these new corrections have on station positioning and satellite POD.

Guilhem Moreaux↗
NASA NTRSDate not supplied

Spectral Analysis of Integrated Pressures on Patches with Unsteady Pressure-Sensitive Paint Measurements

This paper describes the spectral analysis of integrated pressures on patches of the scale model of the Space Launch System (SLS) Block 1B crew vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the ascent transient aerodynamics tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. Recent research has demonstrated that uPSP can be an essential tool for the assessment of the unsteady, aerodynamic phenomena. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. In this paper, 108 patches are defined with x station values and azimuth angles of the scale model. For each patch, the polygons are determined from the surface cells of the grid of the model, clipped with the edges of the patch, and each of the polygons is divided into triangles. The inputs of the pressure integration algorithm are the time series of pressure coefficients on the vertices of the grid of the model, generated by the uPSP data processing program from the videos taken with 4 Phantom high-speed cameras during the wind tunnel tests. The integrated pressure of the patch is determined as the ratio of the sum of the forces on the triangles over the sum of the areas of the triangles. For each of the test cases, the Cross Power Spectral Density (CPSD) and Magnitude-Squared Coherence (msCohere) are computed for frequencies over 1/3 octave bands from the time series of the integrated pressures on the patches. The spectral analysis outputs of different cases of the wind tunnel tests are analyzed and the coherence of patches is investigated over frequencies, x station values, azimuth angles, Mach numbers and SLS vehicle attitudes. The pressure integration and spectral analysis described in this paper were executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aeroscience Evaluation and Test Capabilities (AETC) Project.

Pressure-Sensitive Paint↗
NASA NTRS2010

GOATS - Orbitology Component

The GOATS Orbitology Component software was developed to specifically address the concerns presented by orbit analysis tools that are often written as stand-alone applications. These applications do not easily interface with standard JPL first-principles analysis tools, and have a steep learning curve due to their complicated nature. This toolset is written as a series of MATLAB functions, allowing seamless integration into existing JPL optical systems engineering modeling and analysis modules. The functions are completely open, and allow for advanced users to delve into and modify the underlying physics being modeled. Additionally, this software module fills an analysis gap, allowing for quick, high-level mission analysis trades without the need for detailed and complicated orbit analysis using commercial stand-alone tools. This software consists of a series of MATLAB functions to provide for geometric orbit-related analysis. This includes propagation of orbits to varying levels of generalization. In the simplest case, geosynchronous orbits can be modeled by specifying a subset of three orbit elements. The next case is a circular orbit, which can be specified by a subset of four orbit elements. The most general case is an arbitrary elliptical orbit specified by all six orbit elements. These orbits are all solved geometrically, under the basic problem of an object in circular (or elliptical) orbit around a rotating spheroid. The orbit functions output time series ground tracks, which serve as the basis for more detailed orbit analysis. This software module also includes functions to track the positions of the Sun, Moon, and arbitrary celestial bodies specified by right ascension and declination. Also included are functions to calculate line-of-sight geometries to ground-based targets, angular rotations and decompositions, and other line-of-site calculations. The toolset allows for the rapid execution of orbit trade studies at the level of detail required for the early stage of mission concept development.

Haber, Benjamin M.↗
NASA NTRS2016

SEEDS Direct Imaging of the RV-Detected Companion to V450 Andromedae, and Characterization of the System

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity (RV) method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and RV measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be m(sub 1) = 1.141(sup + 0.037, sub - 0.091) and m(sub 2) = 0.279(sup + 0.023, sub - 0.020) solar mass. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. Hipparcos time-series photometry shows a periodicity of P = 5.743 day, which is also seen in the SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is 380(sup + 220, sub -100) Myr old, consistent with an isochrone analysis (220 (sup +2120, sub -90) Myr). This work may serve as a test case for future studies of low mass stars, brown dwarfs, and exoplanets by combination of RV and direct imaging data.

direct imaging detection of a low-mass companion↗