Search NASA⌕ Search

SEARCH · Search NASA

Results for “trajectory data”

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 181 records · Page 10

A Geometry-Driven Longitudal Topic Model

A simple and scalable framework for longitudinal analysis of Twitter data is developed that combines latent topic models with computational geometric methods. Dimensionality reduction tools from computational geometry are applied to learn the intrinsic manifold on which the latent, temporal topics reside. Then shortest path distances on the manifold are used to link together these topics. The proposed framework permits visualization of the low-dimensional embedding which provides clear interpretation of the complex, high-dimensional trajectories that may exist among latent topics. Practical application of the proposed framework is demonstrated through its ability to capture and effectively visualize natural progression of latent COVID-19 related topics learned from Twitter data. Interpretability of the trajectories is achieved by comparing to real-world events. In addition, the framework permits study of spatial variation in Twitter behavior for learned topics. The analysis demonstrates that the proposed framework is able to capture granular-level impact of COVID-19 on public discussions. We end by arguing that Twitter data, when analyzed within the proposed framework, can serve as a valuable supplementary data stream for COVID-related studies.

97 MATHEMATICS AND COMPUTING↗

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probabilistic error estimation for non-intrusive reduced models learned from data of systems governed by linear parabolic partial differential equations

This work derives a residual-based a posteriori error estimator for reduced models learned with non-intrusive model reduction from data of high-dimensional systems governed by linear parabolic partial differential equations with control inputs. It is shown that quantities that are necessary for the error estimator can be either obtained exactly as the solutions of least-squares problems in a non-intrusive way from data such as initial conditions, control inputs, and high-dimensional solution trajectories or bounded in a probabilistic sense. Here, the computational procedure follows an offline/online decomposition. In the offline (training) phase, the high-dimensional system is judiciously solved in a black-box fashion to generate data and to set up the error estimator. In the online phase, the estimator is used to bound the error of the reduced-model predictions for new initial conditions and new control inputs without recourse to the high-dimensional system. Numerical results demonstrate the workflow of the proposed approach from data to reduced models to certified predictions.

97 MATHEMATICS AND COMPUTING↗

Nonplanar effects in simulations of laser-driven ejecta microjet experiments

Recent experiments of laser-driven ejecta microjets performed at OMEGA 60 reveal tortuous jets whereby the jets appear to deviate from their initial trajectory as they travel in vacuum. To understand these data, we perform two dimensional numerical simulations, considering different target thicknesses, pressures, and models of the drive conditions. In particular, modeling the finite laser spot size appears essential in reproducing qualitatively the non-planar shock observed in the experiment. Simulations capture jet deflection by accounting for a slight misalignment of the laser pointing with respect to the groove axis along with spatial variation of the laser pulse intensity. The principal physical mechanism appears to be that lateral momentum is imparted by release waves arising from the non-planar drive. The induced off-axis velocity is small in comparison to the jet axial velocity but integrates into a pronounced deflection over the course of the experiment. The analysis of jet axial and lateral mass distributions is found to be reproduced reasonably by the simulations. Simulated radiographs are in qualitative agreement with the experiments, though their differences point to potential shortcomings in modeling strictly three-dimensional experiments using two-dimensional hydrodynamic simulations. The simple analysis is able to explain part of the observed discrepancy in simulated vs experimental jet masses.

42 ENGINEERING↗

Time-resolved vacuum-ultraviolet photoelectron spectroscopy of the à 1 A u state of acetylene

Ultrafast time-resolved photoelectron spectra are reported for the vacuum-ultraviolet (VUV) photoionization of acetylene following excitation to the à 1 A u state via UV absorption at 200 nm. The excitation energy lies above the lowest dissociation threshold to C 2 H X̃ 2 Σ + + H, as well as above the threshold for adiabatic dissociation of the à 1 A u state to form C 2 H (à 2 Π) + H. The time-dependent mass spectra and photoelectron spectra provide insight into the intramolecular decay processes of the à 1 A u state. In addition, photoelectron spectra of the à 1 A u state with VUV light access both the X̃ 2 Π u and à 2 Σ g + states of the ion, as well as the predicted, but previously unobserved, 1 2 Π g state, which corresponds to a two-hole, one-particle configuration that lies in close proximity to the à 2 Σ g + state. The 1 2 Π g state is split into 2 A 2 + 2 B 2 and 2 A g + 2 B g states in the cis and trans configurations, respectively. In conclusion, electronic structure calculations, along with trajectory calculations, reproduce the principal features of the experimental data and confirm the assignment of the 1 2 Π g state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimated impacts of forest restoration scenarios on smoke exposures among outdoor agricultural workers in California

Abstract As wildfires continue to worsen across western United States, forest managers are increasingly employing prescribed burns as a way to reduce excess fuels and future wildfire risk. While the ecological benefits of these fuel treatments are clear, little is known about the smoke exposure tradeoffs of using prescribed burns to mitigate wildfires, particularly among at-risk populations. Outdoor agricultural workers are a population at increased risk of smoke exposure because of their time spent outside and the physical demands of their work. Here, we assess the smoke exposure impacts among outdoor agricultural workers resulting from the implementation of six forest management scenarios proposed for a landscape in the Central Sierra, California. We leverage emissions estimates from LANDIS-II to model daily PM 2.5 concentrations with the Hybrid Single Particle Lagrangian Integrated Trajectory Model (HYSPLIT) and link those to agricultural employment data from the Bureau of Labor Statistics. We find a u-shaped result, in that moderate amounts of prescribed burning result in the greatest reduction in total smoke exposure among outdoor agricultural workers, particularly during months of peak agricultural activity due to wildfire-specific smoke reductions. The reduction in total smoke exposure, relative to scenarios with the least amount of management, decreases as more prescribed burning is applied to the landscape due to the contributions of the fuel treatments themselves to overall smoke burden. The results of this analysis may contribute to preparedness efforts aimed at reducing smoke exposures among outdoor agricultural workers, while also informing forest management planning for this specific landscape.

Environmental Sciences & Ecology↗

Advanced Intra-Cycle Detection of Pre-Ignition Events through Phase-Space Transforms of Cylinder Pressure Data

The widespread adoption of boosted, downsized SI engines has brought pre-ignition phenomena into greater focus, as the knock events resulting from pre-ignitions can cause significant hardware damage. Much attention has been given to understanding the causes of pre-ignition and identify lubricant or fuel properties and engine design and calibration considerations that impact its frequency. This helps to shift the pre-ignition limit to higher specific loads and allow further downsizing but does not fundamentally eliminate the problem. Real-time detection and mitigation of pre-ignition would thus be desirable to allow safe engine operation in pre-ignition-prone conditions. This study focuses on advancing the time of detection of pre-ignition in an engine cycle where it occurs. Furthermore, phase space transforms through time-delay embedding of cylinder pressure and principal component analysis were applied to same-cycle detection of pre-ignition and shown to enable detection on the order of a crank degree earlier than deviation in cylinder pressure can be identified through direct statistical observation of the pressure data. Additionally, it appears that the deviation of the trajectory in phase space may offer the opportunity to extend this method to further extend the detection window and allow more time for mitigation actions to occur.

42 ENGINEERING↗

SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 3 (Raw Temperature Replica Exchange Molecular Dynamics Trajectories)

Raw temperature replica-exchange molecular dyanmics trajectories of 23 different SARS CoV-2 Systems including S Protein ACE2-receptor binding-domain, MPro, PLPro, NSP3 ADRP (X-domain/macrodomain/phosphatase), NSP15 (endoribonuclease), NSP9, NSP10, NSP16, and N-protein N-terminus. Systems were prepared with charmm-gui and simulated using the GROMACS simulation software suite. Non-demuxed, i.e. discontinuous/constant temperature window, trajectories are provided in compressed GROMACS xtc and full trr formats. This data supplements the data releases DOI: 10.13139/OLCF/160650 'SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 1 (Simulation Initial Conditions and Parameters)' and DOI: 10.13139/OLCF/1657844 'SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 2 (Extracted Protein Coordinate Trajectories)'

59 BASIC BIOLOGICAL SCIENCES↗

Deriving spatially explicit direct and indirect interaction networks from animal movement data

Abstract Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems [“GPS”]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous‐time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions—individuals occurring at the same location, but at different times—while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease‐relevant interaction networks for two behaviorally differentiated species, wild pigs ( Sus scrofa ) that can host African Swine Fever and mule deer ( Odocoileus hemionus ) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30‐min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM‐Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine‐scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer–resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

59 BASIC BIOLOGICAL SCIENCES↗

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗

Fast 2D Bicephalous Convolutional Autoencoder for Compressing 3D Time Projection Chamber Data

High-energy large-scale particle colliders produce data at high speed in the order of 1 terabytes per second in nuclear physics and petabytes per second in high energy physics. Developing real-time data compression algorithms to reduce such data at high throughput to fit permanent storage has drawn increasing attention. Specifically, at the newly constructed sPHENIX experiment at the Relativistic Heavy Ion Collider (RHIC), a time projection chamber is used as the main tracking detector, which records particle trajectories in a volume of three-dimensional (3D) cylinder. The resulting data are usually very sparse with occupancy around 10.8%. Such sparsity presents a challenge to conventional learning-free lossy compression algorithms, such as SZ, ZFP, and MGARD. The 3D convolutional neural network (CNN)-based approach, Bicephalous Convolutional Autoencoder (BCAE), outperforms traditional methods both in compression rate and reconstruction accuracy. BCAE can also utilize the computation power of graphical processing units suitable for deployment in a modern heterogeneous highperformance computing environment. This work introduces two BCAE variants: BCAE++ and BCAE-2D. BCAE++ achieves a 15% better compression ratio and a 77% better reconstruction accuracy measured in mean absolute error compared with BCAE. BCAE-2D treats the radial direction as the channel dimension of an image, resulting in a 3× speedup in compression throughput. In addition, we demonstrate an unbalanced autoencoder with a larger decoder can improve reconstruction accuracy without significantly sacrificing throughput. Lastly, we observe both the BCAE++ and BCAE-2D can benefit more from using half-precision mode in throughput (76 - 79% increase) without loss in reconstruction accuracy. The source code and links to data and pretrained models can be found at https://github.com/BNL-DAQ-LDRD/NeuralCompression_v2

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Differentially Private Adaptive Noise Injection (DP-ANI) v1.0

Location data is collected from users continuously to understand their mobility patterns. Releasing the user trajectories may compromise user privacy. Therefore, the general practice is to release aggregated location datasets. However, private information may still be inferred from an aggregated version of location trajectories. Differential privacy (DP) protects the query output against inference attacks regardless of background knowledge. This software implements a differential privacy-based privacy model that protects the user's origins and destinations from being inferred from aggregated mobility datasets. This is achieved by injecting Planar Laplace noise to the user origin and destination GPS points. The noisy GPS points are then transformed into a link representation using a link-matching algorithm. Finally, the link trajectories form an aggregated mobility network. The injected noise level is selected using the Sparse Vector Mechanism. This DP selection mechanism considers the link density of the location and the functional category of the localized links. Compared to the different baseline models, including a k-anonymity method, our differential privacy-based aggregation model offers query responses that are close to the raw data in terms of aggregate statistics at both the network and trajectory-levels with maximum 9% deviation from the baseline in terms of network length.

Peisert, Sean [Lawrence Berkeley National Laborato↗

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗

Sea Surface Salinity Reemergence in an Updated North Atlantic In Situ Salinity Dataset

Monthly sea surface salinity (SSS) fields are constructed from observations, using objective mapping on a 1° x 1° grid in the Atlantic between 30°S and 50°N in the 1970-2016 period in an update of the dataset of Reverdin et al. Data coverage is heterogeneous, with increased density in 2002 when Argo floats become available, high density along Voluntary Observing Ship lines, and low density south of 10°S. Using lag correlation, the seasonal reemergence of SSS anomalies is investigated between 20° and 50°N in 5° x 5° boxes during the 1993-2016 period, both locally and remotely following the displacements of the deep mixed-layer waters estimated from virtual float trajectories derived from the daily AVISO surface geostrophic currents. Although SSS data are noisy, local SSS reemergence is detected in about half of the boxes, notably in the northeast and southeast, while little reemergence is seen in the central and part of the eastern subtropical gyre. In the same period, sea surface temperature (SST) reemergence is found only slightly more frequently, reflecting the short data duration. However, taking geostrophic advection into account degrades the detection of remote SSS and even SST reemergence. When anomalies are averaged over broader areas, robust evidence of a second and third SSS reemergence peak is found in the northeastern and southeastern parts of the domain, indicating long cold-season persistence of large-scale SSS anomalies, while only a first SST reemergence is seen. Finally, an oceanic reanalysis is used to confirm that the correlation analysis indeed reflects the reemergence of subsurface salinity anomalies.

54 ENVIRONMENTAL SCIENCES↗

Random Forest Regressor-Based Approach for Detecting Fault Location and Duration in Power Systems

Power system failures or outages due to short-circuits or “faults” can result in long service interruptions leading to significant socio-economic consequences. It is critical for electrical utilities to quickly ascertain fault characteristics, including location, type, and duration, to reduce the service time of an outage. Existing fault detection mechanisms (relays and digital fault recorders) are slow to communicate the fault characteristics upstream to the substations and control centers for action to be taken quickly. Fortunately, due to availability of high-resolution phasor measurement units (PMUs), more event-driven solutions can be captured in real time. In this paper, we propose a data-driven approach for determining fault characteristics using samples of fault trajectories. A random forest regressor (RFR)-based model is used to detect real-time fault location and its duration simultaneously. This model is based on combining multiple uncorrelated trees with state-of-the-art boosting and aggregating techniques in order to obtain robust generalizations and greater accuracy without overfitting or underfitting. Four cases were studied to evaluate the performance of RFR: 1. Detecting fault location (case 1), 2. Predicting fault duration (case 2), 3. Handling missing data (case 3), and 4. Identifying fault location and length in a real-time streaming environment (case 4). A comparative analysis was conducted between the RFR algorithm and state-of-the-art models, including deep neural network, Hoeffding tree, neural network, support vector machine, decision tree, naive Bayesian, and K-nearest neighborhood. Experiments revealed that RFR consistently outperformed the other models in detection accuracy, prediction error, and processing time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗