Plasma Heating Research Articles

Overview of the first Wendelstein 7-X long pulse campaign with fully water-cooled plasma facing components

After a long device enhancement phase, scientific operation resumed in 2022. The main new device components are the water cooling of all plasma facing components and the new water-cooled high heat flux divertor units. Water cooling allowed for the first long-pulse operation campaign. A maximum discharge length of 8 min was achieved with a total heating energy of 1.3 GJ. Safe divertor operation was demonstrated in attached and detached mode. Stable detachment is readily achieved in some magnetic configurations but requires impurity seeding in configurations with small magnetic pitch angle within the edge islands. Progress was made in the characterization of transport mechanisms across edge magnetic islands: Measurement of the potential distribution and flow pattern reveals that the islands are associated with a strong poloidal drift, which leads to rapid convection of energy and particles from the last closed flux surface into the scrape-off layer. Using the upgraded plasma heating systems, advanced heating scenarios were developed, which provide improved energy confinement comparable to the scenario, in which the record triple product for stellarators was achieved in the previous operation campaign. However, a magnetic configuration-dependent critical heating power limit of the electron cyclotron resonance heating was observed. Exceeding the respective power limit leads to a degradation of the confinement.

Ultra-low frequency discharge instability of a barium tungsten hollow cathode at lower self-sustaining margin

Development of low-orbiting constellations stimulated the demand for low-power Hall thrusters, as well as their accompanying low-current hollow cathodes. In this paper, such cathodes were tested on quite low currents to search their self-sustaining limit, and were found to exhibit an ultra-low frequency instability at lower current limit. The instability would grow fast when the cathode current was below 0.3A, and oscillate on 0.1Hz to a few Hz with a visibly blinking plume. At the same time, the emitter temperature was fluctuating with a 15 °C peak-to-peak amplitude. The instability was confirmed to originate from the cathode interior, as a result of depleted plasma heating. Then the thermionic cathode had to rely mainly on neutral gas to heat the emitter. Because thermionic emission was in exponential relation with emitter temperature, the temperature fluctuation was then amplified to cause larger fluctuation in total current as well as in total discharge power, which led to deeper fluctuation in temperature. Regarding the heating of both neutral gas and emitter, an instability model showed that multiple unnoticed factors could determine the onset of instability, including the thermal inertia of neutral gas and emitter, the heat transfer from discharging plasma to neutral gas, the heat loss of cathode and the response time of power supply. As a result of this instability, the plume ion energy was elevated from <50eV to 30∼120eV due to the fluctuating density on cathode exit and fluctuating discharge voltage. It was estimated that the increased ion energy could accelerate the keeper erosion by 4∼22 times.

Estimation of neutral beam injector power transferred to KSTAR plasmas

KSTAR shows the long pulse capability based on superconducting magnet and NBI plays a crucial role in sustaining the plasma via plasma heating and current driving. So, the accurate NBI power measurement transferred to KSTAR plasma is important for the analysis of plasma transport as well as plasma performance parameters. In a long pulse operation, the water coolant temperature is at the steady state condition in longer pulse more than 20s and the coupled neutral beam power to KSTAR plasma during the tokamak experiments is rechecked by using water flow calorimetric method after the experiment. The effect of beam duct scraper loss which was not considered at the neutral beam power calibration process is less than 5 % in terms of neutral beam power. However, in long pulse operation of NBI in KSTAR experiments, high strength of stray PF affects the beam path, the neutral beam power registered on mds + using beam current method is over estimated and it is calculated up to a few percent in terms of neutral beam power using calorimetric method. Therefore, it is necessary to consider the beam power variation by PF effect to interpret the plasma performance degradation in long pulse operation. Lastly. when the ion source tun on or turn off in condition other ion sources are operated, the beam transmission power is also affected because of sharing the beam box for ion sources so that the careful power estimation is necessary for such kind of beam power modulation experiments in KSTAR.

Production of X-ray absorbing ceramics and products based on them for radiotherapy of the periorbital eye region

The paper presents an innovative method for producing high-density ceramics based on bioinert tantalum oxide (Ta2O5) and eye shield products utilizing spark plasma sintering (SPS) technology. Eye shields are in high demand for protecting the cornea and other parts of the eye from ionizing radiation during periorbital eye procedures, such as laser therapy, radiotherapy, and Mohs micrographic surgery (MMS). The dynamics of consolidation of microcrystalline Ta2O5 powder under spark plasma heating conditions in the temperature range of 1000–1400 °C were investigated. The structure, phase, and elemental composition of polycrystalline ceramic samples, including the calculation of crystal structure parameters, were studied using XRD, SEM, and EDS methods. The optimal SPS parameters for achieving maximum relative density (up to 100 %) and Vickers microhardness (up to 510 H V) of the obtained ceramics were determined. The performance properties of the ceramics were examined under conditions of a session of short-focus radiation therapy applied to the periorbital facial area. A dosimetric assessment of the absorbing capacity of the ceramics under X-ray irradiation in modes from 40 to 120 keV was conducted. A prototype of a ceramic eye shield product based on Ta2O5 was developed and manufactured for the first time using SPS technology.

Analysis of spiral antenna for enhancing antenna-plasma coupling impedance for SST-1 tokamak

Abstract A detailed characterization of a high-power radio frequency (RF) broadband circularly polarized two-arm spiral antenna is designed to operate within the frequency range of 0.1–1.0 GHz. The impedance matching network technique is introduced to optimize its performance. The traditional spiral antenna is excited by a vertical or horizontal balun, whereas the proposed design is directly fed by a coaxial cable featuring a planar feeding section specially optimized to achieve broadband input impedance matching. The spiral antenna is designed as per the steady-state superconducting tokamak (SST-1) port space constraints. The simulated efficiency of the RF power coupling with the hydrogen plasma is ∼70 %. Through simulation, it was evident that the proposed antenna exhibited inherent resonance at 0.5 GHz with a reflection coefficient of −27.94 dB and an axial ratio is 3.39 dB respectively. The obtained outcomes unequivocally demonstrate the circular polarization of the designed antenna. Overall, the findings support the enhancement of plasma heating and current drive techniques in fusion research.

Tracing field lines that are reconnecting, or expanding, or both

The explosive release of energy in the solar atmosphere is driven magnetically, but the mechanisms that trigger the onset of the eruption remain controversial. In the case of flares and coronal mass ejections (CMEs), ideal or non-ideal instabilities usually occur in the corona, but it is difficult to obtain direct observations and diagnostics there. To overcome this difficulty, we analyze observational signatures in the upper chromosphere or transition region, particularly brightening and dimming at the base of coronal magnetic structures. In this paper, we examine the time evolution of spatially resolved light curves in two eruptive flares and identify a variety of tempo-spatial sequences of brightening and dimming, such as dimming followed by brightening and dimming preceded by brightening. These brightening–dimming sequences are indicative of the configuration of energy release in the form of plasma heating or bulk motion. We demonstrate the potential of using these analyses to diagnose the properties of magnetic reconnection and plasma expansion in the corona during the early stages of the eruption.

A machine-learning augmented mixed method for simulating α particle transport in inertial confinement fusion

The plasma heating by the α particle transport is a main self-heating source in inertial confinement fusion (ICF) that determines the capsule implosion performance. Due to the high energy of α particle and the high temperature in the ICF capsule hot spot, significant non-equilibrium effect exists and the continuum mechanics breaks down. For the numerical simulation of implosion and charged particle transport, the Boltzmann equation needs to be solved to capture the kinetic effects. However, the 7-dimensional Boltzmann equation, the highly frequent Coulomb collision, and the multi-folded integral stopping power formulation greatly limits the computational efficiency and challenges the computational power. To overcome the high computational cost of high-frequent Coulomb collisions, we propose a mixed collision model according to which the collisions are categorized into low-frequent large-angle collisions and high frequent small-angle grazing collisions. The large-angle collision process is precisely solved based on the Coulomb cross-section. For the highly frequent small-angle grazing, a statistic model is constructed with second-order accuracy in time. The mixed collision model reduces the computational cost of scattering calculation by two orders of magnitude. For the multi-folded integral stopping power formulation, a neural network is used to improve computational efficiency. Based on the proposed algorithm, we develop one-dimensional to three-dimensional module code to directly solve the α particle transport Boltzmann equation. The α transport module code is integrated into the multi-physics LARED-S program. The MC version ICF software is verified by a simulation study of the N210207 and N191110 experiment.

Optimization and preparation for the start-up of the plasma ICR heating system at the KTM tokamak

The paper presents and discusses the results of the research on the preparation of the HF plasma heating system of the KTM tokamak for the start-up. The calculations have been performed to determine the parameters at which the most effective absorption of the HF radiation by the plasma of the KTM tokamak was expected, the results of measurements of the frequency response of the antenna-feeder device of the HF heating system generator of the KTM tokamak plasma have been presented. The result of the carried out work have shown that more efficient plasma heating on a KTM tokamak can be obtained at a frequency of 14 MHz.

Thermodynamics of Alfvénic slow solar wind produced by Alfvénic turbulence

ABSTRACT Alfvén-wave turbulence is known as a plasma heating mechanism associated with the acceleration of fast solar wind, found emanating from open magnetic fields adjacent to coronal holes. In this study, we expand the scope of this mechanism to investigate the thermodynamics of Alfvénic slow solar wind, a phenomenon originating from open fields near a streamer, as observed in recent inner heliospheric missions. We demonstrate a one-dimensional two-fluid model that incorporates three components: (1) low-frequency Alfvén-wave turbulence, serving as the primary dissipation mechanism, (2) a curved magnetic field that reproduces the streamer’s boundary, and (3) the kinetic instabilities to address proton temperature anisotropy. Our findings suggest that this dissipation mechanism can be applied in common to both fast and Alfvénic slow solar winds. We identify the proton-cyclotron instability near the Sun and the oblique and parallel firehose instabilities occurring close to 1 au as crucial factors governing temperature anisotropy. This study contributes to our understanding of the complex thermodynamics of solar winds and provides valuable insights for future space missions.

Study of a degenerate non-elliptic equation to model plasma heating

Study of a degenerate non-elliptic equation to model plasma heating

Observation of alpha-particles in recent D–T experiments on JET

The fusion reaction between deuterium and tritium, D(T,n)4 He is the main source of energy in future thermonuclear reactors. Alpha-particles (4 He-ions) born with an average energy of 3.5 MeV transferring energy to the thermal plasma during their slowing down, should provide the self-sustained D–T plasma burn. The adequate confinement of α-particles is essential to provide efficient heating of the bulk plasma and steady burning of a reactor plasma. That is why the fusion-born α-particle studies have been a priority task in the second D–T experiments (DTE2) on the Joint European Torus (JET) to understand the main mechanisms of their slowing down, redistribution and losses and to develop optimal plasma scenarios. JET with Be-wall and W-divertor, enhanced auxiliary heating systems and improved energetic-particle diagnostic capabilities, producing significant population of α-particles, provided the possibility for comprehensive studying of the α-particle behaviour. Selected results of the confined and lost α-particle measurements, evidence of α-particle self-heating and assessments of the fusion performance are presented in this paper giving an opportunity for further modelling and extrapolation to the International Thermonuclear Experimental Reactor and burning plasma reactors.

1D drift-kinetic numerical model based on semi-implicit particle-in-cell method

The paper presents a new one-dimensional drift-kinetic electrostatic model based on the particle-in-cell method and capable of simulating the processes of plasma heating and confinement in mirror traps. The most of particle and energy losses in these traps occur along the magnetic field lines. The key role in limiting these losses is played by the ambipolar electric potential which creates a potential barrier for electrons and significantly reduces the heat flux that, without this barrier, would go to wall due to the classical electron thermal conductivity. However, modeling the formation of such a potential on real spatial and temporal scales of experiments is a challenging problem, since it requires a detailed description of not only ion, but also electron kinetics. In this work, we propose to solve the problem of taking into account electron kinetic effects on the time scale of plasma confinement in a mirror trap using the particle-in-cell method adapted to the approximate drift-kinetic equations of plasma motion. Unlike other electrostatic particle-in-cell models, which use fully implicit schemes to solve the nonlinear system of Vlasov-Poisson and Vlasov-Ampere equations, we propose a semi-implicit approach. By analogy with the Energy Conserving Semi-Implicit Method (ECSIM), it allows for precise conservation of energy and reduces the procedure for finding the electric field to inverting a tridiagonal matrix without multiple nonlinear iterations. Such a model will be useful for simulating not only collisional losses of hot plasma in fusion experiments, but also for studying the features of creating cold starting plasma in mirror traps using plasma or electron guns.

Investigating crater formation in nanosecond laser ablation of aluminum foils

A nanosecond Nd:YAG laser was used to study the laser ablation of aluminum foil in the phase explosion regime at a laser intensity range of 0.63–3.61 ×1012W/cm2. Laser ablation and plasma characteristics were studied using microscopic ablation crater images, plasma emission spectra, and plasma plume images. Measured plasma density using a Stark width of Al I (396.2 nm) showed a strong linear correlation with crater size, with a Pearson correlation coefficient (r) of 0.97. To understand the origin of this linear correlation, plasma temperature was estimated using Bremsstrahlung emission from 512 to 700 nm. The estimated plasma temperature and aspect ratio of the plasma plume were negatively correlated, having r=−0.76. This negative correlation resulted from a laser-plasma interaction, which heated the plasma and increased its hydrodynamic length. The percentages of laser energy used for plasma heating (Ep/EL) and Al foil ablation (EAl/EL) were estimated from plasma temperature. Increased EAl/EL, such as crater size, with increasing laser intensity, confirms that greater mass ablation is the fundamental reason for the strong linear correlation between crater size and plasma density.

Turbulence generation in the solar corona through nonlinear interaction between 3D kinetic Alfvén wave and magnetic islands

This study investigates the interaction between kinetic Alfvén waves (KAWs) and magnetic islands in solar corona region. KAWs and magnetic islands are believed to be one of the key players in the heating of the coronal plasma and particle acceleration. We have used the dynamical model for KAWs to simulate its propagation through magnetic islands. The equation is solved numerically by using the finite difference and pseudospectral technique for spatial and temporal variations. The simulation results have been studied and signatures of turbulent structures evolving with time have been observed. Turbulence plays a considerable role in tapping out magnetic energy to thermal energy. The plot of the energy spectrum as a function of wavenumber exhibits a power law behavior in the inertial range, characterized with an exponent around −5/3 . Beyond this range, the spectrum becomes steeper. In order to gain a more thorough knowledge of the physical processes that underlie the development of localized structures and to estimate the current sheet scale sizes, a semi-analytically model is also performed. The results indicate that the size of localized structures are comparable to ion gyro-radius, and different intensities of KAW further modify this length scale formed at different times. The purpose of the study is to acquire an understanding of the KAW propagation in magnetic islands and the way it contributes to the formation of turbulence and energy release in the solar corona. A comparative aspect of the non-linear interaction of 3D KAW with magnetic island and null point is also presented which shows that non-linear interaction of 3D KAW with null point, leads to rapid disorganization of the magnetic field and the creation of localized structures.

Plasma heating characterization of the large area inductively coupled plasma etchers with the plasma information for managing the mass production

Meter-scale of the large area inductively coupled plasma etchers with the capacitive power coupling are widely applied for the mass production of OLED (organic light emitting diode) display panels. Because of the large area-to-volume ratio of the etcher, the balance between the power loss and absorption is easily located in the capacitive coupling mode rather than the ideal inductively coupled mode. Therefore, the process results are sensitively governed by the power absorption and plasma heating properties of the reactors. We have introduced a new PI (plasma information) parameter, the ratio of the stochastic heating to Ohmic heating of the plasmas, which is monitorable by using the optical emission spectroscopy data of the processing etchers. With the help of this plasma heating characteristic index, we could optimize the process recipes with the detailed control of the etched hole sidewall passivation and related species generation rate in the plasmas; thus, chamber-to-chamber matching in the huge mass production fab with the higher efficiency was possible. It was demonstrated that the introduced PI index with plasma heating mechanism characterization could be applicable to the VM (virtual metrology) modeling as one of the good information supplying core variables. This PI index has shown a very high correlation with the plasma sheath and ion flux governing phenomena for a large number of mass-produced OLED display glasses. From these results, the introduced plasma heating mechanism-based PI index is expected to be utilized as a good reference index for their performance analysis or PI-VM modelings.

Characterization of water propellant in an electron cyclotron resonance thruster

A coaxial electron cyclotron resonance thruster operating on water vapor propellant is investigated to determine the influence of molecular propellant chemistry on thruster performance. The performance is characterized at different mass flow rates (0.1–0.4 mg/s) and powers (20–200 W) using a thrust stand, ion energy analyzer, and spectrometer. Experimental data are compared to the results from a theoretical model that includes non-equilibrium chemical kinetics. The thruster is observed to transition from a state where power deposition favors propellant dissociation and ionization into a state where it favors electron heating and ion acceleration. The results suggest that the majority of plasma heating occurs in a core region that surrounds the thruster antenna and that propellant flow outside this region is not efficiently utilized for propulsion.

Detection of alpha heating in JET-ILW DT plasmas by a study of the electron temperature response to ICRH modulation

In the JET DTE2 campaign a new method was successfully tested to detect the heating of bulk electrons by α-particles, using the dynamic response of the electron temperature T e to the modulation of ion cyclotron resonance heating (ICRH). A fundamental deuterium (D) ICRH scheme was applied to a tritium-rich hybrid plasma with D-neutral beam injection (NBI). The modulation of the ion temperature T i and of the ICRH accelerated deuterons leads to modulated α-heating with a large delay with respect to other modulated electron heating terms. A significant phase delay of ∼40° is measured between central T e and T i, which can only be explained by α-particle heating. Integrated modelling using different models for ICRH absorption and ICRH/NBI interaction reproduces the effect qualitatively. Best agreement with experiment is obtained with the European Transport Solver/Heating and Current Drive workflow.

Heat and particle exhaust in high-performance plasmas in Wendelstein 7-X

The paper reports for the first time the heat and particle exhaust at the plasma boundary through various edge diagnostics for the high-performance plasma obtained after pellet injection on Wendelstein 7-X. The plasma density at the edge is found to be reduced by a factor of 2 in the high-performance phase, supporting the previously reported density peaking at the plasma centre. The plasma beta effect on the magnetic topology is reflected by the appearance of the second strike line, which is well understood with simulation. However, during the rapid decay phase of the enhanced confinement, a transient localized heat flow of up to 16 MW m−2 is observed at the leading edge of a poorly cooled divertor component, which has not been understood but raises concerns about machine safety.

Does high-latitude ionospheric electrodynamics exhibit hemispheric mirror symmetry?

Abstract. Ionospheric electrodynamics is a problem of mechanical stress balance mediated by electromagnetic forces. Joule heating (the total rate of frictional heating of thermospheric gases and ionospheric plasma) and ionospheric Hall and Pedersen conductances comprise three of the most basic descriptors of this problem. More than half a century after identification of their central role in ionospheric electrodynamics, several important questions about these quantities, including the degree to which they exhibit hemispheric symmetry under reversal of the sign of dipole tilt and the sign of the y component of the interplanetary magnetic field (so-called “mirror symmetry”), remain unanswered. While global estimates of these key parameters can be obtained by combining existing empirical models, one often encounters some frustrating sources of uncertainty: the measurements from which such models are derived, usually magnetic field and electric field or ion drift measurements, are typically measured separately and do not necessarily align. The models to be combined moreover often use different input parameters, different assumptions about hemispheric symmetry, and/or different coordinate systems. We eliminate these sources of uncertainty in model predictions of electromagnetic work J⋅E (in general not equal to Joule heating ηJ2) and ionospheric conductances by combining two new empirical models of the high-latitude ionospheric electric potential and ionospheric currents that are derived in a mutually consistent fashion: these models do not assume any form of symmetry between the two hemispheres; are based on Apex magnetic coordinates (denoted Apex), spherical harmonics, and the same model input parameters; and are derived exclusively from convection and magnetic field measurements made by the Swarm and CHAMP satellites. The model source code is open source and publicly available. Comparison of high-latitude distributions of electromagnetic work in each hemisphere as functions of dipole tilt and interplanetary magnetic field clock angle indicates that the typical assumption of mirror symmetry is largely justified. Model predictions of ionospheric Hall and Pedersen conductances exhibit a degree of symmetry, but clearly asymmetric responses to dipole tilt and solar wind driving conditions are also identified. The distinction between electromagnetic work and Joule heating allows us to identify where and under what conditions the assumption that the neutral wind corotates with the Earth is not likely to be physically consistent with predicted Hall and Pedersen conductances.

Cascade-Dissipation Balance in Astrophysical Plasmas: Insights from the Terrestrial Magnetosheath.

The differential heating of electrons and ions by turbulence in weakly collisional magnetized plasmas and the scales at which such energy dissipation is most effective are still debated. Using a large data sample measured in Earth's magnetosheath by the magnetospheric multiscale mission and the coarse-grained energy equations derived from the Vlasov-Maxwell system, we find evidence of a balance over two decades in scales between the energy cascade and dissipation rates. The decline of the cascade rate at kinetic scales (in contrast with a constant one in the inertial range), is balanced by an increasing ion and electron heating rates, estimated via the pressure strain. Ion scales are found to contribute most effectively to ion heating, while electron heating originates from both ion and electron scales. These results can potentially impact the current understanding of particle heating in turbulent magnetized plasmas as well as their theoretical and numerical modeling.