Doppler Backscattering Research Articles

Study of the effects of the perpendicular velocity gradient on a Doppler backscattering system using a 2D full wave code

Doppler backscattering is a powerful diagnostic technique that is used to perform turbulence measurements in magnetic confinement devices. The effects of the perpendicular velocity distribution on a Doppler backscattering system are studied using a 2D full wave code, in which the finite-difference time-domain method is used to solve Maxwell's equations. A constant velocity and several different velocity gradients are considered, along with a wide range of density fluctuation levels. The numerical results show that the Doppler backscattering mainly occurs around the cut-off layer, and that it has a linear response to turbulence when the fluctuation level is low. However, at higher density fluctuation levels, nonlinear effects become important, and the influence of the area outside the cut-off layer becomes significant. The velocity gradient causes the backscattering spectrum to become non-Gaussian and asymmetrical; this effect has also been noticed in experiments in the Experimental Advanced Superconducting Tokamak.

Preliminary measurements of the edge magnetic field pitch from 2-D Doppler backscattering in MAST and NSTX-U (invited).

The Synthetic Aperture Microwave Imaging (SAMI) system is a novel diagnostic consisting of an array of 8 independently phased antennas. At any one time, SAMI operates at one of the 16 frequencies in the range 10-34.5 GHz. The imaging beam is steered in software post-shot to create a picture of the entire emission surface. In SAMI's active probing mode of operation, the plasma edge is illuminated with a monochromatic source and SAMI reconstructs an image of the Doppler back-scattered (DBS) signal. By assuming that density fluctuations are extended along magnetic field lines, and knowing that the strongest back-scattered signals are directed perpendicular to the density fluctuations, SAMI's 2-D DBS imaging capability can be used to measure the pitch of the edge magnetic field. In this paper, we present preliminary pitch angle measurements obtained by SAMI on the Mega Amp Spherical Tokamak (MAST) at Culham Centre for Fusion Energy and on the National Spherical Torus Experiment Upgrade at Princeton Plasma Physics Laboratory. The results demonstrate encouraging agreement between SAMI and other independent measurements.

Simultaneous measurement of magnetic and density fluctuations via cross-polarization scattering and Doppler backscattering on the DIII-D tokamak.

An upgraded cross-polarization scattering (CPS) system for the simultaneous measurement of internal magnetic fluctuations B̃ and density fluctuations ñ is presented. The system has eight radial quadrature channels acquired simultaneously with an eight-channel Doppler backscattering system (measures density fluctuations ñ and flows). 3-D ray tracing calculations based on the GENRAY ray tracing code are used to illustrate the scattering and geometric considerations involved in the CPS implementation on DIII-D. A unique quasi-optical design and IF electronics system allow direct comparison of B̃ and ñ during dynamic or transient plasma events (e.g., Edge Localized Modes or ELMs, L to H-mode transitions, etc.). The system design allows the interesting possibility of both magnetic-density (B̃-ñ) fluctuation and magnetic-temperature (B̃-T̃) fluctuation cross-phase measurements suitable for detailed tests of turbulence simulations.

Modifications to the synthetic aperture microwave imaging diagnostic.

The synthetic aperture microwave imaging diagnostic has been operating on the MAST experiment since 2011. It has provided the first 2D images of B-X-O mode conversion windows and showed the feasibility of conducting 2D Doppler back-scattering experiments. The diagnostic heavily relies on field programmable gate arrays to conduct its work. Recent successes and newly gained experience with the diagnostic have led us to modify it. The enhancements will enable pitch angle profile measurements, O and X mode separation, and the continuous acquisition of 2D DBS data. The diagnostic has also been installed on the NSTX-U and is acquiring data since May 2016.

Experimental investigation of geodesic acoustic modes on JET using Doppler backscattering

Geodesic acoustic modes (GAMs) have been investigated in JET ohmic discharges using mainly Doppler backscattering. Characteristics and scaling properties of the GAM are studied. Time and spatial resolved measurements of the perpendicular velocity indicate that GAMs are located in a narrow layer at the edge density gradient region with amplitude corresponding to about 50% of the mean local perpendicular velocity. GAMs on JET appear to be regulated by the turbulence drive rather than by their damping rate. It is also shown that the GAM amplitude is ~20% larger in deuterium than in hydrogen plasmas.

Quasi-optical design for systems to diagnose the electron temperature and density fluctuations on EAST.

A system to simultaneously diagnose the electron temperature and density fluctuations is proposed for Experimental Advanced Superconducting Tokamak device. This system includes a common quasi-optical antenna, a correlation electron cyclotron emission (CECE) system that is used to measure the electron temperature fluctuations and a Doppler backscattering (DBS) system that is used to measure the electron density fluctuations. The frequency range of the proposed CECE system is 108-120 GHz, and this corresponds to a radial coverage of normalized radius ((R - R0)/a, R0 = 1850 mm, a = 450 mm) from 0.2 to 0.67 for the plasma operation with a toroidal magnetic field of 2.26 T. This paper focuses on the design of the quasi-optical antenna and aims at optimizing the poloidal resolution for different frequency bands. An optimum result gives the beam radius for the CECE system of 13-15 mm and this corresponds to a wave number range of kθ < 2.4 cm-1. The beam radius is 20-30 mm for V band (50-75 GHz) and 15-20 mm for W band (75-110 GHz).

Role of density gradient driven trapped electron mode turbulence in the H-mode inner core with electron heating

A series of DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] low torque quiescent H-mode experiments show that density gradient driven trapped electron mode (DGTEM) turbulence dominates the inner core of H-mode plasmas during strong electron cyclotron heating (ECH). Adding 3.4 MW ECH doubles Te/Ti from 0.5 to 1.0, which halves the linear DGTEM critical density gradient, locally reducing density peaking, while transport in all channels displays extreme stiffness in the density gradient. This suggests that fusion α-heating may degrade inner core confinement in H-mode plasmas with moderate density peaking and low collisionality, with equal electron and ion temperatures, key conditions expected in burning plasmas. Gyrokinetic simulations using GYRO [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] (and GENE [Jenko et al., Phys. Plasmas 7, 1904 (2000)]) closely match not only particle, energy, and momentum fluxes but also density fluctuation spectra from Doppler backscattering (DBS), with and without ECH. Inner core DBS density fluctuations display discrete frequencies with adjacent toroidal mode numbers, which we identify as DGTEMs. GS2 [Dorland et al., Phys. Rev. Lett. 85, 5579 (2000)] predictions show the DGTEM can be suppressed, to avoid degradation with electron heating, by broadening the current density profile to attain q0&amp;gt;qmin&amp;gt;1.

Gyrokinetic simulation of driftwave instability in field-reversed configuration

Following the recent remarkable progress in magnetohydrodynamic (MHD) stability control in the C-2U advanced beam driven field-reversed configuration (FRC), turbulent transport has become one of the foremost obstacles on the path towards an FRC-based fusion reactor. Significant effort has been made to expand kinetic simulation capabilities in FRC magnetic geometry. The recently upgraded Gyrokinetic Toroidal Code (GTC) now accommodates realistic magnetic geometry from the C-2U experiment at Tri Alpha Energy, Inc. and is optimized to efficiently handle the FRC's magnetic field line orientation. Initial electrostatic GTC simulations find that ion-scale instabilities are linearly stable in the FRC core for realistic pressure gradient drives. Estimated instability thresholds from linear GTC simulations are qualitatively consistent with critical gradients determined from experimental Doppler backscattering fluctuation data, which also find ion scale modes to be depressed in the FRC core. Beyond GTC, A New Code (ANC) has been developed to accurately resolve the magnetic field separatrix and address the interaction between the core and scrape-off layer regions, which ultimately determines global plasma confinement in the FRC. The current status of ANC and future development targets are discussed.

Stationary Zonal Flows during the Formation of the Edge Transport Barrier in the JET Tokamak.

High spatial resolution Doppler backscattering measurements in JET have enabled new insights into the development of the edge Er. We observe fine-scale spatial structures in the edge Er well with a wave number krρi≈0.4-0.8, consistent with stationary zonal flows, the characteristics of which vary with density. The zonal flow amplitude and wavelength both decrease with local collisionality, such that the zonal flow E×B shear increases. Above the minimum of the L-H transition power threshold dependence on density, the zonal flows are present during L mode and disappear following the H-mode transition, while below the minimum they are reduced below measurable amplitude during L mode, before the L-H transition.

GAM observation in the TUMAN-3M tokamak

Results of an experimental study of geodesic acoustic modes (GAM) in the TUMAN-3M tokamak are reported. With Doppler backscattering (DBS) the basic properties of the GAM such as frequency, conditions for the GAM existence and the GAM radial location have been identified. The two-frequency Doppler reflectometer system was employed to reveal an interplay between low frequency sheared poloidal rotation, ambient turbulence level and the GAM intensity. Bicoherence analysis of the DBS data evidences the presence of a nonlinear interaction between the GAM and plasma turbulence.

2D Doppler backscattering using synthetic aperture microwave imaging of MAST edge plasmas

Doppler backscattering (DBS) is already established as a powerful diagnostic; its extension to 2D enables imaging of turbulence characteristics from an extended region of the cut-off surface. The Synthetic Aperture Microwave Imaging (SAMI) diagnostic has conducted proof-of-principle 2D DBS experiments of MAST edge plasma. SAMI actively probes the plasma edge using a wide (±40° vertical and horizontal) and tuneable (10–34.5 GHz) beam. The Doppler backscattered signal is digitised in vector form using an array of eight Vivaldi PCB antennas. This allows the receiving array to be focused in any direction within the field of view simultaneously to an angular range of 6–24° FWHM at 10–34.5 GHz. This capability is unique to SAMI and is a novel way of conducting DBS experiments. In this paper the feasibility of conducting 2D DBS experiments is explored. Initial observations of phenomena previously measured by conventional DBS experiments are presented; such as momentum injection from neutral beams and an abrupt change in power and turbulence velocity coinciding with the onset of H-mode. In addition, being able to carry out 2D DBS imaging allows a measurement of magnetic pitch angle to be made; preliminary results are presented. Capabilities gained through steering a beam using a phased array and the limitations of this technique are discussed.

Comparison of Three Methods in Extracting Coherent Modes from a Doppler Backscatter System**Supported by the National Natural Science Foundation of China under Grant Nos 10990210 and 11475173, and the National Magnetic Confinement Fusion Energy Program of China under Grant Nos 2013GB106002 and 2014GB109002.

We compare three different methods to extract coherent modes from Doppler backscattering (DBS), which are center of gravity (COG) of the complex amplitude spectrum, spectrum of DBS phase derivative (phase derivative method), and phase spectrum, respectively. These three methods are all feasible to extract coherent modes, for example, geodesic acoustic mode oscillation. However, there are still differences between dealing with high frequency modes (several hundred kHz) and low frequency modes (several kHz) hiding in DBS signal. There is a significant amount of power at low frequencies in the phase spectrum, which can be removed by using the phase derivative method and COG. High frequency modes are clearer by using the COG and the phase derivative method than the phase spectrum. The spectrum of DBS amplitude does not show the coherent modes detected by using COG, phase derivative method and phase spectrum. When two Doppler shifted peaks exist, coherent modes and their harmonics appear in the spectrum of DBS amplitude, which are introduced by the DBS phase.

Intermediate-k density and magnetic field fluctuations during inter-ELM pedestal evolution in MAST

Measurements of local density and magnetic field fluctuations near the pedestal top, conditionally averaged over the edge localized mode (ELM) cycle, have been made in Mega Amp Spherical Tokamak (MAST). A Doppler backscattering (DBS) system installed at MAST was used to measure intermediate-k density fluctuations at the top of the pedestal. A novel diagnostic technique combining DBS with cross-polarization scattering (CP-DBS) enabled magnetic field fluctuations to also be locally measured at similar wave numbers. Polarization isolation and other effects for CP-DBS are discussed. Both measurements were used in a series of high-β –4.5) MAST plasmas with large type-I ELMs with an period where microtearing modes (MTMs) had been predicted to be unstable in similar conditions (Dickinson 2012 Phys. Rev. Lett. 108 135002). The measured density fluctuation level increased by a factor of about 4 between 2 and 4 ms after the ELM, which was correlated with the recovery of the density profile while the temperature pedestal height continued to increase slowly. Magnetic field fluctuations showed different temporal behaviors, slowly increasing throughout the ELM cycle as the local β increased. Linear GS2 calculations show both MTM and electron temperature gradient (ETG) modes unstable at similar wave numbers as the measurements (although with more overlap between ETG wave numbers and diagnostic spectral resolution) at the top of the pedestal, along with kinetic ballooning modes are unstable lower in the pedestal (at larger wavelengths). The inferred ratio of fluctuation levels from experiment was . The comparable ratios from GS2 were for the MTM and for the ETG. Both the experimental wave number range and the fluctuation ratio are more similar to the linear characteristics of the ETG than the MTM. These results imply that intermediate-k fluctuations due to the ETG play a role in inter-ELM pedestal evolution.

Influence of lithium coating on the optics of Doppler backscatter system.

This paper presents the first investigation of the effect of lithium coating on the optics of Doppler backscattering. A liquid lithium limiter has been applied in the Experimental Advanced Superconducting Tokamak (EAST), and a Doppler backscattering has been installed in the EAST. A parabolic mirror and a flat mirror located in the vacuum vessel are polluted by lithium. An identical optical system of the Doppler backscattering is set up in laboratory. The power distributions of the emission beam after the two mirrors with and without lithium coating (cleaned before and after), are measured at three different distances under four incident frequencies. The results demonstrate that the influence of the lithium coating on the power distributions are very slight, and the Doppler backscattering can work normally under the dosage of lithium during the 2014 EAST campaign.

Multi-diagnostic approach to geodesic acoustic mode study

Multi-diagnostic approach developed for the GAM research in the spherical tokamak Globus M is described. Doppler backscattering (DBS) method as the tool for the GAM study, together with the diagnostics of plasma density and magnetic field GAM oscillations, were simultaneously used in experiments. The version of the DBS diagnostics with two cut-offs positioned at different poloidal angles of the minor cross-section was employed in Globus-M. For the GAM plasma density oscillation study, the Dα emission was observed at different angles to restore the spatial mode structure of the GAM plasma density oscillations. At the same time, the array of Mirnov coils was used for the GAM-like magnetic oscillation study, and that made it possible to restore the magnetic field perturbation spatial structure. The coherent and cross-bicoherence analyzes were employed to identify the interaction between the GAM velocity oscillation and plasma turbulent fluctuations.A shorter version of this contribution is due to be published in PoS at: 1st EPS conference on Plasma Diagnostics

Ray Tracing for Doppler Backscattering System in the Experimental Advanced Superconducting Tokamak**supported by National Natural Science Foundation of China (Nos. 10990211 and 11105146) and the ITER-CN Project, 973 Program of China (No. 2013GB106002)

The Doppler backscattering system has been widely used for turbulence measurements, and the microwave beam will be backscattered near the cut-off layer when the Brag condition is fulfilled. In tokamak, the ray-tracing code is used to obtain the radial position and perpendicular wave number of the scattering layer for turbulence velocity measurement and the WKB (Wentzel-Kramers-Brillouin) approximation should be satisfied for optical propagation. To calculate the backscattering location and wave number at the cut-off layer only, a single ray tracing in the cross section is enough, while for spatial and wave number resolution calculation, multiple rays reflecting the microwave beam size should be used. Considering the angle between the wave vector and the magnetic field, a three-dimension quasi-optical Gaussian ray tracing is sometimes needed.

Doppler backscattering for spherical tokamaks and measurement of high-k density fluctuation wavenumber spectrum in MAST

The high-k () wavenumber spectrum of density fluctuations has been measured for the first time in MAST (Lloyd et al 2003 Nucl. Fusion 43 1665). This was accomplished with the first implementation of Doppler backscattering (DBS) for core measurements in a spherical tokamak. DBS has become a well-established and versatile diagnostic technique for the measurement of intermediate- k (, and higher) density fluctuations and flows in magnetically confined fusion experiments. Previous implementations of DBS for core measurements have been in standard, large aspect ratio tokamaks. A novel implementation with two-dimensional (2D) steering was necessary to enable DBS measurements in MAST, where the large variation of the magnetic field pitch angle presents a challenge. We report on the scattering considerations and ray tracing calculations used to optimize the design and present data demonstrating measurement capabilities. Initial results confirm the applicability of the design and implementation approaches, showing the strong dependence of scattering alignment on the toroidal launch angle and demonstrating that DBS is sensitive to the local magnetic field pitch angle. We also present comparisons of DBS plasma velocity measurements with charge exchange recombination and beam emission spectroscopy measurements, which show reasonable agreement over most of the minor radius, but imply large poloidal flows approaching the magnetic axis in a discharge with an internal transport barrier. The 2D steering is shown to enable high-k measurements with DBS, at cm () for launch frequencies less than 75 GHz; this capability is used to measure the wavenumber spectrum of turbulence and we find for –11, which is similar to the expectation for the turbulent kinetic cascade of .

Overview of MAST results

The Mega Ampère Spherical Tokamak (MAST) programme is strongly focused on addressing key physics issues in preparation for operation of ITER as well as providing solutions for DEMO design choices. In this regard, MAST has provided key results in understanding and optimizing H-mode confinement, operating with smaller edge localized modes (ELMs), predicting and handling plasma exhaust and tailoring auxiliary current drive. In all cases, the high-resolution diagnostic capability on MAST is complemented by sophisticated numerical modelling to facilitate a deeper understanding. Mitigation of ELMs with resonant magnetic perturbations (RMPs) with toroidal mode number nRMP = 2, 3, 4, 6 has been demonstrated: at high and low collisionality; for the first ELM following the transition to high confinement operation; during the current ramp-up; and with rotating nRMP = 3 RMPs. nRMP = 4, 6 fields cause less rotation braking whilst the power to access H-mode is less with nRMP = 4 than nRMP = 3, 6. Refuelling with gas or pellets gives plasmas with mitigated ELMs and reduced peak heat flux at the same time as achieving good confinement. A synergy exists between pellet fuelling and RMPs, since mitigated ELMs remove fewer particles. Inter-ELM instabilities observed with Doppler backscattering are consistent with gyrokinetic simulations of micro-tearing modes in the pedestal. Meanwhile, ELM precursors have been strikingly observed with beam emission spectroscopy (BES) measurements. A scan in beta at the L–H transition shows that pedestal height scales strongly with core pressure. Gyro-Bohm normalized turbulent ion heat flux (as estimated from the BES data) is observed to decrease with increasing tilt of the turbulent eddies. Fast ion redistribution by energetic particle modes depends on density, and access to a quiescent domain with ‘classical’ fast ion transport is found above a critical density. Highly efficient electron Bernstein wave current drive (1 A W−1) has been achieved in solenoid-free start-up. A new proton detector has characterized escaping fusion products. Langmuir probes and a high-speed camera suggest filaments play a role in particle transport in the private flux region whilst coherence imaging has measured scrape-off layer (SOL) flows. BOUT++ simulations show that fluxes due to filaments are strongly dependent on resistivity and magnetic geometry of the SOL, with higher radial fluxes at higher resistivity. Finally, MAST Upgrade is due to begin operation in 2016 to support ITER preparation and importantly to operate with a Super-X divertor to test extended leg concepts for particle and power exhaust.

Comprehensive comparisons of geodesic acoustic mode characteristics and dynamics between Tore Supra experiments and gyrokinetic simulations

In a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time.

Chaotic density fluctuations in L-mode plasmas of the DIII-D tokamak

Analysis of the time series obtained with the Doppler backscattering system (Hillsheim et al 2009 Rev. Sci. Instrum. 80 0835070) in the DIII-D tokamak (Luxon 2005 Fusion Sci. Technol. 48 828) shows that intermediate wave number plasma density fluctuations in low confinement (L-mode) tokamak plasmas are chaotic. The supporting evidence is based on the shape of the power spectrum; the location of the signal in the complexity-entropy plane (C-H plane) (Rosso et al 2007 Phys. Rev. Lett. 99 154102); and the population of the corresponding Bandt–Pompe (Bandt and Pompe 2002 Phys. Rev. Lett. 88 174102) probability distributions.