Search NASA⌕ Search

SEARCH · Search NASA

Results for “Trajectory Modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Temporal covariation of island arc Sr isotopes and seawater chemistry over the past 2 billion years

The chemical compositions of island arc basalts (IAB) reflect contributions from the mantle as well as fluids and melts from the subducting slab. Addition of radiogenic seawater Sr to oceanic crust through hydrothermal alteration and subsequent subduction is often invoked to explain elevated 87 Sr/ 86 Sr signatures in modern IAB. However, changes in the 87 Sr/ 86 Sr of island arc magmatic rocks through time has not been investigated, limiting our understanding of the factors influencing the Sr budgets of arcs throughout Earth’s history. To address this, we compiled 87 Sr/ 86 Sr values from island arc magmatic rocks ranging in age from modern to Paleoproterozoic, only including data from island arc localities that best preserve initial magmatic 87 Sr/ 86 Sr. Median initial 87 Sr/ 86 Sr values are consistently elevated compared to depleted mantle 87 Sr/ 86 Sr over this period, indicating persistent enrichment in radiogenic Sr in island arcs. Moreover, the elevation in island arc 87 Sr/ 86 Sr relative to the depleted mantle is variable. A notable rise in island arc 87 Sr/ 86 Sr during the late Neoproterozoic coincides with a steep increase in seawater 87 Sr/ 86 Sr and Sr concentration. To investigate this potential connectivity, we modeled the 87 Sr/ 86 Sr of island arc magmas between 0 and 830 Ma with inputs of depleted mantle 87 Sr/ 86 Sr, seawater 87 Sr/ 86 Sr, and seawater Sr concentration. The model reproduces the overall trajectory of the compiled data. We interpret the observed temporal variation in island arc 87 Sr/ 86 Sr values and its close association with fluctuations in seawater chemistry as evidence that changes in marine geochemistry have strongly influenced the Sr isotopic record of island arc magmas over time.

Science & Technology - Other Topics↗
DOE OSTI2024

Command of three-dimensional solitary waves via photopatterning

Multidimensional solitons are prevalent in numerous research fields. In orientationally ordered soft matter system, three-dimensional director solitons exemplify the localized distortion of molecular orientation. However, their precise manipulation remains challenging due to unpredictable and uncontrolled generation. Here, we utilize preimposed programmable photopatterning in nematics to control the kinetics of director solitons. This enables both unidirectional and bidirectional generation at specific locations and times, confinement within micron-scaled patterns of diverse shapes, and directed propagation along predefined trajectories. A focused dynamical model provides insight into the origins of these solitons and aligns closely with experimental observations, underscoring the pivotal role of anchoring conditions in soliton manipulation. Our findings pave the way for diverse fundamental research avenues and promising applications, including microcargo transportation and optical information processing.

Science & Technology - Other Topics↗
DOE OSTI2025

A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations (2006–2011)

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense.

astrometry↗
DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2025

Back trajectories for TBS flights (c1-level)

The ARMTRAJ VAP provides four trajectory datasets initialized at ARM deployment coordinates and configured using ARM datasets. The four trajectory datasets support aerosol, cloud, and planetary boundary layer research. Trajectory calculations use the HYSPLIT model informed by the ERA5 reanalysis dataset at its highest spatial resolution (~31 km). For each sample in each of the four datasets, HYSPLIT will also be run at multiple starting locations surrounding ARM deployments, enabling an ensemble of runs from which the mean and variability (estimated uncertainty) of each sample's trajectory coordinates, thermodynamic properties, or other fields will be reported.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗
DOE OSTI2024

The ab initio non-crystalline structure database: empowering machine learning to decode diffusivity

Non-crystalline materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of non-crystalline materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the largest computed non-crystalline structure database to date, generated from systematic and accurate ab initio molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductivity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching non-crystalline structures provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in universal machine learning potentials, impacting design beyond that of non-crystalline materials. In addition, combining diffusion trajectories from our dataset with models that predict liquidus viscosity and melting temperature could be utilized to develop models for predicting glass-forming ability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗
DOE OSTI2024

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Huge ensembles – Part 2: Properties of a huge ensemble of hindcasts generated with spherical Fourier neural operators

Abstract. In Part 1, we created an ensemble based on spherical Fourier neural operators. As initial condition perturbations, we used bred vectors, and as model perturbations, we used multiple checkpoints trained independently from scratch. Based on diagnostics that assess the ensemble's physical fidelity, our ensemble has comparable performance to operational weather forecasting systems. However, it requires orders-of-magnitude fewer computational resources. Here in Part 2, we generate a huge ensemble (HENS), with 7424 members initialized each day of summer 2023. We enumerate the technical requirements for running huge ensembles at this scale. HENS precisely samples the tails of the forecast distribution and presents a detailed sampling of internal variability. HENS has two primary applications: (1) as a large dataset with which to study the statistics and drivers of extreme weather and (2) as a weather forecasting system. For extreme climate statistics, HENS samples events 4σ away from the ensemble mean. At each grid cell, HENS increases the skill of the most accurate ensemble member and enhances coverage of possible future trajectories. As a weather forecasting model, HENS issues extreme weather forecasts with better uncertainty quantification. It also reduces the probability of outlier events, in which the verification value lies outside the ensemble forecast distribution.

Mahesh, Ankur↗
DOE OSTI2024

Extracting Vehicle Trajectories from Partially Overlapping Roadside Radar

This work presents a methodology for extracting vehicle trajectories from six partially-overlapping roadside radars through a signalized corridor. The methodology incorporates radar calibration, transformation to the Frenet space, Kalman filtering, short-term prediction, lane-classification, trajectory association, and a covariance intersection-based approach to track fusion. The resulting dataset contains 79,000 fused radar trajectories over a 26-h period, capturing diverse driving scenarios including signalized intersections, merging behavior, and a wide range of speeds. Compared to popular trajectory datasets such as NGSIM and highD, this dataset offers extended temporal coverage, a large number of vehicles, and varied driving conditions. The filtered leader–follower pairs from the dataset provide a substantial number of trajectories suitable for car-following model calibration. The framework and dataset presented in this work has the potential to be leveraged broadly in the study of advanced traffic management systems, autonomous vehicle decision-making, and traffic research.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Summertime Carbonaceous Aerosol in Interior Versus Coastal Northern Alaska

Abstract Rapid warming is likely increasing primary production and wildfire occurrence in the Arctic. Projected changes in carbonaceous aerosols during the summer will impact atmospheric chemistry and climate, but our understanding of these processes is limited by sparse observations. Here, we characterize carbonaceous aerosol in Alaska, USA: Toolik Field Station in the Interior and the Atmospheric Radiation Measurement facility at Utqiaġvik on the Arctic coast, during the summers of 2022 and 2023. We estimated PM 2.5 and PM 10 concentrations using laser light scattering (PurpleAir sensors) and examined total carbon (TC) and its organic carbon (OC) and elemental carbon (EC) fractions in total suspended particles (TSP). We investigated the dominant sources of carbonaceous aerosol using air mass backward‐trajectories from the NOAA HYSPLIT model and radiocarbon source apportionment of TC. TC concentrations were about twice as high in the Interior compared to the coast, with contemporary sources dominating at both Toolik (95%–99%) and Utqiaġvik (86%–89%) over minor contributions from fossil sources. Elevated PM, TC, OC, and EC concentrations coincided with major boreal forest fire activity in North America that brought smoke to the region. The radiocarbon signature of EC measured at Toolik during these wildfire events indicated that over 90% of the EC came from contemporary sources. Our measurements demonstrate the potential for Arctic aerosol concentrations to respond significantly to climate warming‐induced changes to the landscape and emphasize the need for continuous atmospheric monitoring to advance our understanding of this rapidly changing environment.

Welch, Allison M. [Department of Earth System Scie↗
DOE OSTI2024

Robust Size-Effects Compensation Through Regularized Lever-Arm Estimation

Size effects are an unavoidable nuisance in inertial navigation using sensors which are not co-located at the navigational point of interest. When estimating transforms between the navigation point and sensor locations, some trajectories preclude observation of all model parameters. Regularization is proposed to avoid over-fitting size-effects models. The result yields robust size effects compensation in other regions of flight.

42 ENGINEERING↗
DOE OSTI2025

Distinguishing Orbiting and Infalling Dark Matter Particles with Machine Learning

Dark matter halos are typically defined as spheres that enclose some overdensity, but these sharp, somewhat arbitrary boundaries introduce nonphysical artifacts such as backsplash halos, pseudo-volution, and an incomplete accounting of halo mass. A more physically motivated alternative is to define halos as the collection of particles that are physically orbiting within their potential well. However, existing methods to classify particles as orbiting or infalling suffer from trade-offs between accuracy, computational cost, and generalizability across cosmologies. We present an efficient, yet accurate, supervised machine learning approach using decision trees. The classification is based on only the particle radii and velocities at two epochs. Compared to detailed analysis of particle trajectories, we find that our model matches the classification of 97% of particles. Consequently, we are able to quickly and accurately reproduce the density profiles of the orbiting and infalling components out to many virial radii. We demonstrate that our model generalizes to a significantly different cosmology that lies outside the training data set. We make publicly available both our final model and the code to train similar models.

79 ASTRONOMY AND ASTROPHYSICS↗