-
Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach
Authors:
Erol Balkovic,
Joaquim Loizu,
Jonathan P. Graves,
Yi-Min Huang,
Christopher B. Smiet
Abstract:
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson e…
▽ More
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson et al. 2012), an equilibrium solver based on the variational principle of the Multi-Region relaxed MHD, featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model $J_\textrm{bs} = C \nabla p$ to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter $Δ'$ and the bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHD code (Huang & Bhattacharjee 2016) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical Modified Rutherford Equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and neoclassical tearing modes, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Efficient single-stage optimization of islands in finite-$β$ stellarator equilibria
Authors:
Christopher Berg Smiet,
Joaquim Loizu,
Erol Balkovic,
Antoine Baillod
Abstract:
We present the first single-stage optimization of islands in finite-$β$ stellarator equilibria. Stellarator optimization is traditionally performed as a two-stage process; in the first stage, an optimal equilibrium is calculated which balances a set of competing constraints, and in the second stage a set of coils is found that supports said equilibrium. Stage one is generally performed using a rep…
▽ More
We present the first single-stage optimization of islands in finite-$β$ stellarator equilibria. Stellarator optimization is traditionally performed as a two-stage process; in the first stage, an optimal equilibrium is calculated which balances a set of competing constraints, and in the second stage a set of coils is found that supports said equilibrium. Stage one is generally performed using a representation for the equilibrium that assumes nestedness of flux surfaces, even though this is not warranted and occasionally undesired. The second stage optimization of coils is never perfect, and the mismatch leads to worse performing equilibria, and further deteriorates if additional constraints such as force minimization, coil torsion or port access are included. The higher fidelity of single-stage optimization is especially important for the optimization of islands as these are incredibly sensitive to changes in the field. In this paper we demonstrate an optimization scheme capable of optimizing islands in finite $β$ stellarator equilibria directly from coils. We furthermore develop and demonstrate a method to reduce the dimensionality of the single-stage optimization problem to that of the first stage in the two-stage approach.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
Design and First Tests of the Trapped Electrons Experiment T-REX
Authors:
Francesco Romano,
Guillaume Le Bars,
Joaquim Loizu,
Marc Nöel,
Jean-Philippe Hogge,
Stefano Alberti,
Jeremy Genoud,
Severino Antonioni,
Lorraine Naux,
Pierrick Giroud-Garampon,
Steve Couturier,
Thierry Leresche,
Damien Fasel
Abstract:
Gyrotrons are the only sources used for electron cyclotron resonance heating (ECRH). Past gyrotron experiments have highlighted some instability issues, leading to restricted operating ranges. One major cause is the possible presence of trapped electrons within the gyrotron's magnetron injection gun (MIG) region, resulting in undesired currents and subsequent operational failures. Hereby, we prese…
▽ More
Gyrotrons are the only sources used for electron cyclotron resonance heating (ECRH). Past gyrotron experiments have highlighted some instability issues, leading to restricted operating ranges. One major cause is the possible presence of trapped electrons within the gyrotron's magnetron injection gun (MIG) region, resulting in undesired currents and subsequent operational failures. Hereby, we present the design and initial experimental findings of the TRapped Electrons eXperiment (T-REX), a novel and unique basic plasma experiment constructed at the Swiss Plasma Center to investigate the underlying physics of trapped electron clouds, therefore entering the realm of non-neutral plasmas. The primary objective is to explore the formation and evolution of electron clouds in gyrotron MIG designs. Typical MIG conditions are replicated in T-REX in terms of geometries, electric and magnetic fields, and supported by kinetic simulations via the 2D Particle-in-Cell (PIC) code FENNECS. T-REX is characterized by two coaxial electrodes, placed in a vacuum chamber on top of a superconducting magnet, with an applied external radial electric field (0.1 - 2 MV/m), and a nearly axial magnetic field (B<0.4T) inducing azimuthal drifts and confining electrons with energies between 0.1-1 keV. The experimental diagnostics include current, voltage, and Faraday probes, and in future plans imaging via phosphor screen, optical emission spectroscopy, streak camera imaging, and, potentially, electric field distribution via the Stark effect. Moreover, T-REX is also designed to be capable of detecting potential growth of diocotron modes via such diagnostics in combination with a segmented outer electrode. Through T-REX, we aim to gain valuable insights into the trapping and dynamics of electrons within MIG regions, paving the way for improved gyrotron performance and reliability in magnetic fusion reactors.
△ Less
Submitted 27 June, 2024;
originally announced June 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.
-
Constructing nested coordinates inside strongly shaped toroids using an action principle
Authors:
Zeno Tecchiolli,
Stuart Hudson,
Joaquim Loizu,
Robert Köberl,
Florian Hindenlang,
Brenno De Lucca
Abstract:
A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0, 2π] x [0, 2π] x [0, 1] to a toroidal domain in R3. This approach employs an action built on the squa…
▽ More
A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0, 2π] x [0, 2π] x [0, 1] to a toroidal domain in R3. This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier-Zernike basis via action minimization. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of 3D MHD equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Parallel flows as a key component to interpret Super-X divertor experiments
Authors:
M. Carpita,
O. Février,
H. Reimerdes,
C. Theiler,
B. P. Duval,
C. Colandrea,
G. Durr-Legoupil-Nicoud,
D. Galassi,
S. Gorno,
E. Huett,
J. Loizu,
L. Martinelli,
A. Perek,
L. Simons,
G. Sun,
E. Tonello,
C. Wüthrich,
the TCV team
Abstract:
The Super-X Divertor (SXD) is an alternative divertor configuration leveraging total flux expansion at the Outer Strike Point (OSP). While the extended 2-Point Model (2PM) predicts facilitated detachment access and control in the SXD configuration, these attractive features are not always retrieved experimentally. These discrepancies are at least partially explained by the effect of parallel flows…
▽ More
The Super-X Divertor (SXD) is an alternative divertor configuration leveraging total flux expansion at the Outer Strike Point (OSP). While the extended 2-Point Model (2PM) predicts facilitated detachment access and control in the SXD configuration, these attractive features are not always retrieved experimentally. These discrepancies are at least partially explained by the effect of parallel flows which, when self-consistently included in the 2PM, reveal the role of total flux expansion on the pressure balance and weaken the total flux expansion effect on detachment access and control, compared to the original predictions. This new model can partially explain the discrepancies between the 2PM and experiments performed on tokamak à configuration variable (TCV), in ohmic L-mode scenarios, which are particularly apparent when scanning the OSP major radius Rt. In core density ramps in lower Single-Null (SN) configuration, the impact of Rt on the CIII emission front movement in the divertor outer leg - used as a proxy for the plasma temperature in the divertor - is substantially weaker than 2PM predictions. Furthermore, in OSP radial sweeps in lower and upper SN configurations, in ohmic L-mode scenarios with a constant core density, the peak parallel particle flux density at the OSP is almost independent of Rt, while the 2PM predicts a linear dependence. Finally, analytical and numerical modeling of parallel flows in the divertor is presented. It is shown that an increase in total flux expansion can favour supersonic flows at the OSP. Parallel flows are also shown to be relevant by analysing SOLPS-ITER simulations of TCV.
△ Less
Submitted 27 February, 2024; v1 submitted 30 June, 2023;
originally announced June 2023.
-
Nonlinear saturation of resistive tearing modes in a cylindrical tokamak with and without solving the dynamics
Authors:
J. Loizu,
D. Bonfiglio
Abstract:
We show that the saturation of resistive tearing modes in a cylindrical tokamak, as well as the corresponding island width, can be directly calculated with an MHD equilibrium code without solving the dynamics and without considering resistivity. The results are compared to initial value resistive MHD simulations and to an analytical nonlinear theory. For small enough islands, the agreement is rema…
▽ More
We show that the saturation of resistive tearing modes in a cylindrical tokamak, as well as the corresponding island width, can be directly calculated with an MHD equilibrium code without solving the dynamics and without considering resistivity. The results are compared to initial value resistive MHD simulations and to an analytical nonlinear theory. For small enough islands, the agreement is remarkable. For sufficiently large islands, the equilibrium calculations, which assume a flat current profile inside the island, overestimate the saturation amplitude. On the other hand, excellent agreement between nonlinear resistive MHD simulations and nonlinear theory is observed for all the considered tearing unstable equilibria.
△ Less
Submitted 21 August, 2023; v1 submitted 26 May, 2023;
originally announced May 2023.
-
Validation of GBS plasma turbulence simulation of the TJ-K stellarator
Authors:
A. J. Coelho,
J. Loizu,
P. Ricci,
M. Ramisch,
A. Köhn-Seemann,
G. Birkenmeier,
K. Rahbarnia
Abstract:
We present a validation of a three-dimensional, two-fluid simulation of plasma turbulence in the TJ-K stellarator, a low temperature plasma experiment ideally suited for turbulence measurements. The simulation is carried out by the GBS code, recently adapted to simulate 3D magnetic fields. The comparison shows that GBS retrieves the main turbulence properties observed in the device, namely the fac…
▽ More
We present a validation of a three-dimensional, two-fluid simulation of plasma turbulence in the TJ-K stellarator, a low temperature plasma experiment ideally suited for turbulence measurements. The simulation is carried out by the GBS code, recently adapted to simulate 3D magnetic fields. The comparison shows that GBS retrieves the main turbulence properties observed in the device, namely the fact that transport is dominated by fluctuations with low poloidal mode number. The poloidal dependence of the radial $\text{E}\times\text{B}$ turbulent flux is compared on a poloidal plane with elliptical flux surfaces, where a very good agreement between experiment and simulation is observed, and on another with triangular flux surfaces, which shows a poorer comparison. The fluctuation levels in both cases are underestimated in the simulations. The equilibrium density profile is well retrieved by the simulation, while the electron temperature and the electrostatic potential profiles, which are very sensitive to the strength and localization of the sources, do not agree well with the experimental measurements.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium
Authors:
Yi-Min Huang,
Yao Zhou,
Joaquim Loizu,
Stuart Hudson,
Amitava Bhattacharjee
Abstract:
Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces and a continuous rotational transform. These currents have two components: a surface current (Dirac $δ$-function in flux surface labeling) that prevents the formation of magnetic islands and an algebraically divergent Pfirsch--Schlüter current de…
▽ More
Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces and a continuous rotational transform. These currents have two components: a surface current (Dirac $δ$-function in flux surface labeling) that prevents the formation of magnetic islands and an algebraically divergent Pfirsch--Schlüter current density when a pressure gradient is present across the rational surface. At flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch--Schlüter current density scaling as $J\sim1/Δι$, where $Δι$ is the difference of the rotational transform relative to the rational surface. If the distance $s$ between flux surfaces is proportional to $Δι$, the scaling relation $J\sim1/Δι\sim1/s$ will lead to a paradox that the Pfirsch--Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm\textendash Kulsrud\textendash Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch--Schlüter current density but also the diamagnetic current density are divergent as $\sim1/Δι$. However, due to the formation of a Dirac $δ$-function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with $s\sim(Δι)^{2}$. Consequently, the singular current density $J\sim1/\sqrt{s}$, making the total current finite, thus resolving the paradox.
△ Less
Submitted 9 February, 2024; v1 submitted 3 March, 2023;
originally announced March 2023.
-
Equilibrium $β$-limits dependence on bootstrap current in classical stellarators
Authors:
A. Baillod,
J. Loizu,
Z. Qu,
H. P. Arbez,
J. P. Graves
Abstract:
While it is important to design stellarators with high magneto-hydrodynamic (MHD) stability $β$-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of $β$ values. As $β$ increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another ki…
▽ More
While it is important to design stellarators with high magneto-hydrodynamic (MHD) stability $β$-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of $β$ values. As $β$ increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another kind of $β$-limit, the so-called equilibrium $β$-limit. In this paper, we explore numerically the dependence of the equilibrium $β$-limit on the bootstrap current strength using the Stepped Pressure Equilibrium Code (SPEC). We develop a diagnostic to determine whether or not magnetic islands are expected to participate significantly to radial transport, and we build an analytical model to predict the expected equilibrium $β$-limit, which recovers the main features of the numerical results. This research opens the possibility to include additional targets in stellarator optimization functions, provides additional understanding on the existence of magnetic surfaces at large $β$, and is a step forward in the understanding of the equilibrium $β$-limit.
△ Less
Submitted 16 May, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
On the relationship between the multi-region relaxed variational principle and resistive inner layer theory
Authors:
A. Kumar,
J. Loizu,
M. J. Hole,
Z. Qu,
S. R. Hudson,
R. L Dewar
Abstract:
We show that the variational energy principle of multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by "ideal interfaces", and in each volume the magnetic field relaxes to a Taylor state where the pressure gradient $\nabla p = 0$. The MRxMHD model is imple…
▽ More
We show that the variational energy principle of multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by "ideal interfaces", and in each volume the magnetic field relaxes to a Taylor state where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the SPEC code so that the equilibrium solution in each region is computed while the preserving force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low beta plasma, we have been able to illustrate a relationship between the resistive singular layer theory [Coppi et al. (1966) Nucl. Fusion 6 101, Glasser et al. The Physics of Fluids 18, 875-888 (1975)], and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner layer theories.
△ Less
Submitted 11 September, 2022;
originally announced September 2022.
-
Self-consistent formation and steady-state characterisation of trapped high energy electron clouds in the presence of a neutral gas background
Authors:
G. Le Bars,
J. -Ph. Hogge,
J. Loizu,
S. Alberti,
F. Romano,
A. Cerfon
Abstract:
This study considers the self-consistent formation and dynamics of electron clouds interacting with a background neutral gas through elastic and inelastic (ionisation) collisions in coaxial geometries similar to gyrotron electron guns. These clouds remain axially trapped as the result of crossed magnetic field lines and electric equipotential lines creating potential wells similar to those used in…
▽ More
This study considers the self-consistent formation and dynamics of electron clouds interacting with a background neutral gas through elastic and inelastic (ionisation) collisions in coaxial geometries similar to gyrotron electron guns. These clouds remain axially trapped as the result of crossed magnetic field lines and electric equipotential lines creating potential wells similar to those used in Penning traps. Contrary to standard Penning traps, in this study we consider a strong externally applied radial electric field which is of the same order as that of the space-charge field. In particular, the combination of coaxial geometry, strong radial electric fields and electron collisions with the residual neutral gas (RNG) present in the chamber induce non-negligible radial particle transport and ionisation. In this paper, the dynamics of the cloud density and currents resulting from electron-neutral collisions are studied using a 2D3V particle-in-cell code. Simulation results and parametric scans are hereby presented. Finally, a fluid model is derived to explain and predict the cloud peak density and peak radial current depending on the externally applied electric and magnetic fields, and on the RNG pressure.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Global fluid simulation of plasma turbulence in a stellarator with an island divertor
Authors:
A. J. Coelho,
J. Loizu,
P. Ricci,
M. Giacomin
Abstract:
Results of a three-dimensional, flux-driven, electrostatic, global, two-fluid turbulence simulation for a 5-field period stellarator with an island divertor are presented. The numerical simulation is carried out with the GBS code, recently extended to simulate plasma turbulence in non-axisymmetric magnetic equilibria. The vacuum magnetic field used in the simulation is generated with the theory of…
▽ More
Results of a three-dimensional, flux-driven, electrostatic, global, two-fluid turbulence simulation for a 5-field period stellarator with an island divertor are presented. The numerical simulation is carried out with the GBS code, recently extended to simulate plasma turbulence in non-axisymmetric magnetic equilibria. The vacuum magnetic field used in the simulation is generated with the theory of Dommaschk potentials, and describes a configuration with a central region of nested flux surfaces, surrounded by a chain of magnetic islands, similarly to the diverted configurations of W7-X. The heat outflowing from the core reaches the island region and is transported along the magnetic islands, striking the vessel walls, which correspond to the boundary of the simulation domain.
The radial transport of particles and heat is found to be mainly driven by a field-aligned coherent mode with poloidal number $m=4$. The analysis of this mode, based on non-local linear theory considerations, shows its ballooning nature. In contrast to tokamak simulations and experiments, where blobs often contribute to transport, we do not observe the presence of intermittent transport events.
△ Less
Submitted 26 January, 2022;
originally announced January 2022.
-
Stellarator optimization for nested magnetic surfaces at finite $β$ and toroidal current
Authors:
A. Baillod,
J. Loizu,
J. P. Graves,
M. Landreman
Abstract:
Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this pap…
▽ More
Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this paper, we build upon their work to show the first finite-$β$, fixed- and free-boundary optimization of SPEC equilibria for good magnetic surfaces. The objective function is constructed with the Greene's residue of selected rational surfaces and the optimization is driven by the SIMSOPT framework [M. Landreman et al., J. Open Source Software 6, 3525 (2021)]. We show that the size of magnetic islands and the consequent regions occupied by chaotic field lines can be minimized in a classical stellarator geometry by optimizing either the injected toroidal current profile, the shape of a perfectly conducting wall surrounding the plasma (fixed-boundary case), or the coils (free-boundary case), in a reasonable amount of computational time. This work shows that SPEC can be used as an equilibrium code both in a two-step or single-step stellarator optimization loop.
△ Less
Submitted 30 November, 2021;
originally announced November 2021.
-
Numerical study of $δ$-function current sheets arising from resonant magnetic perturbations
Authors:
Yi-Min Huang,
Stuart R. Hudson,
Joaquim Loizu,
Yao Zhou,
Amitava Bhattacharjee
Abstract:
General three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially reg…
▽ More
General three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal MHD limit has been a significant challenge. This work presents numerical solutions of the Hahm-Kulsrud-Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad-Shafranov (GS) solver and the SPEC code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), where the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution approaches the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate excellent agreement between the numerical solutions obtained from the two codes through a thorough convergence study.
△ Less
Submitted 1 September, 2023; v1 submitted 20 August, 2021;
originally announced August 2021.
-
On the non-existence of stepped-pressure equilibria far from symmetry
Authors:
Z. S. Qu,
S. R. Hudson,
R. L. Dewar,
J. Loizu,
M. J. Hole
Abstract:
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In th…
▽ More
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. We found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov-Arnold-Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.
△ Less
Submitted 14 September, 2021; v1 submitted 4 August, 2021;
originally announced August 2021.
-
Properties of a new quasi-axisymmetric configuration
Authors:
S. A. Henneberg,
M. Drevlak,
C. Nührenberg,
C. D. Beidler,
Y. Turkin,
J. Loizu,
P. Helander
Abstract:
A novel, compact, quasi-axisymmetric configuration is presented which exhibits low fast-particle losses and is stable to ideal MHD instabilities. The design has fast-particle loss rates below 8\% for flux surfaces within the half-radius, and is shown to have an MHD-stability limit of a normalised pressure of $\langleβ\rangle=3\%$ where $\langleβ\rangle$ is volume averaged. The flux surfaces at var…
▽ More
A novel, compact, quasi-axisymmetric configuration is presented which exhibits low fast-particle losses and is stable to ideal MHD instabilities. The design has fast-particle loss rates below 8\% for flux surfaces within the half-radius, and is shown to have an MHD-stability limit of a normalised pressure of $\langleβ\rangle=3\%$ where $\langleβ\rangle$ is volume averaged. The flux surfaces at various plasma betas and currents as calculated using the SPEC equilibrium code are presented. Neoclassical transport coefficients are shown to be similar to an equivalent tokamak, with a distinct banana regime at half-radius. An initial coil design study is presented to assess the feasibility of this configuration as a fusion-relevant experiment.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.