Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA MACROINSTABILITIES”

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 73 records · Page 4

H-mode pedestal scaling in DIII-D, ASDEX Upgrade, and JET

Multidevice pedestal scaling experiments in the DIII-D, ASDEX Upgrade (AUG), and JET tokamaks are presented in order to test two plasma physics pedestal width models. The first model proposes a scaling of the pedestal width {Delta}/a {proportional_to} {rho}*{sup 1/2} to {rho}* based on the radial extent of the pedestal being set by the point where the linear turbulence growth rate exceeds the ExB velocity. In the multidevice experiment where {rho}* at the pedestal top was varied by a factor of four while other dimensionless parameters where kept fixed, it has been observed that the temperature pedestal width in real space coordinates scales with machine size, and that therefore the gyroradius scaling suggested by the model is not supported by the experiments. The density pedestal width is not invariant with {rho}* which after comparison with a simple neutral fuelling model may be attributed to variations in the neutral fuelling patterns. The second model, EPED1, is based on kinetic ballooning modes setting the limit of the radial extent of the pedestal region and leads to {Delta}{sub {psi} {proportional_to}} {beta}{sub p}{sup 1/2}. All three devices show a scaling of the pedestal width in normalised poloidal flux as {Delta}{sub {psi} {proportional_to}} {beta}{sub p}{sup 1/2}, as described by the kinetic ballooning model; however, on JET and AUG, this could not be distinguished from an interpretation where the pedestal is fixed in real space. Pedestal data from all three devices have been compared with the predictive pedestal model EPED1 and the model produces pedestal height values that match the experimental data well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parasitic signals in the receiving band of the Sub-Harmonic Arc Detection system on JET ICRF Antennas

When testing the SHAD system on JET ICRF antennas, parasitic signals in the detection band (5-20MHz) were detected. We have identified emission from grid breakdown events in the Neutral Beam injectors, and Ion Cyclotron Emission from the plasma. Spurious signals in the band 4-10 MHz are also often observed at the onset of ELM events. Such parasitic signals could complicate the design and operation of SHAD in ICRF systems for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Results of the implementation on a mock-up of the full 3dB hybrid matching option of the ITER ICRH system

Each of the two ICRH antennas for ITER must couple 20MW to the plasma in the 40-55MHz band via an array of 24 radiating shorted straps fed by four generators. The matching system must provide automatic matching control on the mean load provided by the plasma and be resilient (parallel {Gamma}{sub G} parallel <0.2) to a wide range of fast antenna load excursions occurring in ELMy plasmas. Furthermore, good control of the current distribution in the strap array must be possible for the various heating and current drive scenarios. Two load resilient matching options have been considered for ITER: the 4 'Conjugate-T'(CT) and the 4 hybrids ones, the first being presently considered as a back-up option. Automatic control of these 2 options has been developed, and tested for optimization on a low-powered scaled mock-up. Successful implementation of the simultaneous feedback control of 11 actuators for the matching of the 4 CT and for the control of the toroidal phasing has already been achieved. The matching and the array current control of the 3dB hybrid option are provided by simultaneous feedback control of the decouplers and double stub tuners (in total 23 actuators) and this has also been successfully achieved for the full array. The paper discusses the circuit implementation and presents the obtained results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lower hybrid current drive at plasma densities required for thermonuclear reactors

Driving current in high-density plasmas is essential for the progress of thermonuclear fusion energy research based on the tokamak concept. The lower hybrid current drive (LHCD) effect, is potentially the most suitable tool for driving current at large plasma radii, consistent with the needs of ITER steady state scenario. Unfortunately, experiments at reactor grade high plasma densities with kinetic profiles approaching those required for ITER, have shown problems in penetration of the LH power into the core plasma. These plasmas represent a basic reference for designing possible methods useful for assessing the LHCD concept in ITER. On the basis of the phenomenology observed during LHCD experiments carried out in different machines, and model of the spectral broadening effect due to parametric instability, an interpretation and possible solution of the related important problem is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Individual Sawtooth Pacing by Synchronized ECCD in TCV

Previous real-time sawtooth control scenarios using EC actuators have attempted to shorten or lengthen the sawtooth period by optimally positioning the EC absorption near the q = 1 surface. In new experiments we demonstrate for the first time that individual sawtooth crashes can be repetitively induced at predictable times by reducing the stabilizing ECCD power after a predetermined time from the preceding crash. Other stabilizing actuators (e.g. ICRF, NBI) are expected to produce similar effects. Armed with these results, we present a new sawtooth / NTM control paradigm for improved performance in burning plasmas. The potential appearance of neo-classical tearing modes, triggered by long period sawtooth crashes even at low beta, becomes predictable and therefore amenable to preemptive ECCD. The ITER Electron Cyclotron Upper Launcher (EC-UL) design incorporates the needed functionalities for this method to be applied. The methodology and associated TCV experiments will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CIRCULATION AND DISSIPATION ON HOT JUPITERS

Many global circulation models predict supersonic zonal winds and large vertical shears in the atmospheres of short-period Jovian exoplanets. Using linear analysis and nonlinear local simulations, we investigate hydrodynamic dissipation mechanisms to balance the thermal acceleration of these winds. The adiabatic Richardson criterion remains a good guide to linear stability, although thermal diffusion allows some modes to violate it at very long wavelengths and very low growth rates. Nonlinearly, wind speeds saturate at Mach numbers {approx}2 and Richardson numbers {approx}<1/4 for a broad range of plausible diffusivities and forcing strengths. Turbulence and vertical mixing, though accompanied by weak shocks, dominate the dissipation, which appears to be the outcome of a recurrent Kelvin-Helmholtz instability. An explicit shear viscosity, as well as thermal diffusivity, is added to ZEUS to capture dissipation outside of shocks. The wind speed is neither monotonic nor single valued for a range of shear viscosities larger than about 10{sup -3} of the sound speed times the pressure scale height. Coarsening the numerical resolution can also increase the speed. Hence global simulations that are incapable of representing vertical turbulence and shocks, either because of reduced physics or because of limited resolution, may overestimate wind speeds. We recommend that such simulations include artificial dissipation terms to control the Mach and Richardson numbers and to capture mechanical dissipation as heat.

79 ASTRONOMY AND ASTROPHYSICS↗

THREE-DIMENSIONAL RELATIVISTIC MAGNETOHYDRODYNAMIC SIMULATIONS OF CURRENT-DRIVEN INSTABILITY WITH A SUB-ALFVENIC JET: TEMPORAL PROPERTIES

We have investigated the influence of a velocity shear surface on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria in three dimensions. In this study, we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius, the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than that occuring for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and nonlinear development can be similar to the development of a static plasma column.

79 ASTRONOMY AND ASTROPHYSICS↗

USING TWISTED FILAMENTS TO MODEL THE INNER JET IN M 87

Radio and optical images of the M 87 jet show bright filaments, twisted into an apparent double helix, extending from HST-1 to knot A. Proper motions within the jet suggest a decelerating jet flow passing through a slower, accelerating wave pattern. We use these observations to develop a mass and energy flux-conserving model describing the jet flow and conditions along the jet. Our model requires the jet to be an internally hot, but subrelativistic plasma, from HST-1 to knot A. Subsequently, we assume that the jet is in pressure balance with an external cocoon and we determine the cocoon conditions required if the twisted filaments are the result of the Kelvin-Helmholtz (KH) unstable elliptical mode. We find that the cocoon must be cooler than the jet at HST-1 but must be about as hot as the jet at knot A. Under these conditions, we find that the observed filament wavelength is near the elliptical mode maximum growth rate and growth is rapid enough for the filaments to develop and saturate well before HST-1. We generate a pseudo-synchrotron image of a model jet carrying a combination of normal modes of the KH instability. The pseudo-synchrotron image of the jet reveals that (1) a slow decline in the model jet's surface brightness is still about five times faster than the real jet, (2) KH-produced dual helically twisted filaments can appear qualitatively similar to those on the real jet if any helical perturbation to the jet is very small or nonexistent inside knot A, and (3) the knots in the real jet cannot be associated with the twisted filamentary features and are unlikely to be the result of a KH instability. The existence of the knots in the real jet, the limb brightening of the real jet in the radio, and the slower decline of the surface brightness of the real jet indicate that additional processes-such as unsteady jet flow and internal particle acceleration-are occurring within the jet. Disruption of the real jet beyond knot A by KH instability is consistent with the jet and cocoon conditions we find at knot A.

79 ASTRONOMY AND ASTROPHYSICS↗

THE HELICAL JET OF THE VELA PULSAR

We have studied the fascinating dynamics of the nearby Vela pulsar's nebula in a campaign comprising 11 40 ks observations with the Chandra X-Ray Observatory. The deepest images yet revealed the shape, structure, and motion of the 2 arcmin long pulsar jet. We find that the jet's shape and dynamics are remarkably consistent with those of a steadily turning helix projected on the sky. We discuss possible implications of our results, including free precession of the neutron star and MHD instability scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

EQUILIBRIUM DISKS, MAGNETOROTATIONAL INSTABILITY MODE EXCITATION, AND STEADY-STATE TURBULENCE IN GLOBAL ACCRETION DISK SIMULATIONS

Global three-dimensional magnetohydrodynamic (MHD) simulations of turbulent accretion disks are presented which start from fully equilibrium initial conditions in which the magnetic forces are accounted for and the induction equation is satisfied. The local linear theory of the magnetorotational instability (MRI) is used as a predictor of the growth of magnetic field perturbations in the global simulations. The linear growth estimates and global simulations diverge when nonlinear motions-perhaps triggered by the onset of turbulence-upset the velocity perturbations used to excite the MRI. The saturated state is found to be independent of the initially excited MRI mode, showing that once the disk has expelled the initially net flux field and settled into quasi-periodic oscillations in the toroidal magnetic flux, the dynamo cycle regulates the global saturation stress level. Furthermore, time-averaged measures of converged turbulence, such as the ratio of magnetic energies, are found to be in agreement with previous works. In particular, the globally averaged stress normalized to the gas pressure <{alpha}{sub P}>bar = 0.034, with notably higher values achieved for simulations with higher azimuthal resolution. Supplementary tests are performed using different numerical algorithms and resolutions. Convergence with resolution during the initial linear MRI growth phase is found for 23-35 cells per scale height (in the vertical direction).

79 ASTRONOMY AND ASTROPHYSICS↗

EFFECTS OF RESISTIVITY ON MAGNETIZED CORE-COLLAPSE SUPERNOVAE

We studied the role of turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out two-dimensional resistive-MHD simulations. Three cases with different initial strengths of magnetic field and rotation are investigated: (1) a strongly magnetized rotating core, (2) a moderately magnetized rotating core, and (3) a very strongly magnetized non-rotating core. In each case, one ideal-MHD model and two resistive-MHD models are computed. As a result of these computations, each model shows an eruption of matter assisted by magnetic acceleration (and also by centrifugal acceleration in the rotating cases). We found that resistivity attenuates the explosion in cases 1 and 2, while it enhances the explosion in case 3. We also found that in the rotating cases, the main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of the magnetic field and a twisting of poloidal magnetic field lines by differential rotation, which are somewhat dampened down with the presence of resistivity. Although magnetorotational instability seems to occur in the rotating models, it plays only a minor role in magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matter, >2.5, attained in an ideal-MHD model is reduced to some extent in a resistive model. These results indicate that resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

THE STABILITY OF WEAKLY COLLISIONAL PLASMAS WITH THERMAL AND COMPOSITION GRADIENTS

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport phenomena, and long-term evolution of weakly collisional, magnetized plasmas which are stratified in temperature. The insights gained via these studies have led to a significant improvement of our understanding of the processes that determine the physical evolution and observational properties of the intracluster medium (ICM) permeating galaxy clusters. These studies have been carried out under the assumption that the ICM is a homogeneous medium. This, however, might not be a good approximation if heavy elements are able to sediment in the inner region of the galaxy cluster. Motivated by the need to obtain a more complete picture of the dynamical properties of the ICM, we analyze the stability of a weakly collisional, magnetized plane-parallel atmosphere which is stratified in both temperature and composition. This allows us to discuss for the first time the dynamics of weakly collisional environments where heat conduction, momentum transport, and ion-diffusion are anisotropic with respect to the direction of the magnetic field. We show that depending on the relative signs and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous media. We also find that there are new modes which are driven by heat conduction and particle diffusion. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could be alleviated if the mean molecular weight gradient is steep enough, i.e., ({nabla}{mu})/{mu} > ({nabla}T)/T. This study constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of heavy elements, and its consequences for the long-term evolution and observational signatures of the ICM in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DEAD, UNDEAD, AND ZOMBIE ZONES IN PROTOSTELLAR DISKS AS A FUNCTION OF STELLAR MASS

We investigate the viability of the magnetorotational instability (MRI) in X-ray ionized viscous accretion disks around both solar-type stars and very low mass stars. In particular, we determine the disk regions where the MRI can be shut off either by Ohmic resistivity (the so-called dead and undead zones) or by ambipolar diffusion (a region we term the zombie zone). We consider two stellar masses: M {sub *} = 0.7 M {sub Sun} and 0.1 M {sub Sun }. In each case, we assume that: the disk surface density profile is that of a scaled Minimum Mass Solar Nebula, with M {sub disk}/M {sub *} = 0.01 as suggested by current data; disk ionization is driven primarily by stellar X-rays, complemented by cosmic rays and radionuclides; and the stellar X-ray luminosity scales with bolometric luminosity as L{sub X} /L {sub *} Almost-Equal-To 10{sup -3.5}, as observed. Ionization rates are calculated with the MOCCASIN Monte Carlo X-ray transport code, and ionization balance determined using a simplified chemical network, including well-mixed 0.1 {mu}m grains at various levels of depletion. We find that (1) ambipolar diffusion is the primary factor controlling MRI activity in disks around both solar-type and very low mass classical T Tauri stars. Assuming that the MRI yields the maximum possible field strength at each radius, we further find that: (2) the MRI-active layer constitutes only {approx}5%-10% of the total disk mass; (3) the accretion rate ( M-dot ) varies radially in both magnitude and sign (inward or outward), implying time-variable accretion as well as the creation of disk gaps and overdensities, with consequences for planet formation and migration; (4) achieving the empirical accretion rates in solar-type and very low mass stars requires a depletion of well-mixed small grains (via grain growth and/or settling) by a factor of 10-1000 relative to the standard dust-to-gas mass ratio of 10{sup -2}; and (5) the current non-detection of polarized emission from field-aligned grains in the outer disk regions is consistent with active MRI at those radii.

79 ASTRONOMY AND ASTROPHYSICS↗

TURBULENCE IN THE OUTER REGIONS OF PROTOPLANETARY DISKS. I. WEAK ACCRETION WITH NO VERTICAL MAGNETIC FLUX

We use local numerical simulations to investigate the strength and nature of magnetohydrodynamic (MHD) turbulence in the outer regions of protoplanetary disks, where ambipolar diffusion is the dominant non-ideal MHD effect. The simulations include vertical stratification and assume zero net vertical magnetic flux. We employ a super time-stepping technique to ameliorate the Courant restriction on the diffusive time step. We find that in idealized stratified simulations, with a spatially constant ambipolar Elsasser number Am, turbulence driven by the magnetorotational instability (MRI) behaves in a similar manner as in prior unstratified calculations. Turbulence dies away for Am {<=} 1, and becomes progressively more vigorous as ambipolar diffusion is decreased. Near-ideal MHD behavior is recovered for Am {>=} 10{sup 3}. In the intermediate regime (10 {<=} Am {<=} 10{sup 3}) ambipolar diffusion leads to substantial increases in both the period of the MRI dynamo cycle and the characteristic scales of magnetic field structures. To quantify the impact of ambipolar physics on disk accretion, we run simulations at 30 AU and 100 AU that include a vertical Am profile based upon far-ultraviolet (FUV) ionized disk models. These models develop a vertically layered structure analogous to the Ohmic dead zone that is present at smaller radii. We find that, although the levels of surface turbulence can be strong (and consistent with constraints on turbulent line widths at these radii), the inferred accretion rates are at least an order of magnitude smaller than those observed in T Tauri stars. This discrepancy is very likely due to the assumption of zero vertical magnetic field in our simulations and suggests that vertical magnetic fields are essential for MRI-driven accretion in the outer regions of protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗

ANGULAR MOMENTUM TRANSPORT BY ACOUSTIC MODES GENERATED IN THE BOUNDARY LAYER. II. MAGNETOHYDRODYNAMIC SIMULATIONS

We perform global unstratified three-dimensional magnetohydrodynamic simulations of an astrophysical boundary layer (BL)-an interface region between an accretion disk and a weakly magnetized accreting object such as a white dwarf-with the goal of understanding the effects of magnetic field on the BL. We use cylindrical coordinates with an isothermal equation of state and investigate a number of initial field geometries including toroidal, vertical, and vertical with zero net flux. Our initial setup consists of a Keplerian disk attached to a non-rotating star. In a previous work, we found that in hydrodynamical simulations, sound waves excited by shear in the BL were able to efficiently transport angular momentum and drive mass accretion onto the star. Here we confirm that in MHD simulations, waves serve as an efficient means of angular momentum transport in the vicinity of the BL, despite the magnetorotational instability (MRI) operating in the disk. In particular, the angular momentum current due to waves is at times larger than the angular momentum current due to MRI. Our results suggest that angular momentum transport in the BL and its vicinity is a global phenomenon occurring through dissipation of waves and shocks. This point of view is quite different from the standard picture of transport by a local anomalous turbulent viscosity. In addition to angular momentum transport, we also study magnetic field amplification within the BL. We find that the field is indeed amplified in the BL, but only by a factor of a few, and remains subthermal.

79 ASTRONOMY AND ASTROPHYSICS↗

DYNAMO EFFECTS IN MAGNETOROTATIONAL TURBULENCE WITH FINITE THERMAL DIFFUSIVITY

We investigate the saturation level of hydromagnetic turbulence driven by the magnetorotational instability in the case of vanishing net flux. Motivated by a recent paper of Bodo et al., we here focus on the case of a non-isothermal equation of state with constant thermal diffusivity. The central aim of the paper is to complement the previous result with closure parameters for mean-field dynamo models, and to test the hypothesis that the dynamo is affected by the mode of heat transport. We perform computer simulations of local shearing-box models of stratified accretion disks with approximate treatment of radiative heat transport, which is modeled via thermal conduction. We study the effect of varying the (constant) thermal diffusivity, and apply different vertical boundary conditions (BCs). In the case of impenetrable vertical boundaries, we confirm the transition from mainly conductive to mainly convective vertical heat transport below a critical thermal diffusivity. This transition is however much less dramatic when more natural outflow BCs are applied. Similarly, the enhancement of magnetic activity in this case is less pronounced. Nevertheless, heating via turbulent dissipation determines the thermodynamic structure of accretion disks and clearly affects the properties of the related dynamo. This effect, however, may have been overestimated in previous work, and a careful study of the role played by boundaries will be required.

79 ASTRONOMY AND ASTROPHYSICS↗

A SECOND-ORDER GODUNOV METHOD FOR MULTI-DIMENSIONAL RELATIVISTIC MAGNETOHYDRODYNAMICS

We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second-order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the magnetic field. A variety of approximate Riemann solvers are implemented to compute the fluxes of the conserved variables. The methods are tested with a comprehensive suite of multi-dimensional problems. These tests have helped us develop a hierarchy of correction steps that are applied when the integration algorithm predicts unphysical states due to errors in the fluxes, or errors in the inversion between conserved and primitive variables. Although used exceedingly rarely, these corrections dramatically improve the stability of the algorithm. We present preliminary results from the application of these algorithms to two problems in RMHD: the propagation of supersonic magnetized jets and the amplification of magnetic field by turbulence driven by the relativistic Kelvin-Helmholtz instability (KHI). Both of these applications reveal important differences between the results computed with Riemann solvers that adopt different approximations for the fluxes. For example, we show that the use of Riemann solvers that include both contact and rotational discontinuities can increase the strength of the magnetic field within the cocoon by a factor of 10 in simulations of RMHD jets and can increase the spectral resolution of three-dimensional RMHD turbulence driven by the KHI by a factor of two. This increase in accuracy far outweighs the associated increase in computational cost. Our RMHD scheme is publicly available as part of the Athena code.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETIC KELVIN-HELMHOLTZ INSTABILITY AT THE SUN

Flows and instabilities play a major role in the dynamics of magnetized plasmas including the solar corona, magnetospheric and heliospheric boundaries, cometary tails, and astrophysical jets. The nonlinear effects, multi-scale and microphysical interactions inherent to the flow-driven instabilities, are believed to play a role, e.g., in plasma entry across a discontinuity, generation of turbulence, and enhanced drag. However, in order to clarify the efficiency of macroscopic instabilities in these processes, we lack proper knowledge of their overall morphological features. Here we show the first observations of the temporally and spatially resolved evolution of the magnetic Kelvin-Helmholtz instability in the solar corona. Unprecedented high-resolution imaging observations of vortices developing at the surface of a fast coronal mass ejecta are taken by the new Solar Dynamics Observatory, validating theories of the nonlinear dynamics involved. The new findings are a cornerstone for developing a unifying theory on flow-driven instabilities in rarefied magnetized plasmas, which is important for understanding the fundamental processes at work in key regions of the Sun-Earth system.

79 ASTRONOMY AND ASTROPHYSICS↗