Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space distribution”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Quantitative characterization of spatial distribution of particles in materials: Application to materials processing

Most engineering materials contain second phase particles or fibers which serve to reinforce the matrix phase. The effect of reinforcements on material properties is usually analyzed in terms of the average volume fraction and spacing of reinforcements, quantities which are global microstructural characteristics. However, material properties can also depend on local microstructural characteristics; for example, on how uniformly the reinforcing phase is distributed in the material. The analysis method will then be applied to a materials processing problem to discover how processing parameters can be selected to maximize redistribution of the reinforcing phase during processing. Several mathematical analysis methods could be adapted to the problem of characterizing the distribution of particles in materials. A tessellation-based method was selected. In the first phase of the investigation, a software package was written to automate the analysis. Typical results are shown. The analysis technique allows the degree to which particles are clustered together, the size and spacing of particle clusters, and the particle density in clusters to be found. The analysis methods were applied to computer-generated distributions and to a few real particle-containing materials. Methods for analyzing a nonuniform particle distribution in a material can be applied to two broad classes of materials science problems: understanding how the resulting particle distribution affects properties. The analysis method described is applied to a materials processing problem: how to select extrusion conditions to maximize the redistribution of reinforcing particles that are initially nonuniformly distributed. In addition, the tessellation-based method to analyze star distributions in spiral galaxies was adapted, illustrating the diverse types of problems to which the analysis method can be applied.

Parse, J. B.↗
DOE OSTI2023

Visualization of fast ion phase-space flow in plasmas well-below, near and well-above Alfvén eigenmode stability threshold in tokamak

Abstract Transport of fast ions along certain local phase space paths, referred to as fast ion phase-space flow, has been systematically measured by an imaging neutral particle analyzer (INPA) in three plasma regimes, which are well-below, near, and well-above the Alfvén eigenmode stability threshold (Du et al 2021 Phys. Rev. Lett. 127 235002). (1) In plasmas well-below the Alfvén eignenmodes (AE) stability threshold, fast ions are well-confined on passing particle orbits without noticeable transport over the phase space. The observed INPA images agree well with the synthetic INPA images, using the fast ion distribution predicted by neoclassical theory. (2) In plasmas near the AE stability threshold, INPA images in the presence of AE activity moderately deviate from those without AE activity. The image difference can be well interpreted by AE-driven, phase-space fast ion flow. Paths of this flow over the velocity space (or streamlines) are reconstructed by the intersection lines of curved E ʹ and µ surfaces, referred to as E ʹ and µ line (where E ′ ≡ E − ( ω / n ) P ζ ; E , P ζ and µ are the energy, canonical toroidal momentum and magnetic moment of ions; ω and n are the angular frequency and toroidal mode number of AEs, respectively). Resonant fast ions move radially inward by gaining energy and move radially outward by losing energy and the trajectory well aligns with the E ʹ and µ lines that pass through the mode resonances near the injection energy of neutral beams. (3) In plasmas well-above the AE stability threshold, fast ion phase-space dynamics shows additional features. Fast ions are transported out of the birth positions so promptly along the streamlines that the slowing-down process from the injection energy is not observable, exhibiting strong critical gradient behavior at local phase space. As a result, the increase of electron temperature is very small, in spite of an increase of beam power by ∼ 45 % . It should be emphasized that the directions of phase-space transport, induced by AEs with different frequencies, structures and mode numbers, do not largely differ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gluon-mediated inclusive Deep Inelastic Scattering from Regge to Bjorken kinematics

We revisit high energy factorization for gluon mediated inclusive Deep Inelastic Scattering (DIS) for which we propose a new semi-classical approach that accounts systematically for the longitudinal extent of the target in contrast with the shockwave limit. In this framework, based on a partial twist expansion, we derive a factorization formula that involves a new gauge invariant unintegrated gluon distribution which depends explicitly on the Feynman x variable. It is shown that both the Regge and Bjorken limits are recovered in this approach. We reproduce in particular the full one loop inclusive DIS cross-section in the leading twist approximation and the all-twist dipole factorization formula in the strict x = 0 limit. Although quantum evolution is not discussed explicitly in this work, we argue that the proper treatment of the x dependence of the gluon distribution encompasses the kinematic constraint that must be imposed on the phase-space of gluon fluctuations in the target to ensure stability of small-x evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Ulysses field and plasma observations of magnetic holes in the solar wind and their relation to mirror-mode structures

The term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

Winterhalter, Daniel↗
NASA NTRS2002

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗
NASA NTRS1989

Radial diffusion and losses of energetic protons in the 5 to 12 R(S) region of Saturn's magnetosphere

Phase space density profiles for nearly equatorial mirroring energetic protons, derived from Pioneer 11 measurements during the 1979 Saturn encounter, were analyzed using solutions of the time-averaged radial diffusion equation in a dipolar planetary magnetic field. Then, several loss models were considered, ranging from a model that takes into account satellite absorption only to the one that accounts for satellite plus ring-E absorption plus added distributed losses, and the corresponding form of the time-averaged radial diffusion coefficient was determined for each case by a minimum-variance fit to the data-derived phase space densiy profiles. The model comparisons clearly indicated the need for distributed losses, in addition to satellite absorption, in order to explain the smooth decline of the experimental energetic proton phase space density profiles in the 5-12-Saturn radius region of Saturn's magnetosphere.

Hood, L. L.↗
DOE OSTI2023

Bifurcation-like transition of divertor conditions induced by X-point radiation in KSTAR L-mode plasmas *

Abstract Density ramps with ion grad B drift directed into lower single null KSTAR L-mode plasmas are associated with a simultaneous and abrupt reduction of the divertor particle flux on both low- and high-field-side targets when the mid-plane line averaged electron density reaches a given level. Target embedded Langmuir probe signals show a clear ‘cliff edge’ behavior similar to that observed in the divertor target electron temperature in DIII-D H-mode plasmas (Eldon et al 2017 Nucl. Fusion 57 066039; McLean et al 2015 J. Nucl. Mater. 463 533–6). The collapse of the particle flux is observed along the whole divertor target area (from private flux region to the far scrape-off layer (SOL)). The critical upstream density of this target flux cliff is invariant under fuel gas throughput modulation. The transition along the cliff occurs in tens of milliseconds. With the cliff, carbon impurities and deuterium neutrals transported through the X-point to the core produce a strong radiation spot near the X-point, seen on bolometric signals, and increase the upstream density. The experimental observations are consistent with time-dependent SOLPS-ITER simulations, which also demonstrate an abrupt transition of the target flux and upstream density with the increase in X-point radiation. The timescale of the cliff predicted by SOLPS-ITER is consistent with the experiment, although, it is influenced by gas throughput or time-dependent numerical methods. In the L-mode phase space of separatrix electron density and temperature, branches are divided based on target temperature, because the latter is strongly coupled to the radiation front and ionization front due to the monotonic characteristic of the parallel electron temperature distribution. Since the H-mode condition operates at a much higher upstream density and electron temperature in phase space, dissipation from sputtered carbon alone leads to the density limit before reaching the X-point radiation condition. This is therefore consistent with the fact that cliffs have never been observed in H-mode KSTAR experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical investigation of micro-galvanic corrosion in Mg alloys: Role of the cathodic intermetallic phase size and spatial distributions

Magnesium alloys are of increasing interest in structural applications due to their low-density, moderate specific strength and stiffness, recyclability, and high damping among other properties. However, the wide-scale applicability of magnesium alloys in structural applications has been limited due to many factors including its poor corrosion resistance. In this work, a numerical investigation to simulate the micro-galvanic corrosion behavior was performed to examine the influence of size and distribution of cathodic intermetallic phase (ß-Mg17Al12) in a Mg matrix. The ratio of cathodic to anodic surface area was kept constant in each simulation condition to understand the effect of size and spacing distributions. In general, fragmentation of a larger intermetallic particle into smaller ones was found to enhance the localized current density. However, the uniform distribution of this small intermetallic phase throughout the matrix was found to reduce the overall dissolution current density and hence, pitting corrosion severity.

Beura, Vikrant↗
DOE OSTI2025

Correlated X‐Ray and Electron Microscopies of a Single Biphasic GaAs Nanowire

Engineering the properties of semiconductors by changing their crystalline phase is a technologically and economically relevant alternative to doping using foreign elements, with strong potential for photonic and electronic applications. Although major advances have been reported recently for crystal‐phase engineering of III‐V and group IV semiconductor nanowires, interfacing two mismatched crystalline phases in a nanostructure induces several deformation mechanisms, which remain largely unexplored. Here, using state‐of‐the‐art synchrotron X‐ray nanobeam diffraction and transmission electron microscopy, subtle twisting and bending is unveiled within an individual GaAs nanowire containing cubic and hexagonal segments. Their role is discussed in accommodating the inter‐reticular spacing fluctuations, and their variations are correlated to the nanoscale phase distribution and to the effect of the NW support. This study brings direct evidence of a complex combination of deformation mechanisms in biphasic nanowires, which opens a new path to tune the nanowire properties with appealing perspectives for device engineering in nanophotonics and nanomechanics.

III-V semiconductors↗
DOE OSTI2023

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

FEgen (v.1): Field Emission Distribution Generator Freeware Based on Fowler–Nordheim Equation

This work presents an open-source code, FEgen, that generates initial field emission distributions based on the time-dependent Fowler–Nordheim equation. Current in-demand and open-source beam tracking software—such as ASTRA, IMPACT-T, and GPT—do not contain native distribution generator functions suitable for field emitted electron tracking and phase space analysis. Furthermore, the intuitive graphical user interface of FEgen generates input files shown to be fully compatible with both ASTRA and GPT. Here, detailed application examples demonstrating a subset of FEgen’s capabilities are given by creating a customized emitter pattern with application to transverse beam shaping and longitudinal beam dynamic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

Range Systems Simulation for the NASA Shuttle: Emphasis on Disaster and Prevention Management During Lift-Off

This article describes a decision-making system composed of a number of safety and environmental models for the launch phase of a NASA Space Shuttle mission. The components of this distributed simulation environment represent the different systems that must collaborate to establish the Expectation of Casualties (E(sub c)) caused by a failed Space Shuttle launch and subsequent explosion (accidental or instructed) of the spacecraft shortly after liftoff. This decision-making tool employs Space Shuttle reliability models, trajectory models, a blast model, weather dissemination systems, population models, amount and type of toxicants, gas dispersion models, human response functions to toxicants, and a geographical information system. Since one of the important features of this proposed simulation environment is to measure blast, toxic, and debris effects, the clear benefits is that it can help safety managers not only estimate the population at risk, but also to help plan evacuations, make sheltering decisions, establish the resources required to provide aid and comfort, and mitigate damages in case of a disaster.

Rabelo, Lisa↗
NASA NTRS2019

Understanding the Role of 𝛼 Particles in Oblique Heliospheric Shock Oscillations

Recent observations by DSCOVR provide high temporal resolution (50 samples per second) magnetic vector field data that allows investigating the details of oblique heliospheric shock oscillations. It was found that some of these shocks exhibit magnetic oscillations, both downstream and upstream of the shock front. The DSCOVR/MAG magnetic field data are supplemented by an extensive database of low Mach number (M < 3) low‐β (<1) shock data observed by Wind albeit with lower temporal resolution. Motivated by the observations, we use the 2.5D hybrid model of the oblique shocks with α particles in addition to kinetic protons and electron fluid. We model the properties of the oblique shocks for a number of typical parameters found in observations and study the effects of the shock parameters and the relative α particle abundances on the properties of the shock magnetic field, density, and velocity oscillations. We find the α particles “surf” on the shock front and produce a wake of density oscillations. We examine the details of the phase space of the ions as well as the ion velocity distribution functions in various parts of the shock and study their nonthermal properties. We determine the effects of the α particle kinetic properties and abundances on the structure and dynamics of the shock downstream oscillations for a range of parameters relevant to low Mach number low‐β heliospheric shocks.

Ofman, L.↗
DOE OSTI2021

BEAM COLLIMATION IN THE PIP-II LINAC TO BOOSTER TRANSFER LINE

The new PIP-II superconducting linac will deliver a 2 mA average H⁻ beam to the existing Booster synchrotron. The injected beam is accumulated by charge exchange over approximately 300 turns; phase space painting is used to mitigate space charge effects. To limit the power load on the internal waste beam absorber from the transverse tails of the H⁻ distribution missing the foil, the beam will be collimated in both planes in the linac to Booster transfer line using compact collimators of a novel design. Both the number of parasitic hits and the fraction of the beam missing the foil are sensitive functions of the H⁻ beam centroid position with respect to the edge of the foil. The positioning of the collimation is constrained by the availability of suitable space in the transfer line lattice, by specifics of the collimator design, by the phase space orientation at the collimator, and by the betatron phase advance to the foil needed to achieve proper orientation of the spatial distribution at the injection point. In this contribution, we describe the procedure by which collimator positions were optimized. We then discuss the expected performance of the overall system.

Johnson, David E.↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2021

The kinematic properties of Milky Way stellar halo populations

ABSTRACT In the Gaia era, stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern data sets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ‘action diamond’) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially anisotropic Gaia–Enceladus merger remnant. Assuming a 1:1 ratio of Gaia–Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 86 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic data sets can noticeably impact the E−Lz distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗