Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TOKAMAK”

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 37 records · Page 2

A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma

In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – driven by ionization, charge exchange and wall loss – is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A note on the application of the Prigogine theorem to rotation of tokamak-plasmas in absence of external torques

Rotation of tokamak-plasmas, not at the mechanical equilibrium, is investigated using the Prigogine thermodynamic theorem. This theorem establishes that, for systems confined in rectangular boxes, the global motion of the system with barycentric velocity does not contribute to dissipation. This result, suitably applied to toroidally confined plasmas, suggests that the global barycentric rotations of the plasma, in the toroidal and poloidal directions, are pure reversible processes. In case of negligible viscosity and by supposing the validity of the balance equation for the internal forces, we show that the plasma, even not in the mechanical equilibrium, may freely rotate in the toroidal direction with an angular frequency, which may be higher than the neoclassical estimation. In addition, its toroidal rotation may cause the plasma to rotate globally in the poloidal direction at a speed faster than the expression found by the neoclassical theory. The eventual configuration is attained when the toroidal and poloidal angular frequencies reaches the values that minimize dissipation. The physical interpretation able to explain the reason why some layers of plasma may freely rotate in one direction while, at the same time, others may freely rotate in the opposite direction, is also provided. Invariance properties, herein studied, suggest that the dynamic phase equation might be of the second order in time. We then conclude that a deep and exhaustive study of the invariance properties of the dynamical and thermodynamic equations is the most correct and appropriate way for understanding the triggering mechanism leading to intrinsic plasma-rotation in toroidal magnetic configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anomalous fast ion losses at high β on the tokamak fusion test reactor

This paper describes experiments carried out on the Tokamak Fusion Test Reactor (TFTR) [R. J. Hawryluk et al., Plasma Phys. Controlled Fusion 33, 1509 (1991)] to investigate the dependence of β-limiting disruption characteristics on toroidal field strength. The hard disruptions found at the β-limit in high field plasmas were not found at low field, even for β's 50% higher than the empirical β-limit of β{sub n} ≈ 2 at high field. Comparisons of experimentally measured β's to TRANSP simulations suggest anomalous loss of up to half of the beam fast ions in the highest β, low field shots. The anomalous transport responsible for the fast ion losses may at the same time broaden the pressure profile. Toroidal Alfvén eigenmodes, fishbone instabilities, and Geodesic Acoustic Modes are investigated as possible causes of the enhanced losses. Here, we present the first observations of high frequency fishbones [F. Zonca et al., Nucl. Fusion 49, 085009 (2009)] on TFTR. The interpretation of Axi-symmetric Beam-driven Modes as Geodesic Acoustic Modes and their possible correlation with transport barrier formation are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Intrinsic torque reversals induced by magnetic shear effects on the turbulence spectrum in tokamak plasmas

Intrinsic torque, which can be generated by turbulent stresses, can induce toroidal rotation in a tokamak plasma at rest without direct momentum injection. Reversals in intrinsic torque have been inferred from the observation of toroidal velocity changes in recent lower hybrid current drive (LHCD) experiments. This work focuses on understanding the cause of LHCD-induced intrinsic torque reversal using gyrokinetic simulations and theoretical analyses. A new mechanism for the intrinsic torque reversal linked to magnetic shear (s{sup ^}) effects on the turbulence spectrum is identified. This reversal is a consequence of the ballooning structure at weak s{sup ^}. Based on realistic profiles from the Alcator C-Mod LHCD experiments, simulations demonstrate that the intrinsic torque reverses for weak s{sup ^} discharges and that the value of s{sup ^}{sub crit} is consistent with the experimental results s{sup ^}{sub crit}{sup exp}≈0.2∼0.3 [Rice et al., Phys. Rev. Lett. 111, 125003 (2013)]. The consideration of this intrinsic torque feature in our work is important for the understanding of rotation profile generation at weak s{sup ^} and its consequent impact on macro-instability stabilization and micro-turbulence reduction, which is crucial for ITER. It is also relevant to internal transport barrier formation at negative or weakly positive s{sup ^}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuum kinetic modeling of the tokamak plasma edge

The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The impact of collisionality, FLR, and parallel closure effects on instabilities in the tokamak pedestal: Numerical studies with the NIMROD code

The extended-MHD NIMROD code [C. R. Sovinec and J. R. King, J. Comput. Phys. 229, 5803 (2010)] is verified against the ideal-MHD ELITE code [H. R. Wilson et al., Phys. Plasmas 9, 1277 (2002)] on a diverted tokamak discharge. When the NIMROD model complexity is increased incrementally, resistive and first-order finite-Larmour radius effects are destabilizing and stabilizing, respectively. Lastly, the full result is compared to local analytic calculations which are found to overpredict both the resistive destabilization and drift stabilization in comparison to the NIMROD computations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of frequency modulation reflectometer for Korea Superconducting Tokamak Advanced Research tokamak

Frequency modulation reflectometer has been developed to measure the plasma density profile of the Korea Superconducting Tokamak Advanced Research tokamak. Three reflectometers are operating in extraordinary polarization mode in the frequency range of Q band (33.6–54 GHz), V band (48–72 GHz), and W band (72–108 GHz) to measure the density up to 7 × 10{sup 19} m{sup −3} when the toroidal magnetic field is 2 T on axis. The antenna is installed inside of the vacuum vessel. A new vacuum window is developed by using 50 μm thick mica film and 0.1 mm thick gold gasket. The filter bank of low pass filter, notch filter, and Faraday isolator is used to reject the electron cyclotron heating high power at attenuation of 60 dB. The full frequency band is swept in 20 μs. The mixer output is directly digitized with sampling rate of 100 MSamples/s. The phase is obtained by using wavelet transform. The whole hardware and software system is described in detail and the measured density profile is presented as a result.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effects of heating power on divertor in-out asymmetry and scrape-off layer flow in reversed field on Experimental Advanced Superconducting Tokamak

The dependence of divertor asymmetry and scrape-off layer (SOL) flow on heating power has been investigated in the Experimental Advanced Superconducting Tokamak (EAST). Divertor plasma exhibits an outboard-enhanced in-out asymmetry in heat flux in lower single null configuration for in reversed (ion ∇B drift direction toward the upper X-point) field directions. Upper single null exhibits an inboard-favored asymmetry in low heating power condition, while exhibits an outboard-favored asymmetry when increasing the heating power. Double null has the strongest in-out asymmetry in heat flux, favoring the outer divertor. The in-out asymmetry ratios of q{sub t,out}/q{sub t,in} and P{sub out}/P{sub total} increase with the power across the separatrix P{sub loss}, which is probably induced by the enhanced radial particle transport due to a large pressure gradient. The characteristics of the measured SOL parallel flow under various discharge conditions are consistent with the Pfirsch-Schlüter (PS) flow with the parallel Mach number M{sub ∥} decreasing with the line averaged density but increasing with P{sub loss}, in the same direction as the PS flow. The contributions of both poloidal E×B drift and parallel flow on poloidal particle transport in SOL on EAST are also assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparisons between tokamak fueling of gas puffing and supersonic molecular beam injection in 2D simulations

Plasma fueling with high efficiency and deep injection is very important to enable fusion power performance requirements. It is a powerful and efficient way to study neutral transport dynamics and find methods of improving the fueling performance by doing large scale simulations. Two basic fueling methods, gas puffing (GP) and supersonic molecular beam injection (SMBI), are simulated and compared in realistic divertor geometry of the HL-2A tokamak with a newly developed module, named trans-neut, within the framework of BOUT++ boundary plasma turbulence code [Z. H. Wang et al., Nucl. Fusion 54, 043019 (2014)]. The physical model includes plasma density, heat and momentum transport equations along with neutral density, and momentum transport equations. Transport dynamics and profile evolutions of both plasma and neutrals are simulated and compared between GP and SMBI in both poloidal and radial directions, which are quite different from one and the other. It finds that the neutrals can penetrate about four centimeters inside the last closed (magnetic) flux surface during SMBI, while they are all deposited outside of the LCF during GP. It is the radial convection and larger inflowing flux which lead to the deeper penetration depth of SMBI and higher fueling efficiency compared to GP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Active and passive kink mode studies in a tokamak with a movable ferromagnetic wall

High-resolution active and passive kink mode studies are conducted in a tokamak with an adjustable ferromagnetic wall near the plasma surface. Ferritic tiles made from 5.6 mm thick Hiperco{sup ®} 50 alloy have been mounted on the plasma-facing side of half of the in-vessel movable wall segments in the High Beta Tokamak-Extended Pulse device [D. A. Maurer et al., Plasma Phys. Controlled Fusion 53, 074016 (2011)] in order to explore ferritic resistive wall mode stability. Low-activation ferritic steels are a candidate for structural components of a fusion reactor, and these experiments examine MHD stability of plasmas with nearby ferromagnetic material. Plasma-wall separation for alternating ferritic and non-ferritic wall segments is adjusted between discharges without opening the vacuum vessel. Amplification of applied resonant magnetic perturbations and plasma disruptivity are observed to increase when the ferromagnetic wall is close to plasma surface instead of the standard stainless steel wall. Rapidly rotating m/n=3/1 external kink modes have higher growth rates with the nearby ferritic wall. Feedback suppression of kinks is still as effective as before the installation of ferritic material in vessel, in spite of increased mode growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The tokamak density limit: A thermo-resistive disruption mechanism

The behavior of magnetic islands with 3D electron temperature and the corresponding 3D resistivity effects on growth are examined for islands with near-zero net heating in the island interior. We refer to the resulting class of non-linearities as thermo-resistive effects. In particular, the effects of varying impurity mix on the previously proposed local island onset threshold [Gates and Delgado-Aparicio, Phys. Rev. Lett. 108, 165004 (2012)] are examined and shown to be consistent with the well established experimental scalings for tokamaks at the density limit. A surprisingly simple semi-analytic theory is developed which imposes the effects of heating/cooling in the island interior as well as the effects of island geometry. For the class of current profiles considered, it is found that a new term that accounts for the thermal effects of island asymmetry is required in the modified Rutherford equation. The resultant model is shown to exhibit a robust onset of a rapidly growing tearing mode—consistent with the disruption mechanism observed at the density limit in tokamaks. A fully non-linear 3D cylindrical calculation is performed that simulates the effect of net island heating/cooling by raising/suppressing the temperature in the core of the island. In both the analytic theory and the numerical simulation, the sudden threshold for rapid growth is found to be due to an interaction between three distinct thermal non-linearities which affect the island resistivity, thereby modifying the growth dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of three-dimensional fields on bounce averaged particle drifts in a tokamak

The impact of applied 3D magnetic fields on the bounce-averaged precessional drifts in a tokamak plasma are calculated. Local 3D MHD equilibrium theory is used to construct solutions to the equilibrium equations in the vicinity of a magnetic surface for a large aspect ratio circular tokamak perturbed by applied 3D fields. Due to modulations of the local shear caused by near-resonant Pfirsch-Schlüter currents, relatively weak applied 3D fields can have a large effect on trapped particle precessional drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and installation of a ferromagnetic wall in tokamak geometry

Low-activation ferritic steels are leading material candidates for use in next-generation fusion development experiments such as a prospective component test facility and DEMO power reactor. Understanding the interaction of plasmas with a ferromagnetic wall will provide crucial physics for these facilities. In order to study ferromagnetic effects in toroidal geometry, a ferritic wall upgrade was designed and installed in the High Beta Tokamak–Extended Pulse (HBT-EP). Several material options were investigated based on conductivity, magnetic permeability, vacuum compatibility, and other criteria, and the material of choice (high-cobalt steel) is characterized. Installation was accomplished quickly, with minimal impact on existing diagnostics and overall machine performance, and initial results demonstrate the effects of the ferritic wall on plasma stability.

42 ENGINEERING↗

Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions

As a follow-up on the drift-kinetic study of the non-local bootstrap current in the steep edge pedestal of tokamak plasma by Koh et al. [Phys. Plasmas 19, 072505 (2012)], a gyrokinetic neoclassical study is performed with gyrokinetic ions and drift-kinetic electrons. Besides the gyrokinetic improvement of ion physics from the drift-kinetic treatment, a fully non-linear Fokker-Planck collision operator—that conserves mass, momentum, and energy—is used instead of Koh et al.'s linearized collision operator in consideration of the possibility that the ion distribution function is non-Maxwellian in the steep pedestal. An inaccuracy in Koh et al.'s result is found in the steep edge pedestal that originated from a small error in the collisional momentum conservation. The present study concludes that (1) the bootstrap current in the steep edge pedestal is generally smaller than what has been predicted from the small banana-width (local) approximation [e.g., Sauter et al., Phys. Plasmas 6, 2834 (1999) and Belli et al., Plasma Phys. Controlled Fusion 50, 095010 (2008)], (2) the plasma flow evaluated from the local approximation can significantly deviate from the non-local results, and (3) the bootstrap current in the edge pedestal, where the passing particle region is small, can be dominantly carried by the trapped particles in a broad trapped boundary layer. A new analytic formula based on numerous gyrokinetic simulations using various magnetic equilibria and plasma profiles with self-consistent Grad-Shafranov solutions is constructed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rotation profile flattening and toroidal flow shear reversal due to the coupling of magnetic islands in tokamaks

The electromagnetic coupling of helical modes, even those having different toroidal mode numbers, modifies the distribution of toroidal angular momentum in tokamak discharges. This can have deleterious effects on other transport channels as well as on magnetohydrodynamic (MHD) stability and disruptivity. At low levels of externally injected momentum, the coupling of core-localized modes initiates a chain of events, whereby flattening of the core rotation profile inside successive rational surfaces leads to the onset of a large m/n = 2/1 tearing mode and locked-mode disruption. With increased torque from neutral beam injection, neoclassical tearing modes in the core may phase-lock to each other without locking to external fields or structures that are stationary in the laboratory frame. The dynamic processes observed in these cases are in general agreement with theory, and detailed diagnosis allows for momentum transport analysis to be performed, revealing a significant torque density that peaks near the 2/1 rational surface. However, as the coupled rational surfaces are brought closer together by reducing q{sub 95}, additional momentum transport in excess of that required to attain a phase-locked state is sometimes observed. Rather than maintaining zero differential rotation (as is predicted to be dynamically stable by single-fluid, resistive MHD theory), these discharges develop hollow toroidal plasma fluid rotation profiles with reversed plasma flow shear in the region between the m/n = 3/2 and 2/1 islands. The additional forces expressed in this state are not readily accounted for, and therefore, analysis of these data highlights the impact of mode coupling on torque balance and the challenges associated with predicting the rotation dynamics of a fusion reactor—a key issue for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗