Juan Francisco Caneses

Confinement performance predictions for a high field axisymmetric tandem mirror – CORRIGENDUM - Volume 92 Issue 4

This paper presents a Hammir tandem mirror confinement performance analysis based on Realta Fusion’s first-of-a-kind model for axisymmetric magnetic mirror fusion performance. This model uses an integrated end plug simulation model including, heating, equilibrium and transport combined with a new formulation of the plasma operation contours (POPCONs) technique for the tandem mirror central cell. Using this model in concert with machine learning optimization techniques, it is shown that an end plug utilizing high temperature superconducting magnets and modern neutral beams enables a classical tandem mirror pilot plant producing a fusion gain Q > 5. The approach here represents an important advance in tandem mirror design. The high-fidelity end plug model enables calculations of heating and transport in the highly non-Maxwellian end plug to be made more accurately. The detailed end plug modelling performed in this work has highlighted the importance of classical radial transport and neutral beam absorption efficiency on end plug viability. The central cell POPCON technique allows consideration of a wide range of parameters in the relatively simple near-Maxwellian central cell, facilitating the selection of more optimal central cell plasmas. These advances make it possible to find more conservative classical tandem mirror fusion pilot plant operating points with lower temperatures, neutral beam energies and end plug performance requirements than designs in the literature. Despite being more conservative, it is shown that these operating points have sufficient confinement performance to serve as the basis of a viable fusion pilot plant provided that they can be stabilized against magnetohydrodynamic and trapped particle modes.

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

The steady-state linear device “Material Plasma Exposure eXperiment” (MPEX) is currently under construction at Oak Ridge National Laboratory with the goal of enabling Plasma-Material Interaction studies at future fusion reactor relevant plasma conditions. In this work, a newly in-house developed hybrid Particle-In-Cell code-PICOS++ is applied to understand the experimental results obtained from the prototype of MPEX referred to as the “Proto-MPEX” during its helicon-only and helicon with ion cyclotron resonance heating (ICRH) experiments. This study explains the physics of the experimentally observed plasma density-drop at the divertor target in Proto-MPEX device during ICRH. In contrast to previous work on ICRH in MPEX [Kumar et al. Nucl. Fusion, 63, 036004 (2023)], this study demonstrates that the mirror force plays a central role in the Proto-MPEX plasma transport during ICRH, which has new features not previously explored. Force balance analyses reveal that the temperature anisotropy produced by ICRH leads to a significant increase in the mirror force downstream of the resonance where the magnetic field is diverging. This force accelerates ions toward the target and leads to a drop in plasma density to ensure conservation of particle flux. Simulations with ICRH where the magnetic field divergence downstream of the resonance has been removed, do not produce plasma acceleration nor density drop at the target despite efficient ion heating at the resonance. Moreover, simulation results demonstrate that for a given ICRH power, lowering the source rate produces ions with increased perpendicular energy which interact with the mirror force to produce higher plasma acceleration which increases the strength of the density-drop at the target. The strength of the density drop appears to reach an asymptotic limit at a certain threshold ICRH power. Simulations show that this threshold power increases with increasing particle source rate.

The Materials Plasma Exposure eXperiment (MPEX) is a steady-state linear plasma device that will address plasma-material interaction (PMI) science and enable testing of fusion reactor-relevant divertor plasma-facing materials. The MPEX source concept consists of a helicon plasma source to generate the plasma, electron cyclotron heating (ECH) for electron heating, and ion cyclotron heating (ICH) for ion heating. The MPEX source plasma is then transported axially to the PMI material target region to test material samples in fusion reactor-relevant divertor conditions. This paper will summarize the physical and technical basis of MPEX. The paper will first define the MPEX parameters and scenarios at the target relevant to PMI science for various fusion reactor-relevant divertor conditions and show plasma transport modeling results to set the MPEX source parameters. Recent experimental and modeling results from Proto-MPEX, a short-pulse experiment to develop the plasma production, heating, and transport physics for MPEX, will be shown. From these results, it will be shown that MPEX can reach its desired scenarios. The MPEX physical and technical basis will also determine important functional requirements for magnetic field, radiofrequency (RF) power, RF frequency, and neutral pressure in the helicon, ECH, ICH, and PMI regions that are required to achieve the desired MPEX scenarios. The necessity for key in-vessel components such as skimmers, limiters, and microwave absorbers will also be highlighted.

The Material Plasma Exposure eXperiment (MPEX) is a steady state linear device with the goal to perform plasma material interaction studies at future fusion reactor relevant conditions. A prototype of MPEX referred as ‘Proto-MPEX’ is designed to carry out research and development related to source, heating and transport concepts on the planned full MPEX device. The auxiliary heating schemes in MPEX are based on cyclotron resonance heating with radio frequency (RF) waves. Ion cyclotron heating (ICH) and electron cyclotron heating in MPEX are used to independently heat the ions and electrons and provide fusion divertor conditions ranging from sheath-limited to fully detached divertor regimes at a material target. A hybrid particle-in-cell code- PICOS++ is developed and applied to understand the plasma parallel transport during ICH in MPEX/Proto-MPEX to the target. With this tool, evolution of the distribution function of MPEX/Proto-MPEX ions is modeled in the presence of (a) Coulomb collisions, (b) volumetric particle sources and (c) quasi-linear RF-based ICH. The code is benchmarked against experimental data from Proto-MPEX and simulation data from B2.5 EIRENE. The experimental observation of ‘density-drop’ near the target in Proto-MPEX and MPEX during ICH is demonstrated and explained via physics-based arguments using PICOS++ modeling. In fact, the density drops at the target during ICH in Proto-MPEX/MPEX to conserve the flux and to compensate for the increased flow during ICH. Furthermore, sensitivity scans of various plasma parameters with respect to ICH power are performed for MPEX to investigate its role on plasma transport and particle and energy fluxes at the target.

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 Ti and target heat flux are observed to scale linearly with injected ICH power. 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.

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

The principal objective of this work is to report on the power coupled to a tungsten target in the prototype-material plasma exposure experiment device during oblique injection of a microwave beam ( n e > 1 × 10 19   m − 3 ) deuterium helicon plasma column. The experimental setup, electron heating system, electron heating scheme, and IR thermographic diagnostic for quantifying the power transport is described in detail. It is demonstrated that the power transported to the target can be effectively controlled by adjusting the magnetic field profile. Using this method, heat fluxes up to 22 MW m⁻² and power transport efficiencies in the range of 17%–20% have been achieved using 70 kW of microwave power. It is observed that most of the heat flux is confined to a narrow region at the plasma periphery. Ray-tracing calculations are presented which indicate that the power is coupled to the plasma electrons via an O-X-B mode conversion process. Calculations indicate that the microwave power is absorbed in a single pass at the plasma periphery via collisions and in the over-dense region via 2nd harmonic cyclotron resonance of the electron Bernstein wave. The impact of these results is discussed in the context of MPEX.

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.