Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stellar”

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 595 records · Page 33

Collisionless microinstabilities in stellarators. I. Analytical theory of trapped-particle modes

This is the first in a series of papers about collisionless, electrostatic micro-instabilities in stellarators, with an emphasis on trapped-particle modes. It is found that, in so-called maximum-J configurations, trapped-particle instabilities are absent in large regions of parameter space. Quasi-isodynamic stellarators have this property (approximately), and the theory predicts that trapped electrons are stabilizing to all eigenmodes with frequencies below the electron bounce frequency. The physical reason is that the bounce-averaged curvature is favorable for all orbits, and that trapped electrons precess in the direction opposite to that in which drift waves propagate, thus precluding wave-particle resonance. These considerations only depend on the electrostatic energy balance and are independent of all geometric properties of the magnetic field other than the maximum-J condition. However, if the aspect ratio is large and the instability phase velocity differs greatly from the electron and ion thermal speeds, it is possible to derive a variational form for the frequency showing that stability prevails in a yet larger part of parameter space than what follows from the energy argument. Collisionless trapped-electron modes should therefore be more stable in quasi-isodynamic stellarators than in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of continuum damping of Alfvén eigenmodes in tokamak and stellarator equilibria

In an ideal magnetohydrodynamic (MHD) plasma, shear Alfvén eigenmodes may experience dissipationless damping due to resonant interaction with the shear Alfvén continuum. This continuum damping can make a significant contribution to the overall growth/decay rate of shear Alfvén eigenmodes, with consequent implications for fast ion transport. One method for calculating continuum damping is to solve the MHD eigenvalue problem over a suitable contour in the complex plane, thereby satisfying the causality condition. Such an approach can be implemented in three-dimensional ideal MHD codes which use the Galerkin method. Analytic functions can be fitted to numerical data for equilibrium quantities in order to determine the value of these quantities along the complex contour. This approach requires less resolution than the established technique of calculating damping as resistivity vanishes and is thus more computationally efficient. The complex contour method has been applied to the three-dimensional finite element ideal MHD Code for Kinetic Alfvén waves. In this paper, we discuss the application of the complex contour technique to calculate the continuum damping of global modes in tokamak as well as torsatron, W7-X and H-1NF stellarator cases. To the authors' knowledge, these stellarator calculations represent the first calculation of continuum damping for eigenmodes in fully three-dimensional equilibria. The continuum damping of global modes in W7-X and H-1NF stellarator configurations investigated is found to depend sensitively on coupling to numerous poloidal and toroidal harmonics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical investigation of non-perturbative kinetic effects of energetic particles on toroidicity-induced Alfvén eigenmodes in tokamaks and stellarators

The resonant interaction of shear Alfvén waves with energetic particles is investigated numerically in tokamak and stellarator geometry using a non-perturbative MHD-kinetic hybrid approach. The focus lies on toroidicity-induced Alfvén eigenmodes (TAEs), which are most easily destabilized by a fast-particle population in fusion plasmas. While the background plasma is treated within the framework of an ideal-MHD theory, the drive of the fast particles, as well as Landau damping of the background plasma, is modelled using the drift-kinetic Vlasov equation without collisions. Building on analytical theory, a fast numerical tool, STAE-K, has been developed to solve the resulting eigenvalue problem using a Riccati shooting method. The code, which can be used for parameter scans, is applied to tokamaks and the stellarator Wendelstein 7-X. High energetic-ion pressure leads to large growth rates of the TAEs and to their conversion into kinetically modified TAEs and kinetic Alfvén waves via continuum interaction. To better understand the physics of this conversion mechanism, the connections between TAEs and the shear Alfvén wave continuum are examined. It is shown that, when energetic particles are present, the continuum deforms substantially and the TAE frequency can leave the continuum gap. The interaction of the TAE with the continuum leads to singularities in the eigenfunctions. To further advance the physical model and also to eliminate the MHD continuum together with the singularities in the eigenfunctions, a fourth-order term connected to radiative damping has been included. The radiative damping term is connected to non-ideal effects of the bulk plasma and introduces higher-order derivatives to the model. Thus, it has the potential to substantially change the nature of the solution. For the first time, the fast-particle drive, Landau damping, continuum damping, and radiative damping have been modelled together in tokamak- as well as in stellarator geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The shape of LITTLE THINGS dwarf galaxies DDO 46 and DDO 168: Understanding the stellar and gas kinematics

Determining the shape of dwarf irregular (dIrr) galaxies is controversial because if one assumes that these objects are disks and if these disks are randomly distributed over the sky, then their projected minor-to-major axis ratios should follow a particular statistical distribution, which is not observed. Thus, different studies have led to different conclusions. Some believe that the observed distributions can be explained by assuming the dIrrs are thick disks while others have concluded that dIrrs are triaxial. Fortunately, the central stellar velocity dispersion, σ{sub z,0}, combined with maximum rotation speed, V{sub max}, provides a kinematic measure, V{sub max}/σ{sub z,0}, which gives the three-dimensional shape of a system. In this work, we present the stellar and gas kinematics of DDO 46 and DDO 168 from the Local Irregulars That Trace Luminosity Extremes; The H i Nearby Galaxy Survey (LITTLE THINGS) and determine their respective V{sub max}/σ{sub z,0} values. We used the Kitt Peak National Observatory's Mayall 4 m telescope with the Echelle spectrograph as a long-slit spectrograph, which provided a two-dimensional, 3′-long slit. We acquired spectra of DDO 168 along four position angles (PAs) by placing the slit over the morphological major and minor axes and two intermediate PAs. However, due to poor weather conditions during our observing run for DDO 46, we were able to extract only one useful data point from the morphological major axis. We determined a central stellar velocity dispersion perpendicular to the disk, σ{sub z,0}, of 13.5 ± 8 km s{sup −1} for DDO 46 and 〈σ{sub z,0}〉 of 10.7 ± 2.9 km s{sup −1} for DDO 168. We then derived the maximum rotation speed in both galaxies using the LITTLE THINGS H i data. We separated bulk motions from non-circular motions using a double Gaussian decomposition technique and applied a tilted-ring model to the bulk velocity field. We corrected the observed H i rotation speeds for asymmetric drift and found a maximum velocity, V{sub max}, of 77.4 ± 3.7 and 67.4 ± 4.0 for DDO 46 and DDO 168, respectively. Thus, we derived a kinematic measure, V{sub max}/σ{sub z,0}, of 5.7 ± 0.6 for DDO 46 and 6.3 ± 0.3 for DDO 168. Comparing these values to ones determined for spiral galaxies, we find that DDO 46 and DDO 168 have V{sub max}/σ{sub z,0} values indicative of thin disks, which is in contrast to minor-to-major axis ratio studies.

79 ASTRONOMY AND ASTROPHYSICS↗

EVIDENCE FOR TWO DISTINCT STELLAR INITIAL MASS FUNCTIONS: REVISITING THE EFFECTS OF CLUSTER DYNAMICAL EVOLUTION

We measure the velocity dispersions of six galactic globular clusters using spatially integrated spectra, to test for the effects of internal dynamical evolution in the stellar mass-to-light ratios, Y{sub *}, of star clusters. In particular, we revisit whether the low values of Y{sub *} that we found in our previous study, from which we concluded that there are at least two population of stellar clusters with distinct stellar initial mass functions, are artificially depressed by relaxation driven mass loss. The combination of our previous sample of five old clusters and these six now provide an order of magnitude range in cluster mass with which to explore this issue. We find no relationship between cluster mass, or relaxation time, and Y{sub *}. Because relaxation is mass dependent, we conclude that the values of Y{sub *} for these clusters are not strongly affected by dynamical effects, and so confirm the presence of the population of clusters with low Y{sub *}.

79 ASTRONOMY AND ASTROPHYSICS↗

YOUNG STELLAR CLUSTERS WITH A SCHUSTER MASS DISTRIBUTION. I. STATIONARY WINDS

Hydrodynamic models for spherically symmetric winds driven by young stellar clusters with a generalized Schuster stellar density profile are explored. For this we use both semi-analytic models and one-dimensional numerical simulations. We determine the properties of quasi-adiabatic and radiative stationary winds and define the radius at which the flow turns from subsonic to supersonic for all stellar density distributions. Strongly radiative winds significantly diminish their terminal speed and thus their mechanical luminosity is strongly reduced. This also reduces their potential negative feedback into their host galaxy interstellar medium. The critical luminosity above which radiative cooling becomes dominant within the clusters, leading to thermal instabilities which make the winds non-stationary, is determined, and its dependence on the star cluster density profile, core radius, and half-mass radius is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DEVIATION OF STELLAR ORBITS FROM TEST PARTICLE TRAJECTORIES AROUND SGr A* DUE TO TIDES AND WINDS

Monitoring the orbits of stars around Sgr A* offers the possibility of detecting the precession of their orbital planes due to frame dragging, of measuring the spin and quadrupole moment of the black hole, and of testing the no-hair theorem. Here we investigate whether the deviations of stellar orbits from test-particle trajectories due to wind mass loss and tidal dissipation of the orbital energy compromise such measurements. We find that the effects of stellar winds are, in general, negligible. On the other hand, for the most eccentric orbits (e > 0.96) for which an optical interferometer, such as GRAVITY, will detect orbital plane precession due to frame dragging, the tidal dissipation of orbital energy occurs at timescales comparable to the timescale of precession due to the quadrupole moment of the black hole. As a result, this non-conservative effect is a potential source of systematic uncertainty in testing the no-hair theorem with stellar orbits.

79 ASTRONOMY AND ASTROPHYSICS↗

Burying a binary: Dynamical mass loss and a continuous optically thick outflow explain the candidate stellar merger V1309 Scorpii

V1309 Sco was proposed to be a stellar merger and a common envelope transient based on the pre-outburst light curve of a contact eclipsing binary with a rapidly decaying orbital period. Using published data, I show that the period decay timescale P/ P-dot of V1309 Sco decreased from ∼1000 to ∼170 yr in ≲ 6 yr, which implies a very high value of P{sup ¨}. I argue that V1309 Sco experienced an onset of dynamical mass loss through the outer Lagrange point, which eventually obscured the binary. The photosphere of the resulting continuous optically thick outflow expands as the mass-loss rate increases, explaining the ∼200 day rise to optical maximum. The model yields the mass-loss rate of the binary star as a function of time and fits the observed light curve remarkably well. It is also possible to observationally constrain the properties of the surface layers undergoing the dynamical mass loss. V1309 Sco is thus a prototype of a new class of stellar transients distinguished by a slow rise to optical maximum that are driven by dynamical mass loss from a binary. I discuss the implications of these findings for stellar transients and other suggested common envelope events.

79 ASTRONOMY AND ASTROPHYSICS↗

The metastable dynamo model of stellar rotational evolution

This paper introduces a new empirical model for the rotational evolution of Sun-like stars—those with surface convection zones and non-convective interior regions. Previous models do not match the morphology of observed (rotation period)-color diagrams, notably the existence of a relatively long-lived 'C-sequence' of fast rotators first identified by Barnes. This failure motivates the Metastable Dynamo Model (MDM) described here. The MDM posits that stars are born with their magnetic dynamos operating in a mode that couples very weakly to the stellar wind, so their (initially very short) rotation periods at first change little with time. At some point, this mode spontaneously and randomly changes to a strongly coupled mode, the transition occurring with a mass-dependent lifetime that is of the order of 100 Myr. I show that with this assumption, one can obtain good fits to observations of young clusters, particularly for ages of 150-200 Myr. Previous models and the MDM both give qualitative agreement with the morphology of the slower-rotating 'I-sequence' stars, but none of them have been shown to accurately reproduce the stellar-mass-dependent evolution of the I-sequence stars, especially for clusters older than a few hundred million years. I discuss observational experiments that can test aspects of the MDM, and speculate that the physics underlying the MDM may be related to other situations described in the literature, in which stellar dynamos may have a multi-modal character.

79 ASTRONOMY AND ASTROPHYSICS↗

OUTWARD MOTION OF POROUS DUST AGGREGATES BY STELLAR RADIATION PRESSURE IN PROTOPLANETARY DISKS

We study the dust motion at the surface layer of protoplanetary disks. Dust grains in the surface layer migrate outward owing to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study we calculate the mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates, in contrast to homogeneous and spherical compact grains, for which radiation pressure efficiency becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 μm is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum-rich inclusions are hardly transported to the outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles are efficiently transported to the comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from the inner hot region to the outer cold cometary region in the protosolar nebula.

79 ASTRONOMY AND ASTROPHYSICS↗

STELLAR ACTIVITY AND ITS IMPLICATIONS FOR EXOPLANET DETECTION ON GJ 176

We present an in-depth analysis of stellar activity and its effects on radial velocity (RV) for the M2 dwarf GJ 176 based on spectra taken over 10 yr from the High Resolution Spectrograph on the Hobby-Eberly Telescope. These data are supplemented with spectra from previous observations with the HIRES and HARPS spectrographs, and V- and R-band photometry taken over six years at the Dyer and Fairborn observatories. Previous studies of GJ 176 revealed a super-Earth exoplanet in an 8.8-day orbit. However, the velocities of this star are also known to be contaminated by activity, particularly at the 39-day stellar rotation period. We have examined the magnetic activity of GJ 176 using the sodium I D lines, which have been shown to be a sensitive activity tracer in cool stars. In addition to rotational modulation, we see evidence of a long-term trend in our Na I D index, which may be part of a long-period activity cycle. The sodium index is well correlated with our RVs, and we show that this activity trend drives a corresponding slope in RV. Interestingly, the rotation signal remains in phase in photometry, but not in the spectral activity indicators. We interpret this phenomenon as the result of one or more large spot complexes or active regions which dominate the photometric variability, while the spectral indices are driven by the overall magnetic activity across the stellar surface. In light of these results, we discuss the potential for correcting activity signals in the RVs of M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Linking stellar coronal activity and rotation at 500 MYR: A deep Chandra observation of M37

Empirical calibrations of the stellar age-rotation–activity relation (ARAR) rely on observations of the co-eval populations of stars in open clusters. We used the Chandra X-ray Observatory to study M37, a 500-Myr-old open cluster that has been extensively surveyed for rotation periods (P{sub rot}). M37 was observed almost continuously for five days, for a total of 440.5 ks, to measure stellar X-ray luminosities (L{sub X}), a proxy for coronal activity, across a wide range of masses. The cluster’s membership catalog was revisited to calculate updated membership probabilities from photometric data and each star’s distance to the cluster center. The result is a comprehensive sample of 1699 M37 members: 426 with P{sub rot}, 278 with X-ray detections, and 76 with both. We calculate Rossby numbers, R{sub o}= P{sub rot}/τ, where τ is the convective turnover time, and ratios of the X-ray-to-bolometric luminosity, L{sub X}/L{sub bol}, to minimize mass dependencies in our characterization of the rotation-coronal activity relation at 500 Myr. We find that fast rotators, for which R{sub o}<0.09±0.01, show saturated levels of activity, with log(L{sub X}/L{sub bol})= −3.06±0.04. For R{sub o}⩾0.09±0.01, activity is unsaturated and follows a power law of the form R{sub o}{sup β}, where β = −2.03{sub −0.14}{sup +0.17}. This is the largest sample available for analyzing the dependence of coronal emission on rotation for a single-aged population, covering stellar masses in the range 0.4–1.3 M{sub ⊙}, P{sub rot} in the range 0.4–12.8 days, and L{sub X} in the range 10{sup 28.4−30.5} erg s{sup −1}. Our results make M37 a new benchmark open cluster for calibrating the ARAR at ages of ≈500 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating fusion alpha heating in a stellarator reactor

Gyrocenter following simulations of fusion born alpha particles in a stellarator reactor are preformed using the BEAMS3D code. The Wendelstein 7-X high mirror configuration is scaled in geometry and magnetic field to reactor relevant parameters. A $2\,\times\,10^{20}m^{-3}$ density plasma with 20 keV core temperatures is assumed and fusion birth rates calculated for various fusion products assuming a 50/50 deuterium-tritium mixture. It is found that energetic He 4 ions comprise the vast majority of the energetic particle inventory. Slowing down simulations of the He4 population suggest plasma heating consistent with scaled energy confinement times for a stellarator reactor. Losses for this configuration appear large suggesting optimization beyond the scope of the W7-X device is key to a future fusion reactor. These first simulations are designed to demonstrate the capability of the BEAMS3D code to provide fusion alpha birth and heating profiles for stellarator reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Including the vacuum energy in stellarator coil design

Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New method to design stellarator coils without the winding surface

Finding an easy-to-build coils set has been a critical issue for stellarator design for decades. Conventional approaches assume a toroidal 'winding' surface, but a poorly chosen winding surface can unnecessarily constrain the coil optimization algorithm, This article presents a new method to design coils for stellarators. Each discrete coil is represented as an arbitrary, closed, one-dimensional curve embedded in three-dimensional space. A target function to be minimized that includes both physical requirements and engineering constraints is constructed. The derivatives of the target function with respect to the parameters describing the coil geometries and currents are calculated analytically. A numerical code, named flexible optimized coils using space curves (FOCUS), has been developed. Furthermore, applications to a simple stellarator configuration, W7-X and LHD vacuum fields are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Stellarator coil design using cubic splines for improved access on the outboard side

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic particle optimization of quasi-axisymmetric stellarator equilibria

Abstract An important goal of stellarator optimization is to achieve good confinement of energetic particles such as, in the case of a reactor, alphas created by deuterium–tritium fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining this improved QA while minimizing the analytical quantity Γ C , which represents the angle between magnetic flux surfaces and contours of J | | , the second adiabatic invariant. This was performed multiple times, resulting in a group of equilibria with significantly reduced energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator have successfully been reduced by optimizing Γ C . The relationship between energetic particle losses and metrics such as QA error ( E q a ) and Γ C in this set of equilibria were examined via statistical methods and a nearly linear relationship between volume-averaged Γ C and prompt particle losses was found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗