Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle method”

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 1,765 records · Page 98

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗

Transmission measurement based on STM observation to detect the penetration depth of low-energy heavy ions in botanic samples

The penetration depth of low-energy heavy ions in botanic samples was detected with a new transmission measurement. In the measurement, highly oriented pyrolytic graphite (HOPG) pieces were placed behind the botanic samples with certain thickness. During the irradiation of heavy ions with energy of tens of keV, the energetic particles transmitted from those samples were received by the HOPG pieces. After irradiation, scanning tunneling microscope (STM) was applied to observe protrusion-like damage induced by these transmitted ions on the surface of the HOPG. The statistical average number density of protrusions and the minimum transmission rate of the low-energy heavy ions can be obtained. The detection efficiency of the new method for low-energy heavy ions was about 0.1-1 and the background in the measurement can be reduced to as low as 1.0 x 10(8) protrusions/cm2. With this method, the penetration depth of the energetic particles was detected to be no less than 60 micrometers in kidney bean slices when the slices were irradiated by 100 keVAr+ ion at the fluence of 5 x 10(16) ions/cm2. c2002 Elsevier Science Ltd. All rights reserved.

DNA Damage↗

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

Abstract The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s −1 and diffusion coefficient of 0.05 m 2 s −1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Absorption of radiation belt particles by the inner satellites of Jupiter

The study of trapped particle absorption by the inner Jovian satellites is reviewed from the viewpoint of radiation belt physics. Both pre- and post-Pioneer work is discussed but the emphasis is on methods used to deduce radial diffusion coefficients of particle transport from particle data. The phenomenon of particle absorption as observed by experiments on Pioneers 10 and 11 is considered; absorption effects are found to depend on the satellite, and on particle energy and species. Approximate diffusion coefficients derived from the data are found to follow a steeper spatial dependence than previously expected. The assumptions and limitations of absorption analysis and diffusion coefficient estimation are pointed out.

Mogro-Campero, A.↗

Large-scale frictionless jamming with power-law particle size distributions

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Here, invoking a recently developed neighbor binning algorithm, we generate mechanically stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly 4 000 000 particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of nonrattler particles attains the frictionless isostatic value of six in all cases. The size distribution of nonrattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparison of MCNP thermal scattering laws to inelastic neutron scattering data

Neutron transport computations employ thermal scattering laws (TSLs) to account for the influence of material structure and collective excitations at low energy (≲ 10 eV). Recent work has been undertaken to improve TSL accuracy and expand the number of available materials in the Evaluated Nuclear Data File (ENDF). However, comparison of TSL to measured data is often sparse. In this work, we present a method to evaluate TSL accuracy in the Monte Carlo N-Particle (MCNP) code by simulating inelastic neutron scattering (INS) measurement and comparing to data. The method is briefly described with example application to the recently released ENDF-VIII TSL library for polyethylene as well as a cryogenic PE TSL we are separately developing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On the acceleration of charged particles at relativistic shock fronts

The diffusive acceleration of highly relativistic particles at a shock is reconsidered. Using the same physical assumptions as Blandford and Ostriker (1978), but dropping the restriction to nonrelativistic shock velocities, the authors find approximate solutions of the particle kinetic equation by generalizing the diffusion approximation to higher order terms in the anisotropy of the particle distribution. The general solution of the transport equation on either side of the shock is constructed, which involves the solution of an eigenvalue problem. By matching the two solutions at the shock, the spectral index of the resulting power law is found by taking into account a sufficiently large number of eigenfunctions. Low-order truncation corresponds to the standard diffusion approximation and to a somewhat more general method described by Peacock (1981). In addition to the energy spectrum, the method yields the angular distribution of the particles and its spatial dependence.

Kirk, J. G.↗

Particle size and velocity measurement by laser anemometry

The use of a forward scattering Laser Anemometer system for the measurement of particle diameters larger than the fringe spacing is investigated. A systematic calibration experiment demonstrates a direct relationship between signal amplitude and particle diameter for transparent particles between 30 microns and 240 microns in diameter. A light scattering analysis confirms these measurements and indicates how the range of particles which can be measured using the peak (mean) signal depends on the optical arrangement and particle properties. On the basis of this work a real time Laser Anemometer signal processor is modified so that the peak of the mean of each signal can be measured. A method of correcting data for the effect of particles which do not go through the center of the measurement volume is derived. Good agreement is found between particle size distributions measured in sprays by using the Laser Anemometer and size distributions measured by collecting the particles on a slide and using an image analysis computer.

Yule, A. J.↗

Tracking of Marangoni driven motion during laser powder bed fusion

Laser Powder Bed Fusion (LPBF) can be used to create Functionally Graded Materials (FGMs) by selectively distributing property augmenting particles throughout the build matrix. With this method parts can be created with nonhomogeneous properties tuned to spatially varied, pre-determined design constraints. To unlock this capability, the fluid motion that occurs in the molten liquid phase must be thoroughly understood. More importantly, the actual tracking of motion during the LPBF process is needed. Further, a numerical model of heat transfer and fluid flow in LPBF process was created. Using results from the model the motion of particles starting at many different locations was tracked. The average change in position for these particles based on starting location was calculated. The effects of non-dimensional parameters on melt pool size were examined. These results were validated experimentally and compared to what is available in the literature.

42 ENGINEERING↗

Soot Optical Property Study

Recent past studies of soot reaction processes in laminar premixed and nonpremixed flames generally have used the intrusive technique of thermophoretic sampling and analysis by transmission electron microscopy (TEM) to observe soot structure and obtain important fundamental information about soot particle properties, such as soot primary particle diameters, the rate of change of soot primary particle diameter as a function of time (or rate of soot surface growth or oxidation), the amount of soot particle reactive surface area per unit volume, the number of primary soot particles per unit volume, and the rate of formation of primary soot particles (or the rate of soot primary particle nucleation). Given the soot volume per unit volume of the flame (or the soot volume fraction), all these properties are readily found from a measurement of the soot primary particle diameter (which usually is nearly a constant for each location within a laminar flame). This approach is not possible within freely propagating flames, however, because soot properties at given positions in such flames vary relatively rapidly as a function of time in the soot formation and oxidation regions compared to the relatively lengthy sampling times needed to accumulate adequate soot samples and to minimize effects of soot collected on the sampling grid as it moves to and from the sampling position through other portions of the flame. Thus, nonintrusive optical methods must be used to find the soot primary particle diameters needed to define the soot surface reaction properties mentioned earlier. Unfortunately, approximate nonintrusive methods used during early studies of soot reaction properties in flames, found from laser scattering and absorption measurements analyzed assuming either Rayleigh scattering or Mie scattering from polydisperse effective soot particles having the same mass of soot as individual soot aggregates, have not been found to be an effective way to estimate the soot surface reaction area per unit volume. Thus, alternative nonintrusive optical methods of finding these properties must be sought, which was the objective of this phase of the investigation. The alternative method used here involves use of the Rayleigh-Debye-Gans-Polydisperse-Fractal-Aggregate (RDG-PFA) scattering approximation for soot aggregates in flames. Thus, the development of this method will be discussed next before describing its evaluation as a means of nonintrusively measuring soot primary particle diameters in soot-containing flames.

Aung, K. T.↗

Estimation of Sea Spray Aerosol Surface Area Over the Southern Ocean Using Scattering Measurements

This study focuses on methods to estimate dry marine aerosol surface area (SA) from bulk optical measurements. Aerosol SA is used in many models' ice nucleating particle (INP) parameterizations, as well as influencing particle light scattering, hygroscopic growth, and reactivity, but direct observations are scarce in the Southern Ocean (SO). Two campaigns jointly conducted in austral summer 2018 provided co-located measurements of aerosol SA from particle size distributions and lidar to evaluate SA estimation methods in this region. Mie theory calculations based on measured size distributions were used to test a proposed approximation for dry aerosol SA, which relies on estimating effective scattering efficiency (Q) as a function of Ångström exponent (å). For distributions with dry å < 1, Q = 2 was found to be a good approximation within ±50%, but for distributions with dry å > 1, an assumption of Q = 3 as in some prior studies underestimates dry aerosol SA by a factor of 2 or more. We propose a new relationship between dry å and Q, which can be used for –0.2 < å < 2, and suggest å = 0.8 as the cutoff between primary and secondary marine aerosol-dominated distributions. Application of a published methodology to retrieve dry marine aerosol SA from lidar extinction profiles overestimated aerosol SA by a factor of 3–5 during these campaigns. Using Microtops aerosol optical thickness measurements, we derive alternative lidar conversion parameters from our observations, applicable to marine aerosol over the SO.

54 ENVIRONMENTAL SCIENCES↗

Fast nonlinear iterative solver for an implicit, energy-conserving, asymptotic-preserving charged-particle orbit integrator

Here recently, an asymptotic-preserving (AP) particle orbit integrator has been proposed with remarkable properties including exact energy conservation, the ability to capture of all first-order drifts (including the ∇B-drift), the ability to capture trapped-passing boundaries with parallel velocity extremely close to the critical velocity, and the ability to transition from strongly to weakly magnetized spatial regions. The new AP orbit integrator is implicit, employing a Crank-Nicolson (CN) temporal discretization to ensure exact energy conservation. This, in turn, requires a local nonlinear iteration involving particle velocities and positions, and the local electromagnetic fields, to obtain the new-time solution. Ref. [1] did not attempt to provide an efficient solver for this system, and employed a brute-force GMRES-driven Jacobian-free Newton-Krylov (JFNK) solver to invert the particle orbit equations at every timestep for expediency. While JFNK is robust and reliable, it is also expensive and very intrusive for practical implementations of the method (it requires having the JFNK machinery available and solving a 6 x 6 Jacobian system iteratively once per iteration per particle).

97 MATHEMATICS AND COMPUTING↗

Particle resonances in toroidal fusion devices

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through avalanche. The existence of a mode particle resonance depends on properties of the equilibrium, particle trajectories, and perturbation mode harmonic content. Resonance location is a function of particle energy and equilibrium field line helicity. Different methods for finding resonance location and energy dependence are developed. The properties of resonances are discussed using examples from magnetic fusion devices. Here, we show that if mode resonances exist at low particle energy, they very likely also exist at high energy, thus modifying high energy beam particles and fusion products. It is possible for a resonance to appear due to mode induced orbit helicity modification when it is forbidden at low mode amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Seeding techniques used in the Vortex Research Facility

For many years, the Vortex Research Facility used kerosene vapor as seeding particles. The kerosene vapor formed (approx. 1 micron) particles initially. After injection into the test section, the vapor would begin to condense and form larger particles. By switching to solid particles, this has been eliminated. The uniform-sized (approx. 1 micron) solid particle seeding is used both to seed the Laser Velocimeter system and for flow visualization. After injection, the solid particles reach a uniform suspension before the vehicle is launched. It is found that the kerosene vapor poses a problem because it is a heated gas which modified the test section temperature distribution. It also produces high turbulence and secondary air currents because it was continually injected from the time of vehicle launch until model passage. In cases where the effects of either turbulence or the vertical temperature distribution are of interest, this hot gas injection method is unsuitable. Both types of seeding particles provide adequate flow visualization but are sensitive to vertical temperature distributions.

Satran, D. R.↗

Towards Agentic AI on Particle Accelerators

As particle accelerators grow in complexity, traditional control methods face increasing challenges in achieving optimal performance. This paper envisions a paradigm shift: a decentralized multi-agent framework for accelerator control, powered by Large Language Models (LLMs) and distributed among autonomous agents. We present a proposition of a self-improving decentralized system where intelligent agents handle high-level tasks and communication and each agent is specialized control individual accelerator components. This approach raises some questions: What are the future applications of AI in particle accelerators? How can we implement an autonomous complex system such as a particle accelerator where agents gradually improve through experience and human feedback? What are the implications of integrating a human-in-the-loop component for labeling operational data and providing expert guidance? We show two examples, where we demonstrate viability of such architecture.

43 PARTICLE ACCELERATORS↗

Some statistical analysis techniques for the effect of systematic instrument errors on the results of measurements in a space experiment

Mathematical models of certain systematic telemetry signal errors arising during the passage of the signal through the recording and transmitting blocks of the physical instruments aboard the spacecraft are classified and constructed. Statistical analysis methods are used to determine the background function and relative sensitivity coefficients of the channels of a charged particle detector. Results are presented. A block diagram for a method taking into account the effect of distortions is given.

Kondakov, S. P.↗