PPPL

Research Seminars

Upcoming

  • Staircase Safety Factor Profiles in Gyrokinetic Simulations at Low Magnetic Shear (abstract)
    Arnas Volcokas, EPFL
    #s1866, Thursday, 17 Oct 2024, 9:45am

    In this work, we investigate the impact of electromagnetic effects on plasma turbulence self-organization at low magnetic shear using nonlinear gyrokinetics. Our previous electrostatic studies showed that turbulent eddies extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear s is low or zero. Such ``ultra-long’’ eddies have significant consequences on turbulent transport due to parallel self-interaction. At low magnetic shear, parallel self-interaction induces strong corrugations in plasma profiles at low-order rational surfaces, including the formation of stationary current layers. When electromagnetic effects are considered, turbulence-generated currents lead to the development of stationary zonal magnetic potential, locally flattening the safety factor profile to form staircase structures or broaden the safety factor profile minimum. This represents a crucial feedback mechanism between turbulence and the imposed safety factor profile, resulting in a reduction in turbulent transport. We study the corrugated safety factor profiles using both the local flux tube code GENE and the global particle-in-cell code ORB5. To further explore this interaction, we employed a novel extension of the flux tube model, allowing simulations of non-uniform magnetic shear profiles, including minimum-q profiles relevant for Internal Transport Barrier (ITB) formation. Our findings indicate that turbulence-generated current layers can flatten the imposed non-uniformity across the entire domain or substantially widen rational surface regions, consistent with global simulation results. We believe these results are relevant for understanding ITB formation and inform long-standing experimental observations.

  • Two etudes on electromagnetic turbulence in idealised fusion plasmas (abstract)
    Alex Schekochihin, University of Oxford
    #s1867, Tuesday, 15 Oct 2024, 1:45pm

    Schemes involving higher-beta fusion devices have led to intensified interest in electromagnetic turbulence in fusion plasmas. I will discuss this topic in two radically simplified settings: ITG-driven turbulence in a Z-pinch and interchange-driven turbulence in open-field-line systems. Both can be described by minimal models involving just a few “fluid” fields: density and temperature plus velocity and magnetic perturbations perpendicular to the mean magnetic field. Despite a drastic reduction in complexity that many seasoned fusion practitioners might consider excessive and indeed deeply immoral, these simple paradigms yield a wealth of interesting behaviour, some of which seems to capture key physical signatures of electromagnetic effects reported by more realistic numerical explorations — and might perhaps inform future such explorations (failing that, it is interesting enough to pique the curiosity of some theoreticians). Increasing beta in ITG-driven turbulence is shown to move the threshold for the low-to-high-transport (Dimits) transition towards lower temperature gradients. This is because in the struggle between the Reynolds and diamagnetic stresses that determines this transition in an electrostatic plasma [1], the Maxwell stress brought in by the Alfvenic activity associated with higher beta opposes the Reynolds stress and tips the system towards higher-transport states [2]. In the interchange-driven MHD turbulence, increasing beta (which in all such systems is equivalent also to increasing the parallel size of the system, increasing the temperature gradient or the field-line curvature) also pushes the system from an electrostatic low-transport state to an electromagnetic (very)-high-transport one — this is perhaps the simplest conceivable system exhibiting the Dimits transition [3] (simpler even than [1]). It is shown however that when 2D interchange motions are forbidden, the low-beta electrostatic state can be replaced by an electromagnetic one with very high transport and some rather strange properties (while it is interesting theoretically, as well as perhaps in some analogous astrophysical systems, the relevance of this last result to fusion devices is left to post-seminar discussion). Overall, the takeaway message is that electromagnetic turbulence in fusion plasmas is a physics-rich subject that, while challenging practical modellers, should also excite theoreticians.

    [1] P. Ivanov et al. JPP 86, 855860502 (2020) & JPP 88, 905880506 (2022)

    [2] Y. Zhang et al., in preparation (2024)

    [3] W. Clarke et al., in preparation (2024)

Past

  • Resistive hose modes in Tokamaks and stability of sinusoidal flows in the viscous Hasegawa-Mima equation
    Alexandre P. Sainterme, University of Wisconsin-Madison, abstract
    [#s1865, 03 Oct 2024]

    A beam of runaway electrons (REs) within a cold, resistive background plasma can be unstable to resistive hose modes. Helical deflection of the current centroid of the beam away from its initial position is imperfectly canceled by eddy currents induced in the cold background plasma. The resulting macro-scale, over-stable, kink-like oscillations grow on a time scale determined by the resistivity of the cold plasma. Using a fluid model for the RE beam, we find that these high frequency modes are faster growing than resistive MHD instabilities when resistivity is large [Phys. Plasmas 31, 010701 (2024)]. Initial value calculations using the fluid RE model implemented in the extended MHD code NIMROD are compared to a 1D ODE eigenvalue code, and a semi-analytic dispersion relation for a uniform beam current equilibrium. Nonlinear simulations of the hose mode are presented, and the consequences for RE electron confinement in post-disruption tokamaks—where the background thermal plasma is cold and the resistivity is large—are discussed.

    The Hasegawa-Mima equation is a local 2-D model for drift-wave turbulence in a magnetized plasma possessing a density gradient. Steady zonal flows with sinusoidal variation are known to be unstable to Kelvin-Helmholtz instabilities when u₀ / V* > λ² / Lₙ², where u₀ is the amplitude of the flow, V* is the electron diamagnetic drift velocity, λ is the wavelength of the flow, and Lₙ is the density gradient scale length [Zhu, et. al, Phys. Plasmas 25, 082121 (2018)]. With the addition of viscosity, we find that the parameter range for linear instability is modified, and that an additional unstable mode exists at very low kᵧ. The growth rate of the instability with viscosity is shown to go as kᵧ² at low kᵧ' and the frequency corresponds with the drift wave frequency ω*.

  • Magnetorotational Turbulence and Dynamo in Accretion Disks: a Unified Mean-Field Theory
    Tushar Mondal, International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, India, abstract
    [#s1850, 26 Sep 2024]

    Accretion is a fundamental cosmic process involved in phenomena ranging from planet formation to X-ray binaries and the generation of jets from protostars and black holes. Evidence suggests magnetorotational instability (MRI) drives the turbulence that transports angular momentum outwards, allowing gas to spiral inwards in adisk. MRI-induced turbulence depends on coherent magnetic fields, which are susceptible to dissipation unless maintained by a dynamo. The fundamental processes sustaining MRI turbulence and dynamos remain elusive due to the historical separation of mean-field dynamo and angular momentum transport problems.

    We investigate MRI turbulence and dynamo phenomena using direct statistical simulations in a zero net-flux, unstratified shearing box. Our approach begins with the development of a unified mean-field model that integrates the traditionally decoupled problems of large-scale dynamo action and angular momentum transport in accretion disks. The model includes a hierarchical set of equations, capturing up to second-order correlators, with a statistical closure approximation used for third-order correlators. We emphasize the complex web of interactions connecting various components of the stress tensors—Maxwell, Reynolds, and Faraday—through shear, rotation, correlators associated with mean fields, and nonlinear terms. By identifying the dominant interactions, we pinpoint the key mechanisms essential for the generation and maintenance of MRI turbulence. Our general mean-field model for MRI turbulence allows for a self-consistent formulation of the electromotive force, accounting for both inhomogeneities and anisotropies. In the context of large-scale magnetic field generation, we identify two critical mechanisms: the rotation-shear-current effect, which generates radial magnetic fields, and the rotation-shear-vorticity effect, responsible for vertical magnetic fields. We provide explicit, non-perturbative expressions for the transport coefficients associated with these dynamo effects. Importantly, both mechanisms rely on the inherent presence of large-scale vorticity dynamo within MRI turbulence. Lastly, I will discuss the associated magnetic helicity in such unstratified turbulent disks.

  • On the impact of temperature gradient flattening and system size on heat transport in microtearing turbulence (video)
    Ajay Chandrarajan Jayalekshmi, Dutch Institute for Fundamental Energy Research (DIFFER), abstract
    [#s1849, 23 Sep 2024]

    Microtearing instability is one of the major sources of turbulent transport in high-β tokamaks. These modes lead to very localized transport at low-order rational magnetic field lines, and we show that flattening of the local electron temperature gradient at these rational surfaces plays an important role in setting the saturated flux level in microtearing turbulence. This process depends crucially on the density of rational surfaces, and thus the system-size, and gives rise to a worse-than-gyro-Bohm transport scaling for system-sizes typical of existing tokamaks and simulations.

  • Entropy analysis for fusion kinetic plasmas: turbulence and energetic particle instability (video)
    Jungpyo Lee, Hanyang University, Korea, abstract, slides
    [#s1840, 19 Sep 2024]

    Entropy is a good fluid quantity for many quasi-equilibrium systems of physics, chemistry, and information. It has a good property of representing a thermodynamic state, a heat exchange, and a preferential direction of the irreversibility. Although the entropy concept can be useful to analyze the equilibrium (e.g. magnetic equilibrium), it has been hardly understood to be meaningful for the non-equilibrium, opened, and highly nonlinear system like fusion kinetic plasmas. In this presentation, we find the possibility of using entropy concept to interpret the nonlinear unsteady systems by two important examples of the fusion theory: (1) the zonal flow saturation for trapped electron mode (TEM) turbulence, and (2) the frequency chirping of Berk-Breizman (BB) model for energetic particle driven instability.

    Thanks to the entropy property describing the energy exchange by ohmic heating and a triad transfer in the wave-particle interactions, the entropy balance analysis is useful to find the energy flow of a drift-wave turbulence between the driving component by external constraints (temperature and density gradients), the electromagnetic waves, the plasma kinetic free energy, and the collisional dissipation, which giving the free-energy cascade in terms of the spatial wavevector. In this study, we examine the contribution of the zonal flow on the TEM saturation quantitatively by implementing the entropy balance diagnostics in CGYRO simulations. For a TEM case of a high ratio of the electron density gradient to the electron temperature gradients, the effective zonal flow shearing rate (or the advection rate by the zonal flow) is much lower than the growth rate, so the perpendicular diffusion is important as much as the zonal flow saturation, as shown in the entropy transfer evaluation. It could be useful to quantify the role of the zonal flow and the perpendicular diffusion, and improve the existing saturation model of a reduced model (e.g. TGLF) for the TEM mode.

    The BB model is widely used for interpreting the experimental observations of a periodic or a frequency chirping phenomena due to the non-disruptive MHD instability by energetic particles. If the bump-on tail instability is saturated by the external sink with a small collision, the free energy balance equation can be obtained by defining a new effective temperature and using the quasilinear part of the entropy production. The modified free energy balance is still useful to understand the nonlinear phenomena. In the chirping case, we found that the simulation results of the BB model regarding the chirping frequency and the kinetic distribution shape agree with the analysis maximizing the entropy of the background distribution function with the constraints of the external dissipation and Ampere’s law. It implies that the equipartition by entropy increase still sufficiently holds for the nonlinear unsteady process with a small Kubo number.

  • MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas (video)
    Francisco Javier Escoto, CIEMAT (Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas), Spain, abstract
    [#s1838, 14 Aug 2024]

    Stellarator magnetic fields must be optimized to achieve the confinement quality required for fusion reactors. Specifically, in order to exhibit radial neoclassical transport at low collisionality as small as tokamaks, stellarators need to be approximately omnigenous [1]. A magnetic field is omnigenous if, for all particles, the radial component of the drift velocity averages out to zero over the lowest-order orbits. Although radial neoclassical transport has typically been minimized indirectly by means of figures of merit such as the effective ripple [2], recent developments have enabled to carry out this minimization using accurate calculations of radial neoclassical transport at low collisionality [3]. However, not only neoclassical transport across flux-surfaces is important for stellarator optimization. In general, parallel transport produces a net electric current, known as bootstrap current, that can modify the magnetic field and, therefore, needs to be evaluated during the optimization process. This is particularly important if the goal is the design of approximately quasi-isodynamic fields (quasi-isodynamicity is the concept in which the magnetic configuration of Wendelstein 7-X is based on), a subclass of approximately omnigenous fields that have the additional property of giving small bootstrap current [4] and are compatible with island divertors. However, precise calculations of the bootstrap current were too slow to be included in a stellarator optimization loop so far.

    In this seminar I will present MONKES (MONoenergetic Kinetic Equation Solver) [5], a new neoclassical code that solves the same monoenergetic drift-kinetic equation as DKES [6], the workhorse of neoclassical transport calculations in stellarators. MONKES was conceived, among other things, to satisfy the need for fast and accurate calculations of the bootstrap current. By exploiting the tridiagonal structure of the drift-kinetic equation in a Legendre basis, it is possible to obtain accurate results for all monoenergetic transport coefficients (i.e. those giving radial as well as parallel transport) at low collisionality using a single core in approximately one minute. These features make MONKES ideal for its inclusion in stellarator optimization suites for direct optimization of the bootstrap current. Apart from optimization, MONKES can be used for the analysis of experimental discharges and be integrated into predictive transport frameworks.

    References

    [1] J. R. Cary and S. G. Shasharina. Phys. Plasmas 4, 3323 (1997).

    [2] V. V. Nemov et al. Physics of Plasmas 6, 4622 (1999).

    [3] J. L. Velasco et al. Journal of Computational Physics 418, 109512 (2020).

    [4] P. Helander and J. Nührenberg. Plasma Physics and Controlled Fusion 51, 055004 (2009).

    [5] F. J. Escoto et al. Nuclear Fusion 64, 076030 (2024).

    [6] S. P. Hirshman et al. Physics of Fluids 29, 2951 (1986).

  • Gyrokinetic modeling of tungsten transport and radiation in WEST and ASDEX Upgrade plasmas: towards the kinetic modeling of ITER physics (video)
    Julien Dominski, PPPL, abstract
    [#s1823, 20 Jun 2024]

    Understanding the impact of tungsten on tokamak performance and its mitigation by low-Z impurities is key to the success of ITER operations. A novel self-consistent gyrokinetic model, integrated into the XGC code [1,2], has been developed to comprehensively analyze tungsten transport and radiation. In this model, the tungsten ions are represented by a few bundles, with their fractional abundance determined by the atomic balance between ionization and recombination processes derived from ADAS rates [3]. The electron cooling by tungsten radiation is also derived from ADAS rates. This new model enables the study of tungsten radiation, transport dynamics, and interactions with low-Z species. Two complementary tungsten studies will be presented. First, we study the impact of nitrogen on the collisional and turbulent peaking factors of tungsten in a WEST plasma.Analysis of tungsten radiation with synthetic diagnostics will be discussed. Second, we explore the intricate interplay between collisional and turbulent transport in an H-mode plasma scenario of ASDEX Upgrade.

    [1] Dominski et al. J. Plasm. Phys. (2019)

    [2] Dominski et al. Phys Plasmas (2024)

    [3] https://www.adas.ac.uk

  • Integrating Recent Innovations in Single-Stage Stellarator Optimization for the Columbia Stellarator eXperiment
    Antoine Baillod, Columbia University, abstract
    [#s1818, 13 Jun 2024]

    The Columbia Stellarator eXperiment (CSX), currently in the design phase at Columbia University, is focused on investigating quasi-axisymmetric plasma with a small aspect ratio, and on validating recent developments in stellarator technology, theory, and optimization techniques. It is designed to test some of the theoretical predictions of quasi-axisymmetric plasmas, in particular plasma flow damping, MHD stability properties, and the study of trapped particle confinement. The magnetic field is generated by a set of two circular and planar poloidal field coils (PF coils) alongside two shaped interlinked coils (IL coils), with the potential consideration of additional coils to enhance shaping or experimental flexibility. The PF coils and vacuum vessel are repurposed from the former Columbia Non-Neutral Torus (CNT) experiment [1]. The two IL coils will be wound "in-house" at Columbia University, using non-insulated High-Temperature Superconducting (HTS) tapes. These coils undergo shape and strain optimization to produce the desired plasma configuration while adhering to numerous engineering constraints. Discovering a plasma shape that aligns with the physics objectives and can be produced by such a restricted number of coils poses a significant challenge. Indeed, the constrained coil set’s limited capacity to produce varied plasma shapes hinders the application of the traditional two-stage stellarator optimization approach. Instead, novel single-stage optimization techniques are employed, where plasma and coils are optimized concurrently. Despite an increased problem complexity due to the larger number of degrees of freedom, these methods find optimized plasma shapes that can be generated by coils that satisfy engineering constraints. We discuss two single-stage optimization methodologies[2, 3, 4]. We explore their application to the CSX experiment’s design, aiming to identify configurations that fulfill engineering constraints and generate a plasma within a desired regime for the experiment’s physics objectives.

    [1] Pedersen, T. S. et. al. (2006). Construction and Initial Operation of the Columbia Nonneutral Torus. Fusion Science and Technology 50 (3), 372-381

    [2] Jorge, R. et. al. (2023). Single-stage stellarator optimization: combining coils with fixed boundary equi libria. Plasma Physics and Controlled Fusion 67 (7), 074003

    [3] Giuliani, A. et. al. (2022). Direct computation of magnetic surfaces in Boozer coordinates and coil opti mization for quasisymmetry. Journal of Plasma Physics 88 (4), 905880401

    [4] Giuliani, A. et. al. (2023). Direct stellarator coil design using global optimization: application to a compre hensive exploration of quasi-axisymmetric devices. arXiv:2310.19097

  • New kind of ballooning modes with non-separable (in space and time) amplitudes: relevance to astrophysics and non-thermal fusion (video)
    Bruno Coppi, MIT, abstract
    [#s1789, 30 May 2024]

    New tridimensional plasma structures, that are oscillatory and classified as non-separable ballooning modes, can emerge in inhomogeneous plasmas and undergo resonant mode-particle interactions, e.g., with a minority population, that can lead them to modify their spatial profiles. Thus, unlike the case of previously known ballooning modes their amplitudes are not separable functions of space and time. The relevant resonance conditions are intrinsically different from those of the well-known Landau conditions for (ordinary) plasma waves: they involve the mode geometry and affect different regions of the distribution in momentum space at different positions in configuration space. The novel resonant mode-particle interactions constitute a direct (linear) process to exchange energy between different populations without the inefficiencies of nonlinear coupling processes. The new ballooning modes are relevant to circumbinary disks associated with pairs of black holes and to fusion burning plasmas that include an initially thermal population of fusion reacting nuclei and a population of high energy nuclei (reaction products). It is reasonable to expect that the distributions of the reacting nuclei in momentum space will not remain strictly Maxwellian and that the resulting reaction rates will be different from those evaluated for (conventional) thermalized plasmas.

  • Gyrokinetic Exact Fokker-Planck-Landau Collisions in Fusion Plasma Turbulence; and ITER-Relevant Turbulence Broadening of the Divertor Heat Flux Profile in DIII-D QH-Mode Pedestals (video)
    Darin R. Ernst, MIT, abstract, slides
    [#s1788, 16 May 2024]

    In this talk we summarize some of our recent work of interest to PPPL researchers, in two parts:

    Part I: First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models. Previous gyrokinetic simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This work demonstrates significant corrections using the first formulation [1, 2] and implementation [3, 4] of the gyrokinetic Fokker-Planck-Landau collision operator with the exact linearized field-particle terms. Realistic nonlinear gyrokinetic simulations of fusion plasma turbulence show significant corrections relative to the Sugama model collision operator for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes (the exact operator is now released in GENE-3.0). Future work will extend novel spectral methods implemented for the drift-kinetic operator [5] to the gyrokinetic operator.

    Part II: Broadening of the Divertor Heat Flux Profile in DIII-D QH-Modes, Matched by XGC. Multi-machine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. In the experiments presented [6], the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of ELMs or low-n MHD modes. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of (TEM) turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe and infra-red imaging measurements.

    [1] B. Li and D. R. Ernst, Phys. Rev. Lett. 106(19) 195002 (2011). https://doi.org/10.1103/PhysRevLett.106.195002

    [2] Q. Pan and D. R. Ernst, Phys. Rev. E 99, 023201 (2019). https://doi.org/10.1103/PhysRevE.99.023201

    [3] Q. Pan, D. R. Ernst, and D. Hatch, Phys. Rev. E Lett. 103, L051202 (2021). https://doi.org/10.1103/PhysRevE.103.L051202

    [4] Q. Pan, D. R. Ernst and P. Crandall, Physics of Plasmas 27, 042307 (2020). https://doi.org/10.1063/1.5143374

    [5] M. Landreman and D. R. Ernst, J. Comput. Phys. 243, 130 (2013). https://doi.org/10.1016/j.jcp.2013.02.041

    [6] D. R. Ernst, A. Bortolon, C. S. Chang, S. Ku et al., Phys. Rev. Lett. (2024) accepted for publication. https://arxiv.org/abs/2403.00185

  • Kinetic simulations by the DSBGK method with noise reduction (video)
    Jun Li, Agency for Science, Technology and Research (A*STAR), Singapore, abstract, slides
    [#s1786, 25 Apr 2024]

    Particle methods for kinetic simulations are numerically stable, easy to use, convenient for irregular geometries. The main challenge is the stochastic noise that could make the simulations unaffordable in low-speed problems (i.e., low signal-to-noise ratio) as well as transient problems where the time-averaging scheme is invalid. A lot of efforts have been made on the noise reduction by modifying the traditional direct simulation Monte Carlo (DSMC) method in solving the Boltzmann-like equations. At this presentation, I will introduce the direct simulation BGK (DSBGK) method, which can solve BGK-like equations as good approximations to the Boltzmann equation in many problems. As a duality of the DSMC method and the lattice Boltzmann method (LBM), the DSBGK method adopts a large number of simulated particles to represent the distribution function in the phase space, as used in the DSMC method but different from the LBM, while it updates the variables of each particle by integration of the kinetic equation along the corresponding trajectory, as modeled in the LBM but different from the DSMC method. The increments of particles’ variables inside each cell during each time step are obtained by the integration and the corresponding summations are used to regulate (not recompute) the macroscopic variables of the cell concerned, according to the mass, momentum and energy conservation laws. The previous values of cell’s variables are kept as anchors in the auto-regulation scheme to significantly reduce noise associated with the particles’ random movements into and out of each cell. Simulation results in several problems will be presented to show the noise reduction as well as the accuracy validation. Performance comparison with other particle and deterministic methods will be discussed.

  • Modeling interactions between runaway electrons and whistler waves in fusion plasmas (video)
    Yashika Ghai, ORNL, abstract, slides
    [#s1787, 18 Apr 2024]

    Resonant interactions between runaway electrons (REs) and whistler waves in a tokamak may lead to pitch angle scattering of the REs. An increase in RE pitch angles may give rise to the energy dissipation of the runaways via synchrotron radiation. DIII-D experiments on whistler waves have indicated a possibility of intentionally launching whistler waves to mitigate the deleterious effects of REs on the plasma facing components via resonant interactions with whistlers [1,2]. In present work, we have use the coupled KORC-AORSA model to numerically analyze the complex nature of the interactions between whistler waves and runaway electrons in DIII-D. In this framework, we follow full orbit trajectories of large RE ensembles using the Kinetic Orbit Runaway Electron (KORC) code in the presence of whistler wave fields calculated by All Orders Spectral Algorithm (AORSA) code in a DIII-D experimental equilibrium. The nature of RE transport (diffusive/non-diffusive) [3] is analyzed in the presence of whistler fields and the impact of whistler field amplitudes and frequencies is observed on the pitch angle scattering of REs. Our findings indicate a significant increase in RE energy and scattering of the runaways to large pitch angles for whistler fields exceeding a threshold amplitude. The coupled KORC-AORSA simulation model can be further used to get physical insights into tokamak experiments on whistler waves- REs interactions.

    [1] D. A. Spong et al., Phys. Rev. Lett., 120, 155002 (2018).

    [2] W. W. Heidbrink at al., Plasma Phys. Control. Fusion, 61, 14007 (2019).

    [3] D. del-Castillo-Negrete, Phys. Plasmas, 13, 082308 (2006).

  • The thermo-Alfvénic instability: from toy model to torus (video)
    Toby Adkins, University of Otago, abstract, slides
    [#s1799, 16 Apr 2024]

    A comprehensive understanding of electromagnetic effects on the microinstability properties of tokamak plasmas is becoming increasingly important as experimental values of the plasma beta and, therefore, electromagnetic fluctuations will be higher in reactor-relevant tokamak scenarios. Despite significant numerical progress in understanding the behaviour of instabilities such as the micro-tearing mode (MTM) or kinetic ballooning mode (KBM), there is still a lack of clarity about the fundamental physical processes that are responsible for them, owing to the complexity of full toroidal geometry. Constructing simplified models offers a path towards distilling the fundamental physical ingredients behind electromagnetic destabilisation. This talk focuses on electromagnetic instabilities driven by the electron-temperature gradient (ETG) in a local 'toy' model of a tokamak-like plasma. The model has constant equilibrium gradients (including magnetic drifts, but no magnetic shear) and is derived in a low-beta asymptotic limit of gyrokinetics. A new instability is shown to exist in the electromagnetic regime, the so-called 'thermo-Alfvénic instability' (TAI), whose physical mechanism hinges on a competition between diamagnetic drifts (due to the ETG) and rapid parallel streaming along perturbed field lines. Using linear gyrokinetic simulations, the TAI's presence is confirmed in slab geometry. The mapping of the TAI onto a more realistic tokamak equilibrium is considered, demonstrating that it survives aspects of the transition to toroidicity. A comparison is then drawn with the properties of the MTM and KBM, contextualising the TAI within the wider 'zoo' of electromagnetic instabilities commonly observed in tokamak simulations.

  • Resonances and Stellarator Design (video)
    Roscoe White, PPPL, abstract, slides
    [#s1780, 05 Apr 2024]

    High energy particle resonances play an important role in particle confinement in toroidal fusion devices, both tokamaks and stellarators. In stellarators a resonance that matches the periodicity of the equilibrium field produces islands in particle orbits which increase in size with particle energy and can induce loss. As demonstrated in the Japanese stellarator LHD, the presence of a high frequency resonance invariably gives rise to a strong Alfven mode that causes particle loss. In a nonsymmetric stellarator a resonance does not produce a well defined island structure in the orbits, but typically scatters ten percent of orbits of all energies and pitch randomly, modifying mode growth and saturation properties. Avoiding the presence of high energy particle resonances should be a part of device design.

  • Simulating electromagnetic instabilities in tokamaks with gyrokinetic particle-in-cell code GTS (video)
    Edward Startsev, PPPL, abstract, slides
    [#s1785, 28 Mar 2024]
    Recently, the numerical scheme presented in [1] enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma beta. This scheme resolved the long-standing "cancelation problem" that previously hindered gyrokinetic particle-in-cell code simulations of electromagnetic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode (KBM), the internal kink mode, the tearing mode, the micro-tearing mode (MTM) and toroidal alfven eigenmode (TAE) destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We also showcase the preliminary nonlinear simulations of the kinetic ballooning instability and (2,1) island formation due to the tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance. [1] A. Mishchenko, M. Cole, R. Kleiber, A. Konies, Phys. Plasmas 21 (2014) 052113.
  • Isotope Effects in Tokamak Plasmas: Core and Edge (video)
    Wei-li Lee, PPPL, abstract, slides
    [#s1753, 14 Mar 2024]

    Recently, a flurry of activities has been carried out on the isotopic effects in JT-60U [1], JET [2] and DIII-D [3], which has shown favorable confinement trend for heavier hydrogen isotopes. The consensus from these experimental observations was that this is an unsolved puzzle in tokamak plasmas. This is in fact not quite accurate. When the favorable effects was first observed on TFTR [4,5] for hydrogen, deuterium and tritium experiments, a theoretical attempt was indeed made to understand the results by Lee and Santoro [6]. Apparently, this paper has not attracted much attention in the community. Recently, a paper by Lee and White [7] on the H-mode physics has also described the isotope effects at the H-mode pedestal. In this talk, the theoretical interpretations on these isotopic effects based on 1) the resonance broadening theory [8] in the core as well as 2) the force balance equation for the pedestal from the gyrokinetic theory [7] will be described. The implementation of the related physics in an initial value code such as GTC [9] and/or GTS [10] will also be discussed.

    [1] H. Urano and E. Narita , Plasma Phys. Control. Fusion 63, 084003 (2021)

    [2] L. Horvath, C. F. Maggi, A. Chankin et al., Nuclear Fusion 61, 046015 (2021)

    [3] L. Schmitz, Phil. Trans. R. Soc. A381: 20210237 (2022)

    [4] S. D. Scott, M. C. Zarnstorff, C. W. Barnes, R. Bell et al., Phys. Plasmas 2, 2299 (1995)

    [5] S. D. Scott, G. W. Hammett, C. K. Phillips et al., IAEA-CN-64/A6-6 (1997)

    [6] W. W. Lee and R. A. Santoro, Phys. Plasmas 4, 169 (1997)

    [7] W. W. Lee and R. B. White, Phys. Plasmas 26, 040701 (2019)

    [8] T. H. Dupree, Phys. Fluids 11, 2680 (1968)

    [9] Z. Lin, T. S. Hahm, W. W. Lee et al., Science 281, 1835 (1998).

    [10] W. X. Wang, Z. Lin, W. M. Tang, W. W. Lee et al., 13, 092505 (2006)

  • Role of isotopes in microturbulence from linear to saturated Ohmic confinement regimes (video)
    Lei Qi, KFE, abstract, slides
    [#s1754, 07 Mar 2024]

    The first principle gyrokinetic numerical experiments investigating the isotopic dependence of energy confinement achieve a quantitative agreement with experimental empirical scaling, particularly in Ohmic and L-mode tokamak plasmas. Mitigation of turbulence radial electric field intensity |δEr|2 and associated poloidal δE×B fluctuating velocity with the radial correlation length l_cr ∝ Mi^0.11 strongly deviating from the gyro-Bohm scaling is identified as the principal mechanism behind the isotope effects. Three primary contributors are classified, the deviation from gyro-Bohm scaling, zonal flow and trapped electron turbulence stabilization. Zonal flow enhances isotope effects primarily through reinforcing the inverse dependence of turbulence decorrelation rate on isotope mass with ω_c ∝ Mi^-0.76, which markedly differs from the characteristic linear frequency. The findings offer insights into isotope effects, providing critical implications for energy confinement optimization in tokamak plasmas.

  • Turbulent Dynamo action in a 3-dimensional magnetohydrodynamic plasma
    Shishir Biswas, Institute for Plasma Research, abstract
    [#s1781, 27 Feb 2024]

    The generation of small scale, mean or large scale magnetic fields in cosmos and astrophysical bodies is an important problem in astrophysical plasmas. A possible mechanism behind these multi scale magnetic energy growth is explained via dynamo action. Shear flows [1] often coexist in astrophysical conditions and the role of flow shear on the onset of dynamo is only beginning to be investigated. The paradigm of investigation of the exponential growth of magnetic field caused by the interaction of small-scale velocity fluctuations and a flow shear; is commonly referred to as the “shear dynamo problem” [2]. Various laboratory experiments [3], as well as numerical studies have been performed to understand these astrophysical scenarios in detail. According to conventional understanding, for a large scale or mean field dynamo, a lack of reflectional symmetry (e.g., non-zero fluid or kinetic helicity) is required, where as for small scale or fluctuation dynamo it is not. Obviously the role of fluid or kinetic helicity on the onset of dynamo action is a sensible question to ask.

    In this present work we have analyzed kinematic dynamo model i.e, a case wherein (magnetic field does not back-react on velocity field) using a flow recently proposed by Yoshida and Morrison (YM) [4]. An interesting and useful aspect of this flow is that, it is possible to inject finite fluid helicity in the system, by systematically varying certain physically meaningful parameter. Using direct numerical simulation, we demonstrate that by systematically injecting finite fluid helicity, a systematic route emerges that connects “non-dynamo” to “dynamo” regime [5]. Time-averaged magnetic energy spectrum, for various magnitudes of injected fluid helicity is calculated and it is observed that, the spectra contain a visible maxima at a higher mode number, which is the distinguishing feature of small scale dynamo (SSD) [5]. However for a nonlinear dynamo or self-consistent dynamo model, the nonlinear effects start to change the flow (once the magnetic field is large enough) to stop further growth in magnetic field energy, i.e, the flow and magnetic field “back react” on each other. The influence of helical and non-helical drive in such a nonlinear or self-consistent dynamo model is shown to have some crucial dynamics [6]. Evidence of small-scale dynamo (SSD) activity is found for both helical and non-helical drives [6]. The spectrum analysis shows that the kinetic energy evolution adheres to Kolmogorov’s k^−5/3 law, while the magnetic energy evolution follows Kazantsev’s k^3/2 scaling. These scalings are observed to be valid for a range of magnetic Prandtl numbers (Pm) [6]. We have performed the above said studies using an in-house developed, multi-node, multi-card GPU based weakly compressible 3D Magnetohydrodynamic solver (GMHD3D) [7, 8]. Details of this study will be presented.

    References:

    [1] S. Biswas & R. Ganesh, Phys. Fluids 34, 065101 (2022).

    [2] S. Biswas & R. Ganesh, Phys. Plasmas 30, 112902 (2023).

    [3] R. Monchaux, M. Berhanu, et al., Phys. Rev. Lett. 98, 044502 (2007).

    [4] Z. Yoshida & P. J. Morrison, Phys. Rev. Lett. 119, 244501 (2017).

    [5] S. Biswas & R. Ganesh, Physica Scripta, Volume 98, Number 7.

    [6] S. Biswas & R. Ganesh, Manuscript under Preparation (2024).

    [7] S. Biswas, R. Ganesh et al. “GPU Technology Conference (GTC-2022)”,

    https://www.nvidia.com/en-us/on-demand/session/gtcspring22-s41199/.

    [8] S. Biswas & R. Ganesh, Computers and Fluids 272 (2024) 106207.

  • Rigorous Calculation of Transport-Driven Toroidal Rotation with Generally Varying Diffusivity (video)
    Timothy Stoltzfus-Dueck, PPPL, abstract, slides
    [#s1734, 15 Feb 2024]

    Future devices like ITER will have limited capacity to drive toroidal rotation, increasing the risk of instabilities like resistive wall modes. Fortunately, many experiments have found that tokamak plasmas rotate “intrinsically”, that is, without applied torque. The modulated-transport model shows that such rotation may be caused by the interaction of ion drift-orbit excursions with the strong spatial variation of the turbulent momentum diffusivity [1]. The model predicts intriguing qualitative behavior, such as a strong dependence of edge intrinsic toroidal rotation on the major-radial position of the X-point, which was subsequently measured on TCV [2]. The model has also been experimentally validated through further dedicated tests [3, 4], as well as via application in the new European whole-device transport model IMEP [5]. However, certain applications will require a relaxation of the underlying assumptions. In particular, the original model required the turbulent momentum diffusivity to decay exponentially in the radial direction, while experiments often exhibit a more complicated variation. In this work, we generalize the modulated-transport model to allow the turbulent momentum diffusivity to depend on space in an axisymmetric but otherwise arbitrary way. To enable this generality, we assume that the normalized diffusivity is weak, roughly equivalent to assuming that the pedestal-top ion transit time is short compared to the transport time across the pedestal, a condition that is almost always met for experimental applications. Given the increased flexibility, along with a technically much easier calculation, the new approach may serve as a basis for future extensions, including shaped geometry and trapped particles as well as the retention of momentum transport by neutrals.

    [1] T. Stoltzfus-Dueck, Phys. Rev. Lett. 108, 065002 (2012).

    [2] T. Stoltzfus-Dueck et al., Phys. Rev. Lett. 114, 245001 (2015).

    [3] J. A. Boedo et al., Phys. Plasmas 23, 092506 (2016).

    [4] A. Ashourvan, B. A. Grierson, D. J. Battaglia, S. R. Haskey, and T. Stoltzfus-Dueck, Phys. Plasmas 25, 056114 (2018).

    [5] T. Luda et al., Nucl. Fusion 61, 126048 (2021).

  • Fast-ion-driven vertical displacement oscillatory modes in tokamak plasmas (video)
    Tommaso Barberis, Polytechnic University of Turin, abstract, slides
    [#s1745, 09 Feb 2024]

    Axisymmetric modes in elongated plasmas are normally associated with a well-known ideal instability resulting in a vertical shift of the whole plasma column. This vertical instability is stabilized by means of passive feedback consisting of eddy currents induced by the plasma motion in a nearby wall and/or in plasma-facing components. When a thin resistive wall is considered, the n=0 mode dispersion relation can be studied analytically with reduced ideal MHD models and is cubic. Under relevant conditions, two roots are oscillatory and weakly damped. These oscillatory modes present Alfvénic frequency and are dependent on plasma elongation and on the relative position of the plasma boundary and of the wall. The third root is unstable and represents the so- called resistive wall mode (RWM) [1]. We focus on the two oscillatory modes, dubbed Vertical Displacement Oscillatory Modes (VDOM), that can be driven unstable due to their resonant interaction with energetic ions.

    The fast ion drive, involving MeV ions in present days tokamak experiments such as JET, may overcome dissipative and resistive wall damping, setting an instability threshold, as described in Ref. [2]. The effects of energetic particles are added within the framework of the hybrid kinetic-MHD model. An energetic ion distribution function with ∂F/∂E > 0 is required to drive the instability, achievable with pitch angle anisotropy or with an isotropic distribution in velocity space with regions of positive slope as a function of energy. The latter situation can be achieved by considering losses of fast ions or due to fast ion source modulation [3-4]. The theory presented here is partly motivated by the observation of saturated n=0 fluctuations reported in [4,5], which were initially interpreted in terms of a saturated n=0 Global Alfvén Eigenmode (GAE). Modeling of recent JET discharges using the NIMROD [6] extended-MHD code will be presented, focusing on mode structure and frequency dependence. It is early for us to conclude whether the mode observed at JET is a VDOM rather than a GAE, nevertheless, we discuss the main points of distinction between GAE and VDOM that may facilitate their experimental identification.

    References

    [1] T. Barberis, et al. 2022, J.Plasma Phys. 88, 905880511

    [2] T. Barberis, et al 2022 Nucl. Fusion 62 06400

    [3] Ya.I. Kolesnichenko and V.V. Lutsenko 2019 Nucl. Fusion 59 126005

    [4] V. G. Kiptily et al 2022 Plasma Phys. Control. Fusion 64 064001

    [5] H. J. C. Oliver et al. 2017 Phys. Plasmas 24, 122505

    [6] C. Sovinec et al. and the NIMROD Team 2004 J. Comp. Phys. 195 355

  • Transport Physics of the Density Limit (video)
    Patrick Diamond, Isaac Newton Institute, Cambridge U.K. and U.C. San Diego, abstract, slides
    [#s1743, 08 Feb 2024]

    We discuss recent progress in understanding the role of transport physics in density limit phenomena. Our approach is one which combines theory and experiment. Contrary to the conventional wisdom that the density limit is enforced by MHD instability, findings indicate that the L−mode density limit is associated first with the degradation of the edge E × B shear layer. The latter occurs for k‖² vₜₕₑ² / ω νₑᵢ<1. Shear layer decay leads to strongly enhanced turbulence spreading and increased production of density 'blobs'. Interestingly, the spreading flux increases more rapidly with increasing n / n_G than does the particle flux. Shear layer decay is linked to a decline in zonal flow production.

    A simple model for flow, fluctuation and density evolution reveals that the edge density will increase with edge heat flux (power). This favorable trend results from increased Reynolds stress flow drive at higher power. It provides physical insight into the power scaling of density limit, now observed in experiments. A scaling of n ∼ P^(1/3) is suggested for the case of ITG turbulence.

    We briefly discuss recent density limit experiments in negative triangularity plasmas, as well as aspects of the H−mode density limit phenomenon. Implications for burning plasma are discussed.

    Contributions from Ting Long, SWIP ; Rameswar Singh ,UCSD ; Rongjie Hong, UCLA and DIII−D ; Zheng Yan, Univ Wisc and DIII−D ; and George Tynan, UCSD are acknowledged.

  • Suppression of temperature-gradient-driven turbulence by perpendicular flow shear (video)
    Plamen Ivanov, Oxford University, abstract, slides
    [#s1746, 07 Feb 2024]

    Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalized by particle collisions, we formulate a detailed phenomenological theory for the influence of perpendicular flow shear on magnetized-plasma turbulence. Our theory introduces two distinct regimes, called the weakly and strongly sheared regimes, each with its own set of scaling laws for the scale and amplitude of the fluctuations and for the level of turbulent heat transport. We discover that the ratio of the typical radial and poloidal wavenumbers of the fluctuations, i.e., their aspect ratio, plays a central role in determining the dependence of the turbulent transport on the imposed flow shear. Our theoretical predictions are found to be in excellent agreement with numerical simulations of two models of magnetized plasma turbulence: (i) an electrostatic fluid model of slab electron-scale turbulence, and (ii) Cyclone-base-case gyrokinetic ion-scale turbulence.

  • Exploring plasma turbulence with a gyrokinetic moment-based approach (video)
    Antoine Hoffmann, EPFL, abstract, slides
    [#s1750, 05 Feb 2024]

    We present the first local delta-f nonlinear gyrokinetic (GK) simulations based on a gyro-moment (GM) approach, which exploits the projection of the distribution functions onto a Hermite-Laguerre velocity- space basis. We first demonstrate that, in contrast to gyrofluid models, the GM approach reproduces the Dimits shift, notably, with a coarser velocity space resolution than the continuum GK GENE code. In addition, we reveal that the choice of collision operator model (Dougherty, Sugama, Lorentz and Landau) significantly impacts the level of turbulent transport through multi-species zonal flow damping.

    In addition, we show for the first time that the GM approach is able to bridge the gap between GK and reduced fluid modelling by its exact equivalency to the model of Ivanov et al. 2020 when considering the same limits. Leveraging its efficiency and multi-fidelity capability, we finally use the GM approach to explore the impact of triangularity in realistic DIII-D edge conditions across a range of models, spanning from GK electron-ion multi-scale simulations to the reduced fluid limit.

  • Pedestal turbulence in AUG and JET from a global gyrokinetic perspective (video)
    Leonhard Leppin, IPP, Garching, abstract, slides
    [#s1749, 02 Feb 2024]

    Turbulence is one of the key ingredients in shaping H-mode pedestals. Identifying the relevant turbulent transport mechanisms in a pedestal, however, is a great scientific and numerical challenge. Here, we address this challenge by global, nonlinear gyrokinetic simulations of two pedestals: One from ASDEX Upgrade (Type-I ELMy H-mode) and one from JET (hybrid scenario H-mode). The global simulations permit to calculate heat fluxes due to ion-scale turbulence in the steep gradient region encompassing the full pedestal from top to foot. They are supported by detailed characterizations of gyrokinetic instabilities via local, linear simulations at pedestal top, center and foot as well as dedicated nonlinear electron-scale heat flux calculations. Simulations are performed with the gyrokinetic, Eulerian, delta-f code GENE (genecode.org) and employ a new code upgrade of its global, electromagnetic model that enables stable simulations at experimental plasma beta values.

    In both investigated pedestals from AUG and JET, we find turbulent transport to have a complex radial structure that is multi-scale and multi-channel. Electron transport in the AUG pedestal is found to transition in scale. At the pedestal top ion-scale TEM/MTM instabilities fuel electron transport whereas in the pedestal center electron-scale ETG transport takes over. Turbulent ion heat flux is present at the pedestal top and strongly reduces towards the steep gradient region. Magnetic shear is found to locally contribute to the stabilization of microinstabilities and reduction of heat flux. In the JET pedestal, transport due to ITG is found to play a much more important role, particularly on the pedestal top/ outer core. In both pedestals, ExB shear is confirmed to strongly reduce heat fluxes in the global, nonlinear simulations. We discuss implications of our results for the applicability of quasi-linear transport models in the pedestal.

  • Theoretical scaling of the operational density limit in tokamaks and comparison to experimental data (video)
    Paolo Ricci, EPFL, abstract, slides
    [#s1736, 24 Jan 2024]
    A validated and predictive first-principles scaling of the operational density limit is presented. The scaling is based on consideration of plasma turbulence in the tokamak boundary supported by the results of first-principles simulations. The scaling is validated against a multimachine database that includes results from the AUG, JET and TCV tokamaks. By revealing a dependence of the operational density limit on the power crossing the separatrix, the result we obtain has consequences for ITER operation and the design of future fusion reactors.
  • Linear Stability of a weakly magnetized rotating plasma column (video)
    Surabhi Aggarwal, Aix-Marseille University, abstract, slides
    [#s1727, 07 Dec 2023]

    Understanding the formation of large-scale structures in weakly magnetized plasmas represents a crucial step towards developing predictive design capabilities for E×B devices dedicated to investigating fundamental plasma physics phenomena. MISTRAL is such a device based at PIIM laboratory to study plasmas in cross-field configuration (E⊥B). The formation of coherent rotating structures in MISTRAL is supposed to be due to an interplay between various instabilities and the E×B flow. However, a definitive understanding of which instabilities are accountable for their emergence and the specific triggers involved remains elusive. An experimental investigation of MISTRAL plasmas has been performed to lay the basis for the theoretical modeling. A two-fluid model has been developed to discuss the linear stability of rotating plasma columns. Prior works have demonstrated that rotating plasma columns are susceptible to centrifugal flute modes. However, most of the existing models rely on the low-frequency approximation (LFA), which holds true when the instability frequency and equilibrium flow frequency are considerably smaller than the ion-cyclotron frequency. This assumption is challenged in numerous laboratory plasma devices, including weakly magnetized plasma columns like MISTRAL. To address this limitation, a radially global dispersion relation describing the centrifugal instability without the LFA has been derived and linear stability analysis is performed. A comparison has been made between the results obtained using the dispersion relation with the radially local approximation and those obtained using the radially global dispersion relation. This comparison revealed the non-applicability of the local solution to MISTRAL-like plasma systems. Due to the high fraction of neutrals in the present plasma system, the model is further extended to include the effects due to ion-neutral collisions. In this first step, the ion-neutral collision frequency is assumed to be small as compared to the ion-cyclotron frequency. The dispersion relation is then solved with finite ion-neutral collisionality and the linear stability analysis is conducted.

  • Compile-time Symbolic Solver for the Gyrokinetic-MHD Hybrid Code GMEC and 3D Field and Particle code FP3D (video)
    Peiyou Jiang, Zhejiang University, abstract, slides
    [#s1726, 05 Dec 2023]

    In magnetic confinement fusion plasmas, many instabilities have a flute mode character. The field-aligned coordinates bring the benefit of efficient resolution of parallel mode structure along the magnetic field direction. However, the curvilinear coordinates make equations and codes more complex especially in high order PDE.

    The Compile-time Symbolic Solver (CSS) is developed to solve PDEs and ODEs in finite difference method from vector equations directly. CSS is a general-purpose finite difference framework for generating finite difference codes easily and greatly reducing the risk of implementation mistakes.

    For physics model, CSS supports arbitrary equations in arbitrary curvilinear coordinates and multiple boundaries for both PDEs and ODEs. For memory distribution, N-dimension distribution grids with hybrid TBB and MPI parallelization in arbitrary dimensions are implemented. For numerical method, CSS employs Method of Line in numerical difference with arbitrary grid points and offset. The N-dimension B-spline is implemented with arbitrary orders for pushing particles. CSS employs PARDISO to solve matrix problem and Runge–Kutta method for time advance. CSS is a C++20 template metaprogramming code which guarantee zero-overhead at runtime. Furthermore, the instruction optimization makes the codes generated by CSS much faster than usual codes.

    We have used CSS to generate the Gyrokinetic-MHD Hybrid Code GMEC, 3D Field and the Particle calculation code FP3D and a fluid ITG code. For GMEC, we propose a new shifted metric method which is able to stabilize numerical instabilities and avoid the interpolation from MHD field-align grids to particle flux coordinate grids at the same time. The equilibriums can be analytical ones or numerical ones calculated by VMEC or DESC. We have used GMEC to simulate ballooning modes (IBM) with or without the diamagnetic drift term and tearing modes. The simulation results agree well with those of the eigenvalue code MAS. The n=20 IBM costs only 17 seconds using 448 cores. We have also used GMEC to simulate energetic particle-driven TAEs in a circular equilibrium and a CFETR equilibrium. The results of an n=3 TAE agree well with those of M3D-K code.

    We have also used CSS to generate the test particle code FP3D for calculation of magnetic surfaces, rotation transform, particle orbits and neoclassical transport in both tokamaks and stellarators. We have used FP3D to simulate ripple losses in EAST tokamak and neoclassical transport coefficient in NCSX. The results are consistent with previous results. FP3D has been used in design and optimization of stellarators successfully.

    [1] P. Y. Jiang, et al. CSS: Compile-time symbolic solver for finite difference method. To be submitted.

    [2] P. Y. Jiang, et al. GMEC: Gyrokinetic-MHD Hybrid Code. To be submitted.

    [3] P. Y. Jiang, Z. C. Feng, G. D. Yu, and G. Y. Fu, FP3D: A code for calculating 3D magnetic field and particle motion. Submitted to POP.

  • Multiscale simulations of radiofrequency wave scattering in the scrape-off layer (video)
    Bodhi Biswas, University of York, abstract, slides
    [#s1725, 30 Nov 2023]

    Accurately predicting lower hybrid current drive (LHCD) in the weak-damping regime is an outstanding challenge, which suggests important physics is missing in present-day ray-tracing/Fokker-Planck (RTFP) models. In this work, the impact of filamentary scrape-off layer (SOL) turbulence on LH waves is investigated using a new multi-scale scattering model. When coupled to an RTFP code, the resulting simulations of LHCD in Alcator C-Mod show RF power deposition profiles robustly peaked on-axis, leading to good agreement with experimental Motional Stark Effect and hard X-ray measurements. Therefore, it is shown that the rotation of the perpendicular wave-vector due to SOL turbulence is sufficient to bridge the discrepancy between simulation and experiment. Notably, this model predicts an asymmetric broadening of the transmitted wave-spectrum, which is attributed to full-wave scattering effects in the presence of spatially coherent turbulence. This asymmetry leads to rotation of incident power away from the plasma core when SOL densities are sufficiently high. RTFP modeling shows this effect plays a significant role in the anomalous drop in LHCD efficiency observed at high densities.

    The multi-scale scattering model has two steps. (1) Single filament-wave interactions are solved in full-wave formalism using a Mie-scattering technique. (2) Multiple of these filament-wave interactions are modeled using the radiative transfer approximation, in which a photon’s scattering probability depends on the statistical properties of the filament population. The radiative transfer equation (RTE) is then solved using a Monte Carlo scattering term in a ray-tracing model, allowing for self-consistent coupling to RTFP codes. For verification and comparison against other models, the RTE is also solved in a simple slab geometry using a Markov chain. This model shows good agreement with ray-tracing in the Wentzel-Kramer-Brillouin (WKB) limit, and predicts greater, asymmetric scattering beyond the WKB limit. Good agreement is also found with numeric full-wave solutions at sufficiently low filament packing-fraction, which is consistent with the validity limit of the radiative transfer approximation.

    It should be emphasized that this multi-scale scattering model retains many important full-wave effects while remaining computationally inexpensive, allowing fast parameter scans and inter-shot analysis. In addition, this model is highly applicable to the modeling of electron cyclotron wave scattering since the radiative transfer approximation is increasingly valid for waves at higher k.

    B. Biswas et al., “Spectral broadening from turbulence in multiscale lower hybrid current drive simulations,” Nuclear Fusion, 63, 1 (2022).

    B. Biswas et al., “A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments,” Journal of Plasma Physics, 87, 5 (2021).

  • Parity broken fluids: Theory and Applications (video)
    Dylan J Reynolds, CUNY, abstract, slides
    [#s1713, 05 Oct 2023]
    Odd viscosity is a dissipationless transport coefficient that arises is certain classes of parity broken fluids. While typically studied in 2D, the generalization to 3D parity broken flows leads to a much wider class of transport coefficients. In the context of plasmas this parity breaking manifests as gyro viscosity, and is associated with viscous stresses perpendicular to the flow. In this talk I will take a phenomenological approach and outline which transport coefficients are allowed by the symmetry of isotropic 2D and 3D flows. I will then look at some specific systems that highlight a subset of these coefficients, in particular parity odd Hele-Shaw flows and ferrofluids. I will also discuss some microscopic mechanisms that lead to parity broken flows, and look at some applications to active matter and condensed matter systems.
  • Development of the Gyrokinetic-MHD Energetic Particle Code GMEC
    Guoyong Fu, Zhejiang University, abstract
    [#s1712, 01 Sep 2023]
    We report the status of the Gyrokinetic-MHD Energetic particle hybrid Code GMEC being developed for simulations of energetic particle (EP)-driven Alfven instabilities and EP transport in magnetic fusion plasmas such as ITER. In the hybrid model, electrons are treated as a fluid, EPs and thermal ions are described by gyro-kinetic equations. The energetic particle effects enter in the gyrokinetic vorticity equation via pressure terms which are obtained by solving the gyrokinetic equations using PIC method. The field-aligned coordinates and meshes are used to efficiently resolve mode structures of high-n Alfven modes. Five-points 4th order finite differences and 4th order Runge-Kutta method are used for numerical differentiations and time advance respectively. The Compile-time Symbolic Solver (CSS) is developed to generate coding from vector equations directly. CSS is a C++20 template metaprogramming code. It expands vector equations into components scalar equations at compile-time, and greatly simplifies coding of differential equations in toroidal curvilinear coordinates. Both MPI and TBB are used for parallelization. Up to now, a simplified version of GMEC has been developed with initial verifications for ideal ballooning modes and EP-driven TAEs. The alpha particle-driven Alfven eigenmodes in the Chinese Fusion Engineering Test Reactor (CFETR) have also been simulated successfully. Details of GMEC and its applications will be presented.
  • Intrinsic toroidal rotation in tokamak plasmas driven by turbulent and neoclassical processes from global total-f gyrokinetic simulations (video)
    Hongxuan Zhu, abstract
    [#s1680, 27 Jul 2023]
  • Current status of XGC-S and its applications to the Large Helical Device
    Toseo Moritaka, National Institute for Fusion Science, Japan, abstract
    [#s1686, 20 Jul 2023]
    XGC-S is the stellarator version of the global gyrokinetic code XGC, originally developed for whole-volume modeling of tokamaks. We will cover the following topics: 1. The development history and basic code descriptions, 2. A recent application to the Large Helical Device (LHD), and 3. Future capabilities. We consider isotope effects in LHD under the influence of radial electric field and heavy hydrogen components. Both radial electric fields and heavy hydrogen components have similar impacts on thermal conductivity, involving the elongation of the mode structure to increase the heat flux and mode suppression to decrease the heat flux. Quasi-linear estimations indicate that these competitive effects lead to a favorable mass number dependency. Finally, we will present the future capabilities of XGC-S, including new meshing schemes, an electrostatic field solver that could handle complicated magnetic fields in the edge region, extension to a total-f method and electromagnetic physics, and GPU offloading. The electrostatic field structure may explain the up-down asymmetry of divertor particle flux observed in LHD. The advanced XGC features, such as kinetic electrons and multi-species collision, are yet to be implemented in XGC-S. These are also promising to address important issues in LHD.
  • Stellarator profile prediction and evolution using the Trinity3D+GX multi-scale gyrokinetic framework
    Noah Mandell, abstract
    [#s1693, 13 Jul 2023]
    Trinity3D+GX is a framework that leverages multi-scale gyrokinetic theory to model macro-scale profile evolution in fusion plasmas (tokamaks and stellarators) due to micro-scale turbulent processes. In this talk I will first provide a brief background on the multi-scale gyrokinetic theory underpinning the model. I will then discuss the GX gyrokinetic code, which has been developed as a GPU native code that uses an efficient pseudo-spectral discretization scheme to target fast turbulence calculations for fusion reactor design and optimization. This enables GX to be embedded as the micro-turbulence model in the Trinity3D transport solver for tractable fusion profile prediction (and evolution) calculations. I will highlight some preliminary results of modeling W7X plasmas with the Trinity3D+GX system and discuss future plans for using the framework in experimental studies as well as stellarator FPP design and optimization.
  • Recent progress in the influence of energetic particles on tearing modes (remote)
    Huishan Cai, University of Science and Technology of China, abstract
    [#s1687, 06 Jul 2023]
    The dynamics of energetic particles and tearing modes and the interactions between them are of great significance for magnetically confined fusion plasmas. In this review, we focus on the issue: the influence of energetic particles on tearing modes. The influence of energetic particles on tearing modes is described on the basis of a general dispersion relation for tearing modes. The effects of energetic particles are considered separately in the outer region and the island region of a tearing mode.The physics mainly results from themodification of the perturbed parallel current by energetic particles without wave–particle resonance. In addition, the resonance between energetic particles and tearing modes is also reviewed. Our descriptions of physical phenomena here are based on an analytical approach, while the experiments and simulations are used to illustrate and confirm our results. Finally, a number of open issues are discussed.
  • Chaos and Magnetic Reconnection in Low Resistivity Plasmas
    Allen Boozer, Columbia University, abstract
    [#s1679, 22 Jun 2023]
    The mathematics and physics of each of the three aspects of magnetic field evolution—topology, energy, and helicity—is remarkably simple and clear. When the resistivity η is small compared to an imposed evolution, a/v, timescale, which means Rm ≡ μ0va/η >> 1, magnetic field line chaos dominates the evolution of field-line topology in three-dimensional systems. Chaos has no direct role in the dissipation of energy. A large current density, jη ≡ vB/η, is required for energy dissipation to be on a comparable time scale to the topological evolution. Nevertheless, chaos plus Alfv ́en wave damping explain why both timescales tend to be approximately an order of magnitude longer than the evolution timescale a/v. Magnetic helicity is injected onto tubes of field lines when boundary flows have vorticity. Chaos can spread but not destroy magnetic helicity. Resistivity has a negligible effect on helicity accumulation when Rm >> 1. Helicity accumulates within a tube of field lines until the tube erupts and moves far from its original location.
  • Transport-Driven Toroidal Rotation with General Viscosity Profiles
    Timothy Stoltzfus-Dueck, abstract
    [#s1661, 08 Jun 2023]
    Future devices like ITER will have limited capacity to drive toroidal rotation, increasing the risk of instabilities like resistive wall modes. Fortunately, many experiments have found that tokamak plasmas rotate “intrinsically”, that is, without applied torque. The modulated-transport model shows that such rotation may be caused by the interaction of ion drift-orbit excursions with the strong spatial variation of the turbulent momentum diffusivity [1]. The model predicts intriguing qualitative behavior, such as a strong dependence of edge intrinsic toroidal rotation on the major-radial position of the X-point, which was subsequently measured on TCV [2]. The model has also been experimentally validated through further dedicated tests [3, 4], as well as via application in the new European whole-device transport model IMEP [5]. However, certain applications will require a relaxation of the underlying assumptions. In particular, the original model required the turbulent momentum diffusivity to decay exponentially in the radial direction, while experiments often exhibit a more complicated variation. In this work, we generalize the modulated-transport model to allow the turbulent momentum diffusivity to depend on space in an axisymmetric but otherwise arbitrary way. To enable this generality, we assume that the normalized diffusivity is weak, roughly equivalent to assuming that the pedestal-top ion transit time is short compared to the transport time across the pedestal, a condition that is almost always met for experimental applications. Given the increased flexibility, along with a technically much easier calculation, the new approach may serve as a basis for future extensions, including shaped geometry and trapped particles as well as the retention of momentum transport by neutrals. [1] T. Stoltzfus-Dueck, Phys. Rev. Lett. 108, 065002 (2012). [2] T. Stoltzfus-Dueck et al., Phys. Rev. Lett. 114, 245001 (2015). [3] J. A. Boedo et al., Phys. Plasmas 23, 092506 (2016). [4] A. Ashourvan, B. A. Grierson, D. J. Battaglia, S. R. Haskey, and T. Stoltzfus-Dueck, Phys. Plasmas 25, 056114 (2018). [5] T. Luda et al., Nucl. Fusion 61, 126048 (2021).
  • Kinetic theory of performance limiting instabilities in tokamak plasmas Neoclassical Tearing Modes and micro-instabilities (Kink and Kinetic-Ballooning Modes) ,Video
    Alexandra Dudkovskaia , abstract
    [#s1678, 01 Jun 2023]
    Neoclassical tearing modes (NTMs) are identified as one of the main performance limiting, resistive MHD instabilities that exist in tokamak plasmas. They impose a limit on fusion gain, as well as plasma confinement time. Resulting from filamentation of the current density that flows through the tokamak plasma, they modify the equilibrium magnetic topology, breaking down the tokamak toroidal symmetry by forming a chain of magnetic islands. For large magnetic islands (i.e. much larger than the banana orbit width of trapped ions, ρbi), the plasma thermal pressure gradient is removed across them. This therefore reduces the total core plasma pressure, degrading tokamak confinement. The pressure profile flattening across the islands generates a hole in the bootstrap current density close to the island centre (”O-point”), enhancing the current density filamentation even further and amplifying the island. This provides the main drive for magnetic island growth. Fortunately, the NTM behaviour significantly differs from the above when magnetic islands are small, i.e. their width is comparable to ρbi. Experimental observations found that there is some threshold magnetic island width, 2wc ≈ (2 − 3)ρbi [1, 2], below which the pressure gradient is partially restored inside the island, providing its ”self-healing”. This wc is a key parameter of the NTM theory and is responsible for quantifying the NTM control system [3]. A novel drift island formalism is derived in [4] to quantify the NTM threshold in a low beta, large aspect ratio tokamak plasma. Reference [5] improves this theory further by introducing plasma shaping and finite beta effects. In particular, it is found that (1) a higher triangularity plasma is more prone to NTMs (in agreement with tearing mode onset relative frequency measurements in DIII-D, 2022), (2) the conventional (ε 1/2 , where ε is the tokamak inverse aspect ratio) NTM threshold dependence on the tokamak inverse aspect ratio is revisited for finite aspect ratio and (3) the NTM threshold dependence on poloidal beta is obtained and successfully benchmarked against the EAST threshold island width measurements (2022). Effects of the background electric field on the NTM threshold are also investigated [6]. While NTMs are one of the main sources of confinement limit in a tokamak, the core pressure is also significantly influenced by the pedestal physics. A nonlinear electromagnetic global gyrokinetic theory is derived in [7] to ensure that the effects associated with sharp pressure gradients and the conse- quential high bootstrap and Pfirsch-Schluter currents are fully captured in the nonlinear electromagnetic gyrokinetic theory, while allowing arbitrary magnetic field configurations and finite orbit width effects and ensuring consistent ordering. A reduced version of this theory (Bθ ≪ B0, where B0 is the total equi- librium magnetic field and Bθ is its poloidal component) has been implemented in the local turbulence code GS2 (to be referred to as NEO GS2) to quantify the impact of higher order gyrokinetics in sharp pressure gradient regions where the bootstrap current becomes large (such as the pedestal plasma and a spherical tokamak core plasma) [8]. The dominant impact is found to be on kinetic-ballooning modes (KBMs). In particular, it is found that the KBM growth rate is significantly suppressed by inclusion of neoclassical equilibrium effects at large density gradients, representative of the tokamak pedestal val- ues. The latter appears only when the neoclassical electrostatic potential, Φ 1 0 = Φ 1 0 (ψ, θ), dependent on poloidal angle, θ, is calculated consistently with plasma quasi-neutrality. Electrostatic modes are also found to be impacted by the neoclassical equilibrium physics. In contrast, the impact on micro-tearing modes (MTMs) is found to be minimal, based on the test cases considered. References [1] Z. Chang, J. D. Callen, E. D. Fredrickson et al. Phys. Rev. Lett. 74 (1995) 4663 [2] R.J. La Haye, R.J. Buttery, S.P. Gerhardt et al. Phys. Plasmas 19 (2012) 062506 [3] E. Poli, C. Angioni, F.J. Casson et al. Nucl. Fusion 55 (2015) 013023 [4] A V Dudkovskaia, J W Connor, D Dickinson, P Hill, K Imada, S Leigh and H R Wilson Plasma Phys. Control. Fusion 63 (2021) 054001 [5] A V Dudkovskaia, L Bardoczi, J W Connor, D Dickinson, P Hill, K Imada, S Leigh, N Richner, T Shi and H R Wilson Nucl. Fusion 63 (2023) 016020 [6] A V Dudkovskaia, J W Connor, D Dickinson, P Hill, K Imada, S Leigh and H R Wilson Drift kinetic theory of neoclassical tearing modes in tokamak plasmas: polarisation current and its effect on magnetic island threshold physics, to be submitted, (2023) [7] A V Dudkovskaia, H R Wilson, J W Connor, D Dickinson, F I Parra Plasma Phys. Control. Fusion 65 (2023) 045010 [8] A V Dudkovskaia, J W Connor, D Dickinson, H R Wilson Plamsa Phys. Control. Fusion 65 (2023) 054006
  • Towards a reduced transport model for microtearing turbulence in H-mode plasmas
    Myriam Hamed , abstract
    [#s1662, 18 May 2023]
    In magnetic confinement fusion research, predicting turbulent transport in tokamak edge plasma and its effect on fusion device operation is crucial for determining confinement properties. In order to better understand the causes and evolution of electron heat transport in tokamak discharges, a quasilinear transport model has been developed for use in integrated predictive modeling studies. Recent analyses of H-mode plasmas [1] suggest that small-scale instabilities localized near the rational surface, such as microtearing (MT) modes, have a significant effect on con- finement. MT modes draw on the electron temperature gradient as a free-energy source and rearrange magnetic topology through the creation of ion-Larmor-radius-scale magnetic islands, thereby playing a role in determining pedestal characteristics. The stability of MT modes has been extensively studied theoretically, showing that a slab current sheet is stable in the absence of collisions [2]. To evaluate the parametric dependencies of MT and determine new saturation rules, a reduced kinetic transport model for MT has been developed using an electromagnetic quasilinear theory. This reduced model solves the Vlasov and Maxwell equations, and its evaluation inside the resistive layer is obtained from a system of two equations linking the magnetic vector potential and the electric potential. To solve numerically, this system of equations, an eigenvalue code has been developed. The reduced transport model has been tested and compared with gyrokinetic simulations using JET experimental data, showing good agreement. Analysis of nonlinear gyrokinetic simulations shows that this quasilinear transport model for microtearing repro- duces gyrokinetic trends for a variety of parameter regimes[3]. The impact of the electric potential on nonlinear saturation is examined using this model. The electric potential plays a key role in microtearing destabilization by boosting the growth rate of this instability in the presence of collisions. Instability and saturation physics are exam- ined for different pedestal cases and radial positions, with a special focus on the role of electric field fluctuations and the role of zonal flows and fields. In the saturated state, it is found that removing electrostatic fluctuations causes a flux increase, whereas linear stabilization had been observed. This is consistent with a change in saturation mechanism from temperature corrugations to zonal-field and zonal-flow-based energy transfer. References [1] D.R. Hatch et al., Nucl. Fusion 56: 104003 (2016). [2] M. Hamed, et al., Physics of Plasmas, 26(9): 092506 (2019). [3] M. Hamed, et al., Physics of Plasmas 30: 042303 (2023)
  • TBA
    Robert Brzozowski, abstract
    [#s1636, 15 May 2023]
    TBA
  • Confinement, Transport, and Burning Plasma in Magnetic Fusion
    Hyeon K. Park, UNIST Korea, abstract
    [#s1659, 27 Apr 2023]
    Burning plasma in magnetic fusion relies on sufficient a-power and adequate energy confinement time. While D-T experiments with high ion temperature regime (Ti>Te) produced a-power up to ~4 MW and tE~0.7s with an ion heating system, the present ITER, expecting ~5s of tE, has only electron heating systems. No high-power electron heating systems can substitute a-power to identify the threshold electron heating power to achieve discharges with Ti>10 keV due to ion temperature clamping and practical problems (e.g., antennas and narrow resonance layer). In the confinement section, the difference between L and H mode (ETB) is attributed to the configuration difference and the edge density is largely controlled by influx plasmas induced by outflux plasmas from the limiter/divertor plates. The influx plasma is not a quiescent one and the fluctuation level should be high. The ITB position of the various improved confinement regimes such as “Supershot” and “Super H-mode” is highly correlated with the heating profile footprints. Note that ITB positions are formed where the fluctuations are minimum. Transport models that supported the physics of ITB and/or ETB are examined including ExB shear, ITG marginality, etc. A reasonably large size (Vp<200m3 is much smaller than ITER, Vp~800m3) ignition device is feasible with optimized ion heating system and device geometric factors.
  • Plasma flow and acceleration in the magnetic mirror/nozzle configurations
    Andrei Smolyakov, U. Saskatchewan, abstract
    [#s1660, 20 Apr 2023]
    Magnetic nozzle/mirror configurations with converging-diverging magnetic field are used in numerous plasma applications for fusion, electric propulsion and material processing. In fusion devices such as open mirrors and tokamak divertors, the diverging magnetic field is use to reduce the thermal loads to the walls. In electric propulsion, magnetic nozzle is employed to convert the plasma thermal energy into the kinetic energy of the directed flow producing the thrust. Ion acceleration by the magnetic nozzle is used in ion plasma sources for material processing. In many of these applications, plasma flow and acceleration share many common patterns. We present the results of the fluid model taking into account the effects of anisotropic ion pressure. Further generalization includes the role of the induced azimuthal magnetic field and plasma rotation, i.e., coupling with Alfven wave dynamics. It is shown that the inhomogeneous magnetic field couples the axial plasma flow with the evolution of the azimuthal magnetic field and plasma rotation resembling the problem of the magnetically driven flow in astrophysical jets and winds. The kinetic effects have been investigated using the quasineutral hybrid model with kinetic ions and isothermal Boltzmann electrons and full kinetic model including the ions and electrons in quasi-two dimensional (paraxial) model.
  • Sparse Approach for the Particle-In-Cell algorithm acceleration: application to low temperature plasmas (remote)
    Laurent Garrigues, U. Toulouse, abstract
    [#s1629, 20 Apr 2023]
    Explicit Particle-In-Cell (PIC) used to model low temperature plasmas are time consuming. We show that a new Sparse PIC method based on sparse grid approaches combined with the combination technique offers a promising alternative to reduce the computational time maintaining a high accuracy in the modeling results.
  • Shifting and splitting of resonance lines due to dynamical friction in plasmas ,Video
    Vinicius Duarte, abstract
    [#s1628, 06 Apr 2023]
    A quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally-unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics. Reference: V. N. Duarte et al, Phys. Rev. Lett. (2023) https://doi.org/10.1103/PhysRevLett.130.105101.
  • How TRANSP is evolving to meet the needs for ITER and FPP Design and new opportunities for collaboration Video,
    Francesca Poli, abstract
    [#s1635, 30 Mar 2023]
    The TRANSP framework has been a workhorse for the fusion community for over three decades. To meet the needs for HPC capabilities, preparation to ITER operation and FPP design, TRANSP has undergone substantial re-factoring and modernization during the past four years. Currently, TRANSP is being adapted to use the ITER Data Model, which opens the framework to a number of additional physics models. Opportunities for new collaborations are discussed, including ideas for gradually expanding TRANSP to model non-axisymmetric effects.
  • Quantum computing for modeling linear waves in plasmas
    Ivan Novikau, abstract
    [#s1585, 23 Mar 2023]
    Quantum computing (QC) is gaining attention as a potential way to speed up simulations of physical systems. The main efforts are currently focused on modeling quantum systems, e.g. for chemistry and material science. However, applying QC to classical problems, and plasma physics in particular, can also be beneficial. In this talk, I will discuss quantum modeling of linear radiofrequency (RF) waves, which in the future could improve the accuracy and resolution of RF simulations for fusion applications. In the first part of my talk, I will describe an algorithm for solving the initial-value problem for the propagation of RF waves in inhomogeneous cold magnetized plasma using so-called Quantum Signal Processing (QSP) [1]. In the second part of my talk, I will describe an algorithm for solving the boundary-value problem for dissipative linear waves propagating in a medium with a prescribed inhomogeneous dielectric permittivity using so-called Quantum Singular Value Transform (QSVT) [2]. [1] I. Novikau, E. A. Startsev, and I. Y. Dodin, Quantum signal processing for simulating cold plasma waves, Phys. Rev. A 105, 062444 (2022). [2] I. Novikau, I. Y. Dodin, and E. A. Startsev, Simulation of linear non-Hermitian boundary-value problems with quantum singular value transformation, arxiv:2212.09113.
  • (to be rescheduled)
    Priyanjana Sinha, abstract
    [#s1587, 16 Mar 2023]
    TBA
  • (to be rescheduled)
    Zhenyu Wang, abstract
    [#s1586, 09 Mar 2023]
    TBA
  • (to be rescheduled)
    Pallavi Trivedi, abstract
    [#s1584, 02 Mar 2023]
    TBA
  • Reconstruction and interpretation of poloidal fueling asymmetries revealed by line radiation diagnostics in DIII-D video
    George Wilkie, abstract
    [#s1583, 23 Feb 2023]
    Reliable prediction of the turbulent dynamics of the plasma edge remains one of the last frontiers of theory in support of designing a fusion power plant. Previous results from the LLAMA diagnostic on DIII-D revealed a striking asymmetry in the line radiation produced between the high- and low-field sides, and a strong dependence on the direction of the toroidal magnetic field. These results have been unexplained until now. Here we present the only first-principles reproduction of the DIII-D Lyman-alpha signal using kinetic models for the plasma and neutral species with synthetic diagnostics. It is found that the observed change in asymmetry is due to a difference in the primary recycling location which in turn is caused by changes in the plasma flow as the toroidal magnetic field changes direction. The asymmetry in the synthetic signal matches observation when accounting for collisional and drift physics, while turbulence is necessary to capture the order of magnitude.
  • Extended-MHD modeling of transients in spherical tokamaks and SPARC video
    Andreas Kleiner, abstract
    [#s1588, 16 Feb 2023]
    Edge-localized modes (ELMs) and disruptions are two transient phenomena that can cause serious damage to the vessel in reactor-scale tokamaks and thus need to be controlled or mitigated. This talk presents extended-magnetohydrodynamic (MHD) simulations with the goal of accurately predicting ELM stability thresholds in STs and inform the massive gas injection (MGI) layout for disruption mitigation in SPARC. ELMs are typically associated with macroscopic peeling-ballooning (PB) modes in the edge pedestal, which arise due to strong pressure and current density gradients. While in large aspect ratio devices these ideal-MHD modes are well understood, a long-standing problem has been the reliable modeling of such stability boundaries in some ST scenarios, particularly in NSTX. In simulations with the extended-MHD code M3D-C1, it is found that plasma resistivity can significantly alter macroscopic edge-stability in ELMing H-mode discharges in NSTX. These discharges are limited by resistive kink-peeling modes, while the two studied ELM-free scenarios appear limited by ideal ballooning modes. We will also present some extended-MHD analysis of PB stability in MAST and STAR, a preliminary ST-based power plant design. We show how these extended-MHD stability thresholds are incorporated into a higher-fidelity model to predict the pedestal structure in a wider range of tokamak scenarios. The second part of the talk focuses on extended-MHD simulations of disruption mitigation in SPARC via massive gas injection. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures in the wall.
  • A non-local magneto-curvature instability in differentially rotating plasmas. video
    Fatima Ebrahimi, abstract
    [#s1582, 09 Feb 2023]
    Global stability of differentially rotating systems in the presence of magnetic fields is examined. Unlike general global ideal pressure or current-driven instabilities, where an integrated linearized self-adjoint force operator is used to realize the free energy in the flowless MHD, due to the non-self-adjoint property of the MHD force operator a nontrivial modified energy principle with flows is required. Per Frieman & Rotenberg, I will first discuss the modified energy principle, in a differentially rotating system, and then present our results using both global eigenvalue analysis as well as initial value calculations. We find that only global models with spatially varying fields (both magnetic and rotational) can offer the richest mode spectrum, mainly as the result of resonances in the system. A new non-local mode, a magneto-curvature instability (Ebrahimi&Pharr ApJ 2022 https://doi.org/10.3847/1538-4357/ac892d) is obtained. This non-axisymmetric instability is triggered due to Alfven-continuum unstable modes in the presence of non-local effects of the global spatial curvature of flow shear and magnetic field. It will be shown that as the field strength is increased a transition from turbulence to a state dominated by global non-axisymmetric modes is obtained. The implication of this instability for accretion flows as well as laboratory plasmas will be discussed.
  • Modeling the Destabilization of a Network of Nonlinear Interactions during the ELM Onset with Triads of Harmonic Oscillators
    Julien Dominski, abstract
    [#s1555, 02 Feb 2023]
    A mechanism by which a nonlinear perturbation involving many active triads -- a Network of Nonlinear Interactions -- was observed to trigger edge localized modes (ELMs) below the peeling-ballooning (PB) limit [1]. The nonlinear stability of such a network of nonlinear interactions has been modeled with a network of harmonic oscillators coupled together via multiple triads [2]. This model network has been found to transit from a quiet regime of weak nonlinear fluctuations (triads near O-points) towards a regime of strong nonlinear fluctuations (triads near X-points). During this transition, the energy of the dominant waves is transferred to the sub-dominant waves through the strong nonlinear fluctuations, reminiscent of the ELM onset. The rapidity of this transition is found to be inversely proportional to the intensity of the nonlinear coupling between the modes. Moreover, when the nonlinear time scale of the fluctuations is comparable to the time scale of the wave oscillations, it is found that the system is chaotic and that the transition is the most abrupt. [1] J Dominski and A Diallo, Plasma Phys. Control. Fusion 62:095011 (2020) [2] J Dominski and A Diallo, Physics of Plasmas 28:092306 (2021)
  • Hybrid Simulations of Energetic Particle-driven n=1 Modes in Tokamak Plasmas video
    Guoyong Fu, Zhejiang University, abstract
    [#s1623, 27 Jan 2023]
    A systematic simulation study of energetic particle-driven n=1 mode in tokamak plasmas has been carried out using the global kinetic-MHD hybrid code M3D-K [1,2]. This work is focused on the interaction of energetic beam ions and n=1 mode with a monotonic safety factor q profile and q0 < 1. Linear simulations with energetic co-passing particles [1] show excitation of a low-frequency mode of fishbone type with the corresponding resonance of ωφ + ωθ = ω, where ωφ is the energetic ion toroidal transit frequency and ωθ is the energetic ion poloidal transit frequency. The simulated mode frequency is approximately proportional to the energetic ion injection energy and orbit width. The mode structure is similar to that of internal kink mode. These simulation results are similar to the analytic theory of Yu et al. [3]. Furthermore, linear simulations with energetic counter-passing particles [2] show that the instability is either a m/n = 1/1 energetic particle mode (EPM) or a m/n = 1/1 global Alfven eigenmode (GAE) depending on the value of central safety factor. The mode frequencies are close to the tip of Alfven continuum spectrum at the magnetic axis. The excited modes are radially localized near the magnetic axis well within the q = 1 surface. The main wave particle resonance is found to be ωφ + 2ωθ = ω. The nonlinear simulation results show that there is a long period of quasi-steady-state saturation phase with frequency chirping up after initial saturation. Correspondingly, the energetic particle distribution with low energies is flattened in the core of plasma. After this quasi-steady phase, the mode amplitude grows again with frequency jumps down to a low value corresponding to a new mode similar to the energetic co-passing particle-driven low frequency fishbone while the energetic particle distribution is flattened for higher energies in the core of plasma.
  • A Possible Mechanism for the Edge Transport Barrier Formation: The Story of Finite Larmor Radius Effects video
    Wei-li Lee, abstract
    [#s1560, 26 Jan 2023]
    The formation of the H-mode pedestal in magnetized plasmas remains a mystery nearly forty years after it was first experimentally observed [1]. In the ensuing years, it has been observed in nearly every machine experiment. Many different theories have also been proposed to explain its existence. Recently, Burrell [2] argued that sheared E x B flow was most likely responsible for the formation of the pedestal and the improved confinement. This is different from the claim by Lee and White [3] that the formation of the pedestal is the real reason for the improved confinement and, furthermore, its formation subsequently gives rise to the E x B flow. Moreover, they claimed that the H-mode is related to a delicate force balance between the ion pressure gradient and the gyroviscosity arising from the ion Finite Larmor Radius (FLR) effects [3]. The ion FLR effects are also found to be responsible for the creation of the radial electric field well for the H-mode [3,4]. This electric field, together again with the FLR effects associated with the E x B drift, can produce a poloidal current, which modifies the pressure balance. Furthermore, it has been shown that this delicate balance produces a force free field, ∇ x B = (4π /c) J∥, which is the real physics behind the H-mode [3]. A recent experimental paper by Zweben et al. [5] has shown that there is no noticeable change for the poloidal flow in the pedestal region just before the L-H transition, which differs from the claim by Burrell [2]. To answer this chicken-and-egg question, we propose to use the existing fully electromagnetic gyrokinetic codes, e.g., GTC [6] and GTS [7] to simulate the physics inside the separatrix with proper modifications. Specifically, the charge separation due to the finite Larmor radius effects in the regions with steep pressure gradient should be taken into account. Details will be discussed. 1. F. Wagner et al., Phys. Rev. Lett. 53, 1453 (1984) 2. K. H. Burrell, Phys. Plasmas 27, 060501 (2020) 3. W. W. Lee and R. B. White Phys. 26, 040701 Plasmas (2019) 4. W. W. Lee, Phys. Plasmas 23, 070705 (2016) 5. S. J. Zweben, A. Diallo, M. Lampert, T. Stoltzfus-Dueck, and S. Banerjee, Phys. Plasmas 28, 032304 (2021) 6. Z. Lin, T. S. Hahm, W. W. Lee, W. M. Tang and R. B. White, Science 281, 1835 (1998) 7. W. X. Wang, Z. Lin, W. M. Tang, W. W. Lee et al. Phys. Plasmas 13, 092505 (2006)
  • "Entity": general coordinate (QED)(GR)PIC code for astrophysical plasmas video
    Hayk Hakobyan, abstract
    [#s1567, 19 Jan 2023]
    Particle-in-cell has been a go-to approach for modeling plasmas in the environments of compact astrophysical objects for the last decade. Yet, there is no single publicly available code that includes all relevant radation-plasma coupling processes and is capable of modeling global systems. In this talk I will describe development of a new-generation PIC code for extreme astrophysical plasmas, Entity. The code is based on the Kokkos framework, which enables efficient implicit multi-architecture portability including GPUs. The code features algorithms for various radiation-plasma coupling processes, such as Compton scattering, production of electron-positron pairs and their annihilation. The code is designed in general coordinate system, defined by the metric functions; this enables the Entity to also efficiently tackle the global (full-system) models of the magnetospheres of compact objects, which require algorithms on non-cartesian (spherical, cubed sphere) non-uniform grids, and even full general relativity.
  • Electric field screening in pair discharges and generation of pulsar radio emission video
    Elizabeth Tolman, Institute for Advanced Study, abstract
    [#s1554, 15 Dec 2022]
    Pulsar radio emission may be generated in pair discharges which fill the pulsar magnetosphere with plasma as an accelerating electric field is screened by freshly created pairs. In this talk we present a simplified analytic theory for the screening of the electric field in these pair discharges and use it to estimate total radio luminosity and spectrum. The discharge has three stages. First, the electric field is screened for the first time and starts to oscillate. Next, a nonlinear phase occurs. In this phase, the amplitude of the electric field experiences strong damping because the field dramatically changes the momenta of newly created pairs. This strong damping ceases, and the system enters a final linear phase, when the electric field can no longer dramatically change pair momenta. Applied to pulsars, this theory may explain several aspects of radio emission, including the observed luminosity, 10^{28} erg s^{-1}, and the observed spectrum, ω^{-1.4+-1.0}.
  • Collisionless plasma sheaths with a grazing-incidence magnetic field (remote)
    Alessandro Geraldini, EPFL, abstract
    [#s1559, 08 Dec 2022]
    In magnetised plasma sheaths, such as the ones forming next to divertor or limiter targets in a fusion device, very strong electric fields arise on the length scales of the ion gyroradius and Debye length [1]. These reflect electrons and reduce the electron flux so that equal outflow of electrons and ions (ambipolarity) to the targets can be achieved globally. The length and time scales of the sheath are much shorter than those of the bulk plasma (e.g. the Scrape-Off Layer). It is therefore computationally much faster to simulate the plasma with models that average over the shorter scales and thus do not resolve the sheath (e.g. fluid, drift-kinetic, gyrokinetic), although this requires boundary conditions that are consistent with the presence of the sheath. For example, in a kinetic model the electrons must be reflected by the sheath electric field back into the plasma [2]. Hence, the phase space reflection-absorption boundary, or cutoff, of electrons fully specifies the electron boundary conditions. We present computations of the cutoff including the effect of finite electron gyro-orbits, which makes the cutoff parallel velocity a function of magnetic moment. Ions must be pre-accelerated into the sheath and satisfy a constraint known as the Bohm (unmagnetised) [3] or Chodura (magnetised) condition. We derive the kinetic generalisation of the Chodura condition for general magnetic field angles including the effect, which becomes prominent at shallow magnetic field angles, of gradients of the electrostatic potential and distribution functions tangential to the target (e.g. the ExB drifts from such gradients transport ions almost normal to the target, while parallel streaming only has a very small component normal to the target at shallow magnetic field angles). We calculate a critical small magnetic field angle [4] below which a monotonic sheath solution cannot be found, and find that this critical angle increases with electron gyroradius. We also present first preliminary 2-dimensional sheath electrostatic potential solutions including the spatial profile of small-amplitude fluctuations tangential to the target. This work assumes collisionless sheaths, and therefore must be generalised to be applied to colder and more neutral-rich systems such as a detached divertor. [1] R. Chodura, Phys. Fluids (1982). [2] Parker et al. J. Comput. Phys. (1993). [3] K.-U. Riemann, J. Phys. D: Appl. Phys. (1991). [4] R. Ewart, F. Parra, A. Geraldini, PPCF (2021).
  • Modeling turbulence saturation mechanisms for stellarator optimization video
    Benjamin Faber, U. Wisconsin, abstract
    [#s1569, 01 Dec 2022]
    A key challenge to achieving efficient fusion energy in toroidal devices is reducing the particle and energy losses from microinstability-driven turbulence. Efforts to optimize stellarators for reduced turbulence transport have previously been focused on reducing the linear instability drive mechanisms, however this may conflict with other optimization constraints such as macroscopic stability. A different approach to turbulence optimization may be formulated by attempting to reduce the nonlinearly saturated fluctuation amplitudes by increasing coupling to dissipation channels. While multiple dissipation mechanisms exist in plasma turbulence, this talk will focus on ability of fluctuations to nonlinearly couple to stable modes to provide an effective energy sink. This process is largely quantified by a resonant three-wave interaction lifetime between modes. In stellarators, the strength of the three-wave interaction lifetime is dependent on the geometric properties of magnetic field lines and can be estimated effectively from linear eigenmode calculations. This talk will demonstrate the differences in stable mode coupling mechanisms between classes of quasisymmetric geometries and the applicability of the model for stellarator optimization.
  • Integrated Simulation of Energetic Particles in Burning Plasmas video
    Zhihong Lin, UC Irvine, abstract
    [#s1578, 15 Nov 2022]
    Energetic particle transport in burning plasmas depends on cross-scale interactions between macroscopic MHD mode, mesoscale Alfven eigenmodes, and microturbulence, which requires global integrated simulations incorporating multiple physical processes. In this talk, I will first highlight gyrokinetic toroidal code (GTC) simulation finding regulation of reversed shear Alfven eigenmodes (RSAE) by microturbulence, leading to excellent agreement, for the first time, of gyrokinetic simulation results with experimental measurements of RSAE amplitude and mode structure in the DIII-D tokamak [PRL 128, 185001 (2022)]. I will then present results from integrated simulations of energetic particle transport in the ITER operational scenarios carried out by several energetic particle codes in SciDAC and ITPA collaborations. Finally, I will discuss progress on GTC global simulations with kinetic electrons in 3D toroidal geometry including RMP tokamaks and stellarators.
  • Effects of elongation and triangularity on blob dynamics in the SOL video
    Tess Bernard, General Atomics, abstract
    [#s1568, 07 Nov 2022]
    Plasma blobs are coherent turbulent structures of elevated plasma pressure in the scrape-off layer (SOL) of tokamaks that influence radial transport. Previous work has considered the effect of X-point geometry on radial blob velocities [1]. Now a semi-analytic blob velocity scaling has been derived that includes the effects of elongation and triangularity. It predicts that increasing the value of either shaping parameter results in slower radial blob velocities. Using the Gkeyll code, gyrokinetic seeded blob simulations have been carried out assuming DIII-D SOL parameters and inner-wall limited (IWL) geometry, which confirmed scaling predictions. The effect of atomic neutral interactions was also explored. Finally, full flux simulations of positive and negative triangularity IWL DIII-D discharges [2,3] are presented. The negative triangularity case has 50% more blobs than the positive case, and they are slightly faster on average, consistent with the scaling prediction.
  • Problems with Stellarators (remote) video
    Roscoe White, abstract
    [#s1552, 03 Nov 2022]
    Stellarators have an advantage over tokamaks as fusion reactors in that they are disruption free and steady state, whereas tokomaks must be pulsed in order to produce the poloidal confining field. In a stellarator the poloidal field is produced by making the field coils not simply vertical, but with an inclination which is periodically modulated around the torus. But this modulation has also the effect that the magnitude of the toroidal field is not constant, but has a toroidal periodicity. This modulation of the toroidal field has however three deleterious effects. Particle resonances in a plasma are locations where passing particle orbits return to an initial point, they repeat the same motion indefinitely, n periods toroidally and m periods poloidally. The location of a resonance is a function of particle energy and pitch. At very low energy they are located at values of the field line helicity equal to m/n. If the toroidal period of a resonance matches the toroidal period of B, then 1. Unperturbed passing particle orbits form islands, with width increasing with energy, and these islands to some degree impair confinement. 2. Local wells in B are formed along the resonances, producing ripple trapping loss of particles, making early loss of fusion alpha particles a problem for reactor walls. 3. In addition, because of the toroidal dependence of B, Alfven modes do not produce local resonance islands, with the effect of the mode localized to very near the resonance surface. Instead, a broad domain of chaotic orbits is produced, with particle diffusion in all energies large for small mode amplitude.
  • Fundamental Physics Basis for Transport Barriers without Velocity Shear (remote), video
    Mike Kotschenreuther, University of Texas, abstract, slides
    [#s1541, 27 Oct 2022]
    Transport Barriers (TBs) are crucial to magnetic fusion, in the form of edge pedestals (and also Internal Transport Barriers (ITB)). The dominant instability of confined plasmas, the coupled Ion Temperature Gradient and Trapped Electron Mode, must be suppressed for TBs to exist. In a complete departure from previous theoretical analysis, we use statistical mechanical concepts together with extensive gyrokinetic simulations to show that such thermalizing instabilities cannot access the enormous free energy of their gradients due to a fundamental dynamical constraint- the fluctuation induced charge flux must vanish. Together with basic statistical mechanics, the dynamical pathway to instability can be removed by this constraint. Velocity shear, often thought of as crucial for TBs, is thereby rendered unnecessary. In fact, even for well-developed pedestals in present devices, simulations show that this constraint can be even more responsible for sustaining the barrier than velocity shear. On future devices with considerably lower velocity shear, this will undoubtedly have to be true. The constraint becomes restrictive when one plasma species becomes nearly adiabatic ( = Maxwellian response) due to rapid phase space averaging (as in classical statistical mechanics). This averaging depends sensitively upon the magnetic geometry. When phase space averaging rates are faster than the mode dynamics, and when there are also substantial density gradients, the dynamical constraint forces low growth rates independently of the amount of free energy. Nonlinear simulations show that heat fluxes are also commensurately small (reduced by orders of magnitude).
  • A Gyrokinetic Moment-based Method to model the boundary region of fusion devices using advanced collision operators video
    Baptiste Jimmy Frei, EPFL, abstract
    [#s1558, 24 Oct 2022]
    We present a new first-principle model that allows for the proper simulation of the plasma boundary of fusion devices, which encompasses the edge and the scrape-off layer regions [1]. Developed onto a set of fluid-like equations using a Hermite-Laguerre polynomial decomposition of the distribution function that retains the gyrokinetic (GK) Coulomb collision operator [2], the gyro-moment (GM) model offers an ideal analytical and numerical framework to describe the wide range of plasma parameters found in the boundary region. Indeed, the GM model contains the core GK model and the fluid and gyrofluid models used for scrape-off layer simulations as particular limits. We demonstrate that the GM approach can correctly retrieve the properties of microinstabilities that develop at low plasma collisionality, strongly sensitive to kinetic features, in perfect agreement with the GK continuum GENE code. At the same time, we show that the GM model correctly retrieves the fluid limit at high collisionality. A hierarchy of collision operator models with various physics fidelity is developed and numerically implemented using the same GM approach [3,4,5]. This allows us to perform a comparison between collision operator models, revealing large deviations (compared to the Coulomb operator) in linear growth rates, collisional zonal flow damping, and turbulent transport levels. Furthermore, we prove that the GM approach is numerically efficient from the low-collisionality banana regime in H-mode pedestals to the high-collisionality regime of the scrape-off layer. [1] Frei B. J. et al., JPP 82 (2020) [2] Jorge R. et al., JPP 85 (2019) [3] Frei B. J. et al., JPP 87 (2021) [4] Frei B. J. et al., JPP 88 (2022) [5] Frei B. J. et al., PoP 29 (2022)
  • On the Kinetics of Stars and Bars, or The Galactic Tokamak video
    Chris Hamilton, Institute for Advanced Study, abstract
    [#s1536, 13 Oct 2022]
    Many galaxies, including our own galaxy the Milky Way, have a ‘bar’ structure at their center— an elongated collection of millions of stars, that gradually rotates as if it were a solid body. Galaxies are also embedded in massive dark matter haloes. When the rate at which the bar rotates resonates with a dark matter particle’s orbital frequency, the dark matter can suck angular momentum out of the bar, causing it to slow down. Previous theories of the bar-halo interaction calculated this ‘dynamical friction’ on the bar in a manner directly analogous to Landau’s calculation of the collisionless damping of an electric wave (and subsequently to O’Neil’s nonlinear generalisation thereof). This is no coincidence — bar-halo interactions are just one of a plethora of gravitational dynamics problems that have direct plasma-kinetic analogues. In this talk I will introduce the astrophysical context of galactic bars and their host dark matter haloes, and describe some of the aforementioned classic studies of the bar-halo interaction. However, those studies routinely ignored the fact that dark matter particles also experience random ‘diffusive’ kicks from other passing dark matter clumps, gas clouds, and so on. I will describe recent work done in collaboration with Princeton plasma theorists on quantifying the impact of diffusion on bar-halo friction, a problem which turns out to be mathematically identical to that of understanding particle energization in tokamaks. More broadly, I will argue that galactic dynamics has, over the last several decades, largely failed to learn from its more well-developed cousin and therefore has a lot of catching up to do. But I will also argue that in return, we stellar dynamicists can provide you plasma theorists with fresh contexts in which to ply your trade, and an opportunity to work on something both fantastically beautiful and totally useless.
  • Department (re-)introductions - Part II
    Theory Department
    [#s1550, 06 Oct 2022]
  • Toward optimizing the performance of fusion reactors with transport in the loop (remote)
    Noah Mandell, MIT, abstract, slides
    [#s1556, 04 Oct 2022]
    As we approach the breakeven era of fusion, optimizing reactors to make them more efficient and less expensive will be critical to the wide-scale adoption of fusion as a commercial energy source. The main challenge is to achieve high steady-state pressures in the core of the reactor to reach self-sustaining fusion conditions. At the same time, the boundary plasma must be kept sufficiently cool so that the plasma exhausted from the hot core is not dangerous to the device walls. Turbulence is the main source of heat transport from the core to the boundary, which makes understanding how to optimize the reactor design for turbulent transport a key to solving the competing (but coupled) core and boundary challenges. In this talk, I will present a vision for tackling this challenging whole-device-modeling problem in a scalable way. The approach consists of four main modules: (1) fast-but-accurate core turbulence modeling with the GPU-native GX delta-f gyrokinetic code, which leverages pseudo-spectral methods in both configuration (Fourier) and velocity (Hermite-Laguerre) space; (2) multi-scale modeling for predicting and evolving core profiles using a macro-scale transport solver (Trinity) coupled to many radially-local GX micro-turbulence calculations in parallel, leveraging the scale separation between turbulence and transport at reactor scale; (3) kinetic boundary turbulence modeling with the Gkeyll code, a full-f electromagnetic gyrokinetic model for the edge and scrape-off layer; and (4) transport optimization of fusion reactor designs by using (1-3) as a massively-parallel whole-device model inside the optimization loop. For each of these modules I will highlight preliminary results, ongoing work, and future steps.
  • Department (re-)introductions - Part I
    Theory Department
    [#s1549, 29 Sep 2022]
  • Microturbulence in edge of HL-2A tokamak plasma and its interaction with a magnetic island (remote)
    Jingchun Li, SUSTech, abstract, slides
    [#s1553, 27 Sep 2022]
    Transport Barriers (TBs) are crucial to magnetic fusion, in the form of edge pedestals (and also Internal Transport Barriers (ITB)). The dominant instability of confined plasmas, the coupled Ion Temperature Gradient and Trapped Electron Mode, must be suppressed for TBs to exist. In a complete departure from previous theoretical analysis, we use statistical mechanical concepts together with extensive gyrokinetic simulations to show that such thermalizing instabilities cannot access the enormous free energy of their gradients due to a fundamental dynamical constraint- the fluctuation induced charge flux must vanish. Together with basic statistical mechanics, the dynamical pathway to instability can be removed by this constraint. Velocity shear, often thought of as crucial for TBs, is thereby rendered unnecessary. In fact, even for well-developed pedestals in present devices, simulations show that this constraint can be even more responsible for sustaining the barrier than velocity shear. On future devices with considerably lower velocity shear, this will undoubtedly have to be true. The constraint becomes restrictive when one plasma species becomes nearly adiabatic ( = Maxwellian response) due to rapid phase space averaging (as in classical statistical mechanics). This averaging depends sensitively upon the magnetic geometry. When phase space averaging rates are faster than the mode dynamics, and when there are also substantial density gradients, the dynamical constraint forces low growth rates independently of the amount of free energy. Nonlinear simulations show that heat fluxes are also commensurately small (reduced by orders of magnitude).
  • Intrinsic rotation driven by the radial variation of phase velocity and turbulence intensity in tokamaks (remote)
    Denis St-Onge, University of Oxford
    [#s1551, 21 Sep 2022]
  • Preferred Magnetic Axes For Optimal Quasi-Axisymmetry video
    Wrick Sengupta, PPPL, abstract
    [#s1516, 11 Aug 2022]
    Recently, enormous progress has been made in obtaining quasisymmetry (QS) of outstanding precision through numerical optimization. Significant analytical progress has also been made possible thanks to the asymptotic expansions near the magnetic axis (NAE). A critical factor in realizing good QS through the second-order of the NAE is the choice of the magnetic axis. However, because of the complexity of the second-order NAE equations, analytical characterization of these preferred magnetic axes for optimal QS has not been possible so far. In this talk, we shall attempt to answer this question for quasi-axisymmetric (QA) systems. We show that the magnetic axis is well described for small rotational transforms by the same equations that govern Euler-Kirchhoff elastic rod centerlines (Langer and Singer, SIAM review 1996, Pfefferlé et al. PoP 2018). Surprisingly, the connection to these equations can only be made partially within the NAE framework and requires several concepts from the soliton theory. We shall present analytical and numerical evidence supporting our insights for a broad range of QA stellarators.
  • Energetic particle transport in 3D magnetic fields: Loss mechanisms and optimization strategies
    Elizabeth Paul, Princeton University, abstract
    [#s1493, 28 Jul 2022]
    Collisionless physics primarily determines the transport of fusion-born alpha particles in 3D equilibria. Several transport mechanisms have been implicated in stellarator configurations, including stochastic diffusion due to class transitions, ripple trapping, and banana drift-convective orbits. Given the guiding center dynamics in a set of six quasihelical and quasiaxisymmetric equilibria, we perform a classification of trapping states and transport mechanisms. In addition to banana drift convection and ripple transport, diffusive banana tip motion associated with the non-conservation of the parallel adiabatic invariant is substantial among prompt losses, especially in equilibria close to quasiaxisymmetry. Furthermore, many lost trajectories undergo transitions between trapping classes on longer time scales, either with periodic or irregular behavior. We discuss possible optimization strategies for each of the relevant transport mechanisms and perform a comparison between classified guiding center losses and recently-developed metrics for banana drift convection transport. Equilibrium characteristics responsible for distinctions in transport are discussed. Quasihelical configurations are found to have natural protection against both ripple trapping and diffusive banana tip motion leading to a reduction in prompt losses.
  • Portable, Structure Preserving Kinetic Methods in PETSc video
    Mark Adams, Lawrence-Berkeley National Laboratory, abstract
    [#s1515, 21 Jul 2022]
    I will introduce our integrated program for both developing structure-preserving (SP) methods for kinetic applications and deploying these methods as high-performance tools in PETSc (Portable, Extensible, Toolkit for Scientific computing). The metriplectic formalism used to develop these methods is introduced along with a strictly conservative, monotonic entropy particle-in-cell (PIC) application, PETSc-PIC. We briefly discuss several methods developed for this code: a mixed Poisson solver for a C^0 electric field; strictly conservative mapping between particle and finite element bases; monotonic entropy time integrators for collisions; and two new particle based Landau collision operators. A mature, high-order accurate, finite element based Landau collision operator with adaptive mesh refinement, that has been optimized for accelerator architectures using the portable Kokkos programming language, is presented with results using the NVIDIA A100 and AMD MI250X architectures. We present verification studies using plasma resistivity and compare with Spitzer resistivity. New batch GPU linear solvers have been developed for this work. This work is integrated into the PETSc “solver” framework to provide a fully GPU time advance of the Landau collision operator. We show that collision time advance with many species is practical for realistic models of tokamaks on today’s large-scale computers in cylindrical coordinates and that fully 3V models should be feasible in the near future.
  • Ultra long turbulent eddies, magnetic topology, and the triggering of ITBs, (video)
    Justin Ball, EPFL, abstract
    [#s1518, 15 Jul 2022]
    Local nonlinear gyrokinetic simulations of tokamak plasmas demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear is very small. By accurately modeling different field line topologies (e.g. low-order rational, almost rational, or irrational), we show that the parallel self-interaction of such "ultra long" eddies can significantly reduce heat transport. This reveals novel strategies to improve confinement, constitutes experimentally testable predictions, and illuminates past observations of internal transport barriers.
  • Linear stability of ultra-high-beta equilibria video
    Rahul Gaur, University of Maryland, abstract
    [#s1476, 30 Jun 2022]
    The power density of tokamaks scales with the plasma beta as beta^2 which makes high-beta operation an attractive choice for future high-power-density tokamak devices [Menard et al Nucl. Fusion 22]. Ultra-high-beta (beta ~ 1) configurations have previously [Hsu. et al. PoP 96] been explored at the level of asymptotic MHD equilibria by solving the Grad-Shafranov equation in the limit (delta_Hsu)^2 ~ epsilon/(beta q^2) << 1. We extend this by obtaining exact global equilibria numerically. However, various instabilities may limit the utility of such equilibria. To that end, we present an infinite-n ideal-ballooning and linear gyrokinetic analysis of ultra-high-beta (beta~1) equilibria for tokamaks. In the first part, we examine ideal ballooning stability. We find that alpha_MHD ~ 1/(delta_Hsu)^2 >> 1 is large enough to make them "second-stable" to the ideal ballooning mode. Upon ensuring ideal ballooning stability, we examine their stability to the two major sources of electrostatic turbulence: ITG and TEM, using the initial value code GS2. To understand the trend with a changing beta, we compare these equilibria with an intermediate-beta (beta~0.1) and a low-beta (beta~0.01) equilibrium at two different radial locations: the inner core (Normalized radius rho = 0.5) and the outer core (rho = 0.8) for two different triangularities: delta = 0.4 and delta = -0.4. We find that the ultra-high-beta equilibria are stable to both the ITG and TEM over a wide range of gradient scale lengths (R/L_T and R/L_n). Next, we perform a linear electromagnetic study of all the nominal local equilibria to explore the possible effects of Kinetic Ballooning Modes (KBMs). We find that all the high-beta equilibria become more unstable than their low-beta counterparts in the inner core but turn out to be much more stable than both the low or intermediate beta equilibria in the outer core. We also find that the negative-triangularity high-beta equilibria do not show any signs of KBMs. Using a full gyrokinetic code for linear electromagnetic studies at k_perp rho << 1 can be relatively expensive. Therefore, as an alternative, we numerically solve the KBM equations of Tang et al. in the limit omega_bi < omega < omega_be as a reduced KBM model and compare the results with GS2.
  • Multi-channel validation of a quasilinear gyrokinetic transport model in flux-driven integrated modelling (remote) video
    Jonathan Citrin, DIFFER, abstract
    [#s1499, 23 Jun 2022]
    Recent significant progress has been made in the application of first-principle-based reduced turbulent transport models within integrated modelling. We focus on recent developments of the quasilinear gyrokinetic transport model QuaLiKiz. Coupling to tokamak integrated modelling frameworks allows flux-driven core transport modelling for heat, particle, and momentum channels, with JET discharge timescales simulated in ∼100CPUh. We sketch the basis of the QuaLiKiz transport model and its validity in comparison to nonlinear simulations. We then review validation of the model against experimental measurements through flux-driven integrated modelling simulations including multiple transport channels and multiple ions. This capability enables the physics interpretation of present-day experiments – where specific examples of new insight into transport mechanisms will be provided – as well as extrapolation to future machine perf6ormance. We also present a parallel approach for fast integrated modelling based on neural network regression of an extensive QuaLiKiz run database. The neural network transport model is ×10^6 faster than QuaLiKiz itself, opening up new possibilities for first-principle-based scenario optimization and control- oriented applications.
  • Quantitative measurements of ion orbit loss from gyrokinetic simulations video
    Hongxuan Zhu, PPPL, abstract
    [#s1480, 16 Jun 2022]
    Ion-orbit loss is considered important to the radial electric fields Er of tokamak edge plasmas. In neoclassical equilibria, collisions can scatter ions onto the loss orbits and generate a steady-state radial current, which may drive the edge Er away from the confined-region neoclassical value without orbit-loss. To quantitatively measure this effect, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. The validity of the diagnostic is demonstrated by studying the collisional relaxation of Er in the core plasmas. Then, the ion orbit-loss effect is numerically measured in the edge plasmas in the DIII-D geometry. It is found that the effect of the collisional ion orbit loss is more significant for an L-mode plasma compared to an H-mode plasma.
  • A survey of fusion reactor systems codes, and how you can help, video
    Charles Swanson, PPPL, abstract, slides
    [#s1512, 09 Jun 2022]
    In this seminar, I present a survey of the most commonly used fusion reactor systems codes in the literature, and suggest promising avenues toward increased flexibility, fidelity, and utility. Systems codes are the "jack of all trades, master of none" of the fusion reactor modeling world. They aim to model the entire facility at a low fidelity, rather than any one phenomenon or system at high fidelity. They integrate models from plasma physics, engineering, and economics. Systems codes range in complexity from simple spreadsheets to codes which implement numerical optimizers with user-configurable constraints and iteration variables. In this seminar I argue that future systems codes should be modular, flexible, and extensible, allowing the user to implement different workflows and mix-and-match models in their analyses. I identify holes in the landscape of low- and medium-fidelity models, especially those applicable to Stellarators.
  • Theory and modelling of the Centrifugal Mirror Fusion Experiment
    Ian Abel, University of Maryland, abstract, slides
    [#s1388, 10 Feb 2022]
    Due to the success of the Maryland Centrifugal Experiment (MCX) [R. F. Ellis et. al. Phys. Plasmas 8, 2057 (2000)] and initial theoretical analyses, the Centrifugal Mirror concept is being further explored by the construction of the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. This prompts a deeper inquiry into the underlying confinement and stability properties of centrifugal mirrors as a class of devices. We will outline the key advances of this experiment over prior rotating mirrors , and give an updated physics basis for our modelling of the experiment and a future reactor. We will highlight a new semi-analytical calculation for the equilibrium of such a system derived in the asymptotic limit of a rapidly-rotating plasma. This solution is used to shed light on the ultimate limits of rotation in such plasmas. We also use this magnetic equilibrium to compute approximate end loss rates and electrical properties of the plasma. Corrections to the parallel electric field are computed to ensure all losses are ambipolar. This provides a self-consisent basis for a 0D ``systems'' model of a centrifugal mirror machine. As examples we provide design points corresponding to CMFX and an upgraded device capable of achieving breakeven.
  • Geometric integration of Hamiltonian systems on exact symplectic manifolds
    Joshua Burby, LANL, abstract, slides
    [#s1393, 03 Feb 2022]
    Non-dissipative (i.e. Hamiltonian) dynamical systems freeze flux in phase space, just as highly-conductive plasma flows freeze magnetic flux. A time integrator for a non-dissipative system is symplectic when it freezes flux exactly. Symplectic integration is routine in canonical coordinates, where the flux tensor takes the simplest possible form. Much less is understood about symplectic integration in the general non-canonical case, which occurs more frequently in practice. In this talk, I will present a general approach to structure-preserving integration of noncanonical Hamiltonian systems on exact symplectic manifolds. First, the original non-canonical Hamiltonian system is embedded in a larger (essentially) canonical system as a slow manifold. Then a canonical symplectic integrator for the larger system is identified that has approximately the same slow manifold. Provided initial conditions are selected near the slow manifold, the integrator provides a good approximation of the original system. There would be a problem with this approach if the discrete-time slow manifold happened to have any normal instabilities; such instabilities would carry discrete trajectories away from the slow manifold, and the good approximation properties would break down. I will explain how this potential problem is avoided using a newly-developed theory of nearly-periodic maps. By constraining the large system's integrator to be a non-resonant nearly-periodic map, existence of a discrete-time adiabatic invariant is guaranteed. Long-time normal stability of the slow manifold then follows from a Lyapunov-type argument.
  • Global-local gyrokinetic simulations of turbulence in tokamak plasmas using STELLA, video
    Denis St-Onge, University of Oxford, abstract, slides
    [#s1395, 13 Jan 2022]
    We develop a novel approach to gyrokinetics where multiple flux-tube simulations are coupled together in a way that consistently incorporates global profile variation while allowing the use of Fourier basis functions. By doing so, the need for Dirichlet boundary conditions typically employed in global gyrokinetic simulation, where fluctuations are nullified at the simulation boundaries, is obviated. This results in a smooth convergence to the local periodic limit as rho_* -> 0. In addition, our scale-separated approach allows the use of transport-averaged sources and sinks, offering a more physically motivated alternative to the standard sources based on Krook-type operators. Having implemented this approach in the flux-tube code stella, we study the role of transport barriers and avalanche formation in the transition region between the quiescent core and the turbulent pedestal, as well as the efficacy of intrinsic momentum generation by radial profile variation. Finally, we show that near-marginal plasmas can exhibit a radially localized Dimits shift, where strong coherent zonal flows give way to flows which are more turbulent and smaller scale.
  • Three-Dimensional Electron Temperature Gradient Turbulence in the Tokamak Pedestal,video
    Jason Parisi, PPPL, abstract, slides
    [#s1394, 06 Jan 2022]
    Electron temperature gradient (ETG) turbulence in the tokamak edge pedestal has a rich three-dimensional structure, particularly at ion-gyroradius-scales. Nonlinear multiscale gyrokinetic simulations of Joint European Torus (JET) pedestals reveal that ETG pedestal turbulence is highly inhomogeneous in the direction parallel to the magnetic field. Its parallel distribution is determined by the magnetic field geometry, with magnetic drift and finite Larmor radius effects being particularly important. Simulations must run sufficiently long for ion-gyroradius-scale ETG turbulence to saturate and interact with electron-gyroradius-scale ETG turbulence. Simulations without ion-gyroradius-scale ETG turbulence produce at least 65 % higher heat transport, indicating the transport-relevance of ion-gyroradius-scale ETG for multiscale ETG-ETG interactions.
  • TBA
    Fatima Ebrahimi, abstract
    [#s1564, 05 Jan 2022]
    TBA
  • TBA
    Fatima Ebrahimi, abstract
    [#s1581, 05 Jan 2022]
    TBA
  • Plasma simulation with exact Gauss law and without volume mesh, video
    Torsten Keßler, Saarland University, abstract
    [#s1389, 16 Dec 2021]
    I present a solver for the Vlasov-Poisson system that computes the electric field from boundary values of the electric potential only. First, I review the formulation of boundary value problems by integral operators and their discretization by Boundary Element Methods (BEM) with a particular focus on plasma dynamics. The reduction of dimension drastically reduces the number of unknowns while yielding accurate values for the electric field near the boundary. Furthermore, the method exactly preserves the particles as the source of the electric field. I demonstrate the power of the BEM approach with 3+3 dimensional numerical examples such as the formation of sheaths and a particle accelerator with complex geometry and mixed boundary values.
  • Structure-preserving marker-particle discretization of the Landau collision operator, video
    Eero Hirvijoki, Aalto University, abstract, slides
    [#s1387, 09 Dec 2021]
    Particle methods for Fokker-Planck collision operators typically rely on the stochastic approach. Diffusion is interpreted as random kicks and particle motion is described by a stochastic differential equation. This is not the only possibility, though. Diffusion can also be interpreted as a compressible vector field resulting from a gradient of an entropy functional, and the particles pushed along this self-consistently evolving vector field deterministically. This observation puts the collisional motion to an equal footing with the Hamiltonian contribution that moves particles along an incompressible vector field driven by the gradient of the total energy functional. In this seminar, I describe how these ideas can be exploited to discretize the non-linear Landau collision operator, utilizing a collection of marker particles of arbitrary weights while preserving positivity, discrete-time momentum and energy conservation, and monotonic entropy evolution.
  • PTOLEMY: Relic Neutrino Detection with Possible Insights for Plasma Physics, video
    Christopher Tully, Princeton University, abstract
    [#s1373, 21 Oct 2021]
    PTOLEMY is an experimental realization of a detection concept for neutrinos created in the first second following the Big Bang. New initiatives beginning with a setup in the D-Site Test Cell basement in 2011 have yielded transformational developments in collaboration with PPPL in material science [1]. Today, we report on running the PTOLEMY transverse drift filter “in reverse” [2]. This new technique for particle acceleration has remarkable possibilities for a “circular saw” approach to e-beam lithography and a novel tool for high efficiency, charged ion injection through transverse drift into strong magnetic fields. A preliminary study of D+ injection into an NSTX-U-like geometry is explored. [1] High hydrogen coverage on graphene via low temperature plasma with applied magnetic field, F. Zhao, Y. Raitses, X. Yang, A. Tan and C.G. Tully, Carbon 177 (2021) 244-251 [2] https://arxiv.org/abs/2108.10388
  • DREAM: a fluid-kinetic framework for tokamak disruption runaway electron simulations, video
    Mathias Hoppe, Chalmers University, abstract, slides
    [#s1353, 14 Oct 2021]
    Runaway electrons generated during a tokamak disruption pose a severe threat to future reactor-scale devices. Due to the exponential sensitivity of the runaway generation rate to the plasma current, robust avoidance and mitigation schemes cannot be fully validated in the medium-size tokamaks today, making comprehensive and validated runaway electron generation models essential for the development of such schemes. In this contribution we present the Disruption Runaway Electron Analysis Model (DREAM), a new simulation tool specifically designed to study the generation of runaway electrons during tokamak disruptions. The tool combines 1D fluid models for the background plasma (electric field, temperature, poloidal flux, ion charge states) with either fluid or kinetic models for the electrons in tokamak geometry. To enable accurate and efficient simulations of the whole disruption, electrons are separated into three sub-populations based on their energies, allowing different models to be used for thermal, superthermal, and relativistic electrons simultaneously. Notably, the thermal and runaway electrons can be treated using conventional fluid models, while the superthermal electrons are evolved using a reduced kinetic equation, providing precise accounting of the transient---and thus inherently kinetic---hot-tail runaway generation mechanism. In addition to the novel treatment of electrons, DREAM incorporates a number of physical mechanisms which have never before been brought together in a complete, self-consistent disruption simulation, including radial transport of heat and electrons, dynamic evolution of ion charge states, collisions with partially ionized atoms, the effect of passive conducting structures on the electric field, and hyperresistivity. The first studies conducted with DREAM indicate that fast electron radial transport may provide a path to effective runaway electron avoidance in ITER.
  • A discussion around the magnetic differential equation
    Lise-Marie Imbert-Gérard, University of Arizona, abstract
    [#s1351, 07 Oct 2021]
    This talk will focus on mathematical modeling aspects in the context of stellarator design. We will discuss several aspects of the magnetic differential equation, both on a toroidal flux surface and in a 3D volume filled with such surfaces. We will in particular describe some properties related to the periodic setting and the rotational transform, as well as their relation to the existence of solutions. We will also describe how the general Fourier solution's singularities relate to the questions of existence and uniqueness of a solution. We will then open a discussion regarding the choice of assumptions to derive equations of interest to study stellarator design.
  • Explaining the lack of power degradation of the energy confinement in the wide pedestal quiescent H-mode via transport modelling, video
    Saeid Houshmandyar, Institute for Fusion Studies, U. Texas at Austin, abstract
    [#s1363, 27 Sep 2021]
    Wide pedestal quiescent H-mode (WPQH) is an attractive scenario for the future burning plasmas as they operate without ELMs. Unlike the conventional H-modes, WPQH does not show power degradation of the energy (tau_E) confinement, as the heating power is increased. In these experiments, the neutral beam heating power (PNBI) was varied between 3.7 and 5.5 MW, while the net torque from the beams was kept nearly zero. As the PNBI was increased, reduced transport calculated by TRANSP, as well as the increased core ExB shear rate are observed; all suggesting the formation of an ion internal transport barrier (i-ITB) and increased stored energy in the core. In this work, a local quasilinear turbulent transport model enabled by Trapped Gyro Landau Fluid (TGLF), was used to predict the ITB and its stability analysis. tau_E calculated from the new constructed equilibria with the modelled profiles show insensitivity to the increased PNBI; these modelled profiles use TGYRO transport solver with matched energy fluxes between TGLF and TRANSP. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by ExB shear, the Ti/Te ratio and the Shafranov shift. Detailed analysis will be presented in this talk.
  • Generalized delta-f Particle-In-Cell method for arbitrary boundary conditions and forces in the system, video
    Min-Gu Yoo, PPPL, abstract
    [#s1358, 23 Sep 2021]
    The Particle-In-Cell (PIC) method using delta-f markers is widely adopted for simulating various plasma turbulence because delta-f markers can capture small-amplitude perturbations well with low statistical noise. Although the delta-f method can represent arbitrary phenomena without loss of generality, it is typically used to simulate weak turbulence within simple periodic boundaries, narrowing the usefulness of delta-f simulations. This is because the delta-f method has been considered unsuitable for implementing more general physical boundaries, such as absorbing walls or arbitrary sources. Besides, the typical delta-f simulation does not work well when strong forces are applied in the system. In this study, however, we show that the delta-f method can handle all situations in the same way as the full-f method, even for arbitrary boundary conditions and forces in the system. The delta-f method can correctly solve the Vlasov equation in arbitrary situations if the delta-f markers completely fill the entire phase-space of the simulation domain without any vacancy. We propose an efficient and practical way to generate new delta-f markers from the spatial physical boundary and the virtual velocity boundary to fill the phase-space of the simulation domain within some statistical noise. This generalized delta-f PIC method is verified in various 1D-1V physical problems such as Landau damping, two-stream instability, and plasma expansion to the vacuum. The method can be easily extended to multi-dimensional problems. For example, the global gyrokinetic simulation code GTS can successfully simulate thermal quench phenomena with a much lower computational cost by using the novel generalized delta-f PIC method.
  • Low recycling regime for burning plasma, video
    Leonid Zakharov, LiWFusion, abstract
    [#s1352, 16 Sep 2021]
    After two decades since two unlucky attempts to achieve the breakeven Q_DT=1 in tokamaks, this minimum fusion milestone remains unachievable. Although the tokamak program has developed all means for progress, the most critical reserve, i.e. suppression of recycling as the powerful mechanism of plasma cooling, is not yet utilized. For 63 years recycling is close to 100% in tokamaks. This results in excessive heating power, in limited confinement, unsuitable for burning plasma, and in unsolvable problems of power extraction. High recycling makes the entire regime complicated, while plasma unpredictable and disruptive. In 2012 the technology of Continuously Flowing Liquid Lithium (24/78-FLlLi) was invented for reduction recycling to 50% by pumping the escaping plasma particles by a creeping lithium layer. The suppression of the edge plasma cooling leads to a new plasma regime with an order of magnitude better confinement, much simpler plasma physics, insensitive to thermal conduction, and to reducing presently complicated plasma surface interactions to interaction a flow of energetic particles with Li. Simpler plasma control gives hopes to disruption avoidance. The recent assessment of low recycling regime for JET-like parameters predicts Q_DT > 5 in burning plasma at NBI power P_NBI=4 (!) MW. In the author's view, such a regime would be appropriate for the third DT campaign on JET as well as for restoring credibility of fusion.
  • Moment preserving interpolation in velocity-space and resampling in the global total-f gyrokinetic particle-in-cell code XGC, video
    Albert Mollén, PPPL, abstract, slides
    [#s1337, 26 Aug 2021]
    In addition to particle-in-cell methods, the XGC code evaluates dissipative operations such as collisions and heat sources/sinks on a 5-D grid at each fixed time step. This requires a mapping between the marker particles and the 5-D phase-space grid. If the error in the particle energy conservation is undesirably large in this mapping it can cause non-negligible numerical heating in a steep edge pedestal. Here we discuss a novel mapping technique in velocity-space, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. The new interpolation method relies on a particle resampling technique which is used to create, annihilate or redistribute particles in configuration space while preserving a desired number of moments. We will also discuss details of the resampling. [A. Mollén et al. J. Plasma Phys. 87 905870229 (2021)]
  • "Topological waves in magnetized cold plasmas", video
    Y. Fu, Princeton University , abstract, slides
    [#s1318, 15 Jul 2021]
    In the past decade, topological phases of electronic and photonic systems have become a rapidly emerging field of research, which deepened the understanding of the states of matters. One of the essential physical consequences of topological phases is the bulk-edge correspondence, which states that topologically protected edge modes will occur at the interface between topologically different matters. Recently studies have shown that these topological ideas can also be applied to continuous media, such as neutral fluids and plasmas. However, for magnetized plasmas, a comprehensive picture of topological phases and topological phase transitions is still missing. In the present study, we systematically map out all the topological phases and establish the bulk-edge correspondence in cold magnetized plasmas. We find that there exist ten topological phases in the parameter space of density, magnetic field, and parallel wavenumber, and this provides a new classification scheme for magnetized plasmas. A sufficient and necessary condition for the existence of the topologically protected edge modes is given and verified by numerical solutions. We find that topological edge modes exist not only on a plasma-vacuum interface but also on more general plasma-plasma interfaces. These findings broaden the possible applications of these exotic excitations in space and laboratory plasmas.

    Fu, Y., Qin, H. Nat Commun 12, 3924 (2021). https://doi.org/10.1038/s41467-021-24189-3
  • "Axisymmetric modes resonant at the tokamak X-points and a new fast ion instability", video
    F. Porcelli,Polytechnic University of Turin, Italy , abstract
    [#s1311, 24 Jun 2021]
    Axisymmetric MHD perturbations with toroidal mode number n=0 are resonant on the magnetic field lines going through the X-points of the tokamak divertor separatrix. As a consequence, current sheets form along the separatrix, which profoundly affect the stability of these modes. In particular, current sheets at the magnetic separatrix lead to the stabilization of n=0 vertical plasma displacements, at least on the ideal-MHD time scale, adding an important ingredient to the mechanism of passive feedback stabilization. A weakly damped n=0 mode, with a discrete oscillation frequency close to the poloidal Alfvén frequency, is also found. This mode may be driven unstable by the resonant interaction with fast ions, adding a new item to the catalogue of energetic particle driven instabilities.
  • Integration of the guiding-center equations in toroidal fields utilizing a local linearization approach, video
    M. Eder, Fusion@ÖAW, Institut für Theoretische Physik - Computational Physics, Technische Universität Graz, , abstract, slides
    [#s1309, 17 Jun 2021]
    A numerical integration method for guiding-center orbits of charged particles in toroidal fusion devices with three-dimensional field geometry as described in Ref. [1, 2] is presented. Here, high order interpolation of electromagnetic fields in space is replaced by a special linear interpolation, leading to locally linear Hamiltonian equations of motion with piecewise constant coefficients. This approach reduces computational effort and noise sensitivity while the conservation of total energy, magnetic moment and phase space volume is retained. The underlying formulation treats motion in piecewise linear fields exactly and thus preserves the non-canonical symplectic form. The algorithm itself is only quasi-geometric due to a series expansion in the orbit parameter. For practical purposes an expansion to the fourth order retains geometric properties down to computer accuracy in typical examples. When applied to collisionless guiding-center orbits in an axisymmetric tokamak and a realistic three-dimensional stellarator configuration, the method demonstrates correct long-term orbit dynamics. In Monte Carlo evaluation of transport coefficients, the computational efficiency of quasi-geometric integration is an order of magnitude higher than with a standard fourth order Runge-Kutta integrator. Moreover, the integration method is tested for the computation of fusion alpha losses in a realistic stellarator configuration. A Fortran program with the name “Guiding-center ORbit Integration with Local Linearization Approach” (GORILLA) is publicly available as Open Source code on GitHub [3]. References [1] M. Eder et al. 46th EPS Conf. on Plasma Physics, 2019, ECA Vol. 43C, P5.1100. [2] M. Eder et al. Physics of Plasmas 27, 122508 (2020), https://doi.org/10.1063/5.0022117 [3] M. Eder et al. GORILLA GitHub repository: https://github.com/itpplasma/GORILLA
  • Eliminating Finite Grid Instabilities in Gyrokinetic Particle-in-Cell Algorithms
    F. Holderied, IPP Garching, abstract
    [#s1298, 20 May 2021]
    In this talk, I will present a STRUcture-Preserving HYbrid code - STRUPHY - for the simulation of magneto-hydrodynamic (MHD) waves interacting with a small population of energetic particles far from thermal equilibrium (kinetic species). Such configurations appear for instance in deuterium-tritium fusion reactors, where hot α-particles or fast ions coming from external heating devices can resonantly interact with MHD waves and thus compromise confinement time. The implemented model features linear, ideal MHD equations in curved, three-dimensional space, coupled nonlinearly to the full-orbit Vlasov equations via a current coupling scheme (CCS). The implemented algorithm is based on finite element exterior calculus (FEEC) for MHD and particle-in-cell (PIC) methods for the kinetic part; it probably conserves mass, energy, and the divergence-free constraint for the magnetic field, irrespective of metric (= space curvature), mesh parameters and chosen order of the scheme. The motivation for this work stems from the need for reliable long-time simulations of energetic particle physics in complex geometries, covering the whole range of MHD waves. In STRUPHY, the finite element spaces are built from tensor products of univariate B-splines on the logical cuboid and can be made high-order by increasing the polynomial degree. Time-stepping is based on splitting of a skew-symmetric matrix with implicit sub-steps, mitigating CFL conditions from fast magneto-acoustic waves. High-order time splitting schemes can be used in this regard.
  • Eliminating Finite Grid Instabilities in Gyrokinetic Particle-in-Cell Algorithms
    Ben Sturdevant, PPPL, abstract, slides
    [#s1270, 01 Apr 2021]
    We explore the issue of finite grid (aliasing) instabilities in gyrokinetic particle-in-cell (PIC) algorithms. Theory for finite grid instabilities in momentum conserving PIC applied to full particle models including charge separation effects (e.g. Vlasov-Poisson) is well established [1], leading to the requirement of resolving the Debye length. Recent studies with momentum conserving PIC applied to quasi-neutral gyrokinetic models show that the situation can be worse in the sense that the instability is present for arbitrary spatial resolution [2,3]. We show, however, that a simple reformulation of the equations, making use of the continuity equation, eliminates this instability for all practical purposes. Our reformulation has connections to so-called energy-conserving PIC interpolations [4], the split-weight scheme [5], and the vorticity equation. In addition, we explore the effects of finite beta and finite drifts. We demonstrate that our approach is absolutely stable for static plasmas for any spatial resolution, and is stable for drifting plasmas with electron Mach numbers below unity (which is generally ensured  by ambipolarity in the plasmas of interest). [1] A. B. Langdon, J. Comput. Phys. 6 (1970) 247–267. doi:10.1016/0021- 9991(70)90024- 0. [2] G. J. Wilkie, W. Dorland, Phys. Plasmas 23 (2016) 052111. doi:10.1063/1.4948493. [3] B. F. McMillan, Phys. Plasmas 27 (2020) 052106. doi:10.1063/1.5139957. [4] D. C. Barnes, L. Chacon, Comput. Phys. Comm. 258 (2021) 107560. doi:10.1016/j.cpc.2020.107560. [5] I. Manuilskiy, W. W. Lee, Phys. Plasmas 7 (5) (2000) 1381. doi:10.1063/1.873955.
  • Collisionless mechanisms of plasma transport in the presence of stochastic open magnetic field lines
    Min-Gu Yoo, PPPL, abstract, slides
    [#s1269, 19 Mar 2021]
    The collisionless plasma transport in given stochastic magnetic fields has been studied for understanding the mechanisms of the thermal quench in tokamak disruption using a global gyrokinetic simulation code GTS. Previous studies have mostly focused on the dynamics of the passing particles along the open magnetic field lines during the thermal quench. However, we found that a considerable amount of the electrons (<60%) can be trapped in the device due to the magnetic mirror effect although the magnetic field lines are open to the wall. A high-resolution vacuum field analysis of the stochastic layer provides rich information regarding the 3-dimensional magnetic topology for understanding the characteristics of the plasma transport in such systems. In this study, we present a comprehensive picture of the relation between the plasma dynamics and the 3-D topology of the stochastic layer, which is essential to understand thermal quench physics. It was found that the consistent coupling of electron and ion dynamics through the ambipolar electric fields plays a critical role in determining the electron thermal energy transport. The 3-dimensional ambipolar potential builds up in the stochastic layer to keep the quasi-neutrality of the plasma during the thermal quench. The ambipolar potential produces the ExB vortices that mix the plasma across the magnetic field lines. The ExB mixing helps the high-energy trapped electrons to exit to the wall through the favorable open magnetic field lines. As a result, the electron temperature decreases steadily within the time scale of milliseconds.
  • Radar REMPI diagnostic in magnetic fields and low-pressure environments, video
    Christopher Galea, Princeton University, abstract
    [#s1233, 10 Feb 2021]
    Radar resonance-enhanced multi-photon ionization (Radar REMPI) is a remote diagnostic technique which selectively ionizes a species of interest using a tunable laser to create a laser-induced plasma (REMPI) and probes this plasma via coherent microwave scattering (radar) to obtain time-resolved information about the plasma. The high selectivity of the resonance-enhanced ionization allows us to detect a trace species of interest in a gas filled with other species so long as we tune onto an energy level resonance of the trace species. In this talk, we present the results of our theoretical and experimental investigation of the effects of magnetic fields and low pressures (below 1 Torr) on the Radar REMPI diagnostic. In the case of applied magnetic fields, we constructed a toy model which predicts magnetically induced depolarization of the scattered microwaves from the REMPI plasma and observe the effect in experiments. This finding suggests a novel capability of Radar REMPI for the standoff measurement of vector magnetic fields. In the case of low pressures, we observed experimentally that increasing microwave frequency results in a faster decay rate of the scattering signal, which would not be expected from coherent microwave scattering. This discrepancy is resolved by considering the effect of the spatial distribution of plasma on the phase of the scattered microwaves, which we call the decoherence effect. We will delve into the physics of the magnetic field and low-pressure effects and discuss potential applications.
  • Kelvin–Helmholtz waves at Earth’s magnetopause, video
    Shiva Kavosi, University of New Hampshire, abstract, slides
    [#s1234, 10 Feb 2021]
    Establishing the mechanisms by which the solar wind enters Earth's magnetosphere is one of the biggest goals of magnetospheric physics as it forms the basis of space weather phenomena. Magnetic reconnection is believed to be the dominant process by which solar wind plasma enters the magnetosphere. However, for periods of northward interplanetary magnetic field (IMF), reconnection is less likely at the dayside magnetopause, and Kelvin–Helmholtz waves (KHWs) are an important agent for plasma entry and for the excitation of ultra-low-frequency (ULF) waves. Over the past two decades, several space missions have enabled a leap forward in our understanding of this phenomenon at the Earth's magnetopause. Moreover, numerical MHD simulations have been extensively used to study the nonlinear evolution of the KH instability. This talk highlights the contributions to our understanding and the on-going research of KHWs and its role on the plasma transport across the magnetopause using THEMIS (Time History of Events and Macro scale Interactions during Substorms) mission, Magnetospheric Multiscale (MMS) mission and, a global MHD simulation, OpenGGCM.
  • Thermal ion orbit loss in diverted tokamaks and its role in the radial current balance
    Robert Brzozowski, UCLA, abstract
    [#s1213, 29 Jan 2021]
    Ion orbit loss has long been considered a potential player in triggering and/or sustaining the turbulence suppression necessary for the L-H transition. A loss cone model for the X-point mediated loss of thermal ions in diverted tokamaks finds that the loss plays a significant role in establishing the near-edge environment. The effective timescale of any one trajectory’s loss qualitatively changes if transit occurs faster than the ion would be scattered out of the loss region; above this threshold the rate of loss is driven by the rate of transport into the loss cone, while below this threshold the rate of loss is driven by the physical-space drift motions. The consideration of velocity-space relaxation processes along the loss region’s boundary serves both to demarcate the loss cone along this qualitative change in the effective timescale and to estimate the steady-state loss current across the resulting modified boundary. The orbit loss calculations are implemented into the transport code SOLPS as source terms, and the equilibrium radial current balance is investigated in L-mode plasmas near the L-H transition. Within the studied parameter space, the loss current is directly proportional to the ion temperature and exhibits low- and high-density branch behaviors for a given input power. Analysis indicates that the return current primarily consists of two components: an increased inward flow of ions associated with the perpendicular viscosity and a decreased outward diamagnetic flow of ions. The former drives the leading order plasma response, an increase in the shear of the electric field and plasma rotations in the vicinity of the separatrix, while the latter involves lower order poloidal redistributions of the pressure about a flux surface. Following these observations, experimental proposals on ASDEX Upgrade have been made to ascertain the role of ion orbit losses. [1] M. Laishram, D. Sharma, and P. K. Kaw, Phys. Rev. E 91, 063110 (2015). [2] M. Laishram, D. Sharma, P. K. Chattopdhyay, and P. K. Kaw, Phys. Rev. E 95, 033204 (2017). [3] M. Laishram, D. Sharma, and P. Zhu, Journal of Physics D: Applied Physics 53, 025204 (2019). [4] M. Laishram, Driven dust vortex characteristics in plasma with external transverse and weak magnetic field (2020), 2011.03237.
  • Modeling of experimentally observed driven dust vortex characteristics in laboratory plasma, video
    Modhuchandra Laishram, Institute for Plasma Research, India, abstract, slides
    [#s1192, 17 Dec 2020]
    A two-dimensional (2D) hydrodynamic model is developed for characterizing dynamics of driven-dissipative dust cloud confined in an axisymmetric toroidal setup along with an unbounded sheared streaming plasma [1, 2]. This formulation brings out several sources of the dust vorticity due to the background sheared plasma flow fields. The numerical solutions confirm the analytical structure of the driven dust vortex flow in the linear limit [1], but also fairly predict experimentally observed nonlinear characteristics of the dust rotation such as threshold structural bifurcation, the emergence of uniform vorticity core region, recovered scaling law for estimates the kinematic viscosity from the experimentally measurable quantities, and formation of steady-state multiple counter-rotating and co-rotating vortices [2, 3]. The hydrodynamic model is extended for analysis of driven vortex characteristics in presence of external transverse and weak magnetic field (B) in a planner setup and parametric regimes motivated by recent magnetized dusty plasma (MDP) experiments [4]. This analysis has shown that shear in the B can produce a sheared internal field (Ea) in between electrons and ions due to the E × B and ∇B × B-drifts that causes rotation of the dust cloud levitated in the plasma. The flow solution demonstrates that neutral pressure decides the dominance between the ions-drag and the Ea-force. The shear ions-drag generates an anti-clockwise circular vortical structure, whereas the shear Ea-force is very localized and gives rise to a clockwise D-shaped elliptical structure which turns into a meridional structure with decreasing B. In the regime of high pressure and lower B, the Ea-force becomes comparable or dominant over the ion drag and peculiar counter-rotating vortex pairs are developed in the domain. [1] M. Laishram, D. Sharma, and P. K. Kaw, Phys. Rev. E 91, 063110 (2015). [2] M. Laishram, D. Sharma, P. K. Chattopdhyay, and P. K. Kaw, Phys. Rev. E 95, 033204 (2017). [3] M. Laishram, D. Sharma, and P. Zhu, Journal of Physics D: Applied Physics 53, 025204 (2019). [4] M. Laishram, Driven dust vortex characteristics in plasma with external transverse and weak magnetic field (2020), 2011.03237.
  • Extended neoclassical plasma rotation and transport theory, video
    Cheonho Bae, Atlanta International School, abstract, slides
    [#s1191, 15 Dec 2020]
    Extended neoclassical plasma rotation theory has been developed based on the fluid moment equations with collisionality-extended Braginskii’s closure of the viscosity and the first-order poloidal asymmetries in velocities, densities, and electrostatic potential [1-3]. Major recent extensions [3,4] include rotation and transport calculations with generalized D-shaped flux surfaces using the Miller geometry [5] with the Shafranov shifts. Recent rotation calculations of DII-D and KSTAR discharges [3,4], using the nonlinear self-consistent neoclassical GTROTA code [6], indicated agreements to within <15% to the measurements except in the very edge (rho > .90), thus providing confidence on the related transport calculations such as the radial electric field and the poloidal asymmetries. Further extension efforts [7] are in progress with its major extension to Mikhailovskii-Tsypin’s closure of the viscosity and other edge physics to increase the accuracy on the edge rotation and transport calculations, and to include non-axisymmetric toroidal perturbations and develop a neoclassical toroidal viscosity formalism. Capabilities of the GTROTA code have also been extended to allow self-consistent rotation and transport calculations of up to four ion and electron species. Future plans on theory and code development, with its focus on edge rotation and transport, will be discussed. [1] W. M. Stacey, A. W. Bailey, D. J. Sigmar and K. C. Shaing, Nucl. Fusion 25 , 463 (1985). [2] W. M. Stacey and C. Bae, Phys. Plasmas 16, 082501 (2009). [3] C. Bae, W.M. Stacey, W.M. Solomon, Nucl. Fusion, 53 (2013) 043011. [4] C. Bae, W.M. Stacey, S.G. Lee, L. Terzolo, Phys. of Plasmas, 21 (2014). [5] R.L. Miller, M.S. Chu, J.M. Greene, et. al., Phys. of Plasmas, 5 (1998). [6] C. Bae, W.M. Stacey, T.D. Morley, Comp. Phys. Communications, 184 (2013). [7] W. M. Stacey and C. Bae, Phys. Plasmas 22, 062503 (2015).
  • The impact of two-fluid MHD stabilities on the transport of impurities in tokamak plasmas, video
    Jaeheon Ahn, abstract
    [#s1189, 03 Dec 2020]
    The plasma core is likely to be affected by 'sawteeth' playing significant role in core confinement and impurity transport. In particular, heavy impurities accumulated in the core plasma may lead to a radiative collapse. In order to understand and control their dynamics, the characteristics of compound sawteeth and the control of sawteeth relying on current or power deposition in the vicinity of the q=1 surface are explored by means of MHD simulations. In this work, the numerical tool used to simulate sawteeth is the XTOR-2F code, which is a non-linear tridimensional code solving two-fluid MHD equations. Neoclassical transport is shown to be important for heavy impurities in the core region. Meanwhile, ASDEX-U measurements show that the impurity dynamics in presence of sawteeth differ from the predictions made by transport codes. Fluid equations that model the transport of impurities in a highly collisional (Pfirsch-Schlüter) regime are implemented and coupled to the set of two-fluid MHD equations. The simulations show a difference between the impurity profiles with and without sawteeth, consistent with experimental observations. This results from a competition between neoclassical processes and sawtooth relaxations.
  • Numerical studies of sheared flow effects on visco-resistive MHD instabilities and application to ADITYA-U results, video
    Jervis Mendonca, abstract
    [#s1188, 03 Dec 2020]
    The effects of flow shear on the stability of a (2,1) tearing mode are examined using numerical and analytic studies on a number of model systems [1]. For a cylindrical reduced magnetohydrodynamic (MHD) model, linear computations using the CUTIE code show that sheared axial flows have a destabilizing effect, while sheared poloidal flows tend to reduce the growth rate of the mode. These effects are independent of the direction of the flow. For helical flows the sign of the shear in the flow matters. This symmetry breaking is also seen in the nonlinear regime where the island saturation level is found to depend on the sign of the flows. I have subsequently done this study for the visco resistive kink (m= 1, n= 1) mode in the RMHD regime [2].
  • Comparative Heliophysics: Planetary Space Weather in Our Solar System and Beyond
    Chuanfei Dong, abstract
    [#s1190, 02 Dec 2020]
    In the last two decades, the field of heliophysics has witnessed a tremendous surge. Research in heliophysics now encompasses a wide spectrum of fields ranging from space and magnetospheric physics to solar and stellar physics. One of the primary objectives of heliophysics is to understand the Sun and its interactions with planets (including Earth), where space weather plays an important role. In my presentation, I will begin with NASA’s Mars Atmosphere and Volatile EvolutioN (MAVEN) mission and delineate my study on the Martian responses to an interplanetary coronal mass ejection (ICME) - sometimes termed as solar storm - by using a sophisticated multifluid magnetohydrodynamic (MHD) code. I will then discuss how an active young Sun transformed Mars from an early warmer and wetter world to a desiccated and frigid planet with a tenuous atmosphere. In recent several years, it is well recognized that heliophysics plays an increasingly important role in the rapidly growing field of exoplanetary science. I will present my studies on exoplanetary atmospheric losses due to the impact of stellar winds and storms on planets residing within close-in habitable zones of M-type stars, which is a key factor that determines planetary habitability. Last but not least, I will introduce the ten- moment multifluid model developed at Princeton with its applications to the global magnetosphere of Mercury and discuss how this new model may enhance the science return from ESA-JAXA BepiColombo mission to Mercury where I serve as a Co-I.
  • Extending M3D-C1 to stellarator geometry: preliminary results
    Yao Zhou, PPPL, abstract, slides
    [#s1177, 19 Nov 2020]
    Stellarator plasmas have been observed to be nonlinearly stable even when driven beyond linear MHD stability thresholds. Hence, stellarator designs could employ nonlinear stability considerations to relax linear stability constraints, which can often be too restrictive and costly. However, this possibility has not been systematically investigated due to the lack of a state-of-the-art nonlinear initial-value MHD code for stellarators. We aim to fill this gap by extending the M3D-C1 code from tokamak to stellarator geometry. Our approach introduces a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework. The mapping from the logical to the physical (R, phi, Z) coordinates is then used to calculate derivatives in the latter, in terms of which the existing physics equations are written. This way, no significant changes to the extended MHD models within M3D-C1 are required. So far, we have successfully implemented this approach in 2D and verified its results against the existing code. Preliminary results in 3D will also be presented, including proof-of-principle resistive MHD simulations of simple stellarator plasmas such as a rotating ellipse.
  • First Principles Modeling of Point Defects in Bandgap Materials
    Bharat Medasani, PPPL, abstract, slides
    [#s1176, 05 Nov 2020]
    Point defects have a strong impact on the performance of semiconductor and insulator materials used in technological applications, spanning microelectronics to energy conversion and storage. They also play a critical role in the synthesis and growth of oxide films. The nature of the dominant defect types, how they vary with processing conditions, and their impact on materials properties are central aspects that determine the performance of a material in a certain application. This information is, however, difficult to access directly from experimental measurements. Consequently, computational methods, based on electronic density functional theory (DFT), have found widespread use in the calculation of point-defect properties. DFT based computational methods for point defects have inherent errors that require explicit corrections. In this talk I explain the DFT based modeling of point defects and the associated correction schemes using Cr2O3 as an example. Utilizing the points defect data, self-diffusion in Cr2O3 is evaluated.
  • Peeling-ballooning modes in low aspect ratio - towards a predictive edge pedestal model for spherical tokamaks
    Dr. Andreas Kleiner (PPPL), abstract, slides
    [#s1131, 15 Oct 2020]
    Edge-localized modes (ELMs) are a major concern for tokamak devices and their control is crucial for the operation of future reactor-scale machines. ELMs can be triggered when strong gradients are present in the edge transport barrier (edge pedestal). The width and height of the pedestal is often constrained by the occurrence of ideal-MHD peeling-ballooning modes [1] as well as kinetic ballooning modes (KBMs) in the pedestal region; this model has been successfully applied by the EPED model [2] to predict the pedestal height and width in conventional aspect ratio tokamaks. The predictions of the EPED model however often do not accurately describe observations in machines with low aspect ratio (e.g. spherical tokamaks). For ELMing discharges in NSTX the EPED model predicts stability. The reasons for this discrepancy might be associated with the limitation to ideal-MHD computations or the breakdown of the assumption that local ballooning theory well approximates the stability limit of kinetic ballooning modes in the edge of low aspect ratio plasmas. With the goal of obtaining a model to predict ELMs in spherical tokamaks and to find the limiting values for pedestal width and height, we model peeling-ballooning modes including non-ideal effects. The extended-MHD code M3D-C1 [3,4] is applied to determine the stability thresholds of peeling-ballooning modes in NSTX discharges. We mainly focus on the influence of plasma resistivity and rotation on edge stability, but also consider two-fluid effects. By varying pedestal parameters such as the pressure gradient and current density for given (ELMing) discharges, we are able to locate these discharges in parameter space relative to the stability boundary. In this context, we also identify the physics mechanisms that are important to describe these macroscopic edge modes in spherical tokamaks. We find robust resistive peeling-ballooning modes well before the ideal stability threshold is met. These modes thus extend the region of ideal peeling-ballooning instability in the investigated ELMing NSTX discharges. It is found that the actual NSTX plasmas are localized close to, or slightly within, the unstable side of the stability boundary calculated with our model. For large aspect ratio discharges the model is benchmarked with the ELITE code, which is employed in the frame of the EPED model. This study of macroscopic instabilities constitutes a first step towards a model to predict pedestal width and height in H-mode discharges in spherical tokamaks.
  • Learning long division: weak operations in a discontinuous Galerkin kinetic solver
    Mana Francisquez, MIT, abstract, slides
    [#s1120, 08 Jul 2020]
    The Gkeyll computational plasma physics code aims to provide a unified framework for fluid and kinetic studies using state-of-the-art discontinuous Galerkin (DG) methods. In developing Gkeyll we learned a number of lessons in DG theory and application, which are central to Gkeyll's algorithms and perhaps of interest to other DG workers. Some of these lessons stem from encountering operations as simple as division, so in this talk we will motivate DG division and how we have addressed that problem to deliver alias-free and conservative algorithms for, for example, kinetic collision operators. The concept of a weak equality, common to applied math and FEM communities, plays a central role. We will then discuss how weak equality has been leveraged in order to formulate and implement more complex operations, such as spectral transforms for piecewise discontinuous polynomial data and multigrid solutions of DG Poisson equations.
  • An Improved Equation-Free Projective Integration Method for Gyrokinetic Profile Evolution of Tokamak Plasmas
    Benjamin Sturdevant, PPPL, abstract, slides
    [#s1105, 02 Mar 2020]
    Cross-scale interactions between micro-scale transport processes and macro-scale profile evolution are difficult to study using brute-force long-time gyrokinetic simulations due to the large computational resources required and a possible accumulation of error. To address this, a multi-scale time integration method based on the equation-free projective integration method of Keverekidis and Gear [1] has been developed to accelerate kinetic simulations [2]. Previous attempts to use the equation-free method in kinetic plasma simulations were limited due to the appearance of spurious transient oscillations occurring after the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function [3]. In this work, a kinetically-informed lifting algorithm has been added to the equation-free method to mitigate this issue. This scheme has been verified in the gyrokinetic particle-in-cell code XGCa. It has been demonstrated to accurately reproduce profile evolution due to collisional neoclassical orbital dynamics while achieving a computational speed up of over 4x compared to brute force time stepping.

    References:

    [1] I.G. Kevrekidis, C.W. Gear, J. Hyman, P. Kevrekidis, O. Runborg, and C. Theodoropoulos, Comm. Math. Sci. 1 (4) (2003) 715-762.
    [2] B. Sturdevant, S.E. Parker, C.S. Chang, and R. Hager, submitted to Phys. Plasmas (2019).
    [3] M. A. Shay, J.F. Drake, B. Dorland, J. Comput. Phys. 226 (2007) 571-585
  • Global modelling of microinstabilities in stellarators and with electromagnetic effects using XGC
    Michael Cole, PPPL, abstract, slides
    [#s1104, 27 Feb 2020]
    Global gyrokinetic simulations are an increasingly important tool for understanding and designing magnetic confinement fusion devices. The recent operation of the optimized stellarator Wendelstein 7-X in a turbulence-dominated regime gives renewed urgency to the need for stellarator gyrokinetics in general, while electromagnetic effects are likely to play an important role in turbulence transport in both tokamaks and stellarators. In this talk, we report developments and deployment of the full volume gyrokinetic code XGC for stellarator physics and electromagnetic physics. XGC has been extended to treat general 3D toroidal equilibria with modifications to permit efficient electrostatic field solve. The code has then been used to reproduce linear electrostatic ITG calculations in Wendelstein 7-X and PPPL's QUASAR, followed by the first global turbulence transport simulations in QUASAR geometry. In addition, implementation of improved explicit electromagnetic techniques has allowed XGC to reproduce the transition from ion temperature gradient (ITG) to kinetic balloon mode (KBM) dominated microinstability regimes in tokamak cyclone base case geometry. This has then been extended to simulations of KBM in test geometry similar to NSTX-U. Recent flux tube simulations with the GENE code (Aleynikova 2018) show a similar transition from an ITG to a KBM dominated regime at beta values lower than those envisioned for peak Wendelstein 7-X performance and future reactors. One future use of the combined global electromagnetic stellarator gyrokinetic code will be to investigate this and predict turbulence transport in optimized stellarators.
  • Simulations, High Throughput Workflows, and Machine Learning in Materials Science
    Bharat Medasani, University of Delaware, abstract
    [#s1103, 25 Feb 2020]
    In the first part of talk, a classical density functional theory for modeling the ionic environment surrounding macromolecules in aqueous suspensions is presented. Our approach extends the capabilities of conventional approaches by accounting for electrostatic ion correlations, excluded volume correlations, and size asymmetry. I explain the computational schemes for 1d systems (obtained after applying symmetry) and 3d systems. Results for spherical and cylindrical models representing nanoparticles and proteins respectively are presented.
    In the second part, Fireworks, a popular workflow software for running simulations in high throughput at supercomputing centers is presented. Fireworks has been used to run more than 100 million CPU-hours worth of materials science and computational chemistry simulations at National Energy Research Supercomputing Center (NERSC). It supports many features relevant to enabling modern data-driven and high-throughput science applications. I also discuss my work in optimizing the performance of Fireworks.
    In the final part, I present a novel use of machine learning to predict defects in B2 intermetallics, a class of intermetallics containing two elements and with a specific crystal structure. Data on defects in materials is scarce due to the complexities involved in both experiments and simulations. Computational data for defects was generated using Fireworks based workflows. Using the generated data and by applying gradient boosting classification (which was then an innovative machine learning technique) defects in B2 intermetallics were predicted with a high accuracy.
  • Progress in theoretical understanding of the Dimits shift and the tertiary instability in drift-wave turbulence
    Hongxuan Zhu, PPPL, abstract, slides
    [#s1102, 20 Feb 2020]
    A generic understanding of the well-known Dimits shift in electrostatic drift-wave turbulence is obtained, for the first time, by studying the tertiary instability of a zonal flow within reduced turbulence models [Phys. Rev. Lett. 124, 055002 (2020)]. We show that tertiary modes are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. These modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate is derived within the modified Hasegawa-Wakatani model. We show that the growth rate equals the primary-instability growth rate plus a term that depends on the local flow "curvature", i.e., the second radial derivative of U; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This shift is the Dimits shift, which we find explicitly in the Terry-Horton limit, and our analytic predictions agree well with results from numerical simulations. Our theory of the tertiary instability also extends to other turbulence models. For example, the key features of the tertiary instability of ion-temperature-gradient drift waves are reproduced by our theory and verified by gyrokinetic simulations using the code GS2.

    Remote participation: https://zoom.us/j/832071905
  • Towards a quasi-dynamic model for 3D MHD based on energy minimisation and relaxation
    Adelle Wright, Australian National University, abstract, slides
    [#s1100, 13 Feb 2020]
    Recent experiments have found evidence of MHD activity in stellarators, including the observation of sawtooth-like oscillations during electron cyclotron current drive (ECCD) experiments in the second Wendelstein 7-X campaign [1]. In this talk, we present our approach to developing a quasi-dynamical model for 3D MHD that is based on energy minimisation and relaxation, and motivated by the need for predictive and computationally efficient global modelling of macroscopic dynamics in stellarators.

    We propose to approximate plasma evolution by a sequence of MHD equilibria which are connected via re-equilibrating relaxation events and consistent with broader physics setting of extended-MHD. A particular challenge, compared to previous cognate approaches [2], is to accommodate the variety of topological structures (e.g. magnetic islands and stochastic field regions) which can be supported in 3D MHD equilibria.

    For our equilibria we use the Multi-Region Relaxed MHD (MRxMHD) model, proposed by Dewar et al. [3], wherein a plasma is partitioned into a finite number of Taylor-relaxed, force-free (i.e. constant pressure) volumes. The volumes are separated by current sheet interfaces across which a pressure jump is supported leading to equilibria with discontinuous, stepped pressure profiles. We discuss an a priori prescription for determining the sequence of equilibria which is based on the successive break up of MRxMHD interfaces and exploits the Hamiltonian nature of magnetic fields.

    Next, we consider dynamical accessibility, i.e. whether successive equilibria in a quasi-dynamical sequence can be reached by the plasma in a way that is consistent with both MRxMHD and extend-MHD. Since MRxMHD states are hypothesised to form via volume-localised Taylor relaxation, we examine in detail the nonlinear dynamics of relaxation. Using a simple force-free model as a testbed, we discuss the development of diagnostics and analysis tools, with a view to quantitatively constraining the timescale and parameter regime in which such a quasi-dynamical model may be valid.

    [1] M. Zanini et al., EPJ Web of Conferences 203, 02013 (2019).
    [2] R. Clemente et al., Plasma physics and controlled nuclear fusion 1988 V.2 (1989).
    [3] R. L. Dewar et al., Journal of Plasma Physics 81.6 (2015).

    Remote participation: https://zoom.us/j/832071905
  • DCON for SPEC
    Alan H. Glasser, Fusion Theory and Computation , Inc. (FTCI), abstract, slides
    [#s1101, 12 Feb 2020]
    DCON is a widely-used ideal MHD stability code for axisymmetric toroidal plasmas, based on a generalization of Newcomb’s criterion for a cylindrical plasma. This presentation describes a further generalization to non-axisymmetric toroidal plasmas with toroidal periodicity and stellarator symmetry.

    This work was originally formulated for equilibria computed by the VMEC code, but that has encountered several problems. One is that VMEC imposes nested flux surfaces on the solution, which is unrealistic and fails to find a true energy minimum. A second is that it uses radial finite differences, uniformly spaced in toroidal magnetic flux, to solve for the magnetic field B, resulting in very poor resolution near the magnetic. A third is a failure to maintain div B = 0. These problems have proved fatal to the DCON approach.

    The current approach is based on equilibria computed by the more recent SPEC code, developed by Dewar, Hudson, et al. SPEC (Stepwise Pressure Equilibrium Code) partitions the plasma into infinitesimal interfaces with finite plasma pressure discontinuities but [[P + B2/2]] = 0; and volumes between the interfaces where the plasma relaxes to a Taylor state and may have regions of multiple small islands and stochasticity. Radial discretization uses Chebyshev polynomials, with regularization in the innermost volume near the axis of the coordinate system. It solves for the vector potential A rather than the magnetic field B, ensuring div B = 0.

    DCON for SPEC treats the interfaces as having finite width d and then takes the limit as δ -> 0.It derives a sequence of Euler-Lagrange equations for complex vectors of poloidal and toroidal Fourier components in successive interfaces and volumes, coupled together to produce a global solution and a generalized necessary and sufficient condition for ideal MHD stability. The derivation is complete and is being implemented in a new code.

    Remote participation: https://zoom.us/j/832071905
  • Stochastic electron acceleration in laser-plasma interactions
    Yanzeng Zhang, University of California San Diego, USA, abstract, slides
    [#s1099, 11 Feb 2020]
    The generation of energetic electron beams via the interaction of an intense laser pulse with underdense plasma is of great interest for many different applications, where the mechanisms of electron acceleration have been suggested and studied analytically, numerically, and experimentally over many years. It was revealed that the presence of self-generated (or externally applied) quasi-static electric and magnetic fields or a second counter-propagating weak laser pulse could significantly increase the electron energy well beyond the ponderomotive energy scaling (the maximum energy of an electron in the dominant laser only). However, due to the multidimensional spatio-temporal characteristics of the laser waves and strong nonlinearity of relativistic electron in these fields, the analytic investigations of the mechanism of electron acceleration in the earlier studies are quite limited, whereas the numerical simulations, which could also provide information of electron dynamics, are only valid within simulated parameter range. In this talk, I will introduce a novel approach to investigating the electron dynamics, especially paying attention to the stochastic electron acceleration, in different configurations of laser beams and quasi-static EM fields. The main idea of this method is to find proper canonical variables so that the new Hamiltonian describing the electron dynamics is time-independent without the appreciated perturbation. It can significantly simplify the analysis of the electron dynamics and allow us to utilize the fundamental results of previous studies on regular and stochastic motion in Hamiltonian systems. The physics underlying stochastic electron acceleration are revealed and the stochastic conditions are found by deriving the Chirikov-like mappings, which introduce the upper limit of the electron energy gained from the stochastic electron acceleration.

    Remote participation: https://zoom.us/j/832071905
  • Collective Thomson Scattering Diagnostic on Wendelstein 7-X
    Ivana Abramovic, Eindhoven University of Technology, The Netherlands, abstract, slides
    [#s1098, 27 Jan 2020]
    The most advanced and the World’s only optimized stellarator, the Wendelstein 7-X (W7-X) experiment in Greifswald (Germany), came into operation in 2015 with the goal to explore the viability of the optimized stellarator concept as a fusion reactor. Adequate diagnostics of the key plasma parameters are of crucial importance for fulfilling this goal. The research conducted over the course of this project concerns the development of theory and modelling of collective Thomson scattering diagnostic on W7-X. Collective Thomson scattering (CTS) is a powerful microwave diagnostic sensitive primarily to ion dynamics and capable of simultaneously measuring a number of important parameters such as: ion temperature, plasma composition, drift velocities and fast ion population. A CTS measurement results in a spectrum from which values of desired plasma parameters can be inferred. To this end a forward model, eCTS, of the scattering was developed which allows calculation of synthetic CTS spectra based on a set of input parameters. The forward model was integrated into a Bayesian scientific framework which allows one to solve the inverse problem of inferring posterior probability distributions of input parameters of the forward model given the measured spectra. eCTS has successfully been used to analyze the first CTS spectra obtained on W7-X. Additionally neural networks were trained on databases produced by eCTS in a relevant parameter space in order to provide fast inference of bulk ion temperatures as an alternative for Bayesian inference.
    The work also furthers CTS theory by extending the set of parameters measurable by the CTS diagnostic to include the radial electric field and enables exploration of anisotropic temperature effects. Finally, the work explores experimental challenges hindering optimal performance of CTS as well as data analysis related challenges which stand in the way of using this diagnostic to its full potential.

    Remote participation: https://zoom.us/j/832071905
  • Excitation and Damping Mechanisms of Geodesic Acoustic Modes in Tokamaks
    Ivan Novikau, Max-Planck-Institut for Plasma Physics (IPP), Garching, Germany, abstract, slides
    [#s1097, 23 Jan 2020]
    Turbulence in tokamaks generates radially sheared zonal flows (ZFs). Their oscil- latory counterparts, geodesic acoustic modes (GAMs), appear due to the action of the magnetic field curvature. The zonal structures can reduce the radial transport in toka- mak plasma acting as a sink for the turbulence energy through the inverse energy cascad- ing. The GAMs can also be driven by an anisotropic energetic particle (EP) population leading to the formation of global radial structures, called EGAMs. The EGAMs play a role of an intermediate agent between the EPs and thermal plasma, by redistributing EP energy to the bulk plasma through the collisionless wave-particle interaction. In such a way, the EGAMs might contribute to the plasma heating. Thus, investigation of EGAM characteristics, especially in the velocity space, is necessary for precise understanding of the transport phenomena in tokamak plasmas. In addition to ion Landau damping, both the GAMs and EGAMs have been found to be subject to the electron Landau damping. To investigate the influence of the electron damping on the GAMs and EGAMs, a Mode-Particle-Resonance (MPR) diagnostic has been implemented in the global gyrokinetic particle-in-cell code ORB5. This enables to investigate the relative importance and the evolution of wave-particle resonances responsible for the ion and electron Landau damping, and of the EP drive. An ASDEX upgrade (AUG) discharge #31213 is chosen as a reference case for this investigation due to its rich EP nonlinear dynamics.
    Nonlinear GAM excitation by turbulence is also considered. Using realistic temper- ature and density profiles as well as magnetic equilibrium of AUG shot #20787, GAM excitation by ITG instabilities is shown with comparison between numerical and exper- imental GAM frequency spectra. Finally, excitation of a global GAM-like structure in a TCV magnetic configuration is presented as well.

    Remote participation via Zoom: https://zoom.us/j/4634811292
  • Electron Modeling for MHD Simulation in Fusion Plasmas
    Dongjian Liu, Sichuan University, Chengdu, China, abstract, slides
    [#s1096, 22 Jan 2020]
    A series of electron model has been developed to study the low frequency electromagnetic modes in magnetized plasmas[1] [2]. Coupled to the gyrokinetic ions, the global gyrokinetic particle simulation of tearing modes have been developed and verified in the gyrokinetic toroidal code (GTC). GTC linear simulations in the fluid limit of the kink-tearing, resistive tearing modes and collisonless tearing mode in the cylindrical geometry agree well with the magnetohydrodynamic eigenvalue and initial value codes.

    References:
    [1] Dongjian Liu and Liu Chen, Plasma Physics and Controlled Fusion, 53 062002 (2011)
    [2] Dongjian Liu, Wenlu Zhang, Joseph McClenaghan, Jiaqi Wang and Zhihong Lin, Physics of Plasmas 21,122520 (2014).

    Remote participation via Zoom: https://zoom.us/j/4634811292
  • Generation of suprathermal electrons by collective processes in Maxwellian plasmas
    Sabrina F. Tigik, Instituto de Física, Universidade Federal do Rio Grande do Sul,Brazil , abstract
    [#s1095, 06 Jan 2020]
    Weak turbulence theory is the standard formalism for treating nonlinear, low amplitude, kinetic instabilities in collisionless plasmas; therefore, it is concerned only with the description of collective plasma oscillations. However, the long-lasting timescale of nonlinear processes suggests that collisions might affect the late plasma dynamics, acting alongside nonlinear oscillatory effects [1]. Such cases, where collective and collisional processes coexist in the plasma dynamics, have been rigorously addressed only recently [2]. One of the outcomes of this extended formulation is a then-unknown mechanism named electrostatic bremsstrahlung emission. The electrostatic bremsstrahlung is a kind of transient radiation, emitted in all spectrum, caused by continuous interparticle interactions. The portion irradiated in the eigenmode frequency range is capable of altering the wave spectrum, which will then modify the velocity distribution through wave-particle resonance. Considering emissions in the Langmuir wave frequency range, and Maxwellian electrons as the initial state, I have analyzed the time evolution of the system. The result was the scattering off of a small fraction of the thermal electrons in the resonance region, to high velocities, creating a suprathermal population. The tail grows consistently in the beginning, but after long integration time, the growth rate slows down, indicating that the system is arriving at a new asymptotic equilibrium state [3].

    References: [1] S. F. Tigik, L. F. Ziebell, P. H. Yoon, and E. P. Kontar. Two-dimensional time evolution of beam-plasma instability in the presence of binary collisions. Astronomy & Astrophysics, 586:A19, February 2016.
    [2] P. H. Yoon, L. F. Ziebell, E. P. Kontar, and R. Schlickeiser. Weak turbulence theory for collisional plasmas. Physical Review E, 93:033203, Mar 2016.
    [3] S. F. Tigik, L. F. Ziebell, and P. H. Yoon. Generation of suprathermal electrons by collective processes in collisional plasma. The Astrophysical Journal Letters, 849(2):L30, 2017.

    Remote participation via Zoom: https://zoom.us/j/8320719053
  • Flow effects on visco-resistive MHD in a periodic cylindrical tokamak
    Jervis Mendonca, Institute for Plasma Research, HBNI, India , abstract
    [#s1094, 03 Jan 2020]
    Flow and viscosity significantly modify resistive modes in a tokamak, and I have investigated these using the CUTIE code. These studies indicate that flow can be used to improve plasma duration and quality in a tokamak, and this has motivated my investigation. In this presentation, I have begun by studying the (2,1) tearing mode and found several new results, namely, that nature of stabilisation depended on whether axial or poloidal flows were used. I also observed that the sign of shear in helical flow mattered. This symmetry breaking is also seen in the nonlinear regime where the island saturation level is found to depend on the sign of the flows. I have proceeded to study the (1,1) internal kink mode using a Visco-Resistive MHD(V-RMHD) model. I have observed here that stabilisation due to axial flows in particular are affected by the viscosity regime. Symmetry breaking at higher viscosity in linear growth rates and nonlinear saturation levels as well is observed. In summary, for axial, poloidal, and most helical flow cases, there is flow induced stabilisation of the nonlinear saturation level in the high viscosity regime and destabilisation in the low viscosity regime. I have continued my studies in the two fluid regime. In the linear regime, we have studied how the growth rate as well as diamagnetic flow frequency of the modes changes due to fluid effects for a range of viscosity and resistivity values. I have also found diamagnetic drift stabilisation of the (1,1) mode in the two fluid case, that is, we have seen the growth rate of the (1,1) mode reduces with an increase in density gradient. In the nonlinear case, I investigate the evolution of the mode with imposed axial flow, poloidal and helical flows. I find the viscosity regime affects the nonlinear saturation regime.

    Remote participation via Zoom: https://zoom.us/j/8320719053
  • Wide q95 Windows for Edge-Localized-Mode Suppression by Resonant Magnetic Perturbations in the DIII-D Tokamak
    Qiming Hu, PPPL, abstract, slides
    [#s1075, 12 Dec 2019]
    Edge-Localized-Mode (ELM) suppression by resonant magnetic perturbations (RMPs) typically occurs over narrow ranges in the plasma magnetic safety factor q95, however wide q95 windows of ELM suppression are favourable for the robust avoidance of ELMs in ITER and future reactors. Here we show from experiment and nonlinear two-fluid MHD simulation that wide q95 windows for ELM suppression are accessible in low-collisionality DIII-D plasmas relevant to ITER high power mission so long as the applied RMP strength exceeds the threshold for resonant field penetration at the pedestal top by ~1.5X. When the applied RMP is close to the threshold for resonant field penetration at the top of the pedestal, only isolated magnetic islands form near the pedestal top, producing narrow q95 windows of ELM suppression (dq95 ~ 0.1). However, as the threshold for field penetration decreases relative to the RMP amplitude, then multiple magnetic islands can be driven on adjacent rational surfaces near the pedestal top. Multiple magnetic islands lead to q95 window merging and wide regions of ELM suppression (up to dq95 ~ 0.7 seen in DIII-D). Nonlinear MHD simulations are in quantitative agreement with experiment and predict improved access to wide windows of ELM suppression in DIII-D by using n = 4 RMPs, due to the denser rational surfaces and the weak dependence of the penetration threshold and island width on toroidal mode number.
  • Regimes of weak ITG/TEM modes for transport barriers without velocity shear
    Mike Kotschenreuther, University of Texas, Austin , abstract, slides
    [#s1093, 11 Dec 2019]
    TBD
  • A deep dive into the distribution function: understanding phase space dynamics using continuum Vlasov-Maxwell simulations
    J. Juno, University of Maryland , abstract, slides
    [#s1091, 05 Dec 2019]
    In collisionless and weakly collisional plasmas, the particle distribution function is a rich tapestry of the underlying physics. However, actually leveraging the particle distribution function to understand the dynamics of a weakly collisional plasma is challenging. The equation system of relevance, the Vlasov-Maxwell system of equations, is difficult to numerically integrate, and traditional methods such as the particle-in-cell method introduce counting noise into the distribution function.
    But motivated by the physics contained in the distribution function, we have implemented a novel continuum Vlasov-Maxwell method in the Gkeyll framework [1,2]. The algorithm uses the discontinuous Galerkin finite element method to produce a high order accurate solution, and a number of algorithmic breakthroughs have been made to produce a robust, cost efficient solver for the Vlasov-Maxwell system of equations.
    In this talk, I will present both our algorithmic work and the phase space dynamics of a number of plasma processes, focusing on two applications which make manifest the power and utility of a continuum Vlasov-Maxwell method. I will demonstrate how the high fidelity representation of the distribution function, combined with novel diagnostics, permits detailed analysis of the energization processes in a perpendicular collisionless shock. In addition, I will show a set of recent results on the evolution of kinetic instabilities driven by unstable beams of plasma, where the generation of small scale structures in velocity space has dramatic consequences for the overall macroscopic evolution of the plasma[3]. To further motivate the development of the continuum Vlasov-Maxwell method, I will also show how particle noise can modify the dynamics of the latter set of simulations with analogous particle-in-cell simulations.

    [1] J. Juno, A. Hakim, J. TenBarge, E. Shi, W. Dorland (2018). Discontinuous Galerkin algorithms for fully kinetic plasmas. Journal of Computational Physics, 353, 110—147.
    [2] A. Hakim, M. Francisquez, J. Juno, G. Hammett. “Conservative Discontinuous Galerkin Discretizations of Fokker-Planck Operators.” submitted to Journal of Plasma Physics, 2019
    [3] V. Skoutnev, A. Hakim, J. Juno, J. TenBarge. “Temperature-dependent Saturation of Weibel-type Instabilities in Counter-streaming Plasmas.” Astrophysical Journal Letters, 872, 2. (2019)
  • More Intelligent Execution of Ensembles of Nonlinear Kinetic Simulations: learning and improving efficiency and accuracy as we go - BARS, SHAD BARS and mini-BARS
    Bedros Afeyan, Polymath Research Inc.
    [#s1092, 04 Dec 2019]
  • Machine learning of equivariant functions inspired by atomistic modelling with a special view towards fusion materials simulations
    Bastiaan J. Braams, Centrum Wiskunde & Informatica (CWI),The Netherlands , abstract, slides
    [#s1089, 27 Nov 2019]
    Over the past several years big data methods, including but not limited to use of deep convolutional neural networks, have been very successful in computer science applications and there is increasing effort to apply big data or machine learning methods to problems in physical science and engineering. Conversely we are seeing that problems from physical science are influencing machine learning research done in computer science environments. A very important application of big data methods in physical science where we see this mutual influence is the construction of effective interatomic potentials and force fields for atomistic modelling of molecular and condensed phase systems (e.g. [1]). This application shares features with certain applications in three-dimensional image processing in having data associated with point clouds and in seeking to represent functions that are invariant or covariant with respect to a permutation group (applied to the labelling of points in the cloud) and with respect to spatial groups of translations, rotations and inversion. Some by now almost classical big data approaches to the atomistic problem include use of Gaussian process approximation (kernel ridge regression) [2] and use of spherical wavelet expansions [3]. In addition deep neural networks are being applied (e.g. [4], [5]) and here we see the closest link to machine learning research with key words such as Point Cloud Convolutional Networks, Deep Sets, Spherical CNNs, Tensor Field Networks and Gauge Equivariant Neural Networks [6]. The presentation will provide a survey of these machine learning developments in the context of the application in physical science with a special view towards radiation damage simulations for fusion.

    References:
    [1] Ceriotti, Michele. "Atomistic machine learning between predictions and understanding." arXiv preprint arXiv:1902.05158 (2019).
    [2] Bartók, Albert P., and Gábor Csányi. "Gaussian approximation potentials: A brief tutorial introduction." International Journal of Quantum Chemistry 115, no. 16 (2015): 1051-1057.
    [3] Eickenberg, Michael, Georgios Exarchakis, Matthew Hirn, Stéphane Mallat, and Louis Thiry. "Solid harmonic wavelet scattering for predictions of molecule properties." The Journal of chemical physics 148, no. 24 (2018): 241732.
    [4] Schütt, Kristof T., Huziel E. Sauceda, P-J. Kindermans, Alexandre Tkatchenko, and K-R. Müller. "SchNet–A deep learning architecture for molecules and materials." The Journal of Chemical Physics 148, no. 24 (2018): 241722.
    [5] Zhang, Linfeng, Jiequn Han, Han Wang, Wissam Saidi, Roberto Car, and Weinan E. "End-to-end symmetry preserving inter-atomic potential energy model for finite and extended systems." In Advances in Neural Information Processing Systems, pp. 4441-4451. 2018.
    [6] Extensive further references to be provided in the presentation.
  • An adjoint approach for the shape gradients of 3D MHD equilibria
    Elizabeth Paul, University of Maryland , abstract, slides
    [#s1090, 25 Nov 2019]
    The design of modern stellarators often employs gradient-based optimization to navigate the high-dimensional spaces used to describe their geometry. However, computing the numerical gradient of a target function with respect to many parameters can be expensive. The adjoint method allows these gradients to be computed at a much lower cost and without the noise associated with finite differences. In addition to gradient-based optimization, the derivatives obtained from the adjoint method are valuable for local sensitivity analysis and tolerance quantification.
    A continuous adjoint method has been developed for obtaining the derivatives of functions of the MHD equilibrium equations with respect to the shape of the boundary of the domain or the shape of the electro-magnetic coils [1]. This approach is based on the generalization of the self-adjointness of the linearized MHD force operator. The adjoint equation corresponds to a perturbed force balance equation with the addition of a bulk force, rotational transform, or toroidal current perturbation. We numerically demonstrate this approach by adding a small perturbation to the non-linear VMEC [2] solution, obtaining an order 102-103 reduction in cost in comparison with a finite difference approach. Examples are presented for the shape gradient of the rotational transform and vacuum magnetic well, a proxy for MHD stability. The adjoint solution required for the magnetic ripple, a proxy for near-axis quasisymmetry, requires the addition of an anisotropic pressure tensor to the MHD force balance equation. This modification has been implemented in the ANIMEC [3] code. We furthermore demonstrate that this adjoint approach can be applied to compute shape gradients of two important figures of merit [4], the departure from quasisymmetry and the effective ripple in the low-collisionality neoclassical regime, but require the development of new equilibrium solvers. Finally, initial steps toward adjoint solutions with a linearized equilibrium approach will be presented.

    [1] Antonsen, T.M., Paul, E.J. & Landreman, M. 2019 Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Journal of Plasma Physics 85 (2).
    [2] Hirshman, S.P. & Whitson, J.C. 1983 Steepest descent moment method for three-dimensional magnetohydrodynamic equilibria. Physics of Fluids 26 (12), 3553.
    [3] Cooper, W.A., Hirshman, S.P., Merazzi, S. & Gruber, R. 1992 3D magnetohydrodynamic equilibria with anisotropic pressure. Computer Physics Communications 72 (1),1–13.
    [4] Paul, E.J., Antonsen, T.M., Landreman, M., Cooper, W.A. Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria, Part II: Applications. Submitted to Journal of Plasma Physics.
  • Learning data driven discretizations for partial differential equations
    Stephan Hoyer, Google, abstract
    [#s1077, 14 Nov 2019]
    The numerical solution of partial differential equations (PDEs) is challenging because of the need to resolve spatiotemporal features over wide length and timescales. Often, it is computationally intractable to resolve the finest features in the solution. The only recourse is to use approximate coarse-grained representations, which aim to accurately represent long-wavelength dynamics while properly accounting for unresolved small scale physics. Deriving such coarse grained equations is notoriously difficult, and often ad hoc. Here we introduce data driven discretization, a method for learning optimized approximations to PDEs based on actual solutions to the known underlying equations. Our approach uses neural networks to estimate spatial derivatives, which are optimized end-to-end to best satisfy the equations on a low resolution grid. The resulting numerical methods are remarkably accurate, allowing us to integrate in time a collection of nonlinear equations in one spatial dimension at resolutions 4-8x coarser than is possible with standard finite difference methods.
  • Variational Principles and Applications of Local Topological Constants of Motion for Non-Barotropic Magnetohydrodynamics
    Asher Yaholom, Ariel University, Israel, abstract, slides
    [#s1076, 12 Nov 2019]
    Variational principles for magnetohydrodynamics (MHD) were introduced by previous authors both in Lagrangian and Eulerian form. In this presentation we introduce simpler Eulerian variational principles from which all the relevant equations of non-barotropic MHD can be derived for certain field topologies. The variational principle is given in terms of five independent functions for non-stationary non-barotropic flows. This is less than the eight variables which appear in the standard equations of barotropic MHD which are the magnetic field B ⃗ the velocity field v ⃗, the entropy s and the density ρ. The case of non-barotropic MHD in which the internal energy is a function of both entropy and density was not discussed in previous works which were concerned with the simplistic barotropic case. It is important to understand the rule of entropy and temperature for the variational analysis of MHD. Thus we introduce a variational principle of non-barotropic MHD and show that five functions will suffice to describe this physical system.
    We will also discuss the implications of the above analysis for topological constants. It will be shown that while cross helicity is not conserved for non-barotropic MHD a variant of this quantity is. The implications of this to non-barotropic MHD stability is discussed.

    [1] Asher Yahalom, “Simplified Variational Principles for non-Barotropic Magnetohydrodynamics,” J. Plasma Phys. (2016), vol. 82, 905820204.
    [2] Asher Yahalom, “Non-Barotropic Magnetohydrodynamics as a Five Function Field Theory,” International Journal of Geometric Methods in Modern Physics, No. 10 (November 2016), Vol. 13 1650130.
    [3] Asher Yahalom, “A Conserved Local Cross Helicity for Non-Barotropic MHD,” Journal of Geophysical & Astrophysical Fluid Dynamics (2017), Vol. 111, No. 2, 131–137.
    [4] Asher Yahalom, “Non-Barotropic Cross-helicity Conservation Applications in Magnetohydrodynamics and the Aharanov - Bohm effect," Fluid Dynamics Research (2017), Volume 50, Number 1, 011406.
  • First-principle formulation of resonance broadened quasilinear theory near an instability threshold
    Vinícius Duarte, PPPL, abstract, slides
    [#s1069, 19 Sep 2019]
    A method is developed to analytically determine the resonance broadening function in quasilinear theory, due to either Krook or Fokker-Planck scattering collisions of marginally unstable plasma systems where discrete resonance instabilities are excited without any mode overlap. It is demonstrated that a quasilinear system that employs the calculated broadening functions reported here systematically recovers the nonlinear growth rate and mode saturation levels for near-threshold plasmas previously calculated from kinetic theory. The distribution function is also calculated, which enables precise determination of the characteristic collisional resonance width.
    [Based on preprint V. Duarte, N. Gorelenkov, R. White and H. Berk, The collisional resonance function in discrete-resonance quasilinear plasma systems, arXiv:1906.01780v1 (2019)]
  • Global gyrokinetic PIC simulations of electromagnetic perturbations in fusion plasmas
    Alexey Mishchenko, IPP, abstract, slides
    [#s1068, 12 Sep 2019]
    An overview of European gyrokinetic PIC codes will be presented with a focus on numerical schemes used in the electromagnetic regime. Simulations of Alfven modes destabilized by the fast particles and of the zonal flows generated by the unstable Alfven waves will be discussed. Gyrokinetic PIC simulations of the MHD-type instabilities will be addressed on example of the internal kink modes. Simulations of the transition from the ITG-dominated regime to the KBM regime and finally to the microtearing instabilities with plasma beta increasing will be shown.
  • Next-gen distributed computing networks in support of research, discovery and innovation
    Dan Desjardin, Royal Military College of Canada/Distributed Compute Labs, abstract
    [#s1065, 15 Aug 2019]
    Distributed Compute Labs (DCL) is a Canadian educational nonprofit organization responsible for developing and deploying the Distributed Compute Protocol (DCP), a lightweight, easy-to-use, idiomatic, and powerful computing framework built on modern web technology, that allows any device — from smartphones to enterprise web servers — to contribute otherwise idle CPU and GPU capacity to secure and configurable general-purpose computing networks. By leveraging existing devices and infrastructure — a university’s desktop fleet, for example — a large supply of latent computational resources becomes available at no additional capital cost. DCP makes it possible for everyone — from a student in Santa Barbara to a large enterprise in New York City — to have access to large quantities of cost-effective computing resources. In summary, the Distributed Compute Protocol democratizes access to digital infrastructure, reduces barriers, and unleashes innovation.
  • Phase-space theory of the Dimits shift and cross-scale interactions in drift-wave turbulence
    Ilya Dodin, PPPL, abstract, slides
    [#s1066, 08 Aug 2019]
    Interactions of drift-wave turbulence with zonal flows, which are of interest due to their effect on fusion-plasma confinement, can be elucidated by using phase-space methods from quantum theory [1-7]. In this talk, I will show how applying these methods: (i) helps explain the cross-scale interactions between ITG and ETG turbulence seen in gyrokinetic simulations and also (ii) leads to a semi-quantitative analytic theory of the zonal-flow stability and the Dimits shift within the Terry-Horton model of drift-wave turbulence.

    [1] H. Zhu, Y. Zhou, and I. Y. Dodin, Nonlinear saturation and oscillations of collisionless zonal flows, New J. Phys. 21, 063009 (2019).
    [2] Y. Zhou, H. Zhu, and I. Y Dodin, Formation of solitary zonal structures via the modulational instability of drift waves, Plasma Phys. Control. Fusion 61, 075003 (2019).
    [3] D. E. Ruiz, M. E. Glinsky, and I. Y. Dodin, Wave kinetic equation for inhomogeneous drift-wave turbulence beyond the quasilinear approximation, J. Plasma Phys. 85, 905850101 (2019).
    [4] H. Zhu, Y. Zhou, and I. Y. Dodin, On the Rayleigh-Kuo criterion for the tertiary instability of zonal flows, Phys. Plasmas 25, 082121 (2018).
    [5] H. Zhu, Y. Zhou, and I. Y. Dodin, On the structure of the drifton phase space and its relation to the Rayleigh-Kuo criterion of the zonal-flow stability, Phys. Plasmas 25, 072121 (2018).
    [6] H. Zhu, Y. Zhou et al., Wave kinetics of drift-wave turbulence and zonal flows beyond the ray approximation, Phys. Rev. E 97, 053210 (2018).
    [7] D. E. Ruiz, J. B. Parker et al., Zonal-flow dynamics from a phase-space perspective, Phys. Plasmas 23, 122304 (2016).
  • Gyrokinetic study of RMP-driven plasma density and heat transport in tokamak edge plasma using MHD screened RMP field
    Robert Hager, PPPL, abstract, slides
    [#s1064, 01 Aug 2019]
    ITER plans to rely on RMP coils as the primary means for ELM control. However, puzzling observations on present-day experiments complicate understanding the underlying physics: plasma density is pumped out, which can lower the fusion efficiency in ITER, while the electron heat transport is still low in the pedestal. Kinetic level understanding including most of the important physics is needed but has not been available. Gyrokinetic total-f simulation of the plasma transport driven by n=3 resonant magnetic perturbations (RMPs) in a DIII-D H-mode plasma is performed using the gyrokinetic code XGC. The RMP field is calculated in M3D-C1 and coupled into XGC in realistic divertor geometry with neutral particle recycling. The RMP field is stochastic around the pedestal foot but exhibits good KAM surfaces at pedestal top and steep-slope. The simulation qualitatively reproduces the experimental phenomena: plasma density is pumped-out due to enhanced electrostatic turbulence while electron heat transport is low. Different from earlier gyrokinetic studies, the present simulation consistently combines neoclassical and turbulent transport, a fully nonlinear Fokker-Planck collision, neutral particle recycling, and the full 3-D electric field. Density pump-out is not seen without turbulence effects. Reduction of the ExB shearing rate is likely to be responsible, mostly, for the enhanced edge turbulence, which is found to be from trapped electron modes.
  • Velocimetry and the aperture problem for 2D incompressible flows
    Tim Stoltzfus-Dueck, abstract, slides
    [#s1028, 16 May 2019]
    The inference of velocity fields from 2D movies evolving conserved scalars (optical flow) is fundamentally ambiguous due to the well-known “aperture problem”: velocities along isocontours of the scalar are not visible. This may even corrupt the inference of velocity fields averaged at scales longer than the typical length scale of features in the scalar field, as in the drift wave. However, for divergence-free flows, a stream-function formulation allows us to show that the "invisible velocity" vanishes in the surface average inside any closed scalar isocontour. This error-free averaged velocity may be used as an “anchor” for a more reliable inference of the larger-scale velocity field, or to test model-based optical-flow schemes. We have also used the stream-function formulation to derive a new method of optical flow for divergence-free flows. We discuss the new algorithm, including details of discretization, boundary conditions, and image preprocessing that can significantly affect its performance. A simple implementation of the new method is shown to work well for a number of synthetic movies, and is also applied to a GPI movie of edge turbulence in NSTX.
  • Realisability of discontinuous MHD equilibria
    Adelle Wright, Australian National University, Canberra, abstract, slides
    [#s1034, 14 May 2019]
    Smooth 3D MHD equilibria with non-uniform pressure may not exist but, mathematically, there exist 3D MHD equilibria with non-uniform, stepped pressure profiles. The pressure jumps occur at surfaces with highly irrational values of rotational transform and generate singular current sheets. If physically realisable, how such states form dynamically remains to be understood. To be physically realisable states, MHD equilibria must exist for some non-trivial timescale, meaning they must be at least be ideally stable and sufficiently stable to the fastest growing resistive instabilities. This presentation will discuss recent progress towards understanding discontinuous MHD equilibria via a stability analysis of a continuous cylindrical equilibrium model with radially localised pressure gradients, which examines how the resistive stability characteristics of the model change as the localisation of pressure gradients is increased to approach a discontinuous pressure profile in the zero-width limit.
  • Kinetic Effects on Adiabatic Index via Geodesic Acoustic Continuum Calculations
    Fabio Camilo de Souza, University of Sao Paulo, Brazil, abstract, slides
    [#s1027, 08 May 2019]
    Geodesic Acoustic Mode (GAM) is primarily an electrostatic oscillation with the dominant toroidal N=0 and poloidal M=0,1 mode numbers. First predicted within the magnetohydrodynamic (MHD) theory [1], it's frequency is proportional to the square root of the adiabatic index gamma. GAMs can be treated by the kinetic theory, to include such parameters as the plasma rotation of different species [2], fast ions with bump-on-tail like [3] or slowing down [4] distribution function, and other models, to investigate the plasma equilibrium conditions.
    NOVA [5] is an ideal MHD code that computes the Alfvénic and acoustic continua and eigenmodes. The adiabatic index in NOVA is a fixed parameter, typically 5/3. The acoustic oscillations are calculated using the prescribed gamma value and its coupling with others continua.
    As in kinetic calculations it is possible to include other effects for more accuracy of the GAM continuum. We modified NOVA to include a profile for gamma, given as a function of the magnetic surfaces. This modification makes the MHD acoustic continuum to match the Kinetic one. This kinetic gamma allows to compute more accurate eigenmodes which are strongly coupled to GAM continuum structure. This conclusion is similar to all low frequency oscillations including GAMs, the Alfvén-acoustic BAAE modes and others. The understanding of the impact of discharge parameters in these modes can improve the plasma transport control.
    Simulation results for DIII-D will be presented, an Energetic Particle induced GAM appears in the maximum of the modified GAM continuum, the result matches the observed frequency for N=0 oscillation. Implications to the experiments in NSTX will be discussed.

    [1] N. Winsor, et. al., Phys. Fluids 11, 2448 (1968)
    [2] A.G. Elfimov, et. al., Phys. Plasmas 22, 114503 (2015)
    [3] F. Camilo de Souza, et. al, Phys. Lett. A 381, 3066 (2017)
    [4] Z. Qiu, et. al., Plasma Phys. Controlled Fusion 52, 095003 (2010)
    [5] C. Z. Cheng and M. S. Chance, Journal of Computational Phys., 71, 124-146 (1987)
  • Mapping the Sawtooth
    Chris Smiet, PPPL, abstract
    [#s1025, 02 May 2019]
    The magnetic field in a tokamak defines a field line map: a mapping from a poloidal cross-section to that same cross-section by following magnetic field lines. Such a map must necessarily contain fixed points, of which the magnetic axis is an example. The jacobian (the matrix of partial derivatives) $\mathsf{M}$ describes the mapping around such a fixed point to first order, and is part of the Lie group $SL(2,\mathbb{R})$. Different elements of this group act on the euclidean plane as rotations, shear mappings or hyperbolic fixed points. We look at a transition from an ellpitic fixed point (field lines lie on nested surfaces around the fixed point) to an alternating hyperbolic fixed point (field lines map to opposite branches of hyberbolic surfaces) that can occur at $q=2, 2/3, 2/5, 2/7 ...$. Using the NOVA-K code we identify an ideally unstable mode that is localized on the axis and has a high growth rate when the safety factor is close to 2/3. This mode drives the fixed point into the alternating hyperbolic regime. The nonlinear evolution of this instability can lead to complete stochastization of a region near the axis. Though the Sawtooth oscillation has long been observed, there is still disagreement between theoretical models and observations, and no model can reconcile all observations. We speculate that the above transition could explain some of the observations that do not fit other models, such as measurements of q below 1, snakes, and persistent Alfven Eigenmodes.
  • Hyper-Resistive Model of UHE Cosmic Ray Acceleration by AGNs
    Ken Fowler, UC Berkeley, abstract, slides
    [#s1024, 25 Apr 2019]
    Ultra High Energy (UHE) cosmic rays (~ 10^20 eV) may be produced by known processes of acceleration by plasma turbulence in magnetized jets produced by Active Galactic Nuclei (AGNs). A simple model in which turbulence is represented as hyper-resistivity in Ohm’s Law yields several predictions in sufficient agreement with observations to motivate further investigation. Besides jet dimensions, these predictions include the unique extra-galactic cosmic ray energy spectrum (\propto 1/E^3 ) and a different interpretation of the synchrotron radiation by which AGN jets are observed. Crucial to the model is a new theory of jet propagation whereby un-collimated jets generated by General Relativistic MHD simulations evolve to a highly collimated structure, finally evolving at speed 0.01c that explains jet dimensions, while relativistic acceleration parallel to field lines yields both cosmic rays and synchrotron radiation.

    References:
    [1] S. A. Colgate, T. K. Fowler, H. Li & J.Pino, 2014 ApJ 789, 144, on AGN jets
    [2] S. A. Colgate, T. K. Fowler, H. Li et al. 2015 ApJ 813, 136, on jet stability
    [3] T. K. Fowler & H. Li, 2016 J. Plas. Phys. 82, 595820513, on UHE acceleration
    [4] T. K. Fowler, H. Li, R. Anantua, 2019 ArXiV 2615445
  • Realistic 2D quasilinear modeling of fast ion relaxation
    Vinicius Duarte, PPPL, slides
    [#s1071, 19 Apr 2019]
  • Fishbone instability and transport of energetic particles
    Guillaume Brochard, CEA, IRFM, abstract, slides
    [#s1023, 18 Apr 2019]
    Analytical and numerical studies have been carried out to verify linearly the newly implemented 6D PIC module into the nonlinear hybrid code XTOR-K [1],[2]. This code solves the two-fluid extended MHD equations in toroidal geometry while taking into account, self-consistently, kinetic ion populations. The verification has been performed by a linear model [3], developed from [4],[5],[6]. It solves non-perturbatively the kinetic internal kink dispersion relation, with the particularity to take into account non-resonant kinetic terms and passing particles, which have revealed to be crucial features of the fishbone instability. A verification between the model and XTOR-K in its linear phase is presented, regarding the pulsation and growth rate (ω, γ) of the internal kink, and the position in phase space of resonant planes. As expected, the instability is stabilized on the kink branch and then destabilized on the fishbone branch. On the basis of this verification, a series of linear and nonlinear runs have been performed with XTOR-K. Firstly a linear study of the alpha-induced fishbone instability on the ITER 15 MA equilibrium has been done. It highlights the linear thresholds of the instability in the diagram (q0, βh,0), with q0 the on-axis safety factor and βh,0 the on-axis kinetic beta. The fishbone mode is found to be unstable for ITER relevant (q0, βh) couples. Secondly, a first nonlinear simulation has been performed to study the nonlinear evolution of the wave-particle interaction between the kink rotation and the alpha particle precessional motion on the fishbone branch. The equilibrium is taken to have circular ITER-like flux surfaces and alpha particles at peak energy of 1MeV. Such a simulation explores one limit described in [7] where particles non-linearities dominate, far above the linear threshold in βh,0. Results show strong chirping of the fishbone mode associated with the transport of resonant particles. It is found that, for this configuration, the fishbone instability transports 5% of the core’s alphas toward q = 1. Once the resonant particles have reached q = 1, a classical internal kink grows.

    References
    [1] H. L ̈utjens, J-F. Luciani, JCP, (2010), 229, 8130-8143
    [2] D. Leblond, PhD thesis, (2011)
    [3] G. Brochard et al 2018 J. Phys.: Conf. Ser. 1125 012003
    [4] D.Edery, X.Garbet, J-P.Roubin, A.Samain (1992), PPCF, 34, 6, 1089-1112
    [5] F. Porcelli et al, Phys. FLuid B 4 (10), 1992
    [6] F. Zonca, L. Chen, Physics of Plasmas 21, 072121 (2014) 1
    [7] F.Zonca et al, New J. Phys. 17 (2015) 013052
  • Improving computational efficiency of kinetic simulations with physics, mathematics, and machine learning
    George Wilkie, Centrum Wiskunde & Informatica (CWI), abstract, slides
    [#s1022, 11 Apr 2019]
    Kinetic theory has made tremendous progress in recent decades thanks to reduced models and improved computational capacity. Some problems, especially in the non-Maxwellian regime, remain difficult even for large supercomputing clusters. In this talk, I will discuss how such problems can be solved on laptops with the right tools applied. Physical approximations can be made to reduce the burden of predicting the interaction between turbulence and energetic particles. To complement well-established physical reductions of the nonlinear Boltzmann equation, recent advances in applied mathematics are utilized for direct efficient solution. Throughout, I will discuss ongoing and potential applications of machine learning to improve efficiency even further.
  • Simulating Relativistic Astroplasmas from Microphysics to Global Dynamics
    Kyle Parfrey, NASA Goddard Space Flight Center, abstract, slides
    [#s1007, 19 Mar 2019]
    The most extreme and surprising behaviors of black holes and neutron stars are driven by their surrounding magnetic fields and plasmas. Numerical simulations of these systems are complicated by the exotic physical conditions, requiring new approaches. I will present a range of computational methods which are well adapted to challenges such as strongly curved spacetime, energetically dominant electromagnetic fields, and pathological current sheets. In particular, I will describe how a new technique for general-relativistic plasma kinetics will aid in understanding black holes' particle acceleration and jet launching, and in interpreting future observations with the Earth-spanning Event Horizon Telescope.
  • Effects of zonal flows on transport crossphase in Dissipative Trapped-Electron Mode turbulence
    Michael Leconte, National Fusion Research Institute, Korea, abstract, slides
    [#s992, 14 Mar 2019]
    Confinement regimes with edge transport barriers occur through the suppression of turbulent (convective) fluxes in the particle and/or thermal channels, i.e. $Γ = \sum{\sqrt{|n_k|^2}\sqrt{|φ_k|^2} sin δ^{n,φ}_k}$ and $Q = \sum{\sqrt{|T_k|^2}\sqrt{|φ_k|^2} sin δ^{T,φ}_k}$, respectively, for drift-wave turbulence. The quantity $|φ_k|^2$ is the turbulence intensity, while $δ^{n,φ}_k$ and $δ^{T,φ}_k$ are the crossphases. For H-mode, a decorrelation theory predicts that it is the turbulence intensity $|φ_k|^2$ that is mainly affected via flow-induced shearing of turbulent eddies [1]. However, for other regimes (e.g. I-mode), characterized by high energy confinement but low particle confinement, this decrease of turbulence amplitude cannot explain the decoupling of particle v.s. thermal flux, since a suppression of turbulence intensity $|φ_k|^2$ would necessarily affect both fluxes the same way. Here, we explore a possible new stabilizing mechanism: zonal flows may directly affect the transport crossphase. We show the effect of this novel mechanism on the turbulent particle flux, by using a simple fluid model [2,3] for dissipative trapped-electron mode (DTEM), including zonal flows. We first derive the evolution equation for the transport crossphase $δ_k$ between density and potential fluctuations, including contributions from the E × B nonlinearity [4]. By using a parametric interaction analysis including the back-reaction on the pump, we obtain a predator-prey like system of equations for the pump amplitude $φ_p$, the pump crossphase $δ_p$, the zonal amplitude $φ_z$ and the triad phase-mismatch ∆δ. The system displays limit-cycle oscillations where the instantaneous DTEM growth rate - proportional to the crossphase - shows quasi-periodic relaxations where it departs from that predicted by a linear theory. This implies that the crossphase does not respond instantaneously to the driving gradient, instead there is a finite response time, which for DTEM corresponds to the inverse of the de-trapping rate $ν = ν_{ei}/\sqrt{\epsilon}$ with $ν_{ei}$ the electron-ion collisionality and $\epsilon$ the inverse aspect ratio.

    [1] H. Biglari, P.H. Diamond, P.W. Terry, Phys. Fluids B 2, 1 (1990).
    [2] D. A. Baver, P.W. Terry, R. Gatto and E. Fernandez, Phys. Plasmas 9, 3318 (2002).
    [3] F.Y. Gang, P.H. Diamond, J.A. Crotinger and A.E. Koniges, Phys. Fluids 3, 955 (1991).
    [4] C.Y. An, B. Min and C.B. Kim, Plasma Phys. Controlled Fusion 59, 115006 (2017).
  • Astrophysical Collisionless Shock and Current Sheet Instabilities: Particle Modeling and Laboratory Study
    Zhenyu Wang, Princeton University, abstract, slides
    [#s993, 13 Mar 2019]
    In the first part of this talk, I will present the modelling and interpretation of a campaign of laser experiments designed to generate high Mach number magnetized collisionless shocks on OMEGA-EP. We compare the data to the results of 3-D PIC simulations, and describe the signatures of magnetized shock formation, including the early contact discontinuity formation stage, and a later magnetic reflection with magnetic overshoots. We explain the geometrical effects on the radiography introduced by density gradient in expanding plasma and by the curvature of the imposed magnetic field. We conclude that our experiments have reproducibly achieved magnetized shocks with Alfvenic Mach number 3 to 9 in laboratory conditions. In the second part, I will describe the gyrokinetic (GK) electron and fully kinetic (FK) ion particle (GeFi) simulation model and the particle simulation results of waves and current sheet instabilities. In the GeFi model, the GK electron approximation removes the high frequency electron gyromotion and plasma oscillation, but the electron finite Lamor radii effects are retained. For lower-hybrid waves, the GeFi results agree well with the fully kinetic explicit delta-f particle code and the fully kinetic Darwin particle code. Our 3-D GeFi and FK simulation results demonstrate the existence of the lower-hybrid-drift, kink and sausage instability in current sheet under finite guide magnetic fields with the realistic proton-to-electron mass ratio.
  • Kinetic Simulations of Collisionless Plasmas
    Rahul Kumar, Princeton University, abstract, slides
    [#s989, 05 Mar 2019]
    The particle-in-cell method has been remarkably successful in understanding physical processes occurring at the kinetic scales. I will discuss the implementation of the electromagnetic particle-in-cell method for collisionless plasmas in a self-developed code called PICTOR. I will focus on a few techniques employed to improve performance, diagnostic, and scalability of the code. I will then briefly discuss two physics problems to illustrate the efficacy of PIC simulations in addressing a few outstanding problems in plasma physics. First, I will discuss preferential heating and acceleration of heavy ion in Alfvenic turbulence. Then, I'll discuss how self-consistent PIC simulations, combined with the measurements from Voyager spacecraft, could be used to obtain a comprehensive understanding of the dynamics of solar wind termination shock.
  • Available energy of magnetically confined plasmas
    Per Helander, Max Planck Institute for Plasma Physics, Greifswald, Germany, abstract, slides
    [#s990, 28 Feb 2019]
    In this talk, the energy budget of a collisionless plasma subject to electrostatic fluctuations is studied. In particular, the excess of thermal energy over the minimum accessible to it under various constraints that limit the possible forms of plasma motion is considered. This excess measures how much thermal energy is “available” for conversion into plasma instabilities, and therefore constitutes a nonlinear measure of plasma stability. The “available energy” defined in this way becomes an interesting and useful quantity in situations where adiabatic invariants impose non-trivial constraints on the plasma motion. For instance, microstability properties of certain stellarators can be inferred directly from the available energy, without the need to solve the gyrokinetic equation. The technique also suggests that an electron-positron plasma confined by a dipole magnetic field could be entirely free from turbulence.
  • Finite Larmor Radius effects at the high confinement mode pedestal and the related force-free steady state
    Wei-li Lee, PPPL, abstract, slides
    [#s988, 19 Feb 2019]
    For this talk, we will first relate our previous calculations on the radial electric field at the high confinement H-mode pedestal [W. W. Lee and R. B. White, Phys. Plasmas 24, 081204 (2017)] with the actual magnetic fusion experimental measurements. We will then discuss the new pressure balance due to the E x B current, which is induced by the resulting radial electric field, and its impact on the gyrokinetic MHD equations as well as their conservation properties in the force-free steady state.
  • Impurity transport in stellarator plasmas
    Albert Mollén, IPP Greifswald, abstract, slides
    [#s985, 13 Feb 2019]
    In contrast to tokamaks where turbulence typically dominates, a substantial fraction of the radial energy and particle transport in stellarators can often be attributed to collisional processes. The kinetic calculation of collisional transport has for a long time relied on simplified models which use the “mono-energetic” approximation, the simple pitch-angle scattering collision operator and are radially local. But not all experimental observations have been satisfactorily explained, and in recent years more advanced numerical tools have appeared which relax some of the approximations. These improvements in the modelling can be of particular importance when analyzing the impurity transport. I will discuss how the calculation of the impurity transport has advanced in recent years, and what the latest observations in the Wendelstein 7-X stellarator are.
  • Black Aurora, Bohm Diffusions, Quasi-Neutrality Mysteries Explained
    Kwan Chul Lee, NFRI, abstract, slides
    [#s949, 07 Feb 2019]
    Gyrocenter-Shift Analysis has been developed based on the momentum exchange of ion-neutral interaction in the magnetized plasmas. Recent 7 years, after the last PPPL theory seminar held on January 10, 2012, 3 papers on GCS analysis have been published [1-3]. This presentation will focus on the quasi-neutrality controversy after revisiting the basis of GCS analysis. The explanations on the satellite measurement of the black aurora and the transport phenomena named as Bohm diffusions will be presented as well as a related KSTAR experimental proposal.
    [1] K. C. Lee, “Violation of Quasi-neutrality for Ion-neutral Charge-exchange Reactions in Magnetized Plasmas”, J. of Korean Phys. Soc., Vol. 63, No.10, 1944 (2013)
    [2] Kwan Chul Lee, “Analysis of Bohm Diffusions Based on the Ion-Neutral Collisions”, IEEE Trans. on Plasma Science, Vol. 43 No. 2, 494 (2015)
    [3] Kwan Chul Lee, “Electric field formation in three different plasmas: A fusion reactor, arc discharge, and the ionosphere”, Phys. of Plasmas, Vol. 24, 112505 (2017)
  • Bringing the golden standard into the silicon age
    Francesca Poli, PPPL, abstract, slides
    [#s975, 31 Jan 2019]
    A lot happened since TRANSP made its first appearance about thirty years ago. The 'golden standard' for tokamak discharge analysis has evolved to a code that is capable of predicting heating and current drive, thermal and particle transport. Its pool of users has expanded internationally to cover almost every tokamak operating nowadays, demanding modernization, increasing support and additional capabilities. This talk will review the strength and weaknesses of TRANSP, the plans for implementation of new physics modules targeting (as close as possible to) a Whole Device Model. It will discuss areas for partnership with the theory department and ongoing activities and plans in collaboration with the SciDAC projects, for extension of the transport outside the plasma boundary and for self-consistent calculations of transport and stability, including transport induced by energetic particles.
  • Compressional Alfvén eigenmodes excited by runaway electrons
    Chang Liu, PPPL, abstract, slides
    [#s960, 25 Jan 2019]
    Runaway electron phenomena is an important topic in general, and critically important in tokamak disruption studies. Given its high potential for damaging effects, it is critical to understand the physics of their generation and find a mitigation strategy for ITER. Kinetic instabilities associated with high-energy runaway electrons have been shown to play an important role in the behavior of runaway electrons Recently, a new kind of instability with magnetic signals in the Alfvén frequency range have been observed in disruption experiments in DIII-D, which is found to be related to the dissipation of the runaway electron current. In this talk, a candidate explanation for this phenomena is presented, namely resonant interaction with compressional Alfvén eigenmodes (CAE). CAEs driven by energetic ions have been well studied in spherical tokamaks like NSTX. For CAEs driven by resonances with runaway electrons, the damping rate of the modes due to electron-ion collisions are calculated. The model is applied to a time-dependent kinetic simulation of runaway electrons, which includes the bounce-average effect and the enhanced ion pitch-angle scattering due to partial screening. The results match with experiments qualitatively, and provide a promising way to diffuse runaway electrons before their energization. A brief overview of related research into runaway electrons is also given, indicating how this work fits into the wider effort to find mitigation strategies.
  • Magnetohydrodynamical equilibria with current singularities and continuous rotational transform
    Yao Zhou, abstract, slides
    [#s967, 23 Jan 2019]
    We revisit the Hahm-Kulsrud-Taylor (HKT) problem, a classic prototype problem for studying resonant magnetic perturbations and 3D magnetohydrodynamical equilibria. We employ the boundary-layer techniques developed by Rosenbluth, Dagazian, and Rutherford (RDR) for the internal m=1 kink instability, while addressing the subtle difference in the matching procedure for the HKT problem. Pedagogically, the essence of RDR's approach becomes more transparent in the reduced slab geometry of the HKT problem. We then compare the boundary-layer solution, which yields a current singularity at the resonant surface, to the numerical solution obtained using a flux-preserving Grad-Shafranov solver. The remarkable agreement between the solutions demonstrates the validity and universality of RDR's approach. In addition, we show that RDR's approach consistently preserves the rotational transform, which hence stays continuous, contrary to a recent claim that RDR's solution contains a discontinuity in the rotational transform.
  • Gyrokinetic continuum simulations of plasma turbulence in the Texas Helimak
    Tess Bernard, University of Texas - Austin, abstract
    [#s942, 13 Dec 2018]
    The first gyrokinetic simulations of plasma turbulence in the Texas Helimak device are presented. These have been performed using the Gkeyll (http://gkyl.readthedocs.io/) computational framework. The Helimak is a simple magnetized torus with a toroidal and vertical magnetic field and open field lines that terminate on conducting plates at the top and bottom of the device. It has features similar to the scrape-off layer region of tokamaks, such as bad curvature-driven instabilities and sheath boundary conditions on the end plates, which are included in these initial gyrokinetic simulations. A bias voltage can be applied across conducting plates to drive E x B flow and study the effect of velocity shear on turbulence suppression. Comparisons are presented between simulations and measurements from the experiment, showing qualitative similarities, including fluctuation amplitudes and equilibrium profiles that approach experimental values. There are also some important quantitative differences, and I discuss how certain physical and geometric effects may improve agreement in future results.
  • Machine Learning Techniques for Analysis of High-Dimensional Chaotic Spatiotemporal Dynamics
    Jaideep Pathak, University of Maryland, abstract, slides
    [#s927, 06 Dec 2018]
    High-dimensional chaos is a commonly observed feature of many interesting natural systems such as fluid flows or atmospheric dynamics. We formulate a machine learning approach for short-term prediction [1] and for understanding the long-term ergodic properties of the underlying dynamical processes of high-dimensional chaotic spatiotemporal dynamical systems. Specifically, we consider a situation in which limited duration time series data from some dynamical process is available, but a mechanistic, knowledge-based model of how that data is produced is either unavailable or too inaccurate to be useful. Our proposed approach, using a particular kind of recurrent neural network called a Reservoir Computer [2], is computationally efficient and scalable through parallelization to dynamical systems with arbitrarily high-dimensional chaotic attractors. Further, in the context of chaotic dynamical systems, machine learning can also be used in conjunction with an imperfect knowledge-based model for filling in gaps in our underlying mechanistic knowledge that can cause the model to be inaccurate. We demonstrate this technique using simulated data from the Kuramoto-Sivashinsky partial differential equation as well as the Lorenz ‘96 toy model of atmospheric dynamics [3].
    [1] J. Pathak, B. Hunt, M. Girvan, Z. Lu, E. Ott, Phys. Rev. Lett. 120 (2018) 024102.
    [2] H. Jaeger, H. Haa, Science 304 (2004) 78–80.
    [3] E.N. Lorenz, in:, T. Palmer, R. Hagedorn, Predict. Weather Clim., Cambridge University Press, Cambridge, 2006, pp. 40–58.
  • From SOL turbulence to planetary magnetospheres: computational plasma physics at (almost) all scales using the Gkeyll code
    Ammar Hakim, PPPL, abstract, slides
    [#s928, 05 Dec 2018]
    Gkeyll is a computational plasma physics package that aims to simulate plasmas at all scales. At present, the code contains solvers for three major equation systems: Vlasov-Maxwell equations, electromagnetic gyrokinetic equations and multi-fluid moment equations. These span the complete range of plasma physics; electromagnetic shocks, turbulence and first-principles sheath physics, requiring full kinetic treatment; turbulence in tokamak core and SOL, requiring EM gyrokinetics; and planetary magnetospheres, requiring fluid treatment with proper accounting of kinetic effects to capture reconnection and current sheet dynamics. In this talk I will present the status of each of the major solvers implemented in the code, with emphasis on the features of the algorithms as well as the physics being studied. In particular, I will focus on recent progress in implementing a conservative algorithm for collisions; details of a novel algorithm for EM gyrokinetic in the symplectic formulation; and recent progress in developing a robust semi-implicit algorithm for multi-fluid moment equations. I will conclude with the short- and medium-term plans for the project as well as some ideas on strengthening advanced computing at PPPL, leveraging on the work done in Gkeyll as well as other projects.
    https://gkeyll.readthedocs.io/en/latest/
  • Coupled core-edge gyrokinetic simulation and Tungsten impurity transport in JET with XGC
    Julien Dominski, abstract, slides
    [#s929, 04 Dec 2018]
    Incoming Exascale capabilities of super computers will enable whole-device simulations based on first principles plasma physics. To take full advantage of these new capabilities, new numerical schemes and more complete physical models are developed in XGC. XGC is a whole-volume total-f gyrokinetic code optimized for simulation of edge plasma in magnetic fusion devices. One goal of the ECP project is to couple XGC with a core code, such as GENE, to optimize the efficiency of whole-device simulations. The current status of this core-edge coupling project will be presented, including the presentation of the core-edge coupling scheme and of the coupled GENE-XGC simulations. Another goal is to study the influence of tungsten and beryllium on the performance of ITER. The Total-f gyrokinetic XGC is thus being upgraded to simulate the physics of many-species impurities in the whole-volume including SOL. First multi-species simulation of a JET plasma under tungsten contamination will be demonstrated, showing that the lower charge state W can move inward from into core, but the higher charge state W will move outward toward the pedestal top and accumulate, as seen in ASDEX-U.
  • KBM nonlinear dynamics and first-principles-based classifiers with Machine Learning predictions for disruptions
    Ge Dong, Princeton University, abstract, slides
    [#s923, 20 Nov 2018]
    Kinetic ballooning modes (KBM) are widely believed to play a critical role in disruptive dynamics as well as turbulent transport in tokamaks. While the nonlinear evolution of ballooning modes has been proposed as a mechanism for “detonation” in tokamak plasmas, the role of kinetic effects in such nonlinear dynamics remains largely unexplored. In this study saturation mechanism and nonlinear dynamics of KBM is presented, with global gyrokinetic simulation results of KBM nonlinear behavior. Instead of the finite-time singularity predicted by ideal MHD theory, the kinetic instability is shown to develop into an intermediate non-linear regime of exponential growth, followed by a nonlinear saturation regulated by spontaneously generated zonal fields. In the intermediate nonlinear regime, rapid growth of localized current sheet, which can mediate reconnection, is observed. In the future, the linear properties as well as nonlinear mode structures from the simulations could be incorporated into the deep learning models for disruption predictions in the form of a new parameter/channel, as a first-principles physics guide to the AI. The deep learning model could in turn provide feedback on the sensitivity of the parameters, including the linear stability properties of various modes, and nonlinear dynamics of these instabilities, and thus automatically select new inputs for the first-principles codes.
  • Electron-impact excitation of molecular hydrogen: dissociation and vibrationally resolved cross sections
    Dmitry Fursa, Curtin University, abstract
    [#s922, 16 Nov 2018]
    Molecular hydrogen and its isotopologues are present in a range of vibrationally excited states in fusion, atmospheric, and interstellar plasmas. Electron-impact excitation cross sections resolved in both final and initial vibrational levels of the target are required for modeling the properties and dynamics, and controlling the conditions of many low-temperature plasmas. Recently, the convergent close-coupling (CCC) method has been utilized to provide a comprehensive set of accurate excitation, ionization, and grand total cross sections for electrons scattering on H2 in the ground (electronic and vibrational) state, and calculations are being conducted to extend this data set to include cross sections resolved in all initial and final vibrational levels. In this talk I will review the available e-H2 collision data, discuss the resolution of a significant discrepancy between theory and experiment for excitation of the b3Su+ state, and present estimates for dissociation of H2.
  • Strong-flow gyrokinetics with a unified treatment of all length scales
    Amil Sharma, University of Warwick, abstract, slides
    [#s921, 14 Nov 2018]
    Tokamak turbulence exhibits interaction on all length scales, but standard gyrokinetic treatments consider global scale flows and gyroscale flows separately, and assume a separation between these length scales. However, the use of a small-vorticity ordering (Dimits, 2010) allows for the presence of large, time-varying flows on large length scales, whilst providing a unified treatment including shorter length scales near and below the gyroradius. Some examples of strong-flow generalisations of gyrokinetics are presented, followed by a description of the nuances of the equations and numerical implementation that use the ordering of Dimits (2010). Our Euler-Lagrange and Poisson equations contain an implicit dependence that appears as a partial time derivative of the E × B flow. This implicit dependence is analogous to the v||-formulation of gyrokinetics. However, as these implicit terms are small, we use an iterative scheme to resolve this. Additionally, we have developed a stand-alone Poisson solver based on that from the ORB5 code, and use this to simulate certain flow and density gradient driven instabilities in cylindrical geometry.
  • Some Novel Features of Three Dimensional MagnetoHydroDynamic Plasma
    Rupak Mukherjee, Institute for Plasma Research, abstract, slides
    [#s920, 12 Nov 2018]
    Within the framework of MagnetoHydroDynamics, a strong interplay exists between flow and magnetic fields leading to several interesting pathways for energy cascade. In this talk I numerically demonstrate three examples of such interplay using our in-house developed DNS code G-MHD3D which simulates three dimensional single fluid MHD equations. I also suggest analytical arguments in some of our numerical observations. The first problem discusses the phenomena of nonlinear interaction of magnetic and kinetic energies within the premise of two and three dimensional MHD equations leading to periodic exchange of energy. Scaling of the energy exchange frequency with deviation from Alfven resonance and initial wave number of excitation is numerically determined. Results are qualitatively reproduced by analysing the set of single fluid MHD equations through low degrees of freedom coupled ODEs obtained via a Galerkin procedure. Secondly, in three dimensions, at Alfven resonance, for some chaotic flows, the initial flow profile is found to “recur” periodically with the time evolution of the plasma. Such recurrence is unexpected in systems with high degree of freedom (e.g. 3D MHD). The primary cause of such phenomena is analysed using an effective number of active degrees of freedom present in the system. Finally we observe some preliminary results of large and intermediate scale magnetic field generation in plasmas often called as “dynamo” using our code. The growth rate of such ‘fast’ dynamos are compared for different parameters of the system and the fastest dynamos for the parameter set is identified. Physics details and numerical aspects of the development of the code, numerical protocols followed, direct numerical simulations results, numerical tools to diagnose the three dimensional grid data and the analytical arguments in support of the numerical observations will be presented in the talk.
  • A gyrokinetic model for the tokamak periphery
    Rogerio Jorge, EPFL, abstract, slides
    [#s894, 01 Nov 2018]
    We present a new gyrokinetic model that retains the fundamental elements of the plasma dynamics at the tokamak periphery, namely electromagnetic fluctuations at all scales, comparable amplitudes of background and fluctuating components, and a large range of collisionality regimes. Such model is derived within a gyrokinetic full-F approach, describing distribution functions arbitrarily far from equilibrium, and projecting the gyrokinetic equation onto a Hermite-Laguerre velocity space polynomial basis, obtaining a gyrokinetic moment hierarchy. The treatment of arbitrary collisionalities is performed by expressing the full Coulomb collision operator in gyrocentre phase space coordinates, and providing a closed formula for its gyroaverage in terms of the gyrokinetic moments. In the electrostatic regime and long-wavelength limit, the novel gyrokinetic hierarchy reduces to a drift-kinetic moment hierarchy that in the high collisionality regime further reduces to an improved set of drift-reduced Braginskii equations, which are widely used in scrape-off layer simulations. First insights on the linear modes described by our novel gyrokinetic model will be presented.
  • APS Invited dry run - Quasi-linear resonance broadened model for fast ion relaxation in the presence of Alfvénic instabilities
    Nikolai Gorelenkov, PPPL, abstract, slides
    [#s883, 25 Oct 2018]
    The resonance broadened quasi-linear (RBQ) model is developed for the problem of relaxing the fast energetic particle distribution function in constant-of-motion (COM) 3D space [N.N. Gorelenkov et al., Nucl. Fusion 58 (2018) 082016]. The model is generalized by using the QL theory [H. Berk et al., Phys. Plasmas'96] and carefully reexamining the wave particle interaction (WPI) in the presence of realistic AE mode structures and pitch angle scattering with the help of the guiding center code ORBIT. The RBQ model applied in realistic plasma conditions is improved by going beyond the perturbative-pendulum-like approximation for the wave particle dynamics near the resonance. The resonance region is broadened but remains 2-3 times smaller than predicted by an earlier bump-on-tail QL model. In addition the resonance broadening includes the Coulomb collisional or anomalous pitch angle scattering. The RBQ code takes into account the beam ion diffusion in the direction of the canonical toroidal momentum. The wave particle interaction is reduced to one-dimensional dynamics where for the Alfvénic modes typically the particle kinetic energy is nearly constant. The diffusion equation is solved simultaneously for all particles together with the evolution equation for the mode amplitudes. We apply the RBQ code to a DIII-D plasma with elevated q -profile where the beam ion profiles show stiff transport properties [C. Collins et al. PRL'16]. The sources and sinks are included via the Krook operator. The properties of AE driven fast ion distribution relaxation are studied for validations of the applied QL model to DIII-D discharges. Initial results show that the model is robust, numerically efficient, and can predict intermittent fast ion relaxation in present and future burning plasmas.
  • Nonlinear interaction between toroidal Alfvén eigenmodes and tearing modes
    Z. W. Ma, Zhejiang University, abstract, slides
    [#s809, 15 Aug 2018]
    Nonlinear interaction between toroidal Alfvén eigenmodes (TAEs) and the tearing mode is investigated by using the hybrid code CLT-K. It is found that the $n=1$ TAE is first excited by isotropic energetic particles at the linear stage and reaches the first steady state due to wave- particle interaction. After the saturation of the $n=1$ TAE, the $m/n=2/1$ tearing mode grows continuously and reaches its steady state due to nonlinear mode-mode coupling, especially, the $n=0$ component plays a very important role in the tearing mode saturation. The results suggest that the enhancement of the tearing mode activity with increase of the resistivity could weaken the TAE frequency chirping through the interaction between the $p=1$ TAE resonance and the $p=2$ tearing mode resonance for passing particles in the phase space, which is opposite to the classical physical picture of the TAE frequency chirping that is enhanced with increase of dissipation.
  • Understanding coronal structures on the Sun
    Hardi Peter, Max Planck Institute for Solar System Research, Germany, abstract, slides
    [#s799, 17 Jul 2018]
    Since decades coronal heating is a buzzword that is used as a motivation on coronal research. Depending on the level of detail one is interested in, one could define this question anything ranging from answered to not understood at all. 3D MHD models can now produce a corona in a numerical experiment that comes close to the real Sun in complexity. And the fact alone that in these models a three-dimensional loop-dominated time-variable corona is produced could be used as an argument that the problem of coronal heating is solved. However, careful inspection of these model results shows that despite their success they leave many fundamental questions unanswered. In this talk I will address some of these aspects, including the mass and energy exchange between chromosphere and corona, the apparent width of coronal loops, the energy source of hot active region core loops, or the internal structure of loops. In this sense this talk will pose more questions that it provide answers.
  • Machine-learning driven correlation studies: multi-band frequency chirping at NSTX
    Ben Woods, University of York, abstract, slides
    [#s777, 28 Jun 2018]
    Magnetic perturbations in a very broadband range (<30 kHz to >1 GHz) are commonly measured on tokamaks such as NSTX by using Mirnov coils. The spectral behaviour of the perturbations can be categorised as quiescent, fixed-frequency, chirping, or avalanching. Here, ‘chirping’ modes experience a time-dependent frequency shift due to non-linear effects – in some cases, multiple plasma modes chirp in a near-concurrent fashion (mode 'avalanching'). Mode avalanching in the Alfvénic and super-Alfvénic frequency bands is typically correlated with fast-ion loss. However, transition to this phase of mode behaviour is not fully understood. Traditional methods for characterising mode behaviour are highly labour intensive for human characterisation - studying parametric dependences of plasma parameters on mode character proves difficult. Here, preliminary results are presented from machine-learning driven studies of correlations between different mode character (quiescent, fixed-frequency, chirping, avalanching) and weighted averages of plasma parameters obtained from TRANSP (v-fast/v-Alfvén, Q-profile, β-fast/β-Alfvén). The weighted averages form insightful metrics of the stability of plasma modes, allowing for correlations to be drawn (i.e. magnetic shear versus mode character). These results yield similar correlations to previous work by Fredrickson et al. [1] in both the kink/tearing/fishbone frequency band (~1-30 kHz), and the TAE band (~50-200 kHz). An overall framework is presented to utilise this tool for generic tokamaks, for possible future use on MAST-U and DIII-D.
    [1] E. D. Fredrickson, N.N. Gorelenkov et al., Nucl. Fusion 54, 093007 (2014)
  • Conservative Discontinuous Galerkin Discretization for the Landau Collision Operator
    Alexander Frank, Technical University of Munich, abstract, slides
    [#s773, 21 Jun 2018]
    This talk presents a mass-, momentum- and energy conserving discretization of the nonlinear Landau collision integral. The semi-discrete form is achieved using a modal discontinuous Galerkin method on a tensor product mesh and a recovery method to define second derivatives at the element boundaries. Combined with an explicit time stepping scheme this gives a fully discrete conservative form. The conservation properties are proven algebraically and shown numerically for a two dimensional relaxation test problem.
  • Computing local sensitivity and tolerances for stellarators using shape gradients
    Matt Landreman, University of Maryland, abstract, slides
    [#s741, 14 Jun 2018]
    Tight tolerances have been a leading driver of cost in recent stellarator experiments, so improved definition and control of tolerances can have significant impact on progress in the field. Here we relate tolerances to the shape gradient representation that has been useful for shape optimization in industry, used for example to determine which regions of a car or aerofoil most affect drag, and we demonstrate how the shape gradient can be computed for physics properties of toroidal plasmas. The shape gradient gives the local differential contribution to some scalar figure of merit (shape functional) caused by normal displacement of the shape. In contrast to derivatives with respect to quantities parameterizing a shape (e.g. Fourier amplitudes), which have been used previously for optimizing plasma and coil shapes, the shape gradient gives spatially local information and so is more easily related to engineering constraints. We present a method to determine the shape gradient for any figure of merit using the parameter derivatives that are already routinely computed for stellarator optimization, by solving a small linear system relating shape parameter changes to normal displacement. Examples of shape gradients for plasma and electromagnetic coil shapes are given. We also derive and present examples of an analogous representation of the local sensitivity to magnetic field errors; this magnetic sensitivity can be rapidly computed from the shape gradient. The shape gradient and magnetic sensitivity can both be converted into local tolerances, which inform how accurately the coils should be built and positioned, where trim coils and structural supports for coils should be placed, and where magnetic material and current leads can best be located. Both sensitivity measures provide insight into shape optimization, enable systematic calculation of tolerances, and connect physics optimization to engineering criteria that are more easily specified in real space than in Fourier space.
  • Energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator with exact H-theorems
    Eero Hirvijoki, PPPL, abstract, slides
    [#s737, 31 May 2018]
    This talk explores energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator. We discuss two approaches, namely direct Galerkin formulations and discretizations of the underlying infinite-dimensional metriplectic structure of the collision integral. The spatial discretizations are chosen to reproduce the time-continuous conservation laws that correspond to Casimir invariants and to guarantee the positivity of the distribution function. Both the direct and the metriplectic discretization are demonstrated to have exact H-theorems and unique, physically exact equilibrium states. Most importantly, the two approaches are shown to coincide, given the chosen Galerkin method. A temporal discretization, preserving all of the mentioned properties, is achieved with so-called discrete gradients. The proposed algorithm successfully translates all properties of the infinite-dimensional time-continuous Landau collision operator to time- and space-discrete sparse-matrix equations suitable for numerical simulation.
  • Spontaneous reconnection in thin current sheet: the “ideal” tearing mode
    Anna Tenerani, Department of Earth, Planetary and Space Sciences, UCLA, abstract, slides
    [#s734, 24 May 2018]
    Magnetic field reconnection is considered one of the most important mechanisms of magnetic energy conversion and reorganization acting in astrophysical and laboratory plasmas. Although our knowledge has been greatly advanced in the last few decades, the problem of how magnetic reconnection can be triggered explosively in weakly collisional (quasi-ideal) plasmas (as observed e.g in solar flares, geomagnetic substorms and sawtooth crashes in tokamaks) still remains an open field of research. Here we discuss a possible scenario for the triggering of explosive reconnection via the onset of an ‘ideal’ tearing instability within forming current sheets and we show results from MHD numerical simulations that support our scenario. We demonstrate that the same reasoning, if applied recursively, can describe the complete nonlinear disruption of the original current sheet until microscopic marginally-stable current sheets are formed. We show that the ‘ideal’ tearing mode provides a general frame of work that can be extended to include other effects such as kinetic effects and different current profiles.
  • Design of a new quasi-axisymmetric stellarator equilibrium
    Sophia Henneberg, IPP Greifswald, abstract, slides
    [#s730, 17 May 2018]
    A new quasi-axisymmetric, two-field-period stellarator configuration has been designed following a broad study using the optimization code ROSE (Rose Optimizes Stellarator Equilibria). Because of the toroidal symmetry of the magnetic field strength, quasi-axisymmetric stellarators share many neoclassical properties of tokamaks, such as a comparable bootstrap current which can be employed to simplify the coil structure, which is favorable for finding compact equilibria. The ROSE code optimizes the plasma boundary calculated with VMEC based on a set of physical and engineering criteria. Various aspect ratios, number of field periods and iota profiles are investigated. As an evaluation of the design, the bootstrap current, the ideal MHD stability, the fast-particle losses, and the existence of islands are examined. The main result of this extensive study – a compact, MHD-stable, two-field-period stellarator with small fast-particle loss fraction – will be presented.
  • Bringing the GBS plasma turbulent simulation code from limited to diverted configurations
    Paola Paruta, École Polytechnique Fédérale de Lausanne, abstract
    [#s729, 17 May 2018]
    The GBS code has been developed to describe the plasma turbulent behaviour in the SOL,[1] [2], by solving the two-fluids Drift Reduced Braginskii equations.[3] We report on the implementation of diverted magnetic equilibria in GBS: by abandoning flux coordinates systems, which are not defined at the X-point, the model equations are written in toroidal coordinates, and a 4th order finite difference scheme is used for the implementation of the spatial operators on staggered poloidal and toroidal grids. The GBS numerical implementation is verified through the Method of Manufactured Solutions.[4] Its convergence properties are tested. First results of TCV-like simulations are presented.
    [1] F. D. Halpern, P. Ricci et al., J. Comput. Phys. 315, 388 (2016)
    [2] P. Ricci, F. D. Halpern et al., Plasma Phys. Control. Fusion 54, 124047 (2012)
    [3] A. Zeiler, J. F. Drake & B. Rogers, Phys. Plasmas 4, 2134 (1997)
    [4] Patrick J. Roache, J. Fluids Eng. 124, 4 (2002)
  • Modeling Resonant Field Penetration and Its Effects on Transport in the DIII-D Tokamak
    Qiming Hu, PPPL, abstract, slides
    [#s665, 05 Apr 2018]
    We use the nonlinear cylindrical two fluid MHD code (TM1) [1-4] to model the effects of multiple magnetic islands in the DIII-D tokamak. The TM1 code has been previously used to study classical and neoclassical tearing mode (TM) stability, TM stabilization by ECCD, and plasma response (transport and stability) to resonant magnetic perturbations (RMPs). Recently, TM1 has been used to understand the effects of multiple locked magnetic islands on heat transport in DIII-D Ohmic plasmas. It is found that co-existence of 2/1, 3/1 and 4/1 locked islands can produce a large (~50%) reduction in the central electron temperature, even without island overlap. However, the observed reduction in the edge temperature requires island overlap and stochasticity within the TM1. X-ray imaging reveals the appearance of multiple locked islands in the plasma edge at the time of the thermal collapse [5], consistent with TM1 modeling. For ELM suppression studies, we analyze the nonlinear evolution of edge magnetic islands in response to resonant magnetic perturbations [6]. The observed increase in pedestal toroidal rotation with the decrease in the core toroidal rotation is shown to be quantitatively consistent with TM1 modeling of multiple locked non-overlapping magnetic islands in the plasma edge. The TM1 results reveal interesting physics effects of locked modes that motivate further study using more comprehensive physics models.
    [1] Qingquan Yu, Phys. Plasmas 4, 1047 (1997)
    [2] Q. Yu, S. Günter & B. D. Scott, Phys. Plasmas 10, 797 (2003)
    [3] Q. Yu, Nucl. Fusion 50, 025014 (2010)
    [4] S. Günter, Q. Yu et al., J. Comput. Phys. 209, 354 (2005)
    [5] X. D. Du et al., Phys. Rev . Lett. (submitted)
    [6] R. Nazikian, C. Paz-Soldan et al., Phys. Rev. Lett. 114, 105002 (2015)
  • Two etudes on unexpected behaviour of drift-wave turbulence near stability threshold
    Alex Schekochihin, University of Oxford, abstract, slides
    [#s664, 04 Apr 2018]
    I will discuss some recent results — numerical and experimental — on the nature of drift-wave turbulence in MAST, obtained in the doctoral theses of my students Ferdinand van Wyk [1,4], Michael J. Fox [2] and Greg Colyer [3]. At ion scales, in the presence of flow shear, we find numerically [1,4] a type of transition to turbulence that is new (as far as we know) in the tokamaks, but reminiscent of some fluid dynamical phenomena (e.g., pipe flows or accretion discs in astrophysics): close to threshold, the nonlinear saturated state and the associated anomalous heat transport are dominated by long-lived coherent structures, which drift across the domain, have finite amplitudes, but are not volume filling; as the system is taken away from the threshold into the more unstable regime, the number of these structures increases until they overlap and a more conventional chaotic state emerges. Such a transition has its roots in the subcritical nature of the turbulence in the presence of flow shear. It can be diagnosed in terms of the breaking of the statistical up-down symmetry of turbulence: this manifests itself in the form of tilted two-point correlation functions and skewed distributions of the fluctuating density field, found both in simulations and in BES-measured density fields in MAST [2]. The governing (order) parameter in the system is the distance from the threshold, rather than individual values of equilibrium gradients; the symmetries — and drift-wave/zonal-flow turbulence of conventional type — are restored away from the threshold. The experiment appears to lie just at the edge of this latter transition rather than at the exact stability threshold. At electron scales in MAST, the conventional streamer-dominated state of ETG turbulence turns out to be a long-time transient, during which an initially unimportant zonal component continues to grow slowly, eventually leading to a new saturated state dominated by zonal modes, rather similar to ITG turbulence [3]. In this regime, the heat flux turns out to be proportional to the collision rate, in approximate agreement with the experimentally observed collisionality scaling of the energy confinement in MAST. Our explanation of this effect is based on a model of ETG turbulence dominated by zonal–nonzonal interactions and on an analytically derived scaling of the zonal-mode damping rate with the electron–ion collisionality. These developments open some intriguing possibilities both for enterprising theoreticians tired of the V&V routine and for ingenious experimentalists interested in making use of tokamaks to probe transitions to turbulence in nonlinear plasma systems.
    [1] F. van Wyk, E. G. Hichcock et al., J. Plasma Phys. 82, 905820609 (2016)
    [2] M. F. J. Fox, F. van Wyk et al., Plasma Phys. Control. Fusion 59, 034002 (2017)
    [3] G. J. Colyer, A. A. Schekochihin et al., Plasma Phys. Control. Fusion 59, 055002 (2017)
    [4] F. van Wyk, E. G. Highcock et al., Plasma Phys. Control. Fusion 59, 114003 (2017)
  • An adjoint method for gradient-based optimization of stellarator coil shapes
    Elizabeth Joy Paul, University of Maryland, abstract, slides
    [#s659, 15 Mar 2018]
    We present a method for stellarator coil design via gradient-based optimization of the coil-winding surface. The REGCOIL [Landreman, Nucl. Fusion 57, 046003 (2017)] approach is used to obtain the coil shapes on the winding surface using a continuous current potential. We apply the adjoint method to calculate derivatives of the objective function, allowing for efficient computation of analytic gradients while eliminating the numerical noise of approximate derivatives. We are able to improve engineering properties of the coils by targeting the root-mean-squared current density in the objective function. We obtain winding surfaces for W7-X and HSX which simultaneously decrease the normal magnetic field on the plasma surface and increase the surface-averaged distance between the coils and the plasma in comparison with the actual winding surfaces. The coils computed on the optimized surfaces feature a smaller toroidal extent and curvature and increased inter-coil spacing. A technique for computing the local sensitivity of figures of merit to normal displacement of the winding surface is presented, with potential applications for understanding engineering tolerances.
  • Reduced MHD and gyrokinetic studies on auroral plasmas
    Tomo-Hiko Watanabe, Nagoya University, abstract, slides
    [#s658, 14 Mar 2018]
    The reduced MHD and gyrokinetics developed from fusion theory have been applied to a variety of topics in space and astrophysical plasmas. Our theoretical and numerical methods developed from fusion studies at NU have also facilitated to understand key issues in auroral physics. Here, I would like to discuss some of the topics, such as structural formation of auroras and their dynamics, where competitive process of the ballooning, Kelvin-Helmholtz, and feedback instabilities in the magnetosphere-ionosphere (M-I) coupling are investigated [1, 2]. Furthermore, in the nonlinear stage of auroral growth, the M-I coupling system transits into the Alfvénic turbulence state [2]. The auroral turbulence is resulted from interactions of the shear Alfvén waves propagating in the opposite directions, and is an interesting application of the Goldreich-Sridhar theory. I would also like to discuss the gyrokinetic extension of the auroral theory including auroral electron acceleration [3].
    [1] T.-H. Watanabe, Phys. Plasmas 17, 022904 (2010)
    [2] Tomo-Hiko Watanabe, Hiroaki Kurata & Shinya Maeyama, New J. Phys. 18, 125010 (2016)
    [3] T.-H. Watanabe, Geophys. Res. Lett. 41, 6071 (2014)
  • Parametric Instability, Inverse Cascade, and the $1/f$ Range of Solar-Wind Turbulence
    Ben Chandran, University of New Hampshire, abstract, slides
    [#s626, 22 Feb 2018]
    Turbulence likely plays an important role in generating the solar wind, and spacecraft measurements indicate that solar-wind turbulence is largely non-compressive and Alfvén-wave-like. Although compressive fluctuations are sub-dominant, Alfvén waves in the solar wind couple to compressive slow magnetosonic waves (“slow waves”) via the parametric-decay instability. In this instability, an outward-propagating Alfvén wave decays into an outward-propagating slow wave and an inward-propagating Alfvén wave. In this talk, I will describe a weak-turbulence calculation of the nonlinear evolution of the parametric instability in the solar wind at wavelengths much greater than the ion inertial length under the assumption that slow waves, once generated, are rapidly damped. I'll show that the parametric instability leads to an inverse cascade of Alfvén-wave quanta and present several exact solutions to the wave kinetic equations. I will also present a numerical solution to the wave kinetic equations for the solar-wind-relevant case in which most of the Alfvén waves initially propagate away from the Sun in the plasma rest frame. In this case, the outward-propagating Alfvén waves evolve toward a $1/f$ frequency spectrum that shows promising agreement with spacecraft measurements of interplanetary turbulence in the fast solar wind. I will also present predictions that will be tested by NASA's upcoming Solar Probe Plus mission, which will travel much closer to the Sun than any previous spacecraft.
  • On the Global Attractor of 2D Incompressible Turbulence with Random Forcing
    John Bowman, U. Alberta, abstract, slides
    [#s630, 19 Feb 2018]
    We revisit bounds on the projection of the global attractor in the energy-enstrophy plane obtained by Dascaliuc, Foias, and Jolly [1,2]. In addition to providing more elegant proofs of some of the required nonlinear identities, the treatment is extended from the case of constant forcing to the more realistic case of random forcing. Numerical simulations in particular often use a stochastic white-noise forcing to achieve a prescribed mean energy injection rate. The analytical bounds are illustrated numerically for the case of white-noise forcing.
    [1] R. Dascaliuc, C. Foias & M.S. Jolly, J. Dynam. Differential Equations 17, 643 (2005)
    [2] R. Dascaliuc, C. Foias & M.S. Jolly, J. Differential Equations 248, 792 (2010)
  • Magnetic Reconnection in Three Dimensional Space
    Prof. Allen Boozer, Columbia University, abstract, slides
    [#s623, 02 Feb 2018]
    The breaking of magnetic field line connections is of fundamental importance in essentially all applications of plasma physics: laboratory to astrophysics. For sixty years the theory of magnetic reconnection has been focused on two-coordinate models. When dissipative time scales far exceed natural evolution times, such models are not realistic for ordinary three dimensional space. The ideal (dissipationless) evolution of a magnetic field is shown to in general lead to a state in which the magnetic field lines change their connections on an Alfvénic (inertial), not resistive, time scale. Only a finite mass of the lightest current carrier, the electron, is required. During the reconnection, the gradient in $j_\parallel/B$ relaxes while conserving magnetic helicity in the reconnecting region. This implies a definite amount of energy is released from the magnetic field and transferred to shear Alfvén waves, which in turn transfer their energy to the plasma. When there is a strong non-reconnecting component of the magnetic field, called a guide field, $j_\parallel/B$ obeys the same evolution equation as that of an impurity being mixed into a fluid by stirring. Although the enhancement of mixing by stirring has been recognized by every cook for many millennia, the analogous effect in magnetic reconnection is not generally recognized. An interesting mathematical difference is a three-coordinate model is required for the enhancement of magnetic reconnection while only two coordinates are required in fluid mixing. The issue is the number of spatial coordinates required to obtain an exponential spatial separation of magnetic field lines versus streamlines of a fluid flow.
  • A moment approach to plasma fluid/kinetic theory and closures
    Jeong-Young Ji, Utah State University, abstract, slides
    [#s611, 11 Jan 2018]
    A system of exact fluid equations always involves more unknowns than equations. This is called the closure problem. An important aspect of obtaining quantitative closures is an accurate account of collisional effects. Recently, analytical calculations of the Landau (Fokker-Planck) collision operator as well the derivation of an infinite hierarchy of moment equations have been carried out using expansions for the distribution function in terms of irreducible Hermite polynomials. In this talk, I will present solutions to the moment hierarchy that provide closure for the set of five-moment fluid equations. In the collisional limit, improved Braginskii closures are obtained by increasing the number of moments and considering the ion-electron collision effects. For magnetized plasmas, I highlight the effect of long mean free path and derive parallel integral closures for arbitrary collisionality. Finally, I will show how the integral closures can be used to study radial transport due to magnetic field fluctuations and electron parallel transport for arbitrary collisionality.
  • Overview of ENN and Fusion Technology
    Dr. Minsheng Liu, ENN Sci. & Tech. Co. Ltd., China, abstract
    [#s602, 19 Dec 2017]
    The talk includes 3 parts: 1.overview of ENN; 2.ENN research areas and achievements; and 3.ENN fusion technology roadmap.
  • The Development and Applications of CLT
    Wei Zhang, Zhejiang University, abstract
    [#s599, 19 Dec 2017]
    CLT is an explicit, three-dimensional, fully toroidal, non-reduced, Hall-MHD code, which is developed in Zhejiang University. Through CLT, I find that electron diamagnetic rotation, which is well described by the CLT code, can significantly modify the dynamics of the tearing mode. It can also affect the characteristics of the sawtooth oscillations. Besides, I have also studied the influence of driven current on tearing mode instabilities, which can explain some experimental data of EAST. Now CLT is updating to study the influence of RMPs (Resonant Magnetic Perturbations) on tearing mode instabilities. Preliminary results show that the threshold for ‘mode locking’ increases with the frequency of RMPs, which is consistent with theory prediction.
  • The JET 2020 Program: D-T and Alpha physics with the ITER Like Wall
    M. Romanelli, JET, Culham Science Centre, abstract, slides
    [#s560, 18 Dec 2017]
    The 2018-2019 JET program in preparation for DTE2 has been launched in November 2017 and will start in March 2018 with the first of the 2018 analysis and modeling campaigns followed by the first experimental campaign in D lasting until October 2018. The program will feature also an H campaign followed by a full T campaign in 2019 and aims at studying differences in plasma dynamics related to operation with pure and mixed hydrogen isotopes along with the impact of fast particles on heating, transport and confinement followed by a H campaign with NBI in H, and a D campaign for finalising the plasma preparation for DTE2. Extensive modeling and analysis has been devoted to the preparation of successful scenarios that will allow to achieving during D-T operations a fusion power of 15MW continuously for at least 5s and observing clear alpha-particle effects. In this seminar I’ll present the striking results on the isotope (H/D) effect on transport and global confinement from the latest experimental campaigns along with the plans for experiments on transport and confinement in the 2018-2020 campaigns.
    e-mail contact of the main author: michele.romanelli@ukaea.uk
  • Understanding mechanisms underlying ohmic breakdown in a tokamak by considering multi-dimensional plasma responses
    Min-Gu Yoo, Seoul National University, abstract, slides
    [#s598, 18 Dec 2017]
    The ohmic breakdown is generally used to produce initial plasmas in tokamaks. However, the complex electromagnetic structure of tokamaks has obscured the physical mechanism of ohmic breakdown for several decades. Previous studies ignored plasma responses to external electromagnetic fields and adopted only the simplest Townsend avalanche theory. However, we found clear evidence that experimental results cannot be explained by the Townsend theory. Here, we propose a completely new type of breakdown mechanism that systematically considers multi-dimensional plasma responses in the complex electromagnetic topology. As the plasma response, self-electric fields produced by space-charge were found to be crucial for significantly reducing plasma density growth rate and enhancing perpendicular transport via $\small {\bf E}\times{\bf B}$ drifts. A particle simulation code, BREAK, clearly captured these effects and provided a remarkable reproduction of the mysterious experimental results in KSTAR. These new physical insights into complex electromagnetic topology provide general design guideline for a robust breakdown scenario in a tokamak fusion reactor.
  • Investigation of whistler-electron interaction in Earth’s radiation belt
    Lei Zhao, New Mexico Consortium, abstract, slides
    [#s570, 15 Dec 2017]
    In this talk, I will focus on how to describe the energetic electron dynamics through quasilinear whislter-electron interactions in Earth’s magnetosphere. First of all, we explore gyro-resonant wave-particle interaction and quasi-linear diffusion in different magnetic field configurations related to the March 17 2013 storm. We consider the Earth’s magnetic dipole field as a reference, and compare the results against non-dipole field configurations corresponding to quiet and stormy conditions. The latter are obtained with RAM-SCB, a code that models the Earth’s ring current and provides a realistic model of the Earth’s magnetic field. By applying quasi-linear theory, the bounce- and magnetic local time (MLT) averaged electron pitch angle, mixed term and energy diffusion coefficients are calculated for each magnetic field configuration. The results show that the diffusion coefficients become quite independent of the magnetic field configuration for relativistic electrons (∼ 1MeV ), while a realistic model of the magnetic field configuration is necessary to adequately describe the diffusion rates of lower energy electrons (∼ 100keV )
    In the second part of the talk, a Test Particle Model (TPM) is used to explore the limitations of quasilinear theory when applied to whistler wave-electron resonance. We consider the influence of wave amplitude and wave bandwidth on the wave-particle interaction within the quasilinear theory limit. The result implies that quasilinear theory tends to breakdown more easily for energetic particles even at small wave amplitude. While for broad wave bandwidths, it allows for more stochasticity for particles in phase space. On the contrary, electron phase trapping and bunching (nonlinear wave-particle interaction phenomenon) will dominate a narrow bandwidth spectrum that resembles a monochromatic wave.
  • Nonlinear ECDI and anomalous transport in $\small {\bf E} \times {\bf B}$ discharges
    S. Janhunen, U. Saskatchewan, abstract, slides
    [#s559, 14 Dec 2017]
    Cross-field anomalous transport is an important feature affecting the operation and performance of $\small {\bf E} \times {\bf B}$ discharges. Instabilities excited by the $\small {\bf E} \times {\bf B}$ flow cause anomalous current to develop, elicited in the nonlinear regime as a large amplitude coherent wave driven by the energy input from the unstable cyclotron resonances. A persistent train of soliton-like waves characterized by the fundamental harmonic of the electron cyclotron oscillations appears in the ion density. Simultaneously, there is an observable energy cascade toward long wavelength (inverse cascade) which is manifested by the formation of the long wavelength envelope of the wave train in 1D simulations. It is shown that the long wavelength part of the turbulent spectrum provides a dominant contribution to anomalous electron mobility.[1]
    We present results from high fidelity 1D3V and 2D3V particle-in-cell simulations for a simplified Hall-effect thruster like plasma. We describe the non-linear evolution of the system and speculate on mechanisms behind the coherent structures and their interactions. The 2D the picture is complicated by the existence of a simultaneous long-wavelength mode (modified two-stream instability), in addition to the non-linear cascades observable in 1D.
    [1] S. Janhunen, A. Smolyakov et al., Phys. Plasmas 25, 011608 (2018)
    [2] Movie 1) 1D case with $T_e = 10 eV$
    [2] Movie 2) 1D with $T_e = 0 eV$
    [2] Movie 3) Ion density fluctuations in 2D simulation
    [2] Movie 4) Electron temperature in 2D simulation
  • A fast integral equation based solver for the computation of Taylor states in toroidal geometries
    Antoine Cerfon, Courant Institute, New York University , abstract, slides
    [#s531, 08 Nov 2017]
    The stellarator equilibrium code SPEC [S.R. Hudson et al.,, Phys. Plasmas 19, 112502 (2012)] computes 3D equilibria by subdividing the plasma into separate regions assumed to have undergone Taylor relaxation to a minimum energy state subject to conserved fluxes and magnetic helicity, and separated by ideal MHD barriers. In this talk, we present a numerical scheme for the fast and high order computation of Taylor states in toroidal regions based on an integral formulation of the problem. Our formulation offers the advantage that the unknowns are only defined on the boundary of the toroidal regions. As a result, high accuracy is reached with a small number of unknowns, leading to a code which is fast and has low memory requirements. In the context of SPEC, in which the locations of the ideal MHD interfaces are iteratively updated until force balance is satisfied at each interface, our formulation gives the possibility to apply the entire iterative procedure without ever discretizing the plasma volume, only discretizing the ideal interfaces.

    This is joint work with M. O'Neil, L. Greengard, and L.-M. Imbert-Gerard

  • The shearing modes approach to the theory of plasma shear flow
    Vladimir Mikhailenko, Pusan National University, South Korea , abstract, slides
    [#s502, 18 Oct 2017]
    The basic point in the understanding the processes of the evolution of the instabilities and turbulence in the plasma shear flows across the magnetic field is the proper treatment the effects of the persistent distortion of the perturbations by the shearing flow, particularly in the applications of the spectral transforms to the governing equations. The problem of extracting the separate spatial Fourier mode in the stability theory of the plasma shearing flows, and joined with it the problem of the limits of the applicability to such a plasma the spectral methods of the investigation of the stability properties on the base of the dispersion equations, may be resolved by employing the non-modal fluid and kinetic theory, which grounds on the methodology of the shearing modes and completely involves the effect of the persistent deformation of the perturbations by the sheared flows. That theory displays, that the application of the methodology of the shearing modes and convective-shearing coordinates and solution of the initial value problem in the wave vector and time variables, instead of the application of the static spatial Fourier modes and spectral transform in time, has the decisive impact on the understanding the wave-particle interaction in the shearing flow. That theory recovers main linear and non-linear processes and corresponding numerous characteristic time scales, which may be observable in the experiments and numerical simulations, but can’t be distinguished, when the spectral transform in time is applied. A most famous is the ”quench rule”, which was detected in first in the numerical simulations, but was not confirmed analytically in the calculations of the shearing flows stability, grounded on the spectral transformations in time. The primary intent of this report is to show that a non-modal approach is a decisive in the reconciling observational evidence with stability theory of plasma shearing flows and to suggest a more frequent use of that approach.
  • Gyrokinetic simulation of boundary plasma in contact with material wall
    Seung-Hoe Ku, PPPL, abstract, slides
    [#s449, 11 Sep 2017]
    Boundary plasma is in a non-equilibrium statistical state governed by self-organization among multiscale physics, and needs to be modeled with total-f gyrokinetic equations. A unique pariticle-in-cell technique will be introduced that has enabled XGC to be the world’s first, and only, gyrokinetic code that simulates the boundary plasma across the magnetic separatrix into the scrape-off layer in contact with a material wall. Examples of successful boundary physics discoveries enabled by the technique will also be presented.
  • Role of electron physics in 3D two-fluid 10-moment simulations of the Ganymede’s magnetosphere
    Liang Wang, U. New Hampshire, abstract, slides
    [#s445, 31 Aug 2017]
    We studied the role of electron physics in 3D two-fluid 10-moment simulations of the Ganymede’s magnetosphere. The model captures non-ideal physics like the Hall effect, the electron inertia, and anisotropic, non-gyrotropic pressure effects. A series of analyses were carried out: 1) The resulting magnetic field topology and electron and ion convection patterns were investigated. The magnetic fields were shown to agree reasonably well with in-situ measurements by the Galileo satellite. 2) The physics of collisionless magnetic reconnection were carefully examined in terms of the current sheet formation and decomposition of generalized Ohm’s law. The importance of pressure anisotropy and non-gyrotropy in supporting the reconnection electric field is confirmed. 3) We compared surface “brightness” morphology, represented by surface electron and ion pressure contours, with oxygen emission observed by the Hubble Space Telescope (HST). The correlation between the observed emission morphology and spatial variability in electron/ion pressure was demonstrated. We also briefly discussed the relevance of this work to the future JUICE mission.
  • A new coil design code FOCUS for designing stellarator coils without the winding surface
    Caoxiang Zhu, U. Sci. & Tech. China, Hefei, China , abstract, slides
    [#s444, 30 Aug 2017]
    Finding an easy-to-build coils set has been a critical issue for stellarator design for decades. Conventional approaches assume a toroidal “winding” surface and suffer the difficulties of nonlinear optimization. A new coil design method, the FOCUS code, is introduced by representing each discrete coil as an arbitrary, closed space curve. The first and second derivatives of the target function that covers both physical requirements and engineering constraints are calculated analytically. We have employed several advanced nonlinear optimization algorithms, like the nonlinear conjugate gradient and the modified Newton method, for minimizing the target function. Numerical illustrations show that the new method can be applied to different types of coils for various configurations with great flexibilities and robustness. An extended application for analyzing the error field sensitivity is also presented.
  • Stochastic modelling of fluctuations in scrape-off layer plasmas
    Ralph Kube, UIT, Tromso, Norway, abstract, slides
    [#s412, 24 Aug 2017]
    Scrape-off layer plasmas feature intermittent, large-amplitude fluctuations which are attributed to the radial outwards propagation of plasma blobs through this volume. We introduce a stochastic model which describes scrape-off layer time series by superposing uncorrelated pulses with variable amplitude. The resulting time series is Gamma distributed where the lowest order statistical moments are given by the pulse parameters. The power spectral density is governed by the pulse shape - for a double exponential pulse shape presents the power spectral density a flat region for low frequencies and a steep power law scaling for high frequencies. Predictions from this model are compared to fluctuation measurements of electron density, temperature and plasma potential by mirror Langmuir probes in the SOL during an ohmic L-mode discharge in Alcator C-Mod.
  • Simulation of resonant wave-particle interaction in tokamaks
    Meng Li, IFS, The University of Texas, Austin, abstract, slides
    [#s415, 17 Aug 2017]
    We present a numerical procedure for modeling resonant response of energetic particles to waves in tokamaks. With Littlejohn Lagrangian for guiding center motion, we use the action-angle variables to simplify simulations of the fast ion dynamics to one dimensional. The transformation also involves construction of canonical straight field line coordinates, which render a Hamiltonian description of the guiding center motion. This module can be integrated with the modified MHD code AEGIS to simulate wave-particle interactions.
  • Tearing mode dynamics and sawtooth oscillation in Hall-MHD
    Zhiwei Ma, Zhejiang University, Hangzhou, China, abstract, slides
    [#s414, 09 Aug 2017]
    Tearing mode instability is one of the most important dynamic processes in space and laboratory plasmas. Hall effects, resulted from the decoupling of electron and ion motions, could cause the fast development and perturbation structure rotation of the tearing mode and become non-negligible. A high accuracy nonlinear MHD code (CLT) is developed to study Hall effects on the dynamic evolution of tearing modes with Tokamak geometries. It is found that the rotation speed of the mode structure from the simulation is in a good agreement with that from the analytical theory in a single tearing mode. The linear growth rate increases with increase of the ion skin depth. The self-consistent generated rotation largely alters the dynamic behaviors of the double tearing mode and the sawtooth oscillation.
  • Continuum kinetic schemes for Vlasov-Maxwell equations, with applications to laboratory and space plasmas
    Ammar Hakim, PPPL, abstract, slides
    [#s179, 22 Jun 2017]
    Increasingly accurate laboratory experiments and satellite observations have led to a "golden age" in plasma physics. Detailed kinetic features, including distribution functions, can now be measured in-situ, putting severe strain on theory and modeling to explain the experiments/observations. The Particle-In-Cell (PIC) method remains a powerful and widely used tool to study such kinetic physics numerically. Recently, complimenting the PIC approach, significant progress has been made in discretizing the Vlasov equation directly, treating it as a partial differential equation (PDE) in 6D phase-space. In this talk, I present a high-order discontinuous Galerkin (DG) algorithm to solve the Vlasov equation. This continuum scheme leads to noise-free solutions and, with the use of specially optimized computational kernels, can be very efficient, in particular, for problems in which the structure of distribution functions and its higher-moments are required. In addition, with a proper choice of basis functions and numerical fluxes, the scheme conserves energy exactly, while conserving momentum to a high degree of accuracy. Applications of the scheme to (a) kinetic saturation mechanism of Wiebel-like instability and (b) turbulence in the solar wind are presented. We demonstrate a new and novel mechanism for the nonlinear saturation of the Wiebel instability that comes about from a balance between filamentation and a secondary electrostatic two-stream instability. We use 5D simulations of turbulent plasmas to study detailed kinetic physics of magnetized turbulence, showing that the solution contains remarkable amount of detail in the distribution function, leading to new and novel insights the nature of kinetic wave-particle exchange in turbulent plasmas.
  • Nonthermal particle acceleration in magnetic reconnection and turbulence in collisionless relativistic plasmas
    Dmitri Uzdensky, IAS & U. Colorado Boulder, abstract, slides
    [#s177, 12 Jun 2017]
    One of the key recurrent themes in high-energy plasma astrophysics is relativistic nonthermal particle acceleration (NTPA) necessary to explain the bright X-ray and gamma-ray flaring emission with ubiquitous power-law spectra in astrophysical objects such as pulsar wind nebulae, hot accretion flows and coronae of accreting black holes, and black-hole powered relativistic jets in active galactic nuclei and gamma-ray bursts. Two leading physical processes often invoked as possible NTPA mechanisms are collisionless magnetic reconnection and turbulence. In order to understand these processes, as well as their resulting observable radiation signatures, I have recently initiated a broad theoretical and computational research program in kinetic radiative plasma astrophysics. This program employs large-scale first-principles particle-in-cell kinetic simulations (including those that self-consistently incorporate radiation-reaction effects) coupled with analytical theory. In this talk I will review the resulting progress that we have achieved in recent years towards understanding and quantitative characterization of NTPA in reconnection and turbulence over a broad range of physical regimes. I will present 2D and 3D simulation results that demonstrate that both reconnection and turbulence in relativistic collisionless astrophysical plasmas can robustly produce non-thermal energy spectra with power-law indices that show an intriguingly similar characteristic dependence on the plasma magnetization. I will also describe the effects of strong radiative cooling on reconnection and turbulence.
  • Hot Particle Equilibrium code (HPE) with plasma anisotropy and toroidal rotation
    Leonid Zakharov, LiWFusion, abstract, slides
    [#s162, 26 May 2017]
    The HPE code represents the extension of ESC (Equilibrium and Stability Code) to tokamak equilibria with plasma anisotropy and toroidal rotation. In addition to conventional 1-D input profiles of plasma pressure $p(a)$ and safety factor $q(a)$ HPE accepts the poloidal Mach number ${\cal M}(a)$ of plasma rotation and 2-D parallel pressure profile of hot particles $p_\parallel(a,B)$ as the input. Here, $a$ is the normalized radial flux coordinate of magnetic configuration and $B$ is the strength of magnetic field. The HPE code includes the effect of finite width of hot particle orbits and for the case of powerful NBI injection the code can generate the plasma equilibria for theory needs as well as use the experimental (or kinetic simulations) data for interpretation of experiments.
  • Anomalous Diffusion across a Tera Gauss Field in accreting Neutron Stars
    Russell Kulsrud, PPPL, abstract, slides
    [#s161, 25 May 2017]
    A large amount of mass falls on the polar region of neutron star in Xray binaries and the question is, is the mass completely frozen on the field lines or can it diffuse through them? In this talk we present a mechanism for the latter possibility. A strong MHD instability occurs in the top layers of the neutron star driven by the incoming mass. This instability has the same properties as the Schwarzschild instability in the solar convection zone. It gives rise to a turbulent cascade which mixes up the field lines so that lines originally far apart can come with a resistivity diffusion distance and transfer the masses between them. However, the lines of force themselves are not disrupted. This leads to an equilibrium which is marginal with respect to the instability just as happens in the Schwarzschild case.
  • Modeling substrom dipolarizations and particle injections in the terrestrial magnetosphere
    Konstantin Kabin, Royal Military College of Canada, abstract, slides
    [#s156, 11 May 2017]
    Increased fluxes of energetic electrons and ions in the inner magnetosphere of the Earth are often associated with sudden reconfigurations of the magnetotail, often referred to as substorm dipolarizations. We describe a novel model of the magnetotail which is easily controlled by several adjustable parameters, such as the thickness of the tail and the location of transition from dipole-like to tail-like magnetic field lines. This model is fully three-dimensional and includes the day-night asymmetry of the terrestrial magnetosphere, however, the field lines are confined to the meridional planes. Our model is well suited to studies of the magnetotail dipolarizations which we consider to be the tailward movements of the transition between dipole-like and tail-like field lines. We also study the effects of a dipolarizing electromagnetic pulse propagating towards the Earth. The calculated electric and magnetic fields are used to describe the motion of electrons and ions and changes in their energies. In some cases, particle energies increase by a factor of 25 or more. The energized particles are transported earthward where they are often observed by geostationary satellites as substorm injections. The energization level obtained in our model is reasonably consistent with satellite and ground-based observation (e.g. carried out by riometers), and therefore we consider our scenario of the dipolarization process to be feasible.
  • DCON for Stellarators
    Alan Glasser, Fusion Theory & Computation, Inc., abstract, slides
    [#s158, 10 May 2017]
    We report the development of a new version of DCON for nonaxisymmetric toroidal plasmas, e.g. stellarators. The DCON code is widely used for fast and accurate determination of ideal MHD stability of axisymmetric toroidal plasmas. Minimization of the ideal MHD energy principle $\delta W$ is reduced to adaptive numerical integration of a coupled set of ordinary differential equations for the complex poloidal Fourier harmonics of the perturbed normal displacement. For a periodic cylindrical plasma, both the poloidal and toroidal coordinates are ignorable, allowing for treatment of single harmonics $m$ and $n$. For an axisymmetric toroidal plasma, poloidal symmetry is broken, causing different $m$’s to couple and requiring simultaneous treatment of $M$ harmonics. For a nonaxisymmetric plasma, toroidal symmetry is also broken, causing different $n$’s to couple and requiring treatment of $N$ harmonics. For a stellarator with field period $l$, e.g. $l = 5$ for W7X, each toroidal harmonic $n$ is coupled to toroidal harmonics $n+k \, l$, for all integers $k$. Both $M$ and $N$ are truncated on the basis of convergence. Singular surfaces occur at all values of safety factor $q = m/n$ in the plasma. The DCON equations have been generalized to allow for multiple $n$’s with coupling matrices $F$, $K$, and $G$. An interface has been developed for the nonaxisymmetric equilibrium code VMEC, which provides values for these matrices. Their Fourier harmonics are fit to cubic splines as a function of the radial coordinate $s$, allowing for adaptive integration of the ODEs. The status of the code development will be presented.
  • Disintegration threshold of Langmuir solitons in inhomogeneous plasmas
    Yasutaro Nishimura, National Cheng Kung University, Taiwan, abstract, slides
    [#s140, 04 May 2017]
    Dynamics of Langmuir solitons in inhomogeneous plasmas is investigated numerically employing Zakharov equations. The solitons are accelerated toward a lower background density side. With the steep density gradients, balance between the electric field part of the soliton and the density cavity breaks and the solitons disintegrate. The disintegration threshold is given by regarding the electric field part of the soliton as a point mass moving along the self-generated potential well produced by the density cavity. On the other hand, when the density gradient is below the threshold, Langmuir solitons adjust themselves by expelling the imbalanced portion as density cavities at the sound velocity. When the gradient is below the threshold, the electric field part of the soliton bounces back and forth within the potential well. The study is extended to kinetic simulation. Generation mechanism of high energy electron tails in the presence of solitons is discussed. The electron distribution function resembles that of the Lorentzian type. The particle acceleration is explained as a transport process toward high energy side due to overlapping of multiple resonant islands in the phase space.
  • Gyrokinetic simulation of a fast L-H bifurcation dynamics in a realistic diverted tokamak edge geometry
    Seung-Hoe Ku, PPPL, abstract, slides
    [#s141, 27 Apr 2017]
    Despite its critical importance in the fusion program and over 30 years of H-mode operation, there has been no fundamental understanding at the kinetic level on how the H-mode bifurcation occurs. We report the first observation of an edge transport barrier formation event in an electrostatic gyrokinetic simulation carried out in a realistic C-Mod like diverted tokamak edge geometry under strong forcing by a high rate of heat deposition. The results show that the synergistic action between two multiscale dynamics, the turbulent Reynolds-stress driven [1] and the neoclassical X-point orbit loss drive [2] sheared ${\bf E}\times{\bf B}$ flows, works together to quench turbulent transport and form a transport barrier just inside the last closed magnetic flux surface. The synergism helps reconcile experimental reports of the key role of turbulent stress in the bifurcation [3, and references therein] with some other experimental observations that ascribe the bifurcation to X-point orbit loss/neoclassical effects [4,5]. The synergism could also explain other experimental observations that identified a strong correlation between the L-H transition and the orbit loss driven ${\bf E}\times{\bf B}$ shearing rate [6,7]. The synergism is consistent with the general experimental observation that the L-H bifurcation is more difficult with the $\nabla B$-drift away from the single-null X-point, in which the X-point orbit-loss effect is weaker [2].
    [1] P.H. Diamond, S-I Itoh et al., Plasma Phys. Controlled Fusion 47, R35 (2005)
    [2] C.S. Chang, Seunghoe Kue & H. Weitzner, Phys. Plasmas 9, 3884 (2002)
    [3] G.R. Tynan, M. Xu et al., Nucl. Fusion 53, 073053 (2013)
    [4] T. Kobayashi, K. Itoh et al., Phys. Rev. Lett. 111, 035002 (2013)
    [5] M. Cavedon, T. Pütterich et al., Nucl. Fusion 57, 014002 (2017)
    [6] D.J. Battaglia, C.S. Chang et al., Nucl. Fusion 53, 113032 (2017)
    [7] S.M. Kaye, R. Maingi et al., Nucl. Fusion 51, 113109 (2011)
  • Parasitic momentum flux in the tokamak core
    Timothy Stoltzfus-Dueck, PPPL, abstract
    [#s138, 20 Apr 2017]
    Tokamak plasmas rotate spontaneously in the absence of applied torque. This so-called “intrinsic rotation” may be very important for future low-torque devices such as ITER, since rotation can stabilize certain instabilities. In the mid-radius “gradient region”, which reaches from the sawtooth inversion radius out to the pedestal top, intrinsic rotation profiles are sometimes flat and sometimes hollow. Profiles may even transition suddenly between these two states, an unexplained phenomenon referred to as rotation reversal. Theoretical efforts to identify the origin of the mid-radius rotation shear have focused primarily on quasilinear models, in which the phase relationships of some selected instability result in a nondiffusive momentum flux (“residual stress”). In contrast to these efforts, the present work demonstrates the existence of a robust, fully nonlinear symmetry-breaking momentum flux that follows from the free-energy flow in phase space and does not depend on any assumed linear eigenmode structure. The physical origin is an often-neglected portion of the radial ${\bf E}\times {\bf B}$ drift, which is shown to drive a symmetry-breaking outward flux of co-current momentum whenever free energy is transferred from the electrostatic potential to ion parallel flows [1]. The resulting rotation peaking is counter-current and scales as temperature over plasma current. As originally demonstrated by Landau [2], this free-energy transfer (thus also the corresponding residual stress) becomes inactive when frequencies are much higher than the ion transit frequency, which may explain the observed relation of density and counter-current rotation peaking in the core. Simple estimates suggest that this mechanism may be consistent with experimental observations, in both hollow and flat rotation regimes.
    [1] T. Stoltzfus-Dueck, Phys. Plasmas 24, 030702 (2017)
    [2] L. Landau, J. Exp. Theor. Phys. 16, 574 (1946); English translation in J. Phys. (USSR) 10, 25 (1946)
  • Transport in the Coupled Pedestal and Scrape-off layer region of H-mode plasmas
    Michael Churchill, PPPL, abstract, slides
    [#s139, 18 Apr 2017]
    Understanding the multi-scale neoclassical and turbulence physics in the edge region (pedestal + scrape-off layer) is required in order to reliably predict performance in future fusion devices. I will present research exploring characteristics of this transport using the family of X-point Gyrokinetic Codes (XGC). First, the variation of pressure in the scrape-off layer (important to understand in order to avoid divertor wall degradation) is widely believed to follow simple fluid prescriptions, due to high collisionality. However, simulation results in the near-SOL indicate a significant departure from the simple fluid models, even after including additional terms from neutral drag and the Chew-Goldberger-Low form of parallel ion viscosity to the parallel momentum balance. Second, turbulence characteristics in the edge region show nonlocal behavior, including convective transport of turbulent eddies (“blobs”) born just inside the closed field line region out into the SOL. These large intermittent structures can be created even in the absence of collisions in the simulation. Tracking these structures show that on average their radial velocity is very small, even in the near-SOL, while their poloidal velocity is significant. The potential structure within these blobs is monopolar with a peak shifted from the density structure, contrary to the dipolar structure expected by analytical models for blob generation and transport based on interchange turbulence. Finally, coherent phase space structures in blobs are searched for, however only broad regions of velocity space are found to show significant structure.
  • Improved Neoclassical Transport Simulation for Helical Plasmas
    Botsz Huang, The Graduate University for Advanced Studies, Japan, abstract
    [#s135, 30 Mar 2017]
    This work contains the following studies: (1) benchmark of the neoclassical transport models and clarify the impact of their approximations among drift-kinetic equations in helical plasmas; and (2) apply the local models for a quantitative evaluation of bootstrap current in FFHR-d1 DEMO. In part 1, this work reports that the benchmarks of the neoclassical transport codes based on the several local drift-kinetic models. The drift-kinetic models are zero-orbit-width (ZOW), zero-magnetic-drift (ZMD), DKES-like, and global, as classified in [1]. The magnetic geometries of HSX, LHD, and W7-X are employed. It is found that the assumption of ${\bf E} \times {\bf B}$ incompressibility causes discrepancy of neoclassical radial flux and parallel flow among the models when ${\bf E} \times {\bf B}$ is sufficiently large compared to the magnetic drift velocities. For example, ${\cal M}_p ≤ 0.4$ where ${\cal M}_p$ is the poloidal Mach number. On the other hand, when ${\bf E} \times {\bf B}$ and the magnetic drift velocities are comparable, the tangential magnetic drift plays a role in suppressing the unphysical peaking of neoclassical radial-fluxes in the other local models at $E_r ≃ 0$. In low collisional frequency plasmas, the tangential drift suppress such unphysical behavior in the radial transport well. This work demonstrates that the ZOW model not only mitigates the unphysical behavior but also implements evaluation of bootstrap current in LHD with the low computation cost compared to the global model. In part 2, the impact of the parallel momentum conservation on the boot-strap current among the ZOW, DKES, and PENTA models are demonstrated. The ZOW model is extended to include the ion parallel mean flow effect on the electron-ion parallel friction. The DKES model employs the pitch-angle-scattering as the collision operator. The PENTA code employs the Sugama-Nishimura method to correct the momentum balance of DKES results. Therefore, the ZOW model and PENTA codes both conserve the parallel momentum in like-species collisions and include the electron-ion parallel frictions. The work shows that the ZOW and the PENTA code agree with each other well on the calculations of the bootstrap current. Then, this verifies the reliability of the bootstrap current calculation with the ZOW model and the PENTA code for FFHR-d1 DEMO.
    [1] Seikichi Matsuoka, Shinsuke Satake et al., Phys. Plasmas 22, 072511 (2015)
  • Recent progress of understanding 3D magnetic topology in stellarators and tokamaks
    Yasuhiro Suzuki, National Institute for Fusion Science, Japan, abstract, slides
    [#s134, 23 Mar 2017]
    Recent progress of 3D equilibrium calculation will be reported. The 3D equilibrium calculation is a fundamental consideration to understand the magnetic topology. In particular, for stellarators, the topological change due to the beta-sequences has been found and it is an important issue how to make good flux surfaces in general 3D configuration in the presence of the plasma beta. On the other hand, the 3D equilibrium calculation is also an important issue for tokamaks, because the RMP is widely used to control the stability and transport. In this talk, recent results of 3D equilibrium calculations based on the HINT code, which is a 3D equilibrium calculation code without an assumption of perfectly nested flux surfaces. Impacts of the beta-sequences and the plasma rotation will be discussed.
  • Self-organizing knots
    Christopher Smiet, Leiden University, The Netherlands, abstract, slides
    [#s133, 16 Mar 2017]
    Magnetic helicity, a measure for the linking and knotting of magnetic field lines, is a conserved quantity in Ideal MHD. In the presence of resistivity, helicity constrains the rate at which magnetic energy can be dissipated. When a localized, helical magnetic field is set to relax in a low-resistance high-beta plasma, the magnetic pressure drives the plasma to expand whilst the helicity is still approximately conserved. Using numerical simulations I show how this interplay gives rise to a novel MHD equilibrium: the initially linked field lines self-organize to form a structure where field lines lie on nested toroidal surfaces of constant pressure. The Lorentz forces are balanced by the gradient in pressure, with a minimum in pressure on the magnetic axis. Interestingly, the rotational transform is nearly constant on all magnetic surfaces, making the structure topologically nearly identical to a famous knotted structure in Topology: the Hopf fibration. I will explore the nature of this equilibrium, and how it relates geometrically to the structure of the Hopf map. Additional dynamics give rise phenomena that are well known from magnetic confinement devices; magnetic islands can occur at rational surfaces, and in certain regimes the equilibrium becomes nonaxisymmetric, triggering a marginal core-interchange mechanism.
  • Development and application of BOUT++ for large scale turbulence simulation
    Jarrod Leddy - The University of York, UK, abstract, slides
    [#s132, 28 Feb 2017]
    The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the self-consistent evolution of plasma fluctuations, global profiles, and flows, but the numerical tools capable of doing this in realistic (diverted) geometry are only now being developed. BOUT++ is one such tool that has had many recent develops towards this goal. A novel coordinate system has been developed to improve the resolution around the X-point and strike points in the divertor region. A 5-field reduced 2-fluid plasma model for the study of instabilities and turbulence in magnetised plasmas has been built on the BOUT++ framework that allows the evolution of global profiles, electric fields and flows on transport timescales, with flux-driven cross-field transport determined self-consistently by electromagnetic turbulence. Models for neutral evolution have also been included, and the interaction of these neutrals with the plasma is characterised through charge exchange, recombination, ionisation, and radiation. Simulation results for linear devices, MAST-U, and DIII-D are presented that shed light on the nature of plasma-neutral interaction, detachment in the super-X divertor, and turbulence in diverted geometry.
  • Studies on proton effective heating in magnetic reconnection by means of particle simulations
    Shunsuke Usami, National Institute for Fusion Science, Japan, abstract
    [#s130, 23 Feb 2017]
    By means of two-dimensional electromagnetic particle simulations, ion heating mechanism is investigated in magnetic reconnection with a guide magnetic field. These simulations mimic dynamics of two torus plasmas merging through magnetic reconnection in a spherical tokamak (ST) device. It is found that protons are effectively heated in the downstream by the pickup mechanism, since a ring-like structure of proton velocity distribution, which is theoretically predicted to be formed by picked-up ions, is observed at a local region of the downstream. Furthermore, based on the theory by J. F. Drake et al., only heavy ions were believed to be heated by suffering the pickup mechanism, however it is pointed out that the pickup of protons is consistent with the theory in the cases that the plasma beta is much less than 1 in the upstream, which can be satisfied in STs.
  • Low-Frequency $\delta f$ PIC Models with Fully Kinetic Ions
    Benjamin J. Sturdevant, University of Colorado at Boulder, USA, abstract, slides
    [#s131, 22 Feb 2017]
    A fully kinetic ion model is useful for the verification of gyrokinetic turbulence simulations in certain regimes where the gyrokinetic model may break down due to the lack of small ordering parameters. For a fully kinetic ion model to be of value, however, it must first be able to accurately simulate low-frequency drift-type instabilities typically well within the domain of gyrokinetics. In this talk, we present a fully kinetic ion model formulated with weak gradient drive terms and applied to the ion-temperature-gradient (ITG) instability. A $\delta f$ implementation in toroidal geometry is discussed, where orthogonal coordinates are used for the particle dynamics, but field-line-following coordinates are used for the field equation, allowing for high resolution of the field-aligned mode structure. Variational methods are formulated for integrating the particle equations of motion, allowing for accuracy on a long time scale with modest timestep sizes. Finally, an implicit orbit averaging and sub-cycling scheme for the fully kinetic ion model is considered.
  • Exact collisional plasma fluid theories
    Eero Hirvijoki/David Pfefferlé, PPPL, abstract, slides
    [#s129, 10 Feb 2017]
    Following Grad’s procedure, an expansion of the velocity space distribution functions in terms of multi-index Hermite polynomials is carried out to derive a consistent set of collisional fluid equations for plasmas. The velocity-space moments of the often troublesome nonlinear Landau collision operator are evaluated exactly, and to all orders with respect to the expansion. The collisional moments are shown to be generated by applying gradients on two well-known functions, namely the Rosenbluth-MacDonald-Judd-Trubnikov potentials for a Gaussian distribution. The expansion can be truncated at arbitrary order with quantifiable error, providing a consistent and systematic alternative to the Chapman-Enskog procedure which, in plasma physics, amounts to the famous Braginskii equations. To illustrate our approach, we provide the collisional ten-moment equations and prove explicitly that the exact, nonlinear expressions for the momentum- and energy-transfer rate satisfy the correct conservation properties.
  • Recent advances in the variational formulation for reduced Vlasov-Maxwell equations
    Alain Brizard, Saint Michael’s College, Colchester, USA, abstract, slides
    [#s126, 02 Feb 2017]
    The talk presents recent advances in the variational formulation of reduced Vlasov-Maxwell equations. First, the variational formulations of guiding-center Vlasov-Maxwell theory based on Lagrange, Euler, and Euler-Poincaré variational principles are presented. Each variational principle yields a different approach to deriving guiding-center polarization and magnetization effects into the guiding-center Maxwell equations. The conservation laws of energy, momentum, and angular momentum are also derived by Noether method, where the guiding-center stress tensor is now shown to be explicitly symmetric. Next, the Eulerian variational principle for the nonlinear electromagnetic gyrokinetic Vlasov-Maxwell equations is presented in the parallel-symplectic representation, where the gyrocenter Poisson bracket contains contributions from the perturbed magnetic field.
  • Discrete Exterior Calculus Discretization of the Navier-Stokes Equations
    Ravindra Samtaney, King Abdullah University Sci. & Technology, abstract, slides
    [#s127, 10 Jan 2017]
    A conservative discretization of incompressible Navier-Stokes equations over surface simplicial meshes is developed using discrete exterior calculus (DEC). The DEC discretization is carried out for the exterior calculus form of Navier-Stokes equations, where the velocity field is represented by a 1-form. A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme for structured-triangular meshes, and first order accuracy for general unstructured meshes. The mimetic character of many of the DEC operators provides exact conservation of both mass and vorticity, in addition to superior kinetic energy conservation. The employment of various discrete Hodge star definitions based on both circumcentric and barycentric dual meshes is also demonstrated. The barycentric Hodge star allows the discretization to admit arbitrary simplicial meshes instead of being limited only to Delaunay meshes, as in previous DEC-based discretizations. The convergence order attained through the circumcentric Hodge operator is retained when using the barycentric Hodge. The discretization scheme is presented in detail along with numerical test cases demonstrating its numerical convergence and conservation properties. Preliminary results regarding the implementation of hybrid (circumcentric/barycentric) Hodge star operator are also presented. We conclude with some ideas for employing a similar method for magnetohydrodynamics.
  • Global electromagnetic gyrokinetic and hybrid simulations of Alfvén eigenmodes
    Michael Cole, Max Planck Institute for Plasma Physics, Greifswald , abstract, slides
    [#s125, 09 Dec 2016]
    The pursuit of commercial fusion power has driven the development of increasingly complex and complete numerical simulation tools in plasma physics. Recent work with the EUTERPE particle-in-cell code has made possible global, electromagnetic, fully gyrokinetic and fluid-gyrokinetic hybrid simulations in a broad parameter space, where previously global gyrokinetic simulations had been hampered by the so-called ‘cancellation problem’. This has been applied to the simulation of the interaction between Alfvén eigenmodes and energetic particles. In this talk, the range of numerical methods used will be detailed, and it will be shown with practical examples that self-consistent global simulations may be necessary for even a qualitatively accurate prediction of the perturbation of the magnetic field and fast particle transport due to wave-particle interaction. A brief outline will be given of the future direction of this work, such as the possibility of gyrokinetic simulation of the interaction between fine-scale turbulence and MHD modes.
  • Sparse grid techniques for particle-in-cell schemes
    Lee Ricketson, Courant Institute, New York University , abstract, slides
    [#s124, 21 Nov 2016]
    The particle-in-cell (PIC) method has long been the standard technique for kinetic plasma simulation across many applications. The downside, though, is that quantitatively accurate, 3-D simulations require vast computing resources. A prominent reason for this complexity is that the statistical figure of merit is the number of particles per cell. In 3-D, the number of cells grows rapidly with grid resolution, necessitating an astronomical number of particles. To address this challenge, we propose the use of sparse grids: by a clever combination of the results from a variety of grids, each of which is well resolved in at most one coordinate direction, we achieve similar accuracy to that of a full grid, but with far fewer grid cells, thereby dramatically reducing the statistical error. We present results from test cases that demonstrate the new scheme's accuracy and efficiency. We also discuss the limitations of the approach and, in particular, its need for an intelligent choice of coordinate system.
  • Optical collapse and nonlinear laser beam combining
    Pavel Lushnikov, U. New Mexico , abstract, slides
    [#s121, 18 Nov 2016]
    Many nonlinear systems of partial differential equations admit spontaneous formation of singularities in a finite time (blow up). Blow up is often accompanied by a dramatic contraction of the spatial extent of solution, which is called by collapse. A collapse in a nonlinear Schrodinger equation (NLSE) describes the self-focusing of the intense laser beam in the nonlinear Kerr medium (like usual glass) with the propagation distance $z$ playing the role of time. NLSE in the dimension two (two transverse coordinates) corresponds the stationary self-focusing of the laser beam eventually causing optical damage as was routinely observed in experiment since 1960s. NLSE in the dimension three (two transverse coordinates and time) is responsible for the formation of the optical bullet making pulse much shorter in time in addition to the spatial self-focusing. We address the universal self-similar scaling near collapse. In the critical 2D case the collapsing solutions have a form of rescaled soliton such that the $z$-dependence of that scale determines the $z$-dependent collapse width $L(z)$ and amplitude $\sim 1/L(z)$. At the leading order $L(z) \sim (z_c-z)^{1/2}$, where $z_c$ is the collapse time with the required log-log modification of that scaling. Log-log scaling for NLSE was first obtained asymptotically in 1980s and later proven in 2006. However, it remained a puzzle that this scaling was never clearly observed in simulations or experiment. We found that the classical log-log modification NLSE requires double-exponentially large amplitudes of the solution $\sim 10^{10^{100}}$, which is unrealistic to achieve in either physical experiments or numerical simulations. In contrast, we developed a new asymptotic theory which is valid starting from quite moderate (about 3 fold) increase of the solution amplitude compare with the initial conditions. We use that new theory to propose a nonlinear combining of multiple laser beams into a diffraction-limited beam by beam self-focusing in Kerr medium. Multiple beams with total power above critical are combined in near field and propagated through multimode optical fiber. Random fluctuations during propagation first trigger the formation of the strong optical turbulence. During subsequent propagation, the inverse cascade of optical turbulence tends to increase the transverse spatial scale of fluctuation until it efficiently triggers a strong optical collapse event producing diffraction-limited beam with the critical power.
  • On Degenerate Lagrangians, Noncanonical Hamiltonians, Dirac Constraints and their Discretization
    Michael Kraus, Max Planck Institute for Plasma Physics, Garching , abstract, slides
    [#s123, 17 Nov 2016]
    Most systems encountered in plasma physics are Hamiltonian and therefore have a rich geometric structure, most importantly symplecticity and conservation of momentum maps. As most of these systems are formulated in noncanonical coordinates, they are not amenable to standard symplectic discretisation methods, which are popular for the integration of canonical Hamiltonian systems.
    Variational integrators, which can be seen as the Lagrangian equivalent to symplectic methods, seem to provide an alternative route towards the systematic derivation of structure-preserving numerical methods for such systems. However, for noncanonical Hamiltonian systems the corresponding Lagrangian is often found to be degenerate. This degeneracy gives rise to instabilities of the variational integrators which need to be overcome in order to make long-time simulations possible.
    In this talk, recent attempts to devise long-time stable structure-preserving integrators for noncanonical Hamiltonian and degenerate Lagrangian systems will be reviewed. The guiding-centre system will be used to exemplify the problems which arise for such systems and to demonstrate the good long-time fidelity of the newly developed integrators.
  • Development and applications of Verification and Validation procedures
    Fabio Riva, EPFL, Switzerland , abstract
    [#s110, 14 Nov 2016]
    The methodology used to assess the reliability of numerical simulation codes constitutes the Verification and Validation (V&V) procedure. V&V is composed by three separate tasks: the code verification, which is a mathematical issue targeted to assess that the physical model is correctly implemented in a simulation code; the solution verification, which evaluates the numerical errors affecting a simulation; and the validation, which determines the consistency of the code results, and therefore of the physical model, with experimental data.
    To perform a code verification, we propose to use the method of manufactured solutions, a methodology that we have generalized to PIC codes, overcoming the difficulty of dealing with a numerical method intrinsically affected by statistical noise. The solution verification procedure we put forward is based on the Richardson extrapolation, used as higher order estimate of the exact solution. These verification procedures were applied to GBS, a three-dimensional fluid code for SOL plasma turbulence simulation based on a finite difference scheme, and to a unidimensional, electrostatic, collisionless PIC code. To perform a detailed validation of GBS against experimental measurements, we generalized the magnetic geometry of the simulation code to include elongation and non-zero triangularity, and we investigated theoretically the impact of plasma shaping effects on SOL turbulence. An experimental campaign is now planned on TCV, to validate our findings against experimental measurements in tokamak limited configurations.
  • Modeling efforts in hybrid kinetic-MHD and fully kinetic theories
    Cesare Tronci, University of Surrey, UK , abstract, slides
    [#s111, 28 Oct 2016]
    Over the decades, multiscale modeling efforts have resorted to powerful methods, such as asymptotic/perturbative expansions and/or averaging techniques. As a result of these procedures, finer scale terms are typically discarded in the fundamental equations of motion. Although this process has led to well consolidated plasma models, consistency issues may emerge in certain cases especially concerning the energy balance. This may lead to the presence of spurious instabilities that are produced by nonphysical energy sources. The talk proposes alternative techniques based on classical mechanics and its underlying geometric principles. Inspired by Littlejohn's guiding-center theory, the main idea is to apply physical approximations to the action principle (or the Hamiltonian structure) underlying the fundamental system, rather than operating directly on its equations of motion. Here, I will show how this method provides new energy-conserving variants of hybrid kinetic-MHD models, which suppress the spurious instabilities emerging in previous non-conservative schemes. Also, this method allows for quasi-neutral approximations of fully kinetic Vlasov theories, thereby neglecting both radiation and Langmuir oscillations.
  • Extending geometrical optics: A Lagrangian theory for vector waves
    Daniel Ruiz, Princeton University , abstract
    [#s114, 27 Oct 2016]
    Even diffraction aside, the commonly known equations of geometrical optics (GO) are not entirely accurate. GO considers wave rays as classical particles, which are completely described by their coordinates and momenta, but rays have another degree of freedom, namely, polarization. As a result, wave rays can behave as particles with spin. A well-known example of polarization dynamics is wave-mode conversion, which can be interpreted as rotation of the (classical) ``wave spin.'' However, there are other less-known manifestations of the wave spin, such as polarization precession and polarization-driven bending of ray trajectories. This talk presents recent advances in extending and reformulating GO as a first-principle Lagrangian theory, whose effective-gauge Hamiltonian governs both mentioned polarization phenomena simultaneously. Examples and numerical results are presented. When applied to classical waves, the theory correctly predicts the polarization-driven divergence of left- and right- polarized electromagnetic waves in isotropic media, such as dielectrics and nonmagnetized plasmas. In the case of particles with spin, the formalism also yields a point-particle Lagrangian model for the Dirac electron, i.e. the relativistic spin-1/2 electron, which includes both the Stern-Gerlach spin potential and the Bargmann-Michel-Telegdi spin precession. Additionally, the same theory contributes, perhaps unexpectedly, to the understanding of ponderomotive effects in both wave and particle dynamics; e.g., the formalism allows to obtain the ponderomotive Hamiltonian for a Dirac electron interacting with an arbitrarily large electromagnetic laser field with spin effects included.
  • Understanding and Predicting Profile Structure and Parametric Scaling of Intrinsic Rotation
    Weixing Wang, PPPL , abstract
    [#s120, 27 Oct 2016]
    This talk reports on a recent advancement in developing physical understanding and a first-principles-based model for predicting intrinsic rotation profiles in magnetic fusion experiments, including ITER. It is shown for the first time that turbulent fluctuation-driven residual stress (a non-diffusive component of momentum flux) can account for both the shape and magnitude of the observed intrinsic toroidal rotation profile. The model predictions of core rotation based on global gyrokinetic simulations agree well with the experimental measurements for a set of DIII-D ECH discharges. The characteristic dependence of residual stress and intrinsic rotation profile structure on the multi-dimensional parametric space covering turbulence type, q-profile structure, collisionality and up-down asymmetry in magnetic geometry has been studied with the goal of developing physics understanding needed for rotation profile control and optimization. Finally, the first-principles-based model is applied to elucidating the ρ∗-scaling and predicting rotations in ITER regime.
  • Laser-Driven Magnetized Collisionless Shocks
    Derek Schaeffer, Princeton University , abstract
    [#s113, 14 Oct 2016]
    Collisionless shocks -- supersonic plasma flows in which the interaction length scale is much shorter than the collisional mean free path -- are common phenomena in space and astrophysical systems, including the solar wind, coronal mass ejections, supernovae remnants, and the jets of active galactic nuclei. These systems have been studied for decades, and in many the shocks are believed to efficiently accelerate particles to some of the highest observed energies. Only recently, however, have laser and diagnostic capabilities evolved sufficiently to allow the detailed study in the laboratory of the microphysics of collisionless shocks over a large parameter regime. We present experiments that demonstrate the formation of collisionless shocks utilizing the Phoenix laser laboratory and the LArge Plasma Device (LAPD) at UCLA. We also show recent observations of magnetized collisionless shocks on the Omega EP laser facility that extend the LAPD results to higher laser energy, background magnetic field, and ambient plasma density, and that may be relevant to recent experiments on strongly driven magnetic reconnection. Lastly, we discuss a new experimental regime for shocks with results from high-repetition (1 Hz), volumetric laser-driven measurements on the LAPD. These large parameter scales allow us to probe the formation physics of collisionless shocks over several Alfvenic Mach numbers ($M_A$), from shock precursors (magnetosonic solitons with $M_A<1$) to subcritical ($M_A<3$) and supercritical ($M_A>3$) shocks. The results show that collisionless shocks can be generated using a laser-driven magnetic piston, and agree well with both 2D and 3D hybrid and PIC simulations. Additionally, using radiation-hydrodynamic modeling and measurements from multiple diagnostics, the different shock regimes are characterized with dimensionless formation parameters, allowing us to place disparate experiments in a common and predictive framework.
  • Radiation effects on the runaway electron avalanche
    Chang Liu, Princeton University , abstract, slides
    [#s118, 14 Oct 2016]
    Runaway electrons are a critical area of research into tokamak disruptions. A thermal quench on ITER can result in avalanche production of a large amount of runaway electrons and a transfer of the plasma current to be carried by runaway electrons. The potential damage caused by the highly energetic electron beam poses a significant challenge for ITER to achieve its mission. It is therefore extremely important to have a quantitative understanding of the runaway electron avalanche process. It is found that the radiative energy loss and the pitch angle scattering from radiative E&M fields plays an important role in determining the runaway electron distribution in momentum space. In this talk we discuss three kinds of radiation from runaway electrons, synchrotron radiation, Cerenkov radiation, and electron cyclotron emission (ECE) radiation. Synchrotron radiation, which mainly comes from the cyclotron motion of highly relativistic runaway electrons, dominates the energy loss of runaway electrons in the high-energy regime. The Cerenkov radiation from runaway electrons gives an additional correction to the Coulomb logarithm in the collision operator, which changes the avalanche growth rate. The ECE emission mainly comes from electrons in the energy range 1.2<γ<3, and gives an important approach to diagnose the runaway electron distribution in momentum and pitch angle. We developed a novel tool to self-consistently calculate normal mode scattering of runaway electrons using the quasi-linear method, and implement that in the a well-developed runaway electron kinetic simulation code CODE. Using this we successfully reproduce the experimental result of ECE signal qualitatively.
  • Plasmoids formation in a laboratory and large-volume flux closure during simulations of Coaxial Helicity Injection in NSTX-U
    Fatima Ebrahimi, PPPL and Princeton University , abstract
    [#s117, 14 Oct 2016]
    In NSTX-U, transient Coaxial Helicity Injection (CHI) is the primary method for current generation without reliance on the solenoid. A CHI discharge is generated by driving current along open field lines (the injector flux) that connect the inner and outer divertor plates on NSTX/NSTX-U, and has generated over 200 kA of toroidal current on closed flux surfaces in NSTX. Extrapolation of the concept to larger devices requires an improved understanding of the physics of flux closure and the governing parameters that maximizes the fraction of injected flux that is converted to useful closed flux. Here, through comprehensive resistive MHD NIMROD simulations conducted for the NSTX and NSTX-U geometries, two new major findings will be reported. First, formation of an elongated Sweet-Parker current sheet and a transition to plasmoid instability has for the first time been demonstrated by realistic global simulations [1]. This is the first observation of plasmoid instability in a laboratory device configuration predicted by realistic MHD simulations and then supported by experimental camera images from NSTX. Second, simulations have now, for the first time, been able to show large fraction conversion of injected open flux to closed flux in the NSTX-U geometry [2]. Consistent with the experiment, simulations also show that reconnection could occur at every stage of the helicity injection phase. The influence of 3D effects, and the parameter range that supports these important new findings is now being studied to understand the impact of toroidal magnetic field and the electron temperature, both of which are projected to increase in larger ST devices.
    [1] F. Ebrahimi & R. Raman, Phys. Rev. Lett. 114, 205003 (2015)
    [2] F. Ebrahimi & R. Raman, Nucl. Fusion 56, 044002 (2016)
  • Generation of helium and oxygen EMIC waves by the bunch distribution of oxygen ions associated with weak fast magnetosonic shocks in the magnetosphere
    Lou-Chuang Lee, Academia Sinica, Taiwan , abstract, slides
    [#s119, 11 Oct 2016]
    Electromagnetic ion cyclotron (EMIC) waves are often observed in the magnetosphere with frequency usually in the proton and helium cyclotron bands and sometimes in the oxygen band. The temperature anisotropy, caused by injection of energetic ions or by compression of magnetosphere, can efficiently generate proton EMIC waves, but not as efficient for helium or oxygen EMIC waves. Here we propose a new generation mechanism for helium and oxygen EMIC waves associated with weak fast magnetosonic shocks, which are observed in the magnetosphere. These shocks can be associated with either dynamic pressure enhancement or shocks in the solar wind and can lead to the formation of a “bunch” distribution in the perpendicular velocity plane of oxygen ions. The oxygen bunch distribution can excite strong helium EMIC waves and weak oxygen and proton waves. The dominant helium EMIC waves are strong in quasi-perpendicular propagation and show harmonics in frequency spectrum of Fourier analysis. The proposed mechanism can explain the generation and some observed properties of helium and oxygen EMIC waves in the magnetosphere.
  • Penetration and amplification of resonant perturbations in 3D ideal-MHD equilibria
    Stuart Hudson, PPPL , abstract
    [#s112, 10 Oct 2016]
    The nature of ideal-MHD equilibria in three-dimensional geometry is profoundly affected by resonant surfaces, which beget a non-analytic dependence of the equilibrium on the boundary. Furthermore, non-physical currents arise in equilibria with continuously-nested magnetic surfaces and smooth pressure and rotational-transform profiles. We demonstrate that three-dimensional, ideal-MHD equilibria with nested surfaces and δ-function current densities that produce a discontinuous rotational-transform are well defined and can be computed both perturbatively and using fully-nonlinear equilibrium calculations. The results are of direct practical importance: we predict that resonant magnetic perturbations penetrate past the rational surface (i.e. “shielding” is incomplete, even in purely ideal-MHD) and that the perturbation is amplified by plasma pressure, increasingly so as stability limits are approached.
  • Gyrokinetic projection of the divertor heat-flux width from present tokamaks to ITER
    C.S. Chang, PPPL , abstract
    [#s115, 10 Oct 2016]
    The total-f edge gyrokinetic code XGC1 shows that the divertor heat-flux width $λ_q$ in three US tokamaks (DIII-D for conventional aspect ratio, NSTX for tight aspect ratio, and C-Mod for high BP) obeys the experimentally observed $λ_q\propto 1/B_P^\gamma$ scaling in the so called “sheath-dominant regime.” The low-beta edge plasma is non-thermal and approaches the quasi-steady state in a kinetic non-diffusive time scale. Nonlinear Fokker-Planck-Landau collision operator is used. Monte-Carlo neutral atoms are recycled near the material wall. Successful validation of the XGC1 simulation results on three US tokamak devices will be presented in the so called "sheath-limited" regime. It is found that $λ_q$ on DIII-D, NSTX, and lower-$B_P$ C-Mod is dominated by the neoclassical orbit dynamics of the supra-thermal ions. However, C-Mod at higher $B_P$ shows blob-dominance, while still fitting into the $λ_q\propto 1/B_P^\gamma$ graph. Predictive simulation on ITER shows that $λ_q$ is over 5 times greater than that predicted by the empirical $λ_q\propto 1/B_P^\gamma$ scaling.
  • Relativistic Electrons and Magnetic Reconnection in ITER
    Allen Boozer, Columbia University , abstract, slides
    [#s109, 22 Sep 2016]
    ITER, the largest scientific project ever undertaken, “has been designed to prove the feasibility of fusion” energy. That mission could be compromised if the net current in the plasma were transferred to relativistic-electron carriers, which can result in great damage to the device. More than a hundred-fifty papers on the topic have appeared in the twenty years since it was realized that a sudden drop in the plasma temperature could cause such a transfer. The theoretical papers have focused on electron runaway when magnetic field lines remain confined to the plasma volume. Experiments and their simulation imply most magnetic field lines intercept the walls when a sudden drop in the temperature occurs. In simulations not all of the magnetic field lines intercept the walls. If these non-intercepting flux tubes survive until outer magnetic surfaces reform, an especially dangerous situation arises in which the relativistic electrons can be lost in a short pulse along a narrow flux tube.
    Maxwell’s equations imply the conservation of magnetic helicity during fast magnetic reconnections. Helicity conservation clarifies features, such as the spike in the plasma current during a thermal quench, and shows that a rapid, ~1ms, acceleration of electrons can occur. Magnetic reconnection can occur on a time scale of tens of Alfvén transit times either due to the formation of multiple plasmoids or as a result of the exponential sensitivity to non-ideal effects when an ideal evolution causes neighboring magnetic field lines to increase their separation exponentially with distance along the lines. In either case, magnetic flux tubes that had different force-free parallel currents $j_{\parallel}/B$ and different poloidal magnetic fluxes external to their surfaces can be joined in a fast reconnection. The differing parallel currents relax by Alfvén waves, and the change in the poloidal flux associated with a given enclosed toroidal magnetic flux is given by helicity conservation. Based on the observed time scale of current spikes, a few hundred toroidal transits are required to relax the $j_{\parallel}/B$ profile, which is consistent with distance along magnetic field lines implied by the speed of the electron temperature drop. Helicity conservation implies that although changes in the poloidal flux can occur on a time scale determined by the Alfvén transit time, the magnetic helicity and most of the poloidal flux can change only on a long resistive, $L/\mathcal{R}$, time scale. Rapid plasma terminations require plasma cooling.
    The potential for damage, the magnitude of the extrapolation from existing devices, and the importance of the atypical—incidents that occur once in a thousand shots—make theory and simulation essential for ensuring that relativistic runaway electrons will not prevent ITER from achieving its mission. The U.S. DoE Office of Science has established a Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) with a two-year funding of 3.9M$\$$. The danger of runaway of electrons to relativistic energies can be avoided in magnetically confined fusion systems by making the net current the $j_{\parallel}/B$ sufficiently small. This is possible in the non-axisymmetric stellarator geometry though not in axisymmetric tokamaks, such as ITER.
  • Variational current-coupling gyrokinetic-MHD
    Josh Burby, Courant Institute, New York University , abstract, slides
    [#s105, 07 Sep 2016]
    In this talk I will describe the details and derivation of a new current-coupling gyrokinetic-MHD model. In particular, I will show that the model can be derived from a variational principle. Energy, and hot charge are conserved exactly regardless of the form of the background magnetic field. Likewise, when the background field admits a continuous rotation or translation symmetry, the corresponding component of the total momentum is conserved. The theory relies on a new gauge-invariant formulation of the motion of gyrocenters in prescribed electromagnetic fields, and this will be described in detail.
  • Implicit Multiscale Full Kinetics as an Alternative to Gyrokinetics
    Scott E. Parker, University of Colorado , abstract
    [#s103, 30 Aug 2016]
    Recent progress has been made developing full kinetic Lorentz force ion dynamics using implicit multiscale techniques [1]. It is now possible to capture low-frequency physics along with finite Larmor radius (FLR) effects with a fully kinetic multiscale $\delta f$ particle simulation. The utility of such a model is to be able to verify gyrokinetics in situations where the smallness of the ordering parameters is questionable. Additionally, such a model can help identify what higher order terms in gyrokinetics might be important. Orbit averaging and sub-cycling are utilized with an implicit particle time advance based on variational principles. This produces stable and accurate ion trajectories on long time scales. Excellent agreement with the gyrokinetic dispersion relation is obtained including full FLR effects. Ion Bernstein waves and the compressional Alfvén wave are easily suppressed with the implicit time advance. We have developed a global toroidal electrostatic adiabatic electron Lorentz ion code. We will report preliminary linear results benchmarking Lorentz ions with gyrokinetics for the Cyclone base case. We will begin by reviewing recent results from the GEM code simulating electromagnetic gyrokinetic turbulence in the edge pedestal where the timestep required is comparable to the ion cyclotron period.
    [1] Benjamin J. Sturdevant, Scott E. Parker et al., J. Comput. Phys. 316, 519 (2016)
  • Current status of the LHD and the prospect for Deuterium experiment
    Dr. Masaki Osakabe, National Institute for Fusion Science, Japan
    [#s106, 26 Aug 2016]
  • Dissipation and Intermittency in Gyrokinetic Turbulence and Beyond
    Jason TenBarge, University of Maryland , abstract, slides
    [#s102, 25 Aug 2016]
    Turbulence is a ubiquitous process in space and astrophysical plasmas that serves to mediate the transfer of large-scale motions to small scales at which the turbulence can be dissipated and the plasma heated. In situ solar wind observations and direct numerical simulations demonstrate that sub-proton scale turbulence is dominated by highly anisotropic and intermittent, low frequency, kinetic Alfvénic fluctuations. I will review recent work on the dissipation of Alfvénic turbulence observed in gyrokinetic simulations and discuss the coherent structures and intermittency associated with the turbulence, which suggest a non-local and non-self-similar energy cascade. Moving beyond the confines of gyrokinetics, I will also briefly discuss work on a full Eulerian Vlasov-Maxwell code, Gkeyll, being developed at Princeton and the University of Maryland.
  • Effective resistivity in collisionless magnetic reconnection
    Zhi-Wei Ma, Zhejiang University, China , abstract
    [#s80, 18 Aug 2016]
    The well-known, physical mechanism for fast, magnetic reconnection in collisionless plasmas is that the off-diagonal terms of the electron-pressure tensor give rise to a larger electric-fields in the reconnection region. The electron-pressure tensor fully associated with electron kinetic effects is difficultly implemented into the MHD model. In this talk, we try to use a simple equation $E = \eta J$ (where $\eta$ is an effective resistivity) to illustrate the fast reconnection in collisionless, magnetic reconnection. The physical mechanism and formulation of the effective resistivity are addressed.
  • Turbulence in shear MHD flows: Implications for accretion disks
    Farrukh Nauman, Niels Bohr International Academy, Copenhagen , abstract, slides
    [#s63, 04 Aug 2016]
    Accretion flows are found in a large variety of astrophysical systems, from protoplanetary disks to active galactic nuclei. Our present understanding of such flows is severely limited by both observational and numerical resolution. I will discuss some new numerical results on zero magnetic flux shear MHD turbulence and its relation to the magnetic Prandtl number. I will then briefly discuss the effects of rotation on large scale magnetic fields. My talk will end with some speculations about how one might construct a self-consistent model for accretion flows based on our current understanding.
  • Dynamics of the ELMs in pre-crash and crash suppressed period in KSTAR
    Hyeon Park, NFRI/UNIST, Korea , abstract, slides
    [#s77, 28 Jul 2016]
    Following the first operation of H-mode in KSTAR in 2009, study of the edge localized modes (ELM) has been actively conducted. A unique in-vessel control coil (IVCC) set (top, middle and bottom) capable of generating resonant (and non-resonant) magnetic perturbation (RMP) at low n(=1,2) number was successfully utilized to suppress and/or mitigate the ELM-crash in KSTAR. Extensive study of dynamics of the ELMs in both pre-crash and crash suppressed phase under magnetic perturbation with the 2D/3D Electron Cyclotron Emission Imaging (ECEI) system revealed new phenomenology of the ELMs and ELM-crash dynamics that were not available from conventional diagnostics. Since the first 2D images of the ELM time evolution from growth to crash through saturation, the detailed images of the ELMs leading to the crash together with the fast RF emission (<200MHz) signal demonstrated that the pre-crash events are complex. The measured 2D image of the ELM was validated by direct comparison with the synthetic 2D image by the BOUT++ code and non-linear modelling study is in progress. Recently, the observed dynamics of the ELMs at both high and low field sides such as asymmetries in intensity, mode number and rotation direction casted a doubt in peeling-ballooning mode. Response of high field side ELM to the RMP was more pronounced compared to that of the low field side. Other study includes observation of multi-modes and sudden mode number transition. During the ELM-crash suppression experiment, various types of ELM-crash patterns were observed and often the suppression was marginal. The observed semi-coherent turbulence spectra under the RMP provided an evidence of non-linear interaction between the ELMs and turbulence.
  • Statistical origin and properties of kappa distributions
    George Livadiotis, Southwest Research Institute , abstract, slides
    [#s62, 07 Jul 2016]
    Classical particle systems reside at thermal equilibrium with their velocity distribution function stabilized into a Maxwell distribution. On the contrary, collisionless and correlated particle systems, such as space and astrophysical plasmas, are characterized by a non-Maxwellian behavior, typically described by so-called $\kappa$ distributions, or combinations thereof. Empirical $\kappa$ distributions have become increasingly widespread across space and plasma physics. A breakthrough in the field came with the connection of $\kappa$ distributions to non-extensive statistical mechanics. Understanding the statistical origin of $\kappa$ distributions was the cornerstone of further theoretical developments and applications, some of which will be presented in this talk: (i) The physical meaning of thermal parameters, e.g., temperature and kappa index; (ii) the multi-particle description of $\kappa$ distributions; (iii) the generalization to phase-space $\kappa$ distribution of a Hamiltonian with non-zero potential; (iv) the Sackur-Tetrode entropy for $\kappa$ distributions, and (v) the existence of a large-scale phase-space cell, characteristic of collisionless space plasmas, indicating a new quantization constant, $\hbar ^* \sim 10^{-22} Js$.
  • Explosive Solution of a Time Delayed Nonlinear Cubic Equation Derived for Fluids (Hickernell) and Plasmas (Berk-Breizman)
    Herb Berk, Institute for Fusion Studies, U. Texas at Austin, abstract, slides
    [#s61, 27 Jun 2016]
    This presentation will describe new explosive attractor solutions to the universal cubic delay equation found in both the fluid [1] and (for a kinetic system) in the plasma literature [2]. The cubic delay equation describes a system governed by a control parameter $\phi$ (in plasmas its value is determined by the linear properties of the kinetic response). The simulation of the temporal evolution reveals the development of an explosive mode, i.e. a mode growing without bound in a finite time. The two main features of the response are: (1) a well-known explosive envelope $(t_0-t)^{-5/2}$, with $t_0$ the blow-up time of the amplitude; (2) a spectrum with ever-increasing oscillation frequencies whose values depend on the parameter $\phi$. Analytic modeling explains the results and quantitatively nearly replicates the attractor solutions found in the simulations. These analytic attractor solutions are linearly stable except in some cases where the nonlinear solution needs to be corrected to include higher harmonics. Our analysis explains almost all of the rather complicated numerical attractor solutions for the cubic delay equation.
    [1] F.J. Hickernell, J. Fluid Mech. 142, 431 (1984)
    [2] B.N. Breizman, H.L. Berk et al., Phys. Plasmas 4, 1559 (1997)
  • Physics of tokamak flow relaxation to equilibrium
    Bruce Scott, Max-Planck-IPP, EURATOM Association, abstract
    [#s59, 23 Jun 2016]
    The theorem for toroidal angular momentum conservation within gyrokinetic field theory is used as a starting point for consideration of flow equilibration at low frequencies (less than fast-Alfvén or gyrofrequencies). Quasineutrality and perpendicular MHD force balance are inputs to the theory and therefore never violated. However, the gyrocenter densities are not ambipolar in equilibrium, since the flow vorticity is given by their difference. From an arbitrary initial state, flows evolve acoustically and via Landau damping into divergence balance, in which radial force balance of the electric field is a part. On collisional time scales, which in the tokamak core are longer, the neoclassical electric field is brought into balance by collisions, and it is only on these slow time scales that the collisional transport is ambipolar (ie, the time derivative of the vorticity is small). Computations from 2014 showing the relaxation on tokamak core spatial scales are displayed. I will also give relevant cases of edge-layer relaxation and discuss the dependence on the finite poloidal gyroradius. Total-$f$ two-species gyrokinetic relaxation cases from 2009/10 are available to show that the basic processes in fluid and gyrokinetic models are the same for these purposes.
  • Equilibrium Potential Well due to Finite Larmor Radius Effects at the Tokamak Edge
    W. W. Lee, PPPL , abstract, slides
    [#s49, 16 Jun 2016]
    We present a novel mechanism for producing the equilibrium potential well near the edge of a tokamak. Briefly, because of the difference in gyroradii between electrons and ions, an equilibrium electrostatic potential is generated in the presence of spatial inhomogeneity of the background plasma, which, in turn, produces a well associated with the radial electric field, $E_r$, as observed at the edge of many tokamak experiments. We will show that this theoretically predicted $E_r$ field, which can be regarded as producing a long radial wave length zonal flow, agrees well with recent experimental measurements. The work is in collaboration with R. B. White [1].
    [1] W.W. Lee & R.B. White, PPPL Report 5254 (2016)
  • Predicting solar magnetic activity and its implications for global dynamo models
    Nishant Kumar Singh, KTH Royal Institute of Technology, abstract, slides
    [#s51, 09 Jun 2016]
    Using the solar surface mode, i.e. the $f$-mode, we attempt to predict the emergence of active regions (ARs) in the days before they can be seen in magnetograms. Our study is motivated by earlier numerical findings of Singh et al. [1], who showed that, in the presence of a nonuniform magnetic field which is concentrated a few scale heights below the surface, the $f$-mode fans out in the diagnostic $k$-$\omega$ diagram at high wavenumbers. Here we exploit this property using data from the Helioseismic and Magnetic Imager aboard the Solar Dynamics Observatory, and show for about six ARs that at large latitudinal wavenumbers (corresponding to horizontal scales of around 3000 km), the $f$-mode displays strengthening about two days prior to AR formation and thus provides a new precursor for AR formation. I will also discuss ways to isolate signals from newly forming ARs in a crowded environment where existing ones are expected to pollute the neighbouring patches. The idea that the $f$-mode is perturbed days before any visible magnetic activity occurs on the surface can be important in constraining dynamo models aiming at understanding the global magnetic activity of the Sun.
    [1] Nishant K. Singh, Harsha Raichur & Axel Brandenburg, arxiv.org/abs/1601.00629 (2014)
  • Statistical analysis of turbulent transport for flux driven toroidal plasmas
    Johan Anderson, Chalmers University of Technology, abstract, slides
    [#s26, 07 Jun 2016]
    During recent years an overwhelming body of evidence shows that the overall transport of heat and particles is, to a large part, caused by intermittency (or bursty events) related to coherent structures. A crucial question in plasma confinement is thus the prediction of the probability distribution functions (PDFs) of the transport due to these structures and of their formation. This work provides a theoretical interpretation of numerically generated PDFs of intermittent plasma transport events, as well as offering an explanation for elevated PDF tails of heat flux. Specifically, we analyse time traces of heat flux generated by global nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence by the GKNET software [1]. The simulation framework is global, flux-driven and considers adiabatic electrons. In the simulations, SOC type intermittent bursts are frequently observed and transport is often regulated by non-diffusive processes, thus the PDFs of e.g. heat flux are in general non-Gaussian with enhanced tails. A key finding of this study is that the intermittent process in the context of drift-wave turbulence appears to be independent of the specific modelling framework, opening the way to the prediction of its salient features. Although the same PDFs were previously found in local gyrokinetic simulations [2], there are some unique features present inherently coming from the global nature of the physics. The main part of this work consists in providing a theoretical interpretation of the PDFs of radial heat flux. The numerically generated time traces are processed with Box–Jenkins modelling in order to remove deterministic autocorrelations, thus retaining their stochastic parts only. These PDFs have been shown to agree very well with analytical predictions based on a fluid model, on applying the instanton method. In this talk, the theory and comparisons to the numerical work will be presented. The result points to a universality in the modelling of the intermittently stochastic process while the analytical theory offers predictive capability, extending the previous result to be globally applicable.
    [1] K. Imadera et al., 25th FEC, TH/P5-8 (2014)
    [2] Johan Anderson & Pavlos Xanthopoulos, Phys. Plasmas 17, 110702 (2010)
  • Parasitic Momentum Flux in the Tokamak Core
    Tim Stoltzfus-Dueck, PPPL, abstract, slides
    [#s50, 06 Jun 2016]
    A higher-order portion of the ${\bf E}\times {\bf B}$ drift causes an outward flux of co-current momentum when electrostatic potential energy is transferred to ion-parallel flows. The robust symmetry breaking follows from the free energy flow in phase space and does not depend on any assumed linear eigenmode structure. The resulting rotation peaking is counter-current and scales as temperature over plasma current. This peaking mechanism can only act when there are adequate fluctuations at low enough frequencies to excite ion parallel flows, which may explain some experimental observations related to rotation reversals
  • Experimental study of the role of electron pressure in fast magnetic reconnection with a guide field
    Will Fox, PPPL , abstract
    [#s25, 03 Jun 2016]
    Magnetic reconnection, the change of magnetic topology in the presence of plasma, is observed in space, laboratory, and enables the explosive energy release by plasma instabilities, as in solar flares or magnetospheric substorms, and the change in topology allows the rapid heat transport associated with sawtooth relaxation and self-organization in RFPs. In numerous environments, especially in toroidal confinement devices, reconnection proceeds in the presence of a net guide field. We report detailed laboratory observations in MRX of the structure of reconnection current sheets with a guide field regime in a two-fluid plasma regime (ion gyro-radius comparable to the current sheet width) . We observe experimentally for the first time the quadrupolar electron pressure variation in the ion-diffusion region, an analogue of the quadrupolar "Hall" magnetic fields in anti-parallel reconnection. The quadrupolar pressure perturbation was originally predicted by extended MHD simulation as essential to balancing the large parallel reconnection electric fields over the ion-scale current-sheet We observe that electron density variations dominate temperature variations and may provide a new diagnostic of reconnection with finite guide field for fusion experiments and spacecraft missions. We discuss consequences for force balance in the reconnection layer and implications for fast reconnection in fusion devices.
  • 2D Full-wave simulations of plasma waves in space and tokamak plasmas
    Eun-Hwa Kim, PPPL, abstract, slides
    [#s12, 28 Apr 2016]
    A 2D full-wave simulation code (so-called FW2D) has been recently developed. This code currently solves the cold plasma wave equations using the finite element method. One advantage of using the finite element method is that the local basis functions can be readily adapted to boundary shapes and can be packed in such a way as to provide higher resolution in regions where needed. We have constructed a 2D triangular mesh given a specified boundary and a target mesh density function. Moreover, the density of the mesh can be specified based on the expected wavelength obtained from solution of the local dispersion (except close to resonances) so that the most efficient resolution is used. Another advantage of this wave code is short running time. For instance, by using node number of 24,395, the computing time is approximately 300 seconds CPU time to obtain a solution. The wave code has been successfully applied to describe low frequency waves in Earth and Mercury's multi-ion magnetospheres. The results include (a) mode conversion from the incoming fast to the transverse wave modes at the ion-ion hybrid resonance, (b) mode coupling and polarization reversal between left-handed (i.e., electromagnetic ion cyclotron waves: EMIC waves) and right-handed polarized waves (i.e., fast mode), and (c) refraction and reflection of field-aligned propagating EMIC waves near the heavier ion cyclotron frequency. Very recently FW2D code has been improved to adopt tokamak geometry and examine radio frequency (RF) waves in the scape-off layer (SOL) of tokamaks, which is the region of the plasma between last closed flux surface and tokamak vessel. The SOL region is important for RF wave heating of tokamaks because significant wave power loss can occur in this region. This code is ideal for waves in SOL plasma, because realistic boundary shapes and arbitrary density structures can be easily adopted in the code and the SOL plasma can be approximated as cold plasma.
  • Parallel electron force balance and the L-H transition
    Tim Stoltzfus-Dueck, Princeton University, abstract, slides
    [#s13, 14 Apr 2016]
    In a popular description of the L-H transition, energy transfer to the mean flows directly depletes kinetic energy from turbulent fluctuations, resulting in suppression of the turbulence and a corresponding transport bifurcation. However, electron parallel force balance couples non-zonal velocity fluctuations with electron pressure fluctuations on rapid timescales, comparable with the electron transit time. For this reason, energy in the non-zonal velocity stays in a fairly fixed ratio to electron thermal free energy, at least for frequency scales much slower than electron transit. In order for direct depletion of the energy in turbulent fluctuations to cause the L-H transition, energy transfer via Reynolds stress must therefore drain enough energy to significantly reduce the sum of the free energy in non-zonal velocities and electron pressure fluctuations. At low $k$, the electron thermal free energy is much larger than the energy in non-zonal velocities, posing a stark challenge for this model of the L-H transition.
  • Realistic characterizations of chirping instabilities in tokamaks
    Vinícius Duarte, PPPL, abstract, slides
    [#s14, 31 Mar 2016]
    In tokamak plasmas, the dynamics of phase-space structures with their associated frequency chirping is a topic of major interest in connection with mechanisms for fast ion losses. The onset of phase-space holes and clumps which produce chirping phenomena has been theoretically shown to be related to the emergence of an explosive solution of an integro-differential, nonlocal cubic equation [1,2] that governs the early mode amplitude evolution in the nonlinear regime near marginal stability. We have extended the analysis of the model to quantitatively account for multiple resonance surfaces of a given mode in the presence of drag and diffusion (due to collisions and micro-turbulence) operators. Then a more realistic criterion is found that takes into account the details of the mode structure and the variation of transport coefficients in phase space, to determine whether steady-state solutions can or cannot exist. Stable, steady-state solutions indicate that chirping oscillations do not arise, while the lack of steady solutions due to the predominance of drag is indicative that a frequency chirping response is likely in a plasma. Waves measured in experiments have been analyzed using the NOVA and NOVA-K codes, with which we can realistically account for the mode structure and varying resonance responses spread over phase space. In the experiments presently analyzed, compatibility has been found between the theoretical predictions for whether chirping should or should not arise and the experimental observation or lack of observation of toroidicity-induced Alfvén eigenmodes in NSTX, DIII-D and TFTR. We have found that stochastic diffusion due to wave micro-turbulence is the dominant energetic particle transport mechanism in many plasma experiments, and its strength is the key as to whether chirping solutions are likely to arise.
    [1] H.L. Berk, B.N. Breizman & M. Pekker, Phys. Rev. Lett. 76, 1256 (1996)
    [2] M.K. Lilley, B.N. Breizman & S.E. Sharapov, Phys. Rev. Lett. 102, 195003 (2009)
  • Speed-Limited Particle-in-Cell (SLPIC) Method
    John R. Cary, U. Colorado, Boulder, abstract, slides
    [#s15, 25 Mar 2016]
    The Speed-Limited Particle-In-Cell (SLPIC) Method reduces computational requirements for simulations that evolve on ion time scales while keeping appropriate kinetic electron effects. This method works by introducing an ansatz for the distribution function that allows the new unknown phase-space function to be solved by the method of characteristics, where those characteristics move slowly through phase space. Therefore, the solution can be obtained by particle-in-cell (PIC) methods, where the electrons have speeds much smaller than their actual speeds, ultimately leading to a much relaxed numerical (CFL) stability condition on the time step. Ultimately, the time step can be increased by the square root of the ion-electron mass ratio. SLPIC can be easily implemented in existing PIC codes as it requires no changes to deposition and field solution. Its explicit nature makes it ideal for modern computing architectures with vector instruction sets.
  • Free-Boundary Axisymmetric Plasma Equilibria: Computational Methods and Applications
    Holger Heumann, INRIA, France, abstract, slides
    [#s16, 03 Mar 2016]
    We present a comprehensive survey of the various computational methods for finding axisymmetric plasma equilibria. Our focus is on free-boundary plasma equilibria, where either poloidal field coil currents or the temporal evolution of voltages in poloidal field circuit systems are given data. Centered around a piecewise linear finite element representation of the poloidal flux map, our approach allows in large parts the use of established numerical schemes. The coupling of a finite element method and a boundary element method gives consistent numerical solutions for equilibrium problems in unbounded domains. We formulate a Newton-type method for the discretized non-linear problem to tackle the various non-linearities, including the free plasma boundary. The Newton method guarantees fast convergence and is the main building block for the inverse equilibrium problems that we discuss as well. The inverse problems aim at finding either poloidal field coil currents that ensure a desired shape and position of the plasma or at finding the evolution of the voltages in the poloidal field circuit systems that ensure a prescribed evolution of the plasma shape and position. We provide equilibrium simulations for the tokamaks ITER and WEST to illustrate performance and application areas.
  • Mixed finite-element/finite difference method for toroidal field-aligned elliptic electromagnetic equations
    Salomon Janhunen, PPPL, abstract
    [#s17, 29 Feb 2016]
    Gyrokinetic simulations -- such as those performed by the XGC code -- provide a self-consistent framework to investigate a wide range of physics in strongly magnetized high temperature laboratory plasmas, global modes usually considered to be in the realm of MHD simulations. However, the present simulation models generally concentrate on short wavelength electro-magnetic modes mostly to convenience the field solver performance. To incorporate more global fluid-like modes, also non-zonal long wavelength physics needs to be retained. In this work we present development of a fully 3D mixed FEM/FDM electro-magnetic field solver for use in the gyrokinetic code XGC1. We present optimization for use on massively parallel computational platforms, investigation of numerical accuracy characteristics using the method of manufactured solutions, and evaluate importance of field line length calculations on the stability of the discretization. We also invite discussion on the importance of the perpendicular vector potential for pressure driven modes.
  • Gyrokinetic Particle Simulation of Fast Electron Driven Beta-induced Alfvén Eigenmodes
    Wenlu Zhang, Chinese Academy of Science & U.C. Irvine, abstract
    [#s18, 25 Feb 2016]
    The fast electron driven beta induced Alfvén eigenmode (e-BAE) in toroidal plasmas is investigated for the first time using global gyrokinetic particle simulations, where the fast electrons are described by the drift kinetic model. The phase space structure shows that only the processional resonance is responsible for the e-BAE excitations while fast-ion driven BAE can be excited through all the channels such as transit, drift-bounce, and processional resonance. Frequency chirping is observed in nonlinear simulations with both weak and strong drives in the absence of sources and sinks, which provide a complement to the standard `bump-on-tail` paradigm for the frequency chirping of Alfvén eigenmodes. For weakly nonlinear driven case, frequency is observed to be in phase with the particle energy flux, and mode structure is almost the same as linear stage. While in the strongly driven nonlinear case, BAAE is excited along with BAE after the BAE mode saturated. Analysis of nonlinear wave-particle interactions show that the frequency chirping is induced by the nonlinear evolution of the coherent structures in the energetic-particle phase space, where the dynamics of the coherent structure is controlled by the formation and destruction of phrase space islands of energetic particles in the canonical variables. Zonal flow and zonal field are found to affect wave-particle resonance in the nonlinear e-BAE simulations.
  • Theory for Transport Properties of Warm Dense Matter
    Scott Baalrud, University of Iowa, abstract, slides
    [#s5, 18 Feb 2016]
    Progress in a number of research frontiers relies upon an accurate description of the transport coefficients of warm and hot dense matter, characterized by densities near those of solids and temperatures ranging from several to hundreds of eV. Examples include inertial confinement fusion, evolution of giant planets, exoplanets, and other compact astrophysical objects such as white dwarf stars, as well as numerous high energy density laboratory experiments. These conditions are too dense for standard plasma theories to apply and too hot for condensed matter theories to apply. The challenge is to account for the combined effects of strong Coulomb coupling of ions and quantum degeneracy of electrons. This seminar will discuss the first theory to provide fast and accurate predictions of ionic transport coefficients in this regime. The approach combines two recent developments. One is the effective potential theory (EPT), which is a physically motivated approach to extend plasma kinetic theory into the strong coupling regime. The second is a new average atom model, which provides accurate radial density distributions at high-density conditions, accounting for effects such as pressure ionization. Results are compared with state-of-the-art orbital-free density functional theory computations, revealing that the theory is accurate from high temperature through the warm dense matter regime, breaking down when the system exhibits liquid-like behaviors. A number of properties are considered, including diffusion, viscosity and thermal conductivity.
  • Nonlinear Fishbone Dynamics in Spherical Tokamaks with Toroidal Rotation
    Feng Wang, PPPL, abstract, slides
    [#s19, 11 Feb 2016]
    Fishbone is one of the most important energetic particles driven modes in tokamaks. A numerical study of the nonlinear dynamics of fishbone has been carried out in this work. Realistic parameters with finite toroidal plasma rotation are used to understand nonlinear frequency chirping in NSTX. We have carried out a systematic study of nonlinear frequency chirping and energetic particle dynamics. It is found that, linearly, the mode is driven by both trapped particles and passing particles, with resonance condition $\omega_{d} \simeq \omega$ for trapped particles and $\omega_{\phi}+\omega_{\theta}\simeq\omega$ for passing particles, where $\omega_{d}$ is trapped particle toroidal precession frequency, and $\omega_{\phi}$, $\omega_{\theta}$ are passing particle transit frequency in toroidal and poloidal direction. As the mode growing, trapped resonant particles oscillate and move outward radially, which reduces particle precessional frequency. We believe this is the main reason for the mode frequency chirping down. Finally, as the mode frequency chirping down, initially non-resonant particles with lower precessional frequencies become resonant particles in the nonlinear regime. This effect can sustain a quasi-steady state mode amplitude observed in the simulation.
  • Generation of anomalously energetic suprathermal electrons due to collisionless interaction of an electron beam with a nonuniform plasma
    Igor Kaganovich, PPPL, abstract, slides
    [#s4, 04 Feb 2016]
    Electrons emitted by electrodes surrounding or immersed in the plasma are accelerated by the sheath electric field and become the electron beams penetrating the plasma. Recently, it was reported that an electron energy distribution measured in a dc-rf discharge with 800V dc voltage has not only a peak at 800eV corresponding to the electrons emitted from the dc-biased electrode, but also a significant fraction of accelerated electrons with energy up to several hundred eV. The particle-in-cell simulation results show that the acceleration may be caused by the effects related to the plasma nonuniformity. The electron beam excites plasma waves whose wavelength and phase speed gradually decrease towards anode. The short waves near the anode accelerate plasma bulk electrons to suprathermal energies because of multiple interactions of electron with wave region [1]. The two-stream instability of an electron beam propagating in finite-size plasma placed between two electrodes was also studied analytically. It is shown that the growth rate in such a system is much smaller than that of infinite plasma or finite size plasma with periodic boundary conditions. We show that even if width of the plasma matches the resonance condition for standing waves; standing waves do not develop and transform into spatially growing wave, whose growth rate is small compared to that of the standing wave in a system with periodic boundary conditions. The frequency and growth rate as a function of plasma width form a bandwidth structure [2].
    [1] D. Sydorenko, I.D. Kaganovich et al., "Generation of anomalously energetic suprathermal electrons by an electron beam interacting with a nonuniform plasma", arXiv:1503.05048 (2015)
    [2] I.D. Kaganovich & D.Sydorenko, "Band Structure of the Growth Rate of the Two-Stream Instability of an Electron Beam Propagating in a Bounded Plasma" arXiv:1503.04695 (2015).
  • Nonlinear Alfvén waves in generalized magnetohydrodynamics - a creation on Casimir leaves
    Hamdi M. Abdelhamid, University of Tokyo, abstract
    [#s3, 03 Feb 2016]
    Alfvén waves are the most typical electromagnetic phenomena in magnetized plasma. In particular, nonlinear Alfvén waves deeply affect various plasma regimes in laboratory as well as in space, which have a crucial role in plasma heating, turbulence, etc. Large-amplitude Alfvén waves are observed in various systems in space and laboratories, demonstrating an interesting property that the wave shapes are stable even in the nonlinear regime. The ideal magnetohydrodynamics (MHD) model predicts that an Alfvén wave keeps an arbitrary shape constant when it propagates on a homogeneous ambient magnetic field.
    Here we investigate the underlying mechanism invoking a more accurate framework, generalized MHD. When we take into account dispersion effects (we consider both ion and electron inertial effects), the wave forms are no longer arbitrary. Interestingly, these "small-scale effects" change the picture completely; the large-scale component of the wave cannot be independent of the small scale component, and the coexistence of them forbids the large scale component to have a free wave form. This is a manifestation of the nonlinearity-dispersion interplay, which is somewhat different from that of solitons. The Casimir invariants of the system is the root cause of this interesting property.
  • Runaway Mitigation Issues on ITER
    Allen Boozer, Columbia University, abstract, slides
    [#s2, 25 Jan 2016]
    The plasma current in ITER can be transferred from near thermal to relativistic electrons by the runaway phenomenon. If such a current of relativistic electrons were to strike the chamber walls in ITER, the machine could be out of commission for many months. For ITER to be operable as a research device, the shortest credible time between such events must be years. The physics of the runaway process is remarkably simple and clear. The major uncertainty is what range of plasma conditions may arise in post thermal quench ITER plasmas. Consequently, a focused effort that includes theory, experiments, and engineering could relatively quickly clarify whether ITER will be operable with the envisioned mitigation strategy and what mitigation strategies could enhance the credibility that ITER will be operable.
  • Non-twist map bifurcation of drift-lines and drift-island formation in saturated 3D MHD equilibria
    David Pfefferlé, PPPL, abstract
    [#s1, 21 Jan 2016]
    Based on non-canonical perturbation theory of the field-line action [1], guiding-centre drift equations are identified as perturbed magnetic field-line equations. In this context, passing-particle orbits are called drift-lines, and their topology is completely determined by the magnetic configuration. In axisymmetric tokamak fields, drift-lines lie on shifted flux-surfaces, i.e. drift-surfaces, the shift being proportional at lowest order to the parallel gyro-radius and the q-profile [2]. Field-lines as well as drift-lines produce island structures at rational surfaces [3] only when a non-axisymmetric magnetic component is added. The picture is different in the case of 3D saturated MHD equilibrium like the helical core associated with a non-resonant internal kink mode. In assuming nested flux-surfaces, such bifurcated states, expected for a reversed q-profile with qmin close yet above unity [4] and conveniently obtained in VMEC [5], feature integrable field-lines. The helical drift-lines however become resonant with the axisymmetric component in the region of qmin and spontaneously generate drift-islands. Due to the locally reversed sheared q-profile, the drift-island structure follows the bifurcation/reconnection mechanism found in non-twist maps [6, 7]. This result provides a theoretical interpretation of NBI fast ion helical hot-spots in Long-Lived Modes [8] as well as snake-like impurity density accumulation in internal MHD activity [9].
    [1] John R. Cary & Robert G. Littlejohn, Ann. Physics 151, 1 (1983)
    [2] Machiel de Rover, Niek J. Lopes Cardozo & Attila Montvai, Phys. Plasmas 3, 4478 (1996)
    [3] M. Brambilla & A.J. Lichtenberg, Nucl. Fusion 13, 517 (1973)
    [4] F.L. Waelbroeck, Phys. Fluids B 1, 499 (1989)
    [5] W.A. Cooper, J.P. Graves et al., Phys. Rev. Lett. 105, 035003 (2010)
    [6] P.J. Morrison, Phys. Plasmas 7, 2279 (2000)
    [7] R. Balescu, Phys. Rev. E 58, 3781 (1998)
    [8] D. Pfefferlé, J.P. Graves et al., Nucl. Fusion 54, 064020 (2014)
    [9] L. Delgado-Aparicio, L. Sugiyama et al., Phys. Rev. Lett. 110, 065006 (2013)
  • Nonlinear quantum electrodynamics in strong laser fields: From basic concepts to electron-positron photoproduction
    Sebastian Meuren, Max-Planck-Institut für Kernphysik, Germany, abstract, slides
    [#s20, 14 Jan 2016]

    Since the work of Sauter in 1931 it is known that quantum electrodynamics (QED) exhibits a so-called "critical" electromagnetic field scale, at which the quantum interaction between photons and macroscopic electromagnetic fields becomes nonlinear. One prominent example is the importance of light-light interactions in vacuum at this scale, which violates the superposition principle of classical electrodynamics. Furthermore, an electromagnetic field becomes unstable in this regime, as electron-positron pairs can be spontaneously created from the vacuum at the expenses of electromagnetic-field energy (Schwinger mechanism). Unfortunately, the QED critical field scale is so high that experimental investigations are challenging. One promising pathway to explore QED in the nonlinear domain with existing technology consists in the combination of modern (multi) petawatt optical laser systems with highly energetic particles. The suitability of this approach was first demonstrated in the mid-1990s at the seminal SLAC E-144 experiment. Since then, laser technology continuously developed, implying the dawn of a new era of strong-field QED experiments. For instance, the basic processes nonlinear Compton scattering and Breit-Wheeler pair production are expected to influence laser-matter interactions and in particular plasma physics at soon available laser intensities. Therefore, a considerable effort is being undertaken to include these processes into particle-in-cell (PIC) codes used for numerical simulations.

    In the first part of the talk the most prominent nonlinear QED phenomena are presented and discussed on a qualitative level. Afterwards, the mathematical formalism needed for calculations with strong plane-wave background fields is introduced with an emphasize on fundamental concepts. Finally, the nonlinear Breit-Wheeler process is considered more in depth. In particular, the semiclassical approximation is elaborated, which serves as a starting point for the implementation of quantum processes into PIC codes.

  • Attempting a Theoretical Framework for High Energy Density Matter
    Swadesh Mahajan, U. Texas, Austin, abstract, slides
    [#s21, 08 Jan 2016]
    Two distinct approaches to construct a theoretical framework for fully relativistic high energy density (HED) systems- in particular, an assembly of particles in the field of an electromagnetic (EM) wave of arbitrary magnitude- are explored:
    1) The Statistical-Mechanical model for the HED matter is built through the following steps: First, the eigenvalue problem for a quantum relativistic particle immersed in an arbitrary strength EM field is solved; the resulting energy eigenvalue (dependent on $A$, $\omega$ and $k$) is invoked to define the appropriate Boltzmann factor for constructing expressions for physical variables for a weakly interacting system of these field-dressed particles. The fluid equations are the conservation laws,
    2) In the second approach, an equivalent nonlinear quantum mechanics is constructed to represent a hot fluid with and without internal degrees of freedom (like spin). Representative initial results are displayed and discussed. Some notable results are: 1) fundamental changes in the particle energy momentum relationship, 2) The EM wave induces anisotropy in the energy momentum tensor, 3) the EM wave splits the spin-degenerate states, 4) the propagation characteristics of the EM wave are modified by thermal and field effects causing differential self-induced transparency, 5) Particle trapping and ``pushing'' by the high amplitude EM wave. Attempts will be made to highlight testable predictions.
  • The electrostatic response to edge islands induced by Resonant Magnetic Perturbations
    Gianluca Spizzo, Consorzio RFX, Italy, abstract, slides
    [#s22, 03 Dec 2015]
    Measurements of plasma potential have been experimentally determined in great detail in the edge of the RFX reversed-field pinch (RFP) [1], and of the TEXTOR tokamak, with applied magnetic perturbations (MP's) [2]. Generally speaking, the potential has the form $\Phi(r,t,z) = \Phi_0 \sin u$, with $u$ the helical angle: this fact implies a strong correlation between the magnetic field topology and the poloidal/toroidal modulation of the measured plasma potential. In a chaotic tokamak edge, the ion and electron drifts yield a predominantly electron driven radial diffusion when approaching the island X-point, while ion diffusivities are generally an order of magnitude smaller. In the RFP the picture is more complicated, since X-points can act both as drivers of electron diffusion, or dynamical traps (reduced electron diffusion), depending on the helicity of the dominant island [3]. In both devices, this differential electron-to-ion diffusion, causes a strong radial electric field structure pointing outward (inward) from the island O-point. An analytical model for the plasma potential is implemented in the code Orbit [4], and analyses of the ambipolar flow shows that both ion- and electron-dominated transport regimes can exist, which are known as ion and electron roots in stellarators. Moreover, the good agreement found between measured and modeled plasma potential supports that a magnetic island in the plasma edge can act as convective cell. These findings and comparison with stellarator literature, suggests that the role of magnetic islands as convective cells and hence as major radial particle transport drivers could be a generic mechanism in 3D plasma boundary layers.
    [1] P. Scarin, N. Vianello et al., Nucl. Fusion 51, 073002 (2011)
    [1] N. Vianello, C. Rea et al., Plasma Phys. Control. F. 57, 014027 (2015)
    [2] G. Ciaccio, O. Schmitz et al., Phys. Plasmas 22, 102516 (2015)
    [3] G. Spizzo, N. Vianello et al., Phys. Plasmas 21, 056102 (2014)
    [4] R.B. White & M.S. Chance, Phys. Fluids 27, 2455 (1984)