Lower Hybrid Range Research Articles

Investigation of electromagnetic fluctuations in a magnetically screened high beta plasma

We report low-frequency electromagnetic turbulence in large-volume plasma devices in a region that receives magnetically screened plasma. The magnetic screen is produced by the activation of a large solenoid that generates a strong magnetic field ( BEEF ) transverse to the background axial magnetic field ( Bz ). The magnetic field strength variation in BEEF controls the plasma diffusion and is responsible for parametric profile modifications in large volume plasma devices. The turbulence is characterized by the observed features of an electromagnetic mode having frequency ordering in the lower hybrid range (f ci < f < f LH). The free energy source for the excited turbulence is believed to be the prevalent finite radial pressure gradient in the system. The excited electromagnetic mode shows typical scales for the wave number, k⊥ ∼ 0.13 cm−1 and frequency f ∼ 6 kHz and propagates with phase velocity comparable to ion acoustic velocity, C s in electron diamagnetic drift direction. It is believed that the low-frequency pressure gradient-driven electromagnetic mode exhibits coupling with the lower hybrid or whistler wave.

Resonance cones in cold plasma: Origin, singularities, and power flow

In magnetized tenuous plasma, typical at the plasma edge of fusion devices, a nearly electrostatic wave mode with relatively enhanced electric field can propagate along a specific angle with the magnetic field. For this characteristic, it is known as a “resonance cone.” For instance, these waves can be excited by radio frequency antennas in the ion-cyclotron and lower-hybrid range of frequencies. We consider the resonance cones emitted by idealized spatially extended sources. In 2D, we use a novel geometric construction which generalizes the d'Alembert solution to curved boundaries/moving sources, and show, for the first time, that singular electric fields arise under these conditions, thereby bringing the resonance cones in line with the other resonances of the cold plasma theory. Still in 2D, we give an expression for the amount of power radiated by resonance cones in terms of surface quantities on the source, which is finite despite the singular electric field. We generalize the conclusions regarding the presence and location of singular electric fields to the 3D electromagnetic case.

Particle simulations on propagation and resonance of lower hybrid wave launched by phased array antenna in linear devices

In this work, we performed first-principles electromagnetic-kinetic simulations to study a phased antenna array and its interaction with deuterium plasmas within the lower hybrid range of frequency. We first gave wave accessibility and resonance results, which agree well with theoretical prediction. In addition, we further investigated the antenna power spectrum with different antenna phases in the presence of the plasma and compared it with that in a vacuum, which directly indicates wave coupling and plasma absorption. Furthermore, for the case with zero phasing difference, our simulation results show that, albeit the launch is away from the accessibility region, tunneling effect and mode conversion occurred, which enhanced coupling and absorption. Moreover, consistent interactions between the injected wave and the plasma concerning various antenna phase differences are shown. We presented the inchoate response of the plasma in terms of the launching directions. Our results could be favorable for the engineering design of wave heating experiments with a tunable phased antenna array in linear devices, such as simple magnetic mirrors or tandem mirrors.

The high-power helicon program at DIII-D: gearing up for first experiments

Helicon current drive (CD), also called fast wave CD in the lower hybrid range of frequencies, has long been regarded as a promising CD tool for reactor grade plasmas. A newly installed MW-level system at DIII-D will be the first test of this technology in reactor-relevant plasmas, in the sense that full single-pass absorption is expected. A 30-module traveling wave antenna has been installed and optimized in-vessel in early 2020. The linear electromagnetic characteristics of the unloaded module array have been extensively tested both on the bench and in the vessel at instrumentation power levels. Excellent performance has been achieved, ∼2% reflected power and ∼1.5% dissipated power per module in air, in a 10 MHz band around 476 MHz. Stripline feeds on both ends of the antenna allow either co or counter CD. The installation of a 1.2 MW klystron and associated high-power electronics was completed in Fall 2020. Commissioning of the antenna is ongoing. An important goal of this experiment is to validate the helicon CD physics basis using an extensive set of new and upgraded diagnostics.

Investigations into growth of whistlers with energy of energetic electrons

Runaway electrons have great significance in fusion plasmas as they are considered to be responsible for the excitation of instabilities and for damage to the first wall components. The threshold for the generation of the runaway electron population in tokamaks strongly depends on the toroidal magnetic field and bulk electron temperature. The presence of these runaway electrons with a critical density and magnetic field value are correlated with the magnetosonic whistler activity. The Rosenbluth condition applies to the destabilization of upper hybrid modes, that remains mostly stable in large-scale tokamaks like Joint European Tokamakand International Thermonuclear Experimental reactor. On the other hand, whistlers with a frequency in the lower hybrid range are more likely to be destabilized in a lower magnetic field and temperature, along with the whistlers’ excitation threshold. The process of destabilization of whistlers by energetic electrons thus has a finite overlap between fusion and basic plasma devices, and the outcome from the latter can be scalable to fusion machines. We report that in the linear large volume plasma device, the destabilization of quasi-longitudinal whistlers, driven by the reflected particles from the mirror, takes place in the bulk plasma region i.e. away from the mirror in the energetic belt region, and follows the frequency ordering of . Our explanation involves an alternate limit of quasi-linear analysis, commonly employed for recovering the instability threshold of runaways scattering off whistlers in fusion plasmas, and the growth rate due to the reflected particle is proportional to the inverse square root of electron temperature, where the reflected particle fraction is (Sharma and Vlahos 1984 Astrophys. J. 280 405). The measurements required to obtain the energy scaling of runaway electrons with whistler instability are critical for large devices but remain scarce. We have emphasized the analytic correlation between the two regimes in these investigations.

Density dependence of ion cyclotron emission from deuterium plasmas in the large helical device

Ion cyclotron emission (ICE) driven by perpendicular neutral beam-injected (NBI) deuterons, together with the distinctive ICE driven by tangential NBI, have been observed from heliotron–stellarator plasmas in the large helical device (LHD). Radio frequency radiation in the lower hybrid range has also been observed Saito K. et al (2018 Plasma Fusion Res. 13 3402043), with frequency dependent on plasma density. Here we focus on recent measurements of ICE from deuterium plasmas in LHD, which show substantial variation in spectral character, between otherwise similar plasmas that have different local density in the emitting region. We analyse this variation by means of first principles simulations, carried out using a particle-in-cell (PIC) kinetic approach. We show, first, that this ICE is driven by perpendicular NBI deuterons, freshly ionised near their injection point in the outer midplane edge of LHD. We find that these NBI deuterons undergo collective sub-Alfvénic relaxation, which we follow deep into the nonlinear phase of the magnetoacoustic cyclotron instability (MCI). The frequency and wavenumber dependence of the saturated amplitudes of the excited fields determine our simulated ICE spectra, and these spectra are obtained for different local densities corresponding to the different LHD ICE-emitting plasmas. The variation with density of the spectral character of the simulated ICE corresponds well with that of the observed ICE from LHD. These results from heliotron–stellarator plasmas complement recent studies of density-dependent ICE from tokamak plasmas in KSTAR Thatipamula S.G. et al (2016 Plasma Phys. Control. Fusion 58 065003); Chapman B. et al (2017 Nucl. Fusion 57 124004), where the spectra vary on sub-microsecond timescales after an ELM crash. Taken together, these results confirm the strongly spatially localised character of ICE physics, and reinforce the potential of ICE as a diagnostic of energetic ion populations and of the ambient plasma.

Avoidance of impurity-induced current quench using lower hybrid current drive

This work reports observations of a tokamak plasma that experienced a thermal quench due to a large, transient high-Z influx but avoided a current quench. This is argued to be caused by the presence of lower hybrid range of frequency (LHRF) waves that sustain a non-thermal, current-carrying electron population. In Alcator C-Mod L-mode plasmas at kA, m−3, nearly all of the current can be sustained non-inductively by injecting 700 kW of LHRF power at 4.6 GHz and . A sudden influx of a large amount of tungsten, , triggers a cooling wave that propagates at 2–3 m s−1 all the way into the core, dropping on-axis Te from 3 keV to temperatures less than measurement floor of 50 eV. An off-axis reheat begins after 100 ms, but Te profiles remain hollow for 300–350 ms after the injection. Throughout this temperature evolution, the plasma density, current and shape remain unchanged to within 10%. Following the expulsion of the tungsten, the plasma returns to its baseline conditions and the plasma ends as planned with a controlled current ramp-down. Energy balance analysis shows the LHRF power continues to be absorbed in the plasma after the thermal quench, as a significant fraction of it is needed to be consistent with radiated power measurements. Examination of current relaxation time, , and fast-electron slowing down time, , indicate the LHRF must contribute to driving current, despite the low temperatures, as the current remains nominally stationary despite ms and ms for relativistic electrons. These measurements represent an important existence proof of a possible technique for avoidance of disruptions caused by sudden, unplanned influx of impurities in the form of dust or flakes of high-Z wall material. Implications and suggestions for future experimental and modeling and simulation work are summarized.

Lower hybrid range cold magnetized plasma coupling in ANSYS HFSS

The coupling between cold magnetized plasmas with lower hybrid resonance frequency antennas is generally addressed using specifically developed codes. These antenna coupling codes often approximate the plasma to a surface impedance described by a 1D half infinite models and antennas either with simplified 2D dimensions or from 3D CAD models conversion. Such approaches add an additional step of model approximation/conversion which is not convenient to rapidly assess impact of geometry changes on coupling performances. In this work, we assess the ability of using ANSYS HFSS to describe the usual range of experimental magnetized cold plasma inhomogeneous parameters facing LHRF antennas, away from resonances, in order to guide the RF designer during the design phase of an antenna. The coupling calculations of Lower Hybrid Resonance Frequency antennas are performed and successfully benchmarked with the fast coupling codes ALOHA on inhomogeneous cold magnetized plasmas. Good agreements are obtained when the boundary conditions in the full-wave modelling are properly handled. Practical advices and limitations are given in order to define correctly these boundary conditions and insure correct results.

Tests of advanced RF off-axis current drive techniques on DIII-D

The establishment of reactor-relevant radiofrequency heating and current drive techniques is a focus of work on DIII-D in the next five-year period. This paper gives an overview of the planned experimental work in the areas of (1) nearly vertically launched ECCD, (2) ‘helicon’ (whistlers or fast waves in the lower hybrid range of frequencies) current drive, and (3) high-field-side-launch (HFS) lower hybrid (slow wave) current drive. Each of these techniques addresses the need for efficient off-axis current drive for a steady-state tokamak reactor to supplement the bootstrap current and to provide current profile control, and each will be experimentally assessed at a coupled power level of ~1 MW on DIII-D in the next few years.

High field side lower hybrid wave launch for steady state plasma sustainment

This paper explores how the use of waves in the lower hybrid range of frequencies (LHRF) alleviates many unresolved current drive issues when launched from the high field side (HFS) of a tokamak. Efficient, robust, reliable, steady-state current drive is a requirement for a viable tokamak fusion reactor. Launching LHRF waves from the HFS opens the ‘accessibility window’, allowing for use of a lower parallel refractive index, , for a given plasma density and on-axis magnetic field. The lower waves damp in a region of higher temperature, resulting in increased current drive efficiency and a damping profile shifted closer to the mid-radius as compared to low field side launch. Moving the LHRF antenna to the quiescent scrape off layer on the HFS reduces the likelihood of detrimental wave scattering from density fluctuations, as well as ameliorating the plasma-material-interaction and wave coupling challenges. Additional benefits with respect to fast particle and neutron fluxes can be realized by moving to the HFS. An engineering assessment of HFS launch shows that it is feasible to implement an HFS LHRF system on existing tokamaks, with fewer hurdles anticipated for future tokamaks built with HFS antenna systems integrated in the design from the start.

Experiments on helicons in DIII-D—investigation of the physics of a reactor-relevant non-inductive current drive technology**Paper based on contribution EX/P3-22 presented at IAEA FEC, Kyoto, October 2016.

Experiments have begun in DIII-D to evaluate a toroidally-directed spectrum of helicon waves (also known as ‘very high harmonic fast waves’, ‘fast waves in the lower hybrid range of frequencies’, or ‘whistlers’) that can generate efficient non-inductive current drive by Landau absorption in a reactor plasma. Modeling has shown (Prater et al 2014 Nucl. Fusion 54 083024) that non-inductive current drive at mid-radius (ρ ~ 0.5) should be achievable in DIII-D with fast waves at 0.5 GHz in high-beta conditions with an efficiency twice that of other non-inductive current drive tools currently available on DIII-D (neutral beams and electron cyclotron current drive). An innovative traveling wave antenna (TWA) of the ‘comb-line’ type with 12 radiating modules has been constructed, installed in DIII-D, and has been tested at very low power (<0.4 kW) to evaluate the antenna coupling in the linear regime, and to test technological aspects of such structures in the tokamak environment. Results show favorably strong antenna/plasma coupling in several plasma regimes, most significantly in the high-beta ELMing H-mode case with calculated full first-pass absorption of the helicon waves in the desired mid-radius region. No evidence of significant excitation of the undesirable slow wave was found when the total magnetic field adjacent to the structure was approximately aligned with the tilted Faraday screen of the antenna. The implications of these results for a high-power experiment currently under construction are discussed.

Heating and current drive actuators study for FNSF in the ion cyclotron and lower hybrid range of frequency

This paper addresses considerations for the radio frequency (RF) heating and current drive actuators on the Fusion Nuclear Science Facility (FNSF), focusing on the Lower Hybrid Range of Frequency (LHRF) and Ion Cyclotron Range of Frequency (ICRF). FNSF will require long (∼months), steady-state pulses to fulfill its mission of high neutron fluence operation for the qualification of fusion nuclear technologies prototypical for a power plant. RF actuators are uniquely capable of sustaining the plasma without requiring large line-of-sight penetrations through the breeding blanket, which mitigates the impact on tritium breeding and neutron shielding. The technology for generation and transmission of RF waves in the LHRF and ICRF are fully mature, however the lifetime of plasma-facing antenna components in the reactor environment remains a significant technology gap. Moving the location of RF actuators from the low field side (LFS) to the high field side (HFS) of the tokamak reduces fluxes of turbulent plasma and energetic particles to the antenna. Compact HFS antenna designs are able to fit within the blanket sectors with minimal impact on overall dimensions. Material choices for antenna structures represent a significant challenge. Existing materials are unable to simultaneously meet the requirements for high conductivity, high strength, high temperature, and low activation. Pathways forward with respect to materials development are identified.

Characterization of the mutual influence of Ion Cyclotron and Lower Hybrid Range of frequencies systems on EAST

Waves in the Ion Cyclotron (ICRF) and Lower Hybrid (LH) Range of Frequencies are efficient techniques respectively to heat the plasma and drive current. Main difficulties come from a trade-off between good RF coupling and acceptable level of impurities release. The mutual influence of both systems makes such equilibrium often hard to reach [1]. In order to investigate those interactions based on Scrape-Off Layer (SOL) plasma parameters, a new reciprocating probe was designed allying a three tips Langmuir probe with an emissive wire. The emissive filament provides a precise measure of plasma potential [2], which can be used to calibrate Langmuir probe's results. This paper reports on experimental results obtained on EAST, where there are two ICRF antennas and two LH launchers. Among others diagnostics, the new reciprocating probe enabled to evidence the deleterious influence of ICRF power on LHWs coupling in L-mode plasmas. In areas connected with an active ICRF antenna, SOL potentials increase while densities tend to decrease, respectively enhancing impurities release and deteriorating LHWs coupling. This phenomenon has mostly been attributed to RF sheath; the one that forms on top of Plasma Facing Components (PFCs) and causes ExB density convections [3]. From those experiments it seems ICRF has a strong influence on magnetically connected areas, both in the near field – influencing ICRF waves coupling – and in farther locations such as in front of LH grills. Moreover, influence of ICRF on LH system was observed both in L and H modes. Those results are consistent with RF sheath rectification process. Concerning the influence of LHWs on ICRF coupling, nothing was observed in L-mode. Besides during H-mode experiments, LHWs have been identified as having a mitigating effect on ELMs [4], which on average lowers the pedestal, increasing edge densities to the profit of ICRF waves coupling.

Current drive with combined electron cyclotron wave and high harmonic fast wave in tokamak plasmas

The current driven by combined electron cyclotron wave (ECW) and high harmonic fast wave is investigated using the GENRAY/CQL3D package. It is shown that no significant synergetic current is found in a range of cases with a combined ECW and fast wave (FW). This result is consistent with a previous study [Harvey et al., in Proceedings of IAEA TCM on Fast Wave Current Drive in Reactor Scale Tokamaks (Synergy and Complimentarily with LHCD and ECRH), Arles, France, IAEA, Vienna, 1991]. However, a positive synergy effect does appear with the FW in the lower hybrid range of frequencies. This positive synergy effect can be explained using a picture of the electron distribution function induced by the ECW and a very high harmonic fast wave (helicon). The dependence of the synergy effect on the radial position of the power deposition, the wave power, the wave frequency, and the parallel refractive index is also analyzed, both numerically and physically.

Heating and current drive by fast wave in lower hybrid range of frequency on Versatile Experiment Spherical Torus

An efficient heating and current drive scheme in central or off-axis region is required to realize steady state operation of tokamak fusion reactor. And the fast wave in lower hybrid resonance range of frequency could be a candidate for such an efficient scheme in high density and high temperature plasmas. Its propagation and absorption characteristics including current drive and coupling efficiency are analyzed for Versatile Experiment Spherical Torus and it is shown that it is possible to drive current with considerable current drive efficiency in central region. The RF system for the fast wave experiment including klystron, transmission systems, inter-digital antenna, and RF diagnostics are given as well in this paper.

Particle simulations of mode conversion between slow mode and fast mode in lower hybrid range of frequencies

The propagation and mode conversion of lower hybrid waves in an inhomogeneous plasma are investigated by using the nonlinear δf algorithm in a two-dimensional particle-in-cell simulation code based on the gyrokinetic electron and fully kinetic ion (GeFi) scheme [Lin et al., Plasma Phys. Controlled Fusion 47, 657 (2005)]. The characteristics of the simulated waves, such as wavelength, frequency, phase, and group velocities, agree well with the linear theoretical analysis. It is shown that a significant reflection component emerges in the conversion process between the slow mode and the fast mode when the scale length of the density variation is comparable to the local wavelength. The dependences of the reflection coefficient on the scale length of the density variation are compared with the results based on the linear full wave model for cold plasmas. It is indicated that the mode conversion for the waves with a frequency of 2.45 GHz (ω ∼ 3ωLH, where ωLH represents the lower hybrid resonance) and within Tokamak relevant amplitudes can be well described in the linear scheme. As the frequency decreases, the modification due to the nonlinear term becomes important. For the low-frequency waves (ω ∼ 1.3ωLH), the generations of the high harmonic modes and sidebands through nonlinear mode-mode coupling provide new power channels and thus could reduce the reflection significantly.

Whistlers, helicons, and lower hybrid waves: The physics of radio frequency wave propagation and absorption for current drive via Landau damping

This introductory-level tutorial article describes the application of plasma waves in the lower hybrid range of frequencies (LHRF) for current drive in tokamaks. Wave damping mechanisms in a nearly collisionless hot magnetized plasma are briefly described, and the connections between the properties of the damping mechanisms and the optimal choices of wave properties (mode, frequency, wavelength) are explored. The two wave modes available for current drive in the LHRF are described and compared. The terms applied to these waves in different applications of plasma physics are elucidated. The character of the ray paths of these waves in the LHRF is illustrated in slab and toroidal geometries. Applications of these ideas to experiments in the DIII-D tokamak are discussed.

Improved confinement in high-density H-modes via modification of the plasma boundary with lower hybrid wavesa)

Injecting Lower Hybrid Range of Frequency (LHRF) waves into Alcator C-Mod's high-density H-mode plasmas has led to enhanced global energy confinement by increasing pedestal temperature and pressure gradients, decreasing the separatrix density, modifying the pedestal radial electric field and rotation, and decreasing edge turbulence. These experiments indicate that edge LHRF can be used as an actuator to increase energy confinement via modification of boundary quantities. H98-factor increases of up to ∼35% (e.g., H98 from 0.75 to 1.0) are seen when moderate amounts of LH power (PLH/Ptot ∼ 0.15) are applied to H-modes of densities n¯e ∼ 3 × 1020 m−3, corresponding to values ∼0.5 of the Greenwald density. However, the magnitude of the improvement is reduced if the confinement quality of the target H-mode plasma is already good (i.e., H98target ∼ 1). Ray-tracing modeling and accessibility calculations for the LH waves indicate that they do not penetrate to the core. The LHRF power appears to be deposited in p...

Lower hybrid drift instability at a dipolarization front

AbstractWe present observations of a reconnection jet front detected by the Cluster satellites in the magnetotail of Earth, which are commonly referred to as dipolarization fronts. We investigate in detail electric field structures observed at the front which have frequency in the lower hybrid range and amplitudes reaching 40 mV/m. We determine the frequency and phase velocity of these structures in the reference frame of the front and identify them as a manifestation of the lower hybrid drift instability (LHDI) excited at the sharp density gradient at the front. The LHDI is observed in the nonlinear stage of its evolution as the electrostatic potential of the structures is comparable to ∼ 10% of the electron temperature. The front appears to be a coherent structure on ion and MHD scales, suggesting existence of a dynamic equilibrium between excitation of the LHDI and recovery of the steep density gradient at the front.

One dimensional full wave analysis of slow-to-fast mode conversion in lower hybrid frequencies

The linear conversion from the slow wave to the fast wave in the lower hybrid range of frequencies is analyzed numerically by using the set of field equations describing waves in a cold plane-stratified plasma. The equations are solved as a two-point boundary value problem, where the polarizations of each mode are set consistently in the boundary conditions. The scattering coefficients and the field patterns are obtained for various density profiles. It is shown that, for large density scale length, the results agree well with the traditional cognitions. In contrast, the reflected component and the probable transmitted-converted component from the conversion region, which are neglected in the usual calculations, become significant when the scale length is smaller than the wavelength of the mode. The inclusion of these new components will improve the accuracy of the simulated propagation and deposition for the injected rf power when the conversion process is involved within a sharp-varying density profile. Meanwhile, the accessibility of the incident slow wave for the low frequency case is also affected by the scale length of the density profile.