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

An analysis of CCD images of the coma of Comet Halley

The analysis of selected CCD images of the coma of comet P/Halley is presented. The images were taken using specially designed filters that isolate regions of a comet's spectrum such that only sunlight which has been scattered by the dust in the coma is recorded. The modeling analysis objective is to make use of the skills developed in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae and to use those models as a basis for a new dust coma model. This model will include a self-consistant picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. The model is intended to be used as a tool to analyze selected images from the two sets of CCD images with the hope that it will help the understanding of the effects of a number of important processes on the spatial morphology of the observed dust coma. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observe scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; and (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gasdynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology on the outer dust coma.

Combi, Michael↗

Analysis of CCD images of the coma of comet P/Halley

The modeling analysis objective of this project is to make use of the skill acquired in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae as a basis for a new dust coma model. This model will include a self-consistent picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. Our purpose is to use this model as a tool to analyze selected images from two sets of data of the comet P/Halley with the hope that we can help to understand the effects of a number of important processes on the spatial morphology of the observed dust coma. The study will proceed much in the same way as our study of the spatially extended hydrogen coma where we were able to understand the spatial morphology of the Lyman-alpha coma in terms of the partial thermalization of the hot H atoms produced by the photodissociation of cometary H2O and OH. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observed scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gas dynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology of the outer dust coma; and (6) the time and direction dependence of the source of dust.

Combi, Michael R.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SBND Analysis using ML Reconstruction Chain

As part of the Short Baseline Neutrino (SBN) Program at Fermilab, the Short Baseline Near Detector (SBND) is positioned in the Booster Neutrino Beam (BNB) and explores neutrino-argon interactions with unprecedented statistics. SBND is a Liquid Argon Time Projection Chamber (LArTPC). Electrons produced through ionization drift toward three wire planes, providing signals that form 2D images of particle trajectories. I introduce the Scalable Particle Imaging using Neural Embeddings (SPINE) framework, which employs a Machine Learning (ML)-based 3D reconstruction using a series of neural networks. Here, we present SPINE’s reconstruction chain, analysis approaches, and results from our latest simulation samples.

43 PARTICLE ACCELERATORS↗

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗

Effects of solid particles suspended in fluid flow through an axial flow compressor stage

An approximate method for calculating the flow properties of gas-particle mixture flowing over blades in a cascade is studied. Using an analytical method, the solid particle trajectory and location of collisions between the solid particles and the blade surfaces are determined. In addition an experimental investigation of the trajectories and velocities of solid particles suspended in a fluid passing through an axial flow compressor cascade was performed. The cascade blades were made of 2024 aluminum alloy and the solid particles used were quartz sand with average diameter of 165 microns. The following parameters were investigated: blade pressure distribution, and total pressure loss coefficient, both depending on the degree of erosion. In addition, a theoretical estimation of the blade erosion is compared with the experimental data.

Tabakoff, W.↗

Erosion in radial inflow turbines. Volume 3: Trajectories of erosive particles in radial inflow turbines

The theoretical trajectories that erosive particles follow in the gas flow fields of a typical radial inflow turbine were investigated. A discussion of the theoretical trajectories that the particles follow in the scroll, in the nozzles, in the vortex between the nozzles and the rotor, and in the rotor passages is included. The results are presented in terms of the characteristic length, a similarity parameter which relates the particles that follow the same trajectory in equivalent flow fields. For Vol, 1, see N74-19395.

Clevenger, W. B., Jr.↗

Cosmic ray particle dosimetry and trajectory tracing

Five pocket mice (Perognathus longimembris) were flown on Apollo XVII, each with a solid-state (plastic) nuclear track detector implanted beneath its scalp. The subscalp detectors were sensitive to HZE cosmic ray particles with a LET greater than or approximately equal to 0.15 million electron volts per micrometer (MeV/micron). A critical aspect of the dosimetry of the experiment involved tracing individual particle trajectories through each mouse head from particle tracks registered in the individual subscalp detectors, thereby establishing a one-to-one correspondence between a trajectory location in the tissue and the presence or absence of a lesion. The other major aspect was the identification of each registered particle. An average of 16 particles with Z greater than or equal to 6 and 2.2 particles with Z greater than or equal to 20 were found per detector. The track density, 29 tracks/sq cm, when adjusted for detection volume, was in agreement with the photographic emulsion data from an area dosimeter located next to the flight package.

Cruty, M. R.↗

System Configuration Evaluation for Process Settling of Hanford Waste Solid Particles

Direct Feed High-Level Waste (DFHLW) is a potential flowsheet operations approach to initiating high-level waste (HLW) vitrification prior to completion of the Hanford Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. A settle/decant process has been proposed to concentrate solids prior to delivery to the WTP HLW Facility during DFHLW operations, wherein the solids in a settled layer would be remixed with the supernatant liquid remaining after decanting operations to provide the feed at required solids concentrations. Settling would be used in lieu of purpose-built filtration or other solids separation equipment. Pacific Northwest National Laboratory (PNNL) is providing baseline technical support to the Washington River Protection Solutions (WRPS) Flowsheet Integration group. To support planning for DFHLW, WRPS previously requested that PNNL evaluate the current data set available to predict the time needed for HLW solids to settle and the solids concentration and strength of that settled layer, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling times. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers. In addition to the data gaps, an overarching observation was made that there is significant variation in behavior of settling rate and settled layer data. The settling time required to concentrate solids via a settle/decant process was determined from the limited data to have a difference of potentially more than a factor of 5,000 in the estimated settling times, varying from 0.2 to 1,060 days for example depending on process vessel depth and final sediment solids concentration. In contrast, successful processes of liquid-forward output streams resulting from in-tank settling and decanting forward liquid have been reported for operations conducted at the Hanford Site. The purpose of this current report is to further support DFHLW planning by evaluating double-shell tank (DST) and alternate vessel equipment and operational configurations to enable optimization of the settle/decant process to concentrate solids. Hanford waste processing behavior specific to liquid feed availability following a slurry transfer in a DST is summarized, including process stream characteristics and process equipment configurations. The performance of DST process equipment configurations is evaluated for possible improvements using computational fluid dynamics (CFD) and simple analytical models. Potential new vessel design(s) specific to enabling effective settle/decant processes, and cursory summary of other separate and inline solids separations processes, are also provided. The CFD results indicated that improvement in outflow solids concentration was promoted by a reduction in the slurry flow rate, angling the distributor nozzles downward, and lifting the transfer pump. The solid-liquid analysis evaluating particle trajectory confirmed that the potential for particle ingestion (in the transfer pump) was decreased with increased radial separation between the inlet and outlet (transfer pump inlet), decreased inlet flow, and decreased liquid density and viscosity for a neutrally buoyant inlet flow. An assessment was also made of the potential for inflow configuration changes to result in the discrete mounding or piling of solids within the tank. Based on the characterization of the settled waste to date, HLW sediments will be unlikely to sustain a substantial angle of repose to facilitate significant variations in the elevation of the settled solids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultrarelativistic regime with a bulk Lorentz factor of up to γ 0 = 10 6 . Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. Similar to earlier studies, we find that the resulting energy spectra crucially depend on λ/d e (λ is the wavelength of the striped wind and de is the relativistic plasma skin depth). When λ/d e is large (λ ≳ 40d e ), the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ 0 and σ 0 increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor of γ0 ≈ 106 and a magnetization parameter of σ 0 = 10, which can explain the recent observations of high-energy gamma rays from PWNe

79 ASTRONOMY AND ASTROPHYSICS↗