Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stellarator optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Direct optimization of neoclassical ion transport in stellarator reactors

Abstract We directly optimize stellarator neoclassical ion transport while holding neoclassical electron transport at a moderate level, creating a scenario favorable for impurity expulsion and retaining good ion confinement. Traditional neoclassical stellarator optimization has focused on minimizing ϵ eff , the geometric factor that characterizes the amount of radial transport due to particles in the 1 / ν regime. Under expected reactor-relevant conditions, core electrons will be in the 1 / ν regime and core fuel ions will be in the ν regime. Traditional optimizations thus minimize electron transport and rely on the radial electric field ( E r ) that develops to confine the ions. This often results in an inward-pointing E r that drives high- Z impurities into the core, which may be troublesome in future reactors. In this work, we increase the ratio of the thermal transport coefficients L 11 e / L 11 i , which previous research has shown can create an outward-pointing E r . This effect is very beneficial for impurity expulsion. We obtain self-consistent density, temperature, and E r profiles at reactor-relevant conditions for an optimized equilibrium. This equilibrium is expected to enjoy significantly improved impurity transport properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The DESC stellarator code suite Part 3: Quasi-symmetry optimization

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas

Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Optimizing the HSX stellarator for microinstability by coil-current adjustments

Abstract The optimization of helically symmetric experiment (HSX) for reduced microinstability has been achieved by examining a large set of configurations within a neighborhood of the standard operating configuration. This entailed generating a database of more than 10 6 magnetic-field configurations for HSX by varying the currents in external coils. Using a set of volume-averaged metrics and gyrokinetic simulations, this database has helped to identify a set of configurations that can be used to regulate trapped-electron-mode stability in HSX. This set of configurations is also found to correlate flux-surface elongation and triangularity with an increase in magnetic-well depth, an increase in rotational transform, and low neoclassical heat-flux relative to the standard quasi-helically-symmetric configuration. These results demonstrate sensitivity of plasma behavior in response to changes in a 3D magnetic field to both neoclassical and gyrokinetic models, and the experimental potential in HSX to explore turbulence optimization. This perturbative optimization approach is not unique to HSX, and can readily be deployed on existing fusion devices to identify novel magnetic-fields to be used in turbulence-optimization experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Robustness of quasi-symmetry along parametric boundary variation

Quasi-symmetry is a guiding principle to modern stellarator optimization for improved plasma confinement. However, the robustness of optimized configurations, which can be crucial for maintaining performance under diverse engineering constraints and practical limitations, has received relatively little attention. Here, in this work, we present various case studies on this robustness, by investigating variations in 1/ v neoclassical transport when plasma configurations are smoothly altered across distinct optimized targets. These targets, optimized from different families—quasi-axisymmetric, quasi-helical, and quasi-isodynamic—are approximately matched in major radius as part of a flexible stellarator design. Our study shows that an optimized target does not always represent a local minimum in transport and that the robustness of a local minimum when present can vary significantly. Furthermore, there are configurations which belong to no established families but have transport levels as low as those of optimized targets. These results highlight the importance of conducting extended searches with key parametric variations around optimized configurations, to ensure its robustness as well as flexibility if desired.

neoclassical transport↗
DOE OSTI2025

The shear Alfvén continuum of quasisymmetric stellarators

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

fusion plasma↗
DOE OSTI2022

Conceptual design of a heavy ion beam probe diagnostic for the Wendelstein 7-X Stellarator

This article describes the current state of the design of the heavy ion beam probe (HIBP) for Wendelstein 7-X (W7-X). It will be the first HIBP diagnostic on an optimized stellarator and is designed to study electric fields and ion scale turbulence in all W7-X reference magnetic configurations. The use of an existing 2 MV accelerator, located outside of the torus hall, results in the need for a circuitous primary beamline. This increases the complexity of the ion optics design to deliver a focused beam to the plasma. To access most of the magnetic configuration space of W7-X, the secondary beamline and an energy analyzer are designed to pivot, thereby redirecting a wider range of secondary beam trajectories. Signal level estimates indicate that the equilibrium potential can be measured at all radii and that the radial coverage for potential and density fluctuations measurements depends on the plasma density.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Forecasting Dark Matter Subhalo Constraints from Stellar Streams using Implicit Likelihood Inference

The evidence for dark matter (DM) remains compelling, although attempts to understand its particle nature remain inconclusive. One promising method to study DM is detecting DM subhalos through their gravitational interactions with stellar streams. In this study, we apply Neural Posterior Estimation (NPE) to constrain subhalo interaction parameters, including mass, scale radius, velocity, and encounter geometry, from stellar stream kinematics. We generate particle spray simulations based on the Lagrange Cloud stripping technique, focusing on the ATLAS-Aliqa Uma stream as a test case. We train multiple NPE models across multiple observational scenarios, quantifying how kinematic completeness affects inference and forecasting constraints from upcoming surveys including LSST, 4MOST, and 10-year Gaia data. Our results demonstrate that NPE can produce accurate and well-calibrated posteriors. In the idealized case with full 6D coordinates, we achieve subhalo mass uncertainties of 15-20% for a $10^7 \, \mathrm{M_\odot}$ subhalo, with 5D coordinates (excluding radial velocities) achieving similar performance. Under realistic observational conditions, mass uncertainties range from 50% (present-day) to 20-40% (future scenarios), with comparable performance between the photometric-only LSST sample and a smaller sample that includes Gaia proper motions and 4MOST radial velocities. Most notably, we find that velocity bimodality emerges when phase space is poorly sampled, whether due to missing kinematic information or limited stellar tracers. Combining large photometric samples with targeted spectroscopic follow-up can effectively resolves this degeneracy. These results demonstrate the power of implicit likelihood inference for optimizing stellar stream observational strategies and forecasting DM subhalo constraints from upcoming surveys.

Nguyen, Tri [Northwestern U. (main); SkAI, Chicago↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Error field measurements and correction on MUSE permanent magnet stellarator

The magnetic field topology and error fields are measured in the MUSE permanent magnet stellarator with an electron beam and fluorescent rod. The effective magnetization of the magnets is determined by measuring rational flux surface locations. Electron magnetic drift effects must be considered self-consistently due to the small magnetic shear on MUSE. Error fields are characterized by flux surface shape deviation and magnetic island chain width and phase. Error fields are then corrected by adjusting the permanent magnet (PM) holders as rigid bodies. Resonant error fields are on the order of or smaller than 3 x 10 -6 of the toroidal field after correction and surface shape matches well with the design. Resonant components of the quasi-symmetry deviation are estimated, using a novel magnetic drift measurement, to be in the 10 −4 range. The PM approach is validated as a low-cost and precise way to build optimized stellarators.

electron beam mapping↗
DOE OSTI2024

Magnetic fields with general omnigenity

Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has been either axisymmetric, quasi-symmetric or has poloidally closed contours of magnetic field strength $B$ . However, general omnigenous equilibria are a much larger design space than these subsets. A new model is presented and employed in the DESC stellarator optimization suite to represent and discover the full parameter space of omnigenous equilibria. Although exact omnigenity aside from quasi-symmetry is impossible, these results reveal that excellent particle confinement can be achieved in practice. Examples far from quasi-symmetry with poloidally, helically and toroidally closed $B$ contours are attained with DESC and shown to have low neoclassical collisional transport and fast particle losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2023

Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wendelstein 7-X

We assess the turbulent particle transport being responsible for the limitation of the confinement and, thus, the overall performance of the neoclassically optimized stellarator Wendelstein 7-X. The radial particle transport is experimentally inferred from the evaluation of impurity injection into turbulence reduced and enhanced plasma scenarios revealing a completely different confinement behavior. The impact of the density gradient on the turbulent ion transport is theoretically estimated using large-scale non-linear gyro-kinetic simulations enabling, for the first time in Wendelstein 7-X, a quantitative comparison to the experimentally assessed impurity transport properties. We demonstrate that impurity transport in most of the Wendelstein 7-X discharges, up to now impossible to cover only with neoclassical estimations, is dominated by turbulence and can be modelled via gyro-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimization of nonlinear turbulence in stellarators

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 quasi-symmetry 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Optimization of passive superconductors for shaping stellarator magnetic fields

Here, we consider the problem of optimizing a set of passive superconducting coils (PSCs) with currents induced by a background magnetic field rather than power supplies. In the nuclear fusion literature, such coils have been proposed to partially produce the 3D magnetic fields for stellarators and provide passive stabilization. We perform the first optimizations of PSC arrays with respect to the orientation, shape, and location of each coil, jointly minimized with the background fields. We conclude by generating passive coil array solutions for four stellarators.

coil optimization↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On the effect of flux-surface shaping on trapped-electron modes in quasi-helically symmetric stellarators

Using a novel optimization procedure, it has been shown that the Helically Symmetric eXperiment stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [Gerard et al., Nucl. Fusion 63, (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available energy EA, quasi-helical symmetry, and flux-surface shaping parameters as metrics for TEM stabilization. It is found that lower values of EA 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 EA. 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.

Physics↗