Plasma Low Density Research Articles

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

Numerical study of minority ion heating scenarios in CN-H1 stellarator plasma

Abstract A preliminary simulation of ion cyclotron resonance heating (ICRH) in the poloidal cross-section at the antenna location of the CN-H1 (Chinese Heliac 1) stellarator has been conducted for the first time using the full-wave solver TORIC. The heating scheme employed focuses on minority ion heating of He4 (H), with an independent investigation into the effects of wave frequency, toroidal mode number, and minority ion concentration on ICRH. Preliminary simulation results indicate that, under the conditions of wave frequency f = 4.7MHz, central magnetic field B0 = 0.32T and low-temperature, low-density plasma (with density on the order of 1018 m-3 and temperature in the range of tens of eV), more efficient minority-ion heating can be achieved when the toroidal mode number is in the range of 10–20 and the minority ion concentration is between 15% and 30%. This simulation provides theoretical insights and delineates the parameter space for prospective ICRH experiments in CN-H1.

Non-linear dependence of ion heat flux on plasma density at the L–H transition of JET NBI-heated deuterium–tritium plasmas

Abstract Recent JET D–T campaigns opened the possibility of unique isotope studies to investigate the L–H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L–H transitions in the low and high-density branches of JET NBI-heated D–T plasmas. As discussed in the paper, L–H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux (Q i). We present the first power balance analysis of JET NBI-heated D–T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D–T Q i, with a lower Q i in D–T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L–H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Q i deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Q i and the power threshold minimum at JET. At low plasma density, NBI power dominates Q i, while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D–T NBI-heated plasmas at JET.

Suppression of core temperature fluctuations by edge cooling in the J-TEXT tokamak

Abstract The suppression of core temperature fluctuations due to edge cooling in the J-TEXT tokamak is presented. The supersonic molecular beam is injected into the low density plasmas and the core temperature rise appears. The staggered time correlation technique is utilized to infer the temperature fluctuations from an electron cyclotron emission radiometer (ECE). The existing correlation ECE (CECE) also enables the observation of the temperature fluctuations simultaneously. A significant flattening of the density gradients appears almost in the whole plasmas after the cold pulse injection. The behaviors of the temperature fluctuations from both ECE and CECE show a good agreement, except for the fluctuation level, which exhibits a relative variation of approximately 10% − 50%. The relative intensity of temperature fluctuations decreases within the surface of r / a ∼ 0.8 due to the edge cooling. The cross-power analysis and the agreement in measurements suggest that temperature fluctuations from turbulence inferred from a single ECE channel may be possible. Linear simulations suggest the drop in temperature fluctuation intensity may relate to the trapped electron mode turbulence stabilization. This paper also demonstrates the evolution of both gradients and fluctuations before and after supersonic molecular beam injection with density scans.

The slow wave resonance cone in the collisional regime

In the low-density edge plasma of tokamaks, ion cyclotron range of frequencies actuators may parasitically emit slow waves. If the density is sufficiently low, which may be common in large future devices such as international thermonuclear experimental reactor (ITER), these slow waves take the form of so-called resonance cones. The traditional theoretical description of this wave mode relies on formally relating an electrostatic approximation of the frequency-domain wave equation to a time-domain wave equation and relating the cone angle to the wave speed in the time-domain wave equation. In the cold plasma collisional regime, that wave speed is complex. We investigate that scenario in this work.

Impact of sample preparation on bitumen content measurement using laser-induced breakdown spectroscopy

The impact of sample preparation on bitumen content measurement using LIBS was investigated by collecting spectra from wet and dry tailings. A multivariate data analysis model was developed using optimal wavelength selection for bitumen content classification and prediction in tailings. Wet tailings can be classified into three classes (low, medium, and high bitumen) with 12.1 % error, while dry tailings have a classification error of 6.1 %. Quantitative analysis showed a bitumen content prediction error of 4.7 % for wet tailings and 8.9 % for dry tailings. Wet tailings showed a 1.8–2.5 times improvement in the limit of detection range compared to dry tailings. Plasma density and crater size measurements revealed that plasma density fluctuation was 2.7 times lower in wet tailings due to consistent crater formation from laser-tailings interaction. The lower plasma density fluctuation indicates a stable mass ablation for wet samples, which is attributed as the primary reason for significant LIBS performance improvement on wet tailings.

Divertor Tokamak Test: Impact of NBI shine-through and beam-plasma interaction on Divertor Tokamak Test facility

IntroductionIn this work, we aim to explore numerically the behavior of beam energetic particles in the Divertor Tokamak Test (DTT), a superconductive device equipped with a Neutral Beam Injection (NBI) system capable of injecting neutrals up to 510 keV.MethodWe explore beam ionization and beam slowing down for different DTT plasma scenarios. Numerical simulations are performed using the ASCOT suite of codes, including a wide-range scan of plasma density and beam injection energy. For different plasma conditions, we estimate shine-through losses, including the heat fluxes on the first wall thanks to dedicated particle tracing simulations. Orbits of newly-born fast ions are characterized by means of the constant of motion phase space, showing how trapped energetic particles’ population and prompt losses change with plasma density and NBI energy.Results and discussionSlowing down simulations show that NBI injection at 510 keV is well coupled to DTT plasmas. DTT NBI will be one of the sources of auxiliary ion heating, with an absorbed power ratio of up to ∼50% depending on plasma and beam parameters. At low plasma densities, energetic particle confinement is less efficient, and NBI power and/or energy reduction is expected.

Inductively Coupled Plasma Reactive Ion Etching of (−201) β‐Ga2O3 Nano‐Fins

This study reports the effect of inductively coupled plasma reactive ion etching using a BCl3/Ar gas mixture on the sidewall taper angle (STA), etch rate, and surface morphology of (−201) β‐Ga2O3 substrates, as well as (−201)‐oriented β‐Ga2O3 films on sapphire, for nano‐fins of ≈200 nm width. Various etch parameters are systematically investigated for the β‐Ga2O3 films on sapphire, revealing that the STA is lower at high plasma density when chemical component increases, and the ion scattering decreases; achieving a STA of 1.5° ± 0.7°, accompanied by an etch rate of 114 nm min−1 and a reduction in surface roughness compared to before etching. Additionally, a strong positive correlation between the STA and the crystal orientation in the β‐Ga2O3 substrates is observed at lower plasma density, while a weak correlation is found between the STA and the crystal orientation at high plasma density. This study provides valuable insights for the fabrication of β‐Ga2O3 devices.

Modelling study of divertor W leakage for different divertor conditions with Ne seeding in EAST tokamak

Abstract Sequential SOLPS-ITER and DIVIMP simulations are carried out for tungsten (W) edge transport during neon (Ne) injection in EAST, elucidating the pathway of W impurity leakage from the divertor to the core plasma under different dissipative divertor conditions. Divertor conditions from low-recycling to detachment are generated by modulating the Ne injection rate in the SOLPS simulation. With the drift velocities incorporated in DIVIMP, W transport under the drift effect is investigated. For low-density L-mode plasma conditions with favorable Bt direction, the majority of W source originates from the outer divertor target and leaks upstream through the near-SOL EB drift reversed region at both inner and outer divertors. The EB drift is demonstrated to have a significant effect on W leakage and the effect diminishes with the increase of Ne injection rate. Compared to the D2 injection cases, Ne injection helps to achieve a strong dissipative divertor condition at a much lower plasma density, resulting in a higher W leakage ability due to the smaller friction with the background plasma. For Ne injection cases under high-density H-mode conditions, W leaks mostly through the outer divertor region while the inner divertor is fully detached. Compared to the L-mode cases, the shorter power decay lengths in the scrap-off layer (SOL) of the H-mode cases result in smaller parallel EB drift velocities especially in the middle and far SOL region. With the increase of the Ne injection rate, the decrease of the EB drift in the middle SOL leads to a wider near-separatrix W leakage path. Therefore, a worse divertor W screening is expected, which is consistent with the previous published experimental observations [1].

Ion velocity separation mechanism during vacuum spark stage

Abstract Supersonic ion jets produced in vacuum arc discharges have a wide range of applications, where precise control of ion kinetic energy is crucial. However, a comprehensive understanding of the ion acceleration mechanism remains elusive, particularly regarding whether there is ion velocity separation in the vacuum spark stage. In this paper, a 1D spherical implicit particle-in-cell (PIC) with Monte Carlo collision (MCC) model is employed to investigate the ion velocity separation in multi-charged vacuum arc plasma with varying electrode bias voltages and plasma ion densities. The results show that ion kinetic energy can reach hundreds of electron volts due to continuous acceleration by the formed potential valley, which leads to ion velocity separation at low electrode bias voltage or low plasma density. An increasing electrode bias voltage flattens the potential valley, reducing the electric field acceleration. While increasing the plasma density deepens the valley and intensifies Coulomb collisions, resulting in nearly-equal velocities across ions in different charge states. These findings can theoretically explain the discrepancies observed in previous experiments regarding the dependence of the ion velocity on its charge state during the vacuum spark stage.

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρi/a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρi effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

A model for fast electron-driven high-density plasma in the double-cone ignition scheme

Abstract A model for fast electron-driven high-density plasma is proposed to describe the effect of injected fast electrons on the temperature and inner pressure of the plasma in the fast heating process of the double-cone ignition (DCI) scheme. Due to the collision of the two low-density plasmas, the density and volume of the high-density plasma vary. Therefore, the ignition temperature and energy requirement of the high-density plasma vary at different moments, and the required energy for hot electrons to heat the plasma also changes. In practical experiments, the energy input of hot electrons needs to be considered. To reduce the energy input of hot electrons, the optimal moment and the shortest time for injecting hot electrons with minimum energy are analyzed. In this paper, it is proposed to inject hot electrons for a short time to heat the high-density plasma to a relatively high temperature. Then, the alpha particles with the high heating rate and PdV work heat the plasma to the ignition temperature, further reducing the energy required to inject hot electrons. The study of the injection time of fast electrons can reduce the energy requirement of fast electrons for the high-density plasma and increase the probability of successful ignition of the high-density plasma.

Electron density measurement of an ionospheric-like plasma using an impedance probe.

In this work, an impedance probe was designed to measure ionospheric low-density plasma. The probe is made of short cylindrical dipoles with an antenna 20cm in length and 0.25cm in diameter, which can operate in the frequency range of 0-100 MHz. Combined with a vector network analyzer, the performance of the probe was tested in laboratory-created ionospheric-like plasma, and the results were compared with those of the Langmuir probe measurements. The densities measured by the two methods show a consistent trend, and the impedance probe data show much smaller uncertainties, which suggests that the impedance probe can achieve high-precision measurements. Furthermore, by fitting the antenna phase data, the electron-neutral collision frequency can also be obtained. Therefore, the impedance probe provides a feasible method for exploring low-density partially ionized plasma and can be expanded to actual ionospheric exploration in the future.

Micropinch formation dynamics in X pinches.

High temporal resolution x-ray streak camera studies of micropinch formation in Cu hybrid xpinches reveal key plasma conditions. Analysis of Ne-like Cu lines indicate an average electron temperature of about 200eV and 4.5×10^{28}m^{-3} electron density. The spectra suggest that the electron temperature jumps to about 1 keV, inferred from the continuum and the postcontinuum line emission that includes Li-like Cu lines. There is no sign of a rapid temperature change or a substantial surge in radiation emission during the 200 ps precontinuum x-ray burst, suggesting that the radiative collapse process does not play a major role in micropinch formation. Two-dimensional extended Magnetohydrodynamic (MHD) simulations, coupled to a collisional-radiative spectral analysis code, suggest the significance of the rapid radial implosion of high-temperature, low-density plasma, the axial outflow, and the dynamic plasma pressure in micropinch formation.

Toroidal Alfvén eigenmodes excited by energetic electrons in EAST low-density Ohmic plasmas

Abstract Operation in the quiescent regime with abundant trapped energetic electrons (EEs) has been achieved during the current flattop in EAST low-density Ohmic plasmas. This was facilitated by increasing the electron density to a specified level and subsequently reducing it slowly, resulting in the accumulation of a sufficient number of trapped EEs within the energy range of 150–250 keV. During the phase of decreasing electron density, toroidal Alfvén eigenmodes (TAEs) were observed to be excited by these EEs, with frequencies falling within the range of about 100–300 kHz. The experimental parameters were carefully set to satisfy the resonance conditions for TAE excitation by EEs, aligning well with predictions from ideal MHD theory. Statistical analysis indicated different density dependencies between the frequencies of TAEs and the Alfvén frequencies, due to their different radial excitation positions. The radial positions of the TAEs were found to be influenced by the energy distribution and the evolution of trapped EEs, which in turn were affected by the decay rate of electron density and loop voltage. Measurements of Hard X-rays (HXR) confirmed an energy distribution characterized by a “bump-on-tail” shape, with the TAEs observed near the energy bump. Theoretical considerations also demonstrate the possibility that the e-TAE can be driven unstable under this experiment condition even if the mode does not rotate in the electron-diamagnetic drift direction.

Toroidal Alfvén mode instability driven by plasma current in low-density Ohmic plasmas of the spherical tori

Due to intrinsically low magnetic fields, at low density the drift speed of the current-carrying electrons in spherical tokamaks can exceed fraction of the Alfvén speed sufficient for the excitation of the Alfvén gap modes. A particular case of the toroidal mode, observed during minor disruptions in Ohmic shots on SUNIST [Liu et al., Phys. Plasmas 23, 120706 (2016)], is considered in the present communication. Due to the negligible effect of the electron pressure gradient, the growth rate scales linearly with the drift speed, with slope inversely proportional to electron thermal velocity. In the absence of continuum damping, the threshold value of the drift speed for TAE excitation is independent of electron temperature.

Observation of tungsten emission spectra up to W46+ ions in the Large Helical Device and contribution to the study of high-Z impurity transport in fusion plasmas

Spectroscopic studies of emissions released from tungsten ions combined with a pellet injection technique have been conducted in Large Helical Device for contribution to the tungsten transport study in tungsten divertor fusion devices and for expansion of the experimental database of tungsten line emissions. The spectral intensities of W5+, W24+–W28+, W37+, W38+, W41+–W43+, W45+, and W46+ emission lines were measured simultaneously over a wide wavelength range from x-ray to visible. Time evolutions of the various tungsten line spectra indicate that the tungsten confinement time depends on the electron density of the plasma and is long in high density plasmas, on the order of seconds, and short in low density plasmas, on the order of sub-seconds. When the confinement time was long, the tungsten ions remained in the plasma until the end of the discharge, changing their dominant charge with the change in electron temperature. When the confinement time was short, the tungsten ions rapidly decreased in all charge states and disappeared. Space-resolved EUV and visible spectroscopy measurements have revealed that tungsten ions stayed in the core region of the plasma with changing their dominant charge state depending on the electron temperature in the discharges with the long confinement time. Detailed analysis of soft x-ray emission suggested that the confinement time increases with density and becomes saturated when the central electron density exceeds 2 × 1013 cm−3.

A study of arc plasma energy deposition flow control on blunt-headed protuberance induced STBLI

This paper presents experimental results of flow control method using plasma energy deposition on shock wave interactions induced by blunt-headed protuberance on a flat plate at Mach 2.0. Upstream of the leading edge of the protuberance, plasma energy deposition was generated by high-frequency pulsed discharge at three streamwise locations on the flat plate. High-speed Schlieren technique was utilized to investigate the unsteadiness characteristics of shock wave/turbulent boundary layer interaction (STBLI). The visualization results reveal that thermal bubble structures formed by the pulsed discharge could attenuate the separation shock continuously, resulting in a 12.8% decrease in the maximum pressure in the interaction wall. Experimental findings disclosed the dependency of this mitigation efficacy on the average discharge power. The pulsation of the separation shock was also affected by average discharge power, while the pulsation of boundary layer and shear layer was linked to both the discharge power and frequency. The numerical results predicted flow topology under control well with experiment results, low-density plasma layer was formed by high-frequency discharge, which weakening the separation shock substantially, the “λ” shock intensity was also decreased. The flow field of numerical simulation and experiments are in good agreement. Additionally, the control of the thermal bubble structures on STBLI was very effective in terms of reducing the motion and weakening the intensity of the low-frequency separation shock, and augmenting the high frequency component in the turbulent boundary layer and shear layer.

Effect of multi-cusp magnetic fields to generate a high-density hydrogen plasma inside a low pressure H2 cylindrical hollow cathode discharge

Based on probe measurements, the plasma density in a low pressure H2 radio-frequency (RF) magnetized cylindrical hollow-cathode capacitively coupled plasma (CCP) is demonstrated to be enhanced by applying a multi-cusp magnetic field (MCMF). CCPs can be used for fabricating functional thin-films, but are currently limited by low plasma densities of less than or about 1016 m−3. It is found that a maximum plasma density of 2.51, 2.88, and 3.14 × 1016 m−3 is obtained at 3, 4 and 5 Pa, respectively. The employed MCMF arrangement contributes to achieving superior plasma uniformity, both in axial and radial direction, at an RF power of 50 W.

Interactions between MSTIDs and Ionospheric Irregularities in the Equatorial Region Observed on 13–14 May 2013

We investigate the interactions between medium-scale traveling ionospheric disturbances (MSTIDs) and the equatorial ionization anomaly (EIA) as well as between MSTIDs and equatorial plasma bubbles (EPBs) on the night of 13–14 May 2013, based on observations from multiple instruments (an all-sky imager, digisonde, and global positioning system (GPS)). Two dark bands (the low plasma density region) for the MSTIDs were observed moving toward each other, encountering and interacting with the EIA, and subsequently interacting again with the EIA before eventually dissipating. Then, a new dark band of MSTIDs moved in the southwest direction, drifted into the all-sky imager’s field of view (FOV), and interacted with the EIA. Following this interaction, a new dark band split off from the original dark band, slowly moved in the northeast direction, and eventually faded away in a short time. Subsequently, the original southwestward-propagating dark band of the MSTIDs encountered eastward-moving EPBs, leading to an interaction between the MSTIDs and the EPBs. Then, the dark band of the MSTIDs faded away, while the EPBs grew larger with a pronounced westward tilt. The results from various observational instruments indicate the pivotal role played by the high-density region of the EIA in the occurrence of various interaction processes. In addition, this study also revealed that MSTIDs propagating into the equatorial region can significantly impact the morphology and evolution characteristics of EPBs.