Turbulence-driven ion beams in space plasmas

The description of the local turbulent energy transfer and the high-resolution ion distributions measured by the Magnetospheric Multiscale mission together provide a formidable tool to explore the cross-scale connection between the fluid-scale energy cascade and plasma processes at subion scales. When the small-scale energy transfer is dominated by Alfv´enic, correlated velocity, and magnetic field fluctuations, beams of accelerated particles are more likely observed. Both space observations and numerical simulations suggest the nonlinear wave-particle interaction as one possible mechanism for the energy dissipation in space plasmas

Afterlive: A performant code for Vlasov-Hybrid simulations

A parallelized implementation of the Vlasov-Hybrid method [Nunn, 1993] is presented. This method is a hybrid between a gridded Eulerian description and Lagrangian meta-particles. Unlike the Particle-in-Cell method [Dawson, 1983] which simply adds up the contribution of meta-particles, this method does a reconstruction of the distribution function $f$ in every time step for each species. This interpolation method combines meta-particles with different weights in such a way that particles with large weight do not drown out particles that represent small contributions to the phase space density. These core properties allow the use of a much larger range of macro factors and can thus represent a much larger dynamic range in phase space density. The reconstructed phase space density $f$ is used to calculate momenta of the distribution function such as the charge density $\rho$. The charge density $\rho$ is also used as input into a spectral solver that calculates the self-consistent electrostatic field which is used to update the particles for the next time-step. Afterlive (A Fourier-based Tool in the Electrostatic limit for the Rapid Low-noise Integration of the Vlasov Equation) is fully parallelized using MPI and writes output using parallel HDF5. The input to the simulation is read from a JSON description that sets the initial particle distributions as well as domain size and discretization constraints. The implementation presented here is intentionally limited to one spatial dimension and resolves one or three dimensions in velocity space. Additional spatial dimensions can be added in a straight forward way, but make runs computationally even more costly.Comment: Accepted for publication in Computer Physics Communication

Moment-Based Accelerators for Kinetic Problems with Application to Inertial Confinement Fusion

In inertial confinement fusion (ICF), the kinetic ion and charge separation field effects may play a significant role in the difference between the measured neutron yield in experiments and the predicted yield from fluid codes. Two distinct of approaches exists in modeling plasma physics phenomena: fluid and kinetic approaches. While the fluid approach is computationally less expensive, robust closures are difficult to obtain for a wide separation in temperature and density. While the kinetic approach is a closed system, it resolves the full 6D phase space and classic explicit numerical schemes restrict both the spatial and time-step size to a point where the method becomes intractable. Classic implicit system require the storage and inversion of a very large linear system which also becomes intractable. This dissertation will develop a new implicit method based on an emerging moment-based accelerator which allows one to step over stiff kinetic time-scales. The new method converges the solution per time-step stably and efficiently compared to a standard Picard iteration. This new algorithm will be used to investigate mixing in Omega ICF fuel-pusher interface at early time of the implosion process, fully kinetically

Contemporary particle-in-cell approach to laser-plasma modelling

Particle-in-cell (PIC) methods have a long history in the study of laser-plasma interactions. Early electromagnetic codes used the Yee staggered grid for field variables combined with a leapfrog EM-field update and the Boris algorithm for particle pushing. The general properties of such schemes are well documented. Modern PIC codes tend to add to these high-order shape functions for particles, Poisson preserving field updates, collisions, ionisation, a hybrid scheme for solid density and high-field QED effects. In addition to these physics packages, the increase in computing power now allows simulations with real mass ratios, full 3D dynamics and multi-speckle interaction. This paper presents a review of the core algorithms used in current laser-plasma specific PIC codes. Also reported are estimates of self-heating rates, convergence of collisional routines and test of ionisation models which are not readily available elsewhere. Having reviewed the status of PIC algorithms we present a summary of recent applications of such codes in laser-plasma physics, concentrating on SRS, short-pulse laser-solid interactions, fast-electron transport, and QED effects

Kinetic Solvers with Adaptive Mesh in Phase Space

An Adaptive Mesh in Phase Space (AMPS) methodology has been developed for solving multi-dimensional kinetic equations by the discrete velocity method. A Cartesian mesh for both configuration (r) and velocity (v) spaces is produced using a tree of trees data structure. The mesh in r-space is automatically generated around embedded boundaries and dynamically adapted to local solution properties. The mesh in v-space is created on-the-fly for each cell in r-space. Mappings between neighboring v-space trees implemented for the advection operator in configuration space. We have developed new algorithms for solving the full Boltzmann and linear Boltzmann equations with AMPS. Several recent innovations were used to calculate the discrete Boltzmann collision integral with dynamically adaptive mesh in velocity space: importance sampling, multi-point projection method, and the variance reduction method. We have developed an efficient algorithm for calculating the linear Boltzmann collision integral for elastic and inelastic collisions in a Lorentz gas. New AMPS technique has been demonstrated for simulations of hypersonic rarefied gas flows, ion and electron kinetics in weakly ionized plasma, radiation and light particle transport through thin films, and electron streaming in semiconductors. We have shown that AMPS allows minimizing the number of cells in phase space to reduce computational cost and memory usage for solving challenging kinetic problems