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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

The mixed Wentzel-Kramers-Brillouin-full-wave approach and its application to lower hybrid wave propagation and absorption

The mixed Wentzel-Kramers-Brillouin (WKB)-full-wave approach for the calculation of the 2D mode structure in tokamak plasmas is further developed based on our previous work [A. Cardinali and F. Zonca, Phys. Plasmas 10, 4199 (2003) and Z. X. Lu et al., Phys. Plasmas 19, 042104 (2012)]. A new scheme for theoretical analysis and numerical implementation of the mixed WKB-full-wave approach is formulated, based on scale separation and asymptotic analysis. Besides its capability to efficiently investigate the initial value problem for 2D mode structures and linear stability, in this work, the mixed WKB-full-wave approach is extended to the investigation of radio frequency wave propagation and absorption, e.g., lower hybrid waves. As a novel method, its comparison with other approaches, e.g., WKB and beam tracing methods, is discussed. Its application to lower hybrid wave propagation in concentric circular tokamak plasmas using typical FTU discharge parameters is also demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Jump conditions in transonic equilibria

In the present paper, the numerical calculation of transonic equilibria, first introduced with the FLOW code in Guazzotto et al.[Phys. Plasmas 11, 604 (2004)], is critically reviewed. In particular, the necessity and effect of imposing explicit jump conditions at the transonic discontinuity are investigated. It is found that 'standard' (low-{beta}, large aspect ratio) transonic equilibria satisfy the correct jump condition with very good approximation even if the jump condition is not explicitly imposed. On the other hand, it is also found that high-{beta}, low aspect ratio equilibria require the correct jump condition to be explicitly imposed. Various numerical approaches are described to modify FLOW to include the jump condition. It is proved that the new methods converge to the correct solution even in extreme cases of very large {beta}, while they agree with the results obtained with the old implementation of FLOW in lower-{beta} equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Toroidal modeling of penetration of the resonant magnetic perturbation field

A toroidal, quasi-linear model is proposed to study the penetration dynamics of the resonant magnetic perturbation (RMP) field into the plasma. The model couples the linear, fluid plasma response to a toroidal momentum balance equation, which includes torques induced by both fluid electromagnetic force and by (kinetic) neoclassical toroidal viscous (NTV) force. The numerical results for a test toroidal equilibrium quantify the effects of various physical parameters on the field penetration and on the plasma rotation braking. The neoclassical toroidal viscous torque plays a dominant role in certain region of the plasma, for the RMP penetration problem considered in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Toroidal momentum transport in a tokamak caused by symmetry breaking parallel derivatives

A new mechanism for toroidal momentum transport in a tokamak is investigated using the gyro-kinetic model. First, an analytic model is developed through the use of the ballooning transform. The terms that generate the momentum transport are then connected with the poloidal derivative of the ballooning envelope, which are one order smaller in the normalised Larmor radius, compared with the derivative of the eikonal. The mechanism, therefore, does not introduce an inhomogeneity in the radial direction, in contrast with the effect of profile shearing. Numerical simulations of the linear ion temperature gradient mode with adiabatic electrons, retaining the finite {rho}{sub *} effects in the E Multiplication-Sign B velocity, the drift, and the gyro-average, are presented. The momentum flux is found to be linear in the normalised Larmor radius ({rho}{sub *}) but is, nevertheless, generating a sizeable counter-current rotation. The total momentum flux scales linear with the aspect ratio of the considered magnetic surface, and increases with increasing magnetic shear, safety factor, and density and temperature gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Nonlinear dynamics of multiple neoclassical tearing modes in tokamaks

Numerical simulations investigating the interaction of co-existent 2/1 and 3/1 neoclassical tearing modes are presented. The results obtained from an initial value 3D toroidal code that solves a set of generalized reduced MHD equations exhibit a host of complex phenomena arising from the coupling of the two modes. These include a modification of the island saturation widths of the two modes, a significant modification in the perpendicular flow patterns in the vicinity of the islands, and the excitation of geodesic acoustic mode like oscillations that lead to concomitant oscillations in the kinetic and magnetic energies of the islands. These oscillations only occur in the presence of the neoclassical stress tensor contribution and are absent for two coupled classical tearing modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Robustness of predator-prey models for confinement regime transitions in fusion plasmas

Energy transport and confinement in tokamak fusion plasmas is usually determined by the coupled nonlinear interactions of small-scale drift turbulence and larger scale coherent nonlinear structures, such as zonal flows, together with free energy sources such as temperature gradients. Zero-dimensional models, designed to embody plausible physical narratives for these interactions, can help to identify the origin of enhanced energy confinement and of transitions between confinement regimes. A prime zero-dimensional paradigm is predator-prey or Lotka-Volterra. Here, we extend a successful three-variable (temperature gradient; microturbulence level; one class of coherent structure) model in this genre [M. A. Malkov and P. H. Diamond, Phys. Plasmas 16, 012504 (2009)], by adding a fourth variable representing a second class of coherent structure. This requires a fourth coupled nonlinear ordinary differential equation. We investigate the degree of invariance of the phenomenology generated by the model of Malkov and Diamond, given this additional physics. We study and compare the long-time behaviour of the three-equation and four-equation systems, their evolution towards the final state, and their attractive fixed points and limit cycles. We explore the sensitivity of paths to attractors. It is found that, for example, an attractive fixed point of the three-equation system can become a limit cycle of the four-equation system. Addressing these questions which we together refer to as 'robustness' for convenience is particularly important for models which, as here, generate sharp transitions in the values of system variables which may replicate some key features of confinement transitions. Our results help to establish the robustness of the zero-dimensional model approach to capturing observed confinement phenomenology in tokamak fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

On the toroidal current density flowing across a poloidal-magnetic-field null in an axisymmetric plasma

The axisymmetry condition and two of Maxwell's equations are used to show that, in general, there are no nested magnetic surfaces around a poloidal-magnetic-field null for a sufficiently small value of the toroidal current density flowing there. Hence, the toroidal current density at the axis of a magnetic configuration with extreme shear reversal cannot continuously approach zero unless nested surfaces are first broken or particular values are assigned to boundary conditions and other plasma parameters. The threshold of the toroidal current-density at which the topology changes is shown to be set by such parameters, and some examples of the predicted topology transition are presented using analytical solutions of the Grad-Shafranov equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Stability of compressible reduced magnetohydrodynamic equilibria-Analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical {delta}W criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of {delta}W is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Non-linear modulation of short wavelength compressional Alfven eigenmodes

Most Alfvenic activity in the frequency range between toroidal Alfven eigenmodes and roughly one half of the ion cyclotron frequency on National Spherical Torus eXperiment [Ono et al., Nucl. Fusion 40, 557 (2000)], that is, approximately 0.3 MHz up to Almost-Equal-To 1.2 MHz, are modes propagating counter to the neutral beam ions. These have been modeled as Compressional and Global Alfven Eigenmodes (CAE and GAE) and are excited through a Doppler-shifted cyclotron resonance with the beam ions. There is also a class of co-propagating modes at higher frequency than the counter-propagating CAE and GAE. These modes have been identified as CAE, and are seen mostly in the company of a low frequency, n = 1 kink-like mode. In this paper, we present measurements of the spectrum of these high frequency CAE (hfCAE) and their mode structure. We compare those measurements to a simple model of CAE and present a predator-prey type model of the curious non-linear coupling of the hfCAE and the low frequency kink-like mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Role of ion temperature on scrape-off layer plasma turbulence

Turbulence in Scrape-off layer (SOL) of tokamak plasma has been studied numerically using interchange modes with the help of electron continuity, quasineutrality, and ion energy equations. Electron temperature is assumed uniform. We have studied dynamics of seeded plasma blob and plasma turbulence to identify the role of ion temperature and its gradient. The ion temperature elongates the blob poloidally and reduces its radial velocity. Initial dipole nature of the plasma blob potential breaks and generates few more dipoles during its propagation in the SOL. Plasma turbulence simulation shows poloidally elongated density and ion temperature structures that are similar to the seeded blob simulation studies. Fluctuations of the density and ion temperature have been presented as function of scale lengths of the density and ion temperature. Reduction of the SOL width and increase of radial electric field have been measured in the presence of the ion temperature. Particle and energy transports have been also presented as the function of the density and ion temperature scale lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Visible imaging and spectroscopy of disruption runaway electrons in DIII-D

The first visible light images of synchrotron emission from disruption runaway electrons are presented. The forward-detected continuum radiation from runaways is identified as synchrotron emission by comparing two survey spectrometers and two visible fast cameras viewing in opposite toroidal directions. Analysis of the elongation of 2D synchrotron images of oval-shaped runaway beams indicates that the velocity pitch angle v{sub Up-Tack }/v{sub ||} ranges from 0.1 to 0.2 for the detected electrons, with energies above 25 MeV. Analysis of synchrotron intensity from a camera indicates that the tail of the runaway energy distribution reaches energies up to 60 MeV, which agrees with 0D modeling of electron acceleration in the toroidal electric field generated during the current quench. A visible spectrometer provides an independent estimate of the upper limit of runaway electron energy which is roughly consistent with energy determined from camera data. Synchrotron spectra reveal that approximately 1% of the total post-thermal quench plasma current is carried by the detected high-energy runaway population with energies in the range of 25-60 MeV; the bulk of the plasma current thus appears to be carried by relativistic electrons with energy less than 25 MeV. In addition to stable oval shapes, runaway beams with other shapes and internal structure are sometimes observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Particle control and plasma performance in the Lithium Tokamak eXperiment

The Lithium Tokamak eXperiment is a small, low aspect ratio tokamak [Majeski et al., Nucl. Fusion 49, 055014 (2009)], which is fitted with a stainless steel-clad copper liner, conformal to the last closed flux surface. The liner can be heated to 350 °C. Several gas fueling systems, including supersonic gas injection and molecular cluster injection, have been studied and produce fueling efficiencies up to 35%. Discharges are strongly affected by wall conditioning. Discharges without lithium wall coatings are limited to plasma currents of order 10 kA, and discharge durations of order 5 ms. With solid lithium coatings discharge currents exceed 70 kA, and discharge durations exceed 30 ms. Heating the lithium wall coating, however, results in a prompt degradation of the discharge, at the melting point of lithium. These results suggest that the simplest approach to implementing liquid lithium walls in a tokamak—thin, evaporated, liquefied coatings of lithium—does not produce an adequately clean surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Design and implementation of a 150 GHz single-channel millimeter wave interferometer on Joint TEXT tokamak

A simple, single-channel millimeter-wave interferometer system has been designed, fabricated, and installed on the J-TEXT tokamak. For the plasma density anticipated on J-TEXT, a solid-state source operating at 150 GHz has been chosen to minimize errors due to both vibration along the beam path and refraction in the plasma. The new aspect of the interferometer design is to use a subharmonic mixer for detection with a frequency doubled 150 GHz source. It employs a single source which is bias-tuned and modulated with a sawtooth wave form up to 100 kHz in order to generate the intermediate frequency. The 12.5 GHz voltage-controlled oscillator is multiplied to 75 GHz before a final doubler raises it to 150 GHz. A portion of the 75 GHZ power is used for the local oscillator (LO) and is directly connected to the LO input of the subharmonic mixer. The phase is evaluated by a digital phase comparator using a software-based algorithm. Detection noise limits the minimum resolvable phase change with the interferometer to {+-}0.05 fringe, which corresponds to an averaged electron density change along the chord of {+-}1.1 Multiplication-Sign 10{sup 17} m{sup -2}. The maximum measurable electron density is expected to be {approx}9 Multiplication-Sign 10{sup 19} m{sup -3}. A comparison of preliminary results from the millimeter wave interferometer with that from the far-infrared hydrogen cyanide laser (wavelength of 337 {mu}m) interferometer shows good agreement during the pulse flat-top period. The millimeter wave interferometer system will be used as a part of the density feedback control system in the future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Quasilinear model for energetic particle diffusion in radial and velocity space

A quasilinear model for passive energetic particle (EP) turbulent diffusion in radial and velocity space is fitted and tested against nonlinear gyrokinetic tokamak simulations with the GYRO code [J. Candy and R. E. Waltz, Phys. Rev. Lett. 91, 045001 (2003)]. Off diagonal elements of a symmetric positive definite 2 Multiplication-Sign 2 EP diffusion matrix account for fluxes up radial (energy) gradients driven by energy (radial) gradients of the EP velocity space distribution function. The quasilinear ratio kernel of the model is provided by a simple analytic formula for the EP radial and velocity space EP diffusivity relative to radial thermal ion energy diffusivity at each linear mode of the turbulence driven by the thermal plasma. The TGLF [G. M. Staebler, J. E. Kinsey, and R. E. Waltz, Phys. Plasmas 14, 0055909 (2007); ibid. 15, 0055908 (2008)] tokamak transport model provides the linear mode frequency and growth rates to the kernel as well as the nonlinear spectral weight for each mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The wave energy flux of high frequency diffracting beams in complex geometrical optics

We consider the construction of asymptotic solutions of Maxwell's equations for a diffracting wave beam in the high frequency limit and address the description of the wave energy flux transported by the beam. With this aim, the complex eikonal method is applied. That is a generalization of the standard geometrical optics method in which the phase function is assumed to be complex valued, with the non-negative imaginary part accounting for the finite width of the beam cross section. In this framework, we propose an argument which simplifies significantly the analysis of the transport equation for the wave field amplitude and allows us to derive the wave energy flux. The theoretical analysis is illustrated numerically for the case of electron cyclotron beams in tokamak plasmas by using the GRAY code [D. Farina, Fusion Sci. Technol. 52, 154 (2007)], which is based upon the complex eikonal theory. The results are compared to those of the paraxial beam tracing code TORBEAM [E. Poli et al., Comput. Phys. Commun. 136, 90 (2001)], which provides an independent calculation of the energy flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Foreword to Special Issue: Papers from the 54th Annual Meeting of the APS Division of Plasma Physics, Providence, Rhode Island, USA, 2012

Each year, the annual meeting of the APS Division of Plasma Physics (DPP) brings together a broad representation of the many active subfields of plasma physics and enjoys an audience that is equally diverse. The meeting was well attended and largely went as planned despite the interventions of hurricane Sandy which caused the city of Providence to shut-down during the first day of the conference. The meeting began on Monday morning with a review of the physics of cosmic rays, 2012 being the 100th year since their discovery, which illustrated the central importance of plasma physics to astrophysical problems. Subsequent reviews covered the importance of tokamak plasma boundaries, progress towards ignition on the National Ignition Facility (NIF), and magnetized plasma turbulence. The Maxwell prize address, by Professor Liu Chen, covered the field of nonlinear Alfvén wave physics. Tutorial lectures were presented on the verification of gyrokinetics, new capabilities in laboratory astrophysics, magnetic flux compression, and tokamak plasma start-up.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Second harmonic effect on geodesic modes in tokamak plasmas

Results of a kinetic treatment of Geodesic Acoustic Modes (GAMs) that fully takes into account ion parallel dynamics, including the magnetic field component, are presented. The finite-orbit-width (FOW) parameter is considered in the calculation of the second harmonic effect on GAMs. For larger values of the FOW parameter, it is shown that dispersive effects related to the m = 2 harmonics is the cause of the mode frequency splitting and the modes appear due to the interaction with the ion sound mode. Furthermore, the modes may have enhanced damping rates due to second harmonic Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

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↗