Search NASA⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Collisionless Relaxation of the Ion Ring Distribution in Space Plasma

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths' ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

Ofman, Leon↗
NASA NTRS1974

Collisionless relaxation in spiral galaxy models

The increase in random kinetic energy of stars by rapidly fluctuating gravitational fields (collisionless or violent relaxation) in disk galaxy models is investigated for three interaction potentials of the stars corresponding to (1) point stars, (2) rod stars of length 2 kpc, and (3) uniform density spherical stars of radius 2 kpc. To stabilize the galaxy against the large scale bar forming instability, a fixed field corresponding to a central core or halo component of stars was added with the stars containing at most 20 percent of the total mass of the galaxy. Considerable heating occurred for both the point stars and the rod stars, whereas the use of spherical stars resulted in a very low heating rate. The use of spherical stars with the resulting low heating rate will be desirable for the study of large scale galactic stability or density wave propagation, since collective heating effects will no longer mask the phenomena under study.

Hohl, F.↗
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↗
NASA NTRS2014

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗
NASA NTRS1977

Wind enhanced planetary escape - Collisional modifications

Effects of collisions and finite winds characteristic of a highly perturbed atmosphere on the thermal escape of terrestrial hydrogen and helium are investigated using a Monte Carlo approach. The limiting cases of vertical and horizontal winds are considered, and the relaxation layer between the collisionless exosphere and the collision-dominated thermosphere is modeled as a plane-parallel slab of given column density, depth, and atmospheric density. For both gases, the upwardly injected flux at the base of the relaxation layer is compared with the returning downward flux distribution at the same location; the technique is also applied to the atmosphere of Titan. The results show that inclusion of collisions in the escape model for terrestrial hydrogen with winds effectively throttles the escape process, that collisional throttling is negligible for helium when the exobase temperature is at least 5000 K, and that the escape of a planetary-atmosphere constituent depends on the ratio of its gravitational and kinetic energies as well as on the ratio of its mass to that of the background gas.

Curtis, S. A.↗
NASA NTRS1980

Interstellar shock waves

The structure of interstellar shocks driven by supernova remnants and by expanding H II regions around early-type stars is discussed. Jump conditions are examined, along with shock fronts, post-shock relaxation layers, collisional shocks, collisionless shocks, nonradiative shocks, radiative atomic shocks, and shock models of observed nebulae. Effects of shock waves on interstellar molecules are examined, with reference to the chemistry behind shock fronts, infrared and vibrational-rotational cooling by molecules, and observations of shocked molecules. Some current problems and applications of the study of interstellar shocks are summarized, including the initiation of star formation by radiative shock waves, interstellar masers, the stability of shocks, particle acceleration in shocks, and shocks in galactic nuclei.

Mckee, C. F.↗
NASA NTRS1981

On a theory of temporal fluctuations in the electrostatic potential structures associated with auroral arcs

A possible mechanism is presented for the generation of large-amplitude temporal fluctuations in the structure of the electron energization region associated with auroral arcs. The mechanism is based on the observation that the auroral arc system resembles a laboratory circuit consisting of the series connection of battery, resistance and a forward biased diode containing collisionless plasma in which large-amplitude relaxation oscillations are sometimes observed to be superimposed on the steady-state current. It is shown that in both the laboratory and auroral systems, in which a localized auroral arc dynamo, the ionosphere and the electron energization region are involved, the oscillations are controlled by the times for ions and electrons to traverse the acceleration region, which also characterize the low- and high-frequency structure of the fluctuating waveform. It is demonstrated that a plausible one-dimensional double-layer model of the auroral arc acceleration region exhibits the dynamic negative resistance necessary for the generation of oscillations by the present mechanism. Finally, consideration is given to two kinds of auroral phenomena which might be associated with the mechanism: the 10-Hz quasi-periodic flickering aurora and 10-Hz modulations in the intensity of electrostatic hydrogen cyclotron waves.

Silevitch, M. B.↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗
NASA NTRS1995

Heat flux and viscosity of ions in the collisionless solar wind

Between 1 and 2 solar radii, the Coulomb-collision mean free path for thermal ions exceeds the scale height of the solar atmosphere. The expanding solar plasma becomes collisionless and the kinetics of the solar wind are no longer dominated by thermalizing collisions. The usual Braginskii-type expressions for solar wind ion heat flux and viscosity are no longer valid. However, another microscale still exists in the solar wind, dictated by the gyro-radius of ions in the turbulent embedded solar wind magnetic field. Wave-particle interactions will act to isotropize (but not thermalize) particle distributions, and the relevant microscale for this process is the ion gyro-radius. The ion distribution can be modelled as undergoing isotropizing 'collisions,' with the relevant mean free path scaling with gyro-radius. Here, the author presents the heat flux and viscosity expected for solar wind protons which are relaxing to isotropy on a microscale that scales with gyro-radius. The collisionless viscosity and heat flux have a functional dependence different than their collisional analogs. The collisional expressions for ion viscosity and heat flux drastically overestimate the efficiency of diffusive energy and momentum transport actually operative in the solar wind.

Williams, L. L.↗
NASA NTRS1989

Energy loss of fast nonthermal electrons in plasmas

Simple analytical expressions for equilibration times of nonrelativistic monoenergetic electrons in plasmas are evaluated in the 'weak'-beam approximation when the density of the monoenergetic electrons is much smaller than the plasma density. The equilibration time is defined as the time needed by the beam of monenergetic electrons to lose most of its energy as a result of collisions with plasma particles having a Maxwellian energy distribution. The process of energy equilibration is treated as a statistical superposition of both elastic and inelastic collisions in the plasma. The possibility of collisionless equilibration is discussed. Comparison of the equilibration times with the Spitzer relaxation times indicates that the former times are more appropriate for an estimate of the energy loss of the 'weak' electron beams in highly ionized plasmas.

Kunc, Joseph A.↗
NASA NTRS1995

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗
NASA NTRS1975

N-body simulations of disks

The methods used in large-scale n-body simulations are discussed. However, the present review concentrates on results already obtained in n-body simulations using systems containing up to 200,000 simulation stars. Results are presented which show that the stability criterion developed for flattened systems applies only to truly axisymmetric instabilities. Purely stellar disks acquire rather large velocity dispersions, generally two or more times the velocity dispersion required by Toomre (1964) for axisymmetric stability. In computer simulations, the bar-forming instability can be prevented only by comparatively large velocity dispersions. However, simulations including the effects of the galactic halo and core as a fixed background field show that bar formation can be prevented for fixed halo components as large or larger than the self-consistent disk component. Experiments performed to determine the collisional relaxation time for large-scale gravitational n-body calculations show that these models are indeed 'collisionless'.

Hohl, F.↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
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↗
NASA NTRS1993

Discreteness noise versus force errors in N-body simulations

A low accuracy in the force calculation per time step of a few percent for each particle pair is sufficient for collisionless N-body simulations. Higher accuracy is made meaningless by the dominant discreteness noise in the form of two-body relaxation, which can be reduced only by increasing the number of particles. Since an N-body simulation is a Monte Carlo procedure in which each particle-particle force is essentially random, i.e., carries an error of about 1000 percent, the only requirement is a systematic averaging-out of these intrinsic errors. We illustrate these assertions with two specific examples in which individual pairwise forces are deliberately allowed to carry significant errors: tree-codes on supercomputers and algorithms on special-purpose machines with low-precision hardware.

Hernquist, Lars↗
NASA NTRS1973

Relaxation time in disk galaxy simulations.

An initially stationary and stable axisymmetric disk of stars with a mass spectrum of stars is used to determine the collisional relaxation time. The relaxation time as determined from the rate of energy equipartition was found to be 560 rotation periods for stars with 10 times the mean star mass and was 1700 rotation periods for stars with 0.55 times the mean star mass. These times are in general agreement with theoretical predictions for the relaxation time of the two mass groups. The results show that the model used for the large-scale gravitational N-body calculations is indeed 'collisionless.'

Hohl, F.↗