-
Machine learning characterisation of Alfvénic and sub-Alfvénic chirping and correlation with fast ion loss at NSTX
Authors:
Benjamin J. Q. Woods,
Vinícius N. Duarte,
Eric D. Fredrickson,
Nikolai N. Gorelenkov,
Mario Podestà,
Roddy G. L. Vann
Abstract:
Abrupt large events in the Alfvénic and sub-Alfvénic frequency bands in tokamaks are typically correlated with increased fast ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behaviour of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. Analysis allows for comparison between different mode char…
▽ More
Abrupt large events in the Alfvénic and sub-Alfvénic frequency bands in tokamaks are typically correlated with increased fast ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behaviour of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. Analysis allows for comparison between different mode character (such as quiescent, fixed-frequency, chirping, avalanching) and plasma parameters obtained from the TRANSP code such as the ratio of the neutral beam injection (NBI) velocity and the Alfvén velocity ($v_{\textrm{inj.}}/v_{\textrm{A}}$), the $q$-profile, and the ratio of the neutral beam beta and the total plasma beta ($β_{\textrm{beam},i}/β$). In agreement with previous work by Fredrickson et al. [Nucl. Fusion 2014, 54 093007], we find correlation between $β_{\textrm{beam},i}$ and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to non-quiescent behaviour for magnetic fluctuations in the 50 - 200 kHz frequency band is observed along the boundary $v_{\varphi} \lessapprox \frac{1}{4}(v_{\textrm{inj.}} - 3v_{\textrm{A}})$ where $v_{\varphi}$ is the rotation velocity.
△ Less
Submitted 9 January, 2020; v1 submitted 12 March, 2019;
originally announced March 2019.
-
The deteriorating effect of plasma density fluctuations on microwave beam quality
Authors:
Alf Köhn,
Max E. Austin,
Michael W. Brookman,
Kenneth W. Gentle,
Lorenzo Guidi,
Eberhard Holzhauer,
Rob J. La Haye,
Jarrod B. Leddy,
Omar Maj,
Craig C. Petty,
Emanuele Poli,
Terry L. Rhodes,
Antti Snicker,
Matthew B. Thomas,
Roddy G. L. Vann,
Hannes Weber
Abstract:
Turbulent plasma edge density fluctuations can broaden a traversing microwave beam degrading its quality. This can be a problem for scenarios relying on a high spatial localization of the deposition of injected microwave power, like controlling MHD instabilities. Here we present numerical estimations of the scattering of a microwave beam by density fluctuations over a large parameter range, includ…
▽ More
Turbulent plasma edge density fluctuations can broaden a traversing microwave beam degrading its quality. This can be a problem for scenarios relying on a high spatial localization of the deposition of injected microwave power, like controlling MHD instabilities. Here we present numerical estimations of the scattering of a microwave beam by density fluctuations over a large parameter range, including extrapolations to ITER. Two codes are used, the full-wave code IPF-FDMC and the wave kinetic equation solver WKBeam. A successful comparison between beam broadening obtained from DIII-D experiments and corresponding full-wave simulations is shown.
△ Less
Submitted 24 May, 2018;
originally announced May 2018.
-
Positivity-preserving scheme for two-dimensional advection-diffusion equations including mixed derivatives
Authors:
Erasmus J. du Toit,
Martin R. O'Brien,
Roddy G. L. Vann
Abstract:
In this work, we propose a positivity-preserving scheme for solving two-dimensional advection-diffusion equations including mixed derivative terms, in order to improve the accuracy of lower-order methods. The solution to these equations, in the absence of mixed derivatives, has been studied in detail, while positivity-preserving solutions to mixed derivative terms have received much less attention…
▽ More
In this work, we propose a positivity-preserving scheme for solving two-dimensional advection-diffusion equations including mixed derivative terms, in order to improve the accuracy of lower-order methods. The solution to these equations, in the absence of mixed derivatives, has been studied in detail, while positivity-preserving solutions to mixed derivative terms have received much less attention. A two-dimensional diffusion equation, for which the analytical solution is known, is solved numerically to show the applicability of the scheme. It is further applied to the Fokker-Planck collision operator in two-dimensional cylindrical coordinates under the assumption of local thermal equilibrium. For a thermal equilibration problem, it is shown that the scheme conserves particle number and energy, while the preservation of positivity is ensured and the steady-state solution is the Maxwellian distribution.
△ Less
Submitted 3 February, 2018;
originally announced February 2018.
-
Stochastic effects on phase-space holes and clumps in kinetic systems near marginal stability
Authors:
Benjamin J. Q. Woods,
Vinicius N. Duarte,
Anthony P. De-Gol,
Nikolai N. Gorelenkov,
Roddy G. L. Vann
Abstract:
The creation and subsequent evolution of marginally-unstable modes have been observed in a wide range of fusion devices. This behaviour has been successfully explained, for a single frequency shifting mode, in terms of phase-space structures known as a `hole' and `clump'.
Here, we introduce stochasticity into a 1D kinetic model, affecting the formation and evolution of resonant modes in the syst…
▽ More
The creation and subsequent evolution of marginally-unstable modes have been observed in a wide range of fusion devices. This behaviour has been successfully explained, for a single frequency shifting mode, in terms of phase-space structures known as a `hole' and `clump'.
Here, we introduce stochasticity into a 1D kinetic model, affecting the formation and evolution of resonant modes in the system. We find that noise in the fast particle distribution or electric field leads to a shift in the asymptotic behaviour of a chirping resonant mode; this noise heuristically maps onto microturbulence via canonical toroidal momentum scattering, affecting hole and clump formation. The profile of a single bursting event in mode amplitude is shown to be stochastic, with small changes in initial conditions affecting the lifetime of a hole and clump. As an extension to the work of Lang and Fu, we find that an intermediate regime exists where noise serves to decrease the effective collisionality, where microturbulence works against pitch-angle scattering.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
Resolving ECRH deposition broadening due to edge turbulence in DIII-D by heat deposition measurement
Authors:
Michael W Brookman,
Matthew B Thomas,
Jarrod Leddy,
C Craig Petty,
Robert J La Haye,
K Barada,
Terry L Rhodes,
Zheng Yan,
Max E Austin,
Roddy G L Vann
Abstract:
Interaction between microwave power, used for local heating and mode control, and density fluctuations can produce a broadening of the injected beam, as confirmed in experiment and simulation. Increased power deposition width could impact suppression of tearing mode structures on ITER. This work discusses the experimental portion of an effort to understand scattering of injected microwaves by turb…
▽ More
Interaction between microwave power, used for local heating and mode control, and density fluctuations can produce a broadening of the injected beam, as confirmed in experiment and simulation. Increased power deposition width could impact suppression of tearing mode structures on ITER. This work discusses the experimental portion of an effort to understand scattering of injected microwaves by turbulence on the DIII-D tokamak. The corresponding theoretical modeling work can be found in M.B. Thomas et. al.: Submitted to Nuclear Fusion (2017)[Author Note - this paper to be published in same journal]. In a set of perturbative heat transport experiments, tokamak edge millimeter-scale fluctuation levels and microwave heat deposition are measured simultaneously. Beam broadening is separated from heat transport through fitting of modulated fluxes. Electron temperature measurements from a 500 kHz, 48-channel ECE radiometer are Fourier analyzed and used to calculate a deposition-dependent flux. Consistency of this flux with a transport model is evaluated. A diffusive and convective transport solution is linearized and compared with energy conservation-derived fluxes. Comparison between these two forms of heat flux is used to evaluate the quality of ECRF deposition profiles, and fitting finds a significant broadening of 1D equilibrium ray tracing calculations from the benchmarked TORAY-GA ray tracing code is needed. The physical basis, cross-validation, and application of the heat flux method is presented. The method is applied to a range of DIII-D discharges and finds a broadening factor of the deposition profile width which scales linearly with edge density fluctuation level. These experimental results are found to be consistent with the full-wave beam broadening measured by the 3D full wave simulations in the same discharges.
△ Less
Submitted 10 October, 2017;
originally announced October 2017.
-
Resolving ECRH deposition broadening due to edge turbulence in DIII-D by 3D full-wave simulations
Authors:
Matthew B Thomas,
Michael W Brookman,
Max E Austin,
Alf Köhn,
Robert J La Haye,
Jarrod B Leddy,
Roddy G L Vann,
Zheng Yan
Abstract:
Edge plasma density fluctuations are shown to have a significant effect on the electron cyclotron resonance heating (ECRH) beam in the DIII-D tokamak. Experimental measurements of the ECRH deposition profile have been taken in three operating scenarios: L-mode, H-mode and negative triangularity. Each scenario corresponds to distinct turbulence characteristics in the edge region through which the b…
▽ More
Edge plasma density fluctuations are shown to have a significant effect on the electron cyclotron resonance heating (ECRH) beam in the DIII-D tokamak. Experimental measurements of the ECRH deposition profile have been taken in three operating scenarios: L-mode, H-mode and negative triangularity. Each scenario corresponds to distinct turbulence characteristics in the edge region through which the beam must propagate. The measured ECRH deposition profile is significantly broadened by comparison to the profile predicted by the ray tracing code TORAY-GA and has been shown to scale with the severity of edge turbulence. Conventional ray tracing does not include the effects of turbulence and therefore a 3D full-wave cold plasma finite difference time domain code EMIT-3D is presented and used for the simulations. The turbulence is generated through the Hermes model in the BOUT++ framework which takes as input the measured time averaged electron density, temperature and magnetic field profiles for the specific shot in question. The simulated turbulence is constrained to match the experimentally measured (by use of the BES and DBS systems) correlation length and normalised fluctuation levels. The predictions of the beam broadening from the simulations are found to agree very well with the experimentally-observed broadening in all cases: L-mode, H-mode and negative triangularity. Due to the large gradients within the H-mode edge, the resolution uncertainty and error in the measurement from Thomson scattering and BES diagnostics result in a spread in the simulated turbulence amplitude. In light of this a parameter scan through the range in experimental diagnostic measurement uncertainty has been conducted to explore the impact on beam broadening predictions.
△ Less
Submitted 9 October, 2017;
originally announced October 2017.
-
Kinetic study of plasma current start-up under EBW power in tokamak plasmas
Authors:
Erasmus J. du Toit,
Martin R. O'Brien,
Roddy G. L. Vann
Abstract:
Tokamak plasma current start-up assisted by Electron Bernstein waves (EBW) has been demonstrated successfully in a number of experiments. The dynamic start-up phase involves a change in field topology, as the initially open magnetic field lines form closed flux surfaces (CFS) under the initiation of a plasma current. This change in field topology will bring about a change in the current drive (CD)…
▽ More
Tokamak plasma current start-up assisted by Electron Bernstein waves (EBW) has been demonstrated successfully in a number of experiments. The dynamic start-up phase involves a change in field topology, as the initially open magnetic field lines form closed flux surfaces (CFS) under the initiation of a plasma current. This change in field topology will bring about a change in the current drive (CD) mechanism, and, although various mechanisms have been proposed to explain the formation of CFS, no detailed theoretical studies have previously been undertaken. Here, we report on the development of a start-up model for EBW-assisted plasma current start-up in MAST. It is shown that collisions are responsible for only a small part of the CD, while the open magnetic field line configuration leads to an asymmetric confinement of electrons, which is responsible for the greater part of the CD.
△ Less
Submitted 18 May, 2017;
originally announced May 2017.
-
A kinetic study of microwave start-up of tokamak plasmas
Authors:
Erasmus J. du Toit,
Martin R. O'Brien,
Roddy G. L. Vann
Abstract:
A kinetic model for studying the time evolution of the distribution function for microwave start-up is presented. The model for the distribution function is two dimensional in momentum space, but, for simplicity and rapid calculations, has no spatial dependence. Experiments on the Mega Amp Spherical Tokamak have shown that the plasma current is carried mainly by electrons with energies greater tha…
▽ More
A kinetic model for studying the time evolution of the distribution function for microwave start-up is presented. The model for the distribution function is two dimensional in momentum space, but, for simplicity and rapid calculations, has no spatial dependence. Experiments on the Mega Amp Spherical Tokamak have shown that the plasma current is carried mainly by electrons with energies greater than $70$ keV, and effects thought to be important in these experiments are included, i.e. particle sources, orbital losses, the loop voltage and microwave heating, with suitable volume averaging where necessary to give terms independent of spatial dimensions. The model predicts current carried by electrons with the same energies as inferred from the experiments, though the current drive efficiency is smaller.
△ Less
Submitted 3 April, 2017;
originally announced April 2017.
-
Perturbing microwave beams by plasma density fluctuations
Authors:
Alf Köhn,
Eberhard Holzhauer,
Jarrod Leddy,
Matthew B. Thomas,
Roddy G. L. Vann
Abstract:
The propagation of microwaves across a turbulent plasma density layer is investigated with full-wave simulations. To properly represent a fusion edge-plasma, drift-wave turbulence is considered based on the Hasegawa-Wakatani model. Scattering and broadening of a microwave beam whose amplitude distribution is of Gaussian shape is studied in detail as a function of certain turbulence properties. Par…
▽ More
The propagation of microwaves across a turbulent plasma density layer is investigated with full-wave simulations. To properly represent a fusion edge-plasma, drift-wave turbulence is considered based on the Hasegawa-Wakatani model. Scattering and broadening of a microwave beam whose amplitude distribution is of Gaussian shape is studied in detail as a function of certain turbulence properties. Parameters leading to the strongest deterioration of the microwave beam are identified and implications for existing experiments are given.
△ Less
Submitted 15 March, 2017;
originally announced March 2017.
-
Kinetic simulations of X-B and O-X-B mode conversion and its deterioration at high input power
Authors:
A. V. Arefiev,
I. Y. Dodin,
A. Köhn,
E. J. Du Toit,
E. Holzhauer,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
Spherical tokamak plasmas are typically overdense and thus inaccessible to externally-injected microwaves in the electron cyclotron range. The electrostatic electron Bernstein wave (EBW), however, provides a method to access the plasma core for heating and diagnostic purposes. Understanding the details of the coupling process to electromagnetic waves is thus important both for the interpretation o…
▽ More
Spherical tokamak plasmas are typically overdense and thus inaccessible to externally-injected microwaves in the electron cyclotron range. The electrostatic electron Bernstein wave (EBW), however, provides a method to access the plasma core for heating and diagnostic purposes. Understanding the details of the coupling process to electromagnetic waves is thus important both for the interpretation of microwave diagnostic data and for assessing the feasibility of EBW heating and current drive. While the coupling is reasonably well-understood in the linear regime, nonlinear physics arising from high input power has not been previously quantified. To tackle this problem, we have performed one- and two-dimensional fully kinetic particle-in-cell simulations of the two possible coupling mechanisms, namely X-B and O-X-B mode conversion. We find that the ion dynamics has a profound effect on the field structure in the nonlinear regime, as high amplitude short-scale oscillations of the longitudinal electric field are excited in the region below the high-density cut-off prior to the arrival of the EBW. We identify this effect as the instability of the X wave with respect to resonant scattering into an EBW and a lower-hybrid wave. We calculate the instability rate analytically and find this basic theory to be in reasonable agreement with our simulation results.
△ Less
Submitted 22 December, 2016;
originally announced December 2016.
-
Influence of plasma turbulence on microwave propagation
Authors:
Alf Köhn,
Eberhard Holzhauer,
Jarrod Leddy,
Matthew B Thomas,
Roddy G L Vann
Abstract:
It is not fully understood how electromagnetic waves propagate through plasma density fluctuations when the size of the fluctuations is comparable with the wavelength of the incident radiation. In this paper, the perturbing effect of a turbulent plasma density layer on a traversing microwave beam is simulated with full-wave simulations. The deterioration of the microwave beam is calculated as a fu…
▽ More
It is not fully understood how electromagnetic waves propagate through plasma density fluctuations when the size of the fluctuations is comparable with the wavelength of the incident radiation. In this paper, the perturbing effect of a turbulent plasma density layer on a traversing microwave beam is simulated with full-wave simulations. The deterioration of the microwave beam is calculated as a function of the characteristic turbulence structure size, the turbulence amplitude, the depth of the interaction zone and the size of the waist of the incident beam. The maximum scattering is observed for a structure size on the order of half the vacuum wavelength. The scattering and beam broadening was found to increase linearly with the depth of the turbulence layer and quadratically with the fluctuation strength. Consequences for experiments and 3D effects are considered.
△ Less
Submitted 1 June, 2016; v1 submitted 1 April, 2016;
originally announced April 2016.
-
Observations of 2D Doppler backscattering on MAST
Authors:
D. A. Thomas,
K. J. Brunner,
S. J. Freethy,
B. K. Huang,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
The Synthetic Aperture Microwave Imaging (SAMI) diagnostic has conducted proof-of-principle 2D Doppler backscattering (DBS) experiments on MAST. SAMI actively probes the plasma edge using a wide (+-40 degrees vertical and horizontal) and tuneable (10-35.5 GHz) beam. The Doppler backscattered signal is digitised in vector form using an array of eight Vivaldi PCB antennas. This allows the receiving…
▽ More
The Synthetic Aperture Microwave Imaging (SAMI) diagnostic has conducted proof-of-principle 2D Doppler backscattering (DBS) experiments on MAST. SAMI actively probes the plasma edge using a wide (+-40 degrees vertical and horizontal) and tuneable (10-35.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 degrees FWHM at 10-34.5 GHz. This capability is unique to SAMI and is an entirely novel way of conducting DBS experiments. In this paper the feasibility of conducting 2D DBS experiments is explored. Initial measurements of phenomena observed on 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.
△ Less
Submitted 23 July, 2015;
originally announced July 2015.
-
Kinetic simulations of X-B and O-X-B mode conversion
Authors:
A. V. Arefiev,
E. J. Du Toit,
A. Köhn,
E. Holzhauer,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
We have performed fully-kinetic simulations of X-B and O-X-B mode conversion in one and two dimensional setups using the PIC code EPOCH. We have recovered the linear dispersion relation for electron Bernstein waves by employing relatively low amplitude incoming waves. The setups presented here can be used to study non-linear regimes of X-B and O-X-B mode conversion.
We have performed fully-kinetic simulations of X-B and O-X-B mode conversion in one and two dimensional setups using the PIC code EPOCH. We have recovered the linear dispersion relation for electron Bernstein waves by employing relatively low amplitude incoming waves. The setups presented here can be used to study non-linear regimes of X-B and O-X-B mode conversion.
△ Less
Submitted 6 May, 2015;
originally announced May 2015.
-
Electron kinetics inferred from observations of microwave bursts during edge localised modes in the Mega-Amp Spherical Tokamak
Authors:
S. J. Freethy,
K. G. McClements,
S. C. Chapman,
R. O. Dendy,
W. N. Lai,
S. J. P. Pamela,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM c…
▽ More
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM cycle. The observed radiative effects paradoxically imply acceleration of electrons parallel to the magnetic field combined with rapid acquisition of perpendicular momentum. It is shown that this paradox can be resolved by the action of the anomalous Doppler instability which enables fast collective radiative relaxation, in the perpendicular direction, of electrons accelerated in the parallel direction by inductive electric fields generated by the initial ELM instability.
△ Less
Submitted 4 November, 2014;
originally announced November 2014.
-
Propagation in 3D of microwaves through density perturbations
Authors:
T. R. N. Williams,
A. Köhn,
M. R. O'Brien,
R. G. L. Vann
Abstract:
Simulations using 3D and 2D full-wave codes have shown that edge filaments in tokamak plasmas can significantly affect the propagation of microwaves across a broad frequency spectrum, resulting in scattering angles of up to 46 degrees. Parameter scans were carried out for density perturbations comparable in width and amplitude to MAST filaments and the effect on the measured emission was calculate…
▽ More
Simulations using 3D and 2D full-wave codes have shown that edge filaments in tokamak plasmas can significantly affect the propagation of microwaves across a broad frequency spectrum, resulting in scattering angles of up to 46 degrees. Parameter scans were carried out for density perturbations comparable in width and amplitude to MAST filaments and the effect on the measured emission was calculated. 3D effects were discovered in the case of an obliquely incident beam.
In general, the problem of EM propagation past wavelength-sized 3D inhomogeneities is not well understood, yet is of importance for both heating and diagnostic applications in the electron cyclotron frequency range for tokamaks, as well as atmospheric physics. To improve this understanding, a new cold-plasma code, EMIT-3D, was written to extend full-wave microwave simulations in magnetized plasmas to 3D, and make comparisons to the existing 2D code IPF-FDMC. This work supports MAST experiments using the SAMI diagnostic to image microwave emission from the plasma edge due to mode conversion from electron Bernstein waves. Significant fluctuations in the SAMI data mean that detailed modelling is required to improve its interpretation.
△ Less
Submitted 28 November, 2013;
originally announced November 2013.
-
Synthetic aperture microwave imaging with active probing for fusion plasma diagnostics
Authors:
Vladimir F. Shevchenko,
Roddy G. L. Vann,
Simon J. Freethy,
Billy K. Huang
Abstract:
A Synthetic Aperture Microwave Imaging (SAMI) system has been designed and built to obtain 2-D images at several frequencies from fusion plasmas. SAMI uses a phased array of linearly polarised antennas. The array configuration has been optimised to achieve maximum synthetic aperture beam efficiency. The signals received by antennas are down-converted to the intermediate frequency range and then re…
▽ More
A Synthetic Aperture Microwave Imaging (SAMI) system has been designed and built to obtain 2-D images at several frequencies from fusion plasmas. SAMI uses a phased array of linearly polarised antennas. The array configuration has been optimised to achieve maximum synthetic aperture beam efficiency. The signals received by antennas are down-converted to the intermediate frequency range and then recorded in a full vector form. Full vector signals allow beam focusing and image reconstruction in both real time and a post processing mode. SAMI can scan over 16 preprogrammed frequencies in the range of 10-35GHz with a switching time of 300ns. The system operates in 2 different modes simultaneously: both a passive imaging of plasma emission and also an active imaging of the back-scattered signal of the radiation launched by one of the antennas from the same array. This second mode is similar to so-called Doppler backscattering (DBS) reflectometry with 2-D resolution of the propagation velocity of turbulent structures. Both modes of operation show good performance in a real fusion plasma experiments on Mega Amp Spherical Tokamak (MAST). We have obtained the first ever 2-D images of BXO mode conversion windows. With active probing, the first ever turbulence velocity maps have been obtained. In this article we present an overview of the diagnostic and discuss recent results.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.