Markus Held

Interchange Instability and Transport in Matter-Antimatter Plasmas

Symmetric electron-positron plasmas in inhomogeneous magnetic fields are intrinsically subject to interchange instability and transport. Scaling relations for the propagation velocity of density perturbations relevant to transport in isothermal magnetically confined electron-positron plasmas are deduced, including damping effects when Debye lengths are large compared to Larmor radii. The relations are verified by nonlinear full-F gyrofluid computations. Results are analyzed with respect to planned magnetically confined electron-positron plasma experiments. The model is generalized to other matter-antimatter plasmas. Magnetized electron-positron-proton-antiproton plasmas are susceptible to interchange-driven local matter-antimatter separation, which can impede sustained laboratory magnetic confinement.

Streamline integration as a method for two-dimensional elliptic grid generation

Abstract We propose a new numerical algorithm to construct a structured numerical elliptic grid of a doubly connected domain. Our method is applicable to domains with boundaries defined by two contour lines of a two-dimensional function. Furthermore, we can adapt any analytically given boundary aligned structured grid, which specifically includes polar and Cartesian grids. The resulting coordinate lines are orthogonal to the boundary. Grid points as well as the elements of the Jacobian matrix can be computed efficiently and up to machine precision. In the simplest case we construct conformal grids, yet with the help of weight functions and monitor metrics we can control the distribution of cells across the domain. Our algorithm is parallelizable and easy to implement with elementary numerical methods. We assess the quality of grids by considering both the distribution of cell sizes and the accuracy of the solution to elliptic problems. Among the tested grids these key properties are best fulfilled by the grid constructed with the monitor metric approach.

Non-Oberbeck–Boussinesq zonal flow generation

Novel mechanisms for zonal flow (ZF) generation for both large relative density fluctuations and background density gradients are presented. In this non-Oberbeck-Boussinesq (NOB) regime ZFs are driven by the Favre stress, the large fluctuation extension of the Reynolds stress, and by background density gradient and radial particle flux dominated terms. Simulations of a nonlinear full-F gyro-fluid model confirm the predicted mechanism for radial ZF propagation and show the significance of the NOB ZF terms for either large relative density fluctuation levels or steep background density gradients.

Streamline integration as a method for structured grid generation in X-point geometry

Abstract We investigate structured grids aligned to the contours of a two-dimensional flux-function with an X-point (saddle point). Our theoretical analysis finds that orthogonal grids exist if and only if the Laplacian of the flux-function vanishes at the X-point. In general, this condition is sufficient for the existence of a structured aligned grid with an X-point. With the help of streamline integration we then propose a numerical grid construction algorithm. In a suitably chosen monitor metric the Laplacian of the flux-function vanishes at the X-point such that a grid construction is possible. We study the convergence of the solution to elliptic equations on the proposed grid. The diverging volume element and cell sizes at the X-point reduce the convergence rate. As a consequence, the proposed grid should be used with grid refinement around the X-point in practical applications. We show that grid refinement in the cells neighbouring the X-point restores the expected convergence rate.

Unified transport scaling laws for plasma blobs and depletions

We study the dynamics of seeded plasma blobs and depletions in an (effective) gravitational field. For incompressible flows, the radial center of mass velocity of blobs and depletions is proportional to the square root of their initial cross-field size and amplitude. If the flows are compressible, this scaling holds only for ratios of amplitude to size larger than a critical value. Otherwise, the maximum blob and depletion velocity depends linearly on the initial amplitude and is independent of size. In both cases, the acceleration of blobs and depletions depends on their initial amplitude relative to the background plasma density and is proportional to gravity and independent of their cross-field size. Due to their reduced inertia plasma, depletions accelerate more quickly than the corresponding blobs. These scaling laws are derived from the invariants of the governing drift-fluid equations for blobs and agree excellently with numerical simulations over five orders of magnitude for both blobs and depleti...

Lattice Boltzmann model for collisionless electrostatic drift wave turbulence obeying Charney-Hasegawa-Mima dynamics

A lattice Boltzmann method (LBM) approach to the Charney-Hasegawa-Mima (CHM) model for adiabatic drift wave turbulence in magnetised plasmas is implemented. The CHM-LBM model contains a barotropic equation of state for the potential, a force term including a cross-product analogous to the Coriolis force in quasigeostrophic models, and a density gradient source term. Expansion of the resulting lattice Boltzmann model equations leads to cold-ion fluid continuity and momentum equations, which resemble CHM dynamics under drift ordering. The resulting numerical solutions of standard test cases (monopole propagation, stable drift modes and decaying turbulence) are compared to results obtained by a conventional finite difference scheme that directly discretizes the CHM equation. The LB scheme resembles characteristic CHM dynamics apart from an additional shear in the density gradient direction. The occurring shear reduces with the drift ratio and is ascribed to the compressible limit of the underlying LBM. A lattice Boltzmann model for drift wave turbulence in magnetised plasmas is proposed.Typical test cases are compared to simulations of the Charney-Hasegawa-Mima equation.The new model adds an additional shear to the drift wave turbulence behaviour.