Search

Turbulent plasmas notably self-organize to higher energy states upon application of additional free energy sources or modification of edge operating conditions. Mechanisms whereby such bifurcations occur have been actively debated for decades. Enhanced confinement occurs at the plasma edge, where a shortfall of predicted turbulence intensity has been puzzling scientists for decades. We show, from the primitive kinetic equations that both problems are connected and that interplay of confined plasma turbulence with its material boundaries is essential to curing the shortfall of predicted turbulence and to triggering spontaneous transport barrier onset at the plasma edge. Both problems determine access to improved confinement and are central to fusion research. A comprehensive discussion of the underlying mechanisms is proposed. These results, highly relevant to the quest for magnetic fusion may also be generic to many problems in fluids and plasmas where turbulence self-advection is active.

We provide the 21N-moment multi-temperature collision coefficients for the Boltzmann collision operator using the Sonine-Hermite polynomial ansatz in the style of Zhdanov et al. First, we outline the general derivation method. Then, we provide the collision coefficients in the most general form in terms of the Chapman-Cowling integrals for any potential of interaction for the 21N-moment approximation. Then, we provide the the collision coefficients for the specific case of the approximated Coulomb potential cross-section with Debye cutoff. These coefficients help in order to find easy implementation in current SOL/edge fluid packages which currently implement the 21N-moment single-temperature coefficients. This provides the completion of a missing link in the current literature., Comment: 19 pages, 0 figures. Accepted by Contributions to Plasma Physics (Dec 2021). arXiv admin note: substantial text overlap with arXiv:2103.02455

The derivation of the multi-temperature generalized Zhdanov closure is provided starting from the most general form of the left hand side of the moment averaged kinetic equation with the Sonine-Hermite polynomial ansatz for an arbitrary number of moments. The process of arriving at the reduced higher-order moment equations, with its assumptions and approximations, is explicitly outlined. The generalized multi-species, multi-temperature coefficients from the authors' previous article are used to compute values of higher order moments such as heat flux in terms of the lower order moments. Transport coefficients and the friction and thermal forces for magnetic confinement fusion relevant cases with the generalized coefficients are compared to the scheme with the single-temperature coefficients previously provided by Zhdanov et al. It is found that the 21N-moment multi-temperature coefficients are adequate for most cases relevant to fusion. Furthermore, the 21N-moment scheme is also tested against the trace approximation to determine the range of validity of the trace approximation with respect to fusion relevant plasmas. Possible refinements to the closure scheme are illustrated as well, in order to account for quantities which might be significant in certain schemes such as the drift approximation., Comment: Submitted to Plasma Physics and Controlled Fusion, July 15, 2021

Grad's method is used on the linearized Boltzmann collision operator to derive the most general expressions for the collision coefficients for a multi-component, multi-temperature plasma up to rank-2. In doing so, the collision coefficients then get expressed as series sum of pure coefficients of temperature and mass ratios multiplied by the cross-section dependent Chapman-Cowling integrals. These collisional coefficients are compared to previously obtained coefficients by Zhdanov et al [Zhdanov V.M., Transport processes in multi-component plasma, Taylor and Francis (2002)] for 13N-moment multi-temperature scheme. First, the differences in coefficients are compared directly, and then the differences in first approximation to viscosity and friction force are compared. For the 13N-moment multi-temperature coefficients, it is found that they behave reasonably similarly for small temperature differences, but display substantial differences in the coefficients when the temperature differences are high, both for the coefficients and for viscosity and friction force values. Furthermore, the obtained coefficients are compared to the 21N-moment single-temperature approximation provided by Zhdanov et al, and it is seen that the differences are higher than the 13N-moment multi-temperature coefficients, and have substantial differences even in the vicinity of equal temperatures, especially for the viscosity and friction force calculations., Comment: Corrected some typos from the previous version. Accepted by Plasma Physics and Controlled Fusion (April 9,2021)

We present the Multi-Particle-Collision (MPC) dynamics approach to simulate properties of low-dimensional systems. In particular, we illustrate the method for a simple model: a one-dimensional gas of point particles interacting through stochastic collisions and admitting three conservation laws (density, momentum and energy). Motivated from problems in fusion plasma physics, we consider an energy-dependent collision rate that accounts for the lower collisionality of high-energy particles. We study two problems: (i) the collisional relaxation to equilibrium starting from an off-equilibrium state and (ii) the anomalous dynamical scaling of equilibrium time-dependent correlation functions. For problem (i), we demonstrate the existence of long-lived population of suprathermal particles that propagate ballistically over a quasi-thermalized background. For (ii) we compare simulations with the predictions of nonlinear fluctuating hydrodynamics for the structure factors of density fluctuations. Scaling analysis confirms the prediction that such model belong to the Kardar-Parisi-Zhang universality class., Comment: 12 pages, 4 figures. Accepted for publication in Proceedings of the IHP workshop "Numerical aspects of nonequilibrium dynamics", 25-27 April 2017 Paris, France, in "Proceedings in Mathematics and Statistics" (Springer)

In order to improve the presently used ad hoc flux limiter treatment of parallel heat flux transport in edge plasma fluid codes we consider here a generalized version of the Fourier law implementing a non-local kernel for the heat flux computation. The Bohm boundary condition at the wall is recovered introducing a volumetric loss term representing the contribution of suprathermal particles to the energy out flux. As expected, this contribution is negligible in the strongly collisional regime while it becomes more and more dominant for marginally and low collisional regimes. In the second part of the paper, we consider a kinetic approach where collisions are considered using the Multi-Particle-Collision (MPC) algorithm. Kinetic simulation results at medium and low collisionality are also reported., Comment: 13 Pages, 6 figures, Accepted for publication in Contributions to Plasma Physics

The beam-plasma instability, i.e. the response of the plasma bulk to the injection of supra thermal charged-particle beams, results to be appropriately characterized by a long-range interaction system. This physical system hosts a number of very interesting phenomena and, in particular, the emergence of long-lived quasi-stationary states. We characterize the self-consistent distribution functions of such out-of-equilibrium states by means of the Lynden-Bell's theory. The prediction of this theory, based on the statistical mechanics of the Vlasov equation, are checked against the outcomes of numerical simulations of the discrete system. Moreover, a phenomenological study of the effective resonance band for the system response is also addressed. A threshold value has been found in the initial spread of beam-particle momenta. This threshold allows to discriminate between the resonant and non-resonant regimes. The analysis of the thermalization of a few percents of the beam population characterized by large initial momenta (with respect to the main part of the beam itself) is also performed and it confirms and deepens the understanding of the physical meaning of the mentioned threshold., Comment: 21 pages, 14 figures

Inaccurate description of the equilibrium can yield to spurious effects in gyrokinetic turbulence simulations. Also, the Vlasov solver and time integration schemes impact the conservation of physical quantities, especially in long-term simulations. Equilibrium and Vlasov solver have to be tuned in order to preserve constant states (equilibrium) and to provide good conservation property along time (mass to begin with). Several illustrative simple test cases are given to show typical spurious effects that one can observes for poor settings. We explain why Forward Semi-Lagrangian scheme bring us some benefits. Some toroidal and cylindrical GYSELA runs are shown that use FSL., Comment: No. RR-8054 (2012)

We consider the radial transport of test particles due to the ExB drift motion in the guiding center approximation. Using an explicit expression to modify the electrostatic potential, we show that it is possible to construct a transport barrier which suppresses radial transport. We propose an algorithm for the implementation of this local modification computed from an electrostatic potential known on a spatio-temporal grid. The number of particles which escape the inner region defined by the barrier measures the efficiency of the control. We show that the control is robust by showing a significant reduction of radial transport, when applied with a reduced number of probes aligned on a circle., Comment: 13 pages, 9 figures, 1 column

A finite Larmor radius approximation is derived from the classical Vlasov equation, in the limit of large (and uniform) external magnetic field. We also provide an heuristic derivation of the electroneutrality equation in the finite Larmor radius setting. Existence and uniqueness of a solution is proven in the stationary frame for solutions depending only on the direction parallel to the magnetic field and factorizing in the velocity variables.

We consider a 1 1/2degrees of freedom Hamiltonian dynamical system, which models the chaotic dynamics of charged test-particles in a turbulent electric field, across the confining magnetic field in controlled thermonuclear fusion devices. The external electric field E = \nabla\bigvee is modeled by a phenomenological potential V and the magnetic field B is considered uniform. It is shown that, by introducing a small additive control term to the external electric field, it is possible to create a transport barrier for this dynamical system. The robustness of this control method is also investigated. This theoretical study indicates that alternative transport barriers can be triggered without requiring a control action on the device scale as in present Internal Transport Barriers (ITB)., Comment: 18

The ${\bm E}\times{\bm B}$ drift motion of charged test particle dynamics in the Scrape Off Layer (SOL)is analyzed to investigate a transport control strategy based on Hamiltonian dynamics. We model SOL turbulence using a 2D non-linear fluid code based on interchange instability which was found to exhibit intermittent dynamics of the particle flux. The effect of a small and appropriate modification of the turbulent electric potential is studied with respect to the chaotic diffusion of test particle dynamics. Over a significant range in the magnitude of the turbulent electrostatic field, a three-fold reduction of the test particle diffusion coefficient is achieved.

The gyrokinetics formulation of plasmas in strong magnetic fields aims at the elimination of the angle associated with the Larmor rotation of charged particles around the magnetic field lines. In a perturbative treatment or as a time-averaging procedure, gyrokinetics is in general an approximation to the true dynamics. Here we discuss the conditions under which gyrokinetics is either an approximation or an exact operation in the framework of reduction of dynamical systems with symmetry, Comment: 15 pages latex

We present a method of control which is able to create barriers to magnetic field line diffusion by a small modification of the magnetic perturbation. This method of control is based on a localized control of chaos in Hamiltonian systems. The aim is to modify the perturbation locally by a small control term which creates invariant tori acting as barriers to diffusion for Hamiltonian systems with two degrees of freedom. The location of the invariant torus is enforced in the vicinity of the chosen target. Given the importance of confinement in magnetic fusion devices, the method is applied to two examples with a loss of magnetic confinement. In the case of locked tearing modes, an invariant torus can be restored that aims at showing the current quench and therefore the generation of runaway electrons. In the second case, the method is applied to the control of stochastic boundaries allowing one to define a transport barrier within the stochastic boundary and therefore to monitor the volume of closed field lines.

We present a method to control transport in Hamiltonian systems. We provide an algorithm - based on a perturbation of the original Hamiltonian localized in phase space - to design small control terms that are able to create isolated barriers of transport without modifying other parts of phase space. We apply this method of localized control to a forced pendulum model and to a system describing the motion of charged particles in a model of turbulent electric field.

We present a technique to control chaos in Hamiltonian systems which are close to integrable. By adding a small and simple control term to the perturbation, the system becomes more regular than the original one. We apply this technique to a model that reproduces turbulent ExB drift and show numerically that the control is able to drastically reduce chaotic transport.

It is shown that a relevant control of Hamiltonian chaos is possible through suitable small perturbations whose form can be explicitly computed. In particular, it is possible to control (reduce) the chaotic diffusion in the phase space of a Hamiltonian system with 1.5 degrees of freedom which models the diffusion of charged test particles in a ``turbulent'' electric field across the confining magnetic field in controlled thermonuclear fusion devices. Though still far from practical applications, this result suggests that some strategy to control turbulent transport in magnetized plasmas, in particular tokamaks, is conceivable., Comment: Latex2e-Revtex, 4 pages, 6 figures