Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TOKAMAK 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 235 records · Page 13

Effect of electron flow on the ordinary-extraordinary mode conversion

Ordinary-extraordinary mode conversion in the electron cyclotron frequency range is revisited in the presence of a flowing electron component. The analytical expressions of optimal parallel refraction index and conversion efficiency are obtained from a one-dimensional cold plasma model. The presence of flowing electrons leads to an outward shift of the conversion layer and therefore increases the optimal value of parallel refraction index. If this effect is not considered, the efficiency of mode conversion degenerates. In typical tokamak plasmas, this degeneration is about a few percentages, which may induce the reflection of several tens of kilowatts of power from the cutoff layer when injecting megawatts of ECRF power into fusion plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interaction of large scale flow structures with gyrokinetic turbulence

Shear flows have a profound influence on turbulence-driven transport in tokamaks. The introduction of arbitrary initial flow profiles into the code ORB5 [Jolliet et al., Comput. Phys. Commun. 177, 409 (2007)] allows the convenient study of how flows on all length scales both influence transport levels and self-consistently evolve. A formulation is presented which preserves the canonical structure of the background particle distribution when either toroidal or poloidal flows are introduced. Turbulence suppression is possible above a certain shearing rate magnitude for homogeneous shear flows, and little evolution of the shearing rate is seen. However, when a flow with a zone boundary, where the shearing rate reverses at mid-radius, is introduced, the shear flow evolves substantially during the simulation. E x B shear flows with a zone boundary of a positive sign decay to a saturation amplitude, consistent with the well known saturation of turbulently generated zonal flows. Unlike the E x B flow, the parallel flows relax diffusively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic description of rotating Tokamak plasmas with anisotropic temperatures in the collisionless regime

A largely unsolved theoretical issue in controlled fusion research is the consistent kinetic treatment of slowly-time varying plasma states occurring in collisionless and magnetized axisymmetric plasmas. The phenomenology may include finite pressure anisotropies as well as strong toroidal and poloidal differential rotation, characteristic of Tokamak plasmas. Despite the fact that physical phenomena occurring in fusion plasmas depend fundamentally on the microscopic particle phase-space dynamics, their consistent kinetic treatment remains still essentially unchallenged to date. The goal of this paper is to address the problem within the framework of Vlasov-Maxwell description. The gyrokinetic treatment of charged particles dynamics is adopted for the construction of asymptotic solutions for the quasi-stationary species kinetic distribution functions. These are expressed in terms of the particle exact and adiabatic invariants. The theory relies on a perturbative approach, which permits to construct asymptotic analytical solutions of the Vlasov-Maxwell system. In this way, both diamagnetic and energy corrections are included consistently into the theory. In particular, by imposing suitable kinetic constraints, the existence of generalized bi-Maxwellian asymptotic kinetic equilibria is pointed out. The theory applies for toroidal rotation velocity of the order of the ion thermal speed. These solutions satisfy identically also the constraints imposed by the Maxwell equations, i.e., quasi-neutrality and Ampere's law. As a result, it is shown that, in the presence of nonuniform fluid and EM fields, these kinetic equilibria can sustain simultaneously toroidal differential rotation, quasi-stationary finite poloidal flows and temperature anisotropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First neutral beam injection experiments on KSTAR tokamak

The first neutral beam (NB) injection system of the Korea Superconducting Tokamak Advanced Research (KSTAR) tokamak was partially completed in 2010 with only 1/3 of its full design capability, and NB heating experiments were carried out during the 2010 KSTAR operation campaign. The ion source is composed of a JAEA bucket plasma generator and a KAERI large multi-aperture accelerator assembly, which is designed to deliver a 1.5 MW, NB power of deuterium at 95 keV. Before the beam injection experiments, discharge, and beam extraction characteristics of the ion source were investigated. The ion source has good beam optics in a broad range of beam perveance. The optimum perveance is 1.1-1.3 {mu}P, and the minimum beam divergence angle measured by the Doppler shift spectroscopy is 0.8 deg. The ion species ratio is D{sup +}:D{sub 2}{sup +}:D{sub 3}{sup +}= 75:20:5 at beam current density of 85 mA/cm{sup 2}. The arc efficiency is more than 1.0 A/kW. In the 2010 KSTAR campaign, a deuterium NB power of 0.7-1.5 MW was successfully injected into the KSTAR plasma with a beam energy of 70-90 keV. L-H transitions were observed within a wide range of beam powers relative to a threshold value. The edge pedestal formation in the T{sub i} and T{sub e} profiles was verified through CES and electron cyclotron emission diagnostics. In every deuterium NB injection, a burst of D-D neutrons was recorded, and increases in the ion temperature and plasma stored energy were found.

43 PARTICLE ACCELERATORS↗

Radial transport of energetic ions in the presence of trapped electron mode turbulence

The nature of transport of hot ions is studied in the presence of microturbulence generated by the trapped electron mode in a Tokamak using massively parallel, first principle based global nonlinear gyrokinetic simulation, and with the help of a passive tracer method. Passing and trapped hot ions are observed to exhibit inverse and inverse square scaling with energy, while those with isotropic pitch distribution are found to exhibit inverse dependence on energy. For all types of hot ions, namely, isotropic, passing, and trapped, the radial transport appears to be subdiffusive for the parameters considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Convective motion in collisionless trapped electron mode turbulence

Global gyrokinetic particle simulation of collisionless trapped electron mode turbulence in toroidal plasmas finds both diffusive and convective electron motion using a Lagrangian analysis. The convective motion is identified using simulation and analytic theory to arise from the conservation of the second invariant when resonant trapped electrons lose kinetic energy to the drift wave. A resonance broadening model fits well the diffusive and convective electron motion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic energy principle and neoclassical toroidal torque in tokamaks

It is shown that when tokamaks are perturbed, the kinetic energy principle is closely related to the neoclassical toroidal torque by the action invariance of particles. Especially when tokamaks are perturbed from scalar pressure equilibria, the imaginary part of the potential energy in the kinetic energy principle is equivalent to the toroidal torque by the neoclassical toroidal viscosity. A unified description therefore should be made for both physics. It is also shown in this case that the potential energy operator can be self-adjoint and thus the stability calculation can be simplified by minimizing the potential energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct measurement of energetic electron flow in Q-shu University experiment with steady-state spherical tokamak

In magnetized plasmas, the presence of a significant number of energetic electrons has been observed but quantitative characteristics of these electrons are proving difficult to investigate. A Langmuir probe offers a means to provide quantitative measurement of these energetic electrons that takes into account electron emissions (secondary electron emission and electron reflection) from the probe tips and sheath expansion around the probe tips caused by a considerable negative potential. In this paper, these effects are experimentally confirmed and an analytical means to measure energetic electron characteristics are proposed. An analysis of plasmas produced by a high frequency wave is then applied leading to the successful detection of an asymmetric flow of energetic electrons. The estimated electron temperature and current density were approximately 4-5 keV and 2-3 kA/m{sup 2}.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Topology of tokamak plasma equilibria with toroidal current reversal

Some general principles about scalar functions with critical points are used to rigorously ascertain that magnetic equilibria with both toroidal current reversal and nested magnetic surfaces are atypical solutions and highly unstable to arbitrary perturbations of boundary conditions and other parameters. The cause for such is shown to lie in the condition of nested magnetic surfaces and not in the possibility of current reversal and consequent vanishing of the poloidal field inside the plasma. Rather than supporting the claim that instability against experimentally driven perturbations forbids configurations with toroidal current reversal, it is argued that these can be attained if an axisymmetric island system is allowed for in order to break the condition of nested magnetic surfaces. A number of results previously reported in the literature are discussed and reinterpreted under the proposed framework, providing some physical insight on the nature of equilibria with toroidal current reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytic model for coaxial helicity injection in tokamak plasmas

Using a partial differential equation for the time evolution of the mean-field poloidal magnetic flux that incorporates resistivity {eta} and hyper-resistivity {Lambda} terms, an exact analytic solution is obtained for steady-state coaxial helicity injection (CHI) in force-free large aspect ratio tokamaks. The analytic mean-field Ohm's law model allows for calculation of the tokamak CHI current drive efficiency and the plasma inductances at arbitrary levels of magnetic fluctuations, or dynamo activity. The results of the mean-field model suggest that CHI approaching Ohmic efficiency is only possible in tokamaks when the size of the effective current drive boundary layer, {delta}{identical_to}({Lambda}/{eta}){sup 1/2}, becomes greater than half the size of the plasma, {delta}>a/2, with a the plasma minor radius. The electron thermal diffusivity due to magnetic fluctuation induced transport is obtained from the expression {chi}{sub e}={Lambda}/{mu}{sub 0}d{sub e}{sup 2}, with {mu}{sub 0} the permeability of free space and d{sub e} the electron skin depth, which for typical tokamak fusion plasma parameters is on the order of a millimeter. Thus, the ratio of the energy confinement time to the resistive diffusion time in a tokamak plasma driven by steady-state CHI approaching Ohmic efficiency is shown to be constrained by the relation {tau}{sub E}/{tau}{sub {eta}}<(d{sub e}/a){sup 2}{approx_equal}10{sup -6}. The mean-field model suggests that steady-state CHI can be viewed most simply as a boundary layer of stochastically wandering magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New ion source for KSTAR neutral beam injection system

The neutral beam injection system (NBI-1) of the KSTAR tokamak can accommodate three ion sources; however, it is currently equipped with only one prototype ion source. In the 2010 and 2011 KSTAR campaigns, this ion source supplied deuterium neutral beam power of 0.7-1.6 MW to the KSTAR plasma with a beam energy of 70-100 keV. A new ion source will be prepared for the 2012 KSTAR campaign with a much advanced performance compared with the previous one. The newly designed ion source has a very large transparency ({approx}56%) without deteriorating the beam optics, which is designed to deliver a 2 MW injection power of deuterium beams at 100 keV. The plasma generator of the ion source is of a horizontally cusped bucket type, and the whole inner wall, except the cathode filaments and plasma grid side, functions as an anode. The accelerator assembly consists of four multi-circular aperture grids made of copper and four electrode flanges made of aluminum alloy. The electrodes are insulated using PEEK. The ion source will be completed and tested in 2011.

43 PARTICLE ACCELERATORS↗

Gyrokinetic determination of the electrostatic potential of rotating magnetic islands in tokamaks

The electrostatic potential related to a magnetic island structure with imposed width and rotation frequency is studied by means of gyrokinetic simulations, which allow its self-consistent determination via the Poisson equation. An adiabatic response of the trapped ions at the island separatrix leads to a significant smoothing of the potential with respect to analytic calculations based on a complete flattening of the pressure profile inside the island. As a consequence, the magnitude of the polarization current is drastically reduced. When the island size is comparable to the ion banana width, the adiabatic response covers the whole island region, leading to a reduced density flattening for islands rotating in the electron diamagnetic direction. This confirms previous results based on drift-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak L/H mode transition

Through the non field-aligned rotational tokamak equilibrium of a divergence-free plasma flow with a pair of transformed plasma variables w-vector{sub *}=({mu}{rho}){sup 1/2}{nu}-vector and {mu}p{sub *}=({mu}p+w{sub *}{sup 2}/2)[K. H. Tsui, Phys. Plasmas 18, 072502 (2011)], a preliminary understanding of the L/H equilibrium transition is proposed through a feedback cycle, where the higher plasma flux due to external drives enters the rotational Grad-Shafranov equation through the velocity dependent poloidal plasma {beta} to generate the H equilibrium. This H rotational mode has the characteristics of higher normal electric field and plasma pressure. Coupled to the transport properties of E-vector x B-vector drift transport barrier leading to a higher plasma pressure, this makes the H mode a self-sustained equilibrium. The higher plasma {beta} then feeds back to the equilibrium and completes the feedback loop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of non-uniform electron energy distribution function on plasma production in large arc driven negative ion source

Spatially non-uniform electron energy distribution function (EEDF) in an arc driven negative ion source (JAEA 10A negative ion source: 10 A NIS) is calculated numerically by a three-dimensional Monte Carlo kinetic model for electrons to understand spatial distribution of plasma production (such as atomic and ionic hydrogen (H{sup 0}/H{sup +}) production) in source chamber. The local EEDFs were directly calculated from electron orbits including electromagnetic effects and elastic/inelastic collision forces. From the EEDF, spatial distributions of H{sup 0}/H{sup +} production rate were obtained. The results suggest that spatial non-uniformity of H{sup 0}/H{sup +} productions is enhanced by high energy component of EEDF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Symmetries of resistive two-fluid magnetohydrodynamics under reversals of toroidal field, current, and rotation

The discrete symmetries of time-dependent, non-axisymmetric, two-fluid magnetohydrodynamics (MHD) equations are considered. Solutions of the time-independent, axisymmetric, ideal-MHD equations remain solutions under reversals of the toroidal field, current density, or rotation. Introducing non-axisymmetry, resistivity, and two-fluid effects into the equations each break different symmetries. Symmetry groups for solutions possessing up-down spatial symmetry, and the effect of flipping the magnetic geometry across the horizontal midplane, are also considered. It is shown that poloidal velocity may or may not be reversed under a simultaneous reversal of the toroidal field and a vertical reflection of the magnetic geometry, depending on the symmetry of the velocity. Because the symmetry groups of ideal, resistive, and two-fluid MHD are distinct, it may be possible to ascertain the dominant physical mechanism of various phenomena through empirical observations of symmetry-breaking. These results, which hold for nonlinear solutions in arbitrary geometry, should also be of use in testing numerical codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical analysis of radio-frequency sheath-plasma interactions in the ion cyclotron range of frequencies

A new finite element numerical scheme for analyzing self-consistent radio-frequency (RF) sheath-plasma interaction problems in the ion cyclotron range of frequencies is applied to various problems represented by simplified models for the tokamak scrape-off layer. The present code incorporates a modified boundary condition, which is called a sheath boundary condition, that couples the radio-frequency waves and sheaths at the material boundaries by treating the sheath as a thin vacuum layer. A series of numerical analyses in one- and two-dimensional domains show several important physical properties, such as the existence of multiple roots, hysteresis effects, presence and characteristics of the sheath-plasma waves, and the phase shift of a reflected slow wave, some of which are newly identified by introducing a spatially varying plasma density and background magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Omnigenity as generalized quasisymmetry

Any viable stellarator reactor will need to be nearly omnigenous, meaning the radial guiding-center drift velocity averages to zero over time for all particles. While omnigenity is easier to achieve than quasisymmetry, we show here that several properties of quasisymmetric plasmas also apply directly or with only minor modification to the larger class of omnigenous plasmas. For example, concise expressions exist for the flow and current, closely resembling those for a tokamak, and these expressions are explicit in that no magnetic differential equations remain. A helicity (M, N) can be defined for any omnigenous field, based on the topology by which B contours close on a flux surface, generalizing the helicity associated with quasisymmetric fields. For generalized quasi-poloidal symmetry (M = 0), the bootstrap current vanishes, which may yield desirable equilibrium and stability properties. A concise expression is derived for the radial electric field in any omnigenous plasma that is not quasisymmetric. The fact that the tokamak-like analytical calculations are possible in omnigenous plasmas despite their fully 3D magnetic spectrum makes these configurations useful for gaining insight and benchmarking codes. A construction is given to produce omnigenous B({theta},{zeta}) patterns with stellarator symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculations of two-fluid linear response to non-axisymmetric fields in tokamaks

The zero-frequency linear plasma response to static applied non-axisymmetric fields is calculated using a resistive two-fluid model in diverted, toroidal geometry. Within this model, the effects on the plasma response of resistivity, rotation, differential ion and electron velocity, and dissipation are explored. Rotation is generally found to inhibit the formation of islands in the plasma, in qualitative agreement with theoretical results. When two-fluid effects are included, it is found that the penetration of the non-axisymmetric fields is generally greatest when the part of the electron rotation perpendicular to the equilibrium magnetic field is small at the mode-rational surface. Strong rotation shear in the edge is found to enhance the plasma response there. The entire plasma, including the separatrix and scrape-off layer, is included in the computational domain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗