Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tokamak plasma physics”

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

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗
DOE OSTI2022

Data-driven model for divertor plasma detachment prediction

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Halo current rotation scaling in post-disruption plasmas

Abstract Halo current (HC) rotation during disruptions can be potentially dangerous if resonant with the structures surrounding a tokamak plasma. We propose a drift-frequency-based scaling law for the rotation frequency of the asymmetric component of the HC as a function of toroidal field strength and plasma minor radius ( f rot ∝ 1/ B T a 2 ). This scaling law is consistent with results reported for many tokamaks and is motivated by the faster HC rotation observed in the HBT-EP tokamak. Projection of the rotation frequency to ITER and SPARC parameters suggest the asymmetric HC rotation will be on the order of 10 Hz and 60 Hz, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A physics-constrained deep learning surrogate model of the runaway electron avalanche growth rate

A surrogate model of the runaway electron avalanche growth rate in a magnetic fusion plasma is developed. This is accomplished by employing a physics-informed neural network (PINN) to learn the parametric solution of the adjoint to the relativistic Fokker–Planck equation. The resulting PINN is able to evaluate the runaway probability function across a broad range of parameters in the absence of any synthetic or experimental data. This surrogate of the adjoint relativistic Fokker–Planck equation is then used to infer the avalanche growth rate as a function of the electric field, synchrotron radiation and effective charge. Predictions of the avalanche PINN are compared against first principle calculations of the avalanche growth rate with excellent agreement observed across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Role of microtearing mode in DIII-D and future high- β p core plasmas

This paper provides modeling evidence for a new variant of microturbulence instability that regulates the electron thermal fluxes in the transport barrier region of DIII-D "high β p " discharges. This scenario is being investigated as a possible scenario for future steady-state reactors, and thus it is essential to understand the physics processes which govern its performance in order to confidently extrapolate to future devices. Additional analysis examining the extrapolation of this scenario to ITER shows that while the new microtearing mode identified here is unlikely to be unstable in that regime (due to increased temperatures and reduced collisionality), the large shift in magnetic axis (Shafranov shift) that arises in these plasmas may still favor stabilization of common ballooning instabilities which peak on the low field side, but destabilize electrostatic slab modes on the high field side.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Two-fluid model of rf current condensation in magnetic islands

The stabilization of tearing modes with rf waves is subject to a nonlinear effect, termed rf current condensation, that has the potential to greatly enhance and localize current driven within magnetic islands. Here we extend previous investigations of this effect with a two fluid model that captures the balance of diffusive and thermal equilibration processes within the island. We show that the effective power and resulting strength of the condensation effect can be greatly enhanced by avoiding collisional heat loss to the ions. The relative impact of collisions on the overall power balance within the island depends on the ratio of the characteristic diffusion timescale and the electron–ion equilibration time, rather than the latter alone. Although relative heat loss to ions increases with island size, the heating efficiency does as well. In particular, we show that the latter safely dominates for large deposition profiles, as is typically the case for lower hybrid current drive. This supports the possibility of passive stabilization of neoclassical tearing modes without the precise aiming of the rf waves required for electron cyclotron current drive stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The role of toroidal rotation in the very high energy confinement quality observed in super H-mode experiments on DIII-D

In this paper, we report the key role that toroidal rotation and the related E × B shear physics played in the very high energy confinement quality (H 98y2 > 1.5) of super H-mode experiments on DIII-D. Experiments show that the energy confinement quality decreases when toroidal rotation decreases due to the decreased externally controlled torque per particle. Meanwhile, the total pedestal pressure in the experiments remains very high during the rotation and confinement quality change. TGYRO transport modeling suggests the contribution from rotation in the E × B shear is responsible for the confinement quality in excess of standard H-mode (H 98y2 ~ 1). CGYRO gyrokinetic simulations reveal the governing physics in the core plasma of super H-modes: significant up-shift of nonlinear the ITG critical gradient is observed when applying E × B shear physics in the modeling based on experimental data. Furthermore, the effects of other physical parameters and contribution from pedestal height, which may play minor roles in this study, are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Measurements of multiple heat flux components at the divertor target by using surface eroding thermocouples (invited)

The Surface Eroding Thermocouple (SETC) is a robust diagnostic utilized in DIII-D to provide fast, edge-localized modes (ELMs) resolved heat flux measurements, in particular in geometric regions that are too shadowed for traditional infrared thermography. In order to further investigate the power dissipation in the divertor region, a combination of flush-mounted and recessed SETCs was developed to assess the effect on surface heating from non-charged particles at the divertor target. Utilizing the Divertor Materials Evaluation System sample exposure platform, the first demonstration of the feasibility of using this new method to distinguish between the heat flux from charged particles and that from neutrals and radiative heating was achieved. This paper details the process of using the combination of flush SETCs and recessed SETCs to measure the multiple heat flux components at the divertor target and further discusses how to determine two important ratios, α (ratio of heat flux from charged particles deposit on recessed SETC to that deposit on flush SETC) and β (ratio of heat flux from non-charged particles deposit on recessed SETC to that deposit on flush SETC), in the estimation of the heat flux from non-charged particle sources. Using a time dependent ratio α, it was found that ~50% of the total incident heat flux is attributable to the non-charged particles in the fully detached open divertor in DIII-D. Finally, the new application of similar SETC diagnostics in the Small Angle Slot divertor with a V-like configuration and partial tungsten coated surface (SAS-VW) is also introduced.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Understanding the negative triangularity ELM trigger and ELM free state on DIII-D with ECE-imaging

Here, the Electron Cyclotron Emission Imaging (ECEI) diagnostic was used to observe a finite-n interchange mode structure in the edge of negative triangularity shaped plasmas on DIII-D. At a small negative triangularity (δ u = -0.2), the plasma is in the H-mode with ELMs that are triggered by a low-n interchange mode. At a larger negative triangularity (δ u = -0.4) and low NBI power (2 MW), a dithering oscillation is observed that is triggered by a low-n interchange mode, whereas at higher NBI power (>2 MW), the edge reverts to L-mode and the low-n interchange mode is present continuously. In all cases, the edge pressure gradient is clamped when the interchange mode is present. It is concluded that the low-n interchange mode prevents the plasma from transitioning to H-mode at a large negative triangularity. This agrees with linear BOUT++ simulations which suggest that the interchange-type MHD can be a resistive ballooning mode whereby resistivity can significantly increase the finite-n ballooning mode growth rate. The absence of H-mode at large negative triangularity can, thus, be explained by the excitation of low-n pressure driven resistive ballooning modes in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics and applications of dusty plasmas: The Perspectives 2023

Dusty plasmas are electrically quasi-neutral media that, along with electrons, ions, neutral gas, radiation, and electric and/or magnetic fields, also contain solid or liquid particles with sizes ranging from a few nanometers to a few micrometers. These media can be found in many natural environments as well as in various laboratory setups and industrial applications. As a separate branch of plasma physics, the field of dusty plasma physics was born in the beginning of 1990s at the intersection of the interests of the communities investigating astrophysical and technological plasmas. An additional boost to the development of the field was given by the discovery of plasma crystals leading to a series of microgravity experiments of which the purpose was to investigate generic phenomena in condensed matter physics using strongly coupled complex (dusty) plasmas as model systems. Finally, the field has gained an increasing amount of attention due to its inevitable connection to the development of novel applications ranging from the synthesis of functional nanoparticles to nuclear fusion and from particle sensing and diagnostics to nano-contamination control. The purpose of the present perspectives paper is to identify promising new developments and research directions for the field. As such, dusty plasmas are considered in their entire variety: from classical low-pressure noble-gas dusty discharges to atmospheric pressure plasmas with aerosols and from rarefied astrophysical plasmas to dense plasmas in nuclear fusion devices. Both fundamental and application aspects are covered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Performance characterization of x-ray crystal spectroscopy highly oriented pyrolytic graphite reflectors based on x-ray diffractometry experiments

The use of Highly Oriented Pyrolytic Graphite (HOPG) reflectors is often proposed in the design of X-ray Crystal Spectroscopy (XCS) diagnostic systems for the next-generation tokamak devices, including the ITER project. Here, this study introduces an experimental study based on the X-Ray Diffractometry (XRD) method to evaluate the performance of HOPG reflectors. The experimental method provides both the angular responses and the reflectivities of the HOPG reflectors. A demonstrative XRD experiment is conducted, and the details of the experiment are introduced. This method enables precise studies on HOPG reflectors, facilitating the design of XCS diagnostic systems for future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Development of a tomographic reconstruction method for axisymmetric D α emission profiles in the ITER plasma boundary

A tomographic method for reconstructing the axially symmetric D α emission profile in the ITER scrape-off layer and divertor from D α intensity measurements in the fields of view (FoV) of Vis/IR TV and Divertor Impurity Monitor diagnostics is under development. The method takes into account the strong background signal due to stray light from divertor produced by reflections from the metal wall of the vacuum chamber. The method allows to filter out the reflections in D α signals measured in the main chamber. Here, the accuracy of this method is tested within the framework of a synthetic diagnostic, which uses the Raysect and Cherab numerical codes for ray-tracing simulations. The D α emission profiles, simulated with the SOLPS and OEDGE codes, are used as synthetic experimental data. Uncertainty in light reflection properties of the first wall (FW) is the major factor affecting the accuracy of the D α emission reconstruction. The dependence on the uncertainty in the FW surface roughness and on the ratio of diffusive to specular reflection of the error in recovering useful D α signals (the direct light without reflections) in the FoV of the ITER main chamber H-alpha and Visible Spectroscopy Diagnostic is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗