skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A numerical study of Landau damping with PETSc-PIC
We present a study of the standard plasma physics test, Landau damping, using the Particle-In-Cell (PIC) algorithm. The Landau damping phenomenon consists of the damping of small oscillations in plasmas without collisions. In the PIC method, a hybrid discretization is constructed with a grid of finitely supported basis functions to represent the electric, magnetic and/or gravitational fields, and a distribution of delta functions to represent the particle field. Approximations to the dispersion relation are found to be inadequate in accurately calculating values for the electric field frequency and damping rate when parameters of the physical system, such as the plasma frequency or thermal velocity, are varied. We present a full derivation and numerical solution for the dispersion relation, and verify the PETSC-PIC numerical solutions to the Vlasov-Poisson for a large range of wave numbers and charge densities.  more » « less
Award ID(s):
1931524
PAR ID:
10493301
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Mathematical Sciences Publishers
Date Published:
Journal Name:
Communications in Applied Mathematics and Computational Science
Volume:
18
Issue:
1
ISSN:
1559-3940
Page Range / eLocation ID:
135 to 152
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. An accurate description of plasma waves is fundamental for the understanding of many plasma phenomena. It is possible to twist plasma waves such that, in addition to having longitudinal motion, they can possess a quantized orbital angular momentum. One such type of plasma wave is the Laguerre–Gaussian mode. Three-dimensional numerical particle-in-cell simulations demonstrate the existence of stable long-lived plasma waves with orbital angular momentum. These waves can be shown to create large amplitude static magnetic fields with unique twisted longitudinal structures. In this paper, we review the recent progress in studies of helical plasma waves and present a new analytical description of a standing Laguerre–Gaussian plasma wave mode along with 3D particle-in-cell simulation results. The Landau damping of twisted plasma waves shows important differences compared to standard longitudinal plasma wave Landau damping. These effects include an increased damping rate, which is affected by both the focal width and the orbital number of the plasma wave. This increase in the damping rate is of the same order as the thermal correction. Moreover, the direction of momentum picked up by resonant particles from the twisted plasma wave can be significantly altered. By contrast, the radial electric field has a subtle effect on the trajectories of resonant electrons. 
    more » « less
  2. In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. The azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA. 
    more » « less
  3. Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (vsh≈ (5500–22000) km s−1) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (vsh≈ 300 km s−1), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations. 
    more » « less
  4. Abstract. It is generally accepted that modeling Farley–Buneman instabilities requires resolving ion Landau damping to reproduce experimentally observed features. Particle-in-cell (PIC) simulations have been able to reproduce most of these but at a computational cost that severely affects their scalability. This limitation hinders the study of non-local phenomena that require three dimensions or coupling with larger-scale processes. We argue that a form of the five-moment fluid system can recreate several qualitative aspects of Farley–Buneman dynamics such as density and phase speed saturation, wave turning, and heating. Unexpectedly, these features are still reproduced even without using artificial viscosity to capture Landau damping. Comparing the proposed fluid models and a PIC implementation shows good qualitative agreement. 
    more » « less
  5. We present the results of numerical studies of the whistler wave parametric decay instability in the system with the suppressed Landau damping of ion acoustic waves (IAWs) based on the self-consistent Darwin particle-in-cell (PIC) model. It has been demonstrated that a monochromatic whistler wave launched along the background magnetic field couples to a counter-propagating whistler mode and co-propagating ion acoustic mode. The coupling of the electromagnetic mode to the electrostatic mode is guided by a ponderomotive force that forms spatio-temporal beat patterns in the longitudinal electric field generated by the counter-propagating whistler and the pump whistler wave. The threshold amplitude for the instability is determined to be δB w / B 0 = 0.028 and agrees with a prediction for the ion decay instability: δB w / B 0 = 0.042 based on the linear kinetic damping rates, and δB w / B 0 = 0.030 based on the simulation derived damping rates. Increasing the amplitude of the pump whistler wave, the secondary and tertiary decay thresholds are reached, and cascading parametric decay from the daughter whistler modes is observed. At the largest amplitude ( δB w / B 0 ∼ 0.1) the primary IAW evolves into a short-lived and highly nonlinear structure. The observed dependence of the IAW growth rate on the pump wave amplitude agrees with the expected trend; however, quantitatively, the growth rate of the IAW is larger than expected from theoretical predictions. We discuss the relevant space regimes where the instability could be observed and extensions to the parametric coupling of whistler waves with the electron acoustic wave (EAW). 
    more » « less