Hasegawa Wakatani Research Articles

Anomalous transport of multi-species edge plasma with the generalized Hasegawa–Wakatani model and the FLR effects

Anomalous transport of multi-species edge plasma with the generalized Hasegawa–Wakatani model and the FLR effects

StyleGAN as an AI deconvolution operator for large eddy simulations of turbulent plasma equations in BOUT++

We present the use of StyleGAN, a face-synthesis generative adversarial network (GAN) developed by NVidia, as a deconvolution operator for large eddy simulation (LES) of plasma turbulence. The overall methodology, named style eddy simulation, has been integrated into the BOUT++ solver and tested on the original and modified Hasegawa–Wakatani models using different mesh sizes, 2562 and 5122, and different values of the adiabaticity parameter α and background density gradient κ. Using a LES resolution of 32 × 32 and 64 × 64, i.e., 64× smaller resolution than the corresponding direct numerical simulation (DNS), results show convergence toward the ground truth as we tighten the reconstruction tolerance, and an algorithm complexity O(N log N) is compared to the O(N2) of BOUT++. Finally, the trained GAN can be used to create valid initial conditions for a faster DNS by avoiding to start from nonphysical initial perturbations.

Resistive drift wave turbulence and anomalous transport of multi-species plasma

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.

On the Effects of Tokamak Plasma Edge Symmetries on Turbulence Relaxation

The plasma edge of a tokamak configuration is characterized by turbulent dynamics leading to enhanced transport. We construct a simplified 3D Hasegawa–Wakatani model reducing to a single partial differential equation for the turbulent electric potential dynamics. Simulations demonstrate how the 3D turbulence relaxes on a 2D axisymmetric profile, corresponding to the so-called interchange turbulence. The spectral features of this regime are found to be strongly dependent on the initialization pattern. We outline that the emergence of axisymmetric turbulence is also achieved when the corresponding mode amplitude is not initialized. Then, we introduce the symmetries of the magnetic X-point of a tokamak configuration. We linearize the governing equation by treating the poloidal field as a small correction. We show that it is not always possible to solve the electric potential dynamics following a perturbative approach. This finding, which is due to resonance between the modes of the background and the poloidal perturbation, confirms that the X-point symmetries can alter the properties of turbulent transport in the edge region.

Nearly integrable flows and chaotic tangles in the Dimits shift regime of plasma edge turbulence

Transitionally turbulent flows frequently exhibit spatiotemporal intermittency, reflecting a complex interplay between driving forces, dissipation, and transport present in these systems. When this intermittency manifests as observable structures and patterns in the flow, the characterization of turbulence in these systems becomes challenging due to the nontrivial correlations introduced into the statistics of the turbulence by these structures. In this work, we use tools from dynamical systems theory to study intermittency in the Dimits shift regime of the flux-balanced Hasegawa–Wakatani (BHW) equations, which models a transitional regime of resistive drift-wave turbulence relevant to magnetically confined fusion plasmas. First, we show in direct numerical simulations that turbulence in this regime is dominated by strong zonal flows and coherent drift-wave vortex structures, which maintain a strong linear character despite their large amplitude. Using the framework of generalized Liouville integrability, we develop a theory of integrable Lagrangian flows in generic fluid and plasma systems and discuss how the observed zonal flows plus drift waves in the BHW system exhibit a form of “near-integrability” originating from a fluid element relabeling symmetry. We further demonstrate that the BHW flows transition from integrability to chaos via the formation of chaotic tangles in the aperiodic Lagrangian flow, and establish a direct link between the “lobes” associated with these tangles and intermittency in the observed turbulent dissipation. This illustrates how utilizing tools from deterministic dynamical systems theory to study convective nonlinearities can explain aspects of the intermittent spatiotemporal structure exhibited by the statistics of turbulent fields.

Data-driven modal analysis of nonlinear quantities in turbulent plasmas using multi-field singular value decomposition

Nonlinear dynamics in the two-dimensional multi-component plasma turbulence described by the Hasegawa–Wakatani equation is investigated by using a data-driven modal analysis with the singular value decomposition (SVD). The conventional SVD is extended to ‘multi-field SVD’ which can decompose multiple turbulence fields simultaneously by a single set of orthonormal basis functions without imposing a priori scale separations. Then, in addition to the mode amplitude labeled by the singular value, the information on the phase relations in the nonlinear quantities such as a transport flux or a triad energy transfer is extracted in the mode space. Through applications to the two-dimensional plasma turbulence, it is revealed that the multi-field SVD can extract the dominant spatial structures for the turbulent transport and the nonlinear energy transfer, preserving the multi-scale nature of the original turbulent fields. It is also demonstrated that one can reduce the dimensionality or information using the multi-field SVD through comparisons with the conventional Fourier decomposition.

On the 3D turbulence regime in a Tokamak plasma edge

We derive a reduced model for the electrostatic turbulence in a Tokamak edge, when dealing with a resistive plasma and neglecting the spatial gradient of the background density which triggers the linear drift wave response. The obtained dynamics, de facto equivalent to a Hasegawa–Wakatani model, is characterized by a constitutive relation between the electric potential Laplacian and the density, which allows to deal with a single 3D equation governing the electric potential fluctuations. We study the evolution of the model, by separating the n=0 (n indicating toroidal-like number) mode from all the other ones. Then, we linearize the dynamics of the n≠0 modes around the steady 2D spectrum, which describes the spectral features of the electrostatic 2D turbulence. We theoretically and numerically demonstrate the existence of decaying branch of the 3D turbulence, having such a 2D steady spectrum as a natural attractor. This result suggests that the basic constituent of the self-sustained non-linear drift response in the plasma edge has to be individualized in the non-linear 2D electrostatic turbulence.

A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios

We formulate a generalized periodic boundary condition as a limit of the standard twist-and-shift parallel boundary condition that is suitable for simulations of plasmas with low magnetic shear. This is done by applying a phase shift in the binormal direction when crossing the parallel boundary. While this phase shift can be set to zero without loss of generality in the local flux-tube limit when employing the twist-and-shift boundary condition, we show that this is not the most general case when employing periodic parallel boundaries, and may not even be the most desirable. A non-zero phase shift can be used to avoid the convective cells that plague simulations of the three-dimensional Hasegawa–Wakatani system, and is shown to have measurable effects in periodic low-magnetic-shear gyrokinetic simulations. We propose a numerical program where a sampling of periodic simulations at random pseudo-irrational flux surfaces are used to determine physical observables in a statistical sense. This approach can serve as an alternative to applying the twist-and-shift boundary condition to low-magnetic-shear scenarios, which, while more straightforward, can be computationally demanding.

Data-driven model discovery for plasma turbulence modelling

An important problem in nuclear fusion plasmas is the prediction and control of turbulence which drives the cross-field transport, thus leading to energy loss from the system and deteriorating confinement. Turbulence, being a highly nonlinear and multiscale process, is challenging to theoretically describe and computationally model. Most advanced computational models fall into one of the two categories: fluid or gyro-kinetic. They both come at a high computational cost and cannot be applied for routine simulation of plasma discharge evolution and control. Development of reduced models based on (physics informed) artificial neural networks could potentially fulfil the need for affordable simulations of plasma turbulence. However, the training requires an extensive data base and the obtained models lack extrapolation capability to scenarios not originally encountered during training. This leads to reduced models of limited validity which may not prove adequate for predicting scenarios in future machines. In contrast, we explore a data-driven model discovery approach based on sparse regression to infer governing nonlinear partial differential equations directly from the data. Our input data are generated by simulations of drift-wave turbulence according to the Hasegawa–Wakatani and modified Hasegawa–Wakatani models. Balancing model accuracy and complexity enables the reconstruction of the systems of partial differential equations accurately describing the dynamics simulated in the input data sets. Sparse regression is not data hungry and can be extrapolated to unexplored parameter ranges. We explore and demonstrate the potential of this approach for fusion plasma turbulence modelling. The findings show that the methodology is promising for the development of reduced and computationally efficient turbulence models as well as for existing model cross-validation.

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. The density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

Phase and amplitude evolution in the network of triadic interactions of the Hasegawa–Wakatani system

The Hasegawa–Wakatani system, commonly used as a toy model of dissipative drift waves in fusion devices, is revisited with considerations of phase and amplitude dynamics of its triadic interactions. It is observed that a single resonant triad can saturate via three way phase locking, where the phase differences between dominant modes converge to constant values as individual phases increase in time. This allows the system to have approximately constant amplitude solutions. Non-resonant triads show similar behavior only when one of its legs is a zonal wave number. However, when an additional triad, which is a reflection of the original one with respect to the y axis is included, the behavior of the resulting triad pair is shown to be more complex. In particular, it is found that triads involving small radial wave numbers (large scale zonal flows) end up transferring their energy to the subdominant mode which keeps growing exponentially, while those involving larger radial wave numbers (small scale zonal flows) tend to find steady chaotic or limit cycle states (or decay to zero). In order to study the dynamics in a connected network of triads, a network formulation is considered, including a pump mode, and a number of zonal and non-zonal subdominant modes as a dynamical system. It was observed that the zonal modes become clearly dominant only when a large number of triads are connected. When the zonal flow becomes dominant as a “collective mean field,” individual interactions between modes become less important, which is consistent with the inhomogeneous wave-kinetic picture. Finally, the results of direct numerical simulation are discussed for the same parameters, and various forms of the order parameter are computed. It is observed that nonlinear phase dynamics results in a flattening of the large scale phase velocity as a function of scale in direct numerical simulations.

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or −90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s−1. The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. We show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

Evaluation of abrupt energy transfer among turbulent plasma structures using singular value decomposition

A method to quantify the energy transfer among turbulent structures using singular value decomposition (SVD) is presented. We apply the method to numerical turbulence data obtained from a global plasma simulation using the Hasegawa–Wakatani fluid model, in which the Kelvin–Helmholtz instability plays a dominant role. Using the SVD method, the electrostatic potential is decomposed into a background potential deformation, a zonal flow, a coherent mode and an intermittent structure. Thus there are four key structures, as distinct from the three found in conventional theory. The kinetic energy of each structure is evaluated, and the limit cycle among them is obtained. In the limit cycle, an abrupt change of the background is found to be synchronised with the period of the zonal flow. The energy transfer function of each turbulence structure, which is defined on the basis of a vorticity equation, is evaluated. This then provides physical understanding of how the limit cycle is sustained by dynamical changes in the energy transfer among structures over the its period. In addition, it is shown that the abrupt deformation of the background is caused by the non-linear self-coupling of the intermittent structure.

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. In this work, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

Nonlinear interaction and turbulence transition in the limiting regimes of plasma edge turbulence

We study the nonlinear coupling mechanism and turbulent transition in magnetically confined plasma flows based on two representative limiting regime dynamics. The two-field flux-balanced Hasegawa–Wakatani (BHW) model is taken as a simplified approximation to the key physical processes in the energy-conserving nonlinear plasma flows. The limiting regimes separate the effects of finite non-adiabatic resistivity and extreme non-normal dynamics to enable a more detailed investigation on each individual aspect with the help of various mathematical tools. We adopt the strategy from the selective decay theory used for the simpler one-field system to develop new crucial a priori estimations in the two-field model framework. The competing effects from model dissipation, finite particle resistivity, as well as the nonlinear interaction with a zonal mean state to induce dual direction energy transports are characterized from the systematic analysis. Non-normal dynamics with aligned eigendirections is also shown to go through a sharp transition from turbulence to regularized zonal flows. The diverse phenomena implied from the limiting regime analysis are further confirmed from direct numerical simulations of the BHW model.

The influence of magnetic field curvature on intermittency in drift-wave turbulence in the stellarator TJ-K

The influence of magnetic field curvature on the intermittency in density and potential fluctuations in drift-wave turbulence has been investigated in the stellarator TJ-K in the framework of an extended Hasegawa–Wakatani model. A structure function analysis was used to estimate the intermittency level of poloidally resolved drift-wave turbulence measurements of a deuterium plasma in TJ-K. Potential fluctuations were found to be broadly self-similar, whereas density fluctuations were found to be more intermittent in the region with negative normal and positive geodesic curvatures. This behavior could be understood by comparing the data to two-dimensional extended Hasegawa–Wakatani simulations that retain gradients in the magnetic field strength, giving rise to curvature effects. The model is able to reproduce the trends in the experimental data if both normal and geodesic curvature effects are accounted for, as well as the local anisotropy of turbulent length scales. The analysis indicates the importance of local magnetic geometry as a factor in the decoupling of density and potential fluctuations, leading to intermittency in drift-wave turbulence.

Solitary zonal structures in subcritical drift waves: a minimum model

Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same ‘equation of state’, which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.

Flux-balanced two-field plasma edge turbulence in a channel geometry

We investigate the drift wave–zonal flow interaction formulated on a channel domain geometry approximating an isolated plasma edge region with zero net radial transport across the boundary. The recent two-field flux-balanced Hasegawa–Wakatani (BHW) model with improved treatment for a parallel electron response is adapted to the channel geometry configuration, which allows for generalized non-uniform zonal density profiles and is valid for the simulation of turbulent fields near the tokamak or stellarator edge regions. New conserved quantities are constructed based on the channel geometry to help the analysis for the competition between zonal states and non-zonal fluctuations. Effective bounds can be found constraining the maximum growth of total fluctuations and the amplitude of the dominant zonal state based on the conserved quantities. Total statistical variance among all the modes can also be estimated depending on the zonal state strength. The theoretical discoveries are confirmed by detailed numerical experiments from simulations in the channel domain. In addition, the channel geometry provides further support for the important advantage of adopting the balanced flux correction in the BHW model by showing a physically consistent growth rate from a stability analysis for the small-amplitude fluctuation interaction with a prescribed zonal mean profile, in comparison with the persistent instability and strong outward transport found in the modified Hasegawa–Wakatani model even with the increasing zonal density profile. This is again confirmed by direct numerical simulations of the two models. The channel domain BHW model framework with attractive features implies many potential applications in the study of the complex phenomena in plasma edge turbulence.