Search NASA⌕ Search

SEARCH · Search NASA

Results for “Trajectory Reconstruction”

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 91 records · Page 5

Real-Time event reconstruction for Nuclear Physics Experiments using Artificial Intelligence

Charged track reconstruction is a critical task in nuclear physics experiments, enabling the identification and analysis of particles produced in high-energy collisions. Machine learning (ML) has emerged as a powerful tool for this purpose, addressing the challenges posed by complex detector geometries, high event multiplicities, and noisy data. Traditional methods rely on pattern recognition algorithms like the Kalman filter, but ML techniques, such as neural networks, graph neural networks (GNNs), and recurrent neural networks (RNNs), offer improved accuracy and scalability. By learning from simulated and real detector data, ML models can identify and classify tracks, predict trajectories, and handle ambiguities caused by overlapping or missing hits. Moreover, ML-based approaches can process data in near-real-time, enhancing the efficiency of experiments at large-scale facilities like the Large Hadron Collider (LHC) and Jefferson Lab (JLAB). As detector technologies and computational resources evolve, ML-driven charged track reconstruction continues to push the boundaries of precision and discovery in nuclear physics. In these proceedings, we highlight advancements in charged track identification leveraging Artificial Intelligence within the CLAS12 detector, achieving a notable enhancement in experimental statistics compared to traditional methods. Additionally, we showcase real-time event reconstruction capabilities, including the inference of charged particle properties, such as momentum, direction, and species identification, at speeds matching data acquisition rates. These innovations enable the extraction of physics observables directly from the experiment in real-time.

Gavalian, Gagik (ORCID:0000000267385457)↗

Scalable In Situ Computation of Lagrangian Representations via Local Flow Maps

In situ computation of Lagrangian flow maps to enable post hoc time-varying vector field analysis has recently become an active area of research. However, the current literature is largely limited to theoretical settings and lacks a solution to address scalability of the technique in distributed memory. To improve scalability, we propose and evaluate the benefits and limitations of a simple, yet novel, performance optimization. Our proposed optimization is a communication-free model resulting in local Lagrangian flow maps, requiring no message passing or synchronization between processes, intrinsically improving scalability, and thereby reducing overall execution time and alleviating the encumbrance placed on simulation codes from communication overheads. To evaluate our approach, we computed Lagrangian flow maps for four time-varying simulation vector fields and investigated how execution time and reconstruction accuracy are impacted by the number of GPUs per compute node, the total number of compute nodes, particles per rank, and storage intervals. Our study consisted of experiments computing Lagrangian flow maps with up to 67M particle trajectories over 500 cycles and used as many as 2048 GPUs across 512 compute nodes. In all, our study contributes an evaluation of a communication-free model as well as a scalability study of computing distributed Lagrangian flow maps at scale using in situ infrastructure on a modern supercomputer.

Sane, Sudhanshu↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Simultaneous 10 kHz three-dimensional CH 2 O and tomographic PIV measurements in a lifted partially-premixed jet flame

High-speed, three-dimensional (3D) scalar-velocity field measurements are demonstrated in a lifted partially-premixed dimethyl-ether/air jet flame using simultaneous laser-induced fluorescence (LIF) of formaldehyde (CH 2 O) and tomographic particle image velocimetry (TPIV). The 3D LIF measurements are performed by raster scanning the laser beam from a 100 kHz pulse-burst laser across the probe volume using an acousto-optic deflector. The volumetric reconstruction of the LIF signal from ten parallel planes provides quasi-instantaneous 3D LIF measurements that are synchronized with 10 kHz TPIV measurements. Here, the temporally resolved CH 2 O-LIF and velocity field data are used to analyze Lagrangian particle trajectories and displacement speeds at the base of the lifted flame. The particle trajectories reveal flow structures that are difficult to observe in an Eulerian reference frame. Positive and negative displacement speeds are observed at the CH 2 O-LIF surfaces at the inner and outer regions of the jet flame with a maximum displacement speed of approximately eight times the laminar flame speed of a stoichiometric dimethyl-ether/air mixture. The observations in present study support formation of lean premixed flames near the flame base whose role in flame stabilization should be investigated in more detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optics reconstruction in the SCL section of the FNAL Linac

The Side-Coupled Linac (SCL) section of the FNAL linac accelerates the beam from 116 MeV to 401.5 MeV, operating at 22-24 mA beam current. Transverse focusing is performed by 32 quadrupoles, and the beam orbit is guided by 19 dipole correctors and measured by 29 BPMs. The bunch length is measured in a single location by a Bunch Shape Monitor (BSM). This paper presents a three-step reconstruction of the machine optics. First, the transverse and longitudinal Twiss parameters at the start of the SCL section are determined using quadrupole scans and BSM measurements at different settings of an upstream cavity. Second, the quadrupole calibrations are adjusted based on differential-trajectory measurements. Finally, the beam is propagated along the SCL linac using the code TraceWin. A comparison between TraceWin simulations and the beam envelope measured by the 12 wire scanners of the SCL linac was performed. Transverse and longitudinal beam parameters at the entrance and exit of the SCL section will be reported.

Carneiro, Jean-Paul [Fermilab]↗

Optics Reconstruction in the SCL Section of the FNAL Linac

The Side-CoupledLinac (SCL) section oftheFNAL linac acceleratesthe beamfrom 117 MeV to 401.5 MeV, operatingat 22-24mA beam current. Transverse focusing is performedby32 quadrupoles, and the beam orbit is guidedby19dipolecorrectors andmeasuredby29 BPMs. The bunch length is measured in a single location by a Bunch Shape Monitor (BSM). This paper presents a three-step reconstruction of the machine optics. First, the transverse and longitudinal Twiss parameters at the start of the SCL section are determined using quadrupole scans and BSM measurements at different settings o fan upstream cavity. Second, the quadrupole calibrations are adjusted based on differential-trajectory measurements. Finally, the beam is propagated along the SCL linac using the code Trace Win. A comparison between Trace Win simulations and the beam envelope measured bythe12wirescanners of the SCL linac was performed. Transverse and longitudinal beam parameters at the entrance of the transition section will be reported

Carneiro, Jean-Paul [Fermilab]↗

Analysis of tow architecture variability in biaxially braided composite tubes

Spatial deviations in composite braid unit cell dimensions limit mechanical performance and are challenging to model. Full-field measurements of the irregular tow placement in biaxially braided composite tubes were obtained by digital image correlation. Analysis of the reconstructed tow paths showed long-range defects along the hoop direction with short-range defects dominating the biased tow direction. A route toward numerical simulation of the composite response is established by developing a 3D model of the composite using the measured tow trajectories. Finally, a comparison to previously-published triaxial braid results indicates axial tows work to constrain unit cells dimensions along the hoop direction.

3-Dimensional reinforcement↗

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving ICARUS Track Reconstruction Algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino (SBN) program at Fermilab. The main goal of the experiment is to investigate the possibility of sterile neutrinos in the O(1 eV) mass region and provide clarification of the anomaly detected from the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber (LAr-TPC), that can provide excellent 3D imaging and calorimetric reconstruction of any ionizing particles. This detection technique allows a detailed study of neutrino interactions, spanning a wide energy spectrum (from a few keV to several hundreds of GeV). The detector consists of two identical adjacent modules, filled with a total of 760 tons of ultra-pure liquid argon. Each module houses two LAr-TPCs separated by a common cathode with a maximum drift distance of 1.5 m, equivalent to about 1 ms drift time for the nominal $500$ V/m electric drift field. The anode is made of three parallel wire planes positioned 3 mm apart, where the stainless-steel wires are oriented on each plane at a different angle with respect to the horizontal direction ($+60^\degree$,$-60^\degree$,$0^\degree$). The first two planes (Induction 1 and Induction 2) provide a non-destructive charge measurement, whereas the ionization charge is fully collected by the last collection plane. In total, 53248 wires with a 3 mm pitch and length up to 9 m are installed in the detector. In the first stage of the reconstruction, segments of waveforms corresponding to physical signals (hits) are searched for in the deconvolved wire waveform with a threshold-based hit-finding algorithm. Each hit is then fitted with a Gaussian, whose area is proportional to the number of drift electrons generating the signal. In the second stage of the reconstruction, hits are passed as input to Pandora, a framework software composed of different pattern recognition algorithms, that performs a 3D reconstruction of the full image recorded in the collected event, including the identification of interaction vertices and tracks and showers inside the TPC. These are organized into a hierarchical structure (called slice) of particles generated starting from a primary interaction vertex. In some cases, related to the inefficiencies in the hit detection or excessive deflection of the particle trajectory, Pandora breaks the particle's track into two or more smaller pieces and considers each piece as an independent track. We studied this phenomenon focusing on primary muons from ν_μ CC interactions contained in a single module with a track at least 20 cm long, to exclude delta rays. The study determined that about $7-8\%$ of the muon tracks are broken. Approximately $80\%$ of the times, Pandora assigns all segments of the track to the same slice (intra-slice track split), while in the remaining $20\%$ of the cases, one of the segments is associated with another slice (extra-slice track split). To mitigate this phenomenon, we designed an algorithm that detects and stitches the tracks broken by Pandora for the intra-slice split. In Monte Carlo simulations, the algorithm showed an efficiency exceeding $80\%$ and a purity exceeding $93\%$.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗

Topological background discrimination in the PandaX-III neutrinoless double beta decay experiment

The PandaX-III experiment plans to search for neutrinoless double beta decay (0ν ββ ) of 136 Xe in the China JinPing underground Laboratory (CJPL). The experiment will use a high pressure gaseous time projection chamber (TPC) to register both the energy and the electron track topology of an event. This article is devoted to demonstrate our particular detector setup capabilities for the identification of 0ν ββ and the consequent background reduction. As software tool we use REST, a framework developed for the reconstruction and simulation of TPC-based detector systems. We study the potential for background reduction by introducing appropriate parameters based on the properties of 0ν ββ events. We exploit for the first time not only the energy density of the electron track-ends, but also the electron scattering angles produced by an electron near the end of its trajectory. To implement this, we have added new algorithms for detector signal and track processing inside REST. Their assessment shows that background can be reduced by about 7 orders of magnitude while keeping 0ν ββ efficiency above 20% for the PandaX-III baseline readout scheme, a two-dimensional 3 mm pitch stripped readout. More generally, we use the potential of REST to handle 2D/3D data to assess the impact on signal-to-background significance at different detector granularities, and to validate the PandaX-III baseline choice. Finally, we demonstrate the additional potential to discriminate surface background events generated at the readout plane in the absence of to, by making use of event parameters related with the diffusion of electrons.

Majorana↗

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↗

Model and Standard Operating Procedures Supporting Signal Variation Flow Graph Analysis

This document will describe the principles of the Signal Variation Flow Technique, and the uncertainty models generated using it. We focus on the capture of the variation between the ideal signal and the measured signal. The ideal signal is defined to represent the signal output of a system whose full state behavior is known, with no variation in environment or during operation, and whose photons trajectories are not modified by the object. A CT uncertainty analysis is the result of two main steps. First, the radiography regime extends from the source to the collected image, which is 2D in the case of standard CT. This maps all upstream uncertainties into the variation observed on the radiograph and captures all variation in the physical domain. Second, the reconstruction regime extends from the captured images to the reconstructed 3D image. This regime is purely in the mathematical domain and corrects reconstruction algorithm artifacts/anomalies. The present work focuses on building the model through the radiography regime. The reconstruction regime is expected to be largely a study of algorithmic sensitivity, requiring the definition of a range of standardized tests through which the algorithms would be run. Radiographic variation would then be mapped through reconstruction sensitivities to predict the signal variation in the final image. An incomplete list for the reconstructed image variation output basis includes voxel density, edge blur and length variation, in analogous fashion to the basis functions presented in this work. The signal variation occurs in several forms at the radiograph. These forms are gathered into a complete basis set of functions describing all variation on the radiograph. The scale of each basis function is calculated independently via a specific SVFG. This includes 0D pixel noise (0DI), 0D energy noise (0DE), 1D blur (1DB), 1D length (1DL) or 2D position (2DP). These models are orthogonal in that they each explore a space in the signal variation domain that cannot be reached by the other basis functions. The variation basis functions, and associated SVFG models are split by output dimensionality, a term loosely used for categorization, and explained further below.

42 ENGINEERING↗

Time-resolved atomic-resolution Brownian tomography of single nanocrystals reveals size-dependent dynamics

Atomic-resolution structure identification of nanocrystals by graphene liquid cell electron microscopy (GLC-EM) has revealed that small, solubilized platinum nanocrystals consist of an ordered crystalline core surrounded by mobile surface atoms, which dissociate during oxidative etching, resulting in distinct temporal structural states. Requirements imposed by the 3D reconstruction algorithm limit the number of structural states that can be resolved. We introduce a regularized 3D reconstruction algorithm that exploits the redundancy inherent in the experimental data, allowing us to improve the time resolution. Our developments provide a comprehensive molecular picture at unprecedented spatial and temporal resolution of the nonlinear, linear, and fluctuating dynamic phenomena that single nanocrystals undergo during the GLC-EM experiment. We determined atomic structures of 66 temporal structural states, extracted from 15 time trajectories of individual nanocrystals. Large (478 to 698 atoms) and small (<300 atoms) nanocrystals show etching that preserves a stable core, whereas mid-sized (351 to 571 atoms) nanocrystals present dynamics that change the coordination of the core.

Meana-Pañeda, Rubén↗

The influence of recent and future climate change on spring Arctic cyclones

In recent decades, the Arctic has experienced rapid atmospheric warming and sea ice loss, with an ice-free Arctic projected by the end of this century. Cyclones are synoptic weather events that transport heat and moisture into the Arctic, and have complex impacts on sea ice, and the local and global climate. However, the effect of a changing climate on Arctic cyclone behavior remains poorly understood. This study uses high resolution (4 km), regional modeling techniques and downscaled global climate reconstructions and projections to examine how recent and future climatic changes alter cyclone behavior. Results suggest that recent climate change has not yet had an appreciable effect on Arctic cyclone characteristics. However, future sea ice loss and increasing surface temperatures drive large increases in the near-surface temperature gradient, sensible and latent heat fluxes, and convection during cyclones. The future climate can alter cyclone trajectories and increase and prolong intensity with greatly augmented wind speeds, temperatures, and precipitation. Such changes in cyclone characteristics could exacerbate sea ice loss and Arctic warming through positive feedbacks. The increasing extreme nature of these weather events has implications for local ecosystems, communities, and socio-economic activities.

54 ENVIRONMENTAL SCIENCES↗

Modeling root system growth around obstacles

State-of-the-Art models of Root System Architecture (RSA) do not allow simulating root growth around rigid obstacles. Yet, the presence of obstacles can be highly disruptive to the root system. We grew wheat seedlings in sealed petri dishes without obstacle and in custom 3D-printed rhizoboxes containing obstacles. Time-lapse photography was used to reconstruct the wheat root morphology network. We used the reconstructed wheat root network without obstacle to calibrate an RSA model implemented in the R-SWMS software. The root network with obstacles allowed calibrating the parameters of a new function that models the influence of rigid obstacles on wheat root growth. Experimental results show that the presence of a rigid obstacle does not affect the growth rate of the wheat root axes, but that it does influence the root trajectory after the main axis has passed the obstacle. The growth recovery time, i.e. the time for the main root axis to recover its geotropism-driven growth, is proportional to the time during which the main axis grows along the obstacle. Qualitative and quantitative comparisons between experimental and numerical results show that the proposed model successfully simulates wheat RSA growth around obstacles. Our results suggest that wheat roots follow patterns that could inspire the design of adaptive engineering flow networks.

59 BASIC BIOLOGICAL SCIENCES↗

Search in diphoton and dielectron final states for displaced production of Higgs or Z bosons with the ATLAS detector in s = 13 TeV p p collisions

A search is presented for displaced production of Higgs bosons or Z bosons, originating from the decay of a neutral long-lived particle (LLP) and reconstructed in the decay modes H → γ γ and Z → ee . The analysis uses the full Run 2 dataset of proton-proton collisions delivered by the LHC at an energy of $\sqrt{s}$ = 13 TeV between 2015 and 2018 and recorded by the ATLAS detector, corresponding to an integrated luminosity of 139 fb -1 . Exploiting the capabilities of the ATLAS liquid argon calorimeter to precisely measure the arrival times and trajectories of electromagnetic objects, the analysis searches for the signature of pairs of photons or electrons which arise from a common displaced vertex and which arrive after some delay at the calorimeter. The results are interpreted in a gauge-mediated supersymmetry breaking model with pair-produced Higgsinos that decay to LLPs, and each LLP subsequently decays into either a Higgs boson or a Z boson. The final state includes at least two particles that escape direct detection, giving rise to missing transverse momentum. No significant excess is observed above the background expectation. The results are used to set upper limits on the cross section for Higgsino pair production, up to a $\bar{χ}^0_1$ mass of 369 (704) GeV for decays with 100% branching ratio of $\bar{χ}^0_1$ to Higgs (Z) bosons for a $\bar{χ}^0_1$ lifetime of 2 ns. A model-independent limit is also set on the production of pairs of photons or electrons with a significant delay in arrival at the calorimeter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Diagnostic-free onboard battery health assessment

Diverse usage patterns induce complex and variable aging behaviors in lithiumion batteries, complicating accurate health diagnosis and prognosis. Separate diagnostic cycles are often used to untangle the battery’s current state of health from prior complex aging patterns. However, these same diagnostic cycles alter the battery’s degradation trajectory, are time-intensive, and cannot be practically performed in onboard applications. Here, in this work, we leverage portions of operational measurements in combination with an interpretable machine learning model to enable rapid, onboard battery health diagnostics and prognostics without offline diagnostic testing and the requirement of historical data. We integrate mechanistic constraints within an encoder-decoder architecture to extract electrode states in a physically interpretable latent space and enable improved reconstruction of the degradation path. The health diagnosis model framework can be flexibly applied across diverse application interests with slight fine-tuning.

battery aging reconstruction↗

Visualization of fast ion phase-space flow in plasmas well-below, near and well-above Alfvén eigenmode stability threshold in tokamak

Abstract Transport of fast ions along certain local phase space paths, referred to as fast ion phase-space flow, has been systematically measured by an imaging neutral particle analyzer (INPA) in three plasma regimes, which are well-below, near, and well-above the Alfvén eigenmode stability threshold (Du et al 2021 Phys. Rev. Lett. 127 235002). (1) In plasmas well-below the Alfvén eignenmodes (AE) stability threshold, fast ions are well-confined on passing particle orbits without noticeable transport over the phase space. The observed INPA images agree well with the synthetic INPA images, using the fast ion distribution predicted by neoclassical theory. (2) In plasmas near the AE stability threshold, INPA images in the presence of AE activity moderately deviate from those without AE activity. The image difference can be well interpreted by AE-driven, phase-space fast ion flow. Paths of this flow over the velocity space (or streamlines) are reconstructed by the intersection lines of curved E ʹ and µ surfaces, referred to as E ʹ and µ line (where E ′ ≡ E − ( ω / n ) P ζ ; E , P ζ and µ are the energy, canonical toroidal momentum and magnetic moment of ions; ω and n are the angular frequency and toroidal mode number of AEs, respectively). Resonant fast ions move radially inward by gaining energy and move radially outward by losing energy and the trajectory well aligns with the E ʹ and µ lines that pass through the mode resonances near the injection energy of neutral beams. (3) In plasmas well-above the AE stability threshold, fast ion phase-space dynamics shows additional features. Fast ions are transported out of the birth positions so promptly along the streamlines that the slowing-down process from the injection energy is not observable, exhibiting strong critical gradient behavior at local phase space. As a result, the increase of electron temperature is very small, in spite of an increase of beam power by ∼ 45 % . It should be emphasized that the directions of phase-space transport, induced by AEs with different frequencies, structures and mode numbers, do not largely differ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗