Lee A. Berry

Power deposition in high‐density inductively coupled plasma tools for semiconductor processing

A two‐dimensional, computationally efficient numerical model is developed to study power deposition in high‐density inductively coupled plasma sources. Calculations include both inductive coupling, caused by plasma response to external coil currents, and capacitive coupling, caused by plasma response to external voltages on the coils and wafer. Variation in current along the induction coil is determined self‐consistently from the integral constraint of charge conservation. Sheath phenomena are incorporated through previously published analytic models. The system behavior is analogous in some respects to that of a transmission line. Comparison with measurement suggests that this model provides a good description of self‐consistent coil response when the electric field exhibits less than a quarter wavelength per coil turn.

All-orders spectral calculation of radio-frequency heating in two-dimensional toroidal plasmas

Spectral calculations of radio-frequency (rf) heating in tokamak plasmas are extended to two dimensions (2-D) by taking advantage of new computational tools for distributed memory, parallel computers. The integral form of the wave equation is solved in 2-D without any assumption regarding the smallness of the ion Larmor radius (ρ) relative to the perpendicular wavelength (λ⊥). Results are therefore applicable to all orders in k⊥ρ, where k⊥=2π/λ⊥. Previous calculations of rf wave propagation and heating in 2-D magnetized plasmas have relied on finite Larmor radius expansions (k⊥ρ≪1) and are thus limited to relatively long wavelengths. In this paper, no such assumption is made, and we consider short wavelength processes such as the excitation and absorption of ion Bernstein waves in 2-D with k⊥ρ>1. Results show that this phenomenon is far more complex than simple one-dimensional plasma models would suggest. Other applications include fully self-consistent 2-D solutions for high-harmonic fast-wave heating in...

Physics issues of compact drift optimized stellarators

Physics issues are discussed for compact stellarator configurations which achieve good confinement by the fact that the magnetic field modulus |B| in magnetic co-ordinates is dominated by poloidally symmetric components. Two distinct configuration types are considered: (1) those which achieve their drift optimization and rotational transform at low β and low bootstrap current by appropriate plasma shaping; and (2) those which have a greater reliance on plasma β and bootstrap currents for supplying the transform and obtaining quasi-poloidal symmetry. Stability analysis of the latter group of devices against ballooning, kink and vertical displacement modes has indicated that stable β values on the order of 15% are possible. The first class of devices is being considered for a low β near term experiment that could explore some of the confinement features of the high β configurations.

Physics issues in the design of high-beta, low-aspect-ratio stellarator experiments

High-beta, low-aspect-ratio ~‘‘compact’’ ! stellarators are promising solutions to the problem of developing a magnetic plasma configuration for magnetic fusion power plants that can be sustained in steady state without disrupting. These concepts combine features of stellarators and advanced tokamaks and have aspect ratios similar to those of tokamaks ~2‐4!. They are based on computed plasma configurations that are shaped in three dimensions to provide desired stability and transport properties. Experiments are planned as part of a program to develop this concept. A b54% quasi-axisymmetric plasma configuration has been evaluated for the National Compact Stellarator Experiment ~NCSX!. It has a substantial bootstrap current and is shaped to stabilize ballooning, external kink, vertical, and neoclassical tearing modes without feedback or close-fitting conductors. Quasi-omnigeneous plasma configurations stable to ballooning modes at b54% have been evaluated for the Quasi-Omnigeneous Stellarator ~QOS! experiment. These equilibria have relatively low bootstrap currents and are insensitive to changes in beta. Coil configurations have been calculated that reconstruct these plasma configurations, preserving their important physics properties. Theory- and experiment-based confinement analyses are used to evaluate the technical capabilities needed to reach target plasma conditions. The physics basis for these complementary experiments is described.

Advances in full-wave modeling of radio frequency heated, multidimensional plasmas

Previous full-wave models for rf heating in multidimensional plasmas have relied on either cold-plasma or finite Larmor radius approximations. These models assume that the perpendicular wavelength of the rf field is much larger than the ion Larmor radius, and they are therefore limited to relatively long wavelengths and low cyclotron harmonics. Recently, alternate full-wave models have been developed that eliminate these restrictions. These “all orders spectral algorithms” take advantage of new computational techniques for massively parallel computers to solve the integral form of the wave equation in multiple dimensions without any restriction on wavelength relative to orbit size, and with no limit on the number of cyclotron harmonics retained. These new models give high-resolution, two-dimensional solutions for mode conversion and high harmonic fast wave heating in tokamak geometry. In addition, they have been extended to give fully three-dimensional solutions of the integral wave equation for minority ...

Self-consistent full-wave and Fokker-Planck calculations for ion cyclotron heating in non-Maxwellian plasmas

Magnetically confined plasmas can contain significant concentrations of nonthermal plasma particles arising from fusion reactions, neutral beam injection, and wave-driven diffusion in velocity space. Initial studies in one-dimensional and experimental results show that nonthermal energetic ions can significantly affect wave propagation and heating in the ion cyclotron range of frequencies. In addition, these ions can absorb power at high harmonics of the cyclotron frequency where conventional two-dimensional global-wave models are not valid. In this work, the all-orders global-wave solver AORSA [E. F. Jaeger et al., Phys. Rev. Lett. 90, 195001 (2003)] is generalized to treat non-Maxwellian velocity distributions. Quasilinear diffusion coefficients are derived directly from the wave fields and used to calculate energetic ion velocity distributions with the CQL3D Fokker-Planck code [R. W. Harvey and M. G. McCoy, Proceedings of the IAEA Technical Committee Meeting on Simulation and Modeling of Thermonuclear ...

Full wave simulations of fast wave heating losses in the scrape-off layer of NSTX and NSTX-U

Full wave simulations of fusion plasmas show a direct correlation between the location of the fast-wave cut-off, radiofrequency (RF) field amplitude in the scrape-off layer (SOL) and the RF power losses in the SOL observed in the National Spherical Torus eXperiment (NSTX). In particular, the RF power losses in the SOL increase significantly when the launched waves transition from evanescent to propagating in that region. Subsequently, a large amplitude electric field occurs in the SOL, driving RF power losses when a proxy collisional loss term is added. A 3D reconstruction of absorbed power in the SOL is presented showing agreement with the RF experiments in NSTX. Loss predictions for the future experiment NSTX-Upgrade (NSTX-U) are also obtained and discussed.

Full-wave calculation of sheared poloidal flow driven by high-harmonic ion Bernstein waves in tokamak plasmas

A full-wave, one-dimensional spectral model is developed to study sheared poloidal flow driven by high-harmonic ion Bernstein waves (IBWs) in tokamak plasmas. The local plasma conductivity is corrected to lowest order in ρ/L where ρ is the ion Larmor radius and L is the equilibrium scale length. This correction takes into account gradients in equilibrium quantities and is necessary for conservation of energy. It is equivalent to the “odd-order derivative” terms in finite difference models. No assumption is made regarding the smallness of the ion Larmor radius relative to wavelength, and results are applicable to all orders in k⊥ρ where k⊥ is the perpendicular wave number. Previous numerical results for flow drive have relied on expansions in k⊥ρ, and are thus limited to cyclotron harmonics of two and below. In this article, we consider higher-harmonic cases corresponding to recent IBW flow drive experiments on the Tokamak Fusion Test Reactor [B. P. LeBlanc, R. E. Bell, S. Bernabei et al., Phys. Rev. Lett....

Event-based parareal: A data-flow based implementation of parareal

Parareal is an iterative algorithm that, in effect, achieves temporal decomposition for a time-dependent system of differential or partial differential equations. A solution is obtained in a shorter wall-clock time, but at the expense of increased compute cycles. The algorithm combines a fine solver that solves the system to acceptable accuracy with an approximate coarse solver. The critical task for the successful implementation of parareal on any system is the development of a coarse solver that leads to convergence in a small number of iterations compared to the number of time slices in the full time interval, and is, at the same time, much faster than the fine solver. Very fast coarse solvers may not lead to sufficiently rapid convergence, and slow coarse solvers may not lead to significant gains even if the number of iterations to convergence is satisfactory. We find that the difficulty of meeting these conflicting demands can be substantially eased by using a data-driven, event-based implementation of parareal. As a result, tasks for one iteration do not wait for the previous iteration to complete, but are started when the needed data are available. For given convergence properties, the event-based approach relaxes the speed requirements on the coarse solver by a factor of ~K, where K is the number of iterations required for a converged solution. This may, for many problems, lead to an efficient parareal implementation that would otherwise not be possible or would require substantial coarse solver development. In addition, the framework used for this implementation executes a task when the data dependencies are satisfied and computational resources are available. This leads to improved computational efficiency over previous approaches that pipeline or schedule groups of tasks to a particular processor or group of processors.

Mechanisms for the convergence of time-parallelized, parareal turbulent plasma simulations

Parareal is a recent algorithm able to parallelize the time dimension in spite of its sequential nature. It has been applied to several linear and nonlinear problems and, very recently, to a simulation of fully-developed, two-dimensional drift wave turbulence. The mere fact that parareal works in such a turbulent regime is in itself somewhat unexpected, due to the characteristic sensitivity of turbulence to any change in initial conditions. This fundamental property of any turbulent system should render the iterative correction procedure characteristic of the parareal method inoperative, but this seems not to be the case. In addition, the choices that must be made to implement parareal (division of the temporal domain, election of the coarse solver and so on) are currently made using trial-and-error approaches. Here, we identify the mechanisms responsible for the convergence of parareal of these simulations of drift wave turbulence. We also investigate which conditions these mechanisms impose on any successful parareal implementation. The results reported here should be useful to guide future implementations of parareal within the much wider context of fully-developed fluid and plasma turbulent simulations.