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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Issues in Formation of Cryogenic Pellets for Fusion Applications

Cryogenic pellets are used for injection into fusion plasmas to add fuel to build up density and replace the ions lost from fusion reactions and imperfect confinement in the plasma. These pellets are formed at cryogenic temperatures with pure hydrogenic isotopes or mixtures of the isotopes. Technology to make these pellets and inject them into plasmas has been under development for many years, and various methods using freezing or desublimation have been shown to produce high-quality solid pellets suitable for injection. The throughput needed and possible impurity content from the necessary recirculation of fusion exhaust gases are two of the key issues to overcome for fusion pellet fueling systems in long-pulse burning plasmas. Here, we describe the technical challenges associated with these issues and the capability of pellet formation extruders to overcome them.Cryogenic pellets of deuterium, neon, and argon are also used in fusion tokamak devices for disruption mitigation in the form of large pellets that can be injected on demand to quickly dissipate the plasma thermal energy through radiation and add significant density in order to prevent runaway electron formation. Here, the issue is not throughput as with the fueling pellets but rather is the time it takes to form pellets of the size needed and the ability to dislodge them immediately on demand when needed to mitigate a disruption. The method used to make these pellets by desublimation is described, and examples related to how pellet size and input gas parameters affect the formation time are provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2017

Recent progress in understanding electron thermal transport in NSTX

The anomalous level of electron thermal transport inferred in magnetically confined configurations is one of the most challenging problems for the ultimate realization of fusion power using toroidal devices: tokamaks, spherical tori and stellarators. It is generally believed that plasma instabilities driven by the abundant free energy in fusion plasmas are responsible for the electron thermal transport. The National Spherical Torus eXperiment (NSTX) (Ono et al 2000 Nucl. Fusion 40 557) provides a unique laboratory for studying plasma instabilities and their relation to electron thermal transport due to its low toroidal field, high plasma beta, low aspect ratio and large ExB flow shear. Recent findings on NSTX have shown that multiple instabilities are required to explain observed electron thermal transport, given the wide range of equilibrium parameters due to different operational scenarios and radial regions in fusion plasmas. Here we review the recent progresses in understanding anomalous electron thermal transport in NSTX and focus on mechanisms that could drive electron thermal transport in the core region. The synergy between experiment and theoretical/ numerical modeling is essential to achieving these progresses. The plans for newly commissioned NSTX-Upgrade will also be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Shattered pellet injection technology design and characterization for disruption mitigation experiments

The technology to form and shoot high-Z cryogenic solid pellets mixed with deuterium using a gas gun that are shattered upon injection into a plasma has been developed at Oak Ridge National Laboratory for mitigating disruptions. This technology has been selected as the basis for the baseline disruption mitigation system on ITER. The development of shattered pellet injection systems has progressed to be able to accelerate large pellets of pure argon and neon with or without including deuterium. Impact studies have been carried out at shallow angles to determine funnel performance in guiding pellets from multiple barrels into a common injection line and across pumping breaks. The characterization of the shattered spray has also progressed with fragment size measurements as a function of pellet speed showing a strong inverse relationship. Results of these studies are reported with implications for applications on existing and future tokamak devices.

Baylor, Larry↗
DOE OSTI2019

Assessment of equilibrium field coil misalignments on the divertor footprints in NSTX-U

The presence of error fields in a tokamak device is inevitable, so it is of fundamental importance to understand their impact on the operation of the device and therefore up to what level they can be tolerated. In this paper a prediction of the impact on the magnetic footprints on the divertor plates due to the misalignment of NSTX-U equilibrium coils is presented. Resistive MHD solutions are used to predict the magnetic field perturbations experienced by the field lines due to the presence of error fields. A linear relation between the magnitude of the misalignment and the area of the magnetic footprints is found, as well as an inverse proportionality between the size of a magnetic footprint and the plasma penetration of the field lines composing the footprint. This study is intended as a starting point in the process of investigating divertor footprints associated with the plasma response to intrinsic error fields with resistive MHD codes. A basic interpretation of the results suggests that large footprints resulting from equilibrium coils misalignment may be desirable, to some extent, to decrease the peak heat load on the divertor components. Lastly, a more complete interpretation of the results presented here will involve a number of complex physics issues, for example pedestal transport and stability along with the role of various scrape-off layer and divertor processes, as well as experimental data to confirm the predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of negative triangularity shaping on energetic particle driven Alfvén eigenmodes in DIII-D

Shape variations from positive to negative triangularity may lead to improved performance regimes in a tokamak with a reduction in microturbulence as well as improved thermal confinement. The motivation of this investigation is to model and compare the neutral beam driven Alfvén eigenmode (AE) activity in two DIII-D discharges with positive and negative triangularity shaping of plasma. The simulations are performed using the linear version of the FAR3d code, which solves the reduced MHD equations for thermal plasma with addition of moment equations for the energetic ion density and parallel velocity with appropriate Landau closure relations. Overall, our numerical results indicate that for similar physical parameters, the unstable AEs observed in the negative triangularity case have lower growth rates as compared to the positive triangularity regime. Our findings may be useful to analyze the influence of the reverse-D like geometry on the AE instabilities in DIII-D and may lead to better configurations for minimizing fast ion losses in a tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2023

Development of high-performance long-pulse discharge in KSTAR

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impact of error fields and error field correction on heat fluxes in SPARC

Using a single toroidal array of coils to reduce the m,n = 2,1 resonant error field (EF) produced by the misalignment of the axisymmetric coils in SPARC can result in the enhancement of the local divertor heat fluxes. Managing high divertor heat fluxes (q ∥ $\simeq$ 10 GW m -2 ) poses a challenge for compact tokamak devices such as SPARC. The presence of non-axisymmetric magnetic field perturbations adds complexity to the problem by generating intricate 3D edge magnetic topologies that alter the heat flux distributions on the target plates. The aim of this work is to investigate the impact of the EF correction (EFC) on the heat fluxes at the divertor plates in SPARC. The MHD code M3DC1 has been used to simulate the 3D magnetic perturbations generated by the shift and tilt of several axisymmetric coils within specified tolerances, as well as from the array of EFC coils located at the midplane. Using a heuristic model that extends the concept of an axisymmetric heat flux layer to 3D plasmas, the resultant heat flux distributions is derived from magnetic footprints calculated with the MAFOT code. The results show that the EFC could either decrease or further enhance the local heat flux when used to correct the m,n = 2,1 resonant EF to enhance the core plasma performance.

3D fields↗
DOE OSTI2016

Evidence of toroidally localized turbulence with applied 3D fields in the DIII-D tokamak

New evidence indicates that there is significant 3D variation in density fluctuations near the boundary of weakly 3D tokamak plasmas when resonant magnetic perturbations are applied to suppress transient edge instabilities. The increase in fluctuations is concomitant with an increase in the measured density gradient, suggesting that this toroidally localized gradient increase could be a mechanism for turbulence destabilization in localized flux tubes. Two-fluid magnetohydrodynamic simulations find that, although changes to the magnetic field topology are small, there is a significant 3D variation of the density gradient within the flux surfaces that is extended along field lines. This modeling agrees qualitatively with the measurements. The observed gradient and fluctuation asymmetries are proposed as a mechanism by which global profile gradients in the pedestal could be relaxed due to a local change in the 3D equilibrium. In conclusion, these processes may play an important role in pedestal and scrape-off layer transport in ITER and other future tokamak devices with small applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Acoustic MEMS Sensor Array for Quench Detection of CICC Superconducting Cables

Here, a novel quench detection method using micro-electro-mechanical system (MEMS) sensor technology has been investigated in use for high temperature superconducting (HTS) conductors such REBCO tape cables. The sensor array along a superconducting cable, such as a cable-in-conduit-conductor (CICC), is installed in a cooling channel. It will allow sensitive and quick detection for a local quench of a superconducting cable. This work has confirmed that a quench of a single REBCO tape can be detected in liquid nitrogen by a MEMS piezoelectric microphone sensor. The quench detection design utilizing a MEMS sensor array method is discussed for the case of a toroidal field (TF) magnets of a fusion Tokamak device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Acoustic MEMS Sensor Array for Quench Detection of CICC Superconducting Cables

A novel quench detection method using microelectro- mechanical system (MEMS) sensor technology has been investigated in use for high temperature superconducting (HTS) conductors such REBCO tape cables. The sensor array along a superconducting cable, such as a cable-in-conduit-conductor (CICC), is installed in a cooling channel. It will allow sensitive and quick detection for a local quench of a superconducting cable. This work has confirmed that a quench of a single REBCO tape can be detected in liquid nitrogen by a MEMS piezoelectric microphone sensor. The quench detection design utilizing a MEMS sensor array method is discussed for the case of a toroidal field (TF) magnets of a fusion Tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A Drift-Kinetic Analytical Model for SOL Plasma Dynamics at Arbitrary Collisionality

A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously-derived fluid models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Reconnection of Magnetic Fields

Spacecraft observations of steady and nonsteady reconnection at the magnetopause are reviewed. Computer simulations of three-dimensional reconnection in the geomagnetic tail are discussed. Theoretical aspects of the energization of particles in current sheets and of the microprocesses in the diffusion region are presented. Terrella experiments in which magnetospheric reconnection is simulated at both the magnetopause and in the tail are described. The possible role of reconnection in the evolution of solar magnetic fields and solar flares is discussed. A two-dimensional magnetohydrodynamic computer simulation of turbulent reconnection is examined. Results concerning reconnection in Tokamak devices are also presented.

Source record↗
NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗