Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Toroid”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Enhanced plasma performance in C-2W advanced beam-driven field-reversed configuration experiments

TAE Technologies’ fifth-generation fusion device, C-2W (also called ‘Norman’), is the world’s largest compact-toroid device and has made significant progress in field-reversed configuration (FRC) plasma performance. C-2W produces record breaking, macroscopically stable, high-temperature advanced beam-driven FRC plasmas, dominated by injected fast particles and sustained in steady state, which is primarily limited by neutral-beam (NB) pulse duration. The NB power supply system has recently been upgraded to extend the pulse length from 30 ms to 40 ms, which allows for a longer plasma lifetime and thus better characterization and further enhancement of FRC performance. An active plasma control system is routinely used in C-2W to produce consistent FRC performance as well as for reliable machine operations using magnet coils, edge-biasing electrodes, gas injection and tunable-energy NBs. Google’s machine learning framework for experimental optimization has also been routinely used to enhance plasma performance. D

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Observation of a new pedestal stability regime in MAST Upgrade H-mode plasmas

Abstract The first pedestal stability and structure analysis on the new MAST Upgrade (MAST-U) spherical tokamak H-mode plasmas is presented. Our results indicate that MAST-U pedestals are close to the low toroidal mode number ( n ) peeling branch of the peeling-ballooning instability, in contrast with MAST H-mode pedestals which were deeply in the high- n ballooning branch. This offers the possibility of reaching the ELM-free quiescent H-mode (Burrell et al 2005 Plasma Phys. Control. Fusion 47 B37–B52) or high-performance super H-mode (Snyder et al 2015 Nucl. Fusion 55 083026; Snyder et al 2019 Nucl. Fusion 59 086017) regimes. In addition, the coupling between the peeling and ballooning branches is weak in MAST-U, suggesting that a path to very high pedestal pressure gradient at high density may exist with sufficient heating power. A possible explanation for the differences between MAST and MAST-U pedestal stability is given in terms of plasma shaping parameters, in particular squareness and elongation, as well as the pedestal top temperature and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Roles of non-axisymmetric perturbations in free drift vertical displacement events on EAST

Abstract The safe operation of most tokamaks, especially the large ones, relies on the feedback control of vertical displacement events (VDEs). However, most of these feedback control systems are based on axisymmetric VDE models. In this study, we use NIMROD simulations to study the role of non-axisymmetric perturbations in free drift vertical displacement events on EAST. The high- n modes in the non-axisymmetric VDE grow first, which drives the formation of high- n magnetic island chains. Subsequently, the magnetic island chains grow and overlap with each other, leading to the destruction of the magnetic flux surface, which induces a minor disruption and accelerates the start of the following major disruption. The magnetic island and the stochastic magnetic field allow the toroidally non-axisymmetric poloidal plasma current to jet towards the hoop force direction, forming finger-like and filamentary structures. Such a plasma current non-axisymmetry strongly depends on the anisotropy in the thermal transport coefficients.

Physics↗
DOE OSTI2024

Plasma rotation and diamagnetic drift effects on the resistive wall modes in the negative triangularity tokamaks

Abstract It was found previously that the negative triangularity (NT) configuration is more MHD-unstable for low n modes than the positive triangularity (PT) case, although the situation is reversed for intermediate n modes and the NT configuration becomes more stable for intermediate n modes ( n = 3 − 10 ) (Zheng et al 2021 Nucl. Fusion 61 116014). Here, n is the toroidal mode number. In this work, we extend the studies to include the rotation effects, as well as the diamagnetic drift effects, to see how the resistive wall modes (RWMs) in the NT configuration are affected as compared with the PT configuration. This is particularly motivated by noting that the wall interface with the plasma is quite different between the NT and PT configurations. It affects the plasma rotation and diamagnetic drift effects on the low n RWM. We consider the DIII-D-NT-experiment equilibrium reconstructed by the EFIT code. Based on the equilibrium g-file, the extended equilibria are constructed with the VMEC code by varying the beta values while keeping the pressure and poloidal current flux profiles basically unchanged. The bootstrap current contribution to the equilibria is taken into account with the Sauter formula. The MHD stability is then computed using the AEGIS code with the rotation and diamagnetic drift effects taken into account. We found that, although the NT configuration is less stable for n = 1 MHD modes, the rotation and diamagnetic drift stabilization effects on RWMs are more effective in the NT configuration than in the PT one. Note that even in the PT case, the stabilization of RWMs by the rotation and kinetic effects is critical. Because the low-n RWMs in the regular NT case are more unstable, the rotation and diamagnetic drift stabilization effects found in this research are important for the NT tokamak concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of fast particle experiments in the first MAST Upgrade experimental campaigns

Abstract MAST-U is equipped with on-axis and off-axis neutral beam injectors (NBI), and these external sources of super-Alfvénic deuterium fast-ions provide opportunities for studying a wide range of phenomena relevant to the physics of alpha-particles in burning plasmas. The MeV range D-D fusion product ions are also produced but are not confined. Simulations with the ASCOT code show that up to 20% of fast ions produced by NBI can be lost due to charge exchange (CX) with edge neutrals. Dedicated experiments employing low field side (LFS) gas fuelling show a significant drop in the measured neutron fluxes resulting from beam-plasma reactions, providing additional evidence of CX-induced fast-ion losses, similar to the ASCOT findings. Clear evidence of fast-ion redistribution and loss due to sawteeth (ST), fishbones (FB), long-lived modes (LLM), Toroidal Alfvén Eigenmodes (TAE), Edge Localised Modes (ELM) and neoclassical tearing modes (NTM) has been found in measurements with a Neutron Camera (NCU), a scintillator-based Fast-Ion Loss Detector (FILD), a Solid-State Neutral Particle Analyser (SSNPA) and a Fast-Ion Deuterium- α (FIDA) spectrometer. Unprecedented FILD measurements in the range of 1–2 MHz indicate that fast-ion losses can be also induced by the beam ion cyclotron resonance interaction with compressional or global Alfvén eigenmodes (CAEs or GAEs). These results show the wide variety of scenarios and the unique conditions in which fast ions can be studied in MAST-U, under conditions that are relevant for future devices like STEP or ITER.

Physics↗
DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗
DOE OSTI2024

Overview of recent results from the ST40 compact high-field spherical tokamak

ST40 is a compact, high-field ($B_{\mathrm{T0}}\unicode{x2A7D} 2.1\,\,\,\textrm{T}$) spherical tokamak (ST) with a mission to expand the physics and technology basis for the ST route to commercial fusion. The ST40 research programme covers confinement and stability; solenoid-free start-up; high-performance operating scenarios; and plasma exhaust. In 2022, ST40 obtained central deuterium ion temperatures of $9.6 \pm 0.4\ \textrm{keV}$, demonstrating for the first time that pilot plant relevant ion temperatures can be reached in a compact, high-field ST. Analysis of these high-ion temperature plasmas is presented, including a summary of confinement, transport and microstability characteristics, and energetic particle instabilities. Recent scenario development activities have focused on establishing diverted H-mode plasmas across a range of toroidal fields and plasma currents, along with scenarios with high non-inductive current fractions. In future operations, beginning in 2025, a 1 MW dual frequency (104/137 GHz) electron cyclotron (EC) system will be installed to enable the study of EC and electron Bernstein wave plasma start-up and current drive. Predictive modelling of the potential performance of these systems is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recent advance progress of HL-3 experiments

Since the first plasma realized in 2020, a series of key systems on HL-3 (known as HL-2M before) tokamak have been equipped/upgraded, including in-vessel components (the first wall, lower divertor, and toroidal cryogenic/water-cooling/baking/glow discharge systems, etc.), auxiliary heating system of 11 MW, and 28 diagnostic systems (to measure the plasma density, electron temperature, radiation, magnetic field, etc.). Magnet field systems were commissioned firstly for divertor plasma discharges. During the 2nd experimental campaign of HL-3 tokamak, several great progresses have been achieved. Firstly, the successful operation with plasma current larger than 1 MA was achieved under a divertor configuration. Secondly, the advanced divertor concept with two distinct snowflake configurations was realized. It is found that the distribution of ion saturation current and heat flux on bottom plate becomes wide due to magnetic surface expansion, demonstrating the advantage of such configuration in the heat flux mitigation. In addition, using the combination of NBI, ECRH and LHCD, the standard sawtoothing high confinement mode of megampere plasma was firstly accessed on the HL-3. The successful commissioning of HL-3 is beneficial for the initial operation of ITER.

1 MA plasma↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2025

First fluctuation measurements using an Imaging Neutral Particle Analyzer on DIII-D

A recent upgrade to the Imaging Neutral Particle Analyzer (INPA) on DIII-D has allowed for the first fluctuation measurements using an INPA to be taken during a neoclassical tearing mode (NTM). The INPA signal tracked the mode over 150 ms as the mode frequency dropped to zero, capturing both the NTM with poloidal and toroidal mode numbers $m/n=2/1$ and a $3/2$ mode at double the frequency. Analysis shows that the signals originate from charge exchange events near the edge of the plasma, and relative fluctuation amplitudes are greater than 25% for the duration of the NTM. Filtered signals show frequency beating patterns that are phase-space dependent. Simulated signal is dominated by prompt transport from the neutral beams to the INPA sightline, while the contribution from the slowing down distribution is significantly lower. Simulated measurements in the range of pitches that the diagnostic is sensitive to (0.5≤|v∥/v|≤ 0.75) show the signal is dominated by trapped orbits that pass through magnetic islands near the edge of the plasma. Calculations of expected fluctuation levels show that only direct interaction with the NTM can provide the strong relative fluctuation levels seen experimentally. The prompt nature of the orbits and thin radial layer found to contribute to synthetic signals suggest INPA passive data may be used to measure the perturbations of confined orbits on a single pass through a plasma instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗