-
Flexible Stellarator Physics Facility
Authors:
F. I. Parra,
S. -G. Baek,
M. Churchill,
D. R. Demers,
B. Dudson,
N. M. Ferraro,
B. Geiger,
S. Gerhardt,
K. C. Hammond,
S. Hudson,
R. Jorge,
E. Kolemen,
D. M. Kriete,
S. T. A. Kumar,
M. Landreman,
C. Lowe,
D. A. Maurer,
F. Nespoli,
N. Pablant,
M. J. Pueschel,
A. Punjabi,
J. A. Schwartz,
C. P. S. Swanson,
A. M. Wright
Abstract:
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Gyrokinetic investigation of toroidal Alfven eigenmode (TAE) turbulence
Authors:
Ajay C. J.,
Ben McMillan,
Arkaprava Bokshi,
Alessandro di Siena,
M. J. Pueschel,
Juan Ruiz Ruiz
Abstract:
Toroidal Alfvén eigenmodes (TAEs) can transport fusion-born energetic particles out of the plasma volume, thereby decreasing plasma self-heating efficiency and possibly damaging reactor walls. Therefore, understanding TAE destabilisation and identifying saturation mechanisms is crucial to achieving burning plasma. While TAEs have been studies extensively in the past using kinetic-MHD codes, here a…
▽ More
Toroidal Alfvén eigenmodes (TAEs) can transport fusion-born energetic particles out of the plasma volume, thereby decreasing plasma self-heating efficiency and possibly damaging reactor walls. Therefore, understanding TAE destabilisation and identifying saturation mechanisms is crucial to achieving burning plasma. While TAEs have been studies extensively in the past using kinetic-MHD codes, here a fully gyrokinetic study is employed which allows for additional physics. In the case studied, the primary drive mechanism is identified as the resonance between the magnetic drifts and the TAE, and this is seen to be disrupted by equilibrium flow shear which can stabilize the mode by rotating it in the the poloidal plane. It is found that zonal flows do not play a significant role in the saturation of these TAEs, and there are no saturation mechanisms present in the local gyrokinetic picture that are able to saturate the mode at physically relevant transport levels in the case of TAE-only turbulence. Instead, we confirm that the global profile flattening of fast-ion density is the key saturation mechanism. The nonlinear excitation of TAE travelling along the electron diamagnetic direction and its beating with the ion diamagnetic TAE, resulting in large amplitude oscillations that may help detect TAEs more easily in tokamaks, is also reported.
△ Less
Submitted 29 April, 2024;
originally announced April 2024.
-
On the effect of flux-surface shaping on trapped-electron modes in quasi-helically symmetric stellarators
Authors:
M. J. Gerard,
M. J. Pueschel,
B. Geiger,
R. J. J. Mackenbach,
J. M. Duff,
B. J. Faber,
C. C. Hegna,
P. W. Terry
Abstract:
Using a novel optimization procedure it has been shown that the Helically Symmetric eXperiment (HSX) stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [M.J.~Gerard et al., \textit{Nucl.~Fusion} \textbf{63} (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available…
▽ More
Using a novel optimization procedure it has been shown that the Helically Symmetric eXperiment (HSX) stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [M.J.~Gerard et al., \textit{Nucl.~Fusion} \textbf{63} (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available energy $E_\mathrm{A}$, quasi-helical symmetry, and flux-surface shaping parameters as metrics for TEM stabilization. It is found that lower values of $E_\mathrm{A}$ correlate with reduced growth rates, but only when separate flux-surface shaping regimes are considered. Moreover, configurations with improved quasi-helical symmetry demonstrate a similar reduction in growth rates and less scatter compared to $E_\mathrm{A}$. Regarding flux-surface shaping, a set of helical shaping parameters is introduced that show increased elongation is strongly correlated with reduced TEM growth rates, however, only when the quasi-helical symmetry is preserved. Using a newly derived velocity-space-averaged TEM resonance operator, these trends are analyzed to provide insights into the physical mechanism of the observed stabilization. For elongation, stabilization is attributed to geometric effects that reduce the destabilizing particle drifts across the magnetic field. Regarding quasi-helical symmetry, the TEM resonance in the maximally resonant trapping well is shown to increase as the quasi-helical symmetry is broken, and breaking quasi-helical symmetry increases the prevalence of highly resonant trapping wells. While these results demonstrate the limitations of using any single metric as a linear TEM proxy, it is shown that quasi-helical symmetry and plasma elongation are highly effective metrics for reducing TEM growth rates in helical equilibria.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Influence of collisions on trapped-electron modes in tokamaks and low-shear stellarators
Authors:
M. C. L. Morren,
J. H. E. Proll,
J. van Dijk,
M. J. Pueschel
Abstract:
The influence of collisions on the growth rate of trapped-electron modes (TEMs) in core plasmas is assessed through both analytical linear gyrokinetics and linear gyrokinetic simulations. Both methods are applied to the magnetic geometry of the DIII-D tokamak, as well as the Helically Symmetric eXperiment (HSX) and Wendelstein 7-X (W7-X) stellarators, in the absence of temperature gradients. Here…
▽ More
The influence of collisions on the growth rate of trapped-electron modes (TEMs) in core plasmas is assessed through both analytical linear gyrokinetics and linear gyrokinetic simulations. Both methods are applied to the magnetic geometry of the DIII-D tokamak, as well as the Helically Symmetric eXperiment (HSX) and Wendelstein 7-X (W7-X) stellarators, in the absence of temperature gradients. Here we analytically investigate the influence of collisions on the TEM eigenmode frequency by a perturbative approach in the response of trapped particles to the mode, using an energy-dependent Krook operator to model collisions. Although the resulting growth rates exceed perturbative thresholds, they reveal important qualitative dependencies: a geometry-dependent stabilization rate occurs for all wavenumbers at high collisionality, while at low collisionality, a geometry-sensitive mixture of collisionless, resonantly driven, and collisionally destabilized modes is found. Additionally, linear gyrokinetic simulations have been performed with a rigorous pitch-angle scattering operator for the same geometries. In the case of DIII-D and large wavenumber modes in HSX, the trends predicted by analytical theory are reproduced. Dissimilarities are, however, obtained in W7-X geometry and for low wavenumber modes in HSX, which are shown to be due to a collision-induced transition to the Universal Instability as the dominant instability at marginal collisionality.
△ Less
Submitted 24 May, 2024; v1 submitted 22 January, 2024;
originally announced January 2024.
-
Three-Dimensional Shear-Flow Instability Saturation via Stable Modes
Authors:
B. Tripathi,
P. W. Terry,
A. E. Fraser,
E. G. Zweibel,
M. J. Pueschel
Abstract:
Turbulence in three dimensions ($3$D) supports vortex stretching that has long been known to accomplish energy transfer to small scales. Moreover, net energy transfer from large-scale, forced, unstable flow-gradients to smaller scales is achieved by gradient-flattening instability. Despite such enforcement of energy transfer to small scales, it is shown here not only that the shear-flow-instabilit…
▽ More
Turbulence in three dimensions ($3$D) supports vortex stretching that has long been known to accomplish energy transfer to small scales. Moreover, net energy transfer from large-scale, forced, unstable flow-gradients to smaller scales is achieved by gradient-flattening instability. Despite such enforcement of energy transfer to small scales, it is shown here not only that the shear-flow-instability-supplied $3$D-fluctuation energy is largely inverse-transferred from the fluctuation to the mean-flow gradient, but that such inverse transfer is more efficient for turbulent fluctuations in $3$D than in two dimensions ($2$D). The transfer is due to linearly stable eigenmodes that are excited nonlinearly. The stable modes, thus, reduce both the nonlinear energy cascade to small scales and the viscous dissipation rate. The vortex-tube stretching is also suppressed. Up-gradient momentum transport by the stable modes counters the instability-driven down-gradient transport, which also is more effective in $3$D than in $2$D ($\mathrm{\approx} 70\% \mathrm{\,\, vs.\,\,}\mathrm{\approx} 50\%$). From unstable modes, these stable modes nonlinearly receive energy via zero-frequency fluctuations that vary only in the direction orthogonal to the plane of $2$D shear flow. The more widely occurring $3$D turbulence is thus inherently different from the commonly studied $2$D turbulence, despite both saturating via stable modes.
△ Less
Submitted 17 October, 2023; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Global Linear and Nonlinear Gyrokinetic Simulations of Tearing Modes
Authors:
T. Jitsuk,
A. Di Siena,
M. J. Pueschel,
P. W. Terry,
F. Widmer,
E. Poli,
J. S. Sarff
Abstract:
To better understand the interaction of global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code \textsc{Gene} is modified to describe global tearing mode instability via a shifted Maxwellian distribution consistent with experimental equilibria. The implementation of the shifted Maxwellian is tested and benchmarked by com…
▽ More
To better understand the interaction of global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code \textsc{Gene} is modified to describe global tearing mode instability via a shifted Maxwellian distribution consistent with experimental equilibria. The implementation of the shifted Maxwellian is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ($n=5$, $n=6$) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the $m=2$ tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.
△ Less
Submitted 12 October, 2023; v1 submitted 30 August, 2023;
originally announced August 2023.
-
Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence
Authors:
B. Tripathi,
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel,
E. H. Anders
Abstract:
To comprehensively understand saturation of two-dimensional ($2$D) magnetized Kelvin-Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a signif…
▽ More
To comprehensively understand saturation of two-dimensional ($2$D) magnetized Kelvin-Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasilinear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are: an explosive large-scale vortex separation instead of the well-known merger of $2$D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length-scale. These show critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.
△ Less
Submitted 17 July, 2023;
originally announced July 2023.
-
Enhanced transport at high plasma pressure and sub-threshold kinetic ballooning modes in Wendelstein 7-X
Authors:
P. Mulholland,
K. Aleynikova,
B. J. Faber,
M. J. Pueschel,
J. H. E. Proll,
C. C. Hegna,
P. W. Terry,
C. Nührenberg
Abstract:
High-performance fusion plasmas, requiring high pressure $β$, are not well-understood in stellarator-type experiments. Here, the effect of $β$ on ion-temperature-gradient-driven (ITG) turbulence is studied in Wendelstein 7-X (W7-X), showing that subdominant kinetic ballooning modes (KBMs) are unstable well below the ideal MHD threshold and get strongly excited in the turbulence. By zonal-flow eros…
▽ More
High-performance fusion plasmas, requiring high pressure $β$, are not well-understood in stellarator-type experiments. Here, the effect of $β$ on ion-temperature-gradient-driven (ITG) turbulence is studied in Wendelstein 7-X (W7-X), showing that subdominant kinetic ballooning modes (KBMs) are unstable well below the ideal MHD threshold and get strongly excited in the turbulence. By zonal-flow erosion, these sub-threshold KBMs (stKBMs) affect ITG saturation and enable higher heat fluxes. Controlling stKBMs will be essential to allow W7-X and future stellarators to achieve maximum performance.
△ Less
Submitted 25 September, 2023; v1 submitted 6 June, 2023;
originally announced June 2023.
-
Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence
Authors:
B. Tripathi,
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel
Abstract:
Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly…
▽ More
Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified, and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to $98\%$. These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the fall-off rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations.
△ Less
Submitted 5 August, 2022;
originally announced August 2022.
-
Microtearing turbulence saturation via electron temperature flattening at low-order rational surfaces
Authors:
Ajay C. J.,
Ben McMillan,
M. J. Pueschel
Abstract:
Microtearing instability is one of the major sources of turbulent transport in high-$β$ tokamaks. These modes lead to very localized transport at low-order rational magnetic field lines, and we show that flattening of the local electron temperature gradient at these rational surfaces plays an important role in setting the saturated flux level in microtearing turbulence. This process depends crucia…
▽ More
Microtearing instability is one of the major sources of turbulent transport in high-$β$ tokamaks. These modes lead to very localized transport at low-order rational magnetic field lines, and we show that flattening of the local electron temperature gradient at these rational surfaces plays an important role in setting the saturated flux level in microtearing turbulence. This process depends crucially on the density of rational surfaces, and thus the system-size, and gives rise to a worse-than-gyro-Bohm transport scaling for system-sizes typical of existing tokamaks and simulations.
△ Less
Submitted 18 May, 2023; v1 submitted 19 July, 2022;
originally announced July 2022.
-
Mechanism for Sequestering Magnetic Energy at Large Scales in Shear-Flow Turbulence
Authors:
B. Tripathi,
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel
Abstract:
Straining of magnetic fields by large-scale shear flow, generally assumed to lead to intensification and generation of small scales, is re-examined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that, in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales, due to countera…
▽ More
Straining of magnetic fields by large-scale shear flow, generally assumed to lead to intensification and generation of small scales, is re-examined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that, in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales, due to counteracting straining motion of nonlinearly excited large-scale stable eigenmodes. This effect is quantified via dissipation rates, energy transfer rates, and visualizations of magnetic field evolution by artificially removing the stable modes. These analyses show that predictions based upon physics of the linear instability alone miss substantial dynamics, including those of magnetic fluctuations.
△ Less
Submitted 13 June, 2022; v1 submitted 3 May, 2022;
originally announced May 2022.
-
Reduced models for ETG transport in the pedestal
Authors:
David R. Hatch,
Craig Michoski,
Dongyang Kuang,
Ben Chapman-Oplopoiou,
Max Curie,
Michael Halfmoon,
Ehab Hassan,
Mike Kotschenreuther,
Swadesh M. Mahajan,
Gabriele Merlo,
M. J. Pueschel,
Justin Walker,
Cole D. Stephens
Abstract:
This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from the pedestals in a broad range of experimental scenarios. The simulation data has been consolidated in a new database for gyrokinetic simulation data, the Mult…
▽ More
This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from the pedestals in a broad range of experimental scenarios. The simulation data has been consolidated in a new database for gyrokinetic simulation data, the Multiscale Gyrokinetic Database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter eta, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an RMS error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.
△ Less
Submitted 19 January, 2022;
originally announced January 2022.
-
Anisotropic cosmic-ray diffusion in isotropic Kolmogorov turbulence
Authors:
P. Reichherzer,
J. Becker Tjus,
E. G. Zweibel,
L. Merten,
M. J. Pueschel
Abstract:
Understanding the time scales for diffusive processes and their degree of anisotropy is essential for modelling cosmic-ray transport in turbulent magnetic fields. We show that the diffusion time scales are isotropic over a large range of energy and turbulence levels, notwithstanding the high degree of anisotropy exhibited by the components of the diffusion tensor for cases with an ordered magnetic…
▽ More
Understanding the time scales for diffusive processes and their degree of anisotropy is essential for modelling cosmic-ray transport in turbulent magnetic fields. We show that the diffusion time scales are isotropic over a large range of energy and turbulence levels, notwithstanding the high degree of anisotropy exhibited by the components of the diffusion tensor for cases with an ordered magnetic field component. The predictive power of the classical scattering relation as a description for the relation between the parallel and perpendicular diffusion coefficients is discussed and compared to numerical simulations. Very good agreement for a large parameter space is found, transforming classical scattering relation predictions into a computational prescription for the perpendicular component. We discuss and compare these findings and in particular the time scales to become diffusive with the time scales that particles reside in astronomical environments, the so-called escape time scales. The results show that, especially at high energies, the escape times obtained from diffusion coefficients may exceed the time scales required for diffusion. In these cases, the escape time cannot be determined by the diffusion coefficients.
△ Less
Submitted 22 December, 2021;
originally announced December 2021.
-
Implications of turbulence-dependent diffusion on cosmic-ray spectra
Authors:
J. Dörner,
P. Reichherzer,
L. Merten,
J. Becker Tjus,
H. Fichtner,
M. J. Pueschel,
E. G. Zweibel
Abstract:
The propagation of cosmic rays can be described as a diffusive motion in most galactic environments. High-energy gamma-rays measured by Fermi have allowed inference of a gradient in the cosmic-ray density and spectral energy behavior in the Milky Way, which is not predicted by models. Here, a turbulence-dependent diffusion model is used to probe different types of cosmic-ray diffusion tensors. Cru…
▽ More
The propagation of cosmic rays can be described as a diffusive motion in most galactic environments. High-energy gamma-rays measured by Fermi have allowed inference of a gradient in the cosmic-ray density and spectral energy behavior in the Milky Way, which is not predicted by models. Here, a turbulence-dependent diffusion model is used to probe different types of cosmic-ray diffusion tensors. Crucially, it is demonstrated that the observed gradients can be explained through turbulence-dependent energy-scaling of the diffusion tensor.
△ Less
Submitted 13 October, 2021;
originally announced October 2021.
-
Regimes of cosmic-ray diffusion in Galactic turbulence
Authors:
P. Reichherzer,
L. Merten,
J. Dörner,
J. Becker Tjus,
M. J. Pueschel,
E. G. Zweibel
Abstract:
Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process, which is two-fold turbulent mixing in that the particle motion is stochastic with respect to the field lines, needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspect…
▽ More
Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process, which is two-fold turbulent mixing in that the particle motion is stochastic with respect to the field lines, needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at $>5\,$kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.
△ Less
Submitted 12 December, 2021; v1 submitted 27 April, 2021;
originally announced April 2021.
-
Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators
Authors:
J. Smoniewski,
E. Sánchez,
I. Calvo,
M. J. Pueschel,
J. N. Talmadge
Abstract:
The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (N…
▽ More
The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate $k_x$, but the oscillation and damping time scales are longer in local representations, particularly for small $k_x$ approaching the system size.
△ Less
Submitted 13 March, 2021; v1 submitted 22 December, 2020;
originally announced December 2020.
-
Importance of theory, computation and predictive modeling in the US magnetic fusion energy strategic plan
Authors:
Fatima Ebrahimi,
Gary Staebler,
Paul Bonoli,
Francois Waelbroeck,
Chris Hegna,
Lynda LoDestro,
M. J. Pueschel,
I. Joseph
Abstract:
Based on the community input at the National Academy of Sciences (NAS) Madison and Austin workshops in July and December 2017, respectively, this whitepaper was prepared and submitted to the NAS in the category of US fusion theory and computation. This whitepaper was submitted to NAS as one of five community-approved whitepapers. The revised version was also submitted for the Knoxville American Ph…
▽ More
Based on the community input at the National Academy of Sciences (NAS) Madison and Austin workshops in July and December 2017, respectively, this whitepaper was prepared and submitted to the NAS in the category of US fusion theory and computation. This whitepaper was submitted to NAS as one of five community-approved whitepapers. The revised version was also submitted for the Knoxville American Physical Society Division of Plasma Physics Community Planning Process (APS-DPP-CPP) workshop in September 2019.
△ Less
Submitted 26 November, 2020;
originally announced November 2020.
-
The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes
Authors:
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel,
J. M. Schroeder
Abstract:
The Kelvin-Helmholtz (KH) instability of a shear layer with an initially-uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it.…
▽ More
The Kelvin-Helmholtz (KH) instability of a shear layer with an initially-uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with KH-unstable flows in hydrodynamics, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the magnetic field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited and thus transport less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear due to stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.
△ Less
Submitted 21 October, 2020;
originally announced October 2020.
-
Microtearing modes as the source of magnetic fluctuations in the JET pedestal
Authors:
D. R. Hatch,
M. Kotschenreuther,
S. M. Mahajan,
M. J. Pueschel,
C. Michoski,
G. Merlo,
E. Hassan,
A. R. Field,
L. Frassinetti,
C. Giroud,
J. C. Hillesheim,
C. F. Maggi,
C. Perez von Thun,
C. M. Roach,
S. Saarelma,
D. Jarema,
F. Jenko,
JET contributors
Abstract:
We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data, to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands…
▽ More
We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data, to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands with corresponding toroidal mode numbers n=4 and 8. Such disparate fluctuation scales, with substantial gaps between toroidal mode numbers, are commonly observed in pedestal fluctuations. Here we provide a clear explanation, namely the alignment of the relevant rational surfaces (and not others) with the peak in the omega star profile, which is localized in the steep gradient region of the pedestal. We demonstrate that a global treatment is required to capture this effect. Nonlinear simulations suggest that the MTM fluctuations produce experimentally-relevant transport levels and saturate by relaxing the background electron temperature gradient, slightly downshifting the fluctuation frequencies from the linear predictions. Scans in collisionality are compared with simple MTM dispersion relations. At the experimental points considered, MTM growth rates can either increase or decrease with collision frequency depending on the parameters thus defying any simple characterization of collisionality dependence.
△ Less
Submitted 14 July, 2020;
originally announced July 2020.
-
Suppression of turbulence by trapped electrons in optimised stellarators
Authors:
J. H. E. Proll,
P. Xanthopoulos,
P. Helander,
G. G. Plunk,
B. J. Faber,
T. Görler,
H. M. Smith,
M. J. Pueschel
Abstract:
In fusion devices, the geometry of the confining magnetic field has a significant impact on the instabilities that drive turbulent heat loss. This is especially true of stellarators, where the "trapped electron mode" (TEM) is stabilised if specific optimisation criteria are satisfied, as in the Wendelstein 7-X experiment (W7-X). Here we find, by numerical simulation, that W7-X indeed has low TEM-d…
▽ More
In fusion devices, the geometry of the confining magnetic field has a significant impact on the instabilities that drive turbulent heat loss. This is especially true of stellarators, where the "trapped electron mode" (TEM) is stabilised if specific optimisation criteria are satisfied, as in the Wendelstein 7-X experiment (W7-X). Here we find, by numerical simulation, that W7-X indeed has low TEM-driven transport, and also benefits from stabilisation of the ion-temperature-gradient mode, giving theoretical support for the existence of enhanced confinement regimes at finite density gradients.
△ Less
Submitted 24 February, 2020;
originally announced February 2020.
-
Turbulence-Level Dependence of Cosmic-Ray Parallel Diffusion
Authors:
P. Reichherzer,
J. Becker Tjus,
E. G. Zweibel,
L. Merten,
M. J. Pueschel
Abstract:
Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic-ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to tur…
▽ More
Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic-ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient's dependencies on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely-used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in the present work. Consequently, a proper description of cosmic-ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic-ray gradient problem.
△ Less
Submitted 18 August, 2020; v1 submitted 16 October, 2019;
originally announced October 2019.
-
Role of stable modes in driven shear-flow turbulence
Authors:
A. E. Fraser,
M. J. Pueschel,
P. W. Terry,
E. G. Zweibel
Abstract:
A linearly unstable, sinusoidal $E \times B$ shear flow is examined in the gyrokinetic framework in both the linear and nonlinear regimes. In the linear regime, it is shown that the eigenmode spectrum is nearly identical to hydrodynamic shear flows, with a conjugate stable mode found at every unstable wavenumber. In the nonlinear regime, turbulent saturation of the instability is examined with and…
▽ More
A linearly unstable, sinusoidal $E \times B$ shear flow is examined in the gyrokinetic framework in both the linear and nonlinear regimes. In the linear regime, it is shown that the eigenmode spectrum is nearly identical to hydrodynamic shear flows, with a conjugate stable mode found at every unstable wavenumber. In the nonlinear regime, turbulent saturation of the instability is examined with and without the inclusion of a driving term that prevents nonlinear flattening of the mean flow, and a scale-independent radiative damping term that suppresses the excitation of conjugate stable modes. A simple fluid model for how momentum transport and partial flattening of the mean flow scale with the driving term is constructed, from which it is shown that, except at high radiative damping, stable modes play an important role in the turbulent state and yield significantly improved quantitative predictions when compared with corresponding models neglecting stable modes.
△ Less
Submitted 24 July, 2018;
originally announced July 2018.
-
Comparison between measured and predicted turbulence frequency spectra in ITG and TEM regimes
Authors:
J. Citrin,
H. Arnichand,
J. Bernardo,
C. Bourdelle,
X. Garbet,
F. Jenko,
S. Hacquin,
M. J. Pueschel,
R. Sabot
Abstract:
The observation of distinct peaks in tokamak core reflectometry measurements - named quasi-coherent-modes (QCMs) - are identified as a signature of Trapped-Electron-Mode (TEM) turbulence [H. Arnichand et al. 2016 Plasma Phys. Control. Fusion 58 014037]. This phenomenon is investigated with detailed linear and nonlinear gyrokinetic simulations using the \gene code. A Tore-Supra density scan is stud…
▽ More
The observation of distinct peaks in tokamak core reflectometry measurements - named quasi-coherent-modes (QCMs) - are identified as a signature of Trapped-Electron-Mode (TEM) turbulence [H. Arnichand et al. 2016 Plasma Phys. Control. Fusion 58 014037]. This phenomenon is investigated with detailed linear and nonlinear gyrokinetic simulations using the \gene code. A Tore-Supra density scan is studied, which traverses through a Linear (LOC) to Saturated (SOC) Ohmic Confinement transition. The LOC and SOC phases are both simulated separately. In the LOC phase, where QCMs are observed, TEMs are robustly predicted unstable in linear studies. In the later SOC phase, where QCMs are no longer observed, ITG modes are identified. In nonlinear simulations, in the ITG (SOC) phase, a broadband spectrum is seen. In the TEM (LOC) phase, a clear emergence of a peak at the TEM frequencies is seen. This is due to reduced nonlinear frequency broadening of the underlying linear modes in the TEM regime compared with the ITG regime. A synthetic diagnostic of the nonlinearly simulated frequency spectra reproduces the features observed in the reflectometry measurements. These results support the identification of core QCMs as an experimental marker for TEM turbulence
△ Less
Submitted 12 July, 2017;
originally announced July 2017.
-
Coupling of Damped and Growing Modes in Unstable Shear Flow
Authors:
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel
Abstract:
Analysis of the saturation of the Kelvin-Helmholtz (KH) instability is undertaken to determine the extent to which the conjugate linearly stable mode plays a role. For a piecewise-continuous mean flow profile with constant shear in a fixed layer, it is shown that the stable mode is nonlinearly excited, providing an injection-scale sink of the fluctuation energy similar to what has been found for g…
▽ More
Analysis of the saturation of the Kelvin-Helmholtz (KH) instability is undertaken to determine the extent to which the conjugate linearly stable mode plays a role. For a piecewise-continuous mean flow profile with constant shear in a fixed layer, it is shown that the stable mode is nonlinearly excited, providing an injection-scale sink of the fluctuation energy similar to what has been found for gyroradius-scale drift-wave turbulence. Quantitative evaluation of the contribution of the stable mode to the energy balance at the onset of saturation shows that nonlinear energy transfer to the stable mode is as significant as energy transfer to small scales in balancing energy injected into the spectrum by the instability. The effect of the stable mode on momentum transport is quantified by expressing the Reynolds stress in terms of stable and unstable mode amplitudes at saturation, from which it is found that the stable mode can produce a sizable reduction in the momentum flux.
△ Less
Submitted 22 May, 2017; v1 submitted 19 October, 2016;
originally announced October 2016.
-
Core micro-instability analysis of JET hybrid and baseline discharges with carbon wall
Authors:
S. Moradi,
I. Pusztai,
I. Voitsekhovitch,
L. Garzotti,
C. Bourdelle,
M. J. Pueschel,
I. Lupelli,
M. Romanelli,
the JET-EFDA contributors
Abstract:
The core micro-instability characteristics of hybrid and baseline plasmas in a selected set of JET plasmas with carbon wall are investigated through local linear and non-linear and global linear gyro-kinetic simulations with the GYRO code [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. In particular, we study the role of plasma pressure on the micro-instabilities, and scan the pa…
▽ More
The core micro-instability characteristics of hybrid and baseline plasmas in a selected set of JET plasmas with carbon wall are investigated through local linear and non-linear and global linear gyro-kinetic simulations with the GYRO code [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. In particular, we study the role of plasma pressure on the micro-instabilities, and scan the parameter space for the important plasma parameters responsible for the onset and stabilization of the modes under experimental conditions. We find that a good core confinement due to strong stabilization of the micro-turbulence driven transport can be expected in the hybrid plasmas due to the stabilizing effect of the fast ion pressure that is more effective at the low magnetic shear of the hybrid discharges. While parallel velocity gradient destabilization is important for the inner core, at outer radii the hybrid plasmas may benefit from a strong quench of the turbulence transport by $\mathbf{E}\times\mathbf{B}$ rotation shear.
△ Less
Submitted 18 November, 2014;
originally announced November 2014.
-
Electromagnetic stabilization of tokamak microturbulence in a high-$β$ regime
Authors:
J. Citrin,
J. Garcia,
T. Goerler,
F. Jenko,
P. Mantica,
D. Told,
C. Bourdelle,
D. R. Hatch,
G. M. D. Hogeweij,
T. Johnson,
M. J. Pueschel,
M. Schneider
Abstract:
The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$β$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive…
▽ More
The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$β$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive linear and nonlinear gyrokinetic simulations are then carried out for the high-$β$ JET hybrid discharge 75225, at two separate locations at inner and outer radii. It is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment. The enhancement of this effect by suprathermal pressure also remains significant. It is also found that flow shear stabilization is not effective at the inner radii. However, at outer radii the situation is reversed. Electromagnetic stabilization is negligible while the flow shear stabilization is significant. These results constitute the high-$β$ generalization of comparable observations found at low-$β$ at JET. This is encouraging for the extrapolation of electromagnetic ITG stabilization to future devices. An estimation of the impact of this effect on the ITER hybrid scenario leads to a 20% fusion power improvement.
△ Less
Submitted 5 September, 2014;
originally announced September 2014.
-
Nonlinear stabilization of tokamak microturbulence by fast ions
Authors:
J. Citrin,
C. Bourdelle,
J. Garcia,
J. W. Haverkort,
G. M. D. Hogeweij,
F. Jenko,
T. Johnson,
P. Mantica,
M. J. Pueschel,
D. Told,
JET-EFDA contributors
Abstract:
Nonlinear electromagnetic stabilization by suprathermal pressure gradients found in specific regimes is shown to be a key factor in reducing tokamak microturbulence, augmenting significantly the thermal pressure electromagnetic stabilization. Based on nonlinear gyrokinetic simulations investigating a set of ion heat transport experiments on the JET tokamak, described by Mantica et al. [Phys. Rev.…
▽ More
Nonlinear electromagnetic stabilization by suprathermal pressure gradients found in specific regimes is shown to be a key factor in reducing tokamak microturbulence, augmenting significantly the thermal pressure electromagnetic stabilization. Based on nonlinear gyrokinetic simulations investigating a set of ion heat transport experiments on the JET tokamak, described by Mantica et al. [Phys. Rev. Lett. 107 135004 (2011)], this result explains the experimentally observed ion heat flux and stiffness reduction. These findings are expected to improve the extrapolation of advanced tokamak scenarios to reactor relevant regimes.
△ Less
Submitted 26 June, 2013;
originally announced June 2013.
-
Ion temperature profile stiffness: non-linear gyrokinetic simulations and comparison with experiment
Authors:
J. Citrin,
F. Jenko,
P. Mantica,
D. Told,
C. Bourdelle,
R. Dumont,
J. Garcia,
J. W. Haverkort,
G. M. D. Hogeweij,
T. Johnson,
M. J. Pueschel,
JET-EFDA contributors
Abstract:
Recent experimental observations at JET show evidence of reduced ion temperature profile stiffness, hypothesised to be due to concomitant low magnetic shear (s) and significant toroidal rotational flow shear. Non-linear gyrokinetic simulations are performed, aiming to investigate the physical mechanism behind the observations. A comprehensive set of simulations are carried out, comparing the impac…
▽ More
Recent experimental observations at JET show evidence of reduced ion temperature profile stiffness, hypothesised to be due to concomitant low magnetic shear (s) and significant toroidal rotational flow shear. Non-linear gyrokinetic simulations are performed, aiming to investigate the physical mechanism behind the observations. A comprehensive set of simulations are carried out, comparing the impact on the ion heat flux of various parameters that differ within the data-set. These parameters include q, s, rotation, effect of rotation on the magnetohydrodynamic (MHD) equilibrium, R/L_n, beta_e, Z_eff, and the fast particle content. The effect of toroidal flow shear itself is not predicted by the simulations to lead to a significant reduction in ion heat flux, due both to an insufficient magnitude of flow shear and significant parallel velocity gradient destabilisation. It is however found that non-linear electromagnetic effects due to both thermal and fast-particle pressure gradients, even at low beta_e, can significantly reduce the profile stiffness. A total of five discharges are examined, at both inner and outer radii. For all cases studied, the simulated and experimental ion heat flux values agree within reasonable variations of input parameters around the experimental uncertainties.
△ Less
Submitted 16 July, 2013; v1 submitted 9 March, 2013;
originally announced March 2013.
-
Gyrokinetic Studies of Microinstabilities in the RFP
Authors:
Daniel Carmody,
M. J. Pueschel,
P. W. Terry
Abstract:
An analytic equilibrium, the Toroidal Bessel Function Model, is used in conjunction with the gyrokinetic code GYRO to investigate the nature of microinstabilities in a reversed field pinch (RFP) plasma. The effect of the normalized electron plasma pressure (β) on the characteristics of the microinstabilities is studied. A transition between an ion temperature gradient (ITG) driven mode and a micro…
▽ More
An analytic equilibrium, the Toroidal Bessel Function Model, is used in conjunction with the gyrokinetic code GYRO to investigate the nature of microinstabilities in a reversed field pinch (RFP) plasma. The effect of the normalized electron plasma pressure (β) on the characteristics of the microinstabilities is studied. A transition between an ion temperature gradient (ITG) driven mode and a microtearing mode as the dominant instability is found to occur at a β value of approximately 4.5%. Suppression of the ITG mode occurs as in the tokamak, through coupling to shear Alfven waves, with a critical β for stability higher than its tokamak equivalent due to a shorter parallel connection length. There is a steep dependence of the microtearing growth rate on temperature gradient suggesting high profile stiffness. There is evidence for a collisionless microtearing mode. The properties of this mode are investigated, and it is found that curvature drift plays an important role in the instability.
△ Less
Submitted 19 January, 2013;
originally announced January 2013.