Plasma Science and Technology

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s⁻¹, agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

The dielectric barrier discharge (DBD) in air at atmospheric pressure is not suitable for industrial applications due to its randomly distributed discharge filaments. In this paper, the influence of the electric field distribution on the uniformity of DBD is theoretically analyzed and experimentally verified. It is found that a certain degree of uneven electric field distributions can control the development of electron avalanches and regulate their transition to streamers in the gap. The discharge phenomena and electrical characteristics prove that an enhanced Townsend discharge can be formed in atmospheric-pressure air with a curved-plate electrode. The spectral analysis further confirms that the gas temperature of the plasma produced by the curved-plate electrode is close to room temperature, which is beneficial for industrial applications. This paper presents the relationship between the electron avalanche transition and the formation of a uniform DBD, which can provide some references for the development and applications of the DBD in the future.

High fusion triple product has been obtained in the advanced scenarios with high normalized beta (β_N) on the Experimental Advanced Superconducting Tokamak (EAST). A record value of n_i0T_i0τ_E ∼ 1.0 × 10¹⁹ m⁻³ keV s for EAST deuterium plasma has been achieved, which is due to the formation of strong and broad internal transport barriers (ITBs) in n_e, T_e and T_i profiles. Analysis shows that the strong ITB formation could be attributed to the reduction of transport from ITG modes. Based on the analysis, the physical mechanisms and methods to further improve the plasma performance are discussed.

An impurity powder dropper was installed in the 21st campaign of the Large Helical Device experiment (Oct. 2019–Feb. 2020) under a collaboration between the National Institute for Fusion Science and the Princeton Plasma Physics Laboratory for the purposes of real-time wall conditioning and edge plasma control. In order to assess the effective injection of the impurity powders, spectroscopic diagnostics were applied to observe line emission from the injected impurity. Thus, extreme-ultraviolet (EUV) and vacuum-ultraviolet (VUV) emission spectra were analyzed to summarize observable impurity lines with B and BN powder injection. Emission lines released from B and N ions were identified in the EUV wavelength range of 5–300 Å measured using two grazing incidence flat-field EUV spectrometers and in the VUV wavelength range of 300–2400 Å measured using three normal incidence 20 cm VUV spectrometers. BI–BV and NIII–NVII emission lines were identified in the discharges with the B and BN powder injection, respectively. Useful B and N emission lines which have large intensities and are isolated from other lines were successfully identified as follows: BI (1825.89, 1826.40) Å (blended), BII 1362.46 Å, BIII (677.00, 677.14, 677.16) Å (blended), BIV 60.31 Å, BV 48.59 Å, NIII (989.79, 991.51, 991.58) Å (blended), NIV 765.15 Å, NV (209.27, 209.31) Å (blended), NVI 1896.80 Å, and NVII 24.78 Å. Applications of the line identifications to the advanced spectroscopic diagnostics were demonstrated, such as the vertical profile measurements for the BV and NVII lines using a space-resolved EUV spectrometer and the ion temperature measurement for the BII line using a normal incidence 3 m VUV spectrometer.

The physics of electrodeless electric thrusters that use directed plasma to propel spacecraft without employing electrodes subject to plasma erosion is reviewed. Electrodeless plasma thrusters are potentially more durable than presently deployed thrusters that use electrodes such as gridded ion, Hall thrusters, arcjets and resistojets. Like other plasma thrusters, electrodeless thrusters have the advantage of reduced fuel mass compared to chemical thrusters that produce the same thrust. The status of electrodeless plasma thrusters that could be used in communications satellites and in spacecraft for interplanetary missions is examined. Electrodeless thrusters under development or planned for deployment include devices that use a rotating magnetic field; devices that use a rotating electric field; pulsed inductive devices that exploit the Lorentz force on an induced current loop in a plasma; devices that use radiofrequency fields to heat plasmas and have magnetic nozzles to accelerate the hot plasma and other devices that exploit the Lorentz force. Using metrics of specific impulse and thrust efficiency, we find that the most promising designs are those that use Lorentz forces directly to expel plasma and those that use magnetic nozzles to accelerate plasma.

In this paper we present a new experimental observation using a conventional reflectometry technique, poloidal correlation reflectometry (PCR), in the Experimental Advanced Superconducting Tokamak (EAST). The turbulence spectrum detected by the PCR system exhibits an asymmetry and induced Doppler shift ${f}_{\mathrm{D}}$ during the internal kink mode (IKM) rotation phase. This Doppler shift ${f}_{\mathrm{D}}$ is the target measurement of Doppler reflectometry, but captured by conventional reflectometry. Results show that the Doppler shift ${f}_{\mathrm{D}}$ is modulated by the periodic changes in the effective angle between the probing wave and cutoff layer normal, but not by plasma turbulence. The fishbone mode and saturated long-lived mode are typical IKMs, and this modulation phenomenon is observed in both cases. Moreover, the value of the Doppler shift ${f}_{\mathrm{D}}$ is positively correlated with the amplitude of the IKM, even when the latter is small. However, the positive and negative frequency components of the Doppler shift ${f}_{\mathrm{D}}$ can be asymmetric, which is related to the plasma configuration. A simulated analysis is performed by ray tracing to verify these observations. These results establish a clear link between ${f}_{\mathrm{D}}$ and IKM rotation, and are helpful for studying the characteristics of IKM and related physical phenomena.

A gas puff imaging (GPI) diagnostic has been developed and operated on EAST since 2012, and the time-delay estimation (TDE) method is used to derive the propagation velocity of fluctuations from the two-dimensional GPI data. However, with the TDE method it is difficult to analyze the data with fast transient events, such as edge-localized mode (ELM). Consequently, a method called the spatial displacement estimation (SDE) algorithm is developed to estimate the turbulence velocity with high temporal resolution. Based on the SDE algorithm, we make some improvements, including an adaptive median filter and super-resolution technology. After the development of the algorithm, a straight-line movement and a curved-line movement are used to test the accuracy of the algorithm, and the calculated speed agrees well with preset speed. This SDE algorithm is applied to the EAST GPI data analysis, and the derived propagation velocity of turbulence is consistent with that from the TDE method, but with much higher temporal resolution.

A multi-channel polarimeter–interferometer has been developed on the Keda Torus eXperiment (KTX) for the study of equilibrium dynamics and internal magnetic fluctuations. A three-wave technique based on terahertz solid-state sources (~650 GHz) is applied for simultaneous measurements of electron density and Faraday rotation angle. The output power of the microwave source is 2 mW. Faraday rotation effect using a rotating wave plate is tested with phase noise less than 0.8°, and the density phase noise is less than 0.9°. Measurement of Faraday rotation angle and density for discharges on KTX have demonstrated high sensitivity to internal MHD activities.

A high energy and high yield neutron source is a prime requirement for technological studies related to fusion reactor development. It provides a high-energy neutron environment for small-scale fusion reactor components research and testing such as tritium breeding, shielding, plasma-facing materials, reaction cross-section data study for fusion materials, etc. Along with ITER participation, the Institute of Plasma Research, India is developing an accelerator-based 14 MeV neutron source with a yield of 10¹² n s⁻¹. The design of the source is based on the deuterium–tritium fusion reaction. The deuterium beam is accelerated and delivered to the tritium target to generate 14 MeV neutrons. The deuterium beam energy and tritium availability in the tritium target are the base parameters of the accelerator-based neutron source design. The paper gives the physics design of the neutron generator facility of the Institute for Plasma Research. It covers the requirements, design basis, and physics parameters of the neutron generator. As per the analytical results generator can produce more than 1 × 10¹² n s⁻¹ with a 110 keV D⁺ ion beam of 10 mA and a minimum 5 Ci tritium target. However, the detailed simulation with the more realistic conditions of deuteron ion interaction with the tritium titanium target shows that the desired results cannot be achieved with 110 keV. The safe limit of the ion energy should be 300 keV as per the simulation. At 300 keV ion energy and 20 mA current, it reaches 1.6 × 10¹² n s⁻¹. Moreover, it was found that to ensure sufficiently long operation time a tritium target of more than 20 Ci should be used. The scope of the neutron source is not limited to the fusion reactor research studies, it is extended to other areas such as medical radioisotopes research, semiconductor devices irradiations, and many more.

Streamers represent an important stage in the initiation of gap discharge. In this work, we used an eight-frame intensified charge-coupled device camera to capture the streamer development process when a lightning impulse voltage of 95%–100% U_50% was applied in a 3 m rod–plate gap and the streamer velocity was analyzed. Analysis of the observations shows that streamer velocity can be defined by three stages: rapid velocity decline (stage 1), rapid velocity rise (stage 2) and slow velocity decline (stage 3). The effects of electrode shape, applied voltage and gap breakdown or withstanding on streamer velocity were analyzed. The electrode with a larger radius of curvature will result in a higher initial velocity, and a higher voltage amplitude will cause the streamer to propagate faster at stage 3. Gap withstanding or breakdown has no obvious effect on streamer velocity. In addition, the experimental results are compared with previous results and the statistical characteristics of the primary streamer discharge are discussed.

In this study, a pulsed, high voltage driven hollow-cathode electron beam sources through an optical trigger is designed with characteristics of simple structure, low cost, and easy triggering. To validate the new design, the characteristics of hollow-cathode discharge and electron beam characterization under pulsed high voltage drive are studied experimentally and discussed by discharge characteristics and analyses of waveform details, respectively. The validation experiments indicate that the pulsed high voltage supply significantly improves the frequency and stability of the discharge, which provides a new solution for the realization of a high-frequency, high-energy electron beam source. The peak current amplitude in the high-energy electron beam increases from 6.2 A to 79.6 A, which indicates the pulsed power mode significantly improves the electron beam performance. Besides, increasing the capacitance significantly affects the high-current, lower-energy electron beam more than the high-energy electron beam.

To guide the illuminating design to improve the on-state performances of gallium arsenide (GaAs) photoconductive semiconductor switch (PCSS), the effect of spot size on the operation mode of GaAs PCSS based on a semi-insulating wafer with a thickness of 1 mm, triggered by a 1064-nm extrinsic laser beam with the rectangular spot, has been investigated experimentally. It is found that the variation of the spot size in length and width can act on the different parts of the output waveform integrating the characteristics of the linear and nonlinear modes, and then significantly boosts the PCSS toward different operation modes. On this basis, a two-channel model containing the active and passive parts is introduced to interpret the relevant influencing mechanisms. Results indicate that the increased spot length can peak the amplitude of static domains in the active part to enhance the development of the nonlinear switching, while the extended spot width can change the distribution of photogenerated carriers on both parts to facilitate the linear switching and weaken the nonlinear switching, which have been proved by comparing the domain evolutions under different spot sizes.

Diesel particulate matter (DPM) and hydrocarbons (HCs) emitted from diesel engines have a negative affect on air quality and human health. Catalysts for oxidative removal of DPM and HCs are currently used universally but their low removal efficiency at low temperatures is a problem. In this study, Cu-doped CeO₂ loaded on Al₂O₃ coupled with plasma was used to enhance low-temperature oxidation of DPM and HCs. Removals of DPM and HCs at 200 °C using the catalyst were as high as 90% with plasma but below 30% without plasma. Operando plasma diffuse reflectance infrared Fourier transform spectroscopy coupled with mass spectrometry was conducted to reveal the functional mechanism of the oxygen species in the DPM oxidation process. It was found that Cu–CeO₂ can promote the formation of adsorbed oxygen (${\mathrm{M}}^+$–${\mathrm{O}}^-_2$) and terminal oxygen (M=O), which can react with DPM to form carbonates that are easily converted to gaseous CO₂. Our results provide a practical plasma catalysis technology to obtain simultaneous removals of DPM and HCs at low temperatures.

Only a small amount of spectral information is collected because the collection solid angle of the optical fiber probe and lens is very limited when collecting spectral information. To overcome this limitation, this study presents a novel method for acquiring plasma spectral information from various spatial directions. A parabolic-shaped plasma spectral collection device (PSCD) is employed to effectively collect more spectral information into the spectrometer, thereby enhancing the overall spectral intensity. The research objects in this study were soil samples containing different concentrations of heavy metals Pb, Cr, and Cd. The results indicate that the PSCD significantly enhances the spectral signal, with an enhancement rate of up to 45%. Moreover, the signal-to-noise ratio also increases by as much as 36%. Simultaneously, when compared to the absence of a device, it is found that there is no significant variation in plasma temperature when the PSCD is utilized. This observation eliminates the impact of the spatial effect caused by the PSCD on the spectral intensity. Consequently, a concentration-spectral intensity relationship curve is established under the PSCD. The results revealed that the linear fitting R² for Pb, Cr, and Cd increased by 0.011, 0.001, and 0.054, respectively. Additionally, the limit of detection (LOD) decreased by 0.361 ppm, 0.901 ppm, and 0.602 ppm, respectively. These findings indicate that the spectral enhancement rate elevates with the increase in heavy metal concentration. Hence, the PSCD can effectively enhance the spectral intensity and reduce the detection limit of heavy metals in soil.

Recent work has validated a new method for estimating the grain size of microgranular materials in the range of tens to hundreds of micrometers using laser-induced breakdown spectroscopy (LIBS). In this situation, a piecewise univariate model must be constructed to estimate grain size due to the complex dependence of the plasma formation environment on grain size. In the present work, we tentatively construct a unified calibration model suitable for LIBS-based estimation of those grain sizes. Specifically, two unified multivariate calibration models are constructed based on back-propagation neural network (BPNN) algorithms using feature selection strategies with and without considering prior information. By detailed analysis of the performances of the two multivariate models, it was found that a unified calibration model can be successfully constructed based on BPNN algorithms for estimating the grain size in the range of tens to hundreds of micrometers. It was also found that the model constructed with a prior-guided feature selection strategy had better prediction performance. This study has practical significance in developing the technology for material analysis using LIBS, especially when the LIBS signal exhibits a complex dependence on the material parameter to be estimated.

Hydrogels are biomaterials with 3D networks of hydrophilic polymers. The generation of hydrogels is turning to the development of hydrogels with the help of enabling technologies. Plasma can tailor the hydrogels' properties through simultaneous physical and chemical actions, resulting in an emerging technology of plasma-activated hydrogels (PAH). PAH can be divided into functional PAH and biological tissue model PAH. This review systematically introduces the plasma sources, plasma etching polymer surface, and plasma cross-linking involved in the fabrication of PAH. The 'diffusion-drift-reaction model' is used to study the microscopic physicochemical interaction between plasma and biological tissue PAH models. Finally, the main achievements of PAH, including wound treatment, sterilization, 3D tumor model, etc, and their development trends are discussed.

Cold atmospheric plasmas (CAPs) have shown great applicability in agriculture. Many kinds of CAP sources have been studied in agricultural applications to promote plant growth and cure plant diseases. We briefly review the state-of-the-art stimulating effects of atmospheric-pressure dielectric-barrier-discharge (AP-DBD) plasmas, after the direct or indirect treatment of plants for growth promotion and disease control. We then discuss the special demands on the characteristics of the CAP sources for their applications in plant mutation breeding. An atmospheric and room temperature plasma (ARTP) jet generator with a large plasma irradiation area, a high enough concentration of chemically reactive species and a low gas temperature is designed for direct plant mutagenesis. Experimental measurements of the electrical, thermal and optical features of the ARTP generator are conducted. Then, an ARTP-P (ARTP for plant mutagenesis) mutation breeding machine is developed, and a typical case of plant mutation breeding by the ARTP-P mutation machine is presented using Coreopsis tinctoria Nutt. seeds. Physical and agricultural experiments show that the newly-developed ARTP-P mutation breeding machine with a large irradiation area can generate uniform CAP jets with high concentrations of chemically reactive species and mild gas temperatures, and have significant mutagenesis effects on the Coreopsis tinctoria Nutt. seeds. The ARTP-P mutation breeding machine may provide a platform for systematic studies on mutation mechanisms and results for various plant seeds under different operating conditions in future research.

This paper reviews the effects of resonant magnetic perturbation (RMP) on classical tearing modes (TMs) and neoclassical tearing modes (NTMs) from the theory, experimental discovery and numerical results with a focus on four major aspects: (i) mode mitigation, where the TM/NTM is totally suppressed or partly mitigated by the use of RMP; (ii) mode penetration, which means a linearly stable TM/NTM triggered by the externally applied RMP; (iii) mode locking, namely an existing rotating magnetic island braked and finally stopped by the RMP; (iv) mode unlocking, as the name suggests, it is the reverse of the mode locking process. The key mechanism and physical picture of above phenomena are revealed and summarized.

Atmospheric pressure cold plasma, with advantages such as high particle activity, no thermal damage, high efficiency and direct and friendly contact with human tissues, is considered to have great potential in biomedical applications. Therefore, 'plasma medicine' as a new interdiscipline has been developed in the past two decades. This review first briefly describes the development of typical plasma sources suitable for biomedical applications, and those with different discharge forms are simply compared, evaluated and summarized. Subsequently, measurement of the crucial gaseous reactive particles (e.g. OH and O) and their spatio-temporal distributions are introduced. Meanwhile, the generation and variation rules and the related critical macroscopic parameters of the plasma-induced aqueous reactive species are summarized. Finally, related studies in the last ten years on the mechanisms of the plasma-driven microbial inactivation and plasma-induced apoptosis of cancer cells are introduced. Moreover, some scientific problems that need to be urgently solved in the field of plasma medicine are also discussed. This review will provide useful guidance for future related research.

The influence of m/n = 2/1 (m and n are poloidal and toroidal mode numbers) tearing modes on plasma perpendicular flows and micro-fluctuations has been investigated in HL-2A neutral beam injection heated L-mode plasmas. It is found that the local perpendicular rotation velocity and turbulence energy are modulated by the alternation between the island X-point and O-point of the naturally rotating tearing modes. Cross-correlation analysis indicates that the modulation of density fluctuations by the tearing mode is not only limited to the island region, but also occurs in the edge region near the last closed flux surface. The turbulence exhibits distinct spectral characteristics inside and outside the island region. In addition, it is observed that the particle flux near the strike point is also significantly impacted by the tearing modes. The experimental evidence reveals that there are strong core-edge interactions between the core tearing modes and the edge transport.