SEARCH · Search NASA

Results for “Collisionless relaxation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Formation mechanisms of laboratory double layers

The evolution processes of double layers were studied in a series of laboratory experiments. It was found that the existence of virtual cathode-type potential wells at the electron injection boundary was the dominant triggering mechanism. The rapid growth of the potential wells led to collisionless ion trapping and the establishment of the necessary trapped ion population. For double layers with small potential drops, collisionless ion trapping actually induced ion-ion streaming instabilities and the formation of ion phase-space vortices. In this regime, the system often exhibited relaxation-type oscillations which correspond to the disruption and the recovery of the double layers.

Chan, Chung↗

DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

On the collisional damping of plasma velocity space instabilities

For plasma velocity space instabilities driven by particle distributions significantly deviated from a Maxwellian, weak collisions can damp the instabilities by an amount that is significantly beyond the collisional rate itself. This is attributed to the dual role of collisions that tend to relax the plasma distribution toward a Maxwellian and to suppress the linearly perturbed distribution function. The former effect can dominate in cases where the unstable non-Maxwellian distribution is driven by collisionless transport on a timescale much shorter than that of collisions, and the growth rate of the ideal instability has a sensitive dependence on the distribution function. The whistler instability driven by electrostatically trapped electrons is used as an example to elucidate such a strong collisional damping effect of plasma velocity space instabilities, which is confirmed by first-principles kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

Linear and non-linear studies of Alfven waves in space. Stationary and dynamic processes in magnetospheric plasmas

The Final Technical Report on linear and non-linear studies of Alfven waves in space is presented. Areas of research included relaxation of magnetotail plasmas with field-aligned currents; the equilibrium dayside magnetosphere; macroscale particle simulation of kinetic Alfven wave physics; ballooning stability of plasmas with sheared equilibrium flows; theory of the drift-mirror instability; collisionless tearing instability in magnetotail plasmas; and nonadiabatic behavior of the magnetic moment of a charged particle in a dipole magnetic field and the development of stochastic webs.

Bhattacharjee, A.↗

DOE OSTI2023

More relaxed intracluster gas than galaxies in clusters in quasi-equilibrium

ABSTRACT During cluster mergers, the intracluster gas and member galaxies undergo dynamic evolution, but at different time-scales and reach different states. We collect 24 galaxy clusters in quasi-equilibrium state as indicated by the X-ray image, and calculate the cluster orientations and three kinds of dynamical parameters, i.e. the normalized centroid offset, the sphere index, and the ellipticity, for these clusters from the distributions of member galaxies and also the intracluster gas. We find consistent alignments for the orientations estimated from the two components. However, the three kinds of dynamical parameters indicated by member galaxies are systematically larger than those derived from the gas component, suggesting that the gas component is more relaxed than member galaxies. Differences of dynamical features between the intracluster gas and member galaxies are independent of cluster mass and concentration. We conclude that the intracluster gas reaches the dynamic equilibrium state earlier than the almost-collisionless member galaxies.

Astronomy & Astrophysics↗

NASA NTRS1984

Temperature relaxation in supernova remnants, revisited

Some supernova remnants are expanding into a partially neutral medium. The neutral atoms which are engulfed by the fast blast shock are collisionally ionized to eject low-energy secondary electrons. Calculations are conducted of the temperature relaxation through Coulomb collisions among the secondary electrons, the shocked electrons, and the ions, assuming that the three species have Maxwellian velocity distributions. The results are applied to a self-similar blast wave. If the efficiency of collisionless electron heating at the shock front is high in young remnants such as Tycho, the secondary electrons may be much cooler than both the shocked electrons and the ions. In this case, the emergent X-ray continuum spectrum will have a two-temperature, or a power-law, appearance. This effect may have been observed in the bright rim of the remnant of SN 1006.

Itoh, H.↗

NASA NTRS1995

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗

DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗

DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2001

A Semianalytical Ion Current Model for Radio Frequency Driven Collisionless Sheaths

We propose a semianalytical ion dynamics model for a collisionless radio frequency biased sheath. The model uses bulk plasma conditions and electrode boundary condition to predict ion impact energy distribution and electrical properties of the sheath. The proposed model accounts for ion inertia and ion current modulation at bias frequencies that are of the same order of magnitude as the ion plasma frequency. A relaxation equation for ion current oscillations is derived which is coupled with a damped potential equation in order to model ion inertia effects. We find that inclusion of ion current modulation in the sheath model shows marked improvements in the predictions of sheath electrical properties and ion energy distribution function.

Bose, Deepak↗

NASA NTRS1979

Collisionless galaxy simulations

Three-dimensional fully self-consistent computer models were used to determine the evolution of galaxies consisting of 100 000 simulation stars. Comparison of two-dimensional simulations with three-dimensional simulations showed only a very slight stabilizing effect due to the additional degree of freedom. The addition of a fully self-consistent, nonrotating, exponential core/halo component resulted in considerable stabilization. A second series of computer experiments was performed to determine the collapse and relaxation of initially spherical, uniform density and uniform velocity dispersion stellar systems. The evolution of the system was followed for various amounts of angular momentum in solid body rotation. For initally low values of the angular momentum satisfying the Ostriker-Peebles stability criterion, the systems quickly relax to an axisymmetric shape and resemble elliptical galaxies in appearance. For larger values of the initial angular momentum bars develop and the systems undergo a much more drastic evolution.

Hohl, F.↗

DOE OSTI2020

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Examining the self-interaction of dark matter through central cluster galaxy offsets

ABSTRACT While collisionless cold dark matter models have been largely successful in explaining a wide range of observational data, some tensions still exist, and it remains possible that dark matter possesses a non-negligible level of self-interactions. In this paper, we investigate a possible observable consequence of self-interacting dark matter: offsets between the central galaxy and the centre of mass of its parent halo. We examine 23 relaxed galaxy clusters in a redshift range of 0.1–0.3 drawn from clusters in the Dark Energy Survey and the Sloan Digital Sky Survey which have archival Chandra X-ray data of sufficient depth for centre and relaxation determination. We find that most clusters in our sample show non-zero offsets between the X-ray centre, taken to be the centroid within the cluster core, and the central galaxy position. All of the measured offsets are larger, typically by an order of magnitude, than the uncertainty in the X-ray position due to Poisson noise. In all but six clusters, the measured offsets are also larger than the estimated, combined astrometric uncertainties in the X-ray and optical positions. A more conservative cut on concentration to select relaxed clusters marginally reduces but does not eliminate the observed offset. With our more conservative sample, we find an estimated median X-ray to central galaxy offset of $\mu = 6.0 ^{+ 1.4}_{- 1.5}$ kpc. Comparing to recent simulations, this distribution of offsets is consistent with some level of dark matter self-interaction, though further simulation work is needed to place constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1988

Temporal evolution of collisionless sheaths

Experimental measurements have been performed to determine the temporal evolution of Langmuir sheaths near an electrode to which a negative step bias is applied in a collisionless argon plasma. The plasma was produced by a hot-filament discharge in a multidipole device. Plasma potential data were obtained using emissive probes with two different techniques: time-resolved sampling and time-averaged techniques. The sheath is found to initially form close to the electrode, to extend to a maximum separation, and to contract to a steady-state value. The time scale required to reach a steady state is close to the time scale of the presheath relaxation. Characteristics of sheaths in RF plasmas are also measured using a parallel-plate plasma capacitor. It is observed that the plasma potential profile has significant variation with frequency, even for frequencies as low as 1 kHz which are far below the ion plasma frequency (about 1 MHz).

Cho, M. H.↗

NASA NTRS1972

Area influences and floating potentials in Langmuir probe measurements.

Analysis of the influence of a relatively small reference electrode on a Langmuir probe measurement of plasma density and temperature. The analysis employs the collisionless theory of Laframboise for spherical and cylindrical probes and in a completely general way includes a broad range of electron densities and temperatures, as well as all possible values of ion mass and multiplicity of ionization. It is shown that a ratio alpha of reference-electrode area to probe area of 10,000 will guarantee no distortion of the measurement as a result of a shifting reference-electrode potential. It is further shown that the constraint on alpha can be relaxed by approximately two orders of magnitude when the charge-normalized ion mass is decreased from 200 to 1 amu and the ratio of reference-electrode radius to Debye length is decreased from 100 to 0. An additional result of the analysis is the dependence of a probe's floating potential on its geometry and radius, as well as on the properties of the plasma. A discussion of the latter result reveals in a striking manner the significant differences between the ion-current responses of spherical and cylindrical probes.

Szuszczewicz, E. P.↗

NASA NTRS1991

Transport of gyration-dominated space plasmas of thermal origin. I - Generalized transport equations

Grad's 20-moment set of transport equations has been examined in the limit of strong external magnetic fields. This simplified set of equations describes the transport of mass and parallel momentum as well as the transport of parallel and perpendicular energy and heat flow in the magnetic field direction. The effect of collisions was calculated assuming a modified relaxation model. Wave speeds and normal modes of the simplified set of equations were examined for an ion and electron gas. It was found that four of the ten normal modes are electron thermal-heat waves which approximately decouple from the six ion waves in the system. When low-frequency waves are considered (slow-wave approximation), this allows the electron energy and heat flow equations to be solved separately from the ion equations and in a time-independent fashion. When this was done, it was found that, under certain conditions, these equations predict an infinite electron perpendicular temperature, Te-perpendicular, in the collisionless regime. This occurs whenever Te-perpendicular is greater than the parallel temperature, at any point along collisionless and diverging magnetic field lines. When applied appropriately, the significantly simplified set of equations derived here are well suited for application to a large variety of problems in planetary ionospheres and magnetospheres.

Gombosi, Tamas I.↗

DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗