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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗
DOE OSTI2020

Fast Dynamic Aperture Optimization with Reversal Integration

A fast method for dynamic aperture (DA) optimization of storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symme- try, a numerical forward integration di ers from its reversal. For chaotic trajectories, cumulative round-o errors are scaled, which results in an exponential growth in the di erence. The exponential e ect, intrinsically associated with the Lyapunov exponent, is a generic indicator of chaos because it represents the sensitivity of chaotic motion to an initial condition. A chaos indicator of the charged particle motion is then obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. \backward integrations." The indicator was con rmed to be ob- servable through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. The DA of the National Synchrotron Light Source II storage ring, and another test di raction-limited light source ring, were optimized using this method for the purpose of demonstration.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Computer program for prediction of the deposition of material released from fixed and rotary wing aircraft

This is a user manual for the computer code ""AGDISP'' (AGricultural DISPersal) which has been developed to predict the deposition of material released from fixed and rotary wing aircraft in a single-pass, computationally efficient manner. The formulation of the code is novel in that the mean particle trajectory and the variance about the mean resulting from turbulent fluid fluctuations are simultaneously predicted. The code presently includes the capability of assessing the influence of neutral atmospheric conditions, inviscid wake vortices, particle evaporation, plant canopy and terrain on the deposition pattern.

Teske, M. E.↗
NASA NTRS1984

Numerical studies of the deposition of material released from fixed and rotary wing aircraft

The computer code AGDISP (AGricultural DISPersal) has been developed to predict the deposition of material released from fixed and rotary wing aircraft in a single-pass, computationally efficient manner. The formulation of the code is novel in that the mean particle trajectory and the variance about the mean resulting from turbulent fluid fluctuations are simultaneously predicted. The code presently includes the capability of assessing the influence of neutral atmospheric conditions, inviscid wake vortices, particle evaporation, plant canopy and terrain on the deposition pattern. In this report, the equations governing the motion of aerially released particles are developed, including a description of the evaporation model used. A series of case studies, using AGDISP, are included.

Bilanin, A. J.↗
DOE OSTI2022

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗
DOE OSTI2022

Extracting forces from noisy dynamics in dusty plasmas

Extracting environmental forces from noisy data is a common yet challenging task in complex physical systems. Machine learning (ML) represents a robust approach to this problem, yet is mostly tested on simulated data with known parameters. Here we use supervised ML to extract the electrostatic, dissipative, and stochastic forces acting on micron-sized charged particles levitated in an argon plasma (dusty plasma). By tracking the sub-pixel motion of particles in subsequent images, we successfully estimated these forces from their random motion. The experiments contained important sources of non-Gaussian noise, such as drift and pixel-locking, representing a data mismatch from methods used to analyze simulated data with purely Gaussian noise. Our model was trained on simulated particle trajectories that included all of these artifacts, and used more than 100 dynamical and statistical features, resulting in a prediction with 50\% better accuracy than conventional methods. Lastly, in systems with two interacting particles, the model provided non-contact measurements of the particle charge and Debye length in the plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Searching for metastable particles using graph computing

The reconstruction of charged particle trajectories at the Large Hadron Collider and future colliders relies on energy depositions in sensors placed at distances ranging from a centimeter to a meter from the colliding beams. We propose a method of detecting charged particles that decay invisibly after traversing a short distance of about 25 cm inside the experimental apparatus. One of the decay products may constitute the dark matter known to be 84% of all matter at galactic and cosmological distance scales. Our method uses graph computing to cluster spacepoints recorded by two-dimensional silicon pixel sensors into mathematically-defined patterns. The algorithm may be implemented on silicon-based integrated circuits using field-programmable gate array technology to augment or replace traditional computing platforms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Sea Ice Collision Risk Assessment for Tidal Turbine Siting in Cook Inlet, Alaska

Cook Inlet has great potential for tidal stream energy development. However, the presence of drifting sea ice could create hazardous collision risks for tidal turbine farms. Before turbines can be installed in Cook Inlet, sites must be surveyed to determine how often sea ice is present, how fast it will be moved by the current, and where the trajectories of drifting sea ice will be concentrated. In this study, we use remotely sensed data to characterize the seasonal sea ice conditions in Cook Inlet, a hydrodynamic model to map the water velocities, and a particle trajectory tracking model to calculate collision risks. The resulting sea ice coverage and collision risk maps will enable tidal energy developers to choose the best locations for the deployment of tidal turbines and other offshore platforms in Cook Inlet.

16 TIDAL AND WAVE POWER↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
DOE OSTI2024

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

The PIP-II accelerator upgrade at Fermilab represents a groundbreaking leap forward in high-energy physics research. This ambitious initiative involves enhancing Fermilab's accelerator complex by replacing the current linear accelerator with a warm front end (WFE) capable of accelerating H- beams up to 2.1 MeV. Subsequently, a superconducting linac further accelerates these beams up to 800 MeV. To precisely measure the transverse beam profile, a combination of traditional wire scanners at the WFE section and Laser wire scanners along the superconducting linac are planned for implementation. This investigation centers on refining the Faraday cup design for the PIP-II Laser wire scanners by utilizing GEANT4, a Monte Carlo simulation toolkit. Leveraging this method enables a comprehensive analysis of particle trajectories, energy deposition, secondary electron emission, backscattering, etc., facilitating optimization through adjustments to cup geometries, materials, and placement to maximize its efficiency.

43 PARTICLE ACCELERATORS↗
NASA NTRS1985

Monte Carlo calculations of relativistic solar proton propagation in interplanetary space

Particle fluxes and pitch angle distributions of relativistic solar protons at 1 AU were determined by Monte Carlo calculations. The analysis covers two hours after the release of the particles from the Sun and total of eight 100000 particle trajectories were simulated. The pitch angle scattering was assumed to be isotropic ad the scattering mean free path was varied from 0.1 to 4 AU. As an application, the solar injection time and interplanetary scattering mean free path of particles that gave rise to the GLE on May, 1978 were determined. Assuming exponential form, the injection decay time was found to be about 11 minutes. The m.f.p. of pitch angle scattering during the event was about 1 AU.

Lumme, M.↗
NASA NTRS1991

Laser velocimetry measurements in non-isothermal CVD systems

Researchers at the NASA Langley Research Center are applying laser velocimetry (LV) techniques to characterize the fluid dynamics of non-isothermal flows inside fused silica chambers designed for chemical vapor deposition (CVD). Experimental issues involved in the application of LV techniques to this task include thermophoretic effects on the LV seed particles, seeding the hazardous gases, index of refraction gradients in the flow field and surrounding media, optical access, relatively low flow velocities, and analysis and presentation of sparse data. An overview of the practical difficulties these issues represent to the use of laser velocimetry instrumentation for CVD applications is given. A fundamental limitation on the application of LV techniques in non-isothermal systems is addressed which involves a measurement bias due to the presence of thermal gradients. This bias results from thermophoretic effects which cause seed particle trajectories to deviate from gas streamlines. Data from a research CVD reactor are presented which indicate that current models for the interaction of forces such as Stokes drag, inertia, gravity, and thermophoresis are not adequate to predict thermophoretic effects on particle-based velocimetry measurements in arbitrary flow configurations.

Johnson, E. J.↗
NASA NTRS1980

NASCAP modelling computations on large optics spacecraft in geosynchronous substorm environments

The NASA Charging Analyzer Program (NASCAP) is used to evaluate qualitatively the possibility of such enhanced spacecraft contamination on a conceptual version of a large satellite. The evaluation is made by computing surface voltages on the satellite due to encounters with substorm environments and then computing charged particle trajectories in the electric fields around the satellite. Particular attention is paid to the possibility of contaminants reaching a mirror surface inside a dielectric tube because this mirror represents a shielded optical surface in the satellite model used. Deposition of low energy charged particles from other parts of the spacecraft onto the mirror was found to be possible in the assumed moderate substorm environment condition. In the assumed severe substorm environment condition, however, voltage build up on the inside and edges of the dielectric tube in which the mirror is located prevents contaminants from reaching the mirror surface.

Stevens, N. J.↗
NASA NTRS1980

NASCAP modelling computations on large optics spacecraft in geosynchronous substorm environments

Satellites in geosynchronous orbits have been found to be charged to significant negative voltages during encounters with geomagnetic substorms. When satellite surfaces are charged, there is a probability of enhanced contamination from charged particles attracted back to the satellite by electrostatic forces. This could be particularly disturbing to large satellites using sensitive optical systems. In this study the NASA Charging Analyzer Program (NASCAP) is used to evaluate qualitatively the possibility of such enhanced contamination on a conceptual version of a large satellite. The evaluation is made by computing surface voltages on the satellite due to encounters with substorm environments and then computing charged-particle trajectories in the electric fields around the satellite. Particular attention is paid to the possibility of contaminants reaching a mirror surface inside a dielectric tube because this mirror represents a shielded optical surface in the satellite model used. Deposition of low energy charged particles from other parts of the spacecraft onto the mirror was found to be possible in the assumed moderate substorm environment condition. In the assumed severe substorm environment condition, however, voltage build up on the inside and edges of the dielectric tube in which the mirror is located prevents contaminants from reaching the mirror surface.

Stevens, N. J.↗
NASA NTRS1977

The signature of parallel electric fields in a collisionless plasma

The presence of a parallel electric field in the magnetosphere should be apparent through its effect on particles. An analysis is made of particle trajectories in phase space for a steady state collisionless plasma. The analysis utilizes the concepts of distribution function, total energy, and magnetic moment because all three quantities are conserved in a steady state situation in the absence of collisions. Two methods are presented for using particle data to gain information on the potential distribution. The first method involves the use of measurements of both ion and electron distribution functions and invokes the principle of quasi-neutrality. The second method involves the observation of discontinuities in velocity space which can be interpreted in terms of the effect of the electric field on the apparent loss or source cone.

Whipple, E. C., Jr.↗