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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Test particle acceleration in turbulent reconnecting magnetic fields

The effect of turbulence on particle acceleration in a MHD field was investigated by computing test particle trajectories in turbulent MHD reconnecting fields, including reconnection simulations at different magnetic Reynolds numbers. The dynamics of individual particles were investigated making it possible to examine the acceleration mechanism in great detail. It was found that turbulence influences the acceleration in two ways. It enhances the reconnection electric field while producing a stochastic electric field that gives rise to momentum diffusion; and it produces magnetic 'bubbles' and other irregularities that can temporarily trap test particles in the strong reconnection electric field for times comparable to the magnetofluid characteristic time.

Ambrosiano, John↗
NASA NTRS2019

Modeling Heatshield Erosion due to Dust Particle Impacts for a Martian Entry Vehicle

New surface damage models are presented that predict the heat shield erosion due to dust particle impact when a spacecraft enters the Martian atmosphere. Existing models were based on Apollo-era experimental data and approximate methods for tracking the dust particle trajectories through the shock layer. These legacy methods will be compared against new results based on more sophisticated particle tracking methods and recent experimental data.

Planetary entry↗
NASA NTRS2020

Modeling Heatshield Erosion due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗
NASA NTRS2020

Modeling Heatshield Erosion Due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗
NASA NTRS1990

The atmosphere of 2060 Chiron

An explanation for 2060 Chiron's behavior, which focuses on the influence of Chiron's mass on the development of its dust coma, is presented. It is suggested that dust is entrained by the flow of CO or another gas of similar volatility from an active region. It remains gravitationally bound on orbits confined to a region, roughly 5000 km in extent, that lies between the surface and an exopause imposed by radiation pressure forces. The influence of radiation pressure transforms the initial particle trajectories into satellite orbits with a characteristic period of 20 days and orbital residence time of about 25 revolutions. The particle population in the coma slowly increases, explaining Chiron's photometric behavior.

Meech, Karen J.↗
NASA NTRS2001

The Capture of Interstellar Dust: The Pure Poynting-Robertson Case

Ulysses and Galileo spacecraft have discovered interstellar dust particles entering the solar system. In general, particles trajectories not altered by Lorentz forces or radiation pressure should encounter the sun on open orbits. Under Newtonian forces alone these particles return to the interstellar medium. Dissipative forces, such as Poynting Robertson (PR) and corpuscular drag and non-dissipative Lorentz forces can modify open orbits to become closed. In particular, it is possible for the orbits of particles that pass close to the Sun to become closed due to PR drag. Further, solar irradiation will cause modification of the size of the dust particle by evaporation. The combination of these processes gives rise a class of capture orbits and bound orbits with evaporation. Considering only the case of pure PR drag a minimum impact parameter is derived for initial capture by Poynting-Robertson drag. Orbits in the solar radiation field are computed numerically accounting for evaporation with optical and material properties for ideal interstellar particles modeled. The properties of this kind of particle capture are discussed for the Sun but is applicable to other stars.

Jackson, A. A.↗
DOE OSTI2024

A Stochastic Calculus Approach to Boltzmann Transport

Traditional Monte Carlo methods for particle transport utilize source iteration to express the solution, the flux density, of the transport equation as a Neumann series. Our contribution is to show that the particle paths simulated within source iteration are associated with the adjoint flux density and the adjoint particle paths are associated with the flux density. Here, we make our assertion rigorous through the use of stochastic calculus by representing the particle path used in source iteration as a solution to a stochastic differential equation (SDE). The solution to the adjoint Boltzmann equation is then expressed in terms of the same SDE, and the solution to the Boltzmann equation is expressed in terms of the SDE associated with the adjoint particle process. An important consequence is that the particle paths used within source iteration simultaneously provide Monte Carlo samples of the flux density and adjoint flux density in the detector and source regions, respectively. The significant practical implication is that particle trajectories can be reused to obtain both forward and adjoint quantities of interest. To the best our knowledge, the reuse of entire particles paths has not appeared in the literature. Monte Carlo simulations are presented to support the reuse of the particle paths.

Boltzmann transport↗
NASA NTRS2000

Distribution of Dust from Kuiper Belt Objects

Using an efficient computational approach, we have reconstructed the structure of the dust cloud in the Solar system between 0.5 and 100 AU produced by the Kuiper belt objects. Our simulations offer a 3-D physical model of the 'kuiperoidal' dust cloud based on the distribution of 280 dust particle trajectories produced by 100 known Kuiper belt objects; the resulting 3-D grid consists of 1.9 x 10' cells containing 1.2 x 10" particle positions. The following processes that influence the dust particle dynamics are taken into account: 1) gravitational scattering on the eight planets (neglecting Pluto); 2) planetary resonances; 3) radiation pressure; and 4) the Poynting-Robertson (P-R) and solar wind drags. We find the dust distribution highly non-uniform: there is a minimum in the kuiperoidal dust between Mars and Jupiter, after which both the column and number densities of kuiperoidal dust sharply increase with heliocentric distance between 5 and 10 AU, and then form a plateau between 10 and 50 AU. Between 25 and 45 AU, there is an appreciable concentration of kuiperoidal dust in the form of a broad belt of mostly resonant particles associated with Neptune. In fact, each giant planet possesses its own circumsolar dust belt consisting of both resonant and gravitationally scattered particles. As with the cometary belts simulated in our related papers, we reveal a rich and sophisticated resonant structure of the dust belts containing families of resonant peaks and gaps. An important result is that both the column and number dust density are more or less flat between 10 and 50 AU, which might explain the surprising data obtained by Pioneers 10 & 11 and Voyager that the dust number density remains approximately distance-independent in this region. The simulated kuiperoidal dust, in addition to asteroidal and cometary dust, might represent a third possible source of the zodiacal light in the Solar system.

Gorkavyi, Nick N.↗
NASA NTRS2012

Kinetic Description of Ionospheric Outflows Based on the Exact Form of Fokker-Planck Collision Operator: Electrons

We present the results of a finite difference implementation of the kinetic Fokker-Planck model with an exact form of the nonlinear collisional operator, The model is time dependent and three-dimensional; one spatial dimension and two in velocity space. The spatial dimension is aligned with the local magnetic field, and the velocity space is defined by the magnitude of the velocity and the cosine of pitch angle. An important new feature of model, the concept of integration along the particle trajectories, is discussed in detail. Integration along the trajectories combined with the operator time splitting technique results in a solution scheme which accurately accounts for both the fast convection of the particles along the magnetic field lines and relatively slow collisional process. We present several tests of the model's performance and also discuss simulation results of the evolution of the plasma distribution for realistic conditions in Earth's plasmasphere under different scenarios.

Khazanov, George V.↗
NASA NTRS1983

Evaluation of a stochastic model of particle dispersion in a turbulent round jet

A stochastic model of particle dispersion by turbulence, proposed by Gosman and Ioannides (1981), is evaluated. The method employs a k-epsilon model to estimate turbulence properties. Dispersion is determined by computing particle motion, with random sampling to obtain instantaneous flow properties for a statistically significant number of particle trajectories. The stochastic model yields good results particularly when eddy lifetimes are evaluated. The method allows the effects of large relative velocities between the particles and the flow, drag properties at Reynolds numbers greater than the Stokes flow regime, and the variations of local turbulence properties to be readily handled, at least for boundary layer flows.

Shuen, J.-S.↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2020

Fast dynamic aperture optimization with forward-reversal integration

A fast dynamic aperture (DA) optimization method for storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symmetry, a numerical forward integration differs from its reversal. For a chaotic trajectory, cumulative round-off errors are scaled, which results in an exponential growth on the difference. The exponential effect is a generic chaos indicator which represents the sensitivity of the chaotic motion to its initial condition. The chaos indicator of the charged particle motion can be obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. “backward integrations.” The indicator is observable even through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. Finally, the DA of the National Synchrotron Light Source II storage ring, and another test diffraction-limited light source ring, were optimized using this method for the purpose of demonstration.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

A theoretical model of current collection

Numerical simulations were carried out on the high-bias voltage on the conducting satellite for the current collection. Results show that the particle trajectories are dependent on the geomagnetic field orientation and the biased potential on both the satellite and the instrument's detector. It is noted that these particles will not escape from the plasma sheath when the B field is perpendicular to the satellite surface plane. It is also shown that part of the satellite surface will be shielded if the perpendicular magnetic field is greater than 0.35 gauss.

Hwang, Kai-Shen↗
DOE OSTI2025

A data-driven method to estimate the antiproton background in the Mu2e experiment

The Mu2e experiment at Fermilab will search for the Charged Lepton Flavour Violating (CLFV) process of coherent, neutrinoless µ− → e − conversion in the field of an aluminum nucleus. The expected signal is a monochromatic electron with the energy of 104.97 MeV, slightly below the muon rest mass. Observation of a CLFV process would provide unambiguous evidence for Beyond the Standard Model (BSM) physics. Mu2e is sensitive to a wide range of BSM models and has the capability to distinguish between them, guiding us towards the most accurate models. The key features of the Mu2e experiment are: (1) a high intensity pulsed negative muon beam with about 1010 stopped µ −/s, and (2) a sophisticated superconducting solenoid system with a gradient magnetic field to form and guide the intense muon beam to the target. The Mu2e physics data taking is expected to begin in 2027. For Run I, the expected 5σ discovery sensitivity is Rµe = 1.2 × 10−15, with a total expected background of 0.11 ± 0.03 events. In the absence of a signal, the expected upper limit is Rµe < 6.2 × 10−16 at 90% CL. The success of this experiment hinges on the accurate estimation of the background from various SM processes that could provide signal-like electrons. One of the background processes is antiprotons annihilating in the stopping target to produce signal like electrons through π0 → γγ decays followed by γ conversions, and π− → µ−ν¯ decays followed by µ− decay. It is a relatively small background with large uncertainty (100%) due to the lack of antiproton production cross section information for the Mu2e proton beam energy of 8 GeV. We have developed a novel methodology to estimate the antiproton background in-situ. This forms the main theme of the thesis. We observed that at Mu2e energies, antiproton annihilation in the stopping target is the only source of events with multiple, simultaneous particle trajectories. From Geant4 simulations, only about 0.2% of the simulated antiproton annihilation events have a signal-like electron. Meanwhile, ∼ 5% of events have multiple reconstructible particle tracks per event. Therefore, we have devised a methodology to reconstruct the multi-track events and estimate the antiproton background by exploiting the large ratio of the production rates of the two final states.

Chithirasreemadam, Namitha [Pisa U.] (ORCID:000000↗
NASA NTRS1987

Auroral particles

The problems concerning the aurora posed prior to the war are now either solved in principle or were restated in a more fundamental form. The pre-war hypothesis concerning the nature of the auroral particles and their energies was fully confirmed, with the exception that helium and oxygen ions were identified as participating in the auroral particle precipitation in addition to the protons. The nature of the near-Earth energization processes affecting auroral particles was clarified. Charged particle trajectories in various electric field geometries were modeled. The physical problems have now moved from determining the nature and geometry of the electric fields, which accelerate charged particles near the Earth, to accounting for the existence of these electric fields as a natural consequence of the solar wind's interaction with Earth. Ultimately the reward in continuing the work in auroral and magnetospheric particle dynamics will be a deeper understanding of the subtleties of classical electricity and magnetism as applied to situations not blessed with well-defined and invariant geometries.

Evans, David S.↗
NASA NTRSDate not supplied

Characterization of Naturally-Occurring Particles in the Nasa Langley 0.3-m Transonic Cryogenic Tunnel Using PTV

Naturally-occurring particles within the NASA Langley 0.3-m Transonic Cryogenic Tunnel are characterized for their aerodynamic performance using particle tracking velocimetry. Two sets of experiments were conducted to observe different behaviors of the particles. The normal shockwave emanating from the top surface of a supercritical airfoil was used to induce velocity lag in the particles, and the subsequent spatial decay of the velocity was used to estimate the effective diameter of the particles. Mean particle diameters between 1.6 and 1.9 𝝁m were measured, with sizes ranging from 0.2 to 3.5 𝝁m over the entire ensemble. The response of particles to separated flow was investigated. By operating in “high-lift” (low Mach number, high angle of attack) conditions with a semi-span airfoil, the ability of particles to detect separated flow on the upper surface of the airfoil was assessed. Transition from fully attached flow to fully separated flow was observed on the top surface of the airfoil accompanying a variation of angle of attack from 8° to 12°. Examination of velocity distributions indicates less than 10 percent of particle trajectories did not respond to the regions of separated flow. These results are promising, but further facility-specific work is needed to answer the broader question of particle tracking reliability.

Transonic↗
NASA NTRSDate not supplied

Characterization of Naturally-Occurring Particles in the NASA Langley 0.3-m Transonic Cryogenic Tunnel Using PTV

Naturally-occurring particles within the NASA Langley 0.3-m Transonic Cryogenic Tunnel are characterized for their aerodynamic performance using particle tracking velocimetry. Two sets of experiments were conducted to observe different behaviors of the particles. The normal shockwave emanating from the top surface of a supercritical airfoil was used to induce velocity lag in the particles, and the subsequent spatial decay of the velocity was used to estimate the effective diameter of the particles. Mean particle diameters between 1.6 and 1.9 𝝁m were measured, with sizes ranging from 0.2 to 3.5 𝝁m over the entire ensemble. The response of particles to separated flow was investigated. By operating in “high-lift” (low Mach number, high angle of attack) conditions with a semi-span airfoil, the ability of particles to detect separated flow on the upper surface of the airfoil was assessed. Transition from fully attached flow to fully separated flow was observed on the top surface of the airfoil accompanying a variation of angle of attack from 8° to 12°. Examination of velocity distributions indicates less than 10 percent of particle trajectories did not respond to the regions of separated flow. These results are promising, but further facility-specific work is needed to answer the broader question of particle tracking reliability.

Transonic↗
DOE OSTI2021

In-situ magnetic alignment model for additive manufacturing of anisotropic bonded magnets

In this work, we report a mathematical framework which predicts the degree of alignment (DoA) in an in-situ aligned additively manufactured 3D printed bonded magnets. A multiphysics model is developed which couples the harmonious interactions of magnetic particles in a viscous flowing polymer under the presence of an external magnetic field. The hydrodynamic fluid-particle interaction is paired with the magnetophoretic force to predict the particle trajectory and distribution during extrusion through a 3D printer nozzle. Succeeding the force balance, a magnetohydrodynamic torque equilibrium analysis is performed to predict the net-orientation of the magnetic particles as a function of the applied field strength, viscous forces, and particle-to-particle interactions (P2P). Experimental validation of the DoA predictions is performed using 65 vol% Nd-Fe-B+Sm-Fe-N in Nylon12 (DoA exp = 0.620 and DoA theory = 0.686) and 15 vol% Sm-Co in PLA (DoA exp = 0.830 and DoA theory = 0.863). A parametric analysis is performed to analyze the effect of operating and design parameters like alignment field strength, magnetic loading fraction, extrusion load, and particle size. The model predicts a competing behavior between particle-fluid and particle-particle interactions under the presence of an applied field. Lastly, the model provides a framework to efficiently predict the DoA in tandem with a functionalized-magnetic 3D printer and allows the user to adjust the operating parameters according to the desired DoA.

36 MATERIALS SCIENCE↗