Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLOSED PLASMA DEVICES”

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 415 records · Page 23

Ideal ballooning modes in the tokamak scrape-off layer

A drift-reduced Braginskii fluid model is used to carry out a linear and non-linear study of ideal ballooning modes in the tokamak scrape-off layer. First, it is shown that the scrape-off layer finite connection length and boundary conditions modify the ideal stability limit with respect to the closed flux-surface result. Then, in a two-fluid description, it is found that magnetic induction effects can destabilize long wavelength resistive ballooning modes below marginal ideal stability. Non-linear simulations confirm a gradual transition from small scale quasi-electrostatic interchange turbulence to longer wavelength modes as the plasma beta is increased. The transition to global ideal ballooning modes occurs, roughly, at the linearly obtained stability threshold. The transport levels and the pressure gradient as a function of plasma beta obtained in non-linear simulations can be predicted using the non-linear flattening of the pressure profile from the linear modes as a turbulent saturation mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic diagnostics for equilibrium reconstructions with eddy currents on the lithium tokamak experiment

The Lithium Tokamak eXperiment is a spherical tokamak with a close-fitting low-recycling wall composed of thin lithium layers evaporated onto a stainless steel-lined copper shell. Long-lived non-axisymmetric eddy currents are induced in the shell and vacuum vessel by transient plasma and coil currents and these eddy currents influence both the plasma and the magnetic diagnostic signals that are used as constraints for equilibrium reconstruction. A newly installed set of re-entrant magnetic diagnostics and internal saddle flux loops, compatible with high-temperatures and lithium environments, is discussed. Details of the axisymmetric (2D) and non-axisymmetric (3D) treatments of the eddy currents and the equilibrium reconstruction are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics of forced magnetic reconnection in coaxial helicity injection experiments in National Spherical Torus Experiment

We numerically examine the physics of fast flux closure in transient coaxial helicity injection (CHI) experiments in National Spherical Torus Experiment (NSTX). By performing resistive Magnetohydrodynamics (MHD) simulations with poloidal injector coil currents held constant in time, we find that closed flux surfaces are formed through forced magnetic reconnection. Through a local Sweet-Parker type reconnection with an elongated current sheet in the injector region, closed flux surfaces expand in the NSTX global domain. Simulations demonstrate outflows approaching poloidally Alfvénic flows and reconnection times consistent with the Sweet-Parker model. Critical requirements for magnetic reconnection and flux closure are studied in detail. These primary effects, which are magnetic diffusivity, injector flux, injector flux footprint width, and rate of injector voltage reduction, are simulated for transient CHI experiments. The relevant time scales for effective reconnection are τ{sub V}<τ{sub rec}≈τ{sub A}√(S)(1+Pm){sup 1/4}<τ{sub R}, where τ{sub V} is the time for the injector voltage reduction, τ{sub A} is the poloidal Alfvén transit time, τ{sub R} is the global resistive diffusion time, and Pm and S are Prandtl and Lundquist numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic neoclassical transport in the H-mode pedestal

Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrape-off layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. Quantitative agreement between the flux-driven simulation and the experimental electron density, impurity density, and orthogonal measurements of impurity temperature and flow profiles is achieved by adding random-walk particle diffusion to the guiding-center drift motion. The radial electric field (E{sub r}) that maintains ambipolar transport across flux surfaces and to the wall is computed self-consistently on closed and open magnetic field lines and is in excellent agreement with experiment. The E{sub r} inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport are primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-I{sub p} parallel flows in the pedestal, ion temperature anisotropy, and large impurity temperatures in the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detachment scalings derived from 1D scrape-off-layer simulations

Fusion power plants will require detachment to mitigate sputtering and keep divertor heat fluxes at tolerable levels. Controlling detachment on these devices may require the use of real-time scrape-off-layer modeling to complement the limited set of available diagnostics. In this work, we use the configurable Hermes-3 edge modeling framework to perform time-dependent, fixed-fraction-impurity 1D detachment simulations. Although currently far from real-time, these simulations are used to investigate time-dependent effects and the minimum physics set required for control-relevant modeling. We show that these simulations reproduce the expected rollover of the target ion flux — a typical characteristic of detachment onset. We also perform scans of the input heat flux and impurity concentration and show that the steady-state results closely match the scalings predicted by the 0D time-independent Lengyel–Goedheer model. This allows us to indirectly compare to SOLPS simulations, which find a similar scaling but a lower value for the impurity concentration required for detachment for given upstream conditions. We use this result to suggest a series of improvements for the Hermes simulations, and finally show simulations demonstrating the impact of time-dependence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High density LHRF experiments in Alcator C-Mod and implications for reactor scale devices

Parametric decay instabilities (PDI) appear to be an ubiquitous feature of lower hybrid current drive (LHCD) experiments at high density. In density ramp experiments in Alcator C-Mod and other machines the onset of PDI activity has been well correlated with a decrease in current drive efficiency and production of fast electron bremsstrahlung. However whether PDI is the primary cause of the 'density limit', and if so by exactly what mechanism (beyond the obvious one of pump depletion) has not been clearly established. In order to further understand the connection, the frequency spectrum of PDI activity occurring during Alcator C-Mod LHCD experiments has been explored in detail by means of a number of RF probes distributed around the periphery of the C-Mod tokamak including a probe imbedded in the inner wall. The results show that (i) the excited spectra consists mainly of a few discrete ion cyclotron (IC) quasi-modes, which have higher growth than the ion sound branch; (ii) PDI activity can begin either at the inner or outer wall, depending on magnetic configuration; (iii) the frequencies of the IC quasi-modes correspond to the magnetic field strength close to the low-field side (LFS) or high-field side separatrix; and (iv) although PDI activity may initiate near the inner separatrix, the loss in fast electron bremsstrahlung is best correlated with the appearance of IC quasi-modes characteristic of the magnetic field strength near the LFS separatrix. These data, supported by growth rate calculations, point to the importance of the LFS scrape-off layer (SOL) density in determining PDI onset and degradation in current drive efficiency. By minimizing the SOL density it is possible to extend the core density regime over which PDI can be avoided, thus potentially maximizing the effectiveness of LHCD at high density. Increased current drive efficiency at high density has been achieved in FTU and EAST through lithium coating and special fuelling methods, and in recent C-Mod experiments by operating at higher plasma current. Another approach would be to locate the launcher in the inner wall with double null operation. This would reduce the SOL density by an order of magnitude or more and greatly mitigate the effects of PDI as well as other parasitic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of high-performance long-pulse discharge in KSTAR

High-performance long-pulse plasma operation is essential for producing economically viable fusion energy in tokamak devices. To achieve such discharges in KSTAR, firstly, the rapid increase in the temperature of plasma-facing components was mitigated. The temperature increase of the poloidal limiter, especially, was associated with beam-driven fast ion orbit loss and the discrepancy of the equilibrium reconstructed with heated magnetic probes of signal drift. The fast ions lost to the poloidal limiter were reduced by optimizing the plasma shape and the composition of neutral beam injection (NBI). This nonlinear signal drift was successfully reduced by a new thermal shielding protector on the magnetic probes. Secondly, a lower loop voltage approach was implemented to reduce a poloidal flux consumption rate. A plasma current of 400 kA and a line-averaged electron density of ~2.0 × 10 19 m –3 were chosen by considering the L–H power threshold, fast ion orbit loss, and beam shine-through power loss for low loop voltage in KSTAR. In addition, the application of electron cyclotron heating also helped maintain the plasma with low loop voltage (~25 mV) by enhancing the NBI-driven current and achieving a high poloidal beta (β P ) state. KSTAR has achieved a long pulse (~90 s) operation with the high performance of β P $\leqslant$ 2.7, thermal energy confinement enhancement factor (H 98y2 ) ~ 1.1, and fraction of non-inductive current (f NI ) $\leqslant$ 0.96. Still, gradual degradation of the plasma performance has been observed over time in the discharges. In one of the long-pulse discharges, β P reduced by ~18% over the time of ~8τ R (current relaxation time, τ R ~5 s) and ~1067τ E,th (thermal energy confinement time, τ E,th ~45 ms). The degradation may be closely associated with weak, yet growing, and persistent toroidal Alfvén eigenmodes and their effect on fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling, measurement, and 3-D equilibrium reconstruction of the bootstrap current in the Helically Symmetric Experiment

The bootstrap current for three electron cyclotron resonance heated plasma scenarios in a quasihelically symmetric stellarator (the Helically Symmetric Experiment) are analyzed and compared to a neoclassical transport code PENTA. The three conditions correspond to 50 kW input power with a resonance that is off-axis, 50 kW on-axis heating and 100 kW on-axis heating. When the heating location was moved from off-axis to on-axis with 50 kW heating power, the stored energy and the extrapolated steady-state current were both observed to increase. When the on-axis heating power was increased from 50 kW to 100 kW, the stored energy continued to increase while the bootstrap current slightly decreased. This trend is qualitatively in agreement with the calculations which indicate that a large positive electric field for the 100 kW case was driving the current negative in a small region close to the magnetic axis and accounting for the decrease in the total integrated current. This trend in the calculations is only observed to occur when momentum conservation between particle species is included. Without momentum conservation, the calculated bootstrap current increases monotonically. We show that the magnitude of the bootstrap current as calculated by PENTA agrees better with the experiment when momentum conservation between plasma species is included in the calculation. The total current was observed in all cases to flow in a direction to unwind the transform, unlike in a tokamak in which the bootstrap current adds to the transform. The 3-D inductive response of the plasma is simulated to predict the evolution of the current profile during the discharge. The 3-D equilibrium reconstruction code V3FIT is used to reconstruct profiles of the plasma pressure and current constrained by measurements with a set of magnetic diagnostics. The reconstructed profiles are consistent with the measured plasma pressure profile and the simulated current profile when the reconstruction is constrained by the measured data from a diagnostic array that is internal to the vacuum chamber.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisionless microinstabilities in stellarators. III. The ion-temperature-gradient mode

We investigate the linear theory of the ion-temperature-gradient (ITG) mode, with the goal of developing a general understanding that may be applied to stellarators. We highlight the Wendelstein 7X (W7-X) device. Simple fluid and kinetic models that follow closely from existing literature are reviewed and two new first-principle models are presented and compared with results from direct numerical simulation. One model investigates the effect of regions of strong localized shear, which are generic to stellarator equilibria. These “shear spikes” are found to have a potentially significant stabilizing affect on the mode; however, the effect is strongest at short wavelengths perpendicular to the magnetic field, and it is found to be significant only for the fastest growing modes in W7-X. A second model investigates the long-wavelength limit for the case of negligible global magnetic shear. The analytic calculation reveals that the effect of the curvature drive enters at second order in the drift frequency, confirming conventional wisdom that the ITG mode is slab-like at long wavelengths. Using flux tube simulations of a zero-shear W7-X configuration, we observe a close relationship to an axisymmetric configuration at a similar parameter point. It is concluded that scale lengths of the equilibrium gradients constitute a good parameter space to characterize the ITG mode. Thus, to optimize the magnetic geometry for ITG mode stability, it may be fruitful to focus on local parameters, such as the magnitude of bad curvature, connection length, and local shear at locations of bad curvature (where the ITG mode amplitude peaks)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Start current of dielectric-loaded grating in Smith-Purcell radiation

In this paper, a three-dimensional dielectric loaded grating (DLG) is proposed for the Smith-Purcell (SP) device. Taking into the considerations of thickness and width of electron beam, the dispersion equation is derived by using field matches method. The complex frequency is obtained by the numerical solution of dispersion equation, in which the imaginary part represents linear growth rate. The impacts of the electron beam filling factor (EBFF) on growth rate are discussed under the condition that the beam current and beam current density are kept as constants, respectively. In addition, the start current for SP oscillator is obtained by using the dispersion relation combined with boundary conditions. The relationship between the start current and other parameters is discussed and compared with the conventional metal grating. The results show that with the increasing of EBFF, the peak growth rate increases rapidly firstly and then decreases slowly, in which the current and current density are kept as constants, respectively. For the SP oscillator, the start current is increased with the shifting up beam voltage, but it is decreased with the improved EBFF, and only it has a slightly increasing trend when EBFF is close to 1. In addition, the start current is decreased with the increasing of relative dielectric constant, which indicates that by introducing DLG, the start current can be effectively reduced. Theoretical results are in good agreement with that of the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The transition from field emission to collisional space-charge limited current with nonzero initial velocity

Abstract Multiple electron emission mechanisms often contribute in electron devices, motivating theoretical studies characterizing the transitions between them. Previous studies unified thermionic and field emission, defined by the Richardson-Laue-Dushman (RLD) and Fowler–Nordheim (FN) equations, respectively, with the Child-Langmuir (CL) law for vacuum space-charge limited current (SCLC); another study unified FN and CL with the Mott-Gurney (MG) law for collisional SCLC. However, thermionic emission, which introduces a nonzero injection velocity, may also occur in gas, motivating this analysis to unify RLD, FN, CL, and MG. We exactly calculate the current density as a function of applied voltage over a range of injection velocity (i.e., temperature), mobility, and gap distance. This exact solution approaches RLD, FN, and generalized CL (GCL) and MG (GMG) for nonzero injection velocity under appropriate limits. For nonzero initial velocity, GMG approaches zero for sufficiently small applied voltage and mobility, making these gaps always space-charge limited by either GMG at low voltage or GCL at high voltage. The third-order nexus between FN, GMG, and GCL changes negligibly from the zero initial velocity calculation over ten orders of magnitude of applied voltage. These results provide a closed form solution for GMG and guidance on thermionic emission in a collisional gap.

42 ENGINEERING↗

Spontaneous onset of magnetic reconnection in toroidal plasma caused by breaking of 2D symmetry

Magnetic reconnection is studied in the collisionless limit at the Versatile Toroidal Facility (VTF) at MIT. Two distinct magnetic configurations are applied in the experiments; an open magnetic cusp and a closed cusp. In the open cusp configurations, the field lines intersect the the vacuum vessel walls and here axisymmetric oscillatory reconnection is observed. Meanwhile, in the closed cusp configuration, where the field lines are confined inside the experiment, the coupling between global modes and a current sheet leads to powerful bursts of 3D spontaneous reconnection. These spontaneous events start at one toroidal location, and then propagate around the toroidal direction at the Alfven speed (calculated with the strength of the dominant guide field). The three dimensional measurements include the detailed time evolution of the plasma density, current density, the magnetic flux function, the electrostatic potential, and the reconnection rate. The vastly different plasma behavior in the two configurations can be described using a simple theoretical framework, linking together the interdependencies of the reconnection rate, the in-plane electrostatic potential, and the parallel electron currents. We find that it is the breaking of toroidal symmetry by the global mode that allows for a localized disruption of the x-line current and hereby facilitates the onset of spontaneous reconnection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Profile stiffness measurements in the Helically Symmetric experiment and comparison to nonlinear gyrokinetic calculations

Stiffness measurements are presented in the quasi-helically symmetric experiment (HSX), in which the neoclassical transport is comparable to that in a tokamak and turbulent transport dominates throughout the plasma. Electron cyclotron emission is used to measure the local electron temperature response to modulated electron cyclotron resonant heating. The amplitude and phase of the heat wave through the steep electron temperature gradient (ETG) region of the plasma are used to determine a transient electron thermal diffusivity that is close to the steady-state diffusivity. The low stiffness in the region between 0.2 ≤ r/a ≤ 0.4 agrees with the scaling of the steady-state heat flux with temperature gradient in this region. These experimental results are compared to gyrokinetic calculations in a flux-tube geometry using the gyrokinetic electromagnetic numerical experiment code with two kinetic species. Linear simulations show that the ETG mode may be experimentally relevant within r/a ≤ 0.2, while the Trapped Electron Mode (TEM) is the dominant long-wavelength microturbulence instability across most of the plasma. The TEM is primarily driven by the density gradient. Non-linear calculations of the saturated heat flux driven by the TEM and ETG bracket the experimental heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards identifying the mechanisms underlying field-aligned edge-loss of HHFW power on NSTX

Fast-wave heating will be a major heating scheme on ITER, as it can heat ions directly and is relatively unaffected by the large machine size unlike neutral beams. However, fast-wave interactions with the plasma edge can lead to deleterious effects such as, in the case of the high-harmonic fast-wave (HHFW) system on NSTX, large losses of fast-wave power in the scrape off layer (SOL) under certain conditions. In such scenarios, a large fraction of the lost HHFW power is deposited on the upper and lower divertors in bright spiral shapes. The responsible mechanism(s) has not yet been identified but may include fast-wave propagation in the scrape off layer, parametric decay instability, and RF currents driven by the antenna reactive fields. Understanding and mitigating these losses is important not only for improving the heating and current-drive on NSTX-Upgrade but also for understanding fast-wave propagation across the SOL in any fast-wave system. This talk summarizes experimental results demonstrating that the flow of lost HHFW power to the divertor regions largely follows the open SOL magnetic field lines. This lost power flux is relatively large close to both the antenna and the last closed flux surface with a reduced level in between, so the loss mechanism cannot be localized to the antenna. At the same time, significant losses also occur along field lines connected to the inboard edge of the bottom antenna plate. The power lost within the spirals is roughly estimated, showing that these field-aligned losses to the divertor are significant but may not account for the total HHFW loss. To elucidate the role of the onset layer for perpendicular fast-wave propagation with regards to fast-wave propagation in the SOL, a cylindrical cold-plasma model is being developed. This model, in addition to advanced RF codes such as TORIC and AORSA, is aimed at identifying the underlying mechanism(s) behind these SOL losses, to minimize their effects in NSTX-U, and to predict their importance in ITER.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ion cyclotron resonance heating systems upgrade toward high power and CW operations in WEST

The design of the WEST (Tungsten-W Environment in Steady-state Tokamak) Ion cyclotron resonance heating antennas is based on a previously tested conjugate-T Resonant Double Loops prototype equipped with internal vacuum matching capacitors. The design and construction of three new WEST ICRH antennas are being carried out in close collaboration with ASIPP, within the framework of the Associated Laboratory in the fusion field between IRFM and ASIPP. The coupling performance to the plasma and the load-tolerance have been improved, while adding Continuous Wave operation capability by introducing water cooling in the entire antenna. On the generator side, the operation class of the high power tetrodes is changed from AB to B in order to allow high power operation (up to 3 MW per antenna) under higher VSWR (up to 2:1). Reliability of the generators is also improved by increasing the cavity breakdown voltage. The control and data acquisition system is also upgraded in order to resolve and react on fast events, such as ELMs. A new optical arc detection system comes in reinforcement of the V{sub r}/V{sub f} and SHAD systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗