Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distribution functions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Construction of wave dark matter halos: Numerical algorithm and analytical constraints

Here we present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos—such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes and show how it simplifies in the high-energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrödinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗
DOE OSTI2022

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of a drive cylinder for low collisional experiments on magnetic reconnection

For most laboratory plasma experiments, Coulomb collisions between the particle species are sufficiently frequent that the particle distribution functions are relaxed to a near-Maxwellian form. This hampers the applicability of such experiments to phenomena observed in tenuous and near-collisionless space plasma. The Terrestrial Reconnection EXperiment (TREX) at the Wisconsin Plasma Physics Laboratory aims to study collisionless reconnection for parameters relevant to the Earth’s magnetosphere. To reduce the role of collisional effects, a reconnection Drive Cylinder has been developed, which increases both the effective system size of the TREX configuration and the rate at which reconnection can be driven. These two effects now permit TREX to reach a kinetic reconnection regime where collisional effects are minimized. The Drive Cylinder is comprised of 12 single loop drive-coils connected in parallel to a 10 kV capacitor bank. Insulated sheets of aluminum are applied to smooth the magnetic fields and enhance the drive efficiency. Following is a description of the technical details and performance of the Drive Cylinder.

Instruments & Instrumentation↗
DOE OSTI2024

Lower hybrid drift instability in nonthermal plasmas

Lower hybrid drift instability (LHDI) is driven by the cross-field current and operates in the vicinity of the lower-hybrid frequency, between the ion- and electron-gyro frequencies, and with wavelengths between the electron and ion thermal gyro radii. The free energy source that drives this instability resides in the density gradient associated with an inhomogeneous plasma. The existing literature on LHDI assumes that the charged particle distribution function is given by a Maxwellian form, but the space plasma is pervasively observed to feature nonthermal characteristics. This paper extends the theory of LHDI to nonthermal plasmas. The generalized theory of LHDI is, thus, applicable to various space plasma environments characterized by nonthermal plasma velocity distribution functions.

Arya, Neetasha (ORCID:0009000236075605)↗
DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hydrodynamic generators in relativistic kinetic theory

Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Simulating Boundary Region Plasma Instabilities

The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Turbulence and Particle Acceleration in a Relativistic Plasma

Abstract In a collisionless plasma, the energy distribution function of plasma particles can be strongly affected by turbulence. In particular, it can develop a nonthermal power-law tail at high energies. We argue that turbulence with initially relativistically strong magnetic perturbations (magnetization parameter σ ≫ 1) quickly evolves into a state with ultrarelativistic plasma temperature but mildly relativistic turbulent fluctuations. We present a phenomenological and numerical study suggesting that in this case, the exponent α in the power-law particle-energy distribution function, f ( γ ) d γ ∝ γ − α d γ , depends on magnetic compressibility of turbulence. Our analytic prediction for the scaling exponent α is in good agreement with the numerical results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Galactic Bar Resonances with Diffusion: An Analytic Model with Implications for Bar–Dark Matter Halo Dynamical Friction

Abstract The secular evolution of disk galaxies is largely driven by resonances between the orbits of “particles” (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari, and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies—effects that have, with few exceptions, been previously explored only with complex N -body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles’ slow actions. We solve this equation for various values of the dimensionless diffusion strength Δ, and then apply our theory to the calculation of bar–halo dynamical friction. For Δ = 0, we recover the classic result of Tremaine and Weinberg that friction ultimately vanishes, owing to the phase mixing of resonant orbits. However, for Δ > 0, we find that diffusion suppresses phase mixing, leading to a finite torque. Our results suggest that stochasticity—be it physical or numerical—tends to increase bar–halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time.

79 ASTRONOMY AND ASTROPHYSICS↗