Engineering PapersSearch

DOE OSTI · 3017855

Omnigenous umbilic stellarators

Abstract

To better understand the dependence of the magnetic field structure in the plasma edge on the plasma boundary shape, in the context of X-point and island divertor designs, we define and develop a class of stellarators called umbilic stellarators. These equilibria are characterised by a single continuous high-curvature edge on the plasma boundary that goes around multiple times toroidally before meeting itself. We develop a technique that allows us to simultaneously optimise the plasma boundary along with a curve lying on the boundary on which we impose a high curvature while imposing omnigenity – a property of the magnetic field that ensures trapped particle confinement throughout the plasma volume. We find that umbilic stellarators naturally tend to favour piecewise omnigenity instead of omnigenity with a specific helicity. After generating omnigenous umbilic stellarators, we design coil sets for some of them and explore the fieldline structure in the edge and its sensitivity to small fluctuations in the plasma. Finally, using single-stage optimisation, we simultaneously modify the plasma and coil shape and propose an experiment to modify an existing tokamak to a finite-𝛽 stellarator using this technique and explore a potentially simpler way to convert a limited tokamak into a diverted stellarator.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Gaur, Rahul [Princeton Univ., NJ (United States)] (ORCID:0000000343679052), Panici, D. G. [Princeton Univ., NJ (United States)] (ORCID:0000000307364360), Elder, T. M. [Max Planck Institute for Plasma Physics, Greifswald (Germany); Univ. of Maryland, College Park, MD (United States)], Landreman, M. [Univ. of Maryland, College Park, MD (United States)] (ORCID:000000027233577X), Unalmis, K. [Princeton Univ., NJ (United States)] (ORCID:0000000217556764), Elmacioglu, Y. G. [Princeton Univ., NJ (United States)] (ORCID:0009000277472834), Dudt, D. W. [Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)] (ORCID:0000000245573529), Conlin, R. [Univ. of Maryland, College Park, MD (United States)] (ORCID:0000000183662111), Kolemen, E. [Princeton Univ., NJ (United States); Thea Energy, Princeton, NJ (United States)] (ORCID:0000000342123247). 2026-01-26. Omnigenous umbilic stellarators. https://doi.org/10.1017/s0022377825100925

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma

Spectrally accurate, reverse-mode differentiable bounce-averaging algorithm and its applications

We present a fast, spectrally (exponentially) accurate, automatically differentiable bounce-averaging algorithm that is used to simplify kinetic models. Using this algorithm, implemented in the DESC stellarator optimisation suite, we can perform efficient optimisation of many objectives to improve stellarator performance, such as the effective ripple 𝜖 eff metric for the neoclassical transport coefficient in the low collisionality regime and proxies for energetic particle confinement. For the first time, we optimise a finite-beta stellarator to directly reduce neoclassical ripple transport using reverse-mode differentiation. This ensures the computational cost of differentiation is independent of the number of controllable parameters.

fusion plasma

Constructing field-aligned coordinate systems for gyrokinetic simulations of tokamaks in X-point geometries

Structures in tokamak plasmas are elongated along the direction of the magnetic field and short in the directions perpendicular to the magnetic field. Many tokamak simulation codes take advantage of this by using a field-aligned coordinate system. However, field-aligned coordinate systems have a coordinate singularity at magnetic X-points where the poloidal magnetic field vanishes, which makes it difficult to use field-aligned coordinate systems when simulating the core and scrape-off layer simultaneously. Here, we present an algorithm for grid generation and computing geometric quantities in a standard field-aligned coordinate system that avoids the singularity and allows one to conduct two-dimensional gyrokinetic axisymmetric simulations in X-point geometries. Convergence tests of advection, boundary value problems and geometric quantities all show greater than first-order convergence even in the vicinity of the X-point. We also demonstrate the geometric consistency of our algorithm with an example simulation of the spherical tokamak for energy production, which shows machine-precision particle conservation.

fusion plasma