Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BALLISTIC TRAJECTORY”

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.

Machine Learning Analysis of Direct Dynamics Trajectory Outcomes for Thermal Deazetization of 2,3-Diazabicyclo[2.2.1]hept-2-ene

Experimentally, the thermal gas-phase deazetization of 2,3-diazabicyclo[2.2.1]hept-2-ene (1) results in the loss of N2 and the formation of bicyclo products 3 (exo) and 4 (endo) in a nonstatistical ratio, with preference for the exo product. Here, we report unrestricted M06-2X quasiclassical trajectories initialized from the concerted N2 ejection transition state that were able to replicate the experimental preference to form 3. We found that the 3:4 ratio results from the relative amounts of very fast (ballistic) exotype trajectories versus trajectories that lead to the 1,3-diradical intermediate 2. These quasiclassical trajectories provided a set of transition-state vibrational, velocity, momenta, and geometric features for the machine learning analysis. Additionally, a selection of popular supervised classification algorithms (e.g., random forest) provided poor prediction of trajectory outcomes based on only transition-state vibrational quanta and energy features. However, these machine learning models provided more accurate predictions using atomic velocities and atomic positions, attaining ~70% accuracy using initial conditions and between 85 and 95% accuracy at later reaction time steps. This increased accuracy allowed the feature importance analysis to reveal that, at the later-time analysis, the methylene bridge out-of-plane bending is correlated with trajectory outcomes for the formation of either the exo product or toward the diradical intermediate. Possible reasons for the struggle of machine learning algorithms to classify trajectories based on transition-state features is the heavily overlapping feature values, the finite but very large possible vibrational mode combinations, and the possibility of chaos as trajectories propagate. We examined this chaos by comparing a set of nearly identical trajectories that differed by only a very small scaling of the kinetic energies resulting from the transition-state reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gate-tunable Veselago interference in a bipolar graphene microcavity

The relativistic charge carriers in monolayer graphene can be manipulated in manners akin to conventional optics. Klein tunneling and Veselago lensing have been previously demonstrated in ballistic graphene pn-junction devices, but collimation and focusing efficiency remains relatively low, preventing realization of advanced quantum devices and controlled quantum interference. Here, we present a graphene microcavity defined by carefully-engineered local strain and electrostatic fields. Electrons are manipulated to form an interference path inside the cavity at zero magnetic field via consecutive Veselago refractions. The observation of unique Veselago interference peaks via transport measurement and their magnetic field dependence agrees with the theoretical expectation. We further utilize Veselago interference to demonstrate localization of uncollimated electrons and thus improvement in collimation efficiency. Our work sheds new light on relativistic single-particle physics and provide a new device concept toward next-generation quantum devices based on manipulation of ballistic electron trajectory.

36 MATERIALS SCIENCE↗

High-energy ballistic electrons in low-pressure radio-frequency plasmas

This study demonstrates the presence of a small number of high-energy ballistic electrons (HEBEs) that originate from secondary electrons in low-pressure radio-frequency (rf) plasmas. The kinetic behaviors of the HEBEs are illustrated through electron energy probability functions from the fully kinetic particle-in-cell simulations, showing two wavy high-energy tails and two bifurcations during one rf cycle. Test-particle simulations and a semi-analytical method associated with nonlocal electron kinetics are performed to characterize the HEBE trajectories, which reveal the ballistic nature of the HEBEs and their typical bouncing features between the rf sheaths. Parameter dependence of the HEBEs on the discharge conditions (e.g., gas pressure, gap distance, and rf frequency) are identified, which is relevant to the plasma collisionality. With a pronounced presence of HEBEs, the overall impacts of the secondary electron emission on discharge parameters, such as electron power absorption and ionization rate, are also illustrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vehicle Tracking

The detection and tracking of vehicles moving at hypersonic speeds (> Mach 5) at altitudes of ~100 km has been of growing interest due to the ongoing development of hypersonic glide vehicles (HGVs) by the U.S., Russia and China. This challenge is notably different from traditional ballistic missile systems whose re-entry vehicles move at hypersonic velocities but at much higher altitudes on a ballistic course. Hypersonic glide vehicles are also designed to be able to make sudden turns along their flight path making prediction of their trajectories difficult. Their ability to maneuver, combined with their high velocity and altitude make them difficult to track. The combination of these characteristics presents a challenge for anyone attempting to defend against a hypersonic glide vehicle attack by generally shortening the time available to react. That reaction can involve moving potential targets or attempting to detect the hypersonic glide vehicle and intercept it before it reaches the intended target. Figure 1 compares the trajectory of hypersonic vehicles compared to that of a traditional ballistic missile.

42 ENGINEERING↗

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗

Along-Trajectory Acoustic Signal Variations Observed During the Hypersonic Re-Entry of the OSIRIS-REx Sample Return Capsule

The re-entry of the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) sample return capsule (SRC) on 24 September 2023 presented a rare opportunity to study atmospheric entry dynamics through a dense network of ground-based infrasound sensors. As the first interplanetary capsule to re-enter over the United States since Stardust in 2006, this event allowed for unprecedented observations of infrasound signals generated during hypersonic descent. We deployed 39 single-sensor stations across Nevada and Utah, strategically distributed to capture signals from distinct trajectory points. Infrasound data were analyzed to examine how signal amplitude and period vary with altitude and propagation path for a nonablating hypersonic object with well-defined physical and aerodynamic properties. Raytracing simulations incorporated atmospheric specifications from the ground-2-space model to estimate source altitudes for observed signals. Results confirmed ballistic arrivals at all stations, with source altitudes ranging from 44 to 62 km along the trajectory. Signal period and amplitude exhibited strong dependence on source altitude, with higher altitudes corresponding to lower amplitudes, longer periods, and reduced high-frequency content. Regression analysis demonstrated strong correlations between signal characteristics and both altitude and propagation geometry. Our results suggest, when attenuation is considered, the amplitude is primarily determined by the source, with the propagation path playing a secondary role over the distances examined. These findings emphasize the utility of controlled SRC re-entries for advancing our understanding of natural meteoroid dynamics, refining atmospheric entry models, and improving methodologies for planetary defense. The OSIRIS-REx SRC campaign represents the most comprehensive infrasound study of a hypersonic re-entry to date, showcasing the potential of coordinated geophysical observational networks for high-energy atmospheric phenomena, including space debris re-entries.

58 GEOSCIENCES↗

Spin-Vibronic Effect in Photoinduced Electron Transfer

The spin-vibronic effect (SVE) accelerates quantum-mechanically forbidden electronic transitions, but experimental manifestations of this phenomenon remain limited. In this contribution, the role of SVE in photoinduced electron transfer (PET) dynamics is probed by coherent vibrational wavepacket (CVWP) motions in Pt­(II) dimer-naphthalene diimide donor–acceptor complexes. Metal–metal-to-ligand charge-transfer (MMLCT) excitation of the donor triggers ballistic PET, driving the molecule away from equilibrium. The SVE then directs the PET trajectory from the 1 MMLCT excited state to an intermediate ligand-centered triplet state, modulating the CVWP dynamics along the reaction coordinates. Herzberg–Teller type oscillations indicate that the Pt–Pt stretching vibration, arising from the formation of the excited triplet intermediate, serves as the reaction coordinate for charge separation. Finally, our experimental findings advance the understanding of the interplay between electronic, vibrational, and spin degrees of freedom in ultrafast photoactivated processes.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗

Multi‐Sensor Trajectory Reconstruction of the 24 April 2025 Alaska Fireball and Implications for Planetary Defense

On 24 April 2025 at 18:30:57 UTC, a bright daytime fireball over Southcentral Alaska was detected by 37 seismic stations, 16 single infrasound sensors, and four infrasound arrays, yielding 30 ballistic and multiple fragmentation arrivals. Here, the unprecedented density of seismoacoustic coverage enabled detailed reconstruction of the event using acoustic signals, with fragmentation source locations further guiding the identification of Doppler weather radar signatures of a meteorite fall. Incorporation of a radar-derived terminal point yielded a final trajectory solution, which agreed closely with an independent optical trajectory solution from video analysis. The reconstructed entry parameters from seismoacoustic analysis indicate a velocity of 25.3 km/s, an entry angle of 19°, and an energy release of ∼38 t TNT equivalent. Assuming a chondritic composition, the pre-entry object diameter was ∼0.7 m. Using orbital parameters from the optical solution, we estimate meteoroid composition as most likely an L-type ordinary chondrite. The event occurred in the sub-Arctic, where space-based optical systems face challenges in detection, demonstrating the critical role of dense ground-based seismoacoustic networks in characterizing high-latitude atmospheric entries. This uniquely well-recorded event demonstrates the capability of dense seismoacoustic networks to constrain bolide trajectories, energetics, and fragmentation, with radar and optical data providing critical confirmation and complementary perspectives. These results bridge the methodological gap between planetary-defense monitoring of natural impactors and space-traffic analyses of artificial reentries, illustrating how multi-sensor integration can deliver calibration-grade trajectories even for unpredicted events.

Fireball↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗

Imaging the flow of holes from a collimating contact in graphene

A beam of holes formed in graphene by a collimating contact is imaged using a liquid-He cooled scanning probe microscope (SPM). The mean free path of holes is greater than the device dimensions. A zigzag shaped pattern on both sides of the collimating contact absorbs holes that enter at large angles. The image charge beneath the SPM tip defects holes, and the pattern of flow is imaged by displaying the change in conductance between contacts on opposite sides, as the tip is raster scanned across the sample. Collimation is confirmed by bending hole trajectories away from the receiving contact with an applied magnetic field. The SPM images agree well with ray-tracing simulations.

36 MATERIALS SCIENCE↗

Drag Coefficient Analysis of Arbitrary Fragment Geometry with Emphasis on Meshing Sensitivity Studies.

This report characterizes Mach number effects in the hypersonic range from Mach 3 to Mach 7 for an arbitrary, computationally generated fragment. Force and moment coefficients in the x, y, z direction were successfully obtained and resulting trends were compared with theoretical expectations. Three orientations were tested in the form of unit vectors: [1,0,0], [0,1,0], and [0,0,1]; these represent the i, j, k principle axes. The final results ultimately showcased a trend very similar to that shown in previous literature — drag coefficients for a given orientation decreased with increasing Mach number by a very small amount in the hypersonic range. It was therefore concluded that the aerodynamic quantities used to obtain the fragment trajectory at low hypersonic Mach numbers (such as Mach 3) can be used to characterize the aerodynamic qualities at higher hypersonic Mach numbers (Mach 5, 7) with a reasonably small margin. A sensitivity study was also conducted which ultimately defined drag effects as a result of changing grid spacing and normal extrusion parameters for the unstructured tetrahedral mesh used in the simulation. Further work can be done in the form of experimental validation, specifically with regards to wind tunnels and ballistic range testing. Error can be reduced by testing more mesh variabilities and capturing a larger amount of fragment orientations. The results for force and moment characteristics dependent on meshing parameters and high Mach numbers are satisfactory and consistent with expected trends.

42 ENGINEERING↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.

dust↗