Consistent description of kinetics and hydrodynamics of dusty plasma
Markiv, B.; Tokarchuk, M.; National University “Lviv Polytechnic,” 12 Bandera St., 79013 Lviv
2014-02-15
A consistent statistical description of kinetics and hydrodynamics of dusty plasma is proposed based on the Zubarev nonequilibrium statistical operator method. For the case of partial dynamics, the nonequilibrium statistical operator and the generalized transport equations for a consistent description of kinetics of dust particles and hydrodynamics of electrons, ions, and neutral atoms are obtained. In the approximation of weakly nonequilibrium process, a spectrum of collective excitations of dusty plasma is investigated in the hydrodynamic limit.
Generalized Kinetic Description of Steady-State Collisionless Plasmas
NASA Technical Reports Server (NTRS)
Khazanov, G. V.; Liemohn, M. W.; Krivorutsky, E. N.
1997-01-01
We present a general solution to the collisionless Boltzmann (Vlasov) equation for a free-flowing plasma along a magnetic field line using Liouville's theorem, allowing for an arbitrary potential structure including non-monotonicities. The constraints of the existing collisionless kinetic transport models are explored, and the need for a more general approach to the problem of self- consistent potential energy calculations is described. Then a technique that handles an arbitrary potential energy distribution along the field line is presented and discussed. For precipitation of magnetospherically trapped hot plasma, this model yields moment calculations that vary by up to a factor of two for various potential energy structures with the same total potential drop. The differences are much greater for the high-latitude outflow scenario, giving order of magnitude variations depending on the shape of the potential energy distribution.
Cremaschini, Claudio; Tessarotto, Massimo
2011-11-15
A largely unsolved theoretical issue in controlled fusion research is the consistent kinetic treatment of slowly-time varying plasma states occurring in collisionless and magnetized axisymmetric plasmas. The phenomenology may include finite pressure anisotropies as well as strong toroidal and poloidal differential rotation, characteristic of Tokamak plasmas. Despite the fact that physical phenomena occurring in fusion plasmas depend fundamentally on the microscopic particle phase-space dynamics, their consistent kinetic treatment remains still essentially unchallenged to date. The goal of this paper is to address the problem within the framework of Vlasov-Maxwell description. The gyrokinetic treatment of charged particles dynamics is adopted for the construction of asymptotic solutions for the quasi-stationary species kinetic distribution functions. These are expressed in terms of the particle exact and adiabatic invariants. The theory relies on a perturbative approach, which permits to construct asymptotic analytical solutions of the Vlasov-Maxwell system. In this way, both diamagnetic and energy corrections are included consistently into the theory. In particular, by imposing suitable kinetic constraints, the existence of generalized bi-Maxwellian asymptotic kinetic equilibria is pointed out. The theory applies for toroidal rotation velocity of the order of the ion thermal speed. These solutions satisfy identically also the constraints imposed by the Maxwell equations, i.e., quasi-neutrality and Ampere's law. As a result, it is shown that, in the presence of nonuniform fluid and EM fields, these kinetic equilibria can sustain simultaneously toroidal differential rotation, quasi-stationary finite poloidal flows and temperature anisotropy.
Cremaschini, Claudio; Miller, John C.; Tessarotto, Massimo
2011-06-15
A kinetic treatment is developed for collisionless magnetized plasmas occurring in high-temperature, low-density astrophysical accretion disks, such as are thought to be present in some radiatively inefficient accretion flows onto black holes. Quasi-stationary configurations are investigated, within the framework of a Vlasov-Maxwell description. The plasma is taken to be axisymmetric and subject to the action of slowly time-varying gravitational and electromagnetic fields. The magnetic field is assumed to be characterized by a family of locally nested but open magnetic surfaces. The slow collisionless dynamics of these plasmas is investigated, yielding a reduced gyrokinetic Vlasov equation for the kinetic distribution function. For doing this, an asymptotic quasi-stationary solution is first determined, represented by a generalized bi-Maxwellian distribution expressed in terms of the relevant adiabatic invariants. The existence of the solution is shown to depend on having suitable kinetic constraints and conditions leading to particle trapping phenomena. With this solution, one can treat temperature anisotropy, toroidal and poloidal flow velocities, and finite Larmor-radius effects. An asymptotic expansion for the distribution function permits analytic evaluation of all the relevant fluid fields. Basic theoretical features of the solution and their astrophysical implications are discussed. As an application, the possibility of describing the dynamics of slowly time-varying accretion flows and the self-generation of magnetic field by means of a ''kinetic dynamo effect'' are discussed. Both effects are shown to be related to intrinsically kinetic physical mechanisms.
Ramos, J. J.
2010-08-15
A closed theoretical model to describe slow, macroscopic plasma processes in a fusion-relevant collisionality regime is set forward. This formulation is a hybrid one, with fluid conservation equations for particle number, momentum and energy, and drift-kinetic closures. Intended for realistic application to the core of a high-temperature tokamak plasma, the proposed approach is unconventional in that the ion collisionality is ordered lower than in the ion banana regime of neoclassical theory. The present first part of a two-article series concerns the electron system, which is still equivalent to one based on neoclassical electron banana orderings. This system is derived such that it ensures the precise compatibility among the complementary fluid and drift-kinetic equations, and the rigorous treatment of the electric field and the Fokker-Planck-Landau collision operators. As an illustrative application, the special limit of an axisymmetric equilibrium is worked out in detail.
NASA Astrophysics Data System (ADS)
Molinari, V. G.; Rocchi, F.; Sumini, M.
2002-01-01
Aim of this work is to extend the results obtained in a previous study on the magnetic confinement and stability of a quantum degenerate non-neutral fermion plasma. This extension consists in the inclusion in the previously set up model of the effects of the exchange forces, and generalises the Thomas-Fermi (TF) approach used in the referenced work towards a Thomas-Fermi-Dirac (TFD) statistical description. The TF model has not only been used extensively and with success in these years to study atomic, nuclear and molecular properties, or to evaluate features of matter in extreme conditions such as low temperatures and/or high densities typical of astrophysics and inertial confinement fusion experiments, but also to found hydrodynamic theories for the diffusion and stability of fermion plasmas, one component non-neutral degenerate fluids, plasmas etc. In this paper an equation for density profiles in cylindrical symmetry is found, from the semiclassical kinetic theory of quantum gases, which takes into account the effects of temperature, average velocity, external magnetic field and quantum exchange. Numerical solutions of this equation for the case of complete quantum degeneracy are given and comparisons with the previous results are carried out.
NASA Astrophysics Data System (ADS)
Gubchenko, V. M.
2015-12-01
In part I of the work, the physical effects responsible for the formation of low-speed flows in plasma coronas, coupled with formation of coronas magnetosphere-like structures, are described qualitatively. Coronal domain structures form if we neglect scales of spatial plasma dispersion: high-speed flows are accumulated in magnetic tubes of the open domains, while magnetic structures and low-speed flows are concentrated within boundaries of domains. The inductive electromagnetic process occurring in flows of the hot collisionless plasma is shown to underlie the formation of magnetosphere-like structures. Depending on the form of the velocity distribution function of particles (PDF), a hot flow differently reveals its electromagnetic properties, which are expressed by the induction of resistive and diamagnetic scales of spatial dispersion. These determine the magnetic structure scales and structure reconstruction. The inductive electromagnetic process located in lines of the plasma nontransparency and absorption, in which the structures of excited fields are spatially aperiodic and skinned to the magnetic field sources. The toroidal and dipole magnetic sources of different configurations are considered for describing the corona structures during the solar maximum and solar minimum.
Photon kinetic modeling of laser pulse propagation in underdense plasma
Reitsma, A. J. W.; Trines, R. M. G. M.; Bingham, R.; Cairns, R. A.; Mendonca, J. T.; Jaroszynski, D. A.
2006-11-15
This paper discusses photon kinetic theory, which is a description of the electromagnetic field in terms of classical particles in coordinate and wave number phase space. Photon kinetic theory is applied to the interaction of laser pulses with underdense plasma and the transfer of energy and momentum between the laser pulse and the plasma is described in photon kinetic terms. A comparison is made between a one-dimensional full wave and a photon kinetic code for the same laser and plasma parameters. This shows that the photon kinetic simulations accurately reproduce the pulse envelope evolution for photon frequencies down to the plasma frequency.
Kinetic Description of the Impedance Probe
NASA Astrophysics Data System (ADS)
Oberrath, Jens; Lapke, Martin; Mussenbrock, Thomas; Brinkmann, Ralf
2011-10-01
Active plasma resonance spectroscopy is a well known diagnostic method. Many concepts of this method are theoretically investigated and realized as a diagnostic tool, one of which is the impedance probe (IP). The application of such a probe in plasmas with pressures of a few Pa raises the question whether kinetic effects have to be taken into account or not. To address this question a kinetic model is necessary. A general kinetic model for an electrostatic concept of active plasma spectroscopy was presented by R.P. Brinkmann and can be used to describe the multipole resonance probe (MRP). In principle the IP is interpretable as a special case of the MRP in lower order. Thus, we are able to describe the IP by the kinetic model of the MRP. Based on this model we derive a solution to investigate the influence of kinetic effects to the resonance behavior of the IP. Active plasma resonance spectroscopy is a well known diagnostic method. Many concepts of this method are theoretically investigated and realized as a diagnostic tool, one of which is the impedance probe (IP). The application of such a probe in plasmas with pressures of a few Pa raises the question whether kinetic effects have to be taken into account or not. To address this question a kinetic model is necessary. A general kinetic model for an electrostatic concept of active plasma spectroscopy was presented by R.P. Brinkmann and can be used to describe the multipole resonance probe (MRP). In principle the IP is interpretable as a special case of the MRP in lower order. Thus, we are able to describe the IP by the kinetic model of the MRP. Based on this model we derive a solution to investigate the influence of kinetic effects to the resonance behavior of the IP. The authors acknowledge the support by the Deutsche Forschungsgemeinschaft (DFG) via the Ruhr University Research School and the Federal Ministry of Education and Research in frame of the PluTO project.
Population kinetics in dense plasmas
Schlanges, M.; Bornath, T.; Prenzel, R.; Kremp, D.
1996-07-01
Starting from quantum kinetic equations, rate equations for the number densities of the different atomic states and equations for the energy density are derived which are valid for dense nonideal plasmas. Statistical expressions are presented for the rate coefficients taking into account many-body effects as dynamical screening, lowering of the ionization energy and Pauli-blocking. Based on these generalized expressions, the coefficients of impact ionization, three-body recombination, excitation and deexcitation are calculated for nonideal hydrogen and carbon plasmas. As a result, higher ionization and recombination rates are obtained in the dense plasma region. The influence of the many-body effects on the population kinetics, including density and temperature relaxation, is shown then for a dense hydrogen plasma. {copyright} {ital 1996 American Institute of Physics.}
Lagrangian description of warm plasmas
NASA Technical Reports Server (NTRS)
Kim, H.
1970-01-01
Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.
Kinetic Modeling of Divertor Plasma
NASA Astrophysics Data System (ADS)
Ishiguro, Seiji; Hasegawa, Hiroki; Pianpanit, Theerasarn
2015-11-01
Particle-in-Cell (PIC) simulation with the Monte Carlo collisions and the cumulative scattering angle coulomb collision can present kinetic dynamics of divertor plasmas. We are developing two types of PIC codes. The first one is the three dimensional bounded PIC code where three dimensional kinetic dynamics of blob is studied and current flow structures related to sheath formation are unveiled. The second one is the one spatial three velocity space dimensional (1D3V) PIC code with the Monte Carlo collisions where formation of detach plasma is studied. First target of our research is to construct self-consistent full kinetic simulation modeling of the linear divertor simulation experiments. This work is performed with the support and under the auspices of NIFS Collaboration Research program (NIFS15KNSS059, NIFS14KNXN279, and NIFS13KNSS038) and the Research Cooperation Program on Hierarchy and Holism in Natural Science at NINS.
Multiflow approach to plasma kinetics
Ignatov, A. M.
2015-10-15
Instead of the commonly used Vlasov equation, one is able to treat kinetic phenomena in collisionless plasma with the help of the infinite set of hydrodynamic equations. The present paper deals with the linear approximation of multiflow hydrodynamics. It is shown that single-particle and collective excitations analogous to Van Kampen waves are explicitly separated. Expressions for the energy of all eigenmodes are obtained.
Kinetic theory of relativistic plasmas
NASA Technical Reports Server (NTRS)
Gould, R. J.
1981-01-01
The thermalization of particle kinetic motion by binary collisions is considered for a plasma with a Boltzmann constant-temperature product approximately equal to 10 to 100 times the product of the electron mass with the square of the speed of light. At this temperature, the principal mechanism for relaxation of electron motion is via radiationless electron-electron collisions (Moller scattering). Ions are nonrelativistic, but are energetic enough so that their Coulomb scattering can be treated in the Born approximation. Relaxation times are computed and Boltzmann-equation Fokker-Planck operators are derived for the various binary-collision processes. The expression for the rate of kinetic energy exchange between electron and ion gases is derived for the case where the gases are at different temperatures.
Variance Anisotropy in Kinetic Plasmas
NASA Astrophysics Data System (ADS)
Parashar, Tulasi N.; Oughton, Sean; Matthaeus, William H.; Wan, Minping
2016-06-01
Solar wind fluctuations admit well-documented anisotropies of the variance matrix, or polarization, related to the mean magnetic field direction. Typically, one finds a ratio of perpendicular variance to parallel variance of the order of 9:1 for the magnetic field. Here we study the question of whether a kinetic plasma spontaneously generates and sustains parallel variances when initiated with only perpendicular variance. We find that parallel variance grows and saturates at about 5% of the perpendicular variance in a few nonlinear times irrespective of the Reynolds number. For sufficiently large systems (Reynolds numbers) the variance approaches values consistent with the solar wind observations.
Kinetic Profiles in NSTX Plasmas
R.E. Bell; B.P. LeBlanc; C. Bourdelle; D.R. Ernst; E.D. Fredrickson; D.A. Gates; J.C. Hosea; D.W. Johnson; S.M. Kaye; R. Maingi; S. Medley; J.E. Menard; D. Mueller; M. Ono; F. Paoletti; M. Peng; S.A. Sabbagh; D. Stutman; D.W. Swain; E.J. Synakowski; and J.R. Wilson
2001-07-10
The National Spherical Torus Experiment (NSTX) is a low aspect ratio (R/a approximately 1.3) device with auxiliary heating from neutral-beam injection (NBI) and high-harmonic fast-wave heating (HHFW). Typical NSTX parameters are R(subscript ''0'') = 85 cm, a = 67 cm, I(subscript ''p'') = 0.7-1.4 MA, B(subscript ''phi'') = 0.25-0.45 T. Three co-directed deuterium neutral-beam sources have injected P(subscript ''NB'') less than or equal to 4.7 MW. HHFW plasmas typically have delivered P(subscript ''RF'') less than or equal to 3 MW. Important to the understanding of NSTX confinement are the new kinetic profile diagnostics: a multi-pulse Thomson scattering system (MPTS) and a charge-exchange recombination spectroscopy (CHERS) system. The MPTS diagnostic currently measures electron density and temperature profiles at 30 Hz at ten spatial locations. The CHERS system has recently become available to measure carbon ion temperature and toroidal flow at 17 radial positions spanning the outer half of the minor radius with 20 msec time resolution during NBI. Experiments conducted during the last year have produced a wide range of kinetic profiles in NSTX. Some interesting examples are presented below.
Propagation of radiation in fluctuating multiscale plasmas. I. Kinetic theory
Tyshetskiy, Yu.; Pal Singh, Kunwar; Thirunavukarasu, A.; Robinson, P. A.; Cairns, Iver H.
2012-11-15
A theory for propagation of radiation in a large scale plasma with small scale fluctuations is developed using a kinetic description in terms of the probability distribution function of the radiation in space, time, and wavevector space. Large scale effects associated with spatial variations in the plasma density and refractive index of the plasma wave modes and small scale effects such as scattering of radiation by density clumps in fluctuating plasma, spontaneous emission, damping, and mode conversion are included in a multiscale kinetic description of the radiation. Expressions for the Stokes parameters in terms of the probability distribution function of the radiation are used to enable radiation properties such as intensity and polarization to be calculated.
Kinetic equation for spin-polarized plasmas
Cowley, S.C.; Kulsrud, R.M.; Valeo, E.
1984-07-01
The usual kinetic description of a plasma is extended to include variables to describe the spin. The distribution function, over phase-space and the new spin variables, provides a sufficient description of a spin-polarized plasma. The evolution equation for the distribution function is given. The equations derived are used to calculate depolarization due to four processes, inhomogeneous fields, collisions, collisions in inhomogeneous fields, and waves. It is found that depolarization by field inhomogeneity on scales large compared with the gyroradius is totally negligible. The same is true for collisional depolarization. Collisions in inhomogeneous fields yield a depolarization rate of order 10/sup -4/S/sup -1/ for deuterons and a negligible rate for tritons in a typical fusion reactor design. This is still sufficiently small on reactor time scales. However, small amplitude magnetic fluctuations (of order one gauss) resonant with the spin precession frequency can lead to significant depolarization (depolarises triton in ten seconds and deuteron in a hundred seconds.)
NLTE4 Plasma Population Kinetics Database
National Institute of Standards and Technology Data Gateway
SRD 159 NLTE4 Plasma Population Kinetics Database (Web database for purchase) This database contains benchmark results for simulation of plasma population kinetics and emission spectra. The data were contributed by the participants of the 4th Non-LTE Code Comparison Workshop who have unrestricted access to the database. The only limitation for other users is in hidden labeling of the output results. Guest users can proceed to the database entry page without entering userid and password.
Kinetic simulation of a plasma collision experiment
Larroche, O. )
1993-08-01
The ionic Fokker--Planck code which was written for describing plasma shock wave fronts [M. Casanova [ital et] [ital al]. Phys. Rev. Lett. [bold 67], 2143 (1991)] is applied to model the collision of two plasmas in plane geometry. Improvements brought to the code for that purpose are described. The initial phase of the experiment during which the plasmas interpenetrate is accounted for by a simple fluid model, which yields qualitative insight into the phenomena at play as well as an initial condition to start the kinetic simulation. The kinetic results obtained in the stagnation and thermalization phases are discussed with respect to a specific laser-produced plasma collision experiment, as well as to existing fluid and kinetic ( particle-in-cell'') simulations.
Kinetic simulation of a plasma collision experiment
NASA Astrophysics Data System (ADS)
Larroche, Olivier
1993-08-01
The ionic Fokker-Planck code which was written for describing plasma shock wave fronts [M. Casanova et al. Phys. Rev. Lett. 67, 2143 (1991)] is applied to model the collision of two plasmas in plane geometry. Improvements brought to the code for that purpose are described. The initial phase of the experiment during which the plasmas interpenetrate is accounted for by a simple fluid model, which yields qualitative insight into the phenomena at play as well as an initial condition to start the kinetic simulation. The kinetic results obtained in the stagnation and thermalization phases are discussed with respect to a specific laser-produced plasma collision experiment, as well as to existing fluid and kinetic (``particle-in-cell'') simulations.
MHD description of plasma: handbook of plasma physics
Kulsrud, R.M.
1980-10-01
The basic sets of MHD equations for the description of a plasma in various limits are derived and their usefulness and limits of validity are discussed. These limits are: the one fluid collisional plasma, the two fluid collisional plasma, the Chew-Goldberger Low formulation of the guiding center limit of a collisionless plasma and the double-adiabatic limit. Conservation relations are derived from these sets and the mathematics of the concept of flux freezing is given. An example is given illustrating the differences between guiding center theory and double adiabatic theory.
Neutral Vlasov kinetic theory of magnetized plasmas
Tronci, Cesare; Camporeale, Enrico
2015-02-15
The low-frequency limit of Maxwell equations is considered in the Maxwell-Vlasov system. This limit produces a neutral Vlasov system that captures essential features of plasma dynamics, while neglecting radiation effects. Euler-Poincaré reduction theory is used to show that the neutral Vlasov kinetic theory possesses a variational formulation in both Lagrangian and Eulerian coordinates. By construction, the new model recovers all collisionless neutral models employed in plasma simulations. Then, comparisons between the neutral Vlasov system and hybrid kinetic-fluid models are presented in the linear regime.
Inertial range turbulence in kinetic plasmas
Howes, Gregory G.
2008-05-15
The transfer of turbulent energy through an inertial range from the driving scale to dissipative scales in a kinetic plasma followed by the conversion of this energy into heat is a fundamental plasma physics process. A theoretical foundation for the study of this process is constructed, but the details of the kinetic cascade are not well understood. Several important properties are identified: (a) The conservation of a generalized energy by the cascade; (b) the need for collisions to increase entropy and realize irreversible plasma heating; and (c) the key role played by the entropy cascade--a dual cascade of energy to small scales in both physical and velocity space--to convert ultimately the turbulent energy into heat. A strategy for nonlinear numerical simulations of kinetic turbulence is outlined. Initial numerical results are consistent with the operation of the entropy cascade. Inertial range turbulence arises in a broad range of space and astrophysical plasmas and may play an important role in the thermalization of fusion energy in burning plasmas.
Kinetic Approach for Laser-Induced Plasmas
Omar, Banaz; Rethfeld, Baerbel
2008-10-22
Non-equilibrium distribution functions of electron gas and phonon gas excited with ultrashort intense laser pulses are calculated for laser-induced plasmas occurring in solids. The excitation during femtosecond irradiation and the subsequent thermalization of the free electrons, as well as the dynamics of phonons are described by kinetic equations. The microscopic collision processes, such as absorption by inverse bremsstrahlung, electron-electron collisions, and electron-phonon interactions are considered by complete Boltzmann collision integrals. We apply our kinetic approach for gold by taking s-band electron into account and compare it with the case of excitation of d-band electrons.
Linking plasma kinetics to plasma-bio interactions
NASA Astrophysics Data System (ADS)
Bruggeman, Peter
2015-05-01
Cold non-equilibrium atmospheric pressure plasmas have received a lot of attention in the last decade due to their huge potential for biomedical applications. In my group, we have characterized an RF driven APPJ in great detail. The characterization includes electrical measurements, imaging, optical emission spectroscopy, (two photon enhanced) laser induced fluorescence, Thomson scattering, Rayleigh scattering, Raman scattering and mass spectrometry. This led to a detailed knowledge of the electron density, electron temperature, gas temperature, NO, O, OH, O3 densities, ionic species and air concentrations in the plasma effluent. Living organisms for in vitro studies are typically kept in complex solutions or culture media. Plasma-bio interactions involves not only the production of reactive species in the plasma gas phase but also transport to the liquid phase and plasma induced liquid phase chemistry and its impact on the living organisms. Reactive nitrogen and oxygen species have been identified as the key reactive species. Recent results of my group show that controlling the gas phase plasma chemistry can lead to significant different biological responses of the living organisms corresponding to different chemical pathways. The effect of plasma jet interaction with liquids containing mammalian cells, bacteria and virus will be discussed. The outcomes of these studies allow unraveling chemical pathways responsible for plasma-bio interactions and linking plasma kinetics to plasma-bio interactions.
Kinetic equilibria of very high- β plasmas
NASA Astrophysics Data System (ADS)
Steinhauer, Loren; TAE Team
2015-11-01
Plasma equilibria with many large ion orbits, such as an advanced beam-driven field-reversed configuration, are neither static (Grad-Shafranov) nor describable as a flowing, multi-fluid. A fully-kinetic treatment of the ions is essential for such high- β plasmas. A kinetic equilibrium is needed to properly support realistic stability and transport analyses, both of which are strongly affected by large-orbit ions. A hybrid equilibrium model has been developed with a fully-kinetic treatment of both thermal ions and a rapidly-rotating ``beam-ion'' component, such as produced by neutral beam injection, relevant to the C-2U experiments at TAE. It employs analytic Vlasov solutions in that the distribution depends only on the two constants of motion, the Hamiltonian (H) and the canonical angular momentum (Pθ) . Electrons are treated as a pressure-bearing fluid. Since realistic forms of f (H ,Pθ) are affected by collisions, f is limited to solutions of a simplified Fokker-Planck equation. Importantly, a kinetic end-loss condition applies to unconfined ions, using a particle sink at a rate consistent with Monte-Carlo-like simulations of end loss accounting for a strong end mirror.
Kinetic effects in Enceladus plasma environment
NASA Astrophysics Data System (ADS)
Stverak, Stepan; Travnicek, Pavel M.; Sebek, Ondrej; Hellinger, Petr; Khurana, Krishan
2015-04-01
The southern plume of Enceladus represents a significant source of neutrals, ions and dust for the neighboring plasma environment and even for the inner magnetosphere of Saturn. The ion mass loading rate from the plume is not only confirmed by direct plasma measurements but can also be deduced from the strong signatures observed on the Kronian magnetospheric background magnetic field. In view of recent knowledge, namely based on in situ observations provided by numerous Cassini flybys, we try to model numerically and reconstruct the complex plasma environment in the vicinity of the moon by use of a full 3-dimensional hybrid code. With obtained results we attempt to investigate the role of kinetic effects in generation of the ion-cyclotron waves instantaneously observed by the Cassini spacecraft.
On kinetic dissipation in collisionless turbulent plasmas
NASA Astrophysics Data System (ADS)
Parashar, Tulasi Nandan
Plasma turbulence is a phenomenon that is present in astrophysical as well as terrestrial plasmas. The earth is embedded in a turbulent plasma, emitting from the sun, called the solar wind. It is important to understand the nature of this plasma in order to understand space weather. A critical unsolved problem is that of the source of dissipation in turbulent plasmas. It is believed to play a central role in the heating of the solar corona which in turn drives the solar wind. The solar wind itself is observed to be highly turbulent and hotter than predicted through adiabatic expansion models. Turbulence and its associated dissipation have been studied extensively through the use of MHD models. However, the solar wind and large regions of the solar corona have very low collisionality, which calls into question the use of simple viscosity and resistivity in most MHD models. A kinetic treatment is needed for a better understanding of turbulent dissipation. This thesis studies the dissipation of collisionless turbulence using direct numerical hybrid simulations of turbulent plasmas. Hybrid simulations use kinetic ions and fluid electrons. Having full kinetic ion physics, the dissipation in these simulations at the ion scales is self consistent and requires no assumptions. We study decaying as well as quasi steady state systems (driven magnetically). Initial studies of the Orszag-Tang vortex [Orszag, JFM, 1979] (which is an initial condition that quickly generates decaying strong turbulence) showed preferential perpendicular heating of protons (with T_perp /T_|| > 1). An energy budget analysis showed that in the turbulent regime, almost all the dissipation occurs through magnetic interactions. We study the energy budget of waves using the k - o spectra (energy in the wavenumber-frequency space). The k - o spectra of this study and subsequent studies of driven turbulent plasmas do not show any significant power in the linear wave modes of the system. This suggests that
Kinetics of wet sodium vapor complex plasma
Mishra, S. K.; Sodha, M. S.
2014-04-15
In this paper, we have investigated the kinetics of wet (partially condensed) Sodium vapor, which comprises of electrons, ions, neutral atoms, and Sodium droplets (i) in thermal equilibrium and (ii) when irradiated by light. The formulation includes the balance of charge over the droplets, number balance of the plasma constituents, and energy balance of the electrons. In order to evaluate the droplet charge, a phenomenon for de-charging of the droplets, viz., evaporation of positive Sodium ions from the surface has been considered in addition to electron emission and electron/ion accretion. The analysis has been utilized to evaluate the steady state parameters of such complex plasmas (i) in thermal equilibrium and (ii) when irradiated; the results have been graphically illustrated. As a significant outcome irradiated, Sodium droplets are seen to acquire large positive potential, with consequent enhancement in the electron density.
BOOK REVIEW: Kinetic theory of plasma waves, homogeneous plasmas
NASA Astrophysics Data System (ADS)
Porkolab, Miklos
1998-11-01
The linear theory of plasma waves in homogeneous plasma is arguably the most mature and best understood branch of plasma physics. Given the recently revised version of Stix's excellent Waves in Plasmas (1992), one might ask whether another book on this subject is necessary only a few years later. The answer lies in the scope of this volume; it is somewhat more detailed in certain topics than, and complementary in many fusion research relevant areas to, Stix's book. (I am restricting these comments to the homogeneous plasma theory only, since the author promises a second volume on wave propagation in inhomogeneous plasmas.) This book is also much more of a theorist's approach to waves in plasmas, with the aim of developing the subject within the logical framework of kinetic theory. This may indeed be pleasing to the expert and to the specialist, but may be too difficult to the graduate student as an `introduction' to the subject (which the author explicitly states in the Preface). On the other hand, it may be entirely appropriate for a second course on plasma waves, after the student has mastered fluid theory and an introductory kinetic treatment of waves in a hot magnetized `Vlasov' plasma. For teaching purposes, my personal preference is to review the cold plasma wave treatment using the unified Stix formalism and notation (which the author wisely adopts in the present book, but only in Chapter 5). Such an approach allows one to deal with CMA diagrams early on, as well as to provide a framework to discuss electromagnetic wave propagation and accessibility in inhomogeneous plasmas (for which the cold plasma wave treatment is perfectly adequate). Such an approach does lack some of the rigour, however, that the author achieves with the present approach. As the author correctly shows, the fluid theory treatment of waves follows logically from kinetic theory in the cold plasma limit. I only question the pedagogical value of this approach. Otherwise, I welcome this
Fundamental Statistical Descriptions of Plasma Turbulence in Magnetic Fields
John A. Krommes
2001-02-16
A pedagogical review of the historical development and current status (as of early 2000) of systematic statistical theories of plasma turbulence is undertaken. Emphasis is on conceptual foundations and methodology, not practical applications. Particular attention is paid to equations and formalism appropriate to strongly magnetized, fully ionized plasmas. Extensive reference to the literature on neutral-fluid turbulence is made, but the unique properties and problems of plasmas are emphasized throughout. Discussions are given of quasilinear theory, weak-turbulence theory, resonance-broadening theory, and the clump algorithm. Those are developed independently, then shown to be special cases of the direct-interaction approximation (DIA), which provides a central focus for the article. Various methods of renormalized perturbation theory are described, then unified with the aid of the generating-functional formalism of Martin, Siggia, and Rose. A general expression for the renormalized dielectric function is deduced and discussed in detail. Modern approaches such as decimation and PDF methods are described. Derivations of DIA-based Markovian closures are discussed. The eddy-damped quasinormal Markovian closure is shown to be nonrealizable in the presence of waves, and a new realizable Markovian closure is presented. The test-field model and a realizable modification thereof are also summarized. Numerical solutions of various closures for some plasma-physics paradigms are reviewed. The variational approach to bounds on transport is developed. Miscellaneous topics include Onsager symmetries for turbulence, the interpretation of entropy balances for both kinetic and fluid descriptions, self-organized criticality, statistical interactions between disparate scales, and the roles of both mean and random shear. Appendices are provided on Fourier transform conventions, dimensional and scaling analysis, the derivations of nonlinear gyrokinetic and gyrofluid equations
Kinetic Simulations of Dense Plasma Focus Breakdown
NASA Astrophysics Data System (ADS)
Schmidt, A.; Higginson, D. P.; Jiang, S.; Link, A.; Povilus, A.; Sears, J.; Bennett, N.; Rose, D. V.; Welch, D. R.
2015-11-01
A dense plasma focus (DPF) device is a type of plasma gun that drives current through a set of coaxial electrodes to assemble gas inside the device and then implode that gas on axis to form a Z-pinch. This implosion drives hydrodynamic and kinetic instabilities that generate strong electric fields, which produces a short intense pulse of x-rays, high-energy (>100 keV) electrons and ions, and (in deuterium gas) neutrons. A strong factor in pinch performance is the initial breakdown and ionization of the gas along the insulator surface separating the two electrodes. The smoothness and isotropy of this ionized sheath are imprinted on the current sheath that travels along the electrodes, thus making it an important portion of the DPF to both understand and optimize. Here we use kinetic simulations in the Particle-in-cell code LSP to model the breakdown. Simulations are initiated with neutral gas and the breakdown modeled self-consistently as driven by a charged capacitor system. We also investigate novel geometries for the insulator and electrodes to attempt to control the electric field profile. The initial ionization fraction of gas is explored computationally to gauge possible advantages of pre-ionization which could be created experimentally via lasers or a glow-discharge. Prepared by LLNL under Contract DE-AC52-07NA27344.
Kinetic theory of partially ionized complex (dusty) plasmas
Tsytovich, V.N.; De Angelis, U.; Ivlev, A.V.; Morfill, G.E.
2005-08-15
The general approach to the kinetic theory of complex (dusty) plasmas [Tsytovich and de Angelis, Phys. Plasmas 6, 1093 (1999)], which was formulated with the assumption of a regular (nonfluctuating) source of plasma particles, is reformulated to include ionization by electron impact on neutrals as the plasma source and the effects of collisions of ions and dust particles with neutrals.
Nonlinear Kinetic Instabilities in Plasma Wakes
NASA Astrophysics Data System (ADS)
Hutchinson, I. H.; Haakonsen, C. B.
2015-12-01
Relative motion of a plasma and an embedded perturbing solid objectproduces a plasma wake, which is kinetically unstable. For moons,asteroids, spacecraft, probes, and planets without a magnetosphere theresponse is dominantly electrostatic, although generally with abackground magnetic field. Using high-fidelity particle-in-cellsimulations, we have observed the development of kinetic instabilitiesand their non-linear consequences in representative wakes. We havealso explained the observations with semi-analytical non-lineartheory. The ion and electron distribution function shapes are stronglyperturbed in the wake region. The ions form two opposite beamsdirected inward along the guiding magnetic field, in part because ofthe attraction of the wake's electric potential well. The electrondistribution forms a notch or dimple (of reduced phase space density)localized in velocity to orbits that dwell near the wake axis (becauseof repulsion). Those orbits are de-energized by cross-field drift downthe potential-energy ridge. The resulting Langmuir instability spawnselectron holes. The holes that move faster than the ion beams areaccelerated out of the wake by its electrostatic field without growingsubstantially. Some holes, however, remain in the wake at essentiallyzero parallel velocity. They grow, as a result of the same mechanismthat formed the notch: cross-field drift from a lower to a higherdensity. When the density rises by a factor of order two or three,they grow large enough to perturb the ions, tap their free energy, anddisrupt the ion streams well before they would become ion-ionunstable. Crucially, these processes depend strongly on theion/electron mass ratio and require close to physical ratio (1836) insimulations, to reveal their characteristics. Electron holes arisingfrom these processes may be widely present and observable in spaceplasma wakes.
Consistent description of kinetic equation with triangle anomaly
Pu Shi; Gao Jianhua; Wang Qun
2011-05-01
We provide a consistent description of the kinetic equation with a triangle anomaly which is compatible with the entropy principle of the second law of thermodynamics and the charge/energy-momentum conservation equations. In general an anomalous source term is necessary to ensure that the equations for the charge and energy-momentum conservation are satisfied and that the correction terms of distribution functions are compatible to these equations. The constraining equations from the entropy principle are derived for the anomaly-induced leading order corrections to the particle distribution functions. The correction terms can be determined for the minimum number of unknown coefficients in one charge and two charge cases by solving the constraining equations.
Ducted kinetic Alfven waves in plasma with steep density gradients
Houshmandyar, Saeid; Scime, Earl E.
2011-11-15
Given their high plasma density (n {approx} 10{sup 13} cm{sup -3}), it is theoretically possible to excite Alfven waves in a conventional, moderate length (L {approx} 2 m) helicon plasma source. However, helicon plasmas are decidedly inhomogeneous, having a steep radial density gradient, and typically have a significant background neutral pressure. The inhomogeneity introduces regions of kinetic and inertial Alfven wave propagation. Ion-neutral and electron-neutral collisions alter the Alfven wave dispersion characteristics. Here, we present the measurements of propagating kinetic Alfven waves in helium helicon plasma. The measured wave dispersion is well fit with a kinetic model that includes the effects of ion-neutral damping and that assumes the high density plasma core defines the radial extent of the wave propagation region. The measured wave amplitude versus plasma radius is consistent with the pile up of wave magnetic energy at the boundary between the kinetic and inertial regime regions.
Solitary kinetic Alfven waves in dusty plasmas
Li Yangfang; Wu, D. J.; Morfill, G. E.
2008-08-15
Solitary kinetic Alfven waves in dusty plasmas are studied by considering the dust charge variation. The effect of the dust charge-to-mass ratio on the soliton solution is discussed. The Sagdeev potential is derived analytically with constant dust charge and then calculated numerically by taking the dust charge variation into account. We show that the dust charge-to-mass ratio plays an important role in the soliton properties. The soliton solutions are comprised of two branches. One branch is sub-Alfvenic and the soliton velocity is obviously smaller than the Alfven speed. The other branch is super-Alfvenic and the soliton velocity is very close to or greater than the Alfven speed. Both compressive and rarefactive solitons can exist. For the sub-Alfvenic branch, the rarefactive soliton is bell-shaped and it is much narrower than the compressive one. However, for the super-Alfvenic branch, the compressive soliton is bell-shaped and narrower, and the rarefactive one is broadened. When the charge-to-mass ratio of the dust grains is sufficiently high, the width of the rarefactive soliton, in the super-Alfvenic branch, will broaden extremely and a electron depletion will be observed. It is also shown that the bell-shaped soliton can transition to a cusped structure when the velocity is sufficiently high.
Fully implicit kinetic modelling of collisional plasmas
Mousseau, V.A.
1996-05-01
This dissertation describes a numerical technique, Matrix-Free Newton Krylov, for solving a simplified Vlasov-Fokker-Planck equation. This method is both deterministic and fully implicit, and may not have been a viable option before current developments in numerical methods. Results are presented that indicate the efficiency of the Matrix-Free Newton Krylov method for these fully-coupled, nonlinear integro-differential equations. The use and requirement for advanced differencing is also shown. To this end, implementations of Chang-Cooper differencing and flux limited Quadratic Upstream Interpolation for Convective Kinematics (QUICK) are presented. Results are given for a fully kinetic ion-electron problem with a self consistent electric field calculated from the ion and electron distribution functions. This numerical method, including advanced differencing, provides accurate solutions, which quickly converge on workstation class machines. It is demonstrated that efficient steady-state solutions can be achieved to the non-linear integro-differential equation, obtaining quadratic convergence, without incurring the large memory requirements of an integral operator. Model problems are presented which simulate plasma impinging on a plate with both high and low neutral particle recycling typical of a divertor in a Tokamak device. These model problems demonstrate the performance of the new solution method.
Simulations of plasma sheaths using continuum kinetic models
NASA Astrophysics Data System (ADS)
Srinivasan, Bhuvana; Hakim, Ammar
2015-11-01
Understanding plasma sheath physics is important for the performance of devices such as Hall thrusters due to the effect of energetic particles on electrode erosion. Plasma sheath physics is studied using kinetic and multi-fluid models with relevance to secondary electron emissions and plasma-surface interactions. Continuum kinetic models are developed to directly solve the Vlasov-Poisson equation using the discontinuous Galerkin method for each of the ion and electron species. A steady-state sheath is simulated by including a simple model for a neutral fluid. Multi-fluid simulations for the plasma sheath are also performed using the discontinuous Galerkin method to solve a complete set of fluid equations for each of the ion and electron species. The kinetic plasma sheath is compared to a multi-fluid plasma sheath. Supported by Air Force Office of Scientific Research.
Kinetic model for the collisionless sheath of a collisional plasma
NASA Astrophysics Data System (ADS)
Tang, Xian-Zhu; Guo, Zehua
2016-08-01
Collisional plasmas typically have mean-free-path still much greater than the Debye length, so the sheath is mostly collisionless. Once the plasma density, temperature, and flow are specified at the sheath entrance, the profile variation of electron and ion density, temperature, flow speed, and conductive heat fluxes inside the sheath is set by collisionless dynamics, and can be predicted by an analytical kinetic model distribution. These predictions are contrasted here with direct kinetic simulations, showing good agreement.
Kinetic instability of ion acoustic mode in permeating plasmas
Vranjes, J.; Poedts, S.; Ehsan, Zahida
2009-07-15
In plasmas with electron drift (current) relative to static ions, the ion acoustic wave is subject to the kinetic instability which takes place if the directed electron speed exceeds the ion acoustic speed. The instability threshold becomes different in the case of one quasineutral electron-ion plasma propagating through another static quasineutral (target) plasma. The threshold velocity of the propagating plasma may be well below the ion acoustic speed of the static plasma. Such a currentless instability may frequently be expected in space and astrophysical plasmas.
Kinetic signatures and intermittent turbulence in the solar wind plasma.
Osman, K T; Matthaeus, W H; Hnat, B; Chapman, S C
2012-06-29
A connection between kinetic processes and intermittent turbulence is observed in the solar wind plasma using measurements from the Wind spacecraft at 1 A.U. In particular, kinetic effects such as temperature anisotropy and plasma heating are concentrated near coherent structures, such as current sheets, which are nonuniformly distributed in space. Furthermore, these coherent structures are preferentially found in plasma unstable to the mirror and firehose instabilities. The inhomogeneous heating in these regions, which is present in both the magnetic field parallel and perpendicular temperature components, results in protons at least 3-4 times hotter than under typical stable plasma conditions. These results offer a new understanding of kinetic processes in a turbulent regime, where linear Vlasov theory is not sufficient to explain the inhomogeneous plasma dynamics operating near non-Gaussian structures. PMID:23004954
Kinetic description of an electron--LO-phonon system with finite phonon lifetime
Nguyen, V.T.; Mahler, G. )
1992-02-15
We study the cooling of an electron plasma from a kinetic point of view. For this purpose, a quantum theory of fluctuations is applied to derive the kinetic equations for an electron--LO-phonon system from various model Hamiltonians. A polarization approximation is provided that goes beyond perturbation theory of the electron-phonon interaction. The description of electron-phonon energy exchange is shown to be impossible with the interacting Hamiltonian in Froehlich's one-phonon form unless dissipation of the bare LO phonon is included. For a Hamiltonian including effects of the scattering of LO phonons by acoustic phonons, kinetic equations are derived. The equation for LO phonons is shown to describe the collective excitations with finite lifetime, in the limiting case of weak damping of the plasmon-phonon coupled modes. A reduction of the cooling rate similar to the hot-phonon'' effect is shown to occur for the case of weak coupling without assuming a steady state of the LO phonons. Finally, an electron-phonon interaction Hamiltonian in two-phonon form is considered and it is shown that electron-phonon energy exchange may be described in the polarization approximation without introducing a finite phonon lifetime.
Kinetic Modeling of the Moon-Solar Wind Plasma Interaction
NASA Astrophysics Data System (ADS)
Fatemi, S.; Poppe, A. R.; Halekas, J. S.; Delory, G. T.; Holmstrom, M.; Farrell, W. M.
2016-05-01
We use a three-dimensional self-consistent hybrid model of plasma (kinetic ions, fluid electrons) to study solar wind plasma interaction with the Moon. We have studied lunar wake, interaction with crustal fields, and lunar interior with our model.
Intermittent dissipation at kinetic scales in collisionless plasma turbulence.
Wan, M; Matthaeus, W H; Karimabadi, H; Roytershteyn, V; Shay, M; Wu, P; Daughton, W; Loring, B; Chapman, S C
2012-11-01
High resolution kinetic simulations of collisionless plasma driven by shear show the development of turbulence characterized by dynamic coherent sheetlike current density structures spanning a range of scales down to electron scales. We present evidence that these structures are sites for heating and dissipation, and that stronger current structures signify higher dissipation rates. Evidently, kinetic scale plasma, like magnetohydrodynamics, becomes intermittent due to current sheet formation, leading to the expectation that heating and dissipation in astrophysical and space plasmas may be highly nonuniform and patchy. PMID:23215389
Kinetic study of ion-acoustic plasma vortices
Khan, S. A.; Aman-ur-Rehman; Mendonca, J. T.
2014-09-15
The kinetic theory of electron plasma waves with finite orbital angular momentum has recently been introduced by Mendonca. This model shows possibility of new kind of plasma waves and instabilities. We have extended the theory to ion-acoustic plasma vortices carrying orbital angular momentum. The dispersion equation is derived under paraxial approximation which exhibits a kind of linear vortices and their Landau damping. The numerical solutions are obtained and compared with analytical results which are in good agreement. The physical interpretation of the ion-acoustic plasma vortices and their Landau resonance conditions are given for typical case of Maxwellian plasmas.
Kinetic treatment of nonlinear magnetized plasma motions - General geometry and parallel waves
NASA Technical Reports Server (NTRS)
Khabibrakhmanov, I. KH.; Galinskii, V. L.; Verheest, F.
1992-01-01
The expansion of kinetic equations in the limit of a strong magnetic field is presented. This gives a natural description of the motions of magnetized plasmas, which are slow compared to the particle gyroperiods and gyroradii. Although the approach is 3D, this very general result is used only to focus on the parallel propagation of nonlinear Alfven waves. The derivative nonlinear Schroedinger-like equation is obtained. Two new terms occur compared to earlier treatments, a nonlinear term proportional to the heat flux along the magnetic field line and a higher-order dispersive term. It is shown that kinetic description avoids the singularities occurring in magnetohydrodynamic or multifluid approaches, which correspond to the degenerate case of sound speeds equal to the Alfven speed, and that parallel heat fluxes cannot be neglected, not even in the case of low parallel plasma beta. A truly stationary soliton solution is derived.
Kinetic modelling of runaway electron avalanches in tokamak plasmas
NASA Astrophysics Data System (ADS)
Nilsson, E.; Decker, J.; Peysson, Y.; Granetz, R. S.; Saint-Laurent, F.; Vlainic, M.
2015-09-01
Runaway electrons can be generated in tokamak plasmas if the accelerating force from the toroidal electric field exceeds the collisional drag force owing to Coulomb collisions with the background plasma. In ITER, disruptions are expected to generate runaway electrons mainly through knock-on collisions (Hender et al 2007 Nucl. Fusion 47 S128-202), where enough momentum can be transferred from existing runaways to slow electrons to transport the latter beyond a critical momentum, setting off an avalanche of runaway electrons. Since knock-on runaways are usually scattered off with a significant perpendicular component of the momentum with respect to the local magnetic field direction, these particles are highly magnetized. Consequently, the momentum dynamics require a full 3D kinetic description, since these electrons are highly sensitive to the magnetic non-uniformity of a toroidal configuration. For this purpose, a bounce-averaged knock-on source term is derived. The generation of runaway electrons from the combined effect of Dreicer mechanism and knock-on collision process is studied with the code LUKE, a solver of the 3D linearized bounce-averaged relativistic electron Fokker-Planck equation (Decker and Peysson 2004 DKE: a fast numerical solver for the 3D drift kinetic equation Report EUR-CEA-FC-1736, Euratom-CEA), through the calculation of the response of the electron distribution function to a constant parallel electric field. The model, which has been successfully benchmarked against the standard Dreicer runaway theory now describes the runaway generation by knock-on collisions as proposed by Rosenbluth (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355-62). This paper shows that the avalanche effect can be important even in non-disruptive scenarios. Runaway formation through knock-on collisions is found to be strongly reduced when taking place off the magnetic axis, since trapped electrons can not contribute to the runaway electron population. Finally, the
Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions
NASA Technical Reports Server (NTRS)
Aunai, icolas; Belmont, Gerard; Smets, Roch
2012-01-01
Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to
Pedestal Fueling Simulations with a Coupled Kinetic-kinetic Plasma-neutral Transport Code
D.P. Stotler, C.S. Chang, S.H. Ku, J. Lang and G.Y. Park
2012-08-29
A Monte Carlo neutral transport routine, based on DEGAS2, has been coupled to the guiding center ion-electron-neutral neoclassical PIC code XGC0 to provide a realistic treatment of neutral atoms and molecules in the tokamak edge plasma. The DEGAS2 routine allows detailed atomic physics and plasma-material interaction processes to be incorporated into these simulations. The spatial pro le of the neutral particle source used in the DEGAS2 routine is determined from the uxes of XGC0 ions to the material surfaces. The kinetic-kinetic plasma-neutral transport capability is demonstrated with example pedestal fueling simulations.
Energy Conservation Tests of a Coupled Kinetic-kinetic Plasma-neutral Transport Code
Stotler, D. P.; Chang, C. S.; Ku, S. H.; Lang, J.; Park, G.
2012-08-29
A Monte Carlo neutral transport routine, based on DEGAS2, has been coupled to the guiding center ion-electron-neutral neoclassical PIC code XGC0 to provide a realistic treatment of neutral atoms and molecules in the tokamak edge plasma. The DEGAS2 routine allows detailed atomic physics and plasma-material interaction processes to be incorporated into these simulations. The spatial pro le of the neutral particle source used in the DEGAS2 routine is determined from the uxes of XGC0 ions to the material surfaces. The kinetic-kinetic plasma-neutral transport capability is demonstrated with example pedestal fueling simulations.
Theory of spatially non-symmetric kinetic equilibria for collisionless plasmas
Cremaschini, Claudio; Tessarotto, Massimo
2013-01-15
The problem posed by the possible existence/non-existence of spatially non-symmetric kinetic equilibria has remained unsolved in plasma theory. For collisionless magnetized plasmas, this involves the construction of stationary solutions of the Vlasov-Maxwell equations. In this paper, the issue is addressed for non-relativistic plasmas both in astrophysical and laboratory contexts. The treatment is based on a Lagrangian variational description of single-particle dynamics. Starting point is a non-perturbative formulation of gyrokinetic theory, which allows one to construct 'a posteriori' with prescribed order of accuracy an asymptotic representation for the magnetic moment. In terms of the relevant particle adiabatic invariants generalized bi-Maxwellian equilibria are proved to exist. These are shown to recover, under suitable assumptions, a Chapman-Enskog form which permits an analytical treatment of the corresponding fluid moments. In particular, the constrained posed by the Poisson and the Ampere equations are analyzed, both for quasi-neutral and non-neutral plasmas. The conditions of existence of the corresponding non-symmetric kinetic equilibria are investigated. As a notable feature, both astrophysical and laboratory plasmas are shown to exhibit, under suitable conditions, a kinetic dynamo, whereby the equilibrium magnetic field can be self-generated by the equilibrium plasma currents.
Kinetic theory and simulation of multi-species plasmas in tokamaks excited with ICRF microwaves
Kerbel, G.D.; McCoy, M.G.
1984-12-21
This paper presents a description of a bounce-averaged Fokker-Planck quasilinear model for the kinetic description of tokamak plasmas. The non-linear collision and quasilinear resonant diffusion operators are represented in a form conducive to numerical solution with specific attention to the treatment of the boundary layer separating trapped and passing orbit regions of velocity space. The numerical techniques employed are detailed in so far as they constitute significant departure from those used in the conventional uniform magnetic field case. Examples are given to illustrate the combined effects of collisional and resonant diffusion.
Intermittent dissipation and heating in 3D kinetic plasma turbulence
NASA Astrophysics Data System (ADS)
Wan, M.; Matthaeus, W. H.; Roytershteyn, V.; Karimabadi, H.; Parashar, T.; Wu, P.; Shay, M. A.
2014-12-01
The nature of collisionless dissipation has been hotlydebated in recent years, with alternative ideas posed interms of various wave modes, such as kinetic Alfven waves,whistlers, linear Vlasov instabilities, cyclotron resonance,and Landau damping. Here we use large scale, fully kinetic3D simulations of collisionless plasma turbulence which showthe development of turbulence characterized by sheet-likecurrent density structures spanning a range of scales.We present evidence that these structures are sites for heatingand dissipation, and that stronger current structures signifyhigher dissipation rates. The analyses focus on quantities such as J.E, electron and proton temperatures, and PVI of the magnetic field. Evidently, kinetic scale plasma,like magnetohydrodynamics, becomes intermittent due tocurrent sheet formation, leading to the expectationthat heating and dissipation in astrophysical and space plasmasmay be highly nonuniform. Comparison with previousresults from 2D kinetic simulations, as well as high frequencysolar wind observational data will also be discussed.
Kinetic Simulations of Ladder Climbing and Autoresonance of Plasma Waves
NASA Astrophysics Data System (ADS)
Kaminski, Erez; Barth, Ido; Fisch, Nat; Dodin, Ilya
2015-11-01
Quantum like Ladder Climbing and Autoresonance of classical Langmuir waves in bounded plasmas are numerically studied within a kinetic model and compared with earlier fluid model simulations. Both dynamical solutions are excited and controlled via chirped modulations of the background density that preserve the plasma wave quanta. Landau damping determines the system's maximal stable level, imposing a kinetic limit on the maximal level of the Ladder Climbing or Autoresonance dynamics. Vlasov simulations are employed to test the kinetic stability of both dynamics and to find the kinetic limit for different system's parameters. This work was Supported by NNSA grant DE274-FG52-08NA28553, DOE contract DE-AC02-09CH11466, and DTRA grant HDTRA1-11-1-0037.
Suppression of phase mixing in drift-kinetic plasma turbulence
NASA Astrophysics Data System (ADS)
Parker, J. T.; Highcock, E. G.; Schekochihin, A. A.; Dellar, P. J.
2016-07-01
Transfer of free energy from large to small velocity-space scales by phase mixing leads to Landau damping in a linear plasma. In a turbulent drift-kinetic plasma, this transfer is statistically nearly canceled by an inverse transfer from small to large velocity-space scales due to "anti-phase-mixing" modes excited by a stochastic form of plasma echo. Fluid moments (density, velocity, and temperature) are thus approximately energetically isolated from the higher moments of the distribution function, so phase mixing is ineffective as a dissipation mechanism when the plasma collisionality is small.
Intermittent Dissipation and Heating in 3D Kinetic Plasma Turbulence.
Wan, M; Matthaeus, W H; Roytershteyn, V; Karimabadi, H; Parashar, T; Wu, P; Shay, M
2015-05-01
High resolution, fully kinetic, three dimensional (3D) simulation of collisionless plasma turbulence shows the development of turbulence characterized by sheetlike current density structures spanning a range of scales. The nonlinear evolution is initialized with a long wavelength isotropic spectrum of fluctuations having polarizations transverse to an imposed mean magnetic field. We present evidence that these current sheet structures are sites for heating and dissipation, and that stronger currents signify higher dissipation rates. The analyses focus on quantities such as J·E, electron, and proton temperatures, and conditional averages of these quantities based on local electric current density. Evidently, kinetic scale plasma, like magnetohydrodynamics, becomes intermittent due to current sheet formation, leading to the expectation that heating and dissipation in astrophysical and space plasmas may be highly nonuniform. Comparison with previous results from 2D kinetic simulations, as well as high frequency solar wind observational data, are discussed. PMID:25978241
Kinetic theory of electromagnetic ion waves in relativistic plasmas
Marklund, Mattias; Shukla, Padma K.
2006-09-15
A kinetic theory for electromagnetic ion waves in a cold relativistic plasma is derived. The kinetic equation for the broadband electromagnetic ion waves is coupled to the slow density response via an acoustic equation driven by a ponderomotive force-like term linear in the electromagnetic field amplitude. The modulational instability growth rate is derived for an arbitrary spectrum of waves. The monochromatic and random phase cases are studied.
Hamiltonian field description of two-dimensional vortex fluids and guiding center plasmas
Morrison, P.J.
1981-03-01
The equations that describe the motion of two-dimensional vortex fluids and guiding center plasmas are shown to possess underlying field Hamiltonian structure. A Poisson bracket which is given in terms of the vorticity, the physical although noncanonical dynamical variable, casts these equations into Heisenberg form. The Hamiltonian density is the kinetic energy density of the fluid. The well-known conserved quantities are seen to be in involution with respect to this Poisson bracket. Expanding the vorticity in terms of a Fourier-Dirac series transforms the field description given here into the usual canonical equations for discrete vortex motion. A Clebsch potential representation of the vorticity transforms the noncanonical field description into a canonical description.
Collisionless kinetic regimes for quasi-stationary axisymmetric accretion disc plasmas
Cremaschini, C.; Tessarotto, M.
2012-08-15
This paper is concerned with the kinetic treatment of quasi-stationary axisymmetric collisionless accretion disc plasmas. The conditions of validity of the kinetic description for non-relativistic magnetized and gravitationally bound plasmas of this type are discussed. A classification of the possible collisionless plasma regimes which can arise in these systems is proposed, which can apply to accretion discs around both stellar-mass compact objects and galactic-center black holes. Two different classifications are determined, which are referred to, respectively, as energy-based and magnetic field-based classifications. Different regimes are pointed out for each plasma species, depending both on the relative magnitudes of kinetic and potential energies and the magnitude of the magnetic field. It is shown that in all cases, there can be quasi-stationary Maxwellian-like solutions of the Vlasov equation. The perturbative approach outlined here permits unique analytical determination of the functional form for the distribution function consistent, in each kinetic regime, with the explicit inclusion of finite Larmor radius-diamagnetic and/or energy-correction effects.
Fully kinetic simulations of megajoule-scale dense plasma focus
Schmidt, A.; Link, A.; Tang, V.; Halvorson, C.; May, M.; Welch, D.; Meehan, B. T.; Hagen, E. C.
2014-10-15
Dense plasma focus (DPF) Z-pinch devices are sources of copious high energy electrons and ions, x-rays, and neutrons. Megajoule-scale DPFs can generate 10{sup 12} neutrons per pulse in deuterium gas through a combination of thermonuclear and beam-target fusion. However, the details of the neutron production are not fully understood and past optimization efforts of these devices have been largely empirical. Previously, we reported on the first fully kinetic simulations of a kilojoule-scale DPF and demonstrated that both kinetic ions and kinetic electrons are needed to reproduce experimentally observed features, such as charged-particle beam formation and anomalous resistivity. Here, we present the first fully kinetic simulation of a MegaJoule DPF, with predicted ion and neutron spectra, neutron anisotropy, neutron spot size, and time history of neutron production. The total yield predicted by the simulation is in agreement with measured values, validating the kinetic model in a second energy regime.
Kinetics of complex plasma with liquid droplets
Misra, Shikha; Sodha, M. S.; Mishra, S. K.
2013-12-15
This paper provides a theoretical basis for the reduction of electron density by spray of water (or other liquids) in hot plasma. This phenomenon has been observed in a hypersonic flight experiment for relief of radio black out, caused by high ionization in the plasma sheath of a hypersonic vehicle, re-entering the atmosphere. The analysis incorporates a rather little known phenomenon for de-charging of the droplets, viz., evaporation of ions from the surface and includes the charge balance on the droplets and number cum energy balance of electrons, ions, and neutral molecules; the energy balance of the evaporating droplets has also been taken into account. The analysis has been applied to a realistic situation and the transient variations of the charge and radius of water droplets, and other plasma parameters have been obtained and discussed. The analysis through made in the context of water droplets is applicable to all liquids.
Kinetics of complex plasma with liquid droplets
NASA Astrophysics Data System (ADS)
Misra, Shikha; Mishra, S. K.; Sodha, M. S.
2013-12-01
This paper provides a theoretical basis for the reduction of electron density by spray of water (or other liquids) in hot plasma. This phenomenon has been observed in a hypersonic flight experiment for relief of radio black out, caused by high ionization in the plasma sheath of a hypersonic vehicle, re-entering the atmosphere. The analysis incorporates a rather little known phenomenon for de-charging of the droplets, viz., evaporation of ions from the surface and includes the charge balance on the droplets and number cum energy balance of electrons, ions, and neutral molecules; the energy balance of the evaporating droplets has also been taken into account. The analysis has been applied to a realistic situation and the transient variations of the charge and radius of water droplets, and other plasma parameters have been obtained and discussed. The analysis through made in the context of water droplets is applicable to all liquids.
Kinetic simulation of a collisional shock wave in a plasma
Casanova, M.; Larroche, O. ); Matte, J. )
1991-10-14
The ion kinetic structure of a planar collisional shock front in a fully ionized plasma is investigated using a new Vlasov-Fokker-Planck code. The effects of ionic viscosity and ionic thermal conduction are found to be much larger than assumed in usual hydrodynamic plasma simulations with classical transport coefficients. This might have consequences on the numerical modeling of inertial-confinement fusion targets.
Decharging of complex plasmas: first kinetic observations.
Ivlev, A V; Kretschmer, M; Zuzic, M; Morfill, G E; Rothermel, H; Thomas, H M; Fortov, V E; Molotkov, V I; Nefedov, A P; Lipaev, A M; Petrov, O F; Baturin, Yu M; Ivanov, A I; Goree, J
2003-02-01
The first experiment on the decharging of a complex plasma in microgravity conditions was conducted. After switching off the rf power, in the afterglow plasma, ions and electrons rapidly recombine and leave a cloud of charged microparticles. Because of microgravity, the particles remain suspended in the experimental chamber for a sufficiently long time, allowing precise measurements of the rest particle charge. A simple theoretical model for the decharging is proposed which agrees quite well with the experiment results and predicts the rest charge at lower gas pressures. PMID:12633365
Spin Kinetic Models of Plasmas - Semiclassical and Quantum Mechanical Theory
Brodin, Gert; Marklund, Mattias; Zamanian, Jens
2009-11-10
In this work a recently published semiclassical spin kinetic model, generalizing those of previous authors are discussed. Some previously described properties are reviewed, and a new example illustrating the theory is presented. The generalization to a fully quantum mechanical description is discussed, and the main features of such a theory is outlined. Finally, the main conclusions are presented.
Modeling the heating and atomic kinetics of a photoionized neon plasma experiment
NASA Astrophysics Data System (ADS)
Lockard, Tom E.
Motivated by gas cell photoionized plasma experiments performed by our group at the Z facility of Sandia National Laboratories, we discuss in this dissertation a modeling study of the heating and ionization of the plasma for conditions characteristic of these experiments. Photoionized plasmas are non-equilibrium systems driven by a broadband x-ray radiation flux. They are commonly found in astrophysics but rarely seen in the laboratory. Several modeling tools have been employed: (1) a view-factor computer code constrained with side x-ray power and gated monochromatic image measurements of the z-pinch radiation, to model the time-history of the photon-energy resolved x-ray flux driving the photoionized plasma, (2) a Boltzmann self-consistent electron and atomic kinetics model to simulate the electron distribution function and configuration-averaged atomic kinetics, (3) a radiation-hydrodynamics code with inline non-equilibrium atomic kinetics to perform a comprehensive numerical simulation of the experiment and plasma heating, and (4) steady-state and time-dependent collisional-radiative atomic kinetics calculations with fine-structure energy level description to assess transient effects in the ionization and charge state distribution of the plasma. The results indicate that the photon-energy resolved x-ray flux impinging on the front window of the gas cell is very well approximated by a linear combination of three geometrically-diluted Planckian distributions. Knowledge of the spectral details of the x-ray drive turned out to be important for the heating and ionization of the plasma. The free electrons in the plasma thermalize quickly relative to the timescales associated with the time-history of the x-ray drive and the plasma atomic kinetics. Hence, electrons are well described by a Maxwellian energy distribution of a single temperature. This finding is important to support the application of a radiation-hydrodynamic model to simulate the experiment. It is found
Cremaschini, Claudio; Kovář, Jiří; Slaný, Petr; Stuchlík, Zdeněk; Karas, Vladimír
2013-11-01
The possible occurrence of equilibrium off-equatorial tori in the gravitational and electromagnetic fields of astrophysical compact objects has been recently proved based on non-ideal magnetohydrodynamic theory. These stationary structures can represent plausible candidates for the modeling of coronal plasmas expected to arise in association with accretion disks. However, accretion disk coronae are formed by a highly diluted environment, and so the fluid description may be inappropriate. The question is posed of whether similar off-equatorial solutions can also be determined in the case of collisionless plasmas for which treatment based on kinetic theory, rather than a fluid one, is demanded. In this paper the issue is addressed in the framework of the Vlasov-Maxwell description for non-relativistic, multi-species axisymmetric plasmas subject to an external dominant spherical gravitational and dipolar magnetic field. Equilibrium configurations are investigated and explicit solutions for the species kinetic distribution function are constructed, which are expressed in terms of generalized Maxwellian functions characterized by isotropic temperature and non-uniform fluid fields. The conditions for the existence of off-equatorial tori are investigated. It is proved that these levitating systems are admitted under general conditions when both gravitational and magnetic fields contribute to shaping the spatial profiles of equilibrium plasma fluid fields. Then, specifically, kinetic effects carried by the equilibrium solution are explicitly provided and identified here with diamagnetic energy-correction and electrostatic contributions. It is shown that these kinetic terms characterize the plasma equation of state by introducing non-vanishing deviations from the assumption of thermal pressure.
Method of description of the Alfven and magnetosonic branches of inhomogeneous plasma oscillations
Klimushkin, D.Yu.
1994-12-31
A method of description of the Alfven and magnetosonic branches of the 3-D inhomogeneous plasma oscillations is proposed. In the absence of kinetic effects, the electric field of the MHD-wave is two-dimensional and can be split into potential and eddy components. The first component is identified with the Alfven wave, and the second one with the FMS-wave. It was shown that this approach can be applied to analyze the Earth`s magnetosphere oscillations (geomagnetic pulsations). 21 refs.
Kinetic axisymmetric gravitational equilibria in collisionless accretion disk plasmas
Cremaschini, Claudio; Miller, John C.; Tessarotto, Massimo
2010-07-15
A theoretical treatment is presented of kinetic equilibria in accretion disks (AD) around compact objects, for cases where the plasma can be considered as collisionless. The plasma is assumed to be axisymmetric and to be acted on by gravitational and electromagnetic fields; in this paper, the particular case is considered where the magnetic field admits a family of toroidal magnetic surfaces, which are locally mutually nested and closed. It is pointed out that there exist asymptotic kinetic equilibria represented by generalized bi-Maxwellian distribution functions and characterized by primarily toroidal differential rotation and temperature anisotropy. It is conjectured that kinetic equilibria of this type can exist which are able to sustain both toroidal and poloidal electric current densities, the latter being produced via finite Larmor-radius effects associated with the temperature anisotropy. This leads to the possibility of existence of a new kinetic effect - referred to here as a 'kinetic dynamo effect - resulting in the self-generation of toroidal magnetic field even by a stationary plasma, without any net radial accretion flow being required. The conditions for these equilibria to occur, their basic theoretical features, and their physical properties are all discussed in detail.
Emergence of kinetic behavior in streaming ultracold neutral plasmas
McQuillen, P.; Castro, J.; Bradshaw, S. J.; Killian, T. C.
2015-04-15
We create streaming ultracold neutral plasmas by tailoring the photoionizing laser beam that creates the plasma. By varying the electron temperature, we control the relative velocity of the streaming populations, and, in conjunction with variation of the plasma density, this controls the ion collisionality of the colliding streams. Laser-induced fluorescence is used to map the spatially resolved density and velocity distribution function for the ions. We identify the lack of local thermal equilibrium and distinct populations of interpenetrating, counter-streaming ions as signatures of kinetic behavior. Experimental data are compared with results from a one-dimensional, two-fluid numerical simulation.
Kinetic Description of Vacuum Creation of Massive Vector Bosons
Blaschke, D.B.; Prozorkevich, A.V.; Smolyansky, S.A.; Reichel, A.V.
2005-06-01
In the simple model of massive vector field in a flat spacetime, we derive the kinetic equation of non-Markovian type describing the vacuum pair creation under action of external fields of different nature. We use for this aim the nonperturbative methods of kinetic theory in combination with a new element when the transition of the instantaneous quasiparticle representation is realized within the oscillator (holomorphic) representation. We study in detail the process of vacuum creation of vector bosons generated by a time-dependent boson mass in accordance with the framework of a conformal-invariant scalar-tensor gravitational theory and its cosmological application. It is indicated that the choice of the equation of state allows one to obtain a number density of vector bosons that is sufficient to explain the observed number density of photons in the cosmic microwave background radiation.
Kinetic theory of nonlinear transport phenomena in complex plasmas
Mishra, S. K.; Sodha, M. S.
2013-03-15
In contrast to the prevalent use of the phenomenological theory of transport phenomena, a number of transport properties of complex plasmas have been evaluated by using appropriate expressions, available from the kinetic theory, which are based on Boltzmann's transfer equation; in particular, the energy dependence of the electron collision frequency has been taken into account. Following the recent trend, the number and energy balance of all the constituents of the complex plasma and the charge balance on the particles is accounted for; the Ohmic loss has also been included in the energy balance of the electrons. The charging kinetics for the complex plasma comprising of uniformly dispersed dust particles, characterized by (i) uniform size and (ii) the Mathis, Rumpl, and Nordsieck power law of size distribution has been developed. Using appropriate expressions for the transport parameters based on the kinetic theory, the system of equations has been solved to investigate the parametric dependence of the complex plasma transport properties on the applied electric field and other plasma parameters; the results are graphically illustrated.
Igor D. Kaganovich; Oleg Polomarov
2003-05-19
In low-pressure discharges, when the electron mean free path is larger or comparable with the discharge length, the electron dynamics is essentially non-local. Moreover, the electron energy distribution function (EEDF) deviates considerably from a Maxwellian. Therefore, an accurate kinetic description of the low-pressure discharges requires knowledge of the non-local conductivity operator and calculation of the non-Maxwellian EEDF. The previous treatments made use of simplifying assumptions: a uniform density profile and a Maxwellian EEDF. In the present study a self-consistent system of equations for the kinetic description of nonlocal, non-uniform, nearly collisionless plasmas of low-pressure discharges is derived. It consists of the nonlocal conductivity operator and the averaged kinetic equation for calculation of the non-Maxwellian EEDF. The importance of accounting for the non-uniform plasma density profile on both the current density profile and the EEDF is demonstrated.
Non Equilbrium Vibrational Kinetics in Expanding Plasma Flows
Colonna, Gianpiero
2008-12-31
The supersonic expansion of a plasma is a system of interest for aerospace applications, ranging from propulsion to hypersonic wind tunnels. Under these conditions the plasma shows significant departures from chemical and thermal equilibrium, similarly to post-discharge conditions. The multitemperature description is not adequate because the internal level distributions show tails overpopulated with respect to a Boltzmann distribution. The state-to-state approach has to be used, including the interaction with free electrons which follow non-maxwellian distributions.
Kinetic instabilities in a mirror-confined ECR discharge plasma
NASA Astrophysics Data System (ADS)
Mansfeld, Dmitry; Viktorov, Mikhail; Vodopyanov, Alexander; Golubev, Sergey
2015-11-01
Kinetic instabilities of nonequilibrium plasma heated by powerful radiation of gyrotron in electron cyclotron resonance conditions and confined in a mirror magnetic trap are reported. Instabilities are manifested as the generation of short pulses of electromagnetic radiation accompanied by precipitation of hot electrons from magnetic trap. Measuring electromagnetic field with high temporal resolution allowed to observe various dynamic spectra of electromagnetic radiation related to at least five types of kinetic instabilities. The opportunity to recreate different conditions for excitation and amplification of waves in plasma in a single ECR discharge pulse has been demonstrated. This report may be of interest in the context of a laboratory modeling of nonstationary wave-particle interaction processes in nonequilibrium space plasma since the observed phenomena have much in common with similar processes occurring in the magnetosphere of the Earth, planets, and in solar coronal loops. Work was supported by Russian Foundation for Basic Research # 15-32-20770.
A spectral Poisson solver for kinetic plasma simulation
NASA Astrophysics Data System (ADS)
Szeremley, Daniel; Obberath, Jens; Brinkmann, Ralf
2011-10-01
Plasma resonance spectroscopy is a well established plasma diagnostic method, realized in several designs. One of these designs is the multipole resonance probe (MRP). In its idealized - geometrically simplified - version it consists of two dielectrically shielded, hemispherical electrodes to which an RF signal is applied. A numerical tool is under development which is capable of simulating the dynamics of the plasma surrounding the MRP in electrostatic approximation. In this contribution we concentrate on the specialized Poisson solver for that tool. The plasma is represented by an ensemble of point charges. By expanding both the charge density and the potential into spherical harmonics, a largely analytical solution of the Poisson problem can be employed. For a practical implementation, the expansion must be appropriately truncated. With this spectral solver we are able to efficiently solve the Poisson equation in a kinetic plasma simulation without the need of introducing a spatial discretization.
Electron kinetics in a cooling plasma
Helander, P.; Smith, H.; Fueloep, T.; Eriksson, L.-G.
2004-12-01
The distribution function of suprathermal electrons in a slowly cooling plasma is calculated by an asymptotic expansion in the cooling rate divided by the collision frequency. Since the collision frequency decreases with increasing velocity, a high-energy tail forms in the electron distribution function as the bulk population cools down. Under certain simplifying assumptions (slow cooling, constant density, Born approximation of cross sections), the distribution function evolves to a self-similar state where the tail is inversely proportional to the cube of the velocity. Its practical consequences are discussed briefly.
NASA Technical Reports Server (NTRS)
Roth, J. R.
1976-01-01
Parametric variation of independent variables which may affect the characteristics of the NASA Lewis Bumpy Torus plasma have identified those which have a significant effect on the plasma current, ion kinetic temperature, and plasma number density, and those which do not. Empirical power-law correlations of the plasma current, and the ion kinetic temperature and number density were obtained as functions of the potential applied to the midplane electrode rings, the background neutral gas pressure, and the magnetic field strength. Additional parameters studied include the type of gas, the polarity of the midplane electrode rings (and hence the direction of the radial electric field), the mode of plasma operation, and the method of measuring the plasma number density. No significant departures from the scaling laws appear to occur at the highest ion kinetic temperatures or number densities obtained to date.
NASA Technical Reports Server (NTRS)
Roth, J. R.
1976-01-01
Parametric variation of independent variables which may affect the characteristics of bumpy torus plasma have identified those which have a significant effect on the plasma current, ion kinetic temperature, and plasma number density, and those which do not. Empirical power law correlations of the plasma current, and the ion kinetic temperature and number density were obtained as functions of potential applied to the midplane electrode rings, the background neutral gas pressure, and the magnetic field strength. Additional parameters studied included the type of gas, the polarity of the midplane electrode rings, the mode of plasma operation, and the method of measuring the plasma number density. No significant departures from the scaling laws appear to occur at the highest ion kinetic temperatures or number densities obtained to date.
Toward Petaflop First Principles Kinetic Plasma Simulation
NASA Astrophysics Data System (ADS)
Bowers, Kevin
2007-11-01
Due to physical limitations (such as the speed of light), moving data between and even within modern microprocessors is more time consuming than performing computations. As a result, individual processor core performance is stagnant, multicore processors are ubiquitous and traditional programming styles are unable to exploit the potential of modern computers fully. This talk will discuss the architecture and implementation of the 3d electromagnetic relativistic particle-in-cell code VPIC for LANL's Roadrunner supercomputer. Roadrunner is expected to have 13,000 IBM Cell microprocessors (each Cell contains a dual threaded Power core and 8 specialized vector cores) and be capable of over a petaflop (10^15 floating point operations per second). VPIC minimizes data movement and allows vector extensions of modern processors to be utilized portably. This made it possible to port VPIC quickly while achieving unprecedented performance. The initial port performed 0.13 billion particles pushed and accumulated per second per Cell---equivalent to 1.0 billion per second per 8 Cell node or sustaining Roadrunner at 0.4 petaflop. Higher performance is likely as the port is refined. Regardless, already demonstrated performance will enable previously intractable simulations in numerous areas of plasma physics, including magnetic reconnection and laser plasma interactions.
Quantitative description of realistic wealth distributions by kinetic trading models
NASA Astrophysics Data System (ADS)
Lammoglia, Nelson; Muñoz, Víctor; Rogan, José; Toledo, Benjamín; Zarama, Roberto; Valdivia, Juan Alejandro
2008-10-01
Data on wealth distributions in trading markets show a power law behavior x-(1+α) at the high end, where, in general, α is greater than 1 (Pareto’s law). Models based on kinetic theory, where a set of interacting agents trade money, yield power law tails if agents are assigned a saving propensity. In this paper we are solving the inverse problem, that is, in finding the saving propensity distribution which yields a given wealth distribution for all wealth ranges. This is done explicitly for two recently published and comprehensive wealth datasets.
Quantitative description of realistic wealth distributions by kinetic trading models.
Lammoglia, Nelson; Muñoz, Víctor; Rogan, José; Toledo, Benjamín; Zarama, Roberto; Valdivia, Juan Alejandro
2008-10-01
Data on wealth distributions in trading markets show a power law behavior x(-)(1+alpha) at the high end, where, in general, alpha is greater than 1 (Pareto's law). Models based on kinetic theory, where a set of interacting agents trade money, yield power law tails if agents are assigned a saving propensity. In this paper we are solving the inverse problem, that is, in finding the saving propensity distribution which yields a given wealth distribution for all wealth ranges. This is done explicitly for two recently published and comprehensive wealth datasets. PMID:18999570
Cremaschini, Claudio Stuchlík, Zdeněk; Tessarotto, Massimo
2014-03-15
The kinetic description of relativistic plasmas in the presence of time-varying and spatially non-uniform electromagnetic (EM) fields is a fundamental theoretical issue both in astrophysics and plasma physics. This refers, in particular, to the treatment of collisionless and strongly-magnetized plasmas in the presence of intense radiation sources. In this paper, the problem is investigated in the framework of a covariant gyrokinetic treatment for Vlasov–Maxwell equilibria. The existence of a new class of kinetic equilibria is pointed out, which occur for spatially-symmetric systems. These equilibria are shown to exist in the presence of non-uniform background EM fields and curved space-time. In the non-relativistic limit, this feature permits the determination of kinetic equilibria even for plasmas in which particle energy is not conserved due to the occurrence of explicitly time-dependent EM fields. Finally, absolute stability criteria are established which apply in the case of infinitesimal symmetric perturbations that can be either externally or internally produced.
Turbulence and Proton–Electron Heating in Kinetic Plasma
NASA Astrophysics Data System (ADS)
Matthaeus, William H.; Parashar, Tulasi N.; Wan, Minping; Wu, P.
2016-08-01
Analysis of particle-in-cell simulations of kinetic plasma turbulence reveals a connection between the strength of cascade, the total heating rate, and the partitioning of dissipated energy into proton heating and electron heating. A von Karman scaling of the cascade rate explains the total heating across several families of simulations. The proton to electron heating ratio increases in proportion to total heating. We argue that the ratio of gyroperiod to nonlinear turnover time at the ion kinetic scales controls the ratio of proton and electron heating. The proposed scaling is consistent with simulations.
Charging kinetics of dust in interplanetary space plasma
NASA Astrophysics Data System (ADS)
Misra, Shikha; Mishra, S. K.
2013-07-01
A theoretical kinetic model for the physical understanding of the charging of dust particles in the interplanetary space plasma has been developed. In contrast to earlier studies, the present analysis incorporates (i) uniform potential theory for complex plasmas with size distribution of the dust particles, (ii) charge, number and energy balance of the constituents and (iii) appropriate expressions for photoelectric emission from a positively charged particle with inherent charge neutrality of the interplanetary space plasma. Further utilizing the population balance equation (given by Matsoukas and Russel) for the interplanetary dust particles, the fluctuations in steady-state charge (or electric potential) has also been investigated. For the illustration purpose, the computations have been performed for the interplanetary space plasma at 1 au from the sun; for this distance, reasonably good information on the gaseous and dust components are available. As an interesting feature, the theoretical predictions are in reasonably good agreement with observations and earlier estimates.
Wave-kinetic description of finite temperature Bose–Einstein condensates
NASA Astrophysics Data System (ADS)
Mendonça, J. T.
2016-07-01
We derive a system of two coupled wave-kinetic equations for the condensate and for the thermal gas. This provides a simple and elegant model, based on a two fluid interaction, in the spirit of the traditional theory of superfluidity. One of the fluids is the condensed gas, and the other is the fluid of bogolon quasi-particles. As an illustration, we apply this wave-kinetic model to the description of the second-sound.
Cremaschini, Claudio; Stuchlík, Zdeněk; Tessarotto, Massimo
2013-05-15
The problem of formulating a kinetic treatment for quasi-stationary collisionless plasmas in axisymmetric systems subject to the possibly independent presence of local strong velocity-shear and supersonic rotation velocities is posed. The theory is developed in the framework of the Vlasov-Maxwell description for multi-species non-relativistic plasmas. Applications to astrophysical accretion discs arising around compact objects and to plasmas in laboratory devices are considered. Explicit solutions for the equilibrium kinetic distribution function (KDF) are constructed based on the identification of the relevant particle adiabatic invariants. These are shown to be expressed in terms of generalized non-isotropic Gaussian distributions. A suitable perturbative theory is then developed which allows for the treatment of non-uniform strong velocity-shear/supersonic plasmas. This yields a series representation for the equilibrium KDF in which the leading-order term depends on both a finite set of fluid fields as well as on the gradients of an appropriate rotational frequency. Constitutive equations for the fluid number density, flow velocity, and pressure tensor are explicitly calculated. As a notable outcome, the discovery of a new mechanism for generating temperature and pressure anisotropies is pointed out, which represents a characteristic feature of plasmas considered here. This is shown to arise as a consequence of the canonical momentum conservation and to contribute to the occurrence of temperature anisotropy in combination with the adiabatic conservation of the particle magnetic moment. The physical relevance of the result and the implications of the kinetic solution for the self-generation of quasi-stationary electrostatic and magnetic fields through a kinetic dynamo are discussed.
Weakly Ionized Plasmas in Hypersonics: Fundamental Kinetics and Flight Applications
Macheret, Sergey
2005-05-16
The paper reviews some of the recent studies of applications of weakly ionized plasmas to supersonic/hypersonic flight. Plasmas can be used simply as means of delivering energy (heating) to the flow, and also for electromagnetic flow control and magnetohydrodynamic (MHD) power generation. Plasma and MHD control can be especially effective in transient off-design flight regimes. In cold air flow, nonequilibrium plasmas must be created, and the ionization power budget determines design, performance envelope, and the very practicality of plasma/MHD devices. The minimum power budget is provided by electron beams and repetitive high-voltage nanosecond pulses, and the paper describes theoretical and computational modeling of plasmas created by the beams and repetitive pulses. The models include coupled equations for non-local and unsteady electron energy distribution function (modeled in forward-back approximation), plasma kinetics, and electric field. Recent experimental studies at Princeton University have successfully demonstrated stable diffuse plasmas sustained by repetitive nanosecond pulses in supersonic air flow, and for the first time have demonstrated the existence of MHD effects in such plasmas. Cold-air hypersonic MHD devices are shown to permit optimization of scramjet inlets at Mach numbers higher than the design value, while operating in self-powered regime. Plasma energy addition upstream of the inlet throat can increase the thrust by capturing more air (Virtual Cowl), or it can reduce the flow Mach number and thus eliminate the need for an isolator duct. In the latter two cases, the power that needs to be supplied to the plasma would be generated by an MHD generator downstream of the combustor, thus forming the 'reverse energy bypass' scheme. MHD power generation on board reentry vehicles is also discussed.
Weakly Ionized Plasmas in Hypersonics: Fundamental Kinetics and Flight Applications
NASA Astrophysics Data System (ADS)
Macheret, Sergey
2005-05-01
The paper reviews some of the recent studies of applications of weakly ionized plasmas to supersonic/hypersonic flight. Plasmas can be used simply as means of delivering energy (heating) to the flow, and also for electromagnetic flow control and magnetohydrodynamic (MHD) power generation. Plasma and MHD control can be especially effective in transient off-design flight regimes. In cold air flow, nonequilibrium plasmas must be created, and the ionization power budget determines design, performance envelope, and the very practicality of plasma/MHD devices. The minimum power budget is provided by electron beams and repetitive high-voltage nanosecond pulses, and the paper describes theoretical and computational modeling of plasmas created by the beams and repetitive pulses. The models include coupled equations for non-local and unsteady electron energy distribution function (modeled in forward-back approximation), plasma kinetics, and electric field. Recent experimental studies at Princeton University have successfully demonstrated stable diffuse plasmas sustained by repetitive nanosecond pulses in supersonic air flow, and for the first time have demonstrated the existence of MHD effects in such plasmas. Cold-air hypersonic MHD devices are shown to permit optimization of scramjet inlets at Mach numbers higher than the design value, while operating in self-powered regime. Plasma energy addition upstream of the inlet throat can increase the thrust by capturing more air (Virtual Cowl), or it can reduce the flow Mach number and thus eliminate the need for an isolator duct. In the latter two cases, the power that needs to be supplied to the plasma would be generated by an MHD generator downstream of the combustor, thus forming the "reverse energy bypass" scheme. MHD power generation on board reentry vehicles is also discussed.
Kinetic theory of Jeans instability of a dusty plasma.
Pandey, B P; Lakhina, G S; Krishan, V
1999-12-01
A kinetic theory of the Jeans instability of a dusty plasma has been developed in the present work. The effect of grain charge fluctuations due to the attachment of electrons and ions to the grain surface has been considered in the framework of Krook's collisional model. We demonstrate that the grain charge fluctuations alter the growth rate of the gravitational collapse of the dusty plasma. The Jeans length has been derived under limiting cases, and its dependence on the attachment frequency is shown. In the absence of gravity, we see that the damping rate of the dust acoustic mode is proportional to the electron-dust collision frequency. PMID:11970688
Species separation and kinetic effects in collisional plasma shocks
Bellei, C. Wilks, S. C.; Amendt, P. A.; Rinderknecht, H.; Zylstra, A.; Rosenberg, M.; Sio, H.; Li, C. K.; Petrasso, R.
2014-05-15
The properties of collisional shock waves propagating in uniform plasmas are studied with ion-kinetic calculations, in both slab and spherical geometry and for the case of one and two ion species. Despite the presence of an electric field at the shock front—and in contrast to the case where an interface is initially present [C. Bellei et al., Phys. Plasmas 20, 044702 (2013)]—essentially no ion reflection at the shock front is observed due to collisions, with a probability of reflection ≲10{sup −4} for the cases presented. A kinetic two-ion-species spherical convergent shock is studied in detail and compared against an average-species calculation, confirming effects of species separation and differential heating of the ion species at the shock front. The effect of different ion temperatures on the DT and D{sup 3}He fusion reactivity is discussed in the fluid limit and is estimated to be moderately important.
Fluid description of multi-component solar partially ionized plasma
Khomenko, E. Collados, M.; Vitas, N.; Díaz, A.
2014-09-15
We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed, we particularize to some frequently considered cases as for the interaction of matter and radiation.
Transition of electron kinetics in weakly magnetized inductively coupled plasmas
Kim, Jin-Yong; Lee, Hyo-Chang; Kim, Young-Do; Chung, Chin-Wook; Kim, Young-Cheol
2013-10-15
Transition of the electron kinetics from nonlocal to local regime was studied in weakly magnetized solenoidal inductively coupled plasma from the measurement of the electron energy probability function (EEPF). Without DC magnetic field, the discharge property was governed by nonlocal electron kinetics at low gas pressure. The electron temperatures were almost same in radial position, and the EEPFs in total electron energy scale were radially coincided. However, when the DC magnetic field was applied, radial non-coincidence of the EEPFs in total electron energy scale was observed. The electrons were cooled at the discharge center where the electron heating is absent, while the electron temperature was rarely changed at the discharge boundary with the magnetic field. These changes show the transition from nonlocal to local electron kinetics and the transition is occurred when the electron gyration diameter was smaller than the skin depth. The nonlocal to local transition point almost coincided with the calculation results by using nonlocal parameter and collision parameter.
On bias of kinetic temperature measurements in complex plasmas
Kantor, M.; Association Euratom-FOM Institute DIFFER, 3430 BE Nieuwegein; Ioffe Institute, RAS, St. Petersburg 194021 ; Moseev, D.; Salewski, M.
2014-02-15
The kinetic temperature in complex plasmas is often measured using particle tracking velocimetry. Here, we introduce a criterion which minimizes the probability of faulty tracking of particles with normally distributed random displacements in consecutive frames. Faulty particle tracking results in a measurement bias of the deduced velocity distribution function and hence the deduced kinetic temperature. For particles with a normal velocity distribution function, mistracking biases the obtained velocity distribution function towards small velocities at the expense of large velocities, i.e., the inferred velocity distribution is more peaked and its tail is less pronounced. The kinetic temperature is therefore systematically underestimated in measurements. We give a prescription to mitigate this type of error.
A First-Principle Kinetic Theory of Meteor Plasma Formation
NASA Astrophysics Data System (ADS)
Dimant, Yakov; Oppenheim, Meers
2015-11-01
Every second millions of tiny meteoroids hit the Earth from space, vast majority too small to observe visually. However, radars detect the plasma they generate and use the collected data to characterize the incoming meteoroids and the atmosphere in which they disintegrate. This diagnostics requires a detailed quantitative understanding of formation of the meteor plasma. Fast-descending meteoroids become detectable to radars after they heat due to collisions with atmospheric molecules sufficiently and start ablating. The ablated material then collides into atmospheric molecules and forms plasma around the meteoroid. Reflection of radar pulses from this plasma produces a localized signal called a head echo. Using first principles, we have developed a consistent collisional kinetic theory of the near-meteoroid plasma. This theory shows that the meteoroid plasma develops over a length-scale close to the ion mean free path with a non-Maxwellian velocity distribution. The spatial distribution of the plasma density shows significant deviations from a Gaussian law usually employed in head-echo modeling. This analytical model will serve as a basis for more accurate quantitative interpretation of the head echo radar measurements. Work supported by NSF Grant 1244842.
Kinetic modeling of the Saturn ring-ionosphere plasma environment
NASA Technical Reports Server (NTRS)
Wilson, G. R.; Waite, J. H., Jr.
1989-01-01
A time-independent kinetic plasma model was developed on the basis of the Li et al. (1988) semikinetic plasma model and was used to study the interaction of the Saturnian ionosphere and ring plasma. The model includes the gravitational magnetic mirror and centripetal and ambipolar electric forces, and the effect of the mixing of two plasma populations. The results obtained indicate that the density, temperature, and composition of plasma near the rings changing in the direction from the inner C ring to the outer A ring, due to the fact that the predominant source of plasma changes from the ionosphere to the rings. The model results also suggest that the outflow of hydrogen from the ionosphere to the rings may be shut off for field lines passing through the outer B and A ring, due to the ambipolar electric field set up by the warm ring plasma trapped near the ring plane by the centipetal force. In these regions, there will be a net flux of O(+) ions from the rings to the ionosphere.
Plasma transport induced by kinetic Alfven wave turbulence
Izutsu, T.; Hasegawa, H.; Fujimoto, M.; Nakamura, T. K. M.
2012-10-15
At the Earth's magnetopause that separates the hot-tenuous magnetospheric plasma from the cold dense solar wind plasma, often seen is a boundary layer where plasmas of both origins coexist. Plasma diffusions of various forms have been considered as the cause of this plasma mixing. Here, we investigate the plasma transport induced by wave-particle interaction in kinetic Alfven wave (KAW) turbulence, which is one of the candidate processes. We clarify that the physical origin of the KAW-induced cross-field diffusion is the drift motions of those particles that are in Cerenkov resonance with the wave: E Multiplication-Sign B-like drift that emerges in the presence of non-zero parallel electric field component and grad-B drift due to compressional magnetic fluctuations. We find that KAW turbulence, which has a spectral breakpoint at which an MHD inertial range transits to a dissipation range, causes selective transport for particles whose parallel velocities are specified by the local Alfven velocity and the parallel phase velocity at the spectral breakpoint. This finding leads us to propose a new data analysis method for identifying whether or not a mixed plasma in the boundary layer is a consequence of KAW-induced transport across the magnetopause. The method refers to the velocity space distribution function data obtained by a spacecraft that performs in situ observations and, in principle, is applicable to currently available dataset such as that provided by the NASA's THEMIS mission.
Unification of Plasma Fluid and Kinetic Theory via Gaussian Radial Basis Functions
NASA Astrophysics Data System (ADS)
Candy, J. M.
2015-11-01
A fundamental macroscopic description of a magnetized plasma is the Vlasov equation supplemented by the nonlinear inverse-square force Fokker-Planck collision operator [Rosenbluth et al., Phys. Rev. 107, 1957]. The Vlasov part describes advection in a six-dimensional phase space whereas the collision operator contains friction and diffusion coefficients that are weighted velocity-space integrals of the particle distribution function. The Fokker-Planck collision operator is an integro-differential, nonlinear (bilinear) operator. Numerical discretization of the operator, in particular for collisions of unlike species, is extremely challenging. In this work, we describe a new approach to discretize the entire kinetic system based on an expansion in Gaussian Radial Basis functions (RBFs). This approach is particularly well-suited to treat the collision operator because the friction and diffusion coefficients can be analytically calculated. Although the RBF method is known to be a powerful scheme for the interpolation of scattered multidimensional data, Gaussian RBFs also have a deep physical interpretation in statistical mechanics and plasma physics as local thermodynamic equilibria. We outline the general theory, highlight the connection to plasma fluid theories, and also give 2D and 3D numerical solutions of the nonlinear Fokker-Planck equation. A broad spectrum of applications for the new method is anticipated in both astrophysical and laboratory plasmas. In particular, we believe that the RBF method may provide a new bridge between fluid and kinetic descriptions of magnetized plasma. Work supported in part by US DOE under DE-FG02-08ER54963.
Unified description of linear screening in dense plasmas.
Stanton, L G; Murillo, M S
2015-03-01
Electron screening of ions is among the most fundamental properties of plasmas, determining the effective ionic interactions that impact all properties of a plasma. With the development of new experimental facilities that probe high-energy-density physics regimes ranging from warm dense matter to hot dense matter, a unified framework for describing dense plasma screening has become essential. Such a unified framework is presented here based on finite-temperature orbital-free density functional theory, including gradient corrections and exchange-correlation effects. We find a new analytic pair potential for the ion-ion interaction that incorporates moderate electronic coupling, quantum degeneracy, gradient corrections to the free energy, and finite temperatures. This potential can be used in large-scale "classical" molecular dynamics simulations, as well as in simpler theoretical models (e.g., integral equations and Monte Carlo), with no additional computational complexity. The new potential theoretically connects limits of Debye-Hückel-Yukawa, Lindhard, Thomas-Fermi, and Bohmian quantum hydrodynamics descriptions. Based on this new potential, we predict ionic static structure factors that can be validated using x-ray Thomson scattering data. PMID:25871221
NASA Technical Reports Server (NTRS)
Manning, Robert M.
2009-01-01
Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically induced plasma, it has been demonstrated that the classical radiofrequency communications blackout that is experienced during atmospheric reentry can be mitigated through the appropriate control of an external magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and involves an analytical solution for the electric and magnetic fields within the plasma allowing for a description of the attendant transmission, reflection and absorption coefficients. The ability to transmit through the magnetized plasma is due to the magnetic windows that are created within the plasma via the well-known whistler modes of propagation. The case of 2 GHz transmission through a re-entry plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves through the re-entry phase.
Phase-space description of plasma waves: Linear and nonlinear theory
NASA Astrophysics Data System (ADS)
Biro, Thomas
1992-11-01
A (r,k) phase description of waves in plasmas is developed by introducing Gaussian window functions to separate short scale oscillations from long scale modulations of the wave fields and variations in the plasma parameters. To obtain a wave equation that unambiguously separates conservative dynamics from dissipation also in an inhomogeneous and time varying background plasma, the proper form of the current response function, is discussed. On the analogy of the particle distribution function f(v,r,t), a wave density N(k,r,t) is introduced on phase space. This function is proven to satisfy a simple continuity equation. Dissipation is also included, and this allows the damping or growth of wave density along rays to be described. Problems involving geometric optics of continuous media often appear simpler when viewed in phase space, since the flow of N in phase space is incompressible. Within the phase space representation, a very general formula for the second order nonlinear current is obtained in terms of the vector potential. This formula is a convenient starting point for studies of coherent as well as turbulent nonlinear processes. Kinetic equations for weakly inhomogeneous and turbulent plasmas are derived, including the effects of inhomogeneous turbulence, wave convection and refraction.
PDRK: A General Kinetic Dispersion Relation Solver for Magnetized Plasma
NASA Astrophysics Data System (ADS)
Xie, Huasheng; Xiao, Yong
2016-02-01
A general, fast, and effective approach is developed for numerical calculation of kinetic plasma linear dispersion relations. The plasma dispersion function is approximated by J-pole expansion. Subsequently, the dispersion relation is transformed to a standard matrix eigenvalue problem of an equivalent linear system. Numerical solutions for the least damped or fastest growing modes using an 8-pole expansion are generally accurate; more strongly damped modes are less accurate, but are less likely to be of physical interest. In contrast to conventional approaches, such as Newton's iterative method, this approach can give either all the solutions in the system or a few solutions around the initial guess. It is also free from convergence problems. The approach is demonstrated for electrostatic dispersion equations with one-dimensional and two-dimensional wavevectors, and for electromagnetic kinetic magnetized plasma dispersion relation for bi-Maxwellian distribution with relative parallel velocity flows between species. supported by the National Magnetic Confinement Fusion Science Program of China (Nos. 2015GB110003, 2011GB105001, 2013GB111000), National Natural Science Foundation of China (No. 91130031), the Recruitment Program of Global Youth Experts
Kinetic theory of plasma adiabatic major radius compression in tokamaks
NASA Astrophysics Data System (ADS)
Gorelenkova, M. V.; Gorelenkov, N. N.; Azizov, E. A.; Romannikov, A. N.; Herrmann, H. W.
1998-05-01
In order to understand the individual charged particle behavior as well as plasma macroparameters (temperature, density, etc.) during the adiabatic major radius compression (R-compression) in a tokamak, a kinetic approach is used. The perpendicular electric field from the Ohm's law at zero resistivity is made use of in order to describe particle motion during the R-compression. Expressions for both passing and trapped particle energy and pitch angle change are derived for a plasma with high aspect ratio and circular magnetic surfaces. The particle behavior near the passing trapped boundary during the compression is studied to simulate the compression-induced collisional losses of alpha particles. Qualitative agreement is obtained with the alphas loss measurements in deuterium-tritium (D-T) experiments in the Tokamak Fusion Test Reactor (TFTR) [World Survey of Activities in Controlled Fusion Research [Nucl. Fusion special supplement (1991)] (International Atomic Energy Agency, Vienna, 1991)]. The plasma macroparameters evolution at the R-compression is calculated by solving the gyroaveraged drift kinetic equation.
Kinetic phenomena in charged particle transport in gases and plasmas
Petrovic, Zoran Lj.; Dujko, Sasa; Sasic, Olivera; Stojanovic, Vladimir; Malovic, Gordana
2012-05-25
The key difference between equilibrium (thermal) and non-equilibrium (low temperature - a.k.a. cold) plasmas is in the degree in which the shape of the cross sections influences the electron energy distribution function (EEDF). In this paper we will discuss the issue of kinetic phenomena from two different angles. The first will be how to take advantage of the strong influence and use low current data to obtain the cross sections. This is also known as the swarm technique and the product of a ''swarm analysis'' is a set of cross sections giving good number, momentum and energy balances of electrons or other charged particles. At the same time understanding the EEDF is based on the cross section data. Nevertheless sometimes the knowledge of the cross sections and even the behaviour of individual particles are insufficient to explain collective behaviour of the ensemble. The resulting ''kinetic'' effects may be used to favour certain properties of non-equilibrium plasmas and even may be used as the basis of some new plasma applications.
Kinetic Effects in Low Pressure Capacitively Coupled Plasmas
NASA Astrophysics Data System (ADS)
Likhanskii, Alexandre; Roark, Christine; Stoltz, Peter
2011-10-01
We present results of particle-in-cell/Monte Carlo collision simulations of kinetic effects in low pressure capacitively coupled plasma discharge. In particular, we examine discharges of various gases (including Ar, Xe, and others) in the pressure range of 10s of mT and the frequency range of 10s of MHz. We track the formation of high energy electrons (e.g., at the ionization threshold or greater) as a marker for enhanced ionization, and look at the effects of elastic and inelastic collisions on the formation of these high energy electron bunches. We show results for 2D and 3D simulations where we include density gradient effects, and results for plasma chemistry effects on the bulk electron energy distribution function and the ion energy distribution function at a plasma surface interface. We discuss the role of the bunches on electron heating in the plasma bulk and on their presence on how electron heating is treated in fluid simulations of plasma sources.
Propagation of radiation in fluctuating multiscale plasmas. II. Kinetic simulations
Pal Singh, Kunwar; Robinson, P. A.; Cairns, Iver H.; Tyshetskiy, Yu.
2012-11-15
A numerical algorithm is developed and tested that implements the kinetic treatment of electromagnetic radiation propagating through plasmas whose properties have small scale fluctuations, which was developed in a companion paper. This method incorporates the effects of refraction, damping, mode structure, and other aspects of large-scale propagation of electromagnetic waves on the distribution function of quanta in position and wave vector, with small-scale effects of nonuniformities, including scattering and mode conversion approximated as causing drift and diffusion in wave vector. Numerical solution of the kinetic equation yields the distribution function of radiation quanta in space, time, and wave vector. Simulations verify the convergence, accuracy, and speed of the methods used to treat each term in the equation. The simulations also illustrate the main physical effects and place the results in a form that can be used in future applications.
Kinetic simulations of magnetized turbulence in astrophysical plasmas.
Howes, G G; Dorland, W; Cowley, S C; Hammett, G W; Quataert, E; Schekochihin, A A; Tatsuno, T
2008-02-15
This Letter presents the first ab initio, fully electromagnetic, kinetic simulations of magnetized turbulence in a homogeneous, weakly collisional plasma at the scale of the ion Larmor radius (ion gyroscale). Magnetic- and electric-field energy spectra show a break at the ion gyroscale; the spectral slopes are consistent with scaling predictions for critically balanced turbulence of Alfvén waves above the ion gyroscale (spectral index -5/3) and of kinetic Alfvén waves below the ion gyroscale (spectral indices of -7/3 for magnetic and -1/3 for electric fluctuations). This behavior is also qualitatively consistent with in situ measurements of turbulence in the solar wind. Our findings support the hypothesis that the frequencies of turbulent fluctuations in the solar wind remain well below the ion cyclotron frequency both above and below the ion gyroscale. PMID:18352484
NASA Astrophysics Data System (ADS)
Chen, Qiang; Chen, Bin
2012-10-01
In this paper, a hybrid electrodynamics and kinetics numerical model based on the finite-difference time-domain method and lattice Boltzmann method is presented for electromagnetic wave propagation in weakly ionized hydrogen plasmas. In this framework, the multicomponent Bhatnagar-Gross-Krook collision model considering both elastic and Coulomb collisions and the multicomponent force model based on the Guo model are introduced, which supply a hyperfine description on the interaction between electromagnetic wave and weakly ionized plasma. Cubic spline interpolation and mean filtering technique are separately introduced to solve the multiscalar problem and enhance the physical quantities, which are polluted by numerical noise. Several simulations have been implemented to validate our model. The numerical results are consistent with a simplified analytical model, which demonstrates that this model can obtain satisfying numerical solutions successfully.
A linear dispersion relation for the hybrid kinetic-ion/fluid-electron model of plasma physics
NASA Astrophysics Data System (ADS)
Told, D.; Cookmeyer, J.; Astfalk, P.; Jenko, F.
2016-07-01
A dispersion relation for a commonly used hybrid model of plasma physics is developed, which combines fully kinetic ions and a massless-electron fluid description. Although this model and variations of it have been used to describe plasma phenomena for about 40 years, to date there exists no general dispersion relation to describe the linear wave physics contained in the model. Previous efforts along these lines are extended here to retain arbitrary wave propagation angles, temperature anisotropy effects, as well as additional terms in the generalized Ohm’s law which determines the electric field. A numerical solver for the dispersion relation is developed, and linear wave physics is benchmarked against solutions of a full Vlasov–Maxwell dispersion relation solver. This work opens the door to a more accurate interpretation of existing and future wave and turbulence simulations using this type of hybrid model.
Kinetic theory of low-frequency cross-field instability in a weakly ionized plasma. I
Dimant, Y.S.; Sudan, R.N.
1995-04-01
A consistent kinetic theory is developed for the description of electrons under conditions of a low-frequency two-stream {bold E}{times}{bold B} instability in collisionally dominated, weakly ionized plasmas. Starting from the Boltzmann collision integral, a simplified kinetic equation for the electron distribution function has been derived, which takes into account strong pitch-angle scattering of electrons by neutrals, velocity dependence of the electron--neutral collision frequency, etc. Linearized equations describing small oscillations of the electron distribution function and ion density are presented. For the asymptotic case of short waves, the dispersion relation of the {bold E}{times}{bold B} instability has been obtained and analyzed under conditions typical for the lower ionosphere. Under certain conditions, the rigorous kinetic consideration yields substantial changes in results compared to previous theories. The general approach may be applied to other linear and nonlinear low-frequency processes in a weakly ionized plasma. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.
Effects of the g Factor in Semiclassical Kinetic Plasma Theory
Brodin, Gert; Marklund, Mattias; Zamanian, Jens; Ericsson, Aasa; Mana, Piero L.
2008-12-12
A kinetic theory for spin plasmas is put forward, generalizing those of previous authors. In the model, the ordinary phase space is extended to include the spin degrees of freedom. Together with Maxwell's equations, the system is shown to be energy conserving. Analyzing the linear properties, it is found that new types of wave-particle resonances are possible that depend directly on the anomalous magnetic moment of the electron. As a result, new wave modes, not present in the absence of spin, appear. The implications of our results are discussed.
Charged particle dynamics and molecular kinetics in the hydrogen postdischarge plasma
Diomede, P.; Longo, S.; Capitelli, M.
2006-11-15
The afterglow of a parallel plate radio frequency discharge in hydrogen is studied by numerical modelling to compare ion dynamics and chemical effects on the behavior of negative ions. While the ion dynamics requires a kinetic description of space dependent plasma relaxation (at least 1D), chemical effects require a vibrational kinetics of hydrogen molecules. Since previous models did not include both features it has not been possible until now to realize both effects in a single simulation. We apply an updated version of the 1D Bari model which includes a 1.5D (1Dr2Dv) Particle in Cell/Monte Carlo (PIC/MC) multispecies module coupled to the space and time dependent master equation for H{sub 2}(X{sup 1}{sigma}{sub g}{sup +},v=0,...,14) vibrational level population. Negative ion fronts are described in hydrogen for the first time and their impact on the plasma limiting surfaces produces a negative ion current evolution compatible with experimental findings. In the same conditions, the attachment rate overshoot is found to contribute about 7% to the average ion density in the plasma.
NASA Technical Reports Server (NTRS)
Tanaka, Motohiko; Sato, Tetsuya; Hasegawa, A.
1989-01-01
The excitation of the kinetic Alfven wave by resonant mode conversion and longitudinal heating of the plasma by the kinetic Alfven wave were demonstrated on the basis of a macroscale particle simulation. The longitudinal electron current was shown to be cancelled by the ions. The kinetic Alfven wave produced an ordered motion of the plasma particles in the wave propagation direction. The electrons were pushed forward along the ambient magnetic field by absorbing the kinetic Alfven wave through the Landau resonance.
Continuum Kinetic Modeling of the Tokamak Plasma Edge
NASA Astrophysics Data System (ADS)
Dorf, Mikhail
2015-11-01
The problem of edge plasma transport provides substantial challenges for analytical or numerical analysis due to (a) complex magnetic geometry including both open and closed magnetic field lines B, (b) steep radial gradients comparable to ion drift-orbit excursions, and (c) a variation in the collision mean-free path along B from long to short compared to the magnetic connection length. Here, the first 4D continuum drift-kinetic transport simulations that span the magnetic separatrix of a tokamak are presented, motivated in part by the success of continuum kinetic codes for core physics and in part by the potential for high accuracy. The calculations include fully-nonlinear Fokker-Plank collisions and electrostatic potential variations. The problem of intrinsic toroidal rotation driven by ion orbit loss is addressed in detail. The code, COGENT, developed by the Edge Simulation Laboratory collaboration, is distinguished by a fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex magnetic X-point divertor geometry with high accuracy. Previously, successful performance of high-order algorithms has been demonstrated in a simpler closed magnetic-flux-surface geometry for the problems of neoclassical transport and collisionless relaxation of geodesic acoustic modes in a tokamak pedestal, including the effects of a strong radial electric field under H-mode conditions. Work performed for USDOE, at LLNL under contract DE-AC52-07NA27344.
Magnetic Null Points in Kinetic Simulations of Space Plasmas
NASA Astrophysics Data System (ADS)
Olshevsky, Vyacheslav; Deca, Jan; Divin, Andrey; Peng, Ivy Bo; Markidis, Stefano; Innocenti, Maria Elena; Cazzola, Emanuele; Lapenta, Giovanni
2016-03-01
We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic particle-in-cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind, and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly (LMA) and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3-9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and LMA simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.
Kee, R.J.; Rupley, F.M.; Meeks, E.; Miller, J.A.
1996-05-01
This document is the user`s manual for the third-generation CHEMKIN package. CHEMKIN is a software package whose purpose is to facilitate the formation, solution, and interpretation of problems involving elementary gas-phase chemical kinetics. It provides a flexible and powerful tool for incorporating complex chemical kinetics into simulations of fluid dynamics. The package consists of two major software components: an Interpreter and a Gas-Phase Subroutine Library. The Interpreter is a program that reads a symbolic description of an elementary, user-specified chemical reaction mechanism. One output from the Interpreter is a data file that forms a link to the Gas-Phase Subroutine Library. This library is a collection of about 100 highly modular FORTRAN subroutines that may be called to return information on equations of state, thermodynamic properties, and chemical production rates. CHEMKIN-III includes capabilities for treating multi-fluid plasma systems, that are not in thermal equilibrium. These new capabilities allow researchers to describe chemistry systems that are characterized by more than one temperature, in which reactions may depend on temperatures associated with different species; i.e. reactions may be driven by collisions with electrons, ions, or charge-neutral species. These new features have been implemented in such a way as to require little or no changes to CHEMKIN implementation for systems in thermal equilibrium, where all species share the same gas temperature. CHEMKIN-III now has the capability to handle weakly ionized plasma chemistry, especially for application related to advanced semiconductor processing.
[Study on Chemical Kinetic Effect of Dielectric Barrier Discharge Plasma].
Zrang, Peng; Hong, Yan-ji; Shen, Shuang-yan; Ding, Xiao-yu; Ma, Di
2015-03-01
To reveal the mechanism of plasma (assisted the ignition process of methane/air further, schematic of dielectric barrier discharge plasma system with atmospheric air was designed and set up, the emission spectrum of dielectric barrier discharge plasma with atmospheric air was measured, and the active particles produced by the interaction of dielectric barrier discharge plasma with atmospheric air were analyzed with the spectrum technology, the ignition model and calculation methods of sensitivity analysis and reaction path analysis were given, effects of NO and O3 on the ignition delay time were simulated, and the chemical kinetics mechanism of NO and O3 assisted ignition was revealed via sensitivity analysis and reaction path analysis. The results show that main excited particles of N2 and O3 are generated via effect of plasma on the atmospheric air, which are converted into active particles of NO(ξ) and O3 in the end, the life of which are longer than any other active particles, effects of plasma on the ignition is simplified as effects of NO(ξ) and O3 on the ignition; NO and O3 could reduce the ignition delay time significantly, but the amplitude decrease with increase of the initial temperature, this is because the rate of ignition is decided by the oxidation rate of CH3, the oxidized pathway of CH3 is R155 and R156 for auto-ignition and their rates are slower when temperature is low, so the ignition delay time of methane/air is longer; NO could reduce the ignition delay time significantly because of the oxidized pathway of CH3 is changed to R327 CH3O2 + NO = CH3O + NO2, R328 CH3 + NO2 = CH3O + NO for NO(ξ) (assisted ignition process from R155 and R156 for auto-ignition; and the chemical kinetic effect is the dominating factor of O3 on the ignition and which change the reaction path. PMID:26117883
Kinetic studies of microinstabilities in toroidal plasmas: Simulation and theory
Lee, W.W.; Haham, T.S.; Parker, S.E.; Perkins, F.W.; Rath, S.; Rewoldt, G.; Reynders, J.V.W.; Santoro, R.A.; Tang, W.M.
1992-12-01
A comprehensive program for the development and use of particle simulation techniques for solving the gyrokinetic Vlasov-Maxwell equations on massively parallel computers has been carried out at Princeton Plasma Physics Laboratory. This is a key element of our ongoing theoretical efforts to systematically investigate physics issues vital to understanding tokamak plasmas. In this paper, our focus is on spatial-gradient-driven microinstabilities. Their importance is supported by the recent progress in achieving a physics-based understanding of anomalous transport in toroidal systems which has been based on the proposition that these drift-type electrostatic modes dependent on ion temperature gradient (ITG) and trapped particle effects are dominant in the bulk ( confinement'') region. Although their presence is consistent with a number of significant confinement trends, results from high temperature tokamaks such as TFTR have highlighted the need for better insight into the nonlinear properties of such instabilities in long-mean-free-path plasmas. In addressing this general issue, we report important new results including (i) the first fully toroidal 3D gyrokinetic simulation of ITG modes and (ii) realistic toroidal eigenmode calculations demonstrating the unique capability to deal with large scale kinetic behavior extending over many rational surfaces. The effects of ITG modes (iii) on the inward pinch of impurities in 3D slab geometry and (iv) on the existence of microtearing modes in 2D slab are also discussed. Finally, (v) sheared toroidal flow effects on trapped-particle modes are presented.
Kinetic studies of microinstabilities in toroidal plasmas: Simulation and theory
Lee, W.W.; Haham, T.S.; Parker, S.E.; Perkins, F.W.; Rath, S.; Rewoldt, G.; Reynders, J.V.W.; Santoro, R.A.; Tang, W.M.
1992-12-01
A comprehensive program for the development and use of particle simulation techniques for solving the gyrokinetic Vlasov-Maxwell equations on massively parallel computers has been carried out at Princeton Plasma Physics Laboratory. This is a key element of our ongoing theoretical efforts to systematically investigate physics issues vital to understanding tokamak plasmas. In this paper, our focus is on spatial-gradient-driven microinstabilities. Their importance is supported by the recent progress in achieving a physics-based understanding of anomalous transport in toroidal systems which has been based on the proposition that these drift-type electrostatic modes dependent on ion temperature gradient (ITG) and trapped particle effects are dominant in the bulk (``confinement``) region. Although their presence is consistent with a number of significant confinement trends, results from high temperature tokamaks such as TFTR have highlighted the need for better insight into the nonlinear properties of such instabilities in long-mean-free-path plasmas. In addressing this general issue, we report important new results including (i) the first fully toroidal 3D gyrokinetic simulation of ITG modes and (ii) realistic toroidal eigenmode calculations demonstrating the unique capability to deal with large scale kinetic behavior extending over many rational surfaces. The effects of ITG modes (iii) on the inward pinch of impurities in 3D slab geometry and (iv) on the existence of microtearing modes in 2D slab are also discussed. Finally, (v) sheared toroidal flow effects on trapped-particle modes are presented.
Dusty Plasmas - Kinetic Studies of Strong Coupling Phenomena
NASA Astrophysics Data System (ADS)
Morfill, Gregor
2011-10-01
``Dusty plasmas'' can be found almost everywhere - in the interstellar medium, in star and planet formation, in the solar system in the Earth's atmosphere and in the laboratory. In astrophysical plasmas the dust component accounts for only about 1% of the mass, nevertheless this component has a profound influence on the thermodynamics, the chemistry and the dynamics. Important physical processes are charging, sputtering, cooling, light absorption and radiation pressure, connecting electromagnetic forces to gravity. Surface chemistry is another important aspect. In the laboratory there is great interest in industrial processes (e.g. etching, vapor deposition) and at the fundamental physics level - the main topic here - the study of strong coupling phenomena. Here the dust (or microparticles) are the dominant component of the multi-species plasma. The particles can be observed in real time and pace, individually resolved at all relevant length and time scales. This provides an unprecedented means for studying self-organisation processes in many particle systems including the onset of cooperative phenomena. Due to the comparatively large mass of the microparticles (10-12 to 10-9 g) precision experiments are performed on the ISS. The following topics will be discussed: Phase transitions, phase separation, electrorheology, flow phenomena including the onset of turbulence at the kinetic level.
Complex (dusty) plasmas-kinetic studies of strong coupling phenomena
Morfill, Gregor E.; Ivlev, Alexei V.; Thomas, Hubertus M.
2012-05-15
'Dusty plasmas' can be found almost everywhere-in the interstellar medium, in star and planet formation, in the solar system in the Earth's atmosphere, and in the laboratory. In astrophysical plasmas, the dust component accounts for only about 1% of the mass, nevertheless this component has a profound influence on the thermodynamics, the chemistry, and the dynamics. Important physical processes are charging, sputtering, cooling, light absorption, and radiation pressure, connecting electromagnetic forces to gravity. Surface chemistry is another important aspect. In the laboratory, there is great interest in industrial processes (e.g., etching, vapor deposition) and-at the fundamental level-in the physics of strong coupling phenomena. Here, the dust (or microparticles) are the dominant component of the multi-species plasma. The particles can be observed in real time and space, individually resolved at all relevant length and time scales. This provides an unprecedented means for studying self-organisation processes in many-particle systems, including the onset of cooperative phenomena. Due to the comparatively large mass of the microparticles (10{sup -12}to10{sup -9}g), precision experiments are performed on the ISS. The following topics will be discussed: Phase transitions, phase separation, electrorheology, flow phenomena including the onset of turbulence at the kinetic level.
Kinetic and dynamic probability-density-function descriptions of disperse turbulent two-phase flows
NASA Astrophysics Data System (ADS)
Minier, Jean-Pierre; Profeta, Christophe
2015-11-01
This article analyzes the status of two classical one-particle probability density function (PDF) descriptions of the dynamics of discrete particles dispersed in turbulent flows. The first PDF formulation considers only the process made up by particle position and velocity Zp=(xp,Up) and is represented by its PDF p (t ;yp,Vp) which is the solution of a kinetic PDF equation obtained through a flux closure based on the Furutsu-Novikov theorem. The second PDF formulation includes fluid variables into the particle state vector, for example, the fluid velocity seen by particles Zp=(xp,Up,Us) , and, consequently, handles an extended PDF p (t ;yp,Vp,Vs) which is the solution of a dynamic PDF equation. For high-Reynolds-number fluid flows, a typical formulation of the latter category relies on a Langevin model for the trajectories of the fluid seen or, conversely, on a Fokker-Planck equation for the extended PDF. In the present work, a new derivation of the kinetic PDF equation is worked out and new physical expressions of the dispersion tensors entering the kinetic PDF equation are obtained by starting from the extended PDF and integrating over the fluid seen. This demonstrates that, under the same assumption of a Gaussian colored noise and irrespective of the specific stochastic model chosen for the fluid seen, the kinetic PDF description is the marginal of a dynamic PDF one. However, a detailed analysis reveals that kinetic PDF models of particle dynamics in turbulent flows described by statistical correlations constitute incomplete stand-alone PDF descriptions and, moreover, that present kinetic-PDF equations are mathematically ill posed. This is shown to be the consequence of the non-Markovian characteristic of the stochastic process retained to describe the system and the use of an external colored noise. Furthermore, developments bring out that well-posed PDF descriptions are essentially due to a proper choice of the variables selected to describe physical systems
Kinetic and dynamic probability-density-function descriptions of disperse turbulent two-phase flows.
Minier, Jean-Pierre; Profeta, Christophe
2015-11-01
This article analyzes the status of two classical one-particle probability density function (PDF) descriptions of the dynamics of discrete particles dispersed in turbulent flows. The first PDF formulation considers only the process made up by particle position and velocity Z(p)=(x(p),U(p)) and is represented by its PDF p(t; y(p),V(p)) which is the solution of a kinetic PDF equation obtained through a flux closure based on the Furutsu-Novikov theorem. The second PDF formulation includes fluid variables into the particle state vector, for example, the fluid velocity seen by particles Z(p)=(x(p),U(p),U(s)), and, consequently, handles an extended PDF p(t; y(p),V(p),V(s)) which is the solution of a dynamic PDF equation. For high-Reynolds-number fluid flows, a typical formulation of the latter category relies on a Langevin model for the trajectories of the fluid seen or, conversely, on a Fokker-Planck equation for the extended PDF. In the present work, a new derivation of the kinetic PDF equation is worked out and new physical expressions of the dispersion tensors entering the kinetic PDF equation are obtained by starting from the extended PDF and integrating over the fluid seen. This demonstrates that, under the same assumption of a Gaussian colored noise and irrespective of the specific stochastic model chosen for the fluid seen, the kinetic PDF description is the marginal of a dynamic PDF one. However, a detailed analysis reveals that kinetic PDF models of particle dynamics in turbulent flows described by statistical correlations constitute incomplete stand-alone PDF descriptions and, moreover, that present kinetic-PDF equations are mathematically ill posed. This is shown to be the consequence of the non-Markovian characteristic of the stochastic process retained to describe the system and the use of an external colored noise. Furthermore, developments bring out that well-posed PDF descriptions are essentially due to a proper choice of the variables selected to
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
NASA Technical Reports Server (NTRS)
Lipatov, Alexander S.; Combi, Michael R.
2004-01-01
The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.
Advances in petascale kinetic plasma simulation with VPIC and Roadrunner
Bowers, Kevin J; Albright, Brian J; Yin, Lin; Daughton, William S; Roytershteyn, Vadim; Kwan, Thomas J T
2009-01-01
VPIC, a first-principles 3d electromagnetic charge-conserving relativistic kinetic particle-in-cell (PIC) code, was recently adapted to run on Los Alamos's Roadrunner, the first supercomputer to break a petaflop (10{sup 15} floating point operations per second) in the TOP500 supercomputer performance rankings. They give a brief overview of the modeling capabilities and optimization techniques used in VPIC and the computational characteristics of petascale supercomputers like Roadrunner. They then discuss three applications enabled by VPIC's unprecedented performance on Roadrunner: modeling laser plasma interaction in upcoming inertial confinement fusion experiments at the National Ignition Facility (NIF), modeling short pulse laser GeV ion acceleration and modeling reconnection in magnetic confinement fusion experiments.
Kinetic theory of weak turbulence in magnetized plasmas: Perpendicular propagation
Yoon, Peter H.
2015-08-15
The present paper formulates a weak turbulence theory in which electromagnetic perturbations are assumed to propagate in directions perpendicular to the ambient magnetic field. By assuming that all wave vectors lie in one direction transverse to the ambient magnetic field, the linear solution and second-order nonlinear solutions to the equation for the perturbed distribution function are obtained. Nonlinear perturbed current from the second-order nonlinearity is derived in general form, but the limiting situation of cold plasma temperature is taken in order to derive an explicit nonlinear wave kinetic equation that describes three-wave decay/coalescence interactions among X and Z modes. A potential application of the present formalism is also discussed.
Limitation of the ECRIS performance by kinetic plasma instabilities (invited).
Tarvainen, O; Kalvas, T; Koivisto, H; Komppula, J; Kronholm, R; Laulainen, J; Izotov, I; Mansfeld, D; Skalyga, V; Toivanen, V; Machicoane, G
2016-02-01
Electron cyclotron resonance ion source (ECRIS) plasmas are prone to kinetic instabilities due to anisotropic electron velocity distribution. The instabilities are associated with strong microwave emission and periodic bursts of energetic electrons escaping the magnetic confinement. The instabilities explain the periodic ms-scale oscillation of the extracted beam current observed with several high performance ECRISs and restrict the parameter space available for the optimization of extracted beam currents of highly charged ions. Experiments with the JYFL 14 GHz ECRIS have demonstrated that due to the instabilities the optimum Bmin-field is less than 0.8BECR, which is the value suggested by the semiempirical scaling laws guiding the design of ECRISs. PMID:26931921
Wakes in complex plasmas: A self-consistent kinetic theory
NASA Astrophysics Data System (ADS)
Kompaneets, Roman; Morfill, Gregor E.; Ivlev, Alexei V.
2016-06-01
In ground-based experiments with complex (dusty) plasmas, charged microparticles are levitated against gravity by an electric field, which also drives ion flow in the parent gas. Existing analytical approaches to describe the electrostatic interaction between microparticles in such conditions generally ignore the field and ion-neutral collisions, assuming free ion flow with a certain approximation for the ion velocity distribution function (usually a shifted Maxwellian). We provide a comprehensive analysis of our previously proposed self-consistent kinetic theory including the field, ion-neutral collisions, and the corresponding ion velocity distribution. We focus on various limiting cases and demonstrate how the interplay of these factors results in different forms of the shielding potential.
Limitation of the ECRIS performance by kinetic plasma instabilities (invited)
NASA Astrophysics Data System (ADS)
Tarvainen, O.; Kalvas, T.; Koivisto, H.; Komppula, J.; Kronholm, R.; Laulainen, J.; Izotov, I.; Mansfeld, D.; Skalyga, V.; Toivanen, V.; Machicoane, G.
2016-02-01
Electron cyclotron resonance ion source (ECRIS) plasmas are prone to kinetic instabilities due to anisotropic electron velocity distribution. The instabilities are associated with strong microwave emission and periodic bursts of energetic electrons escaping the magnetic confinement. The instabilities explain the periodic ms-scale oscillation of the extracted beam current observed with several high performance ECRISs and restrict the parameter space available for the optimization of extracted beam currents of highly charged ions. Experiments with the JYFL 14 GHz ECRIS have demonstrated that due to the instabilities the optimum Bmin-field is less than 0.8BECR, which is the value suggested by the semiempirical scaling laws guiding the design of ECRISs.
Wakes in complex plasmas: A self-consistent kinetic theory.
Kompaneets, Roman; Morfill, Gregor E; Ivlev, Alexei V
2016-06-01
In ground-based experiments with complex (dusty) plasmas, charged microparticles are levitated against gravity by an electric field, which also drives ion flow in the parent gas. Existing analytical approaches to describe the electrostatic interaction between microparticles in such conditions generally ignore the field and ion-neutral collisions, assuming free ion flow with a certain approximation for the ion velocity distribution function (usually a shifted Maxwellian). We provide a comprehensive analysis of our previously proposed self-consistent kinetic theory including the field, ion-neutral collisions, and the corresponding ion velocity distribution. We focus on various limiting cases and demonstrate how the interplay of these factors results in different forms of the shielding potential. PMID:27415371
Cremaschini, Claudio Stuchlík, Zdeněk; Tessarotto, Massimo
2014-05-15
Astrophysical plasmas in the surrounding of compact objects and subject to intense gravitational and electromagnetic fields are believed to give rise to relativistic regimes. Theoretical and observational evidences suggest that magnetized plasmas of this type are collisionless and can persist for long times (e.g., with respect to a distant observer, coordinate, time), while exhibiting geometrical structures characterized by the absence of well-defined spatial symmetries. In this paper, the problem is posed whether such configurations can correspond to some kind of kinetic equilibrium. The issue is addressed from a theoretical perspective in the framework of a covariant Vlasov statistical description, which relies on the method of invariants. For this purpose, a systematic covariant variational formulation of gyrokinetic theory is developed, which holds without requiring any symmetry condition on the background fields. As a result, an asymptotic representation of the relativistic particle magnetic moment is obtained from its formal exact solution, in terms of a suitably defined invariant series expansion parameter (perturbative representation). On such a basis, it is shown that spatially non-symmetric kinetic equilibria can actually be determined, an example being provided by Gaussian-like distributions. As an application, the physical mechanisms related to the occurrence of a non-vanishing equilibrium fluid 4-flow are investigated.
NASA Astrophysics Data System (ADS)
Chang, Won-Seok; Yu, Dong-Hun; Cho, Deog-Gyun; Yook, Yeong-Geun; Chun, Poo-Reum; Lee, Se-Ah; Kwon, Deuk-Chul; Im, Yeon-Ho
2013-09-01
With continuous decrease of nanoscale design rule, plasma etching processes to form high aspect ratio contact hole still remains a challenge to overcome their inherent drawbacks such as bowing and twisted feature. Due to their complexities there still exist big gaps between current research status and predictable modeling of this process. To address this issue, we proposed a surface kinetic model of silicon nitride etch process under inductively coupled fluorocarbon plasmas. For this work, the cut-off probe and quadrapole mass spectroscopy were used for measuring electrical plasma properties, the ion and neutral radical species. Furthermore, the systematic surface analysis was performed to investigate the thickness and chemical bonding of polymer passivation layer during the etch process. The proposed semi-global surface kinetic model can consider deposition of polymer passivation layer and silicon nitride etching self-consistently. The predicted modeling results showed good agreement with experimental data. We believe that our research will provide valuable information to avoid the empirical development of plasma etching process.
NASA Astrophysics Data System (ADS)
Dodd, E. S.; Barnes, D. C.; Bezzerides, B.; Dubois, D. F.; Vu, H. X.
2003-10-01
RPIC is a reduced-description PIC code designed to investigate laser-plasma instabilities (LPI) in physical systems with vastly-different time scales prevalent under ICF conditions(H.X. Vu, B. Bezzerides, D.F. DuBois, J. Comp. Phys. 156), 12 (1999)., typically studied with the extended Zakharov model. Comparisons between the extended Zakharov model and RPIC were presented in a series of papers(K.Y. Sanbonmatsu, H.X. Vu, D.F. DuBois, and B. Bezzerides, Phys. Rev. Lett. 82), 932 (1999); K.Y. Sanbonmatsu, H.X. Vu, B. Bezzerides, and D.F. DuBois, Phys. Plasmas. 7, 1723,2824 (2000)., where quantitative agreements are obtained in the fluid and quasi-linear regime. In the kinetic regime where particle trapping is important, differences were found. The RPIC model itself is limited, e.g., Langmuir wave frequency harmonics are neglected. Our goal is two fold in comparing RPIC with full PIC in 1-d. First, advantages of RPIC over full PIC will be quantitatively assessed. Second, for strong laser drives, harmonics may be important to LPI physics. We would like to establish the regime of validity for RPIC, and to assess if the regimes where RPIC fails is of interest to ICF indirect drive.
NASA Astrophysics Data System (ADS)
Simakov, Andrei N.; Molvig, Kim
2016-03-01
Paper I [A. N. Simakov and K. Molvig, Phys. Plasmas 23, 032115 (2016)] obtained a fluid description for an unmagnetized collisional plasma with multiple ion species. To evaluate collisional plasma transport fluxes, required for such a description, two linear systems of equations need to be solved to obtain corresponding transport coefficients. In general, this should be done numerically. Herein, the general formalism is used to obtain analytical expressions for such fluxes for several specific cases of interest: a deuterium-tritium plasma; a plasma containing two ion species with strongly disparate masses, which agrees with previously obtained results; and a three ion species plasma made of deuterium, tritium, and gold. These results can be used for understanding the behavior of the aforementioned plasmas, or for verifying a code implementation of the general multi-ion formalism.
Diamagnetic boundary layers - A kinetic theory. [for collisionless magnetized plasmas
NASA Technical Reports Server (NTRS)
Lemaire, J.; Burlaga, L. F.
1976-01-01
A kinetic theory is presented for boundary layers associated with MHD tangential 'discontinuities' in a collisionless magnetized plasma, such as those observed in the solar wind. The theory consists of finding self-consistent solutions of Vlasov's equation and Maxwell's equation for stationary one-dimensional boundary layers separating two Maxwellian plasma states. Layers in which the current is carried by electrons are found to have a thickness of the order of a few electron gyroradii, but the drift speed of the current-carrying electrons is found to exceed the Alfven speed, and accordingly such layers are not stable. Several types of layers in which the current is carried by protons are discussed; in particular, cases are considered in which the magnetic-field intensity, direction, or both, changed across the layer. In every case, the thickness was of the order of a few proton gyroradii, and the field changed smoothly, although the characteristics depended somewhat on the boundary conditions. The drift speed was always less than the Alfven speed, consistent with stability of such structures. These results are consistent with observations of boundary layers in the solar wind near 1 AU.
Effect of grapefruit juice on plasma chloroquine kinetics in mice.
Ali, B H; Al-Qarawi, A; Mousa, H M
2002-08-01
1. It is known that grapefruit juice (GFJ) may interact with drugs concomitantly administered by inhibiting first-pass metabolism during the intestinal absorption phase. However, its interaction with chloroquine has not been studied previously. 2. Grapefruit juice (4 mL/kg) was given orally to mice 1 h prior to oral administration of chloroquine (100 mg/kg) and the concentration of the latter drug was measured fluorometrically in the plasma 0, 0.5, 0.75, 1, 2, 3, 4, 6, 8, 12, 18 and 24 h after its administration. 3. The mean (+/-SEM) of area under the curve values after administration of water +/- control) and GFJ were 5.34 +/- 0.38 and 7.01 +/- 0.66 mg.h/L, respectively. The corresponding mean C(max) values were 763.4 +/- 39.1 and 859.2 +/- 45.2 mg/L and the corresponding T(max) values (median) were 2.65 and 2.95 h. 4. The results suggest that GFJ coingestion increased the plasma concentration of chloroquine and altered some kinetic parameters of chloroquine. The clinical significance of this interaction in patients with malaria needs to be investigated. PMID:12100003
Fundamentals of the Plasma Sail Concept: MHD and Kinetic Studies
NASA Technical Reports Server (NTRS)
Khazanov, G.; Delamere, P.; Kabin, K.; Linde, T. J.; Krivorutsky, E.
2003-01-01
The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000] predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Such velocities will enable travel out of the solar system in period of about seven years-almost an order of magnitude improvement over present chemical based propulsion systems. However, for the parameters of the simulation of Winglee et al. [2000], a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the MHD fluid model, normally applied to planetary magnetospheres, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or, comparable to, the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. We have adopted a two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the momentum transfer to the innermost regions of the plasma sail is negligible.
Diomede, P.; Michau, A.; Redolfi, M.; Hassouni, K.; Morscheidt, W.; Longo, S.; Capitelli, M.
2008-10-15
A comparative study of two models, a kinetic and a fluid one, for the numerical simulation of H{sub 2} plasmas in parallel plate capacitively coupled rf discharges in one space dimension is presented. Both models adopt a multicomponent plasma description including several ionic species, take into account the self-consistent electric field by solving the Poisson equation and include a self-consistent coupling with the vibrational kinetics of the electronic ground state of H{sub 2} molecules. The peculiarities of this particular test case for model comparison with respect to previous studies are highlighted. The merits and the limitations of both approaches are discussed and the results are compared, in particular the steady-state density of charged species, the space-time variation of the electron energy, the vibrational distribution, and the atomic density.
NASA Astrophysics Data System (ADS)
Bénisti, Didier; Morice, Olivier; Gremillet, Laurent; Siminos, Evangelos; Strozzi, David J.
2010-10-01
In this paper, we present our nonlinear kinetic modeling of stimulated Raman scattering in a uniform and collisionless plasma using envelope equations. We recall the derivation of these equations, as well as our theoretical predictions for each of the nonlinear kinetic terms, the precision of which having been carefully checked against Vlasov simulations. We particularly focus here on the numerical resolution of these equations, which requires the additional concept of "self-optimization" that we explain, and we describe the envelope code BRAMA that we used. As an application of our modeling, we present one-dimensional BRAMA simulations of stimulated Raman scattering which predict threshold intensities, as well as time scales for Raman growth above threshold, in very good agreement with those inferred from Vlasov simulations. Finally, we discuss the differences between our modeling and other published ones.
Electrical and kinetic model of an atmospheric rf device for plasma aerodynamics applications
Pinheiro, Mario J.; Martins, Alexandre A.
2010-08-15
The asymmetrically mounted flat plasma actuator is investigated using a self-consistent two-dimensional fluid model at atmospheric pressure. The computational model assumes the drift-diffusion approximation and uses a simple plasma kinetic model. It investigated the electrical and kinetic properties of the plasma, calculated the charged species concentrations, surface charge density, electrohydrodynamic forces, and gas speed. The present computational model contributes to understand the main physical mechanisms, and suggests ways to improve its performance.
The Plasma Interaction Experiment (PIX) description and test program. [electrometers
NASA Technical Reports Server (NTRS)
Ignaczak, L. R.; Haley, F. A.; Domino, E. J.; Culp, D. H.; Shaker, F. J.
1978-01-01
The plasma interaction experiment (PIX) is a battery powered preprogrammed auxiliary payload on the LANDSAT-C launch. This experiment is part of a larger program to investigate space plasma interactions with spacecraft surfaces and components. The varying plasma densities encountered during available telemetry coverage periods are deemed sufficient to determine first order interactions between the space plasma environment and the biased experimental surfaces. The specific objectives of the PIX flight experiment are to measure the plasma coupling current and the negative voltage breakdown characteristics of a solar array segment and a gold plated steel disk. Measurements will be made over a range of surface voltages up to plus or minus kilovolt. The orbital environment will provide a range of plasma densities. The experimental surfaces will be voltage biased in a preprogrammed step sequence to optimize the data returned for each plasma region and for the available telemetry coverage.
A mean-field thermodynamic description of the kinetics of overdriven interfaces
NASA Astrophysics Data System (ADS)
Haxhimali, Tomorr; Belof, Jonathan; Sadigh, Babak
A key aspect of an accurate description of shock-induced structural phase transitions is the rigorous computation of the dynamics of the interfaces between coexisting phases. In the wake of the shock, the system will be exposed to strong gradient fields that give rise to overdriven interfaces during the induced phase transformation. In this work we take a mean-field approach using a time-dependent Ginzburg-Landau formalism to describe the dynamics of such overdriven interfaces. We make a connection of the mean-field result to a quasi-Langevin description, the Kardar-Parisi-Zhang (KPZ) equation, of the kinetics of the interface. Further, larger coarse-grained descriptions of the phase transition such as the Kolmogorov-Johnson-Mehl-Avrami (KJMA) model, which are commonly coupled to hydrodynamic equations that describe the evolution of the temperature and pressure during the shock propagation, ignore the details of the dynamics and structure of the interfacial regions. Overlaying the KPZ description of the interface evolution to these coarse-grained methods will result in physically more accurate multiscale models for shock propagation. We will present results from our efforts in this regard. This work is performed under the auspices of the U. S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Germaschewski, K.; Raeder, J.; Ruhl, H.
2010-12-01
Advances in processor technology provide the opportunity to simulate space plasma dynamics at unprecedented resolution. As processor clock speeds have begun to plateau in recent years, new technologies have emerged that maintain exponential growth in computational capability, in particular multi-core processors and heterogeneous approaches to computing, e.g., the STI Cell processor and general purpose GPUs. We will discuss two projects that aim at porting existing codes to efficiently run on heterogeneous processors. The Particle Simulation Code (PSC) is a 3D fully electromagnetic particle-in-cell code, solving the kinetic plasma equations, including a collision operator. This code is applied to problems requiring a kinetic model, like particle acceleration and modeling the microscopic structure of a reconnecting current sheets. We will discuss the performance gains enabled by porting the code to NVIDIA's GPU CUDA programming environment, as well as the challenges in exploiting the full capabilities of GPUs for the current deposition step. OpenGGCM is a community global magnetosphere model. The main computational challenge is the solution of the 3D MHD equations which are discretized using finite-difference / finite-volume. We ported this code to the Cell processor using a novel code generator. This approach allows us to specify the discretized equations in near-symbolic form as a stencil computation, and then have highly-optimized code be generated automatically. From the same description we are able to generate plain C code, C code with SIMD/SSE2 extensions and code for the Cell processor, yielding significant performance gains. We will also present first results of a new extension to the code generator that creates CUDA code for GPUs.
Galvao, R. A.; Ziebell, L. F.
2012-09-15
In this work, we detail the derivation of a plasma kinetic theory leading to the components of the dielectric tensor for a magnetized dusty plasma with variable charge on the dust particles, considering that the dust component of the plasma contains spherical dust particles with different sizes, which are charged both by inelastic collisions of electrons and ions and by photoionization.
SOLAR WIND TURBULENT SPECTRUM AT PLASMA KINETIC SCALES
Alexandrova, O.; Lacombe, C.; Mangeney, A.; Maksimovic, M.; Grappin, R.
2012-12-01
The description of the turbulent spectrum of magnetic fluctuations in the solar wind in the kinetic range of scales is not yet completely established. Here, we perform a statistical study of 100 spectra measured by the STAFF instrument on the Cluster mission, which allows us to resolve turbulent fluctuations from ion scales down to a fraction of electron scales, i.e., from {approx}10{sup 2} km to {approx}300 m. We show that for k {rho} {sub e} in [0.03, 3] (which corresponds approximately to the frequency in the spacecraft frame f in [3, 300] Hz), all the observed spectra can be described by a general law E(k ){proportional_to}k {sup -8/3} exp (- k {rho} {sub e}), where k is the wavevector component normal to the background magnetic field and {rho} {sub e} the electron Larmor radius. This exponential tail found in the solar wind seems compatible with the Landau damping of magnetic fluctuations onto electrons.
NASA Astrophysics Data System (ADS)
Vu, H. X.; Bezzerides, B.; Dubois, D. F.
1998-11-01
A fully kinetic, reduced-description particle-in-cell (RPIC) model is presented in which deviations from quasineutrality, electron and ion kinetic effects, and nonlinear interactions between low-frequency and high-frequency parametric instabilities are modeled correctly. The model is based on a reduced description where the electromagnetic field is represented by three separate temporal WKB envelopes in order to model low-frequency and high-frequency parametric instabilities. Because temporal WKB approximations are invoked, the simulation can be performed on the electron time scale instead of the time scale of the light waves. The electrons and ions are represented by discrete finite-size particles, permitting electron and ion kinetic effects to be modeled properly. The Poisson equation is utilized to ensure that space-charge effects are included. Although RPIC is fully three dimensional, it has been implemented in only two dimensions on a CRAY-T3D with 512 processors and on the Accelerated Strategic Computing Initiative (ASCI) parallel computer at Los Alamos National Laboratory, and the resulting simulation code has been named ASPEN. Given the current computers available to the authors, one and two dimensional simulations are feasible to, and have been, performed. Three dimensional simulations are much more expensive, and are not feasible at this time. However, with rapidly advancing computer technologies, three dimensional simulations may be feasible in the near future. We believe this code is the first PIC code capable of simulating the interaction between low-frequency and high-frequency parametric instabilites in multiple dimensions. Test simulations of stimulated Raman scattering (SRS), stimulated Brillouin scattering (SBS), and Langmuir decay instability (LDI), are presented.
Numerical modeling of radiation physics in kinetic plasmas [II
NASA Astrophysics Data System (ADS)
Paraschiv, Ioana; Sentoku, Yasuhiko; Mancini, Roberto
2014-10-01
X-ray radiation is an important feature of ultra-intense laser interactions with high Z materials. In order to take into account the radiation effects in the high energy density plasmas created by such interactions, we have modified the collisional particle-in-cell code PICLS to self-consistently model the x-ray radiation transport (RT). Solving the equation of radiation transport requires the creation of a non-LTE database of emissivities and opacities as functions of photon frequency for given densities, bulk electron temperatures, hot electron temperatures, and hot electron fractions. The database was generated using results computed by a non-equilibrium, collisional-radiative atomic kinetics code. Using the two-dimensional RT-PICLS code we have studied the X-ray transport in an ultrafast heated target and the dependence of the emitted K- α radiation on the fast electron dynamics in the solid target. The details of these results obtained from the implementation of the radiation transport model into the PICLS calculations will be reported in this presentation. Work supported by the DOE Office of Science Grant No. DE-SC0008827 and by the NNSA/DOE Grants No. DE-FC52-06NA27616 and DE-NA0002075.
Richardson Extrapolation Based Error Estimation for Stochastic Kinetic Plasma Simulations
NASA Astrophysics Data System (ADS)
Cartwright, Keigh
2014-10-01
To have a high degree of confidence in simulations one needs code verification, validation, solution verification and uncertainty qualification. This talk will focus on numerical error estimation for stochastic kinetic plasma simulations using the Particle-In-Cell (PIC) method and how it impacts the code verification and validation. A technique Is developed to determine the full converged solution with error bounds from the stochastic output of a Particle-In-Cell code with multiple convergence parameters (e.g. ?t, ?x, and macro particle weight). The core of this method is a multi parameter regression based on a second-order error convergence model with arbitrary convergence rates. Stochastic uncertainties in the data set are propagated through the model usin gstandard bootstrapping on a redundant data sets, while a suite of nine regression models introduces uncertainties in the fitting process. These techniques are demonstrated on Flasov-Poisson Child-Langmuir diode, relaxation of an electro distribution to a Maxwellian due to collisions and undriven sheaths and pre-sheaths. Sandia National Laboratories is a multie-program laboratory managed and operated by Sandia Corporation, a wholly owned subisidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under Contract DE-AC04-94AL85000.
Simakov, Andrei Nikolaevich; Molvig, Kim
2016-03-17
Paper I [A. N. Simakov and K. Molvig, Phys. Plasmas23, 032115 (2016)] obtained a fluid description for an unmagnetized collisional plasma with multiple ion species. To evaluate collisional plasmatransport fluxes, required for such a description, two linear systems of equations need to be solved to obtain corresponding transport coefficients. In general, this should be done numerically. Herein, the general formalism is used to obtain analytical expressions for such fluxes for several specific cases of interest: a deuterium-tritium plasma; a plasma containing two ion species with strongly disparate masses, which agrees with previously obtained results; and a three ion species plasmamore » made of deuterium, tritium, and gold. We find that these results can be used for understanding the behavior of the aforementioned plasmas, or for verifying a code implementation of the general multi-ion formalism.« less
Implicit Methods for the Magnetohydrodynamic Description of Magnetically Confined Plasmas
Jardin, S C
2010-09-28
Implicit algorithms are essential for predicting the slow growth and saturation of global instabilities in today’s magnetically confined fusion plasma experiments. Present day algorithms for obtaining implicit solutions to the magnetohydrodynamic (MHD) equations for highly magnetized plasma have their roots in algorithms used in the 1960s and 1970s. However, today’s computers and modern linear and non-linear solver techniques make practical much more comprehensive implicit algorithms than were previously possible. Combining these advanced implicit algorithms with highly accurate spatial representations of the vector fields describing the plasma flow and magnetic fields and with improved methods of calculating anisotropic thermal conduction now makes possible simulations of fusion experiments using realistic values of plasma parameters and actual configuration geometry.
Consistent Holographic Description of Boost-Invariant Plasma
Heller, Michal P.; Surowka, Piotr; Loganayagam, R.; Spalinski, Michal; Vazquez, Samuel E.
2009-01-30
Prior attempts to construct the gravity dual of boost-invariant flow of N=4 supersymmetric Yang-Mills gauge theory plasma suffered from apparent curvature singularities in the late-time expansion. This Letter shows how these problems can be resolved by a different choice of expansion parameter. The calculations presented correctly reproduce the plasma energy-momentum tensor within the framework of second-order viscous hydrodynamics.
Effects of kinetic processes in shaping Io's global plasma environment: A 3D hybrid model
NASA Astrophysics Data System (ADS)
Lipatov, Alexander S.; Combi, Michael R.
2006-02-01
The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the I0 flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io's environment are discussed in this paper. Comparison with the Galileo I0 flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
NASA Technical Reports Server (NTRS)
Lipatov, Alexander S.; Combi, Michael R.
2006-01-01
The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.
Kinetic approach for the ion drag force in a collisional plasma
Ivlev, A.V.; Zhdanov, S.K.; Khrapak, S.A.; Morfill, G.E.
2005-01-01
The linear kinetic approach to calculate the ion drag force in a collisional plasma is generalized. The model collision integral (for ion-neutral collisions) is discussed and employed to calculate the plasma response for arbitrary velocity of the plasma flow and arbitrary frequency of the collisions. The derived plasma response is used to calculate the self-consistent force on the test charged particle. The obtained results are compared to those of the traditional pair collision approach, and the importance of the self-consistent kinetic consideration is highlighted. In conclusion, the applicability of the proposed approach is discussed.
Effects of dust particles in plasma kinetics: Ion dynamics time scales
Angelis, U. de; Tolias, P.; Ratynskaia, S.
2012-07-15
The self-consistent kinetic theory of dusty plasmas [V. N. Tsytovich and U. de Angelis, Phys. Plasmas 6, 1093 (1999)] is extended to frequency regimes relevant for ion dynamics, accounting for both constant and fluctuating plasma sources. In contrast to earlier models, binary plasma collisions are no longer neglected with respect to collisions with dust; hence, the model developed here is also valid for low dust densities. Expressions are found for the system's permittivity, the ion collision integral, and the spectral densities of ion density fluctuations. The structure of the ion kinetic equation is analyzed, and applications of the model for both astrophysical and laboratory environments are discussed.
Effects of dust particles in plasma kinetics: Ion dynamics time scales
NASA Astrophysics Data System (ADS)
de Angelis, U.; Tolias, P.; Ratynskaia, S.
2012-07-01
The self-consistent kinetic theory of dusty plasmas [V. N. Tsytovich and U. de Angelis, Phys. Plasmas 6, 1093 (1999)] is extended to frequency regimes relevant for ion dynamics, accounting for both constant and fluctuating plasma sources. In contrast to earlier models, binary plasma collisions are no longer neglected with respect to collisions with dust; hence, the model developed here is also valid for low dust densities. Expressions are found for the system's permittivity, the ion collision integral, and the spectral densities of ion density fluctuations. The structure of the ion kinetic equation is analyzed, and applications of the model for both astrophysical and laboratory environments are discussed.
Effect of antenna size on electron kinetics in inductively coupled plasmas
Lee, Hyo-Chang; Chung, Chin-Wook
2013-10-15
Spatially resolved measurements of electron energy distribution functions (EEDFs) are investigated in inductively coupled plasmas with two planar antenna coils. When the plasma is sustained by the antenna with a diameter of 18 cm, the nonlocal kinetics is preserved in the argon gas pressure range from 2 mTorr to 20 mTorr. However, electron kinetics transit from nonlocal kinetics to local kinetics in discharge sustained by the antenna coil with diameter 34 cm. The results suggest that antenna size as well as chamber length are important parameters for the transition of the electron kinetics. Spatial variations of plasma potential, effective electron temperature, and EEDF in terms of total electron energy scale are also presented.
THEMIS observation of Kinetic Ballooning/Interchange Waves in the High Bz Plasma Sheet
NASA Astrophysics Data System (ADS)
Panov, Evgeny V.; Nakamura, Rumi; Kubyshkina, Marina V.; Baumjohann, Wolfgang; A, Sergeev, Victor
2015-04-01
Using THEMIS observations of plasma sheet oscillations with kinetic ballooning/interchange instability (BICI) signatures, we investigate the properties of the waves when a high background plasma sheet Bz is seen. We find that such waves are in a better agreement with the existing kinetic simulations. Using adapted Tsyganenko models, we also show conjugate all-sky camera observations in the course of the development of the waves.
Numerical description of discharge characteristics of the plasma needle
Brok, W.J.M.; Bowden, M.D.; Dijk, J. van; Mullen, J.J.A.M. van der; Kroesen, G.M.W.
2005-07-01
The plasma needle is a small atmospheric, nonthermal, radio-frequency discharge, generated at the tip of a needle, which can be used for localized disinfection of biological tissues. Although several experiments have characterized various qualities of the plasma needle, discharge characteristics and electrical properties are still not well known. In order to provide initial estimates on electrical properties and quantities such as particle densities, we employed a two-dimensional, time-dependent fluid model to describe the plasma needle. In this model the balance equation is solved in the drift-diffusion approach for various species and the electron energy, as well as Poisson's equation. We found that the plasma production occurs in the sheath region and results in a steady flux of reactive species outwards. Even at small (<0.1%) admixtures of N{sub 2} to the He background, N{sub 2}{sup +} is the dominant ion. The electron density is typically 10{sup 11} cm{sup -3} and the dissipated power is in the order of 10 mW. These results are consistent with the experimental data available and can give direction to the practical development of the plasma needle.
Note on quantitatively correct simulations of the kinetic beam-plasma instability
Lotov, K. V.; Timofeev, I. V.; Mesyats, E. A.; Snytnikov, A. V.; Vshivkov, V. A.
2015-02-15
A large number of model particles are shown necessary for quantitatively correct simulations of the kinetic beam-plasma instability with the clouds-in-cells method. The required number of particles scales inversely with the expected growth rate, as only a narrow interval of beam velocities is resonant with the wave in the kinetic regime.
Molecular Dynamics Description of Partially Ionized Dense Plasmas
NASA Astrophysics Data System (ADS)
Lagattuta, Ken
2004-11-01
A report on work in progress: the approach to steady-state of partially ionized dense plasmas, containing more than one atomic element, is being simulated with the quasi-classical method known as Fermi Molecular Dynamics (FMD). We recap the FMD method, recalling its several advantages and disadvantages, and present an overview of past work. we have continued to develop the FMD method as a tool for simulating the behaviors of a variety of inhomogeneous, partially ionized, dense plasma systems, in cases for which more rigorous methods are still unavailable. Predictions of the average ionization state Z* of atoms, in a plasma containing more than one atomic element, is complicated by many factors, especially under conditions of high density, and not too high temperature. Average atom models become problematic when two or more atomic elements are present together. In order to address this problem, we have applied the FMD method to plasmas containing selected mixtures of atomic elements, determining Z* for each element over a range of temperatures and densities. LANL archived abstract: LA-UR-04-2186
Jupiter's magnetosphere: Plasma description from the Ulysses flyby
Bame, S.J.; Barraclough, B.L.; Feldman, W.C.; Gisler, G.R.; Gosling, J.T.; McComas, D.J.; Phillips, J.L.; Thomsen, M.F. ); Goldstein, B.E.; Neugebauer, M. )
1992-09-11
Plasma observations at Jupiter show that the outer regions of the Jovian magnetosphere are remarkably similar to those of Earth. Bow-shock precursor electrons and ions were detected in the upstream solar wind, as at Earth. Plasma changes across the bow shock and properties of the magnetosheath electrons were much like those at Earth, indicating that similar processes are operating. A boundary layer populated by a varying mixture of solar wind and magnetospheric plasmas was found inside the magnetopause, again as at Earth. In the middle magnetosphere, large electron density excursions were detected with a 10-hour periodicity as planetary rotation carried the tilted plasma sheet past Ulysses. Deep in the magnetosphere, Ulysses crossed a region, tentatively described as magnetically connected to the Jovian polar cap on one end and to the interplanetary magnetic field on the other. In the inner magnetosphere and Io torus, where corotation plays a dominant role, measurements could not be made because of extreme background rates from penetrating radiation belt particles.
The kinetic response of a stochastic plasma to low frequency perturbations
NASA Technical Reports Server (NTRS)
Hurricane, Omar A.; Pellat, Rene; Coroniti, Ferdinand V.
1994-01-01
Following suggestion that substorm breakup might be caused by an interchange or ballooning instability, several magnetohydrodynamic (MHD) and gyro-kinetic stability analysis have been performed for plasma sheet magnetic field geometries. However, the stochastic ion dynamics in the highly stressed, thin high-beta near earth plasma sheet violate the locality requirements of MHD and the invarience of the magnetic moment, mu, required by gyro-kinetic theories. In this paper we develop a new linear Vlasov kinetic theory (for low frequency modes omega less than omega(sub b), omega(sub b) being the bounce frequency) which includes the dynamics of stochastic ions.
NASA Astrophysics Data System (ADS)
Lie-Svendsen, Øystein; Olsen, Espen Lyngdal
1998-03-01
We have compared kinetic and fluid model descriptions of the proton polar wind outflow, from the collision-dominated, subsonic regime at lower altitudes, through the transition to supersonic flow, and well into the supersonic flow regime. The kinetic model is based on the Fokker-Planck collision operator, and the two fluid models employed are based on the 8-moment expansion and the 16-moment bi-Maxwellian expansions, respectively. We find excellent agreement between the kinetic description and the fluid models for the proton density and flux, even in the transonic and supersonic flow regimes. The models are also in qualitative agreement for the temperature and heat flux moments, although neither fluid model reproduces the negative (downward) kinetic heat fluxes found at high altitudes. The 16-moment fluid model gives a temperature anisotropy similar to the anisotropy derived from the kinetic solution. The assumed forms for the velocity distribution, on which the fluid expansions are based, do not agree with the kinetic velocity distribution, except in the subsonic region where the departure from a Maxwellian distribution is small. Near the fluid critical point the kinetic model develops a double-hump distribution, with an isotropic, low-energy core and an anisotropic, high-energy tail, and at higher altitudes the distribution function develops a ``kidney bean'' shape.
Association between plasma zinc concentration and zinc kinetic parameters in premenopausal women.
Yokoi, Katsuhiko; Egger, Norman G; Ramanujam, V M Sadagopa; Alcock, Nancy W; Dayal, Hari H; Penland, James G; Sandstead, Harold H
2003-11-01
The objective of this study was to measure relationships between plasma zinc (Zn) concentrations and Zn kinetic parameters and to measure relationships of Zn status with taste acuity, food frequency, and hair Zn in humans. The subjects were 33 premenopausal women not taking oral contraceptives and dietary supplements containing iron and Zn. Main outcomes were plasma Zn concentrations, Zn kinetic parameters based on the three-compartment mammillary model using 67Zn as a tracer, electrical taste detection thresholds, and food frequencies. Lower plasma Zn was significantly (P < 0.01) associated with smaller sizes of the central and the lesser peripheral Zn pools, faster disappearance of tracer from plasma, and higher transfer rate constants from the lesser peripheral pool to the central pool and from the central pool to the greater peripheral pool. The break points in the plasma Zn-Zn kinetics relationship were found between 9.94 and 11.5 micromol/l plasma Zn. Smaller size of the lesser peripheral pool was associated with lower frequency of beef consumption and higher frequency of bran breakfast cereal consumption. Hypozincemic women with plasma Zn <10.7 micromol/l or 700 ng/ml had decreased thresholds of electrical stimulation for gustatory nerves. Our results based on Zn kinetics support the conventional cutoff value of plasma Zn (10.7 micromol/l or 700 ng/ml) between normal and low Zn status. PMID:12865259
Turbulence in strongly coupled dusty plasmas using generalized hydrodynamic description
Tiwari, Sanat Kumar; Dharodi, Vikram Singh; Das, Amita; Patel, Bhavesh G.; Kaw, Predhiman
2015-02-15
The properties of decaying turbulence have been studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of strongly coupled dusty plasma medium in two dimensions. The GHD model treats the strongly coupled dusty plasma system as a visco-elastic medium. The incompressible limit of the GHD model is considered here. The studies carried out here are, however, applicable to a wider class of visco-elastic systems, and are not merely confined to the dusty plasma medium. Our simulations studies show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence, where Weissenberg's number also plays an important role on the onset of turbulence. This feature has been observed in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.
NASA Astrophysics Data System (ADS)
Dum, C. T.
1990-06-01
Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).
NASA Technical Reports Server (NTRS)
Dum, C. T.
1990-01-01
Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).
Role of plasma activation in the kinetics of CNT growth in PECVD process
NASA Astrophysics Data System (ADS)
Lebedeva, Irina; Gavrikov, Alexey; Baranov, Alexey; Belov, Maxim; Knizhnik, Andrey; Potapkin, Boris; Sommerer, Timothy
2009-10-01
The work presents kinetic modeling of the effect of acceleration for the growth kinetics of carbon nanotubes by hydrocarbon gas mixture modification with plasma discharge. The plasma activation creates active species in hydrocarbon gas mixture, which can easily adsorb and dissociate on the catalyst surface. So plasma treatment of the gas mixture in the CVD process allows to increase the carbon supply rate by a few orders of magnitude compared to that in thermal CVD process. On the other hand, plasma can also provide etching of carbon species from the catalyst surface. To correctly reproduce both of these effects of plasma, the kinetic model of growth of carbon nanotubes is developed based on first-principles analysis of heterogeneous processes on the catalyst surface and detailed kinetics of gas phase chemistry. The model is used to compare the growth rates of carbon nanotubes in thermal and plasma-enhanced CVD processes and to determine critical gas pressures, at which CNT growth kinetics switches from the adsorption limitation to the limitation by reaction and diffusion on the catalyst.
Kinetic effects on streaming instabilities in electron-positron-ion plasmas
NASA Astrophysics Data System (ADS)
Shan, S. Ali; Saleem, H.
2009-02-01
Streaming instabilities in electron-positron-ion plasmas are investigated using kinetic approach in several different limits. The effects of the variation of background temperatures of electrons Teo and positrons Tpo on the growth rates are also presented for the case of ion beam streaming into electron-positron plasmas and positrons beam streaming into electron-ion plasmas. It is noticed that the increase of number density of positrons gives a destabilizing trend to the electrostatic perturbations in the system.
BRIEF COMMUNICATION: On the drift kinetic equation driven by plasma flows
NASA Astrophysics Data System (ADS)
Shaing, K. C.
2010-07-01
A drift kinetic equation that is driven by plasma flows has previously been derived by Shaing and Spong 1990 (Phys. Fluids B 2 1190). The terms that are driven by particle speed that is parallel to the magnetic field B have been neglected. Here, such terms are discussed to examine their importance to the equation and to show that these terms do not contribute to the calculations of plasma viscosity in large aspect ratio toroidal plasmas, e.g. tokamaks and stellarators.
Plasma Boundaries and Kinetic-Scale Electric Field Structures in the Inner Magnetosphere
NASA Astrophysics Data System (ADS)
Malaspina, David; Larsen, Brian; Ergun, R. E.; Skoug, Ruth; Wygant, John; Reeves, Geoffrey; Jaynes, Allison
2016-07-01
Recent advances in spacecraft instrumentation have enabled fresh examination of coupling between macro-scale and micro-scale physics in the terrestrial magnetosphere, demonstrating not only that cross-scale interactions are a key component of magnetospheric dynamics, but also that plasma boundaries play a crucial role in mediating cross-scale coupling. We use Van Allen Probe observations to study the cross-scale interaction between inner magnetospheric plasma boundaries (including the plasmapause and injection fronts) and kinetic-scale electric field structures including kinetic Alfven waves, double layers, phase space holes, and nonlinear whistler mode waves. We focus on the spatial distribution of these kinetic structures in the inner magnetosphere and their interaction with plasma boundaries. We demonstrate that both the occurrence probability and amplitude of these structures peak at plasma boundaries. Further, it is found that regions of kinetic-scale electric field structure activity travel with plasma boundaries. These observations imply that kinetic-scale electric field structures are continually generated by instabilities localized to these boundaries, constraining their ability to energize radiation belt particles over large spatial regions.
Simulating Electron Cyclotron Resonance Heating in Kinetic and Dielectric Plasma Models with VORPAL
NASA Astrophysics Data System (ADS)
Roark, Christine; Smithe, David; Stoltz, Peter; Tech-X Corporation Team
2011-10-01
We present results of electron cyclotron resonance heating (ECRH) in a plasma sustained by microwaves using VORPAL. Specifically, we look at the electron temperature, sheath size, rate of plasma formation and power absorbed for simulations with an argon gas at 10s of mTorr pressure and 2.45 GHz. We look at the effects of including elastic, inelastic and ionizing Monte Carlo collisions on the formation of the kinetic plasma. We also discuss the use of higher-order particle algorithms for smoothing out the particle current and charge which can help reduce unphysical heating in PIC simulations of high pressure, low temperature plasmas and the effect this has on sheath size and electron temperature. We then compare these simulations to a method replacing the kinetic particles with an equivalent plasma dielectric model.
NASA Astrophysics Data System (ADS)
Qi, Lei; Kwon, Jaemin; Hahm, T. S.; Jo, Gahyung
2016-06-01
Nonlinear bounce-averaged kinetic theory [B. H. Fong and T. S. Hahm, Phys. Plasmas 6, 188 (1999)] is used for magnetically trapped electron dynamics for the purpose of achieving efficient gyrokinetic simulations of Trapped Electron Mode (TEM) and Ion Temperature Gradient mode with trapped electrons (ITG-TEM) in shaped tokamak plasmas. The bounce-averaged kinetic equations are explicitly extended to shaped plasma equilibria from the previous ones for concentric circular plasmas, and implemented to a global nonlinear gyrokinetic code, Gyro-Kinetic Plasma Simulation Program (gKPSP) [J. M. Kwon et al., Nucl. Fusion 52, 013004 (2012)]. Verification of gKPSP with the bounce-averaged kinetic trapped electrons in shaped plasmas is successfully carried out for linear properties of the ITG-TEM mode and Rosenbluth-Hinton residual zonal flow [M. N. Rosenbluth and F. L. Hinton, Phys. Rev. Lett. 80, 724 (1998)]. Physics responsible for stabilizing effects of elongation on both ITG mode and TEM is identified using global gKPSP simulations. These can be understood in terms of magnetic flux expansion, leading to the effective temperature gradient R / L T ( 1 - E ') [P. Angelino et al., Phys. Rev. Lett. 102, 195002 (2009)] and poloidal wave length contraction at low field side, resulting in the effective poloidal wave number kθρi/κ.
Analytical description of a neutral-induced tripole vortex in a plasma.
Vranjes, J; Okamoto, A; Yoshimura, S; Poedts, S; Kono, M; Tanaka, M Y
2002-12-23
An analytical description of a stationary triple vortex, observed in a cylindrical plasma, is presented. The concentration of neutrals, which is rather high in the experiment, turns out to be of crucial importance due to a spatially dependent distribution. In the radial direction the neutral concentration is paraboliclike, yielding an effective radial force directed towards the axis of the system. This neutral force causes the rotation of the plasma in the direction which is opposite to the E-->xB--> drift. The stationary triple vortex develops for a starting Gaussian-density distribution and a rigid-body rotation of the plasma column. PMID:12484828
Ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves
NASA Technical Reports Server (NTRS)
Moghaddam-Taaheri, E.; Goertz, C. K.; Smith, R. A.
1989-01-01
A two-dimensional quasi-linear numerical code was developed for studying ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves. The model assumes that the central plasma sheet is the particle source, and that the last magnetic field lines on which kinetic Alfven waves exist and diffusion occurs can be either open or closed. As the possible source for the excitement of the kinetic Alfven waves responsible for ion diffusion, the resonant mode conversion of the surface waves to kinetic Alfven waves is considered. It is shown that, depending on the topology of the magnetic field at the lobe side of the simulation system, i.e., on whether field lines are open or closed, the ion distribution function may or may not reach a steady state.
Sozarukova, M M; Polimova, A M; Proskurnina, E V; Vladimirov, Yu A
2016-01-01
Oxidative stress is a pathogenetic factor of many diseases. The control of its level is important for early diagnosis and therapy adjustment. In this work, antioxidant status was estimated in blood plasma. In the system of 2,2'-azo-bis(2-amidinopropane)dihydrochloride-luminol a set of chemiluminescence kinetic curve parameters is proposed for oxidative stress level estimation (the latent period τ(lat) and the increasing of analytical signal ΔI(CL)). Uric acid and albumin were shown as the main components that responsible for changes in chemiluminescence kinetic curve of plasma. Serum albumin undergoes oxidative modification in dose-depend manner under the action of UV irradiation, it causes the enhancement of antioxidant properties. Changes in plasma chemiluminescence kinetics are proposed as a measure of oxidative stress in human body. PMID:27192837
Modeling of High Kinetic Energy Plasma Jets for Fusion Applications
NASA Astrophysics Data System (ADS)
Bogatu, I. N.; Galkin, S. A.; Kim, J. S.
2006-10-01
We used semi-analytical models for high velocity (>200 km/s) and density (>10^17 cm-3) plasma jets to describe the acceleration in coaxial electrodes geometry, the collision, and plasma liner implosion, assuming that jets have merged into a spherical or cylindrical shell. The results are compared with experimental data and are being used for guiding LSP and MACH2 codes simulation and for optimization. The simplest model which uses the adiabatic invariant for oscillator revealed the basic relation between the velocity and the parameters of the plasma accelerator. Plasma slug model was extended for including friction and mass addition by electrode erosion. A simple model of blow-by instability by using the canting angle of the plasma current was formulated. As plasma jets collision at high interfacial Mach number generates shock fronts, we analyzed their possible consequences on the merging process and liner formation. The structure of the spherical shell liner during adiabatic implosion and the effect of the shock wave generated at void closure on the confinement time were also investigated.
NASA Astrophysics Data System (ADS)
Mulec, M.; Ivanov, I. B.; Heyn, M. F.; Kernbichler, W.
2012-03-01
Resistive wall modes (RWMs) are studied within the kinetic model proposed by Heyn et al. [Nucl. Fusion 46, S159 (2006); Phys. Plasmas 18, 022501 (2011)], which accounts for Landau damping, transit-time magnetic pumping, and Coulomb collisions in cylindrical geometry. Results for the reversed field pinch plasma are compared to the magnetohydrodynamic results obtained by Guo et al., [Phys. Plasmas 6, 3868 (1999)]. Stabilization of the external kink mode by an ideal wall as well as stabilization of the resistive wall mode by toroidal plasma rotation is obtained. In contrast to MHD modelling, which predicts a stability window for the resistive wall position, kinetic modelling predicts a one sided window only, i.e., the resistive wall must be sufficiently close to plasma to achieve rotational stabilization of the mode but there is no lower limit on the wall position. Stabilizing rotation speeds are found somewhat smaller when compared to MHD results. In addition, for the present plasma configuration, the kinetic model predicts resistive wall mode stabilization only in one direction of toroidal rotation. In the opposite direction, a destabilizing effect is observed. This is in contrast to MHD where mode stabilization is symmetric with respect to the direction of the toroidal plasma rotation.
Time-dependent gas phase kinetics in a hydrogen diluted silane plasma
NASA Astrophysics Data System (ADS)
Nunomura, S.; Yoshida, I.; Kondo, M.
2009-02-01
The gas phase kinetics in a high-pressure hydrogen diluted silane plasma has been studied at time scales of 10-2-6×102 s. The time-resolved gas phase composition shows the following kinetics at different time scales: silane decomposition and polysilane generation in ≲2×10-1 s, nanoparticle formation and plasma density reduction in 10-1-100 s, polysilane accumulation in 100-102 s, and silane depletion and electrode heating in ≳101 s. Disilane radicals are implied to be the dominant film precursors in addition to silyl radicals.
Time-dependent gas phase kinetics in a hydrogen diluted silane plasma
Nunomura, S.; Kondo, M.; Yoshida, I.
2009-02-16
The gas phase kinetics in a high-pressure hydrogen diluted silane plasma has been studied at time scales of 10{sup -2}-6x10{sup 2} s. The time-resolved gas phase composition shows the following kinetics at different time scales: silane decomposition and polysilane generation in < or approx. 2x10{sup -1} s, nanoparticle formation and plasma density reduction in 10{sup -1}-10{sup 0} s, polysilane accumulation in 10{sup 0}-10{sup 2} s, and silane depletion and electrode heating in > or approx. 10{sup 1} s. Disilane radicals are implied to be the dominant film precursors in addition to silyl radicals.
Description of plasma focus current sheath as the Turner relaxed state of a Hall magnetofluid
Auluck, S. K. H.
2009-12-15
The central mystery of plasma focus research is the two orders-of-magnitude-higher-than-thermal fusion reaction rate and the fact that both the space-resolved neutron spectra and space-resolved reaction proton spectra show features which can be ascribed only to a rotational motion of the center-of-mass of the reacting deuteron population. It has been suggested earlier [S. K. H. Auluck, IEEE Trans. Plasma Sci. 25, 37 (1997)] that this and other experimental observations can be consistently explained in terms of a hypothesis involving rotation of the current carrying plasma annulus behind the imploding gas-dynamic shock. Such rotation (more generally, mass flow) is an in-built feature of relaxed state of a two-fluid plasma [R. N. Sudan, Phys. Rev. Lett. 42, 1277 (1979)]. Relaxation in the 'Hall magnetofluid' approximation, in which the generalized Ohm's law includes the Hall effect term and the magnetic convection term but omits the contributions to the electric field from resistive dissipation, electron pressure gradient, thermoelectric effect, electron inertia, etc., has been extensively studied by many authors. In the present paper, Turner's [IEEE Trans. Plasma Sci. PS-14, 849 (1986)] degenerate solution for the relaxed state of the Hall magnetohydrodynamic plasma has been adapted to the case of an infinitely long annular current carrying plasma, a tractable idealization of the current sheath of a plasma focus. The resulting model is consistent with experimental values of ion kinetic energy and observation of predominantly radially directed neutron emission in good shots.
NASA Astrophysics Data System (ADS)
Wilczek, Sebastian; Trieschmann, Jan; Eremin, Denis; Brinkmann, Ralf Peter; Schulze, Julian; Schuengel, Edmund; Derzsi, Aranka; Korolov, Ihor; Hartmann, Peter; Donkó, Zoltán; Mussenbrock, Thomas
2016-06-01
Low pressure capacitive radio frequency (RF) plasmas are often described by equivalent circuit models based on fluid approaches that predict the self-excitation of resonances, e.g., high frequency oscillations of the total current in asymmetric discharges, but do not provide a kinetic interpretation of these effects. In fact, they leave important questions open: How is current continuity ensured in the presence of energetic electron beams generated by the expanding sheaths that lead to a local enhancement of the conduction current propagating through the bulk? How do the beam electrons interact with cold bulk electrons? What is the kinetic origin of resonance phenomena? Based on kinetic simulations, we find that the energetic beam electrons interact with cold bulk electrons (modulated on a timescale of the inverse local electron plasma frequency) via a time dependent electric field outside the sheaths. This electric field is caused by the electron beam itself, which leaves behind a positive space charge, that attracts cold bulk electrons towards the expanding sheath. The resulting displacement current ensures current continuity by locally compensating the enhancement of the conduction current. The backflow of cold electrons and their interaction with the nonlinear plasma sheath cause the generation of multiple electron beams during one phase of sheath expansion and contribute to a strongly non-sinusoidal RF current. These kinetic mechanisms are the basis for a fundamental understanding of the electron power absorption dynamics and resonance phenomena in such plasmas, which are found to occur in discharges of different symmetries including perfectly symmetric plasmas.
Kinetic dissipation and anisotropic heating in a turbulent collisionless plasma
Parashar, T. N.; Shay, M. A.; Cassak, P. A.; Matthaeus, W. H.
2009-03-15
The kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations. In magnetohydrodynamics (MHD) this configuration leads rapidly to broadband turbulence. At large length scales, the evolution of the hybrid simulations is very similar to MHD, with magnetic power spectra displaying scaling similar to a Kolmogorov scaling of -5/3. At small scales, differences from MHD arise, as energy dissipates into heat almost exclusively through the magnetic field. The magnetic energy spectrum of the hybrid simulation shows a break where linear theory predicts that the Hall term in Ohm's law becomes significant, leading to dispersive kinetic Alfven waves. A key result is that protons are heated preferentially in the plane perpendicular to the mean magnetic field, creating a proton temperature anisotropy of the type observed in the corona and solar wind.
Plasma kinetics of ethanol conversion in a glow discharge
NASA Astrophysics Data System (ADS)
Levko, D. S.; Tsymbalyuk, A. N.; Shchedrin, A. I.
2012-11-01
The mechanism of ethanol conversion in a nonequilibrium glow discharge has been studied. It is shown that molecular hydrogen is produced in reactions between ethanol molecules and hydrogen atoms in the initial stage and in reactions involving active H, CH2OH, CH3CHOH, and formaldehyde in the final stage. Comparison with experimental data shows that the kinetic mechanism used in these calculations correctly predicts the concentrations of the main components of the gas mixture.
Analytical solution of the kinetic equation for a uniform plasma in a magnetic field
Ji, Jeong-Young; Held, Eric D.
2010-07-15
The kinetic equation for a single-component uniform plasma in a magnetic field is analytically solved by the moment method. The linear system of ordinary differential equations for the moments is decomposed into subsystems of lower dimensions by a geometric method. The eigensystem of each subsystem shows that parallel moments decay monotonically, but perpendicular lth harmonic moments decay while oscillating with the l,l-2,...,-th harmonics of gyrofrequency. A generalization to a multicomponent plasma is discussed.
Continuum kinetic modeling of the tokamak plasma edge
NASA Astrophysics Data System (ADS)
Dorf, M. A.; Dorr, M. R.; Hittinger, J. A.; Cohen, R. H.; Rognlien, T. D.
2016-05-01
The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.
Lauber, Ph. Guenter, S.; Koenies, A.; Pinches, S.D.
2007-09-10
In a plasma with a population of super-thermal particles generated by heating or fusion processes, kinetic effects can lead to the additional destabilisation of MHD modes or even to additional energetic particle modes. In order to describe these modes, a new linear gyrokinetic MHD code has been developed and tested, LIGKA (linear gyrokinetic shear Alfven physics) [Ph. Lauber, Linear gyrokinetic description of fast particle effects on the MHD stability in tokamaks, Ph.D. Thesis, TU Muenchen, 2003; Ph. Lauber, S. Guenter, S.D. Pinches, Phys. Plasmas 12 (2005) 122501], based on a gyrokinetic model [H. Qin, Gyrokinetic theory and computational methods for electromagnetic perturbations in tokamaks, Ph.D. Thesis, Princeton University, 1998]. A finite Larmor radius expansion together with the construction of some fluid moments and specification to the shear Alfven regime results in a self-consistent, electromagnetic, non-perturbative model, that allows not only for growing or damped eigenvalues but also for a change in mode-structure of the magnetic perturbation due to the energetic particles and background kinetic effects. Compared to previous implementations [H. Qin, mentioned above], this model is coded in a more general and comprehensive way. LIGKA uses a Fourier decomposition in the poloidal coordinate and a finite element discretisation in the radial direction. Both analytical and numerical equilibria can be treated. Integration over the unperturbed particle orbits is performed with the drift-kinetic HAGIS code [S.D. Pinches, Ph.D. Thesis, The University of Nottingham, 1996; S.D. Pinches et al., CPC 111 (1998) 131] which accurately describes the particles' trajectories. This allows finite-banana-width effects to be implemented in a rigorous way since the linear formulation of the model allows the exchange of the unperturbed orbit integration and the discretisation of the perturbed potentials in the radial direction. Successful benchmarks for toroidal Alfven
Exact kinetic theory for the instability of an electron beam in a hot magnetized plasma
Timofeev, I. V.; Annenkov, V. V.
2013-09-15
Efficiency of collective beam-plasma interaction strongly depends on the growth rates of dominant instabilities excited in the system. Nevertheless, exact calculations of the full unstable spectrum in the framework of relativistic kinetic theory for arbitrary magnetic fields and particle distributions were unknown until now. In this paper, we give an example of such a calculation answering the question whether the finite thermal spreads of plasma electrons are able to suppress the fastest growing modes in the beam-plasma system. It is shown that nonrelativistic temperatures of Maxwellian plasmas can stabilize only the oblique instabilities of relativistic beam. On the contrary, non-Maxwellian tails typically found in laboratory beam-plasma experiments are able to substantially reduce the growth rate of the dominant longitudinal modes affecting the efficiency of turbulent plasma heating.
Test-particle method in kinetic theory of a plasma.
NASA Technical Reports Server (NTRS)
Matsuda, K.
1971-01-01
The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.
NASA Astrophysics Data System (ADS)
Kim, June Young; Cho, Won-Hwi; Dang, Jeong-Jeung; Chung, Kyoung-Jae; Hwang, Y. S.
2016-02-01
Electron kinetics regime is characterized with the evolution of electron energy probability functions (EEPFs) in inductively coupled hydrogen plasmas. Measurements on EEPFs are carried out with a radio-frequency-compensated single Langmuir probe at the center of a planar-type hydrogen plasma driven by 13.56 MHz wave frequency. Measured EEPFs deviate considerably from the Maxwellian distribution only at relatively high pressures (15-40 mTorr), and the effective electron temperature steeply decreases as the gas pressure increases. Such evolution of the EEPF shapes with pressures is discussed in the consideration of the electron energy relaxation length and various characteristic frequencies. It is found that the EEPFs show locally depleted electron energy distribution where the electron-molecule vibrational collision frequency exceeds the electron-electron collision frequency at the local kinetics regime, while the measured EEPF is not dependent on the vibrational collision frequency at the non-local kinetics regime. Variation of the EEPF shape with distance from the heating region at the local kinetics regime is also well explained in the context of the energy relaxation length and electron-molecule collision frequencies. This study indicates that the control of electron energy distribution should be carried out in the consideration of electron kinetic regime depending on the energy relaxation length for various hydrogen plasma sources.
Diomede, P.; Longo, S.; Capitelli, M.
2005-05-16
We present a 1D(r)2D(v) particle code for capacitively coupled radio frequency discharge plasmas in hydrogen, which includes a rigorous kinetic modeling of ion transport and several solutions to speed up the convergence. In a test case the effect of surface atom recombination and molecule vibrational deactivation on H- concentration is investigated.
Comparisons of dense-plasma-focus kinetic simulations with experimental measurements
Schmidt, A.; Link, A.; Welch, D.; Ellsworth, J.; Falabella, S.; Tang, V.
2014-06-01
Dense-plasma-focus (DPF) Z-pinch devices are sources of copious high-energy electrons and ions, x rays, and neutrons. The mechanisms through which these physically simple devices generate such high-energy beams in a relatively short distance are not fully understood and past optimization efforts of these devices have been largely empirical. Previously we reported on fully kinetic simulations of a DPF and compared them with hybrid and fluid simulations of the same device. Here we present detailed comparisons between fully kinetic simulations and experimental data on a 1.2 kJ DPF with two electrode geometries, including neutron yield and ion beam energy distributions. A more intensive third calculation is presented which examines the effects of a fully detailed pulsed power driver model. We also compare simulated electromagnetic fluctuations with direct measurement of radiofrequency electromagnetic fluctuations in a DPF plasma. These comparisons indicate that the fully kinetic model captures the essential physics of these plasmas with high fidelity, and provide further evidence that anomalous resistivity in the plasma arises due to a kinetic instability near the lower hybrid frequency.
Non-thermal plasma destruction of allyl alcohol in waste gas: kinetics and modelling
NASA Astrophysics Data System (ADS)
DeVisscher, A.; Dewulf, J.; Van Durme, J.; Leys, C.; Morent, R.; Van Langenhove, H.
2008-02-01
Non-thermal plasma treatment is a promising technique for the destruction of volatile organic compounds in waste gas. A relatively unexplored technique is the atmospheric negative dc multi-pin-to-plate glow discharge. This paper reports experimental results of allyl alcohol degradation and ozone production in this type of plasma. A new model was developed to describe these processes quantitatively. The model contains a detailed chemical degradation scheme, and describes the physics of the plasma by assuming that the fraction of electrons that takes part in chemical reactions is an exponential function of the reduced field. The model captured the experimental kinetic data to less than 2 ppm standard deviation.
High-order continuum kinetic method for modeling plasma dynamics in phase space
Vogman, G. V.; Colella, P.; Shumlak, U.
2014-12-15
Continuum methods offer a high-fidelity means of simulating plasma kinetics. While computationally intensive, these methods are advantageous because they can be cast in conservation-law form, are not susceptible to noise, and can be implemented using high-order numerical methods. Advances in continuum method capabilities for modeling kinetic phenomena in plasmas require the development of validation tools in higher dimensional phase space and an ability to handle non-cartesian geometries. To that end, a new benchmark for validating Vlasov-Poisson simulations in 3D (x,vx,vy) is presented. The benchmark is based on the Dory-Guest-Harris instability and is successfully used to validate a continuum finite volumemore » algorithm. To address challenges associated with non-cartesian geometries, unique features of cylindrical phase space coordinates are described. Preliminary results of continuum kinetic simulations in 4D (r,z,vr,vz) phase space are presented.« less
Comparative study of gyrokinetic, hybrid-kinetic and fully kinetic wave physics for space plasmas
NASA Astrophysics Data System (ADS)
Told, D.; Cookmeyer, J.; Muller, F.; Astfalk, P.; Jenko, F.
2016-06-01
A set of numerical solvers for the linear dispersion relations of the gyrokinetic (GK), the hybrid-kinetic (HK), and the fully kinetic (FK) model is employed to study the physics of the KAW and the fast magnetosonic mode in these models. In particular, we focus on parameters that are relevant for solar wind oriented applications (using a homogeneous, isotropic background), which are characterized by wave propagation angles averaging close to 90°. It is found that the GK model, while lacking high-frequency solutions and cyclotron effects, faithfully reproduces the FK {{Alfv\\acute{e}n}} wave physics close to, and sometimes significantly beyond, the boundaries of its range of validity. The HK model, on the other hand, is much more complete in terms of high-frequency waves, but owing to its simple electron model it is found to severely underpredict wave damping rates even on ion spatial scales across a large range of parameters, despite containing full kinetic ion physics.
NASA Astrophysics Data System (ADS)
Mašek, Martin; Rohlena, Karel
2015-05-01
Influence of kinetic effects on 3-wave interaction was examined within the frame of stimulated Raman backward scattering (SRBS) in a rarefied laser corona. The plasma is supposed to be weakly collisional with a negligible density gradient. The model is centred on the physical situation of shock ignition at a large scale direct drive compression experiments. The modelling uses a 1D geometry in a Maxwell-Vlasov model. The method used is a truncated Fourier-Hermite expansion numerically stabilized by a model collisional term with a realistic value of the collision frequency. In parallel, besides the linear theory of SRBS, a coupled mode 3-wave equation system (laser driving wave, Raman back-scattered wave and the daughter forward scattered plasma wave) is solved to demonstrate the correspondence between the full kinetic model and 3-wave interaction with no electron kinetics involved to identify the differences between both the solutions arising due to the electron kinetic effects. We concentrated mainly on the Raman reflectivity, which is one of the important parameters controlling the efficiency of the shock ignition scheme. It was found that the onset of the kinetic effects has a distinct intensity threshold, above which the Raman reflectivity may go down due to the electron kinetics. In addition, we were trying to identify the most important features of the electron phase space behaviour, such as particle trapping in potential minima of the generated plasma wave and its consequences for the 3-wave interaction. The role of the trapped electrons seems to be crucial for a deformation of the plasma wave dispersion curve, as indicated in some earlier work.
One-dimensional hybrid-direct kinetic simulation of the discharge plasma in a Hall thruster
Hara, Kentaro; Boyd, Iain D.; Kolobov, Vladimir I.
2012-11-15
In order to model the non-equilibrium plasma within the discharge region of a Hall thruster, the velocity distribution functions (VDFs) must be obtained accurately. A direct kinetic (DK) simulation method that directly solves the plasma Boltzmann equation can achieve better resolution of VDFs in comparison to particle simulations, such as the particle-in-cell (PIC) method that inherently include statistical noise. In this paper, a one-dimensional hybrid-DK simulation, which uses a DK simulation for heavy species and a fluid model for electrons, is developed and compared to a hybrid-PIC simulation. Time-averaged results obtained from the hybrid-DK simulation are in good agreement with hybrid-PIC results and experimental data. It is shown from a comparison of using a kinetic simulation and solving the continuity equation that modeling of the neutral atoms plays an important role for simulations of the Hall thruster discharge plasma. In addition, low and high frequency plasma oscillations are observed. Although the kinetic nature of electrons is not resolved due to the use of a fluid model, the hybrid-DK model provides spatially and temporally well-resolved plasma properties and an improved resolution of VDFs for heavy species with less statistical noise in comparison to the hybrid-PIC method.
Transition from gas to plasma kinetic equilibria in gravitating axisymmetric structures
Cremaschini, Claudio; Stuchlík, Zdeněk
2014-04-15
The problem of the transition from gas to plasma in gravitating axisymmetric structures is addressed under the assumption of having initial and final states realized by kinetic Maxwellian-like equilibria. In astrophysics, the theory applies to accretion-disc scenarios around compact objects. A formulation based on non-relativistic kinetic theory for collisionless systems is adopted. Equilibrium solutions for the kinetic distribution functions describing the initial neutral matter and the resulting plasma state are constructed in terms of single-particle invariants and expressed by generalized Maxwellian distributions. The final plasma configuration is related to the initial gas distribution by the introduction of appropriate functional constraints. Qualitative aspects of the solution are investigated and physical properties of the system are pointed out. In particular, the admitted functional dependences of the fluid fields carried by the corresponding equilibrium distributions are determined. Then, the plasma is proved to violate the condition of quasi-neutrality, implying a net charge separation between ions and electrons. This result is shown to be independent of the precise realization of the plasma distribution function, while a physical mechanism able to support a non-neutral equilibrium state is proposed.
NASA Astrophysics Data System (ADS)
Kubes, P.; Prykarpatsky, A. K.; Zagrodzinski, J.; Prykarpatsky, Y. A.
In this article we will follow the approach developed in articles N.~N.~Bogoliubov, V.~Hr.~Samoilenko, Ukr. Fiz. Zh., 37, 147 (1992); J.~Gibbon, Physica D, 3, 503 (1981) using modern Lie--algebraic and symplectic geometry methods. It is devoted to the description of Boltzman--Vlasov type kinetic equations and some two--dimensional hydrodynamic Benney type flows associated with them. In our case of the cylindrical symmetry taking place at the interrupted magnetic z--pinch in plasma we used intensively the corresponding two--dimensionality of the plasma flow under consideration which made it possible to build a kinetic model of the plasmoid vortex structure with a conserved number of linkages of vortex lines. The latter can be used to explain the observed earlier stability of the plasmoid structure at the magnetic z--pinch.
Plasma-resistivity-induced strong damping of the kinetic resistive wall mode.
He, Yuling; Liu, Yueqiang; Liu, Yue; Hao, Guangzhou; Wang, Aike
2014-10-24
An energy-principle-based dispersion relation is derived for the resistive wall mode, which incorporates both the drift kinetic resonance between the mode and energetic particles and the resistive layer physics. The equivalence between the energy-principle approach and the resistive layer matching approach is first demonstrated for the resistive plasma resistive wall mode. As a key new result, it is found that the resistive wall mode, coupled to the favorable average curvature stabilization inside the resistive layer (as well as the toroidal plasma flow), can be substantially more stable than that predicted by drift kinetic theory with fast ion stabilization, but with the ideal fluid assumption. Since the layer stabilization becomes stronger with decreasing plasma resistivity, this regime is favorable for reactor scale, high-temperature fusion devices. PMID:25379920
Kinetic Alfven solitary waves in a magnetized plasma with superthermal electrons
Panwar, A. E-mail: ryu201@postech.ac.kr Ryu, C. M. E-mail: ryu201@postech.ac.kr; Bains, A. S. E-mail: ryu201@postech.ac.kr
2015-09-15
A study of the ion Larmor radius effects on the solitary kinetic Alfven waves (SKAWs) in a magnetized plasma with superthermal electrons is presented by employing the kinetic theory. The linear dispersion relation of SKAW is shown to depend on the superthermal parameter κ, ion to electron temperature ratio, and the angle of wave propagation. Using the Sagdeev potential approach, the energy balance equation has been derived to study the dynamics of SKAWs. The effects of various plasma parameters are investigated for the propagation of SKAWs. It is shown that only compressive solitons can exist and in the Maxwellian limit our results are in good agreement with previous studies. Further, the characteristics of small amplitude SKAWs are investigated. Present study could be useful for the understanding of SKAWs in a low β plasma in astrophysical environment, where particle distributions are superthermal in nature.
Kinetic shear Alfvén instability in the presence of impurity ions in tokamak plasmas
Lu, Gaimin; Shen, Y.; Xie, T.; He, Zhixiong; He, Hongda; Qi, Longyu; Cui, Shaoyan
2013-10-15
The effects of impurity ions on the kinetic shear Alfvén (KSA) instability in tokamak plasmas are investigated by numerically solving the integral equations for the KSA eigenmode in the toroidal geometry. The kinetic effects of hydrogen and impurity ions, including transit motion, finite ion Larmor radius, and finite-orbit-width, are taken into account. Toroidicity induced linear mode coupling is included through the ballooning-mode representation. Here, the effects of carbon, oxygen, and tungsten ions on the KSA instability in toroidal plasmas are investigated. It is found that, depending on the concentration and density profile of the impurity ions, the latter can be either stabilizing or destabilizing for the KSA modes. The results here confirm the importance of impurity ions in tokamak experiments and should be useful for analyzing experimental data as well as for understanding anomalous transport and control of tokamak plasmas.
Orszag Tang vortex - Kinetic study of a turbulent plasma
Parashar, T. N.; Servidio, S.; Shay, M. A.; Matthaeus, W. H.; Cassak, P. A.
2010-03-25
Kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations based on particle in cell ions and fluid electrons. In magnetohydrodynamics (MHD) this configuration leads rapidly to broadband turbulence. An earlier study estimated the dissipation in the system. A comparison of MHD and hybrid simulations showed similar behavior at large scales but substantial differences at small scales. The hybrid magnetic energy spectrum shows a break at the scale where Hall term in the Ohm's law becomes important. The protons heat perpendicularly and most of the energy is dissipated through magnetic interactions. Here, the space time structure of the system is studied using frequency-wavenumber (k-omega) decomposition. No clear resonances appear, ruling out the cyclotron resonances as a likely candidate for the perpendicular heating. The only distinguishable wave modes present, which constitute a small percentage of total energy, are magnetosonic modes.
Kinetic theory of a two-dimensional magnetized plasma. II - Balescu-Lenard limit.
NASA Technical Reports Server (NTRS)
Vahala, G.
1972-01-01
The kinetic theory of a two-dimensional one-species plasma in a uniform dc magnetic field is investigated in the small plasma parameter limit. The plasma consists of charged rods interacting through the logarithmic Coulomb potential. Vahala and Montgomery earlier (1971) derived a Fokker-Planck equation for this system, but it contained a divergent integral, which had to be cut off on physical grounds. This cutoff is compared to the standard cutoff introduced in the two-dimensional unmagnetized Fokker-Planck equation. In the small plasma parameter limit, it is shown that the Balescu-Lenard collision term is zero in the long time average limit if only two-body interactions are considered. The energy transfer from a test particle to an equilibrium plasma is discussed and is also shown to be zero in the long time average limit. This supports the unexpected result of zero Balescu-Lenard collision term.
Non-LTE Steady-State Kinetics of He-Air Atmospheric Pressure Plasmas
NASA Astrophysics Data System (ADS)
Petrova, Tzvetelina; Petrov, George; Gillman, Eric; Boris, David; Hernández, Sandra; Walton, Scott
2015-11-01
A non-LTE, steady-state collisional-radiative kinetics model is developed to study discharges produced in mixtures of He, N2 and O2 (He-Air) at atmospheric pressures. The model is based on a self-consistent solution of coupled Boltzmann equation for the electron energy distribution function, electron energy balance equation, gas thermal balance equation, and a system of non-linear equations for species that govern plasma chemistry (electrons, ions, radicals, atoms and molecules in ground and excited states). The model and results can be applied to study a variety of atmospheric pressure plasmas generated in He-Air mixtures, such as plasma jets, dielectric barrier discharges, laser-induced plasmas, microwave plasmas, etc. In this talk, collisional rates and species densities are obtained as a function of He-to-air ratio and the results are benchmarked against available experimental data. Work supported by the NRL Base Program.
Liu, Chang; Dodin, Ilya Y.
2015-08-15
The nonlinear frequency shift is derived in a transparent asymptotic form for intense Langmuir waves in general collisionless plasma. The formula describes both fluid and kinetic effects simultaneously. The fluid nonlinearity is expressed, for the first time, through the plasma dielectric function, and the kinetic nonlinearity accounts for both smooth distributions and trapped-particle beams. Various known limiting scalings are reproduced as special cases. The calculation avoids differential equations and can be extended straightforwardly to other nonlinear plasma waves.
Energy conserving continuum algorithms for kinetic & gyrokinetic simulations of plasmas
NASA Astrophysics Data System (ADS)
Hakim, A.; Hammett, G. W.; Shi, E.; Stoltzfus-Dueck, T.
2015-11-01
We present high-order, energy conserving, continuum algorithms for the solution of gyrokinetic equations for use in edge turbulence simulations. The distribution function is evolved with a discontinuous Galerkin scheme, while the fields are evolved with a continuous finite-element method. These algorithms work for a general, possibly non-canonical, Poisson bracket operator and conserve energy exactly. Benchmark simulations with ETG turbulence in 3X/2V are shown, as well as initial applications of the algorithms to turbulence in a simplified SOL geometry. Sheath boundary conditions with recycling and secondary electron emission are implemented, and a Lenard-Bernstein collision operator is included. Extension of these algorithms to full Vlasov-Maxwell equations are presented. It is shown that with a particular choice of numerical fluxes the total (particle+field) energy is conserved. Algorithms are implemented in a flexible and open-source framework, Gkeyll, which also includes fluid models, allowing potential hybrid simulations of various plasma problems. Supported by the Max-Planck/Princeton Center for Plasma Physics, and DOE Contract DE-AC02-09CH11466.
Two dimensional kinetic analysis of electrostatic harmonic plasma waves
NASA Astrophysics Data System (ADS)
Fonseca-Pongutá, E. C.; Ziebell, L. F.; Gaelzer, R.; Yoon, P. H.
2016-06-01
Electrostatic harmonic Langmuir waves are virtual modes excited in weakly turbulent plasmas, first observed in early laboratory beam-plasma experiments as well as in rocket-borne active experiments in space. However, their unequivocal presence was confirmed through computer simulated experiments and subsequently theoretically explained. The peculiarity of harmonic Langmuir waves is that while their existence requires nonlinear response, their excitation mechanism and subsequent early time evolution are governed by essentially linear process. One of the unresolved theoretical issues regards the role of nonlinear wave-particle interaction process over longer evolution time period. Another outstanding issue is that existing theories for these modes are limited to one-dimensional space. The present paper carries out two dimensional theoretical analysis of fundamental and (first) harmonic Langmuir waves for the first time. The result shows that harmonic Langmuir wave is essentially governed by (quasi)linear process and that nonlinear wave-particle interaction plays no significant role in the time evolution of the wave spectrum. The numerical solutions of the two-dimensional wave spectra for fundamental and harmonic Langmuir waves are also found to be consistent with those obtained by direct particle-in-cell simulation method reported in the literature.
Electromagnetic fluctuations in magnetized plasmas. I. The rigorous relativistic kinetic theory
Schlickeiser, R. E-mail: yoonp@umd.edu; Yoon, P. H. E-mail: yoonp@umd.edu
2015-07-15
Using the system of the Klimontovich and Maxwell equations, the general linear fluctuation theory for magnetized plasmas is developed. General expressions for the electromagnetic fluctuation spectra (electric and magnetic fields) from uncorrelated plasma particles in plasmas with a uniform magnetic field are derived, which are covariantly correct within the theory of special relativity. The general fluctuation spectra hold for plasmas of arbitrary composition, arbitrary momentum dependences of the plasma particle distribution functions, and arbitrary orientations of the wave vector with respect to the uniform magnetic field. Moreover, no restrictions on the values of the real and the imaginary parts of the frequency are made. The derived fluctuation spectra apply to both non-collective fluctuations and collective plasma eigenmodes in magnetized plasmas. In the latter case, kinetic equations for the components of fluctuating electric and magnetic fields in magnetized plasmas are derived that include the effect of spontaneous emission and absorption. In the limiting case of an unmagnetized plasmas, the general fluctuation spectra correctly reduce to the unmagnetized fluctuation spectra derived before.
Bistoletti, M; Alvarez, L; Lanusse, C; Moreno, L
2014-01-01
1. To optimise the use of albendazole (ABZ) as an anthelmintic in hens, the effects of fasting and type of diet on the plasma kinetics of ABZ and its metabolites were evaluated. 2. Twenty-four hens were distributed into 4 groups: In experiment I the Fed group were fed ad libitum, while the Fasted group was fasted over a 12-h period. In experiment II the Pelleted group was fed with pelleted commercial food, while the Grain group was fed with cereal grains. All the groups were treated with ABZ by oral route. Blood samples were taken and plasma analysed by HPLC. 3. ABZ and its metabolites albendazole-sulphoxide (ABZSO) and albendazole-sulphone (ABZSO2) were recovered in plasma in all the groups. The 12-h fasting period did not modify the disposition kinetics of ABZ in hens. The type of feed affected ABZ kinetics. ABZSO concentration profile was higher and detected for longer in the Grain group compared to the Pelleted group. Statistical differences were not found for AUC0-∞ values, whereas the T1/2for and T1/2el were different between groups. 4. Factors affecting ABZ kinetic behaviour should be taken into account to optimise its use to ensure the sustainability of the limited available anthelmintic therapeutic tools in avian parasite control. PMID:25159169
Bendib, A.
2008-09-23
The conference is devoted to the study of systems consisting of a large number of particles by using the kinetic theory. In a first part, we present a general overview of the kinetic theory. In particular, the role of the correlations between particles is shown and discussed through the main models reported in the literature. In a second part, we present three applications to the transport properties in plasmas and neutral gases. The first application is devoted to the transport in hot plasmas perturbed with respect to the global equilibrium. The quasi-static and collisionless distribution function and transport coefficients are established. The influence of relativistic effects is also discussed. The second application deals with strongly inhomogeneous magnetized plasmas. The transport coefficients of Braginskii are calculated numerically in the local and the weakly nonlocal approximations. New nonlocal transport coefficients are emphasized. Finally, we apply the kinetic theory to the neutral gases by calculating the semi-collisional dispersion relation of acoustic waves. In particular, the dispersion and the damping of these waves in rarefied gases are highlighted. The method used to solve the kinetic equations is compared with the conventional method of Chapman-Enskog.
Kinetic theory of a two-dimensional magnetized plasma.
NASA Technical Reports Server (NTRS)
Vahala, G.; Montgomery, D.
1971-01-01
Several features of the equilibrium and nonequilibrium statistical mechanics of a two-dimensional plasma in a uniform dc magnetic field are investigated. The charges are assumed to interact only through electrostatic potentials. The problem is considered both with and without the guiding-center approximation. With the guiding-center approximation, an appropriate Liouville equation and BBGKY hierarchy predict no approach to thermal equilibrium for the spatially uniform case. For the spatially nonuniform situation, a guiding-center Vlasov equation is discussed and solved in special cases. For the nonequilibrium, nonguiding-center case, a Boltzmann equation, and a Fokker-Planck equation are derived in the appropriate limits. The latter is more tractable than the former, and can be shown to obey conservation laws and an H-theorem, but contains a divergent integral which must be cut off on physical grounds. Several unsolved problems are posed.
Continuum kinetic modeling of the tokamak plasma edge
Dorf, M. A.; Dorr, M.; Rognlien, T.; Hittinger, J.; Cohen, R.
2016-03-10
In this study, the first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasmatransport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalousmore » radial transport.« less
Dust kinetic Alfvén solitary and rogue waves in a superthermal dusty plasma
NASA Astrophysics Data System (ADS)
Saini, N. S.; Singh, Manpreet; Bains, A. S.
2015-11-01
Dust kinetic Alfvén solitary waves (DKASWs) have been examined in a low-β dusty plasma comprising of negatively charged dust grains, superthermal electrons, and ions. A nonlinear Korteweg-de Vries (KdV) equation has been derived using the reductive perturbation method. The combined effects of superthermality of charged particles (via κ), plasma β, obliqueness of propagation (θ), and dust concentration (via f) on the shape and size of the DKASWs have been examined. Only negative potential (rarefactive) structures are observed. Further, characteristics of dust kinetic Alfvén rogue waves (DKARWs), by deriving the non-linear Schrödinger equation (NLSE) from the KdV equation, are studied. Rational solutions of NLSE show that rogue wave envelopes are supported by this plasma model. It is observed that the influence of various plasma parameters (superthermality, plasma β, obliqueness, and dust concentration) on the characteristics of the DKARWs is very significant. This fundamental study may be helpful in understanding the formation of coherent nonlinear structures in space and astrophysical plasma environments where superthermal particles are present.
Dust kinetic Alfvén solitary and rogue waves in a superthermal dusty plasma
Saini, N. S. Singh, Manpreet; Bains, A. S.
2015-11-15
Dust kinetic Alfvén solitary waves (DKASWs) have been examined in a low-β dusty plasma comprising of negatively charged dust grains, superthermal electrons, and ions. A nonlinear Korteweg-de Vries (KdV) equation has been derived using the reductive perturbation method. The combined effects of superthermality of charged particles (via κ), plasma β, obliqueness of propagation (θ), and dust concentration (via f) on the shape and size of the DKASWs have been examined. Only negative potential (rarefactive) structures are observed. Further, characteristics of dust kinetic Alfvén rogue waves (DKARWs), by deriving the non-linear Schrödinger equation (NLSE) from the KdV equation, are studied. Rational solutions of NLSE show that rogue wave envelopes are supported by this plasma model. It is observed that the influence of various plasma parameters (superthermality, plasma β, obliqueness, and dust concentration) on the characteristics of the DKARWs is very significant. This fundamental study may be helpful in understanding the formation of coherent nonlinear structures in space and astrophysical plasma environments where superthermal particles are present.
Tarvainen, O. Laulainen, J.; Komppula, J.; Kronholm, R.; Kalvas, T.; Koivisto, H.; Izotov, I.; Mansfeld, D.; Skalyga, V.
2015-02-15
Electron cyclotron resonance ion source (ECRIS) plasmas are prone to kinetic instabilities due to anisotropy of the electron energy distribution function stemming from the resonant nature of the electron heating process. Electron cyclotron plasma instabilities are related to non-linear interaction between plasma waves and energetic electrons resulting to strong microwave emission and a burst of energetic electrons escaping the plasma, and explain the periodic oscillations of the extracted beam currents observed in several laboratories. It is demonstrated with a minimum-B 14 GHz ECRIS operating on helium, oxygen, and argon plasmas that kinetic instabilities restrict the parameter space available for the optimization of high charge state ion currents. The most critical parameter in terms of plasma stability is the strength of the solenoid magnetic field. It is demonstrated that due to the instabilities the optimum B{sub min}-field in single frequency heating mode is often ≤0.8B{sub ECR}, which is the value suggested by the semiempirical scaling laws guiding the design of modern ECRISs. It is argued that the effect can be attributed not only to the absolute magnitude of the magnetic field but also to the variation of the average magnetic field gradient on the resonance surface.
Arbitrary amplitude kinetic Alfven solitary waves in two temperature electron superthermal plasma
NASA Astrophysics Data System (ADS)
Singh, Manpreet; Singh Saini, Nareshpal; Ghai, Yashika
2016-07-01
Through various satellite missions it is observed that superthermal velocity distribution for particles is more appropriate for describing space and astrophysical plasmas. So it is appropriate to use superthermal distribution, which in the limiting case when spectral index κ is very large ( i.e. κ→∞), shifts to Maxwellian distribution. Two temperature electron plasmas have been observed in auroral regions by FAST satellite mission, and also by GEOTAIL and POLAR satellite in the magnetosphere. Kinetic Alfven waves arise when finite Larmor radius effect modifies the dispersion relation or characteristic perpendicular wavelength is comparable to electron inertial length. We have studied the kinetic Alfven waves (KAWs) in a plasma comprising of positively charged ions, superthermal hot electrons and Maxwellian distributed cold electrons. Sagdeev pseudo-potential has been employed to derive an energy balance equation. The critical Mach number has been determined from the expression of Sagdeev pseudo-potential to see the existence of solitary structures. It is observed that sub-Alfvenic compressive solitons and super-Alfvenic rarefactive solitons exist in this plasma model. It is also observed that various parameters such as superthermality of hot electrons, relative concentration of cold and hot electron species, Mach number, plasma beta, ion to cold electron temperature ratio and ion to hot electron temperature ratio have significant effect on the amplitude and width of the KAWs. Findings of this investigation may be useful to understand the dynamics of coherent non-linear structures (i.e. KAWs) in space and astrophysical plasmas.
A Kinetic Plasma-Pumped Rare Gas Laser
NASA Astrophysics Data System (ADS)
Parsey, Guy; Güçlü, Yaman; Verboncoeur, John; Christlieb, Andrew
2015-09-01
Extending from diode-pumped alkali vapor lasers (DPAL), Han and Heaven have shown that rare gas metastable states, np5 (n + 1) s[ 3 / 2 ] 2 , can operate as the base of a three-level laser with excition of the (n + 1) s --> (n + 1) p transitions. Though both the rare gas lasers (RGL) and DPALs can be excited with incoherent optical pumping, RGLs do not suffer from the highly reactive behavior of alkali metals. Since metastable populations are maintained via electric discharge, we propose using a tuned electron energy distribution function (EEDF) to modify RGL efficiencies and drive the population inversion. The EEDF is maintained by the discharge along with the introduction of electron sources. Using our kinetic global modeling framework (KGMf) and three gas systems (helium buffered argon and krypton along with pure argon), we first validate the intracavity intensity laser model and then generate gain and energy efficiency baselines for each system. Parameter scanning methods are then used to find optimized EEDFs and system parameters for metastable production, generation of a lasing population inversion, and increasing RGL operation efficiencies. Finally, we determine if an RGL can operate without optical pumping. Supported by AFOSR and a MSU Strategic Partnership Grant
Phase-space description of plasma waves. Part 1. Linear theory
NASA Astrophysics Data System (ADS)
Biro, T.; Rönnmark, K.
1992-06-01
We develop an (r, k) phase-space description of waves in plasmas by introducing Gaussian window functions to separate short-scale oscillations from long-scale modulations of the wave fields and variations in the plasma parameters. To obtain a wave equation that unambiguously separates conservative dynamics from dissipation in an inhomogeneous and time-varying background plasma, we first discuss the proper form of the current response function. In analogy with the particle distribution function f(v, r, t), we introduce a wave density N(k, r, t) on phase space. This function is proved to satisfy a simple continuity equation. Dissipation is also included, and this allows us to describe the damping or growth of wave density along rays. Problems involving geometric optics of continuous media often appear simpler when viewed in phase space, since the flow of N in phase space is incompressible.
Quantification of Transthyretin Kinetic Stability in Human Plasma Using Subunit Exchange
2015-01-01
The transthyretin (TTR) amyloidoses are a group of degenerative diseases caused by TTR aggregation, requiring rate-limiting tetramer dissociation. Kinetic stabilization of TTR, by preferential binding of a drug to the native tetramer over the dissociative transition state, dramatically slows the progression of familial amyloid polyneuropathy. An established method for quantifying the kinetic stability of recombinant TTR tetramers in buffer is subunit exchange, in which tagged TTR homotetramers are added to untagged homotetramers at equal concentrations to measure the rate at which the subunits exchange. Herein, we report a subunit exchange method for quantifying the kinetic stability of endogenous TTR in human plasma. The subunit exchange reaction is initiated by the addition of a substoichiometric quantity of FLAG-tagged TTR homotetramers to endogenous TTR in plasma. Aliquots of the subunit exchange reaction, taken as a function of time, are then added to an excess of a fluorogenic small molecule, which immediately arrests further subunit exchange. After binding, the small molecule reacts with the TTR tetramers, rendering them fluorescent and detectable in human plasma after subsequent ion exchange chromatography. The ability to report on the extent of TTR kinetic stabilization resulting from treatment with oral tafamidis is important, especially for selection of the appropriate dose for patients carrying rare mutations. This method could also serve as a surrogate biomarker for the prediction of the clinical outcome. Subunit exchange was used to quantify the stabilization of WT TTR from senile systemic amyloidosis patients currently being treated with tafamidis (20 mg orally, once daily). TTR kinetic stability correlated with the tafamidis plasma concentration. PMID:24661308
Hydrodynamic description of an unmagnetized plasma with multiple ion species. I. General formulation
Simakov, Andrei Nikolaevich; Molvig, Kim
2016-03-17
A generalization of the Braginskii ion fluid description [S. I. Braginskii, Sov. Phys. JETP 6, 358 (1958)] to the case of an unmagnetized collisional plasma with multiple ion species is presented. An asymptotic expansion in the ion Knudsen number is used to derive the individual ion species continuity, as well as the total ion mass density, momentum, and energy evolution equations accurate through the second order. Expressions for the individual ion species drift velocities with respect to the center of mass reference frame, as well as for the total ion heat flux and viscosity, which are required to close themore » fluid equations, are evaluated in terms of the first-order corrections to the lowest order Maxwellian ion velocity distribution functions. A variational formulation for evaluating such corrections and its relation to the plasma entropy are presented. Employing trial functions for the corrections, written in terms of expansions in generalized Laguerre polynomials, and maximizing the resulting functionals produces two systems of linear equations (for “vector” and “tensor” portions of the corrections) for the expansion coefficients. A general matrix formulation of the linear systems as well as expressions for the resulting transport fluxes are presented in forms convenient for numerical implementation. The general formulation is employed in the companion paper [A. N. Simakov and K. Molvig, Hydrodynamic description of an unmagnetized plasma with multiple ion species. II. Two and three ion species plasmas, submitted to Phys. Plasmas (2015)] to evaluate the individual ion drift velocities and the total ion heat flux and viscosity for specific cases of two and three ion species plasmas.« less
Capitelli, M.; De Pascale, O.; Diomede, P.; Gorse, C.; Longo, S.; Pagano, D.; Gicquel, A.; Hassouni, K.
2005-04-06
Different approaches to study vibrational kinetics coupled to electron one for modeling different kinds of negative ion sources are presented. In particular two types of sources are investigated. The first one is a classical negative ion source in which the plasma is generated by thermoemitted electrons; in the second one, electrons already present in the mixture are accelerated by an RF field to sufficiently high energy to ionize the gas molecules. For the first kind of ion source a new computational scheme is presented to couple heavy particle and electron kinetics. Moreover models developed for an RF inductive discharge and for a parallel plate discharge are described.
Kinetics and Statistical Behaviour of Iron Recovery from Red Mud using Plasma Arc Furnace
NASA Astrophysics Data System (ADS)
Rath, Swagat S.; Jayasankar, K.; Satapathy, Bijoy K.; Mishra, Barada K.; Mukherjee, Partha S.
2011-06-01
Study of the recovery of pig iron from dry red mud through plasma smelting has been done. The reductant used was graphite. Various fluxes were used for the formation of slag. The parameters varied were smelting time, basicity and concentration of reductant and fluxes. It was observed that the reduction kinetics improved progressively with the increase of all the variables up to certain level and on further increase the metal recovery decreased. The results were interpreted kinetically and rate determining step was evaluated. Multi Linear Regression Analysis studies were carried out. The variables were subjected to Principal Component Analysis. Four factors could explain 75.4% of the variables.
Kinetic simulation of edge instability in fusion plasmas
NASA Astrophysics Data System (ADS)
Fulton, Daniel Patrick
In this work, gyrokinetic simulations in edge plasmas of both tokamaks and field reversed. configurations (FRC) have been carried out using the Gyrokinetic Toroidal Code (GTC) and A New Code (ANC) has been formulated for cross-separatrix FRC simulation. In the tokamak edge, turbulent transport in the pedestal of an H-mode DIII-D plasma is. studied via simulations of electrostatic driftwaves. Annulus geometry is used and simulations focus on two radial locations corresponding to the pedestal top with mild pressure gradient and steep pressure gradient. A reactive trapped electron instability with typical ballooning mode structure is excited in the pedestal top. At the steep gradient, the electrostatic instability exhibits unusual mode structure, peaking at poloidal angles theta=+- pi/2. Simulations find this unusual mode structure is due to steep pressure gradients in the pedestal but not due to the particular DIII-D magnetic geometry. Realistic DIII-D geometry has a stabilizing effect compared to a simple circular tokamak geometry. Driftwave instability in FRC is studied for the first time using gyrokinetic simulation. GTC. is upgraded to treat realistic equilibrium calculated by an MHD equilibrium code. Electrostatic local simulations in outer closed flux surfaces find ion-scale modes are stable due to the large ion gyroradius and that electron drift-interchange modes are excited by electron temperature gradient and bad magnetic curvature. In the scrape-off layer (SOL) ion-scale modes are excited by density gradient and bad curvature. Collisions have weak effects on instabilities both in the core and SOL. Simulation results are consistent with density fluctuation measurements in the C-2 experiment using Doppler backscattering (DBS). The critical density gradients measured by the DBS qualitatively agree with the linear instability threshold calculated by GTC simulations. One outstanding critical issue in the FRC is the interplay between turbulence in the FRC. core
Kinetic study of the secondary plasma created in the ITER neutraliser
NASA Astrophysics Data System (ADS)
Dure, F.; Lifschitz, A.; Bretagne, J.; Maynard, G.; Katsonis, K.; Simonin, A.; Minea, T.
2009-03-01
The properties of the secondary plasma created inside the ITER Neutral Beam Injector (NBI) neutraliser, through the interaction of the high energetic hydrogen beam with the molecular hydrogen gas, have been analysed. Starting from the results of our OBI-2 PIC Monte-Carlo numerical code, detailed kinetic of the hydrogen plasma has been studied using a Collisional-Radiative model. In this model, the electron distribution function is determined by solving a Boltzmann equation, whereas main plasma species are derived from balance equations. This paper presents preliminary results obtained in a 0D geometry, boundary conditions bing introduced through effective rates for gain and loss of particles at the neutraliser walls. It has been found that the main ion specie is H2+, essentially coming from the ionisation of the target gas. The electron energy distribution function is not maxwellian and its mean energy is about 5 eV. The plasma-wall interactions yield a strong contribution, in particular regarding the density of molecular ion H3+. Assuming several independent slices of plasma along the negative ions beam axis, the axial profile of the secondary plasma has been analysed. It has been found that the density and mean energy profiles of the plasma electrons are directly related to the plasma potential profile, which in turn closely follows the gas density one.
Spectral evolution of two-dimensional kinetic plasma turbulence in the wavenumber-frequency domain
Comişel, H.; Verscharen, D.; Narita, Y.; Motschmann, U.
2013-09-15
We present a method for studying the evolution of plasma turbulence by tracking dispersion relations in the energy spectrum in the wavenumber-frequency domain. We apply hybrid plasma simulations in a simplified two-dimensional geometry to demonstrate our method and its applicability to plasma turbulence in the ion kinetic regime. We identify four dispersion relations: ion-Bernstein waves, oblique whistler waves, oblique Alfvén/ion-cyclotron waves, and a zero-frequency mode. The energy partition and frequency broadening are evaluated for these modes. The method allows us to determine the evolution of decaying plasma turbulence in our restricted geometry and shows that it cascades along the dispersion relations during the early phase with an increasing broadening around the dispersion relations.
Plasma surface functionalization and dyeing kinetics of Pan-Pmma copolymers
NASA Astrophysics Data System (ADS)
Labay, C.; Canal, C.; Rodríguez, C.; Caballero, G.; Canal, J. M.
2013-10-01
Fiber surface modification with air corona plasma has been studied through dyeing kinetics under isothermal conditions at 30 °C on an acrylic-fiber fabric with a cationic dye (CI Basic Blue 3) analyzing the absorption, desorption and fixing on the surface of molecules having defined cationic character. The initial dyeing rate in the first 60 s indicates an increase of 58.3% in the dyeing rate due to the effect of corona plasma on the acrylic fiber surface. At the end of the dyeing process, the plasma-treated fabrics absorb 24.7% more dye, and the K/S value of the acrylic fabric increases by 8.8%. With selected dyestuff molecules, new techniques can be designed to amplify the knowledge about plasma-treated surface modifications of macromolecules.
A coarse-grained kinetic equation for neutral particles in turbulent fusion plasmas
Mekkaoui, A.; Marandet, Y.; Genesio, P.; Rosato, J.; Stamm, R.; Capes, H.; Koubiti, M.; Godbert-Mouret, L.; Catoire, F.
2012-06-15
A coarse-grained kinetic equation for neutral particles (atoms, molecules) in magnetized fusion plasmas, valid on time scales large compared to the turbulence correlation time, is presented. This equation includes the effects of plasma density fluctuations, described by gamma statistics, on the transport of neutral particles. These effects have so far been neglected in plasma edge modeling, in spite of the fact that the amplitude of fluctuations can be of order unity. Density fluctuations are shown to have a marked effect on the screening of neutrals and on the spatial localization of the ionization source, in particular at high density. The coarse-grained equations obtained in this work are readily implemented in edge code suites currently used for fusion plasma analysis and future divertor design (ITER, DEMO).
Dimension reduction of non-equilibrium plasma kinetic models using principal component analysis
NASA Astrophysics Data System (ADS)
Peerenboom, Kim; Parente, Alessandro; Kozák, Tomáš; Bogaerts, Annemie; Degrez, Gérard
2015-04-01
The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.
High-Performance Kinetic Plasma Simulations with GPUs and load balancing
NASA Astrophysics Data System (ADS)
Germaschewski, Kai; Ahmadi, Narges; Abbott, Stephen; Lin, Liwei; Wang, Liang; Bhattacharjee, Amitava; Fox, Will
2014-10-01
We will describe the Plasma Simulation Code (PSC), a modern particle-in-cell code with GPU support and dynamic load balancing capabilities. For 2-d problems, we achieve a speed-up of up to 6 × on the Cray XK7 ``Titan'' using its GPUs over the well-known VPIC code, which has been optimized for conventional CPUs with SIMD support. Our load-balancing algorithm employs a space-filling Hilbert-Peano curve to maintain locality and has shown to keep the load balanced within approximately 10% in production runs which otherwise slow down up to 5 × with only static load balancing. PSC is based on the
Kinetic plasma processes occurring in the outer plasmasphere
NASA Technical Reports Server (NTRS)
Wilson, Gordon R.
1992-01-01
One area of data analysis work that was begun under this contract is the fitting of the perpendicular velocity distributions of equatorially trapped ions with a Kappa function. This type of characterization of the trapped ions will be very useful for comparison with velocity distributions produced by the model. A second area of data analysis is to study data from consecutive passes when DE 1's apogee was near the magnetic equator and the spacecraft was often skimming along nearly the same L shell. In 1982 three such periods occurred in May, June, and July. For these consecutive events we have Kp histories, density measurements from a number of sources (Whistler data, DE SFR, ISEE SFR) and consecutive samples of ion pitch angle distributions along field lines. It is clear from this data how the pitch angle distributions evolve during a flux tube refilling event. Our modeling of the flow of plasma along closed field lines is following two basic tracks. The first is a study of the basic refilling process without the effect of wave-particle heating near the equator or the effect of large or abrupt field-aligned electric potential drops. This model includes the effects of Coulomb self-collisions and collisions with the O+ ions in the topside ionosphere. The second track is a study of the effects of wave produced pitch-angle scattering and perpendicular heating occurring near the magnetic equator, in connection with the development of large potential drops that result from electron heating and the development of density gradients.
Kinetic modeling of E-to-H mode transition in inductively coupled hydrogen plasmas
NASA Astrophysics Data System (ADS)
Nishida, K.; Mattei, S.; Mochizuki, S.; Lettry, J.; Hatayama, A.
2016-06-01
Radio Frequency (RF) Inductively Coupled Plasmas (ICPs) are widely known for their two discharge modes, i.e., H-mode and E-mode, where the dynamics of the plasmas are completely different from each other. We have performed a kinetic simulation of a hydrogen plasma discharge in order to clarify the discharge mechanism and the E-to-H transition of the RF ICPs. The numerical simulation results, such as the time variations of spatial distribution of electron density and the power dissipated in the plasma, show the characteristic changes of the plasma dynamics due to E-to-H mode transition. Especially, the drastic change during the mode transition has been observed in the time evolution of the electron energy distribution function (EEDF). The EEDF deviates from a Maxwellian distribution before/after the transition and the deviation is more significant in the E-mode phase. These results indicate the importance of kinetic modeling for the physical understanding of E-to-H transition.
Kinetic Alfvén solitary and rogue waves in superthermal plasmas
Bains, A. S.; Li, Bo Xia, Li-Dong
2014-03-15
We investigate the small but finite amplitude solitary Kinetic Alfvén waves (KAWs) in low β plasmas with superthermal electrons modeled by a kappa-type distribution. A nonlinear Korteweg-de Vries (KdV) equation describing the evolution of KAWs is derived by using the standard reductive perturbation method. Examining the dependence of the nonlinear and dispersion coefficients of the KdV equation on the superthermal parameter κ, plasma β, and obliqueness of propagation, we show that these parameters may change substantially the shape and size of solitary KAW pulses. Only sub-Alfvénic, compressive solitons are supported. We then extend the study to examine kinetic Alfvén rogue waves by deriving a nonlinear Schrödinger equation from the KdV equation. Rational solutions that form rogue wave envelopes are obtained. We examine how the behavior of rogue waves depends on the plasma parameters in question, finding that the rogue envelopes are lowered with increasing electron superthermality whereas the opposite is true when the plasma β increases. The findings of this study may find applications to low β plasmas in astrophysical environments where particles are superthermally distributed.
First steps towards the reaction kinetics of HMDSO in an atmospheric pressure plasma jet in argon
NASA Astrophysics Data System (ADS)
Loffhagen, Detlef; Becker, Markus M.; Foest, Rüdiger; Schäfer, Jan; Sigeneger, Florian
2014-10-01
Hexamethyldisiloxane (HMDSO) is a silicon-organic compound which is often used as precursor for thin-film deposition by means of plasma polymerization because of its high deposition rate and low toxicity. To improve the physical understanding of the deposition processes, fundamental investigations have been performed to clarify the plasma-chemical reaction pathways of HMDSO and their effect on the composition and structure of the deposited film. The current contribution represents the main primary and secondary plasma-chemical processes and their reaction products in the effluent region of an argon plasma jet at atmospheric pressure. The importance of the different collision processes of electrons and heavy particles are discussed. Results of numerical modelling of the plasma jet and the Ar-HMDSO reaction kinetics indicate that the fragmentation of HMDSO is mainly initiated by collisions with molecular argon ions, while Penning ionization processes play a minor role for the reaction kinetics in the effluent region of the jet. The work has been supported by the German Research Foundation (DFG) under Grant LO 623/3-1.
Weakly relativistic quantum kinetic theory for electrostatic wave modes in magnetized plasmas
Hussain, Azhar; Stefan, Martin; Brodin, Gert
2014-03-15
We have derived the electrostatic dispersion relation in a magnetized plasma using a recently developed quantum kinetic model based on the Dirac equation. The model contains weakly relativistic spin effects such as Thomas precession, the polarization currents associated with the spin and the spin-orbit coupling. It turns out that for strictly electrostatic perturbations the non-relativistic spin effects vanish, and the modification of the classical dispersion relation is solely associated with the relativistic terms. Several new wave modes appear due the electron spin effects, and an example for astrophysical plasmas are given.
Guan, G.; Mauel, M.E. ); Holber, W.M.; Caughman, J.B.O. )
1992-12-01
A fluid description of the presheath of a magnetized plasma is used to model a divergent electron cyclotron resonance (ECR) plasma source. The fluid equations are moments of the time-independent Boltzman equation when cross-field particle motion occurs only through a static {bold E}{times}{bold B} drift. Closure is obtained by neglecting third-order moments. The electrons are assumed to have constant temperature along the magnetic field, to obey a Maxwell--Boltzmann potential-density relationship, and to be warmer than the ions. Interactions between plasma and neutral gas are included by specifying the profile of the gas density along the magnetic field and collision cross sections. A form of the equations is derived that can be used to study ions with anisotropic temperatures. The model is used to estimate the axial profiles of the density, ion flow, and ion temperature of an ECR plasma source. The calculated global relationships between (1) the electron temperature and the equilibrium neutral gas density, and (2) the absorbed microwave power and the equilibrium plasma density are comparable with experimental measurements. Furthermore, the calculated ion temperature is comparable to recently reported measurements (Appl. Phys. Lett. {bold 57}, 661 (1990) and Appl. Phys. Lett. {bold 58}, 458 (1991)).
Kinetic theory of low-frequency cross-field instability in a weakly ionized plasma. II
Dimant, Y.S.; Sudan, R.N.
1995-04-01
The consistent kinetic approach developed in Paper I [Ya. S. Dimant and R. N. Sudan, Phys. Plasmas {bold 2}, 1157 (1995)] is applied to obtain the general dispersion relation of the two-stream {bold E}{times}{bold B} instability in collisionally dominated weakly ionized plasmas for wave frequencies small compared to the ion--neutral collision frequency. This dispersion relation covers the whole low-frequency band from the asymptotic short-wave limit studied in Paper I to the long-wave limit. Previous theories employing simplified kinetic theory or fluid equations for electron behavior are only correct in the long-wave limit. The principal new results are that the threshold conditions for this instability and the growth rates are altered from those predicted by earlier simplified theories. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.
NASA Astrophysics Data System (ADS)
Yin, L.; Albright, B. J.; Bergen, B.; Bowers, K. J.; Vold, E. L.; Molvig, K.; Fernández, J. C.; Bang, W.; Bradley, P. A.; Gautier, D. C.; Hamilton, C. E.; Palaniyappan, S.; Santiago Cordoba, M. A.; Hegelich, B. M.; Dyer, G.; Roycroft, R.
2015-11-01
Mixing of high-Z/low-Z interfaces in dense plasma media is a problem of importance for understanding mix in inertial confinement fusion experiments and recent experiments at the LANL Trident facility. In this presentation, we apply the VPIC particle-in-cell code with a binary collision model to explore kinetic effects of the atomic mixing. Comparisons are made to published analytic theory and hybrid modeling results and conditions are identified under which plasma kinetic behavior may lead to anomalously rapid atomic mixing. Work performed under the auspices of the U.S. DOE by the LANS, LLC, Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396. Funding provided by the Los Alamos National Laboratory Directed Research and Development Program.
AEGIS-K code for linear kinetic analysis of toroidally axisymmetric plasma stability
NASA Astrophysics Data System (ADS)
Zheng, L. J.; Kotschenreuther, M. T.; Van Dam, J. W.
2010-05-01
A linear kinetic stability code for tokamak plasmas: AEGIS-K (Adaptive EiGenfunction Independent Solutions-Kinetic), is described. The AEGIS-K code is based on the newly developed gyrokinetic theory [L.J. Zheng, M.T. Kotschenreuther, J.W. Van Dam, Phys. Plasmas 14 (2007) 072505]. The success in recovering the ideal magnetohydrodynamics (MHD) from this newly developed gyrokinetic theory in the proper limit leads the AEGIS-K code to be featured by being fully kinetic in essence but hybrid in appearance. The radial adaptive shooting scheme based on the method of the independent solution decomposition in the MHD AEGIS code [L.J. Zheng, M.T. Kotschenreuther, J. Comp. Phys. 211 (2006) 748] is extended to the kinetic calculation. A numerical method is developed to solve the gyrokinetic equation of lowest order for the response to the independent solutions of the electromagnetic perturbations, with the quasineutrality condition taken into account. A transform method is implemented to allow the pre-computed Z-function (i.e., the plasma dispersion function) to be used to reduce the integration dimension in the moment calculation and to assure the numerical accuracy in determining the wave-particle resonance effects. Periodic boundary condition along the whole banana orbit is introduced to treat the trapped particles, in contrast to the usual reflection symmetry conditions at the banana tips. Due to the adaptive feature, the AEGIS-K code is able to resolve the coupling between the kinetic resonances and the shear Alfvén continuum damping. Application of the AEGIS-K code to compute the resistive wall modes in ITER is discussed.
Fokker–Planck kinetic modeling of suprathermal α-particles in a fusion plasma
Peigney, B.E.
2014-12-01
We present an ion kinetic model describing the transport of suprathermal α-particles in inertial fusion targets. The analysis of the underlying physical model enables us to develop efficient numerical methods to simulate the creation, transport and collisional relaxation of fusion reaction products (α-particles) at a kinetic level. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A two-energy-scale approach leads to a self-consistent modeling of the coupling between suprathermal α-particles and the thermal bulk of the imploding plasma. This method provides an accurate numerical treatment of energy deposition and transport processes involving suprathermal particles. The numerical tools presented here are then validated against known analytical results. This enables us to investigate the potential role of ion kinetic effects on the physics of ignition and thermonuclear burn in inertial confinement fusion schemes.
Geodesic acoustic mode in anisotropic plasmas using double adiabatic model and gyro-kinetic equation
Ren, Haijun; Cao, Jintao
2014-12-15
Geodesic acoustic mode in anisotropic tokamak plasmas is theoretically analyzed by using double adiabatic model and gyro-kinetic equation. The bi-Maxwellian distribution function for guiding-center ions is assumed to obtain a self-consistent form, yielding pressures satisfying the magnetohydrodynamic (MHD) anisotropic equilibrium condition. The double adiabatic model gives the dispersion relation of geodesic acoustic mode (GAM), which agrees well with the one derived from gyro-kinetic equation. The GAM frequency increases with the ratio of pressures, p{sub ⊥}/p{sub ∥}, and the Landau damping rate is dramatically decreased by p{sub ⊥}/p{sub ∥}. MHD result shows a low-frequency zonal flow existing for all p{sub ⊥}/p{sub ∥}, while according to the kinetic dispersion relation, no low-frequency branch exists for p{sub ⊥}/p{sub ∥}≳ 2.
Plummer, Hillary A; Oliver, Gretchen D
2016-08-01
In order to decrease the amount of time that it takes the catcher to throw the ball, a catcher may chose to throw from the knees. Upper extremity kinematics may play a significant role in the kinetics about the elbow observed in catchers throwing from the knees. If relationships between kinematics and kinetics exist then the development of training and coaching instruction may help in reduced upper extremity injury risk. Twenty-two baseball and softball catchers (14.36±3.86years; 165.11±17.54cm; 65.67±20.60kg) volunteered. The catchers exhibited a less trunk rotation (5.6±16.2°), greater elbow flexion (87.9±21.4°) and decreased humeral elevation (71.1±12.3°) at the event of maximum shoulder external rotation as compared to what has previously reported in catchers. These variables are important, as they have previously been established as potential injury risk factors in pitchers, however it is not yet clear the role these variables play in catchers' risk of injury. A positive relationship between elbow varus torque during the deceleration phase and elbow flexion at MIR was observed (r=0.609; p=0.003). Throwing from the knees reduces a catcher's ability to utilize the proximal kinetic chain and this may help to explain why their kinematics and kinetics differ from what has previously been presented in the literature. PMID:26360828
A new dynamic fluid-kinetic model for plasma transport within the plasmaspheric plume
NASA Astrophysics Data System (ADS)
Wang, Y.; Tu, J.; Song, P.
2011-12-01
A new dynamic fluid-kinetic (DyFk) model is proposed and developed for investigating the plasma transport from the plasmasphere to the dayside magnetopause through the plasmaspheric plume. This model treats a closed flux tube in a local sense, in contrast to the global sense. The flux tube is allowed to move both radially from near the Earth to the magnetopause, which may result in expansion in its volume, and azimuthally around the Earth. Plasma may flow along the flux tube. The numerical simulation model couples a truncated version of the field line interhemispheric plasma (FLIP) model at altitudes below 800 km and a generalized semi-kinetic (GSK) model above it with an overlapped boundary region in each of the hemispheres. A self-consistently treatment of the ionospheric losses and production with possible heat sinks couples to a kinetic treatment of the multiple ion species (O+/ H+/ He+) and electrons in the plasmasphere. This model includes the effects of the convection of the plasmaspheric flux tube, parallel electric field, magnetic mirror force, centrifugal force, changing ionospheric conditions, Coulomb and ion-neutral collisions, and anisotropic temperatures, as well as the wave-particle interaction. The preliminary simulation results of the multi-species ion transport within a plasmaspheric plume will be presented.
Kolobov, Vladimir; Arslanbekov, Robert; Frolova, Anna
2014-12-09
The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.
NASA Astrophysics Data System (ADS)
Kolobov, Vladimir; Arslanbekov, Robert; Frolova, Anna
2014-12-01
The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.
Inertial-range kinetic turbulence in pressure-anisotropic astrophysical plasmas
NASA Astrophysics Data System (ADS)
Kunz, M. W.; Schekochihin, A. A.; Chen, C. H. K.; Abel, I. G.; Cowley, S. C.
2015-10-01
> A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al., Astrophys. J. Suppl. Ser., vol. 182, 2009, pp. 310-377) to the case where the mean distribution function of the plasma is pressure-anisotropic and different ion species are allowed to drift with respect to each other - a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas such as the intracluster medium. Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g. the Alfvén ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. Beyond these order-unity corrections, the main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvénic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvénic cascade is fluid, satisfying RMHD equations (with the Alfvén speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping (and for a bi-Maxwellian plasma splits into three independent collisionless cascades). Secondly, the organising principle of this turbulence is elucidated in the form of a conservation law for the appropriately generalised kinetic free energy. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses, and that these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics
Plasma arginine and leucine kinetics and urea production rates in burn patients.
Yu, Y M; Young, V R; Castillo, L; Chapman, T E; Tompkins, R G; Ryan, C M; Burke, J F
1995-05-01
We measured plasma arginine and leucine kinetics and rates of urea production (appearance) in 12 severely burned patients (mean body surface burn area, 48%) during a basal state (low-dose intravenous glucose) and while receiving routine, total parenteral nutrition ([TPN] fed state) including an L-amino acid mixture, supplying a generous level of nitrogen (mean, 0.36 g N.kg-1.d-1). The two nutritional states were studied in random order using a primed 4-hour constant intravenous tracer infusion protocol. Stable-nuclide-labeled tracers were L-[guanidino-13C]arginine, L-[1-13C]leucine, [18O]urea, and NaH13CO3 (prime only), with blood and expired air samples drawn at intervals to determine isotopic abundance of arginine, citrulline, ornithine, alpha-ketoisocaproate ([KIC] for leucine), and urea in plasma and 13CO2 in breath. Results are compared with data obtained in these laboratories in healthy adults. Leucine kinetics (flux and disappearance into protein synthesis) indicated the expected higher turnover in burn patients than in healthy controls. Mean leucine oxidation rates are also higher and compared well with values predicted from urea production rates, provided that urea nitrogen recycling via intestinal hydrolysis is taken into account. The plasma urea flux was also higher than for normal subjects. Arginine fluxes as measured in the systemic whole body, via the plasma pool, were correspondingly higher in burned patients than in healthy controls and were in good agreement with values predicted from leucine-KIC kinetics. However, systemic whole-body arginine flux measured via the plasma pool was only 20% of the arginine flux estimated from the urea flux plus the rate of protein synthesis.(ABSTRACT TRUNCATED AT 250 WORDS) PMID:7752916