HifiStell: High Fidelity Simulations of Stellarators
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Until now, quasi-isodynamic magnetic fields have been the only known stellarator configurations that, at low collisionality, give small radial neoclassical transport and zero bootstrap current for arbitrary plasma profiles, the latter facilitating control of the magnetic configuration. The recently introduced notion of piecewise omnigenous fields has enormously broadened the space of stellarator configurations with small radial neoclassical transport. Here, in this Letter, the existence of piecewise omnigenous fields that give zero bootstrap current is proven analytically and confirmed numerically. These results establish piecewise omnigenity as an alternative approach to stellarator reactor design.
A three-dimensional computer code has been developed to study the magnetohydrodynamic equilibrium and stability of a diffuse or sharp boundary plasma in toroidal geometry. It is shown how equilibria with net toroidal current identically zero can be determined and how growth rates of instabilities can be calculated. Applications are made to an l = 2, 3 stellarator configuration that offers the possibility of achieving a critical value as high as 10% for the plasma parameter
Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the particle distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through an avalanche. The existence of a mode particle resonance depends on the properties of the equilibrium and particle parameters, and their number, location, and density can vary with device design. Recently, the advent of more powerful computing capabilities and advanced theoretical understanding has led to the design of non-axisymmetric devices or stellarators, which could prove to be more advantageous than tokamaks. Stellarators have the advantage of being immune to major disruptions because of the very low plasma current. One of the problems shared by both types of devices is the existence of resonances in particle orbits, which can lead to large amplitude high frequency instabilities and subsequent induced particle loss. We examine the number of resonances, their location, and dependence on particle energy for some stellarator designs.
Fast magnetic reconnection is defined by the topology of the magnetic field lines changing on a timescale that is approximately an order of magnitude longer than the topology-conserving ideal-evolution timescale. Fast reconnection is an intrinsic property of Faraday's law when the evolving magnetic field depends non-trivially on all three spatial coordinates and is commonly observed—even when the effects that allow topology breaking are arbitrarily small. The associated current density need only be enhanced by a factor of approximately ten and flows in thin but broad ribbons along the magnetic field. These results follow from the variation in the separation of neighboring pairs of magnetic field lines, which in an ideal evolution typically increases exponentially with time, and the existence of a spatial scale below which magnetic field lines freely change their identities due to non-ideal effects, such as resistivity. Traditional reconnection theory ignores exponentially large variations and relies on the current density reaching a magnitude that is exponentially larger than is actually required. Here, an analysis of the behavior of magnetic field lines in the neighborhood of an arbitrarily chosen line is used to obtain more precise and rigorous results on intrinsic reconnection. The maximum parallel kinetic energy of collisionless charged particles is shown to have an exponential increase in time during a generic magnetic evolution.
I am a staff scientist at Sandia National Laboratories (SNL). I work on multiple fundamental-science projects and lead modeling/data-analysis for the stellar opacity experiments. I also work with theorists to refine plasma material-property calculations to make NNSA simulations more predictable. These are challenging and important problems for national security. In fact, our experiments raised questions about “opacity”, the property of matter that controls energy transport inside stars. This work was published by Nature due to its serious implication over broad applications, and our team received NNSA Defense Program Award of Excellence, both in 2015. Work at SNL is collegial and full of learning through interactions with the world’s finest experimentalists and theorists. When I started my Ph.D. program, working at a national lab was a dream position. Today, I collaborate with top scientists to solve mission-critical problems. The Stewardship Science Academic Programs (SSAP) and National Laser Users’ Facility (NLUF) played a vital role in the efficient learning and seamless transition to SNL.
Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.
The deposit of nonradiative heat and momentum in the mantle of a hot star affects the interpretation of the stellar spectrum in two ways. First, a superheated and moving plasma should be considered when doing the analysis, and second, a model atom which is appropriate for the physical state of the line forming regions. Some examples are presented for H and He showing how the changes in the electron temperature affect the solution of the equations of statistical equilibrium. The observed spectra of the Wolf-Raynet stars HD 191765, HD 192103, and HD 192163 are compatible with a normal H/He abundance ratio.
The radio emission from dMe flare stars is discussed using Very Large Array and Arecibo observations as examples. Active flare stars emit weak, unpolarized, quiescent radio radiation that may be always present. Although thermal bremsstrahlung and/or thermal gyroresonance radiation account for the slowly-varying, quiescent radio radiation of solar active regions, these processes cannot account for the long-wavelength quiescent radiation observed from nearby dMe flare stars. It has been attributed to nonthermal gyrosynchrotron radiation, but some as yet unexplained mechanism must be continually producing the energetic electrons. Long duration, narrow-band radiation is also emitted from some nearby dMe stars at 20 cm wavelength. Such radiation may be attributed to coherent plasma radiation or to coherent electron-cyclotron masers. Impulsive stellar flares exhibit rapid variations that require radio sources that are smaller than the star in size, and high brightness temperatures greater than 10(exp 15) K that are also explained by coherent radiation processes. Quasi-periodic temporal fluctuations suggest pulsations during some radio flares. Evidence for frequency structure and positive or negative frequency drifts during radio flares from dMe stars is also presented.
The radio emission from dMe flare stars is discussed using Very Large Array and Arecibo observations as examples. Active flare stars emit weak, unpolarized, quiescent radio radiation that may be always present. Although thermal bremsstrahlung and/or thermal gyroresonance radiation account for the slowly-varying, quiescent radio radiation of solar active regions, these processes cannot account for the long-wavelength quiescent radiation observed from nearby dMe flare stars. It has been attributed to nonthermal gyrosynchrotron radiation, but some as yet unexplained mechanism must be continually producing the energetic electrons. Long duration, narrow-band radiation is also emitted from some nearby dMe stars at 20 cm wavelength. Such radiation may be attributed to coherent plasma radiation or to coherent electron-cyclotron masers. Impulsive stellar flares exhibit rapid variations that require radio sources that are smaller than the star in size, and high brightness temperatures greater than 10(exp 15) K that are also explained by coherent radiation processes. Quasi-periodic temporal fluctuations suggest pulsations during some radio flares. Evidence for frequency structure and positive or negative frequency drifts during radio flares from dMe stars is also presented.
A numerical integration is undertaken of the integral expression for the thermal conductivity due to collisions of relativistic degenerate electrons and compared to the limiting-case analytic formula. The integration is designed to encompass all temperature/electron-plasma-frequency temperature ratios between the melting temperature and the Fermi temperature. High accuracy fits are demonstrated by interpolating the values of the integrals of the function and by using an asymptotic function by Urpin and Yakovlev (1980). The numerically integrated expression compares favorably to the limiting-case analytic asymptotic formula by Urpin and Yakovlev, and the results are valid for temperatures above and below the electron-plasma-frequency temperature. The present techniques can be used in stellar opacity calculations and in the study of the propagation of deflagration fronts of compact stellar remnants.
We report on variations in important X-ray emission lines in a series of Chandra grating spectra of the supermassive colliding wind binary star eta Car, including key phases around the X-ray minimum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of eta Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission lines provide the only direct measure of the flow dynamics of the companion's wind along the wind-wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with symmetric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant velocity shifts and correlated increases in line widths through the observations. The R = forbidden-to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit in all observations, suggesting that the line-forming region is > 1.6 stellar radii from the companion star, and that the emitting plasma may be in a non-equilibrium state. We show that simple geometrical models cannot simultaneously fit both the observed centroid variations and changes in line width as a function of phase. We show that the observed profiles can be fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision boundary. We use this analysis to help constrain the line formation region as a function of orbital phase, and the orbital geometry. Subject headings: X-rays: stars -stars: early-type-stars: individual (q Car)
The HST Treasury Program Advanced Spectral Library Project: Cool Stars was designed to collect representative, high-quality UV spectra of eight evolved FM type cool stars. The Space Telescope Imaging Spectrograph (STIS)echelle spectra of these objects enable investigations of a broad range of topics, including stellar and interstellar astrophysics. This paper provides a guide to the spectra of the two evolved M stars, the M2 Iab supergiant Oriand the M3.4 giant Cru, with comparisons to the prototypical K1.5 giant Boo. It includes identifications of the significant atomic and molecular emission and absorption features and discusses the character of the photospheric and chromospheric continua and line spectra. The fluorescent processes responsible for a large portion of the emission-line spectrum, the characteristics of the stellar winds, and the available diagnostics for hot and cool plasmas are also summarized. This analysis will facilitate the future study of the spectra, outer atmospheres, and winds, not only of these objects but of numerous other cool, low-gravity stars, for years to come.
The HST (Hubble Space Telescope) Treasury Program Advanced Spectral Library Project: Cool Stars was designed to collect representative, high-quality ultraviolet spectra of eight evolved F-M type cool stars. The Space Telescope Imaging Spectrograph (STIS) echelle spectra of these objects enable investigations of a broad range of topics, including stellar and interstellar astrophysics. This paper provides a guide to the spectra of the two evolved M stars, the M2 Iab super giant alpha Ori and the M3.4 giant gamma Cru, with comparisons to the prototypical K1.5 giant alpha Boo. It includes identifications of the significant atomic and molecular emission and absorption features and discusses the character of the photospheric and chromospheric continua and line spectra. The fluorescent processes responsible for a large portion of the emission-line spectrum, the characteristics of the stellar winds, and the available diagnostics for hot and cool plasmas are also summarized. This analysis will facilitate the future study of the spectra, outer atmospheres, and winds, not only of these objects but of numerous other cool, low-gravity stars, for years to come.
In recent years many efforts have been undertaken to simplify coil designs for stellarators due to the difficulties in fabricating non-planar coils. Additionally, the FOCUS code removes the need for a winding surface and represents the coils as arbitrary curves in 3D. In the following work, the implementation of a spline representation for the coils in FOCUS is described, along with the implementation of a new engineering constraint to design coils with a straighter outer section. The new capabilities of the code are shown as an example on HSX, NCSX, and a prototype quasi-axisymmetric reactor-sized stellarator. The flexibility granted by splines along with the new constraint will allow for stellarator coil designs with improved accessibility and simplified maintenance.
Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.
We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.
This project supported a U.S. contribution to the TJ-II stellarator research program focused on distinguishing radially directed particle transport from energy transport in magnetized fusion plasmas. The motivating physics issue is that turbulent particle flux and turbulent heat or energy flux are not necessarily locked together: changes in density, electron temperature, plasma potential, and electric-field fluctuations can produce different phase relationships and therefore different radial fluxes. Resolving these relationships is important for understanding improved-confinement behavior, transport-barrier-like regimes, and the broader ability to predict and control confinement in stellarators and related toroidal devices.