Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA RADIAL PROFILES”

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 703 records · Page 39

Direct optimization of neoclassical ion transport in stellarator reactors

Abstract We directly optimize stellarator neoclassical ion transport while holding neoclassical electron transport at a moderate level, creating a scenario favorable for impurity expulsion and retaining good ion confinement. Traditional neoclassical stellarator optimization has focused on minimizing ϵ eff , the geometric factor that characterizes the amount of radial transport due to particles in the 1 / ν regime. Under expected reactor-relevant conditions, core electrons will be in the 1 / ν regime and core fuel ions will be in the ν regime. Traditional optimizations thus minimize electron transport and rely on the radial electric field ( E r ) that develops to confine the ions. This often results in an inward-pointing E r that drives high- Z impurities into the core, which may be troublesome in future reactors. In this work, we increase the ratio of the thermal transport coefficients L 11 e / L 11 i , which previous research has shown can create an outward-pointing E r . This effect is very beneficial for impurity expulsion. We obtain self-consistent density, temperature, and E r profiles at reactor-relevant conditions for an optimized equilibrium. This equilibrium is expected to enjoy significantly improved impurity transport properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe CMEs in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approx. 1.5Rs (solar radii), which coincides with the distance at which the Alfven speed profile has a minimum value.We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfven speed peaks (approx. 3Rs - 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approx 1.5Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Solar Energetic Particles↗

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe coronal mass ejections (CMEs) a in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approximately 1.5Rs (solar radii), which coincides with the distance at which the Alfv?n speed profile has a minimum value. We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfv?n speed peaks (?3Rs ? 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approximately 1.5 Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2 Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Gopalswamy, N.↗

Digital holographic interferometry employing Fresnel transform reconstruction for the study of flow shear stabilized Z-pinch plasmas

The ZaP-HD Flow Z-Pinch project provides a platform to explore how shear flow stabilized Z-pinches could scale to high-energy-density plasma (HEDP) and fusion reactor conditions. ZaP-HD generates shear stabilized, axisymmetric Z-pinches with stable lifetimes approaching $60$$\mu$s. The goal of the project is to increase the plasma density and temperature compared to the previous ZaP project by compressing the plasma to smaller radii ($\approx 1$ mm). Radial and axial plasma electron density structure is measured using digital holographic interferometry (DHI), which provides the necessary fine spatial resolution. ZaP-HD's DHI system uses a 2 ns Nd:YAG laser pulse with a second harmonic generator ($\lambda = 532$ nm) to produce holograms recorded by a Nikon D3200 digital camera. The holograms are numerically reconstructed with the Fresnel transform reconstruction method to obtain the phase shift caused by the interaction of the laser beam with the plasma. This provides a two-dimensional map of line-integrated electron density, which can be Abel inverted to determine the local number density. The DHI resolves line-integrated densities down to $3\times 10^{20}$ m$^{-2}$ with spatial resolution near $10 \mu$m. This work presents the first application of Fresnel transform reconstruction as an analysis technique for a plasma diagnostic, and it analyzes the method's accuracy through study of synthetic data. It then presents an Abel inversion procedure that utilizes data on both sides of a Z-pinch local number density profile to maximize profile symmetry. Error estimation and Abel inversion are applied to measured data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High performance double-null plasmas under radiating divertor and mantle scenarios on DIII-D

Experiments to combine a radiating divertor or mantle with high power, high βN plasma operation have successfully reduced the divertor heat flux by 40% but have also encountered several challenging problems. For example, injecting either neon or argon seed impurities have led both to significant fuel dilution of the core plasma and to the emergence of harmful tearing mode activity which compromised plasma beta and energy confinement time. Increased divertor closure ended in more effective control over the injected seed impurity inventory, although it was also observed that small radial variations in the placement of the outer separatrix strike point within the DIII-D ‘closed’ divertor could lead to clear differences in the impurity build-up inside the core plasma. Active particle pumping of argon seed impurities through the divertor leg on the high-field side of a double-null divertor configuration showed little effect on seed impurity inventory and only modest control over the deuterium fuel inventory. Normally, 80% or greater of injected argon was removed by the corresponding cryo-pump on the low-field (outboard) side of that divertor, with the remaining argon being removed by the outboard cryo-pump from the opposite divertor. For single-null divertor cases, the contribution of the inner (high-field) cryo-pump to impurity and main-ion control appears to be complicated by the levels of divertor recycling and degree of divertor detachment. In general, divertor design that impedes the escape of seed impurities, plasma shaping that allows a higher stability limit, and electron cyclotron deposition profiles that prevent tearing mode onset and impede impurity accumulation in the core are under consideration for addressing the issues presented in our work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation into the formation of the scrape-off layer density shoulder in JET ITER-like wall L-mode and H-mode plasmas

The low temperature boundary layer plasma (scrape-off layer or SOL) between the hot core and the surrounding vessel determines the level of power loading, erosion and implantation of material surfaces, and thus the viability of tokamak-based fusion as an energy source. Here, this study explores mechanisms affecting the formation of flattened density profiles, so-called 'density shoulders', in the low-field side (LFS) SOL, which modify ion and neutral fluxes to surfaces—and subsequent erosion. There is evidence against local enhancement of ionization inducing shoulder formation. We find that increases in SOL parallel resistivity, Λ div (=[ L || ν ei Ω i ]/ c s Ω e ), postulated to lead to shoulder growth through changes in SOL turbulence characteristics, correlates with increases in SOL shoulder amplitude, A s , but only under a subset of conditions (D 2 -fuelled L-mode density scans with outer strike point on the horizontal target). Λ div fails to correlate with A s for cases of N 2 seeding or during sweeping of the strike point across the horizontal target. The limited correlation of Λ div and A s is also found for H-mode discharges. Thus, while it may be necessary for Λ div to be above a threshold of ~1 for shoulder formation and/or growth, another mechanism is required. More significantly, we find that in contrast to parallel resistivity, outer divertor recycling, as quantified by the total outer divertor Balmer D α emission, I – D α , does scale with A s where Λ div does and even where Λ div does not. Divertor recycling could lead to SOL density shoulder formation through: (a) reducing the parallel to the field flow (loss) of ions out of the SOL to the divertor; and (b) changes in radial electric fields which lead to E × B poloidal flows as well as potentially affecting SOL turbulence birth characteristics. Therefore, changes in divertor recycling may be the sole process involved in bringing about SOL density shoulders or it may be that it acts in tandem with parallel resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TGLF Recalibration for ITER Standard Case Parameters FY2015: Theory and Simulation Performance Target Final Report

This work was motivated by the observation, as early as 2008, that GYRO simulations of some ITER operating scenarios exhibited nonlinear zonal-flow generation large enough to effectively quench turbulence inside r /a ~ 0.5. This observation of flow-dominated, low-transport states persisted even as more accurate and comprehensive predictions of ITER profiles were made using the state-of-the-art TGLF transport model. This core stabilization is in stark contrast to GYRO-TGLF comparisons for modern-day tokamaks, for which GYRO and TGLF are typically in very close agreement. So, we began to suspect that TGLF needed to be generalized to include the effect of zonal-flow stabilization in order to be more accurate for the conditions of reactor simulations. While the precise cause of the GYRO-TGLF discrepancy for ITER parameters was not known, it was speculated that closeness to threshold in the absence of driven rotation, as well as electromagnetic stabilization, created conditions more sensitive the self-generated zonal-flow stabilization than in modern tokamaks. Need for nonlinear zonal-flow stabilization: To explore the inclusion of a zonal-flow stabilization mechanism in TGLF, we started with a nominal ITER profile predicted by TGLF, and then performed linear and nonlinear GYRO simulations to characterize the behavior at and slightly above the nominal temperature gradients for finite levels of energy transport. Then, we ran TGLF on these cases to see where the discrepancies were largest. The predicted ITER profiles were indeed near to the TGLF threshold over most of the plasma core in the hybrid discharge studied (weak magnetic shear, q > 1). Scanning temperature gradients above the TGLF power balance values also showed that TGLF overpredicted the electron energy transport in the low-collisionality ITER plasma. At first (in Q3), a model of only the zonal-flow stabilization (Dimits shift) was attempted. Although we were able to construct an ad hoc model of the zonal flows that fit the GYRO simulations, the parameters of the model had to be tuned to each case. A physics basis for the zonal flow model was lacking. Electron energy transport at short wavelength: A secondary issue – the high-k electron energy flux – was initially assumed to be independent of the zonal flow effect. However, detailed studies of the fluctuation spectra from recent multiscale (electron and ion scale) GYRO simulations provided a critical new insight into the role of zonal flows. The multiscale simulations suggested that advection by the zonal flows strongly suppressed electron-scale turbulence. Radial shear of the zonal E×B fluctuation could not compete with the large electron-scale linear growth rate, but the k x -mixing rate of the E×B advection could. This insight led to a preliminary new model for the way zonal flows saturate both electron- and ion-scale turbulence. It was also discovered that the strength of the zonal E×B velocity could be computed from the linear growth rate spectrum. The new saturation model (SAT1), which replaces the original model (SAT0), was fit to the multiscale GYRO simulations as well as the ion-scale GYRO simulations used to calibrate the original SAT0 model. Thus, SAT1 captures the physics of both multiscale electron transport and zonal-flow stabilization. In future work, the SAT1 model will require significant further testing and (expensive) calibration with nonlinear multiscale gyrokinetic simulations over a wider variety of plasma conditions – certainly more than the small set of scans about a single C-Mod L-mode discharge. We believe the SAT1 model holds great promise as a physics-based model of the multiscale turbulent transport in fusion devices. Correction to ITER performance predictions: Finally, the impact of the SAT1model on the ITER hybrid case is mixed. Without the electron-scale contribution to the fluxes, the Dimits shift makes a significant improvement in the predicted fusion power as originally posited. Alas, including the high-k electron transport reduces the improvement, yielding a modest net increase in predicted fusion power compared to the TGLF prediction with the original SAT0 model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observational goals for Max '91 to identify the causative agent for impulsive bursts

Recent studies of impulsive hard x ray and microwave bursts suggest that a propagating causative agent with a characteristic velocity of the order of 1000 km/s is responsible for these bursts. The results of these studies are summarized and observable distinguishing characteristics of the various possible agents are highlighted, with emphasis on key observational goals for the Max '91 campaigns. The most likely causative agents suggested by the evidence are shocks, thermal conduction fronts, and propagating modes of magnetic reconnection in flare plasmas. With new instrumentation planned for Max '91, high spatial resolution observations of hard x ray sources have the potential to identify the agent by revealing detailed features of source spatial evolution. Observations with the Very Large Array and other radio imaging instruments are of great importance, as well as detailed modeling of coronal loop structures to place limits on their density and temperature profiles. With the combined hard x ray and microwave imaging observations, aided by loop model results, the simplest causative agent to rule out would be the propagating modes of magnetic reconnection. To fit the observational evidence, reconnection modes would need to travel at approximately the same velocity (the Alfven velocity) in different coronal structures that vary in length by a factor of 10(exp 3). Over such a vast range in loop lengths, it is difficult to believe that the Alfven velocity is constant. Thermal conduction fronts would be suggested by sources that expand along the direction of B and exhibit relatively little particle precipitation. Particle acceleration due to shocks could produce more diverse radially expanding source geometries with precipitation at loop footprints.

Batchelor, D. A.↗

METAL-MASS-TO-LIGHT RATIOS OF THE PERSEUS CLUSTER OUT TO THE VIRIAL RADIUS

We analyzed XMM-Newton data of the Perseus cluster out to {approx}1 Mpc, or approximately half the virial radius. Using the flux ratios of Ly{alpha} lines of H-like Si and S to K{alpha} line of He-like Fe, the abundance ratios of Si/Fe and S/Fe of the intracluster medium (ICM) were derived using the APEC plasma code v2.0.1. The temperature dependence of the line ratio limits the systematic uncertainty in the derived abundance ratio. The Si/Fe and S/Fe in the ICM of the Perseus cluster show no radial gradient. The emission-weighted averages of the Si/Fe and S/Fe ratios outside the cool core are 0.91 {+-} 0.08 and 0.93 {+-} 0.10, respectively, in solar units according to the solar abundance table of Lodders. These ratios indicate that most Fe was synthesized by supernovae Ia. We collected K-band luminosities of galaxies and calculated the ratio of Fe and Si mass in the ICM to K-band luminosity, iron-mass-to-light ratio (IMLR), and silicon-mass-to-light ratio (SMLR). Within {approx}1 Mpc, the cumulative IMLR and SMLR increase with radius. Using Suzaku data for the northwest and east directions, we also calculated the IMLR out to {approx}1.8 Mpc, or about the virial radius. We constrained the SMLR out to this radius and discussed the slope of the initial mass function of stars in the cluster. Using the cumulative IMLR profile, we discuss the past supernova Ia rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Initial operation of a newly developed multichord motional Stark effect diagnostic in KSTAR

A photo-elastic modulator based 25-chord motional Stark effect (MSE) diagnostic has been successfully developed and commissioned in Korea Superconducting Tokamak Advanced Research. The diagnostic measures the radial magnetic pitch angle profile of the Stark splitting of a D-alpha line at 656.1 nm by the electric field associated with the neutral deuterium heating beam. A tangential view of the neutral beam provides a good spatial resolution of 1–3 cm for covering the major radius from 1.74 m to 2.28 m, and the time resolution is achieved at 10 ms. An in-vessel calibration before the vacuum closing as well as an in situ calibration during the tokamak operation was performed by means of specially designed polarized lighting sources. In this work, we present the final design of the installed MSE diagnostic and the first results of the commissioning.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Fast and Accurate Predictions of Radio Frequency Power Deposition and Current Profile via Data-driven Modeling

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. The machine learning models use a database of 16,000+ GENRAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods ensure that the database covers the range of 9 input parameters ($n_{e0}$, $T_{e0}$, $I_p$, $B_t$, $R_0$, $n_{||}$, $Z_{eff}$, $V_{loop}$, $P_{LHCD}$) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to ~ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating Laser Preheat and Applied Magnetic Fields Relevant to the MagLIF Fusion Scheme

The MAGnetized Liner Inertial Fusion (MagLIF) scheme has achieved thermonuclear fusion yields on the Z Facility by imploding a cylindrical liner filled with D 2 fuel that is preheated with a multi-kJ laser and pre-magnetized with an axial Bz=10 T magnetic field. Preheating (T e = 100-200 eV) and pre-magnetizing (10-30 T) the fuel serves to reduce the implosion velocity required to achieve multi-keV fusion-relevant temperatures at stagnation with a modest radial convergence. The challenge of fuel preheat in MagLIF is to deposit multiple-kJ of energy into the underdense (n e /n c <0.1) fuel over ~10 mm target length efficiently and without introducing contaminants. Once the fuel is heated the applied axial magnetic field (ω ce τ e ~ 10) needs to suppress electron thermal conduction sufficiently to prevent unacceptable heat losses to the liner walls. In this LDRD we investigated laser energy deposition at two facilities: The OMEGA-EP laser at the Laboratory for Laser Energetics and the Z-beamlet laser at Sandia National Labs utilizing the PECOS chamber. Multiple experiments were carried out investigating laser transmission through LEH foils, laser heating of underdense gasses and the effects of magnetization on laser preheat. The studies find that magneto-hydrodynamic simulations are able to reproduce energy deposition at MagLIF-like conditions but that at the intensities currently used to preheat MagLIF significant laser plasma instabilities (LPI) occur which partly explain the inability of codes to reproduce previous MagLIF preheat studies. The experiments find that reducing the intensity and smoothing the beam dramatically reduces the amount of stimulated Brillouin backscatter and produces deposition profiles more similar to those produced in simulations. The experiments have provided a large and varied dataset that can be compared to simulations. As part of the LDRD new experimental capabilities have also been developed that will be used to design future MagLIF integrated experiments and investigate fuel magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rotational stability of a long field-reversed configuration

Rotationally driven modes of long systems with dominantly axial magnetic field are considered. We apply the incompressible model and order axial wavenumber small. A recently developed gyro-viscous model is incorporated. A one-dimensional equilibrium is assumed, but radial profiles are arbitrary. The dominant toroidal (azimuthal) mode numbers ℓ=1 and ℓ=2 modes are examined for a variety of non-reversed (B) and reversed profiles. Previous results for both systems with rigid rotor equilibria are reproduced. New results are obtained by incorporation of finite axial wavenumber and by relaxing the assumption of rigid electron and ion rotation. It is shown that the frequently troublesome ℓ=2 field reversed configuration (FRC) mode is not strongly affected by ion kinetic effects (in contrast to non-reversed cases) and is likely stabilized experimentally only by finite length effects. It is also shown that the ℓ=1 wobble mode has a complicated behavior and is affected by a variety of configuration and profile effects. The rotationally driven ℓ=1 wobble is completely stabilized by strong rotational shear, which is anticipated to be active in high performance FRC experiments. Thus, observed wobble modes in these systems are likely not driven by rotation alone.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

X-point position dependence of edge intrinsic toroidal rotation on the Tokamak à Configuration Variable

Recent theoretical work predicts intrinsic toroidal rotation in the tokamak edge to depend strongly on the normalized major radial position of the X-point. With this motivation, we conducted a series of Ohmic L-mode shots on the Tokamak à Configuration Variable, moving the X-point from the inboard to the outboard edge of the last closed flux surface in both lower and upper single null configurations. The edge toroidal rotation evolved from strongly co-current for an inboard X-point to either vanishing or counter-current for an outboard X-point, in agreement with the theoretical expectations. The whole rotation profile shifted roughly rigidly with the edge rotation, resulting in variation of the peak core rotation by more than a factor of two. Core rotation reversals had little effect on the edge rotation. Edge rotation was slightly more counter-current for unfavorable than favorable ∇B drift discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ORIGIN OF THE DENSE CORE MASS FUNCTION IN CONTRACTING FILAMENTS

Mass functions of starless dense cores (CMFs) may arise from contraction and dispersal of core-forming filaments. In an illustrative model, a filament contracts radially by self-gravity, increasing the mass of its cores. During this contraction, FUV photoevaporation and ablation by shocks and winds disperse filament gas and limit core growth. The stopping times of core growth are described by a waiting-time distribution. The initial filament column density profile and the resulting CMF each match recent Herschel observations in detail. Then low-mass cores have short growth ages and arise from the innermost filament gas, while massive cores have long growth ages and draw from more extended filament gas. The model fits the initial density profile and CMF best for mean core density 2 Multiplication-Sign 10{sup 4} cm{sup -3} and filament dispersal timescale 0.5 Myr. Then the typical core mass, radius, mean column density, and contraction speed are respectively 0.8 solar masses, 0.06 pc, 6 Multiplication-Sign 10{sup 21} cm{sup -2}, and 0.07 km s{sup -1}, also in accord with observed values.

79 ASTRONOMY AND ASTROPHYSICS↗