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 163 records · Page 9

Ripple modifications to alpha transport in tokamaks

We report magnetic field ripple is inherent in tokamaks since the toroidal magnetic field is generated by a finite number of toroidal field coils. The field ripple results in departures from axisymmetry that cause radial transport losses of particles and heat. These ripple losses are a serious concern for alphas near their birth speed $v_0$ since alpha heating of the background plasma is required to make fusion reactors into economical power plants. Ripple in tokamaks gives rise to at least two alpha transport regimes of concern. As the slowing down time $τ_s$ is much larger than the time for an alpha just born to make a toroidal transit, a regime referred to as the 1/$ν∝τ_s$ regime can be encountered, with ν the appropriate alpha collision frequency. In this regime the radial transport losses increase as $v_0τ_s/R$, with R the major radius of the tokamak. The deleterious effect of ripple transport is mitigated by electric and magnetic drifts within the flux surface. When drift tangent to the flux surface becomes significant another ripple regime, referred to as the $\sqrt{ν}$ regime, is encountered where a collisional boundary layer due to the drift plays a key role. We evaluate the alpha transport in both regimes, taking account of the alphas having a slowing down rather than a Maxwellian distribution function and their being collisionally scattered by a collision operator appropriate for alphas. Alpha ripple transport is found to be in the $\sqrt{ν}$ regime where it will be a serious issue for typical tokamak reactors as it will be well above the axisymmetric neoclassical level and can be large enough to deplete the alpha slowing down distribution function unless toroidal rotation is strong.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A high-density and high-confinement tokamak plasma regime for fusion energy

The tokamak approach, utilizing a toroidal magnetic field configuration to confine a hot plasma, is one of the most promising designs for developing reactors that can exploit nuclear fusion to generate electrical energy. To reach the goal of an economical reactor, most tokamak reactor designs simultaneously require reaching a plasma line-averaged density above an empirical limit—the so-called Greenwald density and attaining an energy confinement quality better than the standard high-confinement mode. However, such an operating regime has never been verified in experiments. In addition, a long-standing challenge in the high-confinement mode has been the compatibility between a high-performance core and avoiding large, transient edge perturbations that can cause very high heat loads on the plasma-facing-components in tokamaks. Here we report the demonstration of stable tokamak plasmas with a line-averaged density approximately 20% above the Greenwald density and an energy confinement quality of approximately 50% better than the standard high-confinement mode, which was realized by taking advantage of the enhanced suppression of turbulent transport granted by high density-gradients in the high-poloidal-beta scenario. Furthermore, our experimental results show an integration of very low edge transient perturbations with the high normalized density and confinement core. The operating regime we report supports some critical requirements in many fusion reactor designs all over the world and opens a potential avenue to an operating point for producing economically attractive fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gamma ray imager on the DIII-D tokamak

A gamma ray camera is built for the DIII-D tokamak [J. Luxon, Nucl. Fusion 42, 614 (2002)] that provides spatial localization and energy resolution of gamma flux by combining a lead pinhole camera with custom-built detectors and optimized viewing geometry. This diagnostic system is installed on the outer midplane of the tokamak such that its 123 collimated sightlines extend across the tokamak radius while also covering most of the vertical extent of the plasma volume. A set of 30 bismuth germanate detectors can be secured in any of the available sightlines, allowing for customizable coverage in experiments with runaway electrons in the energy range of 1–60 MeV. Commissioning of the gamma ray imager includes the quantification of electromagnetic noise sources in the tokamak machine hall and a measurement of the energy spectrum of background gamma radiation. First measurements of gamma rays coming from the plasma provide a suitable testbed for implementing pulse height analysis that provides the energy of detected gamma photons.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Destabilizing effects of edge infernal components on resistive wall modes in advanced tokamak scenarios

The stability of the n = 1 resistive wall modes (RWMs) dominated by the pressure-driven infernal components is investigated using the ideal magnetohydrodynamics (MHD) code AEGIS for the advanced tokamak scenarios. Here, n is the toroidal mode number. In the advanced tokamak scenarios, due to the large fraction of bootstrap current contribution, the profile of safety factor q is deeply reversed in magnetic shear in the central core region and locally flattened within the edge pedestal. Consequently, the pressure-driven infernal components develop in the corresponding flat-q regions of both core and edge. However, the edge infernal components dominate the n = 1 RWM structure and lead to lower β N limits for the advanced tokamak scenarios. In the framework of ideal MHD, the edge rotation is found the most critical to the stabilization due to the dominant influence of the edge infernal components, which should be maintained sufficiently large in magnitude and range in order for the rotation alone to fully suppress the n = 1 RWM in typical advanced tokamak scenarios.

36 MATERIALS SCIENCE↗

Numerical modelling of geodesic acoustic mode relaxation in a tokamak edge

Here, the edge of a tokamak in a high confinement (H mode) regime is characterized by steep density gradients and a large radial electric field. Recent analytical studies demonstrated that the presence of a strong radial electric field consistent with a subsonic pedestal equilibrium modifies the conventional results of the neoclassical formalism developed for the core region. In the present work we make use of the recently developed gyrokinetic code COGENT to numerically investigate neoclassical transport in a tokamak edge including the effects of a strong radial electric field. The results of numerical simulations are found to be in good qualitative agreement with the theoretical predictions and the quantitative discrepancy is discussed. In addition, the present work investigates the effects of a strong radial electric field on the relaxation of geodesic acoustic modes (GAMs) in a tokamak edge. Numerical simulations demonstrate that the presence of a strong radial electric field characteristic of a tokamak pedestal can enhance the GAM decay rate, and heuristic arguments elucidating this finding are provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope effects on L-H threshold and confinement in tokamak plasmas

The dependence of plasma transport and confinement on the main hydrogenic ion isotope mass is of fundamental importance for understanding turbulent transport and, therefore, for accurate extrapolations of confinement from present tokamak experiments, which typically use a single hydrogen isotope, to burning plasmas such as ITER, which will operate in deuterium–tritium mixtures. Knowledge of the dependence of plasma properties and edge transport barrier formation on main ion species is critical in view of the initial, low-activation phase of ITER operations in hydrogen or helium and of its implications on the subsequent operation in deuterium–tritium. The favourable scaling of global energy confinement time with isotope mass, which has been observed in many tokamak experiments, remains largely unexplained theoretically. Moreover, the mass scaling observed in experiments varies depending on the plasma edge conditions. In preparation for upcoming deuterium–tritium experiments in the JET tokamak with the ITER-like Be/W Wall (JET-ILW), a thorough experimental investigation of isotope effects in hydrogen, deuterium and tritium plasmas is being carried out, in order to provide stringent tests of plasma energy, particle and momentum transport models. Recent hydrogen and deuterium isotope experiments in JET-ILW on L-H power threshold, L-mode and H-mode confinement are reviewed and discussed in the context of past and more recent isotope experiments in tokamak plasmas, highlighting common elements as well as contrasting observations that have been reported. The experimental findings are discussed in the context of fundamental aspects of plasma transport models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integral analysis of the effect of material dimension and composition on tokamak neutronics *

Abstract The neutronics performance of a tokamak has been identified as an important factor in designing a fusion power plant. The design of the tokamak should not only meet operational parameters such as sufficient tritium breeding, but also safety parameters such as low structural material activation. This paper investigates the impacts of the neutronics metrics for the ARC-class tokamak, a compact tokamak with an immersion blanket, by perturbing the first five layers of structural material—first wall, inner vacuum vessel, coolant salt channel, neutron multiplier, and outer vacuum vessel. The goal of this work is to provide insight into shaping and scaling the flux on each layer to obtain optimized operational and safety metrics through quantification of the responses from each perturbation. Results show that increased first wall thickness can increase the tritium breeding ratio (TBR) in specific configurations with high 6 Li enrichments and that vacuum vessels decrease TBR for low- 6 Li enrichment configurations. It was also found that the neutron multiplier can either increase or decrease TBR depending on the configuration. The response of metrics to the change in layer thickness and enrichment also varies depending on the vacuum vessel material. The integral impacts of 6 Li enrichment, layer thicknesses, and vacuum vessel material choice are investigated and presented in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analyses of Extreme Ultraviolet Spectra from Tokamak Plasmas

The most recent work was a one-year extension of work started in 2015, for a total duration of three years. During this time we performed detailed spectral analyses of extreme ultraviolet spectra taken by Lawrence Livermore National Laboratory on the National Spherical Torus Experiment and the Lithium Tokamak Experiment at Princeton Plasma Physics Laboratory, and the Alcator C-Mod tokamak at the Massachusetts Institute of Technology. Our analyses focused on the K-shell emission of lithium and the L-shell emission of oxygen, fluorine, carbon, aluminum, silicon, sulfur, and calcium, measured with the Long-Wavelength Extreme Ultraviolet Spectrometer from LLNL. Data analyses are complete, as documented by poster and oral presentations and by published papers. The most recent and final year focused on finishing and publishing the analysis of data from the NSTX tokamak and from the Alcator tokamak, which included documenting charge exchange-produced carbon lines in the EUV spectral region and emission lines from fluorinelike ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Magnetic fusion experiments have the largest gradients in the world that are in steady state. These gradients lead to instabilities, and thus, most of plasma theory over the past four decades have been devoted to the development of instability theory and their use in interpreting plasma results. Instabilities in tokamaks lead to 3 different phenomenological outcomes: coherent saturation, turbulent saturation, and sudden relaxation (which may include disruptions). For long-wavelength instabilities, the NIMROD extended MHD code has emerged as an important tool for understanding tokamak instabilities. Because of the long history of NIMROD, a new version of it has been started at Tech-X to be able to address multi-species capabilities as well as exploit modern GPU systems. The new version builds on the previous version and improves the workflow by enabling new equilibria from experiments to be resolved. The fundamental formulation of the equations underlying NIMROD has also been explored in multiples ways. First, understanding instabilities in tokamaks through numerical methods is aided by being able to understand the rich history of analytic studies. A paper was published which aid in understanding the literature by simplifying some of the analytic machinery inherent in these studies. Extended MHD, like gyrokinetics, are quasineutral models. From a theoretical point of view, plasma quasineutrality approximation is best understood as the same as the magneto-quasistatic approximation of the Maxwell equations. Finally, the standard model of tokamak theory is that of instability theory. A simple dynamical systems model has been developed to better illustrate the strengths and weaknesses of this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ripple modifications to alpha transport in tokamaks

Magnetic field ripple is inherent in tokamaks since the toroidal magnetic field is generated by a finite number of toroidal field coils. The field ripple results in departures from axisymmetry that cause radial transport losses of particles and heat. These ripple losses are a serious concern for alphas near their birth speed v0 since alpha heating of the background plasma is required to make fusion reactors into economical power plants. Ripple in tokamaks gives rise to two alpha transport regimes of concern. As the slowing down time τs is much larger than the time for an alpha just born to make a toroidal transit, a regime referred to as the 1/ ν∝ τs regime is first encountered, with ν the appropriate alpha collision frequency. In this regime the radial transport losses increase as v0 τs / R , with R the major radius of the tokamak. The deleterious effect of ripple transport is mitigated by electric and magnetic drifts within the flux surface. When drift is strong enough a second regime referred to as the ν is encountered where a collisional boundary layer due to the drift plays a key role. We evaluate the alpha transport in both regimes, taking account of the alphas having a slowing down rather than a Maxwellian distribution function and their being collisionally scattered by a collision operator appropriate for alphas. Alpha ripple transport is found to be in the ν regime where it will be a serious issue for typical tokamak reactors as it will be well above the axisymmetric neoclassical level and can be large enough to deplete the alpha slowing down distribution function unless toroidal rotation is strong.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Runaway electrons in tokamaks

The generation of runaway electrons is a complex and important phenomenon that impacts many areas of plasma physics. Due to the decrease of electron collision frequency with increasing velocity, electrons under strong electric field can experience unlimited “runaway” acceleration. In tokamaks, runaway electrons can be produced in disruptions, due to the strong inductive electric field formed as the thermal energy of plasma gets rapidly lost. This population of runaway electrons can undergo an exponential growth, denoted the runaway electron avalanche, due to hard collisions between relativistic runaway electrons and low energy electrons. It is predicted that in a large tokamak device like the International Thermonuclear Experimental Reactor (ITER), a runway electron beam generated in a disruption event can potentially cause severe damage to the device, which poses a significant challenge for ITER to achieve its mission. It is therefore extremely important to seek an effective mitigation mechanism for runaway electrons. Experimental efforts have been made to study the properties of runaway electrons in tokamaks, including their generation, diffusion, and radiation. In order to understand these experimental results, extensive theoretical and simulation studies of runaway electron physics are required. The main topic of this thesis is to study the wave particle interaction associated with runaway electron beams in tokamaks. The runaway electrons can emit and absorb electromagnetic waves through resonances, and can be diffused in momentum space by the waves. Initially, we address the Cherenkov radiation of runaway electrons, which originates from the polarization of the plasma medium. The energy and momentum loss of the Cherenkov radiation can be modeled by adding a correction to the Coulomb logarithm in the collisional drag force. Subsequently, we address pitch angle scattering caused by normal modes in the plasma, which are driven unstable by the anisotropicity of the runaway electron beam. The fluctuating electromagnetic fields are found to act as a seed for the unstable normal modes. Numerical simulations show that the pitch angle scattering effect from the normal modes, mainly whistler waves, can be significantly larger than that from collisional pitch angle scattering. Finally, we present a synthetic diagnostic tool we developed to calculate the electron cyclotron emission (ECE) from the runaway electrons, and successfully reproduce the prompt growth of the ECE signal observed in DIII-D quiescent runaway electron (QRE) experiments. Within the thesis, we also present the application of the adjoint method to runaway electron research, and show the calculations of the runaway probability function (RPF) and the expected loss time (ELT). These calculations not only help depict the dynamics of runaway electrons in momentum space, but also can be used to efficiently calculate experimentally relevant quantities such as the critical electric field for runaway electron avalanche and the avalanche growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretation of the effects of electron cyclotron power absorption in pre-disruptive tokamak discharges in ASDEX Upgrade

Tokamak disruptions are events of fatal collapse of the magnetohydrodynamic (MHD) confinement configuration, which cause a rapid loss of the plasma thermal energy and the impulsive release of magnetic energy and heat on the tokamak first wall components. The physics of the disruptions is very complex and non-linear, strictly associated with the dynamics of magnetic tearing perturbations. The crucial problem of the response to the effects of localized heat deposition and current driven by external (rf) sources to avoid or quench the MHD tearing instabilities has been investigated both experimentally and theoretically on the ASDEX Upgrade tokamak. The analysis of the conditions under which a disruption can be prevented by injection of electron cyclotron (EC) rf power, or, alternatively, may be caused by it, shows that the local EC heating can be more significant than EC current drive in ensuring neoclassical tearing modes (NTMs) stability, due to two main reasons: first, the drop of temperature associated with the island thermal short circuit tends to reduce the neoclassical character of the instability and to limit the EC current drive generation; second, the different effects on the mode evolution of both the location of the power deposition relative to the island separatrix and the island shape deformation lead to less strict requirements of precise power deposition focussing. A contribution to the validation of theoretical models of the events associated with NTM is given and can be used to develop concepts for their control, relevant also for ITER-like scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A flowing liquid lithium limiter for the Experimental Advanced Superconducting Tokamak

A program involving the extensive and systematic use of lithium (Li) as a “first,” or plasma-facing, surface in Tokamak fusion research devices located at Institute of Plasma Physics, Chinese Academy of Sciences, was started in 2009. Many remarkable results have been obtained by the application of Li coatings in Experimental Advanced Superconducting Tokamak (EAST) and liquid Li limiters in the HT-7 Tokamak—both located at the institute. In furtherance of the lithium program, a flowing liquid lithium (FLiLi) limiter system has been designed and manufactured for EAST. The design of the FLiLi limiter is based on the concept of a thin flowing film which was previously tested in HT-7. Exploiting the capabilities of the existing material and plasma evaluation system on EAST, the limiter will be pre-wetted with Li and mechanically translated to the edge of EAST during plasma discharges. The limiter will employ a novel electro-magnetic pump which is designed to drive liquid Li flow from a collector at the bottom of limiter into a distributor at its top, and thus supply a continuously flowing liquid Li film to the wetted plasma-facing surface. This paper focuses on the major design elements of the FLiLi limiter. In addition, a simulation of incoming heat flux has shown that the distribution of heat flux on the limiter surface is acceptable for a future test of power extraction on EAST.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗

Continuous, edge localized ion heating during non-solenoidal plasma startup and sustainment in a low aspect ratio tokamak

Plasmas in the Pegasus spherical tokamak are initiated and grown by the non-solenoidal local helicity injection (LHI) current drive technique. The LHI system consists of three adjacent electron current sources that inject multiple helical current filaments that can reconnect with each other. Anomalously high impurity ion temperatures are observed during LHI with T i,OV ≤ 650 eV, which is in contrast to T i,OV ≤ 70 eV from Ohmic heating alone. Spatial profiles of T i,OV indicate an edge localized heating source, with T i,OV ~ 650 eV near the outboard major radius of the injectors and dropping to ~150 eV near the plasma magnetic axis. Experiments without a background tokamak plasma indicate the ion heating results from magnetic reconnection between adjacent injected current filaments. In these experiments, the HeII T i perpendicular to the magnetic field is found to scale with the reconnecting field strength, local density, and guide field, while ${{T}_{\text{i},\parallel}}$ experiences little change, in agreement with two-fluid reconnection theory. In conclusion, this ion heating is not expected to significantly impact the LHI plasma performance in Pegasus, as it does not contribute significantly to the electron heating. However, estimates of the power transfer to the bulk ion are quite large, and thus LHI current drive provides an auxiliary ion heating mechanism to the tokamak plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗