Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear plasma device”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ion cyclotron heating at high plasma density in Proto-MPEX

The physics of ion cyclotron heating (ICH) relevant to the steady-state linear machine MPEX (Material Plasma Exposure eXperiment) has been explored in its predecessor, short-pulse device: Proto-MPEX. MPEX will utilize fundamental ICH to increase heat flux at the target and produce ion temperatures and velocity distributions with improved fidelity to those found in a tokamak divertor region, in comparison to those produced by substrate biasing. In the experiments on Proto-MPEX described here, bulk ion temperatures up to ~15 eV have been achieved with 20 kW net ICH power at 6.5 MHz, using ICH heating of a deuterium plasma produced by a helicon plasma source. The heat flux at the target has been observed to increase throughout the plasma cross section, including in the core region. Core $T_i$ and target heat flux are observed to scale linearly with injected ICH power. Further, measurements of plasma loading and target heat flux as a function of the magnetic field strength at the antenna, together with modeling of the wave propagation from the antenna to the ion cyclotron resonance using the ANTENA and COMSOL codes with a warm plasma dielectric tensor, indicate that power is coupled to the core plasma via fast wave excitation of a kinetic Alfvén wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Data for Construction of a Compact Array of Microplasma Jet Devices and Its Application for Random Mutagenesis of Rhodosporidium toruloides

A small and efficient DNA mutation-inducing machine was constructed with an array of microplasma jet devices (7 × 1) that can be operated at atmospheric pressure for microbial mutagenesis. Using this machine, we report disruption of a plasmid DNA and generation of mutants of an oleaginous yeast Rhodosporidium toruloides . Specifically, a compact-sized microplasma channel (25 × 20 × 2 mm3) capable of generating an electron density of greater than 1013 cm–3 was constructed to produce reactive species (N2*, N2+, O, OH, and Hα) under helium atmospheric conditions to induce DNA mutagenesis. The length of microplasma channels in the device played a critical role in augmenting both the volume of plasma and the concentration of reactive species. First, we confirmed that microplasma treatment can linearize a plasmid by creating nicks in vitro. Second, we treated R. toruloides cells with a jet device containing 7 microchannels for 5 min; 94.8% of the treated cells were killed, and 0.44% of surviving cells showed different colony colors as compared to their parental colony. Microplasma-based DNA mutation is energy-efficient and can be a safe alternative for inducing mutations compared to conventional methods using toxic mutagens. This compact and scalable device is amenable for industrial strain improvement involving large-scale mutagenesis.

Conversion↗
DOE OSTI2021

Integrated model predictions on the impact of substrate damage on gas dynamics during ITER burning-plasma operations

Divertor design and choice of plasma-facing materials (PFM) will be essential to the success of next-generation fusion reactors as they operate under more powerful scenarios. Understanding and controlling interactions between the plasma and PFM is essential to making these choices. Within these plasma–material interactions and especially in tungsten (W), the interplay between the most abundant plasma species (hydrogen isotopes and helium, He) with the wall material alters fuel retention. However, this interplay is yet to be sufficiently understood to confidently project fuel retention levels to future fusion devices. The paper presents a series of integrated simulations of fusion plasmas and their interaction with tungsten. Specifically, this study assesses the impact of He plasma pre-exposure on hydrogenic species retention during 100 s of burning plasma operations (BPO) in ITER. Multiple pre-exposure scenarios are considered, including sub-surface damage resulting from exposures in the linear device PISCES and from early ITER He-operation. The predictions from these consecutive He-BPO exposures show that fuel content and spatial distribution in the material are largely determined by the He-induced damage, as manifest in: (i) changes in surface temperature expected during BPO have little effect on fuel retention in the presence of He-induced damage; (ii) gas content stabilizes quickly in substrates pre-exposed in PISCES, at levels set by the concentration of pre-existing vacancies, while it continues to increase in substrates initially pristine or pre-exposed to ITER He plasmas; (iii) the presence of He and He–V clusters in the near-surface region locally increases hydrogenic retention, but decreases its permeation; this results in hydrogenic species that remain closer to the surface in pre-damaged substrates, while the bulk content is higher for initially pristine cases. In summary, the interaction and binding of D and T with the pre-existing He–V clusters modifies retention and permeation of hydrogen species during ITER BPO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

New linear stability parameter to describe low-β electromagnetic microinstabilities driven by passing electrons in axisymmetric toroidal geometry

In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electromagnetic turbulence in increased β plasmas in the Large Plasma Device

The variation of pressure-gradient-driven turbulence with plasma $\beta$ (up to $\beta \approx 15\,\%$ ) is investigated in linear, magnetized plasma. The magnitude of magnetic fluctuations is observed to increase substantially with increasing $\beta$ . More importantly, parallel magnetic fluctuations are observed to dominate at higher $\beta$ values, with $\delta B_\parallel / \delta B_\perp \approx 2$ and $\delta B / B_0 \approx 1\,\%$ . Parallel magnetic fluctuations are strongly correlated with density fluctuations and the two are observed to be out of phase. The relative magnitude of and cross-phase between density and parallel magnetic field fluctuations are consistent with the dynamic pressure balance ( $P+{B_{0}^2}/{2\mu _0} = \textrm {constant}$ ). A local slab model theory for electromagnetic, modified drift Alfvén waves, including parallel magnetic fluctuations, shows partial agreement with experimental observations.

Physics↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2021

Predicting the critical gradient of ITG turbulence in fusion plasmas

The quasilinear mixing-length approach to efficient prediction of transport in fusion devices is improved to account for the ‘Dimits’ upshift between linear and nonlinear critical pressure gradients in zonal-flow-saturated turbulence regimes. This modification uses the frequency mismatch between modes interacting turbulently to track changes in saturation efficiency. Near criticality, energy is transferred exclusively to stable eigenmodes, rapidly increasing the efficacy of the nonlinearity. Finally, the modified quasilinear model is able to predict below-threshold turbulent ion-temperature-gradient-driven transport accurately and also yields significantly improved predictions for trapped-electron-mode transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FAR3d

The FAR3d model calculates the linear and nonlinear stability properties of energetic particle driven Alfven instabilities for both tokamak and stellarator plasma confinement devices using gyro-landau closure methods. This is an important fundamental physics problem for existing fusion energy experiments and for future fusion reactors.

Varela, Jacobo↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Disentanglement of density and rotation dependences of the field penetration threshold on the J-TEXT tokamak

The field penetration threshold of magnetic perturbations has been observed to vary non-monotonically with an increase of density in ohmic plasmas on the J-TEXT tokamak. This observation appears contradicting the previous empirical density scaling law. Disentanglement of plasma density and rotation dependences of the field penetration threshold has been carried out. It shows that the field penetration threshold depends only weakly on the density but linearly on the plasma rotation. This result is not only important for the prediction of error field tolerance in fusion devices, but also opens a question on the role of density in the forced magnetic reconnection process in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Laboratory study of the PFRC-2's initial plasma densification stages

Initial plasma densification by odd-parity rotating magnetic fields (RMF o ) applied to the linear magnetized Princeton field-reversed configuration (PFRC-2) device with fill gases at pressures near 1 mTorr proceeds through two phases: a slow one, characterized by a rise time $τ_{slow}$ ~ 100 $μ$s, followed by a fast one, characterized by $τ_{fast}$ ~ 10 $μ$s. The transition from slow to fast occurs at a line-integral-averaged electron density, t n e , near 2$\times$ 10 11 cm –3 , independent of magnetic field. Here, over most of the range of experimental parameters investigated, as the PFRC-2 axial magnetic field strength was increased, RMF o power decreased, gas fill pressure lowered, or lower atomic mass unit (AMU) fill gas used, the duration of the slow phase lengthened from 50 $μ$s to longer than 10 ms after the RMF o power began. The post-fast-phase maximum n e increases with the fill-gas AMU, exceeding 5 × 10 13 cm –3 for Ar. The slow phase is consistent with atomic physics processes and field-parallel sound-speed losses. The fast phase may be explained by improved axial confinement, possibly augmented by radial or axial contraction of the plasma. Another possible explanation, a large increase in electron temperature, is inconsistent with x-ray emission. The n e behavior is discussed in relation to the E to H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Using the Stix finite element RF code to investigate operation optimization of the ICRF antenna on Alcator C-Mod

Abstract As the Ion Cyclotron Radio Frequency range (ICRF) heating becomes more favorable in fusion devices, the urgency of predicting and mitigating impurity generation that arises from it becomes more pressing. In the ICRF regime, rectified Radio Frequency (RF) sheaths are known to form at antenna and material edges that influence negative effects like sputtering and a decrease in heating efficiency. Methods to mitigate the formation of these RF sheaths through RF image currents cancellation have been experimentally studied. A power-phasing scan done on Alcator C-Mod in which the amount of power on the two inner straps ( P in ) versus the total 4 straps ( P tot ) was varied showed a minimization of enhanced potentials between P in / P tot ∼ 0.7–0.9 while impurities were minimized for P in / P tot ∼ 0.5–0.8. New capabilities in the realm of representing the RF sheath numerically now allow for these experiments to be simulated. Given the size of the sheath relative to the scale of the device, it can be approximated as a Boundary Condition (BC). A new parallelized cold-plasma wave equation solver called Stix implements a non-linear sheath impedance model BC formulated by Myra et al (2015 Phys. Plasmas 22 062507) through the method of finite elements using the MFEM library [ http://mfem.org ]. It is seen that Stix shows qualitative agreement with the measured C-Mod enhanced potentials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗