Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stellar plasmas”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enhancement of mid-/high-Z impurity transport by continuous Li-granule dropping in a stellarator plasma

An enhancement of core impurity transport is observed in high-density plasmas of the stellarator large helical device, heated by neutral beam injection, when continuous lithium (Li) granule dropping is performed. In the reported experiments, in which the tracer-encapsulated solid pellet is employed to inject trace amounts of titanium (Ti) and molybdenum (Mo) into the plasma core, confinement times for these impurities are seen to reduce significantly when Li dropping is applied, this reduction being more notable for Mo. To gain some initial insight into these observations, simulations are performed using the drift-kinetic transport code SFINCS for the Mo case. These simulations indicate that, while neoclassical transport prevails for the main plasma components (electrons, majority ions, and low-Z impurities), the classical contribution appears dominant for transporting Mo impurities. In summary, this work reports the first experimental observation of the degradation of mid-Z and high-Z impurity confinement induced by the continuous dropping of Li granules into a high-density stellarator plasma. In the case of the Mo impurity, simulations suggest that classical transport is the key mechanism underlying the enhanced impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗
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 OSTI2022

First lithium experiments in HIDRA and evidence of helium retention during quasi-steady-state stellarator plasma operations

Recent experiments in Hybrid Illinois Device for Research and Applications (HIDRA) have had operational discharges between t discharge = 60 and 1000 s using electron cyclotron resonant heating (ECRH) of the plasma. This means that quasi-steady-state plasma discharges reach conditions to study long-pulse plasma material interactions (PMIs). The newly commissioned HIDRA-Material Analysis Test-stand PMI diagnostic is used to place a drop of lithium onto a heated tungsten surface, transfer the sample in-vacuo and expose it in a helium plasma. Helium is of interest as there is an open question to whether lithium will be able to remove helium ash in real fusion devices. The introduction of the W-Li sample in HIDRA resulted in evaporation of lithium into the helium plasma during a 600 s pulse and caused a reduction of over 90% in neutral pressure during the discharge. It was also observed that the plasma density and temperature increased by over 2.5 times. Using spectroscopy and a helium collisional radiative model, the peak temperature and density of the helium plasma can be monitored during the discharge. During lithium evaporation, as significant lithium ionization occurs, there is a 85% drop in the HIDRA vessel neutral pressure, despite a constant flow rate of He gas. This reduction in neutral pressure is supported by spectroscopy data with corresponding reductions in He I line intensities (587 nm, 667 nm, 706 nm, and 728 nm), as well as those of other impurities. At one point in the discharge a lithium plasma is created, as indicated by an increase in Li + emission and a complete reduction in He + emission, but the electron density jumps from ne = 3 × 10 18 m –3 to over ne = 8 × 10 18 m –3 while the core temperature stays relatively constant between Te = 16 eV and 20 eV. Once lithium has completely evaporated from the sample and the majority of the ionized lithium has diffused from the plasma to the vessel walls, pressure and spectroscopy data paired with He collisional radiative model calculations shows a re-establishment of a helium plasma in a low recycling regime. In this regime, the density drops down to ne = 2 × 10 18 m –3 and the electron temperature increases from T e = 20 eV to over T e = 50 eV indicating an increase in helium heating efficiency. This is also indicated by the He+ emission re-establishing and having a higher intensity. Here, we show the results from the first lithium campaign in HIDRA. In the presence of lithium, and in particular when lithium ions are present, the helium disappears from the plasma via an as of yet unknown complex relationship that needs to be further studied. The most likely explanation is that the lithium ions are distributed around the vessel and able to trap helium to the surface turning HIDRA into a large gettering surface. These results have potential implications on future plasma facing component design using liquid lithium for impurity and recycling control using limiters and divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-C 1 code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-C 1 . A C 1 coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-C 1 extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of a stellarator plasma with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spectroscopic diagnostics of high-temperature plasma in stellar coronae using Fe XXIV-XXIV K-shell lines with XRISM

The RS CVn type binary star GT Muscae was observed during its quiescence using the Resolve X-ray microcalorimeter spectrometer onboard XRISM. The main and satellite lines of the Fe XXIV-XXIV K-shell transitions were resolved for the first time from stellar sources. We conducted line ratio analysis to investigate any deviations from collisional ionization equilibrium and Maxwell electron energy distribution with a single temperature. By using five combinations of direct excitation lines and dielectronic recombination satellite lines in three line complexes (Fe He$\alpha$, Ly$\alpha$, and He$\beta$), we found that the plasma is well characterized by two-temperature thermal plasmas with temperatures of 1.7 and 4.3 keV, which is consistent with thermal broadening of Fe xxv, and the broad-band fitting results in the 1.7–10 keV band. Other forms of deviation from a single-temperature plasma, such as different ionization and electron temperatures or the $\kappa$ distribution for the electron energy distributions, are not favored, which is reasonable for stellar coronae at quiescence. This study demonstrates the utility of the Fe K-shell line ratio diagnostics to probe plasma conditions using X-ray microcalorimeters.

X-rays: stars↗
DOE OSTI2025

Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted‐Chain Approach

Under the extreme conditions found in small stars, where electron degeneracy and Coulomb coupling are significant, accurate modeling of Thomson scattering is crucial for determining opacity, a primary quantity for stellar energy transport. We use hypernetted-chain calculations, incorporating quantum pseudopotentials and electron-exchange effects to obtain the electron–electron static structure factor to calculate the Thomson scattering transport cross-section for conditions prevailing in the interior of small stars. These results are compared to those from average-atom simulations and analytical calculations. Our findings support laboratory astrophysics experiments aimed at benchmarking opacity models for stellar interiors, particularly for red dwarf stars, and help to bridge theoretical models with observations.

ICF↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantum physics of stars

Stars are slowly developing objects; the lifetimes of the different burning phases are determined by the strength of nuclear reactions, which in turn are defined by the quantum structure of the associated nuclei at the threshold and the respective reaction mechanisms. Stars, from the nuclear physics perspective, are cold environments where only a few of the key nuclear reactions have been measured at the actual stellar plasma temperatures. This is also the case for more dynamic astrophysical phenomena from the big bang to stellar explosions. Most of the nuclear reaction rates are therefore based on theoretical extrapolations. A number of discrepancies between these predictions and the associated stellar signatures have been observed, and many may be due to low-energy or near-threshold quantum effects. These effects need to be understood in order to reliably model nuclear reaction processes, not only for stars but also for low-temperature plasma environments such as controlled magnetic or inertial confinement fusion systems, which operate in similar temperature regimes. This review summarizes the various theoretical techniques presently used for deriving reaction rates and discusses possible quantum effects that may impact the reaction cross section near the reaction threshold. These resemble enhanced single-particle and cluster structures near threshold and associated interference effects. New experimental techniques such as deep-underground accelerators or the study of transfer reactions to mimic the quantum-mechanical transition strength, the so-called Trojan horse method, provide ways to directly or indirectly probe the reaction features that determine the reaction rates at stellar energies. Furthermore, this is demonstrated on a number of key nuclear reactions for different nucleosynthesis environments. Finally, current inconsistencies between experimental predictions and observations are discussed.

Models & methods for nuclear reactions↗
DOE OSTI2022

Design of an arrangement of cubic magnets for a quasi-axisymmetric stellarator experiment

The usage of permanent magnets to shape the confining magnetic field of a stellarator has the potential to reduce or eliminate the need for non-planar coils. As a proof-of-concept for this idea, here we have developed a procedure for designing an array of cubic permanent magnets that works in tandem with a set of toroidal-field (TF) coils to confine a stellarator plasma. All of the magnets in the design are constrained to have identical geometry and one of three polarization types in order to simplify fabrication while still producing sufficient field accuracy. We present some of the key steps leading to the design, including the geometric arrangement of the magnets around the device, the procedure for optimizing the polarizations according to the three allowable magnet types, and the choice of magnet types to be used. We apply these methods to design an array of rare-Earth permanent magnets that can be paired with a set of planar TF coils to confine a quasi-axisymmetric plasma with a toroidal magnetic field strength of about 0.5 T on axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Needed computations and computational capabilities for stellarators

Stellarator plasmas are externally controlled to a degree unparalleled by any other fusion concept, magnetic or inertial. This control is largely through the magnetic fields produced by external coils. The development of fusion energy could be expedited by carrying out remarkably straight-forward computations to define strategies for exploiting this external control. In addition to these computations, which have a reliability limited only by competence, certain physics areas that affect the development of stellarator power plants should have a more intense study. The low cost and speed with which computations can be carried out relative to experiments have implications for the development of fusion. Computations should be used to develop a strategy that to the extent possible allows major issues to be circumvented. Required computations for this strategy are the subject of this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Light Curves and Event Rates of Axion Instability Supernovae

Abstract It was recently proposed that exotic particles can trigger a new stellar instability that is analogous to the e − e + pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ∼400 keV–2 MeV and coupling with photons of g a γ ∼ 10 −5 GeV −1 . It is found that the 56 Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10–20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7–2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function.

Mori, Kanji (ORCID:0000000325951657)↗
DOE OSTI2026

High resolution measurements of core density collapse dynamics in LHD

The core density collapse (CDC) is an abrupt instability that degrades the performance of super dense core stellarator plasmas on the Large Helical Device (LHD). In LHD, CDCs have proved to be an operational limit in so-called outward-shifted configurations. New measurements with high time resolution Thomson scattering characterize dynamics of the instability in the core. The density profile is observed to collapse inwards, narrowing the density peak till the core itself collapses. A hot, radially localized feature is also observed in the electron temperature profile during the collapse and disappears within a few hundred microseconds of CDC onset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Greedy permanent magnet optimization

Abstract A number of scientific fields rely on placing permanent magnets in order to produce a desired magnetic field. We have shown in recent work that the placement process can be formulated as sparse regression. However, binary, grid-aligned solutions are desired for realistic engineering designs. We now show that the binary permanent magnet problem can be formulated as a quadratic program with quadratic equality constraints, the binary, grid-aligned problem is equivalent to the quadratic knapsack problem with multiple knapsack constraints, and the single-orientation-only problem is equivalent to the unconstrained quadratic binary problem. We then provide a set of simple greedy algorithms for solving variants of permanent magnet optimization, and demonstrate their capabilities by designing magnets for stellarator plasmas. The algorithms can a-priori produce sparse, grid-aligned, binary solutions. Despite its simple design and greedy nature, we provide an algorithm that compares with or even outperforms the state-of-the-art algorithms while being substantially faster, more flexible, and easier to use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗