Collisional implicit particle-in-cell methods for fully kinetic simulations of dynamically compressing plasmas
Explore the source record for details and available documents.
Physics - Plasma physics↗
SEARCH · Engineering Papers
Results for “collisional plasma”
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.
Using Kinetic Entropy to Understand Dissipation in Re-connection and Turbulence
This report summarizes results obtained at Space Science Institute (SSI) under a collaborative project funded by the NSF-DOE partnership in Basic Plasma Physics. The overall goal of the project was to perform a study of kinetic entropy as a diagnostic of dissipation in several archetypal models, such as magnetic reconnection and decaying plasma turbulence. The specific task of the SSI group was to perform collisionless and weakly collisional Particle-in-Cell (PIC) simulations modeling these systems. Overall, the project has demonstrated that, as expected, kinetic entropy carries an imprint of weakly collisional processes commonly thought to be associated with dissipation and thus provides an important diagnostic of dissipation. Specific entropy-related metrics have been proposed that could be used as quantitative measures of dissipation. These metrics have been implemented as diagnostics in PIC simulations and their utility, as well as relation to other, previously proposed, measures of dissipation have been demonstrated. Further, PIC simulations were used to investigate a general question of how energy is dissipated in magnetized weakly collisional turbulence in both non-relativistic and relativistic plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
New model for the ion collection by cylindrical probes over a wide range of collisionality
Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Global Core Plasma Model
Over 40 years of ground and spacecraft plasmaspheric measurements have resulted in many statistical descriptions of plasmaspheric properties. In some cases, these properties have been represented as analytical descriptions that are valid for specific regions or conditions. For the most part, what has not been done is to extend regional empirical descriptions or models to the plasmasphere as a whole. In contrast, many related investigations depend on the use of representative plasmaspheric conditions throughout the inner magnetosphere. Wave propagation, involving the transport of energy through the magnetosphere, is strongly affected by thermal plasma density and its composition. Ring current collisional and wave particle losses also strongly depend on these quantities. Plasmaspheric also plays a secondary role in influencing radio signals from the Global Positioning System satellites. The Global Core Plasma Model (GCPM) is an attempt to assimilate previous empirical evidence and regional models for plasmaspheric density into a continuous, smooth model of thermal plasma density in the inner magnetosphere. In that spirit, the International Reference Ionosphere is currently used to complete the low altitude description of density and composition in the model. The models and measurements on which the GCPM is currently based and its relationship to IRI will be discussed.
Gallagher, Dennis L.↗
Empirical Modeling of the Plasmasphere
Over 40 years of ground and spacecraft plasmaspheric measurements have resulted in many statistical descriptions of plasmaspheric properties. In some cases, these properties have been represented as analytical descriptions that are valid for specific regions or conditions. For the most part, what has not been done is to extend regional empirical descriptions or models to the plasmasphere as a whole. In contrast, many related investigations depend on the use of representative plasmaspheric conditions throughout the inner magnetosphere. Wave propagation, involving the transport of energy through the magnetosphere, is strongly affected by thermal plasma density and its composition. Ring current collisional and wave particle losses also strongly depend on these quantities, Plasmaspheric also plays a secondary role in influencing radio signals from the Global Positioning System satellites. The Global Core Plasma Model (GCPM) is an attempt to assimilate previous empirical evidence and regional models for plasmaspheric density into a continuous, smooth model of then-nal plasma density in the inner magnetosphere. In that spirit, the International Reference Ionosphere is currently used to complete the low altitude description of density and composition in the model. The models and measurements on which the GCPM is currently based and its relationship to IRI will be discussed.
Gallagher, Dennis L.↗
Quasilinear modelling of collisional trapped electron modes
The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.
fusion plasma↗
Demonstrating the impact of ro-vibrationally excited H 2 on divertor detachment via population modelling
EDGE2D-EIRENE profiles characteristic for low recycling, high recycling and detached plasmas are postprocessed applying a vibrationally resolved collisional radiative (CR) model based on the Yacora solver to predict excited state densities and thus the molecular Lyman-, Werner- and Fulcher band emission intensities. These emission intensities are compared to predictions from the AMJUEL database standardly used in EIRENE. Lower Fulcher band and higher Lyman- and Werner band emission is predicted by the Yacora model in comparison to AMJUEL suggesting that inconsistent results from model based evaluation of measured spectra are expected depending on which model and molecular emission band is considered. The implications on detachment rate predictions are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of Shocks on Emission from Central Engines of Active Galactic Nuclei. I
In this paper we show that perturbations of the accretion flow within the central engines of some active galactic nuclei (AGNS) are likely to form shock waves in the accreting plasma. Such shocks, which may be either collisional or collisionless, can contribute to the observed high-energy temporal and spectral variability. Our rationale is the following: Observations show that the continuum emission probably originates in an optically thin, hot plasma in the AGN central engine. The flux and spectrum from this hot plasma varies significantly over light crossing timescales. Several authors have suggested that macroscopic perturbations contained within this plasma are the sources of this variability. In order to produce the observed emission the perturbations must be radiatively coupled with the optically thin hot matter and must also move with high velocities. We suggest that shocks, which can be very effective in randomizing the bulk motion of the perturbations, are responsible for this coupling. Shocks should form in the central engine, because the temperatures and magnetic fields are probably reduced below their virial values by radiative dissipation. Perturbations moving at Keplerian speeds, or strong non-linear excitations, result in supersonic and super-Alfvenic velocities leading to shock waves within the hot plasma. We show that even a perturbation smaller than the emitting region can form a shock that significantly modifies the continuum emission in an AGN, and that the spectral and temporal variability from such a shock generally resembles those of radio-quiet AGNS. As an example, the shock inducing perturbation in our model is a small main-sequence star, the capturing and eventual accretion of which are known to be a plausible process. We argue that shocks in the central engine may also provide a natural triggering mechanism for the "cold" component of Guilbert & Rees two-phase medium and an efficient mecha- nism for angular momentum transfer. Current and future missions, such as ASCA, XTE, XMM, AXAF, and ASTRO-E may determine the importance of shock-related emission from the central engines of AGNS.
Sivron, R.↗
The effects of laser pulse length and collisional ionization on the acceleration of titanium ions
Abstract The interaction of a relativistic laser pulse ( > 10 18 W cm −2 µ m 2 ) with foil targets can accelerate ions to energies of tens of MeV u −1 with optimized laser and target parameters. We report the results on simulations of the interaction of a 6.0 × 10 20 W cm −2 laser pulse incident on ultrathin (10–500 nm) titanium foils to investigate the roles of laser pulse duration and ionization mechanisms in the acceleration of titanium ions. While holding peak intensity constant, two laser pulse durations were investigated, 140 and 650 fs. The optimum thickness is dependent on pulse duration, as it requires the concurrence of target transparency with the incidence of the peak laser intensity. The collisional processes do not play a significant role in Ti ion beam generation from the 140 fs laser pulse duration at the optimum thickness (30 nm). However, for the 650 fs laser optimum, collisions improve the conversion efficiency of highly energetic ( > 10 MeV u −1 ), high charge titanium (Ti 21 − 22 + ) by a factor of 20, and the titanium ion cutoff energy by ∼15%. This improvement is due to the fact that collisional ionization increases the electron density of the plasma, which delays the time of relativistic transparency, causing collisions to decrease the optimum foil thickness from 150 to 100 nm. Additionally, collisional ionization increases the charge-to-mass ratio of the titanium, and injects more electrons into the accelerating sheath field. At the optimum thickness, target normal sheath acceleration is the dominant mechanism of acceleration, with additional contributions from radiation pressure and shock wave acceleration.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas
A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.
ETG↗
Electron temperature gradient driven transport model for tokamak plasmas
A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Thermal quench induced by a composite pellet-produced plasmoid
Injecting shattered pellets is the critical concept of the envisaged ITER disruption mitigation system (DMS). Rapid deposition of large amounts of material should presumably result in controlled cooling of the entire plasma. A considerable transfer of thermal energy from the electrons of the background plasma to the ions accompanies a localized material injection due to the ambipolar expansion along the magnetic field line of the cold and dense plasmoid produced by the ablated pellet. Radiation initially plays the dominant role in the energy balance of a composite plasmoid containing high-Z impurities. A competition between the ambipolar expansion and the radiative losses defines the Thermal Quench scenario, including the amount of pre-quench thermal energy radiated on a short collisional timescale—possibly detrimental for the plasma-facing components. The present work quantifies plasmoid energy balance for disruption mitigation parameters. For pure hydrogen injection, up to 90% of the pre-pellet electron thermal energy may go to the newly injected ions. We also demonstrate that a moderate high-Z impurity content within the plasmoid can reduce highly localized radiation at the beginning of the expansion. The thermal energy will then dissipate on the much longer ion collisional timescale, which would be attractive for ITER DMS.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Radiation from a high-temperature, low-density plasma - The X-ray spectrum of the solar corona
Solar coronal X ray spectrum calculation of high- temperature low-density plasma, considering line emission from electron collisional excitation and radiation
Koren, M.↗
X-ray line emission from the Tycho supernova remnant
The observation of the X-ray spectrum of the Tycho supernova remnant in the energy range 0.5 to 20 keV is discussed. Four significant line features in the spectrum: The K alpha lines of silicon, sulphur, and iron; and the L lines of iron are examined. Comparisons between the silicon and sulphur equivalent widths and K alpha iron line energies of Tycho and Cas A are discussed. Suggest that the X-ray emitting plasma in Tycho is further from collisional ionization equilibrium than that of Cas A.
Pravdo, S. H.↗
X-ray line emission from the Tycho supernova remnant
The A-2 experiment onboard HEAO 1 has observed the X-ray spectrum of the Tycho supernova remnant in the energy range 0.5-20 keV. The spectrum reveals four significant line features: the K(alpha) lines of silicon, sulfur, and iron; and the L lines of iron. Comparisons between the silicon and sulfur equivalent widths and K(alpha) iron line energies of Tycho and Cas A suggest that the X-ray emitting plasma in Tycho is further from collisional ionization equilibrium than that in Cas A.
Pravdo, S. H.↗
On the dynamo generation of flux ropes in the Venus ionosphere
Small scale magnetic field structures or 'flux ropes' observed in the ionosphere of Venus can be interpreted as the result of a kinematic dynamo process acting on weak seed fields. The seed fields result from the prevailing downward convection of magnetic flux from the vicinity of the ionopause, while small scale fluctuations in the velocity of the ionospheric plasma, which can be caused by collisional coupling to gravity waves in the neutral atmosphere, provide the mechanism by which the field is twisted and redistributed into features of similar scale. This mechanism naturally explains some of the average properties of flux ropes such as the variation of their characteristics with altitude and solar zenith angle. It also elucidates the relationship between the large scale and small scale ionospheric magnetic fields.
Luhmann, J. G.↗
Dielectronic recombination as a function of electric field strength
Dielectronic recombination (DR) is the dominant recombination mechanism at coronal temperatures and densities. We present a procedure for calculating DR rate coefficients as a function of electric field strength and apply this method to carbon ions. We focus on the competing effects of enhancement by plasma microfields and rate decrease through collisional excitation and ionization. We find that, in the case of C(3+), a significant rate enhancement results, leading to a reinterpretation of C IV emission-line intensities in the sun and late-type stars. We further consider how macroscopic electric fields, in particular motional electric fields, can affect DR rate coefficients, demonstrating dramatic rate increases for a number of the carbon ions.
Reisenfeld, Daniel B.↗
Understanding JET-C quiescent phases with edge harmonic magnetohydrodynamic activity and comparison with behaviour under ITER-like wall conditioning
An analysis of edge localised mode-free (quiescent) H-mode discharges exhibiting edge harmonic magnetoydrodynamic activity in the JET-carbon wall machine is presented. It is observed that the otherwise quiescent pulses with multiple-n harmonic oscillations are sustained until a threshold in pedestal electron density and collisionality is crossed. The macroscopic pedestal parameters associated with the quiescent phase are compared with those of a database of JET-ELMy discharges with both carbon and ITER-like wall (ILW). This comparison provides the identification of the existence regions in the relevant pedestal and global plasma parameters for edge harmonic oscillations (EHOs) in JET plasmas. Although the ELMy database scans pedestal collisionality and β values typical of ET-carbon quiescent operation, shaping and current are not simultaneously compatible with EHO existence. Nevertheless, ILW operation with JET-carbon quiescent-like parameters could in principle be achieved, and improved pedestal performance could be observed in more recent JET-ILW pulses.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗