Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stellarator”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Improved performance of stellarator coil design optimization

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic well and Mercier stability of stellarators near the magnetic axis

We have recently demonstrated that by expanding in small distance from the magnetic axis compared with 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 the 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 and Johnson's stability criterion for resistive modes approximately coincides with Mercier's ideal condition.

Physics↗
DOE OSTI2021

Constrained stellarator coil curvature optimization with FOCUS

Finding less complicated coils that have adequately low field errors is a crucial step in stellarator development. One coil metric that is of high importance is the maximum curvature of the coil centreline, or coil single filament. Conductors cannot be bent below some threshold minimum radius of curvature. High coil curvatures can cause strains to exceed acceptable levels, especially in superconducting coils. We investigate three ways to optimize coil curvature and find that applying penalty functions to the coil curvature solves for coils that have a constrained maximum curvature and low field error. Penalty functions are implemented in FOCUS and coil solutions optimized for an HSX-like ‘plasma boundary’ are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗
DOE OSTI2021

The on-axis magnetic well and Mercier's criterion for arbitrary stellarator geometries

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.

Physics↗
DOE OSTI2022

Measures of quasisymmetry for stellarators

Quasisymmetric stellarators are an attractive class of optimised magnetic confinement configurations. The property of quasisymmetry (QS) is in practice limited to be approximate, and thus the construction requires measures that quantify the deviation from the exact property. In this paper we study three measure candidates used in the literature, placing the focus on their origin and a comparison of their forms. The analysis shows clearly the lack of universality in these measures. As these metrics do not directly correspond to any physical property (except when exactly quasisymmetric), optimisation should employ additional physical metrics for guidance. It is suggested that close to QS minima, one should treat QS metrics through inequality constraints so that additional physics metrics dominate optimisation. The impact of different quasisymmetric measures on optimisation is presented through an example, for which the standard metric that weights the asymmetric Fourier modes of the field magnitude appears to perform best.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

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 OSTI2024

An asymptotic Grad–Shafranov equation for quasisymmetric stellarators

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 $\beta$ of the order of the ratio 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 is the set of 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.

Nikulsin, Nikita (ORCID:0000000318611777)↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2020

High-resolution X-ray spectroscopy of the stellar wind in Vela X-1 during a flare

In this work, we present a ~130 ks observation of the prototypical wind-accreting, high-mass X-ray binary Vela X-1 collected with XMM-Newton at orbital phases between 0.12 and 0.28. A strong flare took place during the observation that allows us to investigate the reaction of the clumpy stellar wind to the increased X-ray irradiation. To examine the wind’s reaction to the flare, we performed both time-averaged and time-resolved analyses of the RGS spectrum and examined potential spectral changes. We focused on the high-resolution XMM-Newton RGS spectra and divided the observation into pre-flare, flare, and post-flare phases. We modeled the time-averaged and time-resolved spectra with phenomenological components and with the self-consistent photoionization models calculated via CLOUDY and XSTAR in the pre-flare phase, where strong emission lines due to resonant transitions of highly ionized ions are seen. In the spectra, we find emission lines corresponding to K-shell transitions in highly charged ions of oxygen, neon, magnesium, and silicon as well as radiative recombination continua (RRC) of oxygen. Additionally, we observe potential absorption lines of magnesium at a lower ionization stage and features identified as iron L lines. The, CLOUDY and XSTAR photoionization models provide contradictory results, either pointing towards uncertainties in theory or possibly a more complex multi-phase plasma, or both. We are able to demonstrate the existence of a plethora of variable narrow features, including the firm detection of oxygen lines and RRC that RGS enables to observe in this source for the first time. We show that Vela X-1 is an ideal source for future high-resolution missions, such as XRISM and Athena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Once in a blue stream: Detection of recent star formation in the NGC 7241 stellar stream with MEGARA

Aims.In this work we study the striking case of a narrow blue stream with a possible globular cluster-like progenitor around the NGC 7241 galaxy and its foreground dwarf companion. We want to figure out if the stream was generated by tidal interaction with NGC 7241 or if it first interacted with the foreground dwarf companion and later both fell together toward NGC 7241. Methods.We used four sets of observations, including a follow-up spectroscopic study of this stream based on data taken with the MEGARA instrument at the 10.4 m Gran Telescopio Canarias using the integral field spectroscopy mode, the Mount Lemmon 0.80 m telescope, the Telescopio NazionaleGalileo, the DESI Imaging Legacy surveys, and GALEX archival data. We also used high-resolution zoomed-in cosmological simulations. Results.Our data suggest that the compact object we detected in the stream is a foreground Milky Way halo star. Near this compact object we detect emission lines overlapping a less compact, bluer, and fainter blob of the stream that is clearly visible in both ultraviolet and optical deep images. From its heliocentric systemic radial velocity derived from the [O III]λ5007 Å lines (V syst = 1548.58 ± 1.80 km s −1 ) and new UV and optical broadband photometry, we conclude that this overdensity could be the actual core of the stream, with an absolute magnitude ofM g ∼ −10 and ag − r = 0.08 ± 0.11, consistent with a remnant of a low-mass dwarf satellite undergoing a current episode of star formation. From the width of the stream and assuming a circular orbit, we calculate that the progenitor mass can be typical of a dwarf galaxy, but it could also be substantially lower if the stream is on a very radial orbit or if it was created by tidal interaction with the companion dwarf instead of with NGC 7241. These estimates also suggest that this is one of the lowest mass streams detected to date beyond the Local Group. Finally, we find that blue stellar streams containing star formation regions are commonly predicted by high-resolution cosmological simulations of galaxies lighter than the Milky Way. This scenario is consistent with the processes explaining the bursty star formation history of some dwarf satellites, which are followed by a gas depletion and a fast quenching once they enter within the virial radius of their host galaxies for the first time. Thus, it is likely that the stream’s progenitor is undergoing a star formation burst comparable to those that have shaped the star formation history of several Local Group dwarfs in the last few gigayears.

Astronomy & Astrophysics↗
DOE OSTI2024

Re-evaluation of the sequential 3 α reaction in stellar conditions

The formation of carbon via the triple alpha process is a key step in stellar production reactions of light ions and heavily impacts the life-cycle of red giant stars. Here, we propose a novel theoretical approach for sequential mechanism of the helium burning process, based on the Imaginary Time Method (ITM) and semiclassical models. We compare our results to those of the NACRE collaboration, which has set the standard for the past decades. Our calculations differ by several orders of magnitude in the low temperature region, which may seriously affect the early evolution of helium burning stars.

Depastas, Theodoros↗
DOE OSTI2021

Characterization of injection and confinement improvement through impurity induced profile modifications on the Wendelstein 7-X stellarator

Pulsed injections of boron carbide granules into Wendelstein 7-X stellarator (W7-X) plasmas transiently increase the plasma stored energy and core ion temperatures above the reference W7-X experimental programs by up to 30%. In a series of 4 MW electron cyclotron resonance heating experiments, the PPPL Probe Mounted Powder Injector provided 50 ms bursts of 100 μm granules every 350 ms at estimated quantities ranging from approximately 1 mg/pulse to over 30 mg/pulse. For each injection, the stored energy was observed to initially drop and the radiated power transiently increased, while the radial electron density profile rose at the edge as material was assimilated. Once the injected boron carbide was fully absorbed, the density rise transitioned to the core while the stored energy increased above the previous baseline level by an amount linearly correlated with the injection quantity. During the injection, the ion temperature gradient steepened with peak core ion temperatures observed to increase from a nominal 1.7 keV to over 2.6 keV for the largest injection amounts. Enhanced performance is accompanied by a reversal of the radial electric field at ρ < 0.3, indicating that the core transport has switched to the ion root. These observations are suggestive of a change in transport and provide further evidence that externally induced profile modifications provide a possible path to enhanced W7-X performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Generalized Mercier stability criterion for stellarators

The Mercier criterion is a well-known stability criterion for tokamaks. It is derived from a 2 × 2 matrix eigenvalue problem arising from the expansion of resonant solutions about a singular surface where m−nq=0, with m and n being the poloidal and toroidal mode numbers, respectively, and q being the safety factor. The stability criterion is that the eigenvalues must be real, otherwise, the solution oscillates, violating the Newcomb crossing criterion. Because of the non-axisymmetry of stellarators, different toroidal as well as poloidal harmonics couple to each other. It follows that each singular surface can have multiple resonant harmonics, with multiplicity M≥1. The corresponding matrix eigenvalue problem involves a 2M×2M matrix, resulting in M pairs of positive and negative eigenvalues. The generalized stability criterion is that all eigenvalues must be real. While the original Mercier criterion can be expressed in terms of quadratures of equilibrium quantities over the singular surface, which can be evaluated anywhere, the generalized Mercier criterion can only be evaluated on rational q surfaces with a given set of resonant harmonics.

Physics↗
DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗
DOE OSTI2024

Localized phase contrast imaging at the Wendelstein 7-X stellarator

In its basic form, phase contrast imaging (PCI) provides line-integrated measurements of electron density fluctuations in plasmas. As turbulent fluctuations in magnetically confined plasmas have wave vectors almost perpendicular to the background magnetic field, the signals scattered by fluctuations from different parts of the PCI line-of-sight (LoS) are spatially separated in focal planes of the plasma. This allows localized PCI measurements by placing a mask in such a plane, to only permit signals from specific parts of the LoS to reach the PCI detectors. The present paper describes modeling and design of localization masks for the PCI system at the Wendelstein 7-X (W7-X) stellarator as well as the first results obtained using the masks in the recent long-pulse W7-X experimental campaign. During this project, we have extended the theory describing the mask response within the Fraunhofer diffraction model. As a novel development, we show from first principles that the mask response is determined by the fraction of power of the scattered beam spots that passes the mask. These insights have been used to select the W7-X mask design, consisting of a circular cutout, allowing the unscattered beam spot to pass the mask, with wedges covering a fixed angular range outside the central cutout. In the recent W7-X experimental campaign, the masks have verified the location of the main turbulence features observed by the PCI system and provided new information about the location of short-wavelength magnetohydrodynamic modes.

Physics↗
DOE OSTI2024

A Spitzer Model-C Stellarator fusion reactor

Tokamaks are affected by two serious problems. The first problem is the radiation damage to the internal reactor structures by energetic neutrons carrying 80% of the total fusion energy that are generated by the deuterium–tritium fusion cycle that tokamaks are planning to use. The second problem is a loss of plasma control because of the onset of instabilities, causing the sudden shutdown of the tokamak itself. Certainly, no electric utility knowing about these disruptive events would dare to use a tokamak reactor for powering the turbines of its power plant! This paper describes a new fusion reactor based on the Model-C Stellarator that would not be affected by these fatal phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗