-
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.
-
Simplified and Flexible Coils for Stellarators using Single-Stage Optimization
Authors:
R. Jorge,
A. Giuliani,
J. Loizu
Abstract:
Single-stage optimization, also known as combined plasma-coil algorithms or direct coil optimization, has recently emerged as a possible method to expedite the design of stellarator devices by including, in a single step, confinement, stability, and engineering constraints. In this work, we show how such frameworks allow us to find new designs in a streamlined manner, yielding a broad range of new…
▽ More
Single-stage optimization, also known as combined plasma-coil algorithms or direct coil optimization, has recently emerged as a possible method to expedite the design of stellarator devices by including, in a single step, confinement, stability, and engineering constraints. In this work, we show how such frameworks allow us to find new designs in a streamlined manner, yielding a broad range of new configurations. Examples are shown for stellarators with a small number of coils and quasisymmetric stellarators with only one to three coils per half field period, with external trim coils, helical coils, and a single set of coils generating both a quasi-axisymmetric and a quasi-helical equilibrium.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Electron root optimisation for stellarator reactor designs
Authors:
E. Lascas Neto,
R. Jorge,
C. D. Beidler,
J. Lion
Abstract:
In this work, we propose a method of optimising stellarator devices to favour the presence of an electron root solution of the radial electric field. Such a solution can help avoid heavy impurity accumulation, improve neoclassical thermal ion confinement and helium ash exhaust and possibly reduce turbulence. This study shows that an optimisation for such a root is possible in quasi-isodynamic stel…
▽ More
In this work, we propose a method of optimising stellarator devices to favour the presence of an electron root solution of the radial electric field. Such a solution can help avoid heavy impurity accumulation, improve neoclassical thermal ion confinement and helium ash exhaust and possibly reduce turbulence. This study shows that an optimisation for such a root is possible in quasi-isodynamic stellarators. Examples are shown for both vacuum and finite plasma pressure configurations.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
An asymptotic Grad-Shafranov equation for quasisymmetric stellarators
Authors:
Nikita Nikulsin,
Wrick Sengupta,
Rogerio Jorge,
Amitava Bhattacharjee
Abstract:
A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $β$ of the order of the ratio of induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for…
▽ More
A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $β$ of the order of the ratio of induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad-Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions are the first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement.
△ Less
Submitted 21 June, 2024; v1 submitted 12 March, 2024;
originally announced March 2024.
-
Tokamak to Stellarator Conversion using Permanent Magnets
Authors:
M. Madeira,
R. Jorge
Abstract:
With the advances in the optimization of magnetic field equilibria, stellarators have become a serious alternative to the tokamak, bringing this concept to the forefront of the pursuit of fusion energy. In order to be successful in experimentally demonstrating the viability of optimized stellarators, we must overcome any potential hurdles in the construction of its electromagnetic coils. Finding c…
▽ More
With the advances in the optimization of magnetic field equilibria, stellarators have become a serious alternative to the tokamak, bringing this concept to the forefront of the pursuit of fusion energy. In order to be successful in experimentally demonstrating the viability of optimized stellarators, we must overcome any potential hurdles in the construction of its electromagnetic coils. Finding cost-effective ways of increasing the number of operating optimized stellarators could be key in cementing this magnetic confinement concept as a contender for a reactor. In this work, an alternative to modular coils, permanent magnets, are studied and are shown to enable the possibility of converting a tokamak into a stellarator. This is then applied to the case of ISTTOK tokamak, where an engineering design study is conducted.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta
Authors:
W. Sengupta,
E. Rodriguez,
R. Jorge,
M. Landreman A. Bhattacharjee
Abstract:
A systematic theory of the asymptotic expansion of the magnetohydrodynamic (MHD) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared to the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in wh…
▽ More
A systematic theory of the asymptotic expansion of the magnetohydrodynamic (MHD) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared to the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHD system closely mimics the vacuum system. Thus, a unified treatment of MHD fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free, and MHD equilibrium in different geometries are presented.
△ Less
Submitted 26 February, 2024;
originally announced February 2024.
-
Energetic Particle Tracing in Optimized Quasisymmetric Stellarator Equilibria
Authors:
P. A. Figueiredo,
R. Jorge,
J. Ferreira,
P. Rodrigues
Abstract:
Recent developments in the design of magnetic confinement fusion devices have allowed the construction of exceptionally optimized stellarator configurations. The near-axis expansion in particular has proven to enable the construction of magnetic configurations with good confinement properties while taking only a fraction of the usual computation time to generate optimized magnetic equilibria. Howe…
▽ More
Recent developments in the design of magnetic confinement fusion devices have allowed the construction of exceptionally optimized stellarator configurations. The near-axis expansion in particular has proven to enable the construction of magnetic configurations with good confinement properties while taking only a fraction of the usual computation time to generate optimized magnetic equilibria. However, not much is known about the overall features of fast-particle orbits computed in such analytical, yet simplified, equilibria when compared to those originating from accurate equilibrium solutions. This work aims to assess and demonstrate the potential of the near-axis expansion to provide accurate information on particle orbits and to compute loss fractions in moderate to high aspect ratios. The configurations used here are all scaled to fusion-relevant parameters and approximate quasisymmetry in various degrees. This allows us to understand how deviations from quasisymmetry affect particle orbits and what are their effects on the estimation of the loss fraction. Guiding-center trajectories of fusion-born alpha particles are traced using gyronimo and SIMPLE codes under the NEAT framework, showing good numerical agreement. Discrepancies between near-axis and MHD fields have minor effects on passing particles but significant effects on trapped particles, especially in quasihelically symmetric magnetic fields. Effective expressions were found for estimating orbit widths and passing-trapped separatrix in quasisymmetric near-axis fields. Loss fractions agree in the prompt losses regime but diverge afterward.
△ Less
Submitted 15 April, 2024; v1 submitted 13 November, 2023;
originally announced November 2023.
-
Optimization of Nonlinear Turbulence in Stellarators
Authors:
Patrick Kim,
Stefan Buller,
Rory Conlin,
William Dorland,
Daniel W. Dudt,
Rahul Gaur,
Rogerio Jorge,
Egemen Kolemen,
Matt Landreman,
Noah R. Mandell,
Dario Panici
Abstract:
We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC. Since using GX allows for fast nonlinear simulations, we directly optimize for r…
▽ More
We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC. Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasisymmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.
△ Less
Submitted 22 February, 2024; v1 submitted 28 October, 2023;
originally announced October 2023.
-
Single-Stage Stellarator Optimization: Combining Coils with Fixed Boundary Equilibria
Authors:
R. Jorge,
A. Goodman,
M. Landreman,
J. Rodrigues,
F. Wechsung
Abstract:
We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using…
▽ More
We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using a quadratic flux term in the objective function. Four use cases are presented to demonstrate the effectiveness of the approach, including simple and complex stellarator geometries. As shown here, this method outperforms previous 2-stage approaches, achieving smaller plasma objective function values when coils are taken into account.
△ Less
Submitted 25 April, 2023; v1 submitted 21 February, 2023;
originally announced February 2023.
-
Direct Microstability Optimization of Stellarator Devices
Authors:
R. Jorge,
W. Dorland,
P. Kim,
M. Landreman,
N. R. Mandell,
G. Merlo,
T. Qian
Abstract:
Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks in stellarators. In this letter, we show that a significant decrease in the turbulent heat flux can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasi-linear heat flux at each step of…
▽ More
Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks in stellarators. In this letter, we show that a significant decrease in the turbulent heat flux can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasi-linear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and turbulent transport.
△ Less
Submitted 23 January, 2023;
originally announced January 2023.
-
Constructing precisely quasi-isodynamic magnetic fields
Authors:
Alan Goodman,
Katia Camacho Mata,
Sophia A Henneberg,
Rogerio Jorge,
Matt Landreman,
Gabriel Plunk,
Hakan Smith,
Ralf Mackenbach,
Per Helander
Abstract:
We present a novel method for numerically finding quasi-isodynamic stellarator magnetic fields with excellent fast-particle confinement and extremely small neoclassical transport. The method works particularly well in configurations with only one field period. We examine the properties of these newfound quasi-isodynamic configurations, including their bootstrap currents, particle confinement, and…
▽ More
We present a novel method for numerically finding quasi-isodynamic stellarator magnetic fields with excellent fast-particle confinement and extremely small neoclassical transport. The method works particularly well in configurations with only one field period. We examine the properties of these newfound quasi-isodynamic configurations, including their bootstrap currents, particle confinement, and available energy for trapped-electron driven turbulence, as well as the degree to which they change when a finite pressure profile is added. We finally discuss the differences between the magnetic axes of the optimized solutions and their respective initial conditions, and conclude with the prospects for future quasi-isodynamic optimization.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
Direct construction of stellarator-symmetric quasi-isodynamic magnetic configurations
Authors:
Katia Camacho Mata,
Gabriel G. Plunk,
Rogerio Jorge
Abstract:
We develop the formalism of the first order near-axis expansion of the MHD equilibrium equations described in Garren & Boozer (1991), Plunk et al. (2019) and Plunk et al. (2021), for the case of a quasi-isodynamic, N-field period, stellarator symmetric, single-well magnetic field equilibrium. The importance of the magnetic axis shape is investigated, and we conclude that control of the curvature a…
▽ More
We develop the formalism of the first order near-axis expansion of the MHD equilibrium equations described in Garren & Boozer (1991), Plunk et al. (2019) and Plunk et al. (2021), for the case of a quasi-isodynamic, N-field period, stellarator symmetric, single-well magnetic field equilibrium. The importance of the magnetic axis shape is investigated, and we conclude that control of the curvature and torsion is crucial to obtain omnigenous configurations with finite aspect ratio and low effective ripple, especially for a higher number of field periods. For this reason a method is derived to construct classes of axis shapes with favourable curvature and torsion. Solutions are presented, including a three-field-period configuration constructed at an aspect ratio of A=20, with a maximum elongation of e=3.2 and an effective ripple under 1%, which demonstrates that high elongation is not a necessary feature of QI stellarators.
△ Less
Submitted 16 June, 2022;
originally announced June 2022.
-
A single-field-period quasi-isodynamic stellarator
Authors:
R. Jorge,
G. G. Plunk,
M. Drevlak,
M. Landreman,
J. -F. Lobsien,
K. Camacho Mata,
P. Helander
Abstract:
A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configuration…
▽ More
A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configurations defined by a set of adjustable magnetic-field parameters. This optimization, completed in 30 seconds on a single cpu core using the SIMSOPT code, yields a solution with excellent confinement, as measured by the conventional figure of merit for neoclassical transport, effective ripple, at a modest aspect ratio of eight. The optimization parameters that led to this configuration are described, its confinement properties are assessed, and a set of magnetic-field coils is found. The resulting transport at low collisionality is much smaller than that of W7-X, and the device needs significantly fewer coils thanks to the reduced number of field periods.
△ Less
Submitted 1 September, 2022; v1 submitted 11 May, 2022;
originally announced May 2022.
-
Development of Advanced Linearized Gyrokinetic Collision Operators Using a Moment Approach
Authors:
B. J. Frei,
J. Ball,
A. C. D. Hoffmann,
R. Jorge,
P. Ricci,
L. Stenger
Abstract:
The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge R. et al., J. Plasma Phys. 85, 905850604 (2019)) and of the widely-used linearized GK Sugama collision operator (Sugama H. et al., Phys. Plasmas 16, 112503 (2009)) is reported. An approach based on a Hermite-Laguerre moment expansion of the perturbed gyrocenter distribution…
▽ More
The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge R. et al., J. Plasma Phys. 85, 905850604 (2019)) and of the widely-used linearized GK Sugama collision operator (Sugama H. et al., Phys. Plasmas 16, 112503 (2009)) is reported. An approach based on a Hermite-Laguerre moment expansion of the perturbed gyrocenter distribution function is used, referred to as gyro-moment expansion. This approach allows considering arbitrary perpendicular wavenumber and expressing the two linearized GK operators as a linear combination of gyro-moments where the expansion coefficients are given by closed analytical expressions that depend on the perpendicular wavenumber and on the temperature and mass ratios of the colliding species. The drift-kinetic (DK) limits of the GK linearized Coulomb and Sugama operators are also obtained. Comparisons between the gyro-moment approach with the GK continuum code GENE are reported focusing on the ion-temperature-gradient (ITG) instability and zonal flow (ZF) damping, finding an excellent agreement. In particular, we demonstrate that the GK linearized Sugama yields a stronger collisional damping of the ZF residual compared to the GK linearized Coulomb. Finally, we show that the numerical efficiency of the gyro-moment approach increases with collisionality, a desired property for boundary plasma applications.
△ Less
Submitted 13 February, 2022; v1 submitted 23 April, 2021;
originally announced April 2021.
-
Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators
Authors:
R. Jorge,
M. Landreman
Abstract:
The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using…
▽ More
The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations.
△ Less
Submitted 24 February, 2021;
originally announced February 2021.
-
The On-Axis Magnetic Well and Mercier's Criterion for Arbitrary Stellarator Geometries
Authors:
Patrick Kim,
Rogerio Jorge,
William Dorland
Abstract:
A simplified analytical form of the on-axis magnetic well and Mercier's criterion for interchange instabilities for arbitrary three-dimensional magnetic field geometries is derived. For this purpose, a near-axis expansion based on a direct coordinate approach is used by expressing the toroidal magnetic flux in terms of powers of the radial distance to the magnetic axis. For the first time, the mag…
▽ More
A simplified analytical form of the on-axis magnetic well and Mercier's criterion for interchange instabilities for arbitrary three-dimensional magnetic field geometries is derived. For this purpose, a near-axis expansion based on a direct coordinate approach is used by expressing the toroidal magnetic flux in terms of powers of the radial distance to the magnetic axis. For the first time, the magnetic well and Mercier's criterion are then written as a one-dimensional integral with respect to the axis arclength. When compared with the original work of Mercier, the derivation here is presented using modern notation and in a more streamlined manner that highlights essential steps. Finally, these expressions are verified numerically using several quasisymmetric and non-quasisymmetric stellarator configurations including Wendelstein 7-X.
△ Less
Submitted 7 April, 2021; v1 submitted 14 November, 2020;
originally announced November 2020.
-
Four-dimensional drift-kinetic model for scrape-off layer plasmas
Authors:
L. M. Perrone,
R. Jorge,
P. Ricci
Abstract:
A four-dimensional plasma model able to describe the scrape-off layer region of tokamak devices at arbitrary collisionality is derived in the drift-reduced limit. The basis of the model is provided by a drift-kinetic equation that retains the full non-linear Coulomb collision operator and describes arbitrarily far from equilibrium distribution functions. By expanding the dependence of distribution…
▽ More
A four-dimensional plasma model able to describe the scrape-off layer region of tokamak devices at arbitrary collisionality is derived in the drift-reduced limit. The basis of the model is provided by a drift-kinetic equation that retains the full non-linear Coulomb collision operator and describes arbitrarily far from equilibrium distribution functions. By expanding the dependence of distribution function over the perpendicular velocity in a Laguerre polynomial basis and integrating over the perpendicular velocity, a set of four-dimensional moment equations for the expansion coefficients of the distribution function is obtained. The Coulomb collision operator, as well as Poisson's equation, are evaluated explicitly in terms of perpendicular velocity moments of the distribution function.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
The Use of Near-Axis Magnetic Fields for Stellarator Turbulence Simulations
Authors:
R. Jorge,
M. Landreman
Abstract:
The design of turbulence optimized stellarators has so far relied on three-dimensional equilibrium codes such as VMEC in order to find the minimum of a given objective function. In this work, we propose a complimentary approach based on the near-axis expansion to compute the geometry parameters of neoclassicaly optimized stellarators used in turbulence studies. As shown here, the near-axis expansi…
▽ More
The design of turbulence optimized stellarators has so far relied on three-dimensional equilibrium codes such as VMEC in order to find the minimum of a given objective function. In this work, we propose a complimentary approach based on the near-axis expansion to compute the geometry parameters of neoclassicaly optimized stellarators used in turbulence studies. As shown here, the near-axis expansion can be a reasonable approximation of the geometric parameters relevant for turbulence and stability simulations of the core of existing optimized stellarator designs. In particular, we examine the geometry coefficients that appear in the gyrokinetic equation, the drift-reduced fluid equations and the ideal ballooning equation. This approach may allow for the development of new stellarator optimization techniques significantly faster than conventional methods.
△ Less
Submitted 20 August, 2020;
originally announced August 2020.
-
Magnetic well and Mercier stability of stellarators near the magnetic axis
Authors:
Matt Landreman,
Rogerio Jorge
Abstract:
We have recently demonstrated that by expanding in small distance from the magnetic axis compared to the major radius, stellarator shapes with low neoclassical transport can be generated efficiently. To extend the utility of this new design approach, here we evaluate measures of magnetohydrodynamic interchange stability within the same expansion. In particular, we evaluate magnetic well, Mercier's…
▽ More
We have recently demonstrated that by expanding in small distance from the magnetic axis compared to the major radius, stellarator shapes with low neoclassical transport can be generated efficiently. To extend the utility of this new design approach, here we evaluate measures of magnetohydrodynamic interchange stability within the same expansion. In particular, we evaluate magnetic well, Mercier's criterion, and resistive interchange stability near a magnetic axis of arbitrary shape. In contrast to previous work on interchange stability near the magnetic axis, which used an expansion of the flux coordinates, here we use the inverse expansion in which the flux coordinates are the independent variables. Reduced expressions are presented for the magnetic well and stability criterion in the case of quasisymmetry. The analytic results are shown to agree with calculations from the VMEC equilibrium code. Finally, we show that near the axis, Glasser, Greene, & Johnson's stability criterion for resistive modes approximately coincides with Mercier's ideal condition.
△ Less
Submitted 2 September, 2020; v1 submitted 26 June, 2020;
originally announced June 2020.
-
Construction of Quasisymmetric Stellarators Using a Direct Coordinate Approach
Authors:
R. Jorge,
W. Sengupta,
M. Landreman
Abstract:
Optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic fie…
▽ More
Optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface.
△ Less
Submitted 13 March, 2020;
originally announced March 2020.
-
Near-Axis Expansion of Stellarator Equilibrium at Arbitrary Order in the Distance to the Axis
Authors:
R. Jorge,
W. Sengupta,
M. Landreman
Abstract:
A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight on the character of high-aspect-ratio stellarator shapes, such as…
▽ More
A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight on the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta-limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect-ratio of the ideal MHD equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet-Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren-Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest order solution and the W7-X magnetic field is performed.
△ Less
Submitted 28 January, 2020; v1 submitted 6 November, 2019;
originally announced November 2019.
-
Nonlinear Gyrokinetic Coulomb Collision Operator
Authors:
R. Jorge,
B. J. Frei,
P. Ricci
Abstract:
A gyrokinetic Coulomb collision operator is derived, which is particularly useful to describe the plasma dynamics at the periphery region of magnetic confinement fusion devices. The derived operator is able to describe collisions occurring in distribution functions arbitrarily far from equilibrium with variations on spatial scales at and below the particle Larmor radius. A multipole expansion of t…
▽ More
A gyrokinetic Coulomb collision operator is derived, which is particularly useful to describe the plasma dynamics at the periphery region of magnetic confinement fusion devices. The derived operator is able to describe collisions occurring in distribution functions arbitrarily far from equilibrium with variations on spatial scales at and below the particle Larmor radius. A multipole expansion of the Rosenbluth potentials is used in order to derive the dependence of the full Coulomb collision operator on the particle gyroangle. The full Coulomb collision operator is then expressed in gyrocentre phase-space coordinates, and a closed formula for its gyroaverage in terms of the moments of the gyrocenter distribution function in a form ready to be numerically implemented is provided. Furthermore, the collision operator is projected onto a Hermite-Laguerre velocity space polynomial basis and expansions in the small electron-to-ion mass ratio are provided.
△ Less
Submitted 4 November, 2019; v1 submitted 22 May, 2019;
originally announced June 2019.
-
A gyrokinetic model for the plasma periphery of tokamak devices
Authors:
B. J. Frei,
R. Jorge,
P. Ricci
Abstract:
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order description of the singl…
▽ More
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts, assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behavior of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function and includes a non-linear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a velocity-space Hermite-Laguerre polynomial basis and the gyrokinetic equation is projected onto the same basis, obtaining the spatial and temporal evolution of the Hermite-Laguerre expansion coefficients. The Hermite-Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwell's equations derived from a variational principle, with the velocity integrals of the gyroaveraged distribution function explicitly evaluated.
△ Less
Submitted 22 April, 2020; v1 submitted 15 April, 2019;
originally announced April 2019.
-
Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality
Authors:
R. Jorge,
P. Ricci,
S. Brunner,
S. Gamba,
V. Konovets,
N. F. Loureiro,
L. M. Perrone,
N. Teixeira
Abstract:
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength…
▽ More
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard-Bernstein and Dougherty collision operators, and the electron-ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode.
△ Less
Submitted 5 March, 2019; v1 submitted 26 November, 2018;
originally announced November 2018.
-
Theory of the Drift-Wave Instability at Arbitrary Collisionality
Authors:
R. Jorge,
P. Ricci,
N. F. Loureiro
Abstract:
A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis, to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show…
▽ More
A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis, to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show that our framework allows retrieving established collisional and collisionless limits. At the intermediate collisionalities relevant for present and future magnetic nuclear fusion devices, deviations with respect to collision operators used in state-of-the-art turbulence simulation codes show the need for retaining the full Coulomb operator in order to obtain both the correct instability growth rate and eigenmode spectrum, which, for example, may significantly impact quantitative predictions of transport. The exponential convergence of the spectral representation that we propose makes the representation of the velocity space dependence, including the full collision operator, more efficient than standard finite difference methods.
△ Less
Submitted 18 September, 2018; v1 submitted 18 June, 2018;
originally announced June 2018.
-
Electroluminescence TPCs at the thermal diffusion limit
Authors:
C. A. O. Henriques,
C. M. B. Monteiro,
D. González-Díaz,
C. D. R Azevedo,
E. D. C. Freitas,
R. D. P. Mano,
M. R. Jorge,
A. F. M. Fernandes,
J. J. Gómez-Cadenas,
L. M. P. Fernandes,
C. Adams,
V. Álvarez,
L. Arazi,
K. Bailey,
F. Ballester,
J. M. Benlloch-Rodríguez,
F. I. G. M. Borges,
A. Botas,
S. Cárcel,
J. V. Carrión,
S. Cebrián,
C. A. N. Conde,
J. Díaz,
M. Diesburg,
J. Escada
, et al. (56 additional authors not shown)
Abstract:
The NEXT experiment aims at searching for the hypothetical neutrinoless double-beta decay from the ${}^{136}$Xe isotope using a high-purity xenon TPC. Efficient discrimination of the events through pattern recognition of the topology of primary ionisation tracks is a major requirement for the experiment. However, it is limited by the diffusion of electrons. It is known that the addition of a small…
▽ More
The NEXT experiment aims at searching for the hypothetical neutrinoless double-beta decay from the ${}^{136}$Xe isotope using a high-purity xenon TPC. Efficient discrimination of the events through pattern recognition of the topology of primary ionisation tracks is a major requirement for the experiment. However, it is limited by the diffusion of electrons. It is known that the addition of a small fraction of a molecular gas to xenon reduces electron diffusion. On the other hand, the electroluminescence (EL) yield drops and the achievable energy resolution may be compromised. We have studied the effect of adding several molecular gases to xenon (CO${}_{2}$, CH${}_{4}$ and CF${}_{4}$) on the EL yield and energy resolution obtained in a small prototype of driftless gas proportional scintillation counter. We have compared our results on the scintillation characteristics (EL yield and energy resolution) with a microscopic simulation, obtaining the diffusion coefficients in those conditions as well. Accordingly, electron diffusion may be reduced from about 10 mm/$\sqrt{\mathrm{m}}$ for pure xenon down to 2.5 mm/$\sqrt{\mathrm{m}}$ using additive concentrations of about 0.05%, 0.2% and 0.02% for CO${}_{2}$, CH${}_{4}$ and CF${}_{4}$, respectively. Our results show that CF${}_{4}$ admixtures present the highest EL yield in those conditions, but very poor energy resolution as a result of huge fluctuations observed in the EL formation. CH${}_{4}$ presents the best energy resolution despite the EL yield being the lowest. The results obtained with xenon admixtures are extrapolated to the operational conditions of the NEXT-100 TPC. CO${}_{2}$ and CH${}_{4}$ show potential as molecular additives in a large xenon TPC, CH${}_{4}$ showing the best performance and stability to be used in the NEXT-100 TPC, with an extrapolated energy resolution of 0.4% at 2.45 MeV for concentrations below 0.4%.
△ Less
Submitted 30 October, 2018; v1 submitted 15 June, 2018;
originally announced June 2018.
-
A Drift-Kinetic Analytical Model for SOL Plasma Dynamics at Arbitrary Collisionality
Authors:
R. Jorge,
P. Ricci,
N. F. Loureiro
Abstract:
A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition…
▽ More
A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously-derived fluid models.
△ Less
Submitted 27 November, 2017; v1 submitted 4 September, 2017;
originally announced September 2017.
-
Secondary scintillation yield of Xenon with sub-percent levels of CO2 additive: efficiently reducing electron diffusion in HPXe optical TPCs for rare-event detection
Authors:
C. A. O. Henriques,
E. D. C. Freitas,
C. D. R. Azevedo,
D. González-Díaz,
R. D. P. Mano,
M. R. Jorge,
L. M. P. Fernandes,
C. M. B. Monteiro,
J. J. Gómez-Cadenas,
V. Álvarez,
J. M. Benlloch-Rodríguez,
F. I. G. M. Borges,
A. Botas,
S. Cárcel,
J. V. Carrión,
S. Cebrián,
C. A. N. Conde,
J. Díaz,
M. Diesburg,
J. Escada,
R. Esteve,
R. Felkai,
P. Ferrario,
A. L. Ferreira,
A. Goldschmidt
, et al. (45 additional authors not shown)
Abstract:
We have measured the electroluminescence (EL) yield of Xe-CO2 mixtures, with sub-percent CO2 concentrations. We demonstrate that the EL production is still high in these mixtures, 70% and 35% relative to that produced in pure xenon, for CO2 concentrations around 0.05% and 0.1%, respectively. The contribution of the statistical fluctuations in EL production to the energy resolution increases with i…
▽ More
We have measured the electroluminescence (EL) yield of Xe-CO2 mixtures, with sub-percent CO2 concentrations. We demonstrate that the EL production is still high in these mixtures, 70% and 35% relative to that produced in pure xenon, for CO2 concentrations around 0.05% and 0.1%, respectively. The contribution of the statistical fluctuations in EL production to the energy resolution increases with increasing CO2 concentration and, for our gas proportional scintillation counter, it is smaller than the contribution of the Fano factor for concentrations below 0.1% CO2. Xe-CO2 mixtures are important alternatives to pure xenon in TPCs based on EL signal amplification with applications in the important field of rare event detection such as directional dark matter, double electron capture and double beta decay detection. The addition of CO2 to pure xenon at the level of 0.05-0.1% can reduce significantly the scale of electron diffusion from 10 mm/sqrt(m) to 2.5 mm/sqrt(m), with high impact on the HPXe TPC discrimination efficiency of the events through pattern recognition of the topology of primary ionisation trails.
△ Less
Submitted 12 April, 2017; v1 submitted 5 April, 2017;
originally announced April 2017.
-
Plasma Turbulence in the Scrape-off Layer of the ISTTOK Tokamak
Authors:
Rogerio Jorge,
Paolo Ricci,
Federico D. Halpern,
Nuno F. Loureiro,
Carlos Silva
Abstract:
The properties of plasma turbulence in a poloidally limited scrape-off layer (SOL) are addressed, with focus on ISTTOK, a large aspect ratio tokamak with a circular cross section. Theoretical investigations based on the drift-reduced Braginskii equations are carried out through linear calculations and non-linear simulations, in two- and three-dimensional geometries. The linear instabilities drivin…
▽ More
The properties of plasma turbulence in a poloidally limited scrape-off layer (SOL) are addressed, with focus on ISTTOK, a large aspect ratio tokamak with a circular cross section. Theoretical investigations based on the drift-reduced Braginskii equations are carried out through linear calculations and non-linear simulations, in two- and three-dimensional geometries. The linear instabilities driving turbulence and the mechanisms that set the amplitude of turbulence as well as the SOL width are identified. A clear asymmetry is shown to exist between the low-field and the high-field sides of the machine. A comparison between experimental measurements and simulation results is presented.
△ Less
Submitted 2 October, 2016; v1 submitted 30 June, 2016;
originally announced June 2016.
-
First in-beam studies of a Resistive-Plate WELL gaseous multiplier
Authors:
S. Bressler,
L. Moleri,
M. Pitt,
S. Kudella,
D. C. R. Azevedo,
F. D. Amaro,
M. R. Jorge,
J. M. F. dos Santos,
J. F. C. A. Veloso,
H. Natal da Luz,
L. Arazi,
E. Olivieri,
A. Breskin
Abstract:
We present the results of the first in-beam studies of a medium size (10$\times$10 cm$^2$) Resistive-Plate WELL (RPWELL): a single-sided THGEM coupled to a pad anode through a resistive layer of high bulk resistivity ($\sim$10$^9 Ω$cm). The 6.2~mm thick (excluding readout electronics) single-stage detector was studied with 150~GeV muons and pions. Signals were recorded from 1$\times$1 cm$^2$ squar…
▽ More
We present the results of the first in-beam studies of a medium size (10$\times$10 cm$^2$) Resistive-Plate WELL (RPWELL): a single-sided THGEM coupled to a pad anode through a resistive layer of high bulk resistivity ($\sim$10$^9 Ω$cm). The 6.2~mm thick (excluding readout electronics) single-stage detector was studied with 150~GeV muons and pions. Signals were recorded from 1$\times$1 cm$^2$ square copper pads with APV25-SRS readout electronics. The single-element detector was operated in Ne\(5% $\mathrm{CH_{4}}$) at a gas gain of a few times 10$^4$, reaching 99$\%$ detection efficiency at average pad multiplicity of $\sim$1.2. Operation at particle fluxes up to $\sim$10$^4$ Hz/cm$^2$ resulted in $\sim$23$\%$ gain drop leading to $\sim$5$\%$ efficiency loss. The striking feature was the discharge-free operation, also in intense pion beams. These results pave the way towards robust, efficient large-scale detectors for applications requiring economic solutions at moderate spatial and energy resolutions.
△ Less
Submitted 20 January, 2016; v1 submitted 11 October, 2015;
originally announced October 2015.
-
Nonlinear Acoustics -- Perturbation Theory and Webster's Equation
Authors:
Rogério Jorge
Abstract:
Webster's horn equation (1919) offers a one-dimensional approximation for low-frequency sound waves along a rigid tube with a variable cross-sectional area. It can be thought as a wave equation with a source term that takes into account the nonlinear geometry of the tube. In this document we derive this equation using a simplified fluid model of an ideal gas. By a simple change of variables, we co…
▽ More
Webster's horn equation (1919) offers a one-dimensional approximation for low-frequency sound waves along a rigid tube with a variable cross-sectional area. It can be thought as a wave equation with a source term that takes into account the nonlinear geometry of the tube. In this document we derive this equation using a simplified fluid model of an ideal gas. By a simple change of variables, we convert it to a Schrödinger equation and use the well-known variational and perturbative methods to seek perturbative solutions. As an example, we apply these methods to the Gabriel's Horn geometry, deriving the first order corrections to the linear frequency. An algorithm to the harmonic modes in any order for a general horn geometry is derived.
△ Less
Submitted 17 November, 2013;
originally announced November 2013.
