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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

The Soft State of Cygnus X-1 Observed with NuSTAR: A Variable Corona and a Stable Inner Disk

We present a multi-epoch hard X-ray analysis of Cygnus X-1 in its soft state based on four observations with the Nuclear Spectroscopic Telescope Array (NuSTAR). Despite the basic similarity of the observed spectra, there is clear spectral variability between epochs. To investigate this variability, we construct a model incorporating both the standard disk-corona continuum and relativistic reflection from the accretion disk, based on prior work on Cygnus X-1, and apply this model to each epoch independently. We find excellent consistency for the black hole spin and the iron abundance of the accretion disk, which are expected to remain constant on observational timescales. In particular, we confirm that Cygnus X-1 hosts a rapidly rotating black hole, 0.93 < approx. a* < approx. 0.96, in broad agreement with the majority of prior studies of the relativistic disk reflection and constraints on the spin obtained through studies of the thermal accretion disk continuum. Our work also confirms the apparent misalignment between the inner disk and the orbital plane of the binary system reported previously, finding the magnitude of this warp to be approx.10deg-15deg. This level of misalignment does not significantly change (and may even improve) the agreement between our reflection results and the thermal continuum results regarding the black hole spin. The spectral variability observed by NuSTAR is dominated by the primary continuum, implying variability in the temperature of the scattering electron plasma. Finally, we consistently observe absorption from ionized iron at approx. 6.7 keV, which varies in strength as a function of orbital phase in a manner consistent with the absorbing material being an ionized phase of the focused stellar wind from the supergiant companion star.

black hole physics↗
DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Magnetic Reconnection as the Driver of the Solar Wind

We present EUV solar observations showing evidence for omnipresent jetting activity driven by small-scale magnetic reconnection at the base of the solar corona. We argue that the physical mechanism that heats and drives the solar wind at its source is ubiquitous magnetic reconnection in the form of small-scale jetting activity (i.e., a.k.a. jetlets). This jetting activity, like the solar wind and the heating of the coronal plasma, are ubiquitous regardless of the solar cycle phase. Each event arises from small-scale reconnection of opposite polarity magnetic fields producing a short-lived jet of hot plasma and Alfv´en waves into the corona. The discrete nature of these jetlet events leads to intermittent outflows from the corona, which homogenize as they propagate away from the Sun and form the solar wind. This discovery establishes the importance of small-scale magnetic reconnection in solar and stellar atmospheres in understanding ubiquitous phenomena such as coronal heating and solar wind acceleration. Based on previous analyses linking the switchbacks to the magnetic network, we also argue that these new observations might provide the link between the magnetic activity at the base of the corona and the switchback solar wind phenomenon. These new observations need to be put in the bigger picture of the role of magnetic reconnection and the diverse form of jetting in the solar atmosphere.

magnetic reconnection↗
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 OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Ground-based observations of comets, the Jupiter plasma Torus, and Io

Aspects of cometary and magnetospheric physics were investigated by means of ground-based astronomical spectroscopy. High-throughput, dual-etalon Fabry-Perot spectrometers were used to obtain very high resolution spectra of atomic, molecular, and ionic emission lines from the diffuse gases and plasmas associated with comets and the Jupiter plasma torus. The Fabry-Perot spectrometers were also used with a charge coupled device (CCD) camera to obtain images of these extended emission sources in individual spectral lines at high spectral resolution. A new program using the McMath solar-stellar spectrograph to observe emission lines from Io was recently initiated. The McMath spectrograph has a high resolution mode which allows the detection of narrow, relatively faint emission lines superimposed on Io's reflected solar spectrum.

Scherb, Frank↗
NASA NTRS2012

X-Ray Emission from an Asymmetric Blast Wave and a Massive White Dwarf in the Gamma Ray Emitting Nova V407 CYG

Classical nova events in symbiotic stars, although rare, offer a unique opportunity to probe the interaction between ejecta and a dense environment in stellar explosions. In this work, we use X-ray data obtained with Swift and Suzaku during the recent classical nova outburst in V407 Cyg to explore such an interaction. We find evidence of both equilibrium and non-equilibrium ionization plasmas at the time of peak X-ray brightness, indicating a strong asymmetry in the density of the emitting region. Comparing a simple model to the data, we find that the X-ray evolution is broadly consistent with nova ejecta driving a forward shock into the dense wind of the Mira companion. We detect a highly absorbed soft X-ray component in the spectrum during the first 50 days of the outburst that is consistent with supersoft emission from the nuclear burning white dwarf. The high temperature and short turn off time of this emission component, in addition to the observed breaks in the optical and UV lightcurves, indicate that the white dwarf in the binary is extremely massive. Finally, we explore the connections between the X-ray and GeV-ray evolution, and propose that the gamma ray turn-off is due to the stalling of the forward shock as the ejecta reach the red giant surface.

Nelson, Thomas↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
NASA NTRS2000

Emission Line Spectra from Low-Density Laboratory Plasmas

Using spectroscopic equipment optimized for laboratory astrophysics, we are performing systematic measurements of the line emission from astrophysically relevant ions in the wavelength band between 1 and 400 Angstroms important to X-ray missions such as Chandra, XMM, Astro-E, and EUVE. Obtained in a controlled laboratory setting at electron densities similar to those found in stellar coronae, the data are used to test spectral modeling codes for accuracy and completeness. Our effort includes the compilation of the iron L-shell emission lines from 6-18 Angstroms and the iron M-shell emission lines from 50-200 Angstroms. Many lines have been identified for the first time, and the fluxes from lines missing in the spectral modeling codes are assessed. Our measurements also assess the accuracy of line excitation calculations, including direct electron-impact excitation, dielectronic recombination, and resonance excitation. These measurements yield a calibration of specific diagnostic line ratios. Examples of our current measurements are given.

Beiersdorfer, P.↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗
DOE OSTI2021

Laser blow-off with photodiode detection system at the helically symmetric experiment

Abstract A laser blow-off (LBO) system has been installed in the HSX stellarator to investigate impurity transport. A Nd:YAG laser is used to ablate a target area on a glass slide that has been coated with a thin film of aluminum and rapidly deposits a small, controlled quantity of aluminum atoms into HSX plasma. To study the radial propagation and confinement of the injected aluminum impurities, time-resolved emission measurements are made using a 20 channel AXUV photodiode detection system. The LBO system has achieved its design purpose, most notably the ability to make controlled impurity injection without significant perturbation to the background plasma. A detailed description of the setup as well as initial results of operation are presented.

Instruments & Instrumentation↗
NASA NTRS1991

On two-dimensional relativistic stellar winds

Stellar wind is of great interest, because many of the astrophysical systems possess wind like structures. For a long time only nonrelativistic winds have been studied. However, recently it was proposed by Kennel et al. that the relativistic pulsar wind with the plasma, consisting of electrons and positrons, can be responsible for the observed features of the Crab nebula. The study of Kennel et al. revealed the inconsistency of the assumption of the wind zero temperature with the observational data. It has been shown that only for high relativistic temperatures can high Mach numbers be reached, which allows the possibility of a shock formation. The analysis of Kennel et al. was extended into the 2-D case by using conventional magnetohydrodynamic equations for a relativistic plasma with an isotropic relativistic temperature. The state equation is assumed polytropic.

Gedalin, M. E.↗
NASA NTRS1992

Scientific objectives of a Solar Probe mission

The Solar Probe is a mission of exploration and discovery in an uninvestigated region near the sun. The Solar Probe can enhance general understanding of the corona and the solar wind. The first in situ measurements of coronal particles and fields are expected to provide definitive discriminators among many currently proposed models of processes that shape the structure and dynamics of the outer corona and solar wind. In addition, observations of the development of turbulence, nonlinear wave processes, plasma heating, and particle acceleration in the corona can advance these areas of plasma physics in a regime neither duplicated by earth-based laboratories nor by previous space exploration. It also fills the last gap between knowledge of the terrestrial response to solar output as well as provide insight for interpreting observations of general stellar envelopes.

Feldman, W.↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Computer simulation of plasma and N-body problems

Any spiral structure in computer-generated galaxies is generally short lived and the final state is a rotating bar. The bar thus obtained rotates more slowly than the stars. It has been argued that core/halo components have a stabilizing effect on galaxies and result in longer lived spiral structure. However, numerical experiments with large fixed stellar components representing the core/halo component show that multiarmed spiral structure develops and persists for many rotations but only in an evolving manner. That is, the spiral structure is either wound up into a tight pattern or it is wound up and then reappears again. A recent study of the effect of fixed core/halo components does show that the bar instability is indeed inhibited by a sufficiently large fixed component. The present study determines the effect of a self-consistent (rather than fixed) core/halo component in order to show whether there are any instabilities (such as two-stream) or other important interactions present that may be suppressed with a fixed core. Also studied were the effects of finite thickness of the disk and of three-dimensional essentially spherical core/halo components.

Harries, W. L.↗
DOE OSTI2026

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Orbital and short-term variability of iron absorption

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe K band on second timescales. Using data from the Performance Verification phase in 2024 April, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser, and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve’s capability to probe plasma environments in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2023

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗