-
Using Convolutional Neural Networks to detect Edge Localized Modes in DIII-D from Doppler Backscattering measurements
Authors:
N. Q. X. Teo,
V. H. Hall-Chen,
K. Barada,
R. J. H. Ng,
L. Gu,
A. K. Yeoh,
Q. T. Pratt,
X. Garbet,
T. L. Rhodes
Abstract:
In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods $\unicode{x2013}$ D-alpha spectr…
▽ More
In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods $\unicode{x2013}$ D-alpha spectroscopy and Doppler backscattering (DBS). The former detects ELMs by measuring Balmer alpha emission while the latter uses microwave radiation to probe the plasma. DBS has the advantage of having higher temporal resolution and robustness to damage. These advantages of DBS diagnostics may be beneficial for future operational tokamaks and thus data processing techniques for DBS should be developed in preparation. In sight of this, we explore the training of neural networks to detect ELMs from DBS data, using D-alpha data as the ground truth. With shots found in the DIII-D database, the model is trained to classify each time step based on the occurrence of an ELM event. The results are promising. When tested on shots similar to those used for training, the model is capable of consistently achieving a high f1-score of 0.93. This score is a performance metric for imbalanced datasets that ranges between 0 and 1. We evaluate the performance of our neural network on a variety of ELMs $\unicode{x2013}$ grasssy, suppressed, and wide pedestal $\unicode{x2013}$ finding broad applicability. Beyond ELMs, our work demonstrates the wider feasibility of applying neural networks to data from DBS diagnostics.
△ Less
Submitted 3 July, 2024; v1 submitted 3 June, 2024;
originally announced June 2024.
-
ELM-free Enhanced Dα H-mode with Near Zero NBI Torque Injection in DIII-D Tokamak
Authors:
T. Macwan,
K. Barada,
J. F. Parisi,
R. Groebner,
T. L. Rhodes,
S. Banerjee,
C. Chrystal,
Q. Pratt,
Z. Yan,
H. Wang,
L. Zeng,
M. E. Austin,
N. A. Crocker,
W. A. Peebles
Abstract:
Enhanced $D_α$ H-mode (EDA H-mode), an ELM-free H-mode regime, is explored in neutral beam heated, lower single null plasmas with near zero torque injection. This regime exhibits a good energy confinement ($\mathrm{H}_{\mathrm{98y2}}$ $\sim 1$) with $β_N \sim 2$, high density, regime access at low input power, and no ELMs. This paper further presents the time-resolved measurements of electron and…
▽ More
Enhanced $D_α$ H-mode (EDA H-mode), an ELM-free H-mode regime, is explored in neutral beam heated, lower single null plasmas with near zero torque injection. This regime exhibits a good energy confinement ($\mathrm{H}_{\mathrm{98y2}}$ $\sim 1$) with $β_N \sim 2$, high density, regime access at low input power, and no ELMs. This paper further presents the time-resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field during the EDA H-mode phase and examines the dynamics of the edge quasi-coherent mode (QCM). Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of $k_{\perp}ρ_s \sim 0.1-1.2$ with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where $k_{\perp}$ is the binormal wavenumber and $ρ_s$ is the ion sound radius. Linear gyrokinetic simulations of an EDA H-mode discharge with CGYRO indicates that the trapped electron mode (TEM) and electron temperature gradient (ETG) are dominant instabilities in the region where QCM is unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimental observed characteristics of the QCM. These similarities suggest that the QCM might be a TEM instability existing in the edge region of the EDA H-mode plasmas.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
Stability and Transport of Gyrokinetic Critical Pedestals
Authors:
J. F. Parisi,
A. O. Nelson,
W. Guttenfelder,
R. Gaur,
J. W. Berkery,
S. M. Kaye,
K. Barada,
C. Clauser,
A. Diallo,
D. R. Hatch,
A. Kleiner,
M. Lampert,
T. Macwan,
J. E. Menard
Abstract:
A gyrokinetic threshold model for pedestal width-height scaling prediction is applied to multiple devices and to a shaping and aspect-ratio scan giving $Δ_{\mathrm{ped}} = 0.92 A^{1.04} κ^{-1.24} 0.38^δ β_{θ,\mathrm{ped}}^{1.05}$ for pedestal width $Δ_{\mathrm{ped}}$, aspect-ratio $A$, elongation $κ$, triangularity $δ$, and normalized pedestal height $β_{θ,\mathrm{ped}}$. We also find a width-tran…
▽ More
A gyrokinetic threshold model for pedestal width-height scaling prediction is applied to multiple devices and to a shaping and aspect-ratio scan giving $Δ_{\mathrm{ped}} = 0.92 A^{1.04} κ^{-1.24} 0.38^δ β_{θ,\mathrm{ped}}^{1.05}$ for pedestal width $Δ_{\mathrm{ped}}$, aspect-ratio $A$, elongation $κ$, triangularity $δ$, and normalized pedestal height $β_{θ,\mathrm{ped}}$. We also find a width-transport scaling $Δ_{\mathrm{ped} } = 0.028 \left(q_e/Γ_e - 1.7 \right)^{1.5} \sim η_e ^{1.5}$ where $q_e$ and $Γ_e$ are turbulent electron heat and particle fluxes and $η_e = \nabla \ln T_e / \nabla \ln n_e$ for electron temperature $T_e$ and density $n_e$. Pedestals close to those limited by kinetic-ballooning-modes (KBMs) have modified turbulent transport properties compared to strongly driven KBMs. The role of flow shear is studied as a width-height scaling constraint and pedestal saturation mechanism for a standard and wide pedestal discharge.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Kinetic-Ballooning-Bifurcation in Tokamak Pedestals Across Shaping and Aspect-Ratio
Authors:
J. F. Parisi,
A. O. Nelson,
R. Gaur,
S. M. Kaye,
F. I. Parra,
J. W. Berkery,
K. Barada,
C. Clauser,
A. J. Creely,
A. Diallo,
W. Guttenfelder,
J. W. Hughes,
L. A. Kogan,
A. Kleiner,
A. Q. Kuang,
M. Lampert,
T. Macwan,
J. E. Menard,
M. A. Miller
Abstract:
We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers pot…
▽ More
We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.
△ Less
Submitted 7 April, 2024; v1 submitted 8 December, 2023;
originally announced December 2023.
-
Characterization of mesoscopic turbulent transport events with long-radial-range correlation in DIII-D H-mode plasmas
Authors:
R. Hong,
T. L. Rhodes,
Y. Ren,
P. H. Diamond,
X. Jian,
L. Zeng,
K. Barada,
Z. Yan,
G. R. McKee
Abstract:
A dimensionless collisionality scan has been performed in H-mode plasmas on DIII-D tokamak, with detailed measurements of intermediate-to-high wavenumber turbulence using Doppler backscattering systems. It is found that the shorter wavelength turbulence develops into spatially asymmetric turbulent structures with a long-radial-range correlation (LRRC) in the mid-radius region of high-collisionalit…
▽ More
A dimensionless collisionality scan has been performed in H-mode plasmas on DIII-D tokamak, with detailed measurements of intermediate-to-high wavenumber turbulence using Doppler backscattering systems. It is found that the shorter wavelength turbulence develops into spatially asymmetric turbulent structures with a long-radial-range correlation (LRRC) in the mid-radius region of high-collisionality discharges. Linear \textsc{cgyro} simulations indicate that the underlying turbulence is likely driven by the electron-temperature-gradient (ETG) mode. The LRRC transport events are highly intermittent and show a power spectrum of \(S_{\tilde{n}}(k_\perp) \propto k^{-1}_\perp\) for density fluctuations, which is often associated with self-organized criticality. The magnitude and the radial scale of those turbulent structures increase significantly when the $E_{r}\times B$ mean flow shearing rate decreases. The enhanced LRRC transport events appear to be correlated with the degraded energy confinement time. The emergence of such LRRC transport events may serve as a candidate explanation for the degrading nature of \emph{H}-mode core plasma confinement at high collisionality.
△ Less
Submitted 28 February, 2023;
originally announced March 2023.
-
Parallel Vertex Cover Algorithms on GPUs
Authors:
Peter Yamout,
Karim Barada,
Adnan Jaljuli,
Amer E. Mouawad,
Izzat El Hajj
Abstract:
Finding small vertex covers in a graph has applications in numerous domains. Two common formulations of the problem include: Minimum Vertex Cover, which finds the smallest vertex cover in a graph, and Parameterized Vertex Cover, which finds a vertex cover whose size is less than or equal to some parameter $k$. Algorithms for both formulations traverse a search tree, which grows exponentially with…
▽ More
Finding small vertex covers in a graph has applications in numerous domains. Two common formulations of the problem include: Minimum Vertex Cover, which finds the smallest vertex cover in a graph, and Parameterized Vertex Cover, which finds a vertex cover whose size is less than or equal to some parameter $k$. Algorithms for both formulations traverse a search tree, which grows exponentially with the size of the graph or the value of $k$.
Parallelizing the traversal of the vertex cover search tree on GPUs is challenging for multiple reasons. First, the search tree is a narrow binary tree which makes it difficult to extract enough sub-trees to process in parallel to fully utilize the GPU's resources. Second, the search tree is highly imbalanced which makes load balancing across a massive number of parallel GPU workers challenging. Third, keeping around all the intermediate state needed to traverse many sub-trees in parallel puts high pressure on the GPU's memory resources and may act as a limiting factor to parallelism.
To address these challenges, we propose an approach to traverse the vertex cover search tree in parallel using GPUs while handling dynamic load balancing. Each thread block traverses a different sub-tree using a local stack, however, we also use a global worklist to balance load. Blocks contribute branches of their sub-trees to the global worklist on an as-needed basis, while blocks that finish their sub-trees get new ones from the global worklist. We use degree arrays to represent intermediate graphs so that the representation is compact in memory to avoid limiting parallelism, but self-contained which is necessary for load balancing. Our evaluation shows that compared to prior work, our hybrid approach of using local stacks and a global worklist substantially improves performance and reduces load imbalance, especially on difficult instances of the problem.
△ Less
Submitted 21 April, 2022;
originally announced April 2022.
-
Observation of Long-Radial-Range-Correlation in Turbulence in High-Collisionality High-Confinement Fusion Plasmas
Authors:
R. Hong,
T. L. Rhodes,
P. H. Diamond,
Y. Ren,
L. Zeng,
X. Jian,
K. Barada,
G. Wang,
W. A. Peebles
Abstract:
We report on the observation of spatially asymmetric turbulent structures with a long radial correlation length in the core of high-collisionality H-mode plasmas on DIII-D tokamak. These turbulent structures develop from shorter wavelength turbulence and have a radially elongated structure. The envelope of turbulence spans a broad radial range in the mid-radius region, leading to streamer-like tra…
▽ More
We report on the observation of spatially asymmetric turbulent structures with a long radial correlation length in the core of high-collisionality H-mode plasmas on DIII-D tokamak. These turbulent structures develop from shorter wavelength turbulence and have a radially elongated structure. The envelope of turbulence spans a broad radial range in the mid-radius region, leading to streamer-like transport events. The underlying turbulence is featured by intermittency, long-term memory effect, and the characteristic spectrum of self-organized criticality. The amplitude and the radial scale increase substantially when the shearing rate of the mean flow is reduced below the turbulent scattering rate. The enhanced LRRC transport events are accompanied by apparent degradation of normalized energy confinement time. These findings constitute the first experimental observation of long-radial-range turbulent transport events in high-collisionality H-mode plasmas, and demonstrate the role of mean shear flows in the formation and propagation of turbulence with long-radial-range correlation.
△ Less
Submitted 15 February, 2023; v1 submitted 20 February, 2022;
originally announced February 2022.
-
Observation of quasi-coherent density fluctuation in scrape-off layer enhancing boundary transport in high-betaN hybrid plasmas on DIII-D
Authors:
R. Hong,
T. L. Rhodes,
Z. Y. Li,
H. Wang,
L. Zeng,
K. Barada,
G. Wang,
J. G. Watkins,
W. A. Peebles
Abstract:
We report the observation of a quasi-coherent density fluctuation (QCF) by the Doppler backscattering system in the scrape-off layer (SOL) region of the DIII-D tokamak. This QCF is observed in high-power, high-performance hybrid plasmas with near double-null divertor (DND) shape during the electron cyclotron heating period. This mode is correlated with a steepened SOL density profile and leads to…
▽ More
We report the observation of a quasi-coherent density fluctuation (QCF) by the Doppler backscattering system in the scrape-off layer (SOL) region of the DIII-D tokamak. This QCF is observed in high-power, high-performance hybrid plasmas with near double-null divertor (DND) shape during the electron cyclotron heating period. This mode is correlated with a steepened SOL density profile and leads to significantly elevated particle and heat fluxes between ELMs. The SOL QCF is a long-wavelength ion-scale fluctuation and propagates in the ion diamagnetic direction in the plasma frame. Its radial expanse is about 1.5-2 cm, well beyond the typical width of heat flux on DIII-D. Also, the SOL QCF does not show any clear dependence on the effective SOL collisionality and thus may raise issues on the control of plasma-material interactions in low collisionality plasmas in which the blob-induced transport is reduced. A linear simulation using BOUT++ with a 5-field reduced model is performed and compared with experimental observations. In simulation results, an interchange-like density perturbation can be driven by the SOL density gradient, and its peak location and the radial width of the density perturbation are in agreement with the experimental observations.
△ Less
Submitted 19 March, 2021; v1 submitted 10 November, 2020;
originally announced November 2020.
-
Understanding the working of a B-dot probe
Authors:
Sayak Bose,
Manjit Kaur,
Kshitish K Barada,
Joydeep Ghosh,
Prabal K Chattopadhyay,
Rabindranath Pal
Abstract:
Magnetic pickup loops or B-dot probes are one of the oldest known sensors of time-varying magnetic fields. The operating principle is based on Faraday's law of electromagnetic induction. However, obtaining accurate measurements of time-varying magnetic fields using these kinds of probes is a challenging task. A B-dot probe and its associated circuit are prone to electrical oscillations. As a resul…
▽ More
Magnetic pickup loops or B-dot probes are one of the oldest known sensors of time-varying magnetic fields. The operating principle is based on Faraday's law of electromagnetic induction. However, obtaining accurate measurements of time-varying magnetic fields using these kinds of probes is a challenging task. A B-dot probe and its associated circuit are prone to electrical oscillations. As a result, the measured signal may not faithfully represent the magnetic field sampled by the B-dot probe. In this paper, we have studied the transient response of a B-dot probe and its associated circuit to a time-varying magnetic field. Methods of removing the oscillations pertaining to the detector structure are described. After removing the source of the oscillatory signal, we have shown that the time-integrated induced emf measured by the digitiser is linearly proportional to the magnetic field sampled by the B-dot probe, thus verifying the faithfulness of the measured signal.
△ Less
Submitted 22 September, 2020;
originally announced September 2020.
-
Effect of Inhomogeneous magnetic field on Plasma generation in a low magnetic field helicon discharge
Authors:
Sonu Yadav,
Prabal K Chattopadhyay,
Kshitish K. Barada,
Soumen Ghosh,
Joydeep Ghosh
Abstract:
The ionization efficiency of helicon plasma discharge is explored by changing the low axial magnetic field gradients near the helicon antenna. The highest plasma density is found for a most possible diverging field near the antenna by keeping the other operating condition constant. Measurement of axial wave number together with estimated radial wavenumber suggests the oblique mode propagation of h…
▽ More
The ionization efficiency of helicon plasma discharge is explored by changing the low axial magnetic field gradients near the helicon antenna. The highest plasma density is found for a most possible diverging field near the antenna by keeping the other operating condition constant. Measurement of axial wave number together with estimated radial wavenumber suggests the oblique mode propagation of helicon wave along the resonance cone boundary. Propagation of helicon wave near the resonance cone angle boundary can excite electrostatic fluctuations which subsequently can deposit energy in the plasma. This process has been shown to be responsible for peaking in density in low field helicon discharges, where the helicon wave propagates at an angle with respect to the applied uniform magnetic field. The increased efficiency can be explained on the basis of multiple resonances for multimode excitation by the helicon antenna due to the availability of a broad range of magnetic field values in the near field of the antenna when a diverging magnetic field is applied in the source.
△ Less
Submitted 10 January, 2019;
originally announced January 2019.
-
Hollow density formation in magnetically expanding helicon plasma
Authors:
Sonu Yadav,
Soumen Ghosh,
Sayak Bose,
K. K. Barada,
R. pal,
P. K. Chattopadhyay
Abstract:
Measurement of radial density profile in both the source and expansion chambers of a helicon plasma device have revealed that it is always centrally peaked in the source chamber, whereas in the expansion chamber near the diverging magnetic field it becomes hollow above a critical value of the magnetic field. This value corresponds to that above which both electrons and ions become magnetized. The…
▽ More
Measurement of radial density profile in both the source and expansion chambers of a helicon plasma device have revealed that it is always centrally peaked in the source chamber, whereas in the expansion chamber near the diverging magnetic field it becomes hollow above a critical value of the magnetic field. This value corresponds to that above which both electrons and ions become magnetized. The temperature profile is always peaked off- axis and tail electrons are found at the peak location in both the source and expansion chambers. Rotation of the tail electrons in the azimuthal direction in the expansion chamber due to gradient-B drift produces more ionization off-axis and creates a hollow density profile; however, if the ions are not magnetized, the additional ionization does not cause hollowness.
△ Less
Submitted 20 February, 2018;
originally announced February 2018.
-
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.