Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PARTICLE 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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Ice Particle Transport Analysis With Phase Change for the E(sup 3) Turbofan Engine Using LEWICE3D Version 3.2

Ice Particle trajectory calculations with phase change were made for the Energy Efficient Engine (E(sup 3)) using the LEWICE3D Version 3.2 software. The particle trajectory computations were performed using the new Glenn Ice Particle Phase Change Model which has been incorporated into the LEWICE3D Version 3.2 software. The E(sup 3) was developed by NASA and GE in the early 1980 s as a technology demonstrator and is representative of a modern high bypass turbofan engine. The E(sup 3) flow field was calculated using the NASA Glenn ADPAC turbomachinery flow solver. Computations were performed for the low pressure compressor of the E(sup 3) for a Mach 0.8 cruise condition at 11,887 m assuming a standard warm day for ice particle sizes of 5, 20, and 100 microns and a free stream particle concentration of 0.3 g/cu m. The impingement efficiency results showed that as particle size increased average impingement efficiencies and scoop factors increased for the various components. The particle analysis also showed that the amount of mass entering the inner core decreased with increased particle size because the larger particles were less able to negotiate the turn into the inner core due to particle inertia. The particle phase change analysis results showed that the larger particles warmed less as they were transported through the low pressure compressor. Only the smallest 5 micron particles were warmed enough to produce melting and the amount of melting was relatively small with a maximum average melting fraction of 0.836. The results also showed an appreciable amount of particle sublimation and evaporation for the 5 micron particles entering the engine core (22 percent).

Bidwell, Colin, S.↗
DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗
DOE OSTI2021

Charged Particle Tracking via Edge-Classifying Interaction Networks

Recent work has demonstrated that geometric deep learning methods such as graph neural networks (GNNs) are well suited to address a variety of reconstruction problems in high-energy particle physics. In particular, particle tracking data are naturally represented as a graph by identifying silicon tracker hits as nodes and particle trajectories as edges, given a set of hypothesized edges, edge-classifying GNNs identify those corresponding to real particle trajectories. In this work, we adapt the physics-motivated interaction network (IN) GNN toward the problem of particle tracking in pileup conditions similar to those expected at the high-luminosity Large Hadron Collider. Assuming idealized hit filtering at various particle momenta thresholds, we demonstrate the IN’s excellent edge-classification accuracy and tracking efficiency through a suite of measurements at each stage of GNN-based tracking: graph construction, edge classification, and track building. The proposed IN architecture is substantially smaller than previously studied GNN tracking architectures; this is particularly promising as a reduction in size is critical for enabling GNN-based tracking in constrained computing environments. Furthermore, the IN may be represented as either a set of explicit matrix operations or a message passing GNN. Efforts are underway to accelerate each representation via heterogeneous computing resources towards both high-level and low-latency triggering applications.

accelerator physics↗
NASA NTRS1982

Solid particle dynamic behavior through twisted blade rows

The particle trajectory calculations provide the essential information which is required for predicting the pattern and intensity of turbomachinery erosion. Consequently, the evaluation of the machine performance deterioration due to erosion is extremely sensitive to the accuracy of the flow field and blade geometry representation in the trajectory computational model. A model is presented that is simple and efficient yet versatile and general to be applicable to axial, radial and mixed flow machines, and to inlets, nozzles, return passages and separators. The results of the computations are presented for the particle trajectories through a row of twisted vanes in the inlet flow field. The effect of the particle size on their trajectories, blade impacts, and on their redistribution and separation are discussed.

Hamed, A.↗
DOE OSTI2022

Ultracold neutron storage simulation using the Kassiopeia software package

The Kassiopeia software package was originally developed to simulate electromagnetic fields and charged particle trajectories for neutrino mass measurement experiments. Recent additions to Kassiopeia also allow it to simulate neutral particle trajectories in magnetic fields based on their magnetic moments. Two different methods were implemented: an exact method that can work for arbitrary fields and an adiabatic method that is limited to slowly-varying fields but is much faster for large precession frequencies. Additional interactions to simulate reflection of ultracold neutrons (UCNs) from material walls and to allow spin–flip pulses were also added. These tools were used to simulate neutron precession in a room temperature neutron electric dipole moment experiment and predict the values of the longitudinal and transverse relaxation times as well as the trapping lifetime. All three parameters are found to closely match the experimentally determined values when simulated with both the exact and adiabatic methods, confirming that Kassiopeia is able to accurately simulate neutral particles. This opens the door for future uses of Kassiopeia to prototype the next generation of atomic traps and UCN experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Calorimetric classification of track-like signatures in liquid argon TPCs using MicroBooNE data

The MicroBooNE liquid argon time projection chamber located at Fermilab is a neutrino experiment dedicated to the study of short-baseline oscillations, the measurements of neutrino cross sections in liquid argon, and to the research and development of this novel detector technology. Accurate and precise measurements of calorimetry are essential to the event reconstruction and are achieved by leveraging the TPC to measure deposited energy per unit length along the particle trajectory, with mm resolution. We describe the non-uniform calorimetric reconstruction performance in the detector, showing dependence on the angle of the particle trajectory. Such non-uniform reconstruction directly affects the performance of the particle identification algorithms which infer particle type from calorimetric measurements. This work presents a new particle identification method which accounts for and effectively addresses such non-uniformity. The newly developed method shows improved performance compared to previous algorithms, illustrated by a 93.7% proton selection efficiency and a 10% muon mis-identification rate, with a fairly loose selection of tracks performed on beam data. The performance is further demonstrated by identifying exclusive final states in ν μ CC interactions. While developed using MicroBooNE data and simulation, this method is easily applicable to future LArTPC experiments, such as SBND, ICARUS, and DUNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Tracer Particles for Core-collapse Supernova Nucleosynthesis: The Advantages of Moving Backward

Abstract After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multidimensional simulations. However, obtaining nucleosynthesis yields based on such multidimensional core-collapse supernova simulations is not straightforward. Frequently, tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamic simulation. Reconstruction can be done in a number of ways, and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction to the point when nuclear statistical equilibrium was last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Three-dimensional water droplet trajectory code validation using an ECS inlet geometry

A task was completed under NASA contract, the purpose of which was to validate a three-dimensional particle trajectory code with existing test data obtained from the Icing Research Tunnel at NASA-LeRC. The geometry analyzed was a flush-mounted environmental control system (ECS) inlet. Results of the study indicated good overall agreement between analytical predictions and wind tunnel test results at most flight conditions. Difficulties were encountered when predicting impingement characteristics of the droplets less than or equal to 13.5 microns in diameter. This difficulty was corrected to some degree by modifications to a module of the particle trajectory code; however, additional modifications will be required to accurately predict impingement characteristics of smaller droplets.

Breer, Marlin D.↗
NASA NTRS1989

Trajectory-capture cell instrumentation for measurement of dust particle mass, velocity and trajectory, and particle capture

The development of the polyvinylidene fluoride (PVDF) dust detector for space missions--such as the Halley Comet Missions where the impact velocity was very high as well as for missions where the impact velocity is low was extended to include: (1) the capability for impact position determination - i.e., x,y coordinate of impact; and (2) the capability for particle velocity determination using two thin PVDF sensors spaced a given distance apart - i.e., by time-of-flight. These developments have led to space flight instrumentation for recovery-type missions, which will measure the masses (sizes), fluxes and trajectories of incoming dust particles and will capture the dust material in a form suitable for later Earth-based laboratory measurements. These laboratory measurements would determine the elemental, isotopic and mineralogical properties of the captured dust and relate these to possible sources of the dust material (i.e., comets, asteroids), using the trajectory information. The instrumentation described here has the unique advantages of providing both orbital characteristics and physical and chemical properties--as well as possible origin--of incoming dust.

Simpson, J. A.↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Two refinements to the quasi-linear theory of ion radial diffusion are proposed and examined analytically with simulations of particle trajectories. The resonance-broadening correction by Dungey (1965) is applied to the quasi-linear diffusion theory by Faelthammar (1965) for an individual model storm. Quasi-linear theory is then applied to the mean diffusion coefficients resulting from simulations of particle trajectories in 20 model storms. The correction for drift-resonance broadening results in quasi-linear diffusion coefficients with discrepancies from the corresponding simulated values that are reduced by a factor of about 3. Further reductions in the discrepancies are noted following the averaging of the quasi-linear diffusion coefficients, the simulated coefficients, and the resonance-broadened coefficients for the 20 storms. Quasi-linear theory provides good descriptions of particle transport for a single storm but performs even better in conjunction with the present ensemble-averaging.

Chen, Margaret W.↗
DOE OSTI2021

Improvement in the spatial resolution of heavy ion beam probe measurements through application of ion optics

A technique for more accurately modeling and improving the spatial resolution of heavy ion beam probe measurements is described. We use a set of particle trajectories to numerically determine the focusing properties of a complicated three-dimensional magnetic field and characterize those properties with a transfer matrix. We then modify the transfer matrix approach of traditional ion optics to include a parameter that describes the ionization location of the detected ions. The ion optics model calculated using this technique enables a more accurate description of the particle trajectories than previously feasible. The model also allows one to easily determine an initial beam focus that could be used during experimental operation to optimize the spatial resolution of measurements. The technique has been applied to the design of a heavy ion beam probe diagnostic for the Wendelstein 7-X stellarator and improvements in the modeled spatial resolution by a factor of about 2 over previous estimates are possible. As a result, the improved spatial resolution will enable measurements of plasma fluctuations with smaller wavelengths than would otherwise be possible.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional ice accretion modeling method is initiated. This first step toward creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flowfields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is intended as a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3D method for a MS-317 swept wing geometry are projected onto a 2D plane normal to the wing leading edge and compared to 2D results for the same geometry. These results indicate that the flowfield over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3D calculation.

Potapczuk, M. G.↗
NASA NTRS2006

Enhancing the ATIC Charge Resolution

The Advanced Thin Ionization Calorimeter (ATIC) experiment measures the energy spectra of elements, from H to Fe, in the energy region from about 100 GeV to tens of TeV. The ATIC instrument was flown twice in long-duration balloon flights around the South Pole; the ATIC-1 test flight during Dec. 2000 - Jan. 2001 and the ATIC-2 science flight during Dec. 2002 - Jan. 2003. Analyses of both datasets have, to date, relied upon the highly segmented Silicon Matrix (SiM) detector to separate the incident cosmic ray from the calorimeter backscatter and to identify the charge. This method has worked well, enabling ATIC to separate protons from helium and to resolve all the major species up through iron. This charge resolution can be significantly improved by restricting the analysis to particle trajectories that pass through two SiM pixels at the cost of using only a fraction of the potential instrument geometry. However, immediately below the SiM is the two layer SI hodoscope (x, y) consisting of Bicron BC-408 plastic scintillator 2 cm wide, 1 cm thick, 88.2 cm long strips viewed by Hamamatsu R5611 photomultiplier tubes on each end of each strip. The primary purpose of the ATIC hodoscopes is to provide a fast trigger, and each hodoscope includes two crossed layers of strips (42 per layer in the case of Sl) providing supplemental particle trajectory information. The hodoscope readout electronics were designed to provide reasonable charge resolution over the dynamic range from protons through iron. This presentation discusses the S 1 hodoscope energy deposit calibrations, examines the charge resolution possible with this detector and investigates combining the S1 and SiM charge measurements to improve the overall ATIC charge resolution while minimizing degradation of the instrument geometry.

Guzik, T. G.↗
NASA NTRS1985

Multiple Scattering Effects in de/dx-E Instruments for Isotopic Composition Studies

The development of cosmic ray telescopes to separate individual isotopes of heavy elements using the dE/dx-vs.-E technique depends on the incorporation of precise trajectory sensing elements into these systems. In typical implementations the particle trajectory is derived from a set of position measurements made prior to a particle's entering the first energy loss detector. The use of the trajectory obtained in this way to correct energy loss signals for the actual pathlengths through the detectors depends on the assumption that the particle trajectory is a straight line. In order to resolve iron isotopes the angles of incidence, theta the particle track must be known rather accurately.

Wiedenbeck, M. E.↗
NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flow fields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane normal to the wing leading edge and compared to 2-D results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that the flow field over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3-D calculation.

Potapczuk, Mark G.↗
DOE OSTI2024

Geant4 Simulations on Faraday Cup Design for PIP-II Laser Wire Scanner System

The Proton Improvement Plan-II (PIP-II) accelerator upgrade at Fermilab represents a groundbreaking leap forward in high-energy physics research. This ambitious initiative involves enhancing Fermilab's accelerator complexby replacing the current linear accelerator (linac) with awarm front end (WFE) capable of accelerating H- beamsup to 2.1 MeV. Subsequently, a superconducting linac further accelerates these beams up to 800 MeV. To preciselymeasure the transverse beam profile, a combination of traditional wire scanners at the WFE section and Laser wirescanners along the superconducting linac are planned forimplementation. This investigation is centered on refiningthe Faraday cup design for the PIP-II Laser wire scannersby utilizing GEANT4, a Monte Carlo simulation toolkit.Leveraging this method enables a comprehensive analysisof particle trajectories, energy deposition, secondary particle emission, backscattering, etc., facilitating optimizationthrough adjustments to cup geometries, materials, andplacement to maximize its efficacy in beam diagnostics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

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↗