SEARCH · Search NASA
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.
We developed, DESC, a new pseudo-spectral stellarator equilibrium and optimization code. The algorithm was specifically designed to be fast (runs on GPU) and easy use. The code is now being used many stellarator fusion companies and institutions.
A new coherence imaging spectroscopy (CIS) diagnostic optimized to measure the C 2+ impurity ion temperature T i spatial distribution in the divertor plasma of the W7-X stellarator is designed, tested, and validated. Using CIS to obtain T i in the edge of magnetically confined plasmas has historically been challenging because Doppler broadening and Zeeman splitting have comparable effects on the shape of spectral emission lines. To distinguish between these two mechanisms, a novel approach to birefringent crystal design is employed to minimize the diagnostic’s sensitivity to Zeeman splitting. The recently developed pixelated multi-delay CIS approach is also used to obtain four times as much spectral information as traditional CIS approaches. The T i -optimized CIS diagnostic is validated in a long-pulse W7-X plasma by comparison with a high-resolution spectrometer whose sightlines overlap with the CIS field of view. The CIS and spectrometer T i profiles have the same shape and agree to within 10% on average and 25% in the worst case. Images of the T i distribution near the divertor show toroidally elongated bands aligned with the magnetic field, with T i ranging between 10 and 40 eV.
Predictive 3D optimization reveals a novel approach to modify a nonaxisymmetric magnetic perturbation to be entirely harmless for tokamaks, by essentially restoring quasisymmetry in perturbed particle orbits as much as possible. Such a quasisymmetric magnetic perturbation (QSMP) has been designed and successfully tested in the KSTAR and DIII-D tokamaks, demonstrating no performance degradation despite the large overall amplitudes of nonaxisymmetric fields and strong response otherwise expected in the tested plasmas. The results indicate that a quasisymmetric optimization is a robust path of error field correction across the resonant and nonresonant field spectrum in a tokamak, leveraging the prevailing concept of quasisymmetry for general 3D plasma confinement systems such as stellarators. Finally, the optimization becomes, in fact, a simple eigenvalue problem to the so-called torque response matrices if a perturbed equilibrium is calculated consistent with nonaxisymmetric neoclassical transport.
Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.
Deuterium–tritium (DT) burning requires a long energy confinement times compared to collision times, so the particle distribution functions must approximate local Maxwellians. Non-equilibrium thermodynamics is applicable, which gives relations among transport, entropy production, the collision frequency, and the deviation from a Maxwellian. The distribution functions are given by the Fokker–Planck equation, which is an advection–diffusion equation. A large hyperbolic operator, the Vlasov operator with the particle trajectories as its characteristics, equals a small diffusive operator, the collision operator. The collisionless particle trajectories would be chaotic in stellarators without careful optimization. This would lead to rapid entropy production and transport—far beyond what is consistent with a self-sustaining DT burn. Omnigeneity is the weakest general condition that is consistent with a sufficiently small entropy production associated with the thermal particle trajectories. Omnigeneity requires that the contours of constant magnetic field strength be unbounded in at least one of the two angular coordinates in magnetic surfaces and that there be a symmetry in the field-strength wells along the field lines. Even in omnigenous plasmas, fluctuations due to microturbulence can produce chaotic particle trajectories and the gyro-Bohm transport is seen in many stellarator and tokamak experiments. The higher the plasma temperature above 10 keV, the smaller the transport must be compared to gyro-Bohm for a self-sustaining DT burn. The hot alphas of DT fusion heat the electrons. When the ion–electron equilibration time is long compared to the ion energy confinement time, a self-sustaining DT burn is not possible, which sets a limit on the electron temperature.
A common optimization problem in the areas of magnetized plasmas and fusion energy is the design of magnets to produce a given three-dimensional magnetic field distribution to high precision. When designing arrays of permanent magnets for stellarator plasma confinement, such problems have tens of thousands of degrees of freedom whose solutions, for practical reasons, should be constrained to discrete spaces. We perform a direct comparison between two algorithms that have been developed previously for this purpose, and demonstrate that composite procedures that apply both algorithms in sequence can produce substantially improved results. One approach uses a continuous, quasi-Newton procedure to optimize the dipole moments of a set of magnets and then projects the solution onto a discrete space. The second uses an inherently discrete greedy optimization procedure that has been enhanced and generalized for this work. Further, the approaches are both applied to design arrays cubic rare-Earth permanent magnets to confine a quasi-axisymmetric plasma with a magnetic field on axis of 0.5 T. The first approach tends to find solutions with higher field accuracy, whereas the second can find solutions with substantially (up to 30%) fewer magnets. When the approaches are combined, they can obtain solutions with magnet quantities comparable to the second approach while matching the field accuracy of the first.
The Southern Photometric Local Universe Survey (S-PLUS) is a novel project that aims to map the Southern Hemisphere using a twelve filter system, comprising five broad-band SDSS-like filters and seven narrow-band filters optimized for important stellar features in the local universe. In this paper we use the photometry and morphological information from the first S-PLUS data release (S-PLUS DR1) cross-matched to unWISE data and spectroscopic redshifts from Sloan Digital Sky Survey DR15. We explore three different machine learning methods (Gaussian Processes with GPz and two Deep Learning models made with TensorFlow) and compare them with the currently used template-fitting method in the S-PLUS DR1 to address whether machine learning methods can take advantage of the twelve filter system for photometric redshift prediction. Using tests for accuracy for both single-point estimates such as the calculation of the scatter, bias, and outlier fraction, and probability distribution functions (PDFs) such as the Probability Integral Transform (PIT), the Continuous Ranked Probability Score (CRPS) and the Odds distribution, we conclude that a deep-learning method using a combination of a Bayesian Neural Network and a Mixture Density Network offers the most accurate photometric redshifts for the current test sample. In conclusion, it achieves single-point photometric redshifts with scatter (σ NMAD ) of 0.023, normalized bias of -0.001, and outlier fraction of 0.64% for galaxies with r_auto magnitudes between 16 and 21.
A Grad-Shafranov equation (GSE) valid for compact quasisymmetric stellarators is derived by an asymptotic expansion around a vacuum field carried to first order. We obtain an equation for the existence of flux surfaces leading up to the GSE. The flux surface label must simultaneously satisfy the existence equation and the GSE, which generally leads to an overdetermined problem. We show how the overdetermined problem can be resolved within our model for a class of hybrid devices similar to that studied by Henneberg and Plunk [Phys. Rev. Res. 6, L022052 (2024)]. We are also able to solve the existence equation for flux surfaces analytically in the most general case by introducing a special coordinate system. This will enable us to carry out an optimization seeking to minimize the error in our GSE while obeying the flux surface existence equation, which will allow us to find solutions outside the class of hybrid devices. This will allow for a coarse-grained approximate search in the space of quasisymmetric equilibria that should be faster than a conventional stellarator optimization. Nevertheless, it would still be necessary to fine-tune the approximate solutions using conventional tools to obtain a more precise optimized equilibrium.
A heavy ion beam probe (HIBP) diagnostic is being developed for studies of plasma equilibrium and turbulence in the optimized Wendelstein 7-X (W7-X) stellarator. Operation of W7-X has experimentally demonstrated that its optimized magnetic field results in improved neoclassical particle confinement and, as a result, turbulence is the predominant cause of energy transport. The HIBP will have the unique ability to provide experimental data needed to complement models of both neoclassical and turbulent transport. It will acquire direct measurements in the W7-X plasma interior of the electric potential (needed for understanding ambipolar particle flux) and fluctuations of electron density and potential (needed for understanding turbulence). The HIBP for W7-X will inject singly charged ion beams with energies of up to 2 MeV and is designed to access the upper cross section of the W7-X plasma. We use trajectory simulations to illustrate the plasma coverage that the diagnostic can achieve in the reference magnetic configurations of W7-X. We calculate beam signal levels, discuss anticipated measurement sensitivity of broadband fluctuations of electron density and plasma potential, and show how they depend on plasma density. We also discuss the diagnostic sensitivity to equilibrium plasma potential.
The transport of carbon impurities has been studied in the helically symmetric stellarator experiment (HSX) using active and passive charge exchange recombination spectroscopy (CHERS). For the analysis of the CHERS signals, the STRAHL impurity transport code has been re-written in the python programming language and optimized for the application in stellarators. In addition, neutral hydrogen densities both along the NBI line of sight as well as for the background plasma have been calculated using the FIDASIM code. By using the basinhopping algorithm to minimize the difference between experimental and predicted active and passive signals, significant levels of impurity diffusion are observed. In this work, comparisons with neoclassical calculations from DKES/PENTA show that the inferred levels exceed the neoclassical transport by about a factor of four in the core and more than 100 times towards the plasma edge, thus indicating a high level of anomalous transport. This observation is in agreement with experimental heat diffusivites determined from a power balance analysis which exhibits strong anomalous transport as well.
An adjoint method to calculate the gradient of island width in stellarators is presented and applied to a set of magnetic field configurations. The underlying method for calculation of the island width is that of Cary & Hanson ( Phys. Fluids B, vol. 3, issue 4, 1991, pp. 1006–1014) (with a minor modification), and requires that the residue of the island centre be small. Therefore, the gradient of the residue is calculated in addition. Both the island width and the gradient calculations are verified using an analytical magnetic field configuration introduced by Reiman & Greenside ( Comput. Phys. Commun. , vol. 43, issue 1, 1986, pp. 157–167). The method is also applied to the calculation of the shape gradient of the width of a magnetic island in a National Compact Stellarator Experiment (NCSX) vacuum configuration with respect to positions on a coil. A gradient-based optimization is applied to a magnetic field configuration studied by Hanson & Cary ( Phys. Fluids , vol. 27, issue 4, 1984, pp. 767–769) to minimize stochasticity by adding perturbations to a pair of helical coils. Although only vacuum magnetic fields and an analytical magnetic field model are considered in this work, the adjoint calculation of the island width gradient could also be applied to a magnetohydrodynamic (MHD) equilibrium if the derivative of the magnetic field, with respect to the equilibrium parameters, is known. Using the island width gradient calculation presented here, more general gradient-based optimization methods can be applied to design stellarators with small magnetic islands. Moreover, the sensitivity of the island size may itself be optimized to ensure that coil tolerances, with respect to island size, are kept as high as possible.
Not Available
Identification of chemically similar stars using elemental abundances is core to many pursuits within Galactic archeology. However, measuring the chemical likeness of stars using abundances directly is limited by systematic imprints of imperfect synthetic spectra in abundance derivation. We present a novel data-driven model that is capable of identifying chemically similar stars from spectra alone. We call this relevant scaled component analysis (RSCA). RSCA finds a mapping from stellar spectra to a representation that optimizes recovery of known open clusters. By design, RSCA amplifies factors of chemical abundance variation and minimizes those of nonchemical parameters, such as instrument systematics. The resultant representation of stellar spectra can therefore be used for precise measurements of chemical similarity between stars. We validate RSCA using 185 cluster stars in 22 open clusters in the Apache Point Observatory Galactic Evolution Experiment survey. We quantify our performance in measuring chemical similarity using a reference set of 151,145 field stars. We find that our representation identifies known stellar siblings more effectively than stellar-abundance measurements. Using RSCA, 1.8% of pairs of field stars are as similar as birth siblings, compared to 2.3% when using stellar-abundance labels. We find that almost all of the information within spectra leveraged by RSCA fits into a two-dimensional basis, which we link to [Fe/H] and α-element abundances. We conclude that chemical tagging of stars to their birth clusters remains prohibitive. However, using the spectra has noticeable gain, and our approach is poised to benefit from larger data sets and improved algorithm designs.
First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, power balance analysis has shown that turbulent transport poses a serious limitation to the global performance of the machine. Furthermore, several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of T i ≃ 1.7 keV, except in the few scenarios in which turbulence can be suppressed. In order to understand the precise turbulent mechanisms at play and thus design improved performance scenarios, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. As a first step in this direction, in this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ( k ⊥ ρ i ≃ 1) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients and . Furthermore, we discuss the importance of other potentially relevant parameters such as T e / T i , E r or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio η i = , as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which T i,core > 1.7 keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.
Abstract A simplified view of the space of optimised stellarators has the potential to guide and aid the design efforts of magnetic confinement configurations suitable for future fusion reactors. We present one such view for the class of quasisymmetric stellarators based on their approximate description near their centre (magnetic axis). The result is a space that captures existing designs and presents new ones, providing a common framework to study them. Such a simplified construction offers a basic topological approach, guided by certain theoretical and physical choices, which this paper presents in detail.
A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. Furthermore, the method exploits the Biot–Savart Law's exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value.
The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.
The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.