skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A method for obtaining three-dimensional computation equilibrium of non-neutral plasmas using WARP

Abstract

Computer simulation studies of the stability and transport properties of trapped non-neutral plasmas require the numerical realization of a three-dimensional plasma distribution. This paper presents a new numerical method for obtaining, without an explicit model for physical collisions in the code, a low noise three-dimensional computational equilibrium distribution. This requires both the loading of particles into an idealized distribution and the relaxation from that distribution toward an approximate numerical equilibrium. The equilibrium can then be modified through a slow change of system parameters, to generate other equilibria. In the present work we apply this method to a UC Berkeley experiment on electron confinement in magnetic geometries appropriate for the ALPHA anti-hydrogen experiment, using the three-dimensional Particle-In-Cell code WARP. WARP's guiding center mover and its option to switch between different solvers during a simulation are highly valuable because they speed up the simulations; they enable the practical use of the new technique for generating numerical equilibrium states of trapped nonneutral plasmas.

Authors:
; ; ; ; ;
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
Accelerator& Fusion Research Division
OSTI Identifier:
928303
Report Number(s):
LBNL-256E
Journal ID: ISSN 0021-9991; JCTPAH; TRN: US0804322
DOE Contract Number:
DE-AC02-05CH11231
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics
Country of Publication:
United States
Language:
English
Subject:
70; COMPUTERIZED SIMULATION; CONFINEMENT; DISTRIBUTION; ELECTRONS; PLASMA; RELAXATION; SIMULATION; STABILITY; TRANSPORT; VELOCITY; WARP Non-Neutral Plasmas

Citation Formats

Wurtele, J., Wurtele, J., Friedman, A., Grote, D.P., Vay, J-L., and Gomberoff, K. A method for obtaining three-dimensional computation equilibrium of non-neutral plasmas using WARP. United States: N. p., 2006. Web.
Wurtele, J., Wurtele, J., Friedman, A., Grote, D.P., Vay, J-L., & Gomberoff, K. A method for obtaining three-dimensional computation equilibrium of non-neutral plasmas using WARP. United States.
Wurtele, J., Wurtele, J., Friedman, A., Grote, D.P., Vay, J-L., and Gomberoff, K. Sat . "A method for obtaining three-dimensional computation equilibrium of non-neutral plasmas using WARP". United States. doi:. https://www.osti.gov/servlets/purl/928303.
@article{osti_928303,
title = {A method for obtaining three-dimensional computation equilibrium of non-neutral plasmas using WARP},
author = {Wurtele, J. and Wurtele, J. and Friedman, A. and Grote, D.P. and Vay, J-L. and Gomberoff, K.},
abstractNote = {Computer simulation studies of the stability and transport properties of trapped non-neutral plasmas require the numerical realization of a three-dimensional plasma distribution. This paper presents a new numerical method for obtaining, without an explicit model for physical collisions in the code, a low noise three-dimensional computational equilibrium distribution. This requires both the loading of particles into an idealized distribution and the relaxation from that distribution toward an approximate numerical equilibrium. The equilibrium can then be modified through a slow change of system parameters, to generate other equilibria. In the present work we apply this method to a UC Berkeley experiment on electron confinement in magnetic geometries appropriate for the ALPHA anti-hydrogen experiment, using the three-dimensional Particle-In-Cell code WARP. WARP's guiding center mover and its option to switch between different solvers during a simulation are highly valuable because they speed up the simulations; they enable the practical use of the new technique for generating numerical equilibrium states of trapped nonneutral plasmas.},
doi = {},
journal = {Journal of Computational Physics},
number = ,
volume = ,
place = {United States},
year = {Sat Mar 25 00:00:00 EST 2006},
month = {Sat Mar 25 00:00:00 EST 2006}
}
  • Computer simulation studies of the stability and transport properties of trapped non-neutral plasmas require the numerical realization of a three-dimensional plasma distribution. This paper presents a new numerical method for obtaining, without an explicit model for physical collisions in the code, a low noise three-dimensional computational equilibrium distribution. This requires both the loading of particles into an idealized distribution and the relaxation from that distribution toward an approximate numerical equilibrium. The equilibrium can then be modified through a slow change of system parameters, to generate other equilibria. In the present, work we apply this method to a UC Berkeley experimentmore » on electron confinement in magnetic geometries appropriate for the ALPHA anti-hydrogen experiment, using the three-dimensional particle-in-cell code WARP. WARP's guiding center mover and its option to switch between different solvers during a simulation are highly valuable because they speed up the simulations; they enable the practical use of the new technique for generating numerical equilibrium states of trapped non-neutral plasmas.« less
  • It is known that magnetic perturbations can mitigate edge localized modes (ELMs) in experiments, for example, MAST [Kirk et al., Nucl. Fusion 53, 043007 (2013)]. One hypothesis is that the magnetic perturbations cause a three dimensional corrugation of the plasma and this corrugated plasma has different stability properties to peeling-ballooning modes compared to an axisymmetric plasma. It has been shown in an up-down symmetric plasma that magnetic perturbations in tokamaks will break the usual axisymmetry of the plasma causing three dimensional displacements [Chapman et al., Plasma Phys. Controlled Fusion 54, 105013 (2012)]. We produce a free boundary three-dimensional equilibrium ofmore » a lower single null MAST relevant plasma using VMEC [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)]. The safety factor and pressure profiles used for the modelling are similar to those deduced from axisymmetric analysis of experimental data with ELMs. We focus on the effect of applying n = 3 and n = 6 magnetic perturbations using the resonant magnetic perturbation (RMP) coils. A midplane displacement of over ±1 cm is seen when the full current is applied. The current in the coils is scanned and a linear relationship between coil current and midplane displacement is found. The pressure gradient in real space in different toroidal locations is shown to change when RMPs are applied. This effect should be taken into account when diagnosing plasmas with RMPs applied. The helical Pfirsch-Schlüter currents which arise as a result of the assumption of nested flux surfaces are estimated for this equilibrium. The effect of this non-axisymmetric equilibrium on infinite n ballooning stability is investigated using COBRA [Sanchez et al., J. Comput. Phys. 161, 576–588 (2000)]. The infinite n ballooning stability is analysed for two reasons; it may give an indication of the effect of non-axisymmetry on finite n peeling-ballooning modes, responsible for ELMs; and infinite n ballooning modes are correlated to kinetic ballooning modes which are thought to limit the pressure gradient of the pedestal [Snyder et al., Phys. Plasmas 16, 056118 (2009)]. The ballooning mode growth rate gains a variation in toroidal angle. The equilibria with midplane displacements due to RMP coils have a higher ballooning mode growth rate than the axisymmetric case and the possible implications are discussed.« less
  • When a single-species plasma is confined in a harmonic Penning trap at cryogenic temperature, the thermal equilibrium is approximately a uniform density spheroid (ellipsoid of revolution). Normal modes corresponding to quadrupole excitations of this plasma have recently been measured. In this paper, nonlinear equations of motion are derived for these quadrupole oscillations. For large amplitudes, the oscillations deform a spheroidal plasma into a triaxial ellipsoid with time-dependent shape and orientation. The integrals of the motion are found and the cylindrically symmetric finite-amplitude oscillations of a spheroid are studied. An analysis of all possible ellipsoidal equilibria is also carried out. Newmore » equilibria are discovered which correspond to finite-amplitude versions of the noncylindrically symmetric linear quadrupole oscillations. The equilibria are shown to fall into two classes in which the ellipsoids are either tilted or aligned with respect to the magnetic field. Some of these equilibria have densities well above the Brillouin limit.« less
  • The effect of weak external axisymmetric magnetic and electrostatic perturbations on the equilibrium of a non-neutral plasma confined in a Malmberg-Penning trap is analyzed. A semi-analytic solution for the potential variations inside the trap is found in a paraxial limit of the perturbations for the case of global thermal equilibrium. The populations of magnetically and electrostatically trapped particles created by the external perturbations are characterized, and their fractions are calculated explicitly for a bi-Maxwellian distribution function. 2D numerical simulations of the thermal equilibrium of a pure electron plasma in the presence of axial magnetic field perturbations are performed to checkmore » the limits of validity of the analytical 1D approximation.« less