SEARCH · Engineering Papers

Results for “Surface fluxes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Impact of β n and spectrum of n = 1 applied fields on fast ion losses in DIII-D

Abstract The effect of n = 1 magnetic perturbations (MPs) on prompt fast ion losses in DIII-D has been investigated using the light ion beam probe technique. The effect of normalized beta, β n , and the plasma response to the MPs is studied for several MP spectra. Magnetics data show a strong dependence of plasma response to β n that depends on the phase difference in perturbation coils, Δ ϕ U L . For most phases, the response increases with β n , however, the response is suppressed for phases between 180 ∘ and 240 ∘ . Experimental data from scintillator based fast ion loss detectors shows that relative fluctuation of lost ions (20%–30% of the steady signal) does not diminish across the L- to H-mode transition for Δ ϕ U L = 240 ∘ , despite the 34% decrease in radial response field and 50% decrease in poloidal response field. Simulations of the DIII-D discharges using M3D-C1 and ASCOT5 find that prompt losses from neutral beam injection in each investigated case hit the first wall at several concentrated locations: the vessel floor, midplane diagnostic ports, vessel ceiling, and inner wall. Simulated losses from co-injection neutral beams are born outside the last closed flux surface, corresponding to radial orbit displacements of 3–6 cm due to the MPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Comparison of MHD stability properties between QH-mode and ELMy H-mode plasmas by considering plasma rotation and ion diamagnetic drift effects

Magnetohydrodynamic (MHD) stability at tokamak edge pedestal in a quiescent H-mode (QH-mode) and type-I ELMy H-mode plasmas in DIII-D experiment was analyzed by considering plasma rotation and ion diamagnetic drift effects. QH-mode plasma is marginally stable to kink/peeling mode (K/PM), but ELMy H-mode one is almost unstable to peeling-ballooning mode (PBM). It was identified that there are three physics features responsible for the difference in the MHD stability properties between QH-mode plasma and ELMy H-mode one. These are the distance of pedestal foot from the last closed flux surface (LCFS), the amount of the ion diamagnetic drift frequency at pedestal, and impact of coupled rotation and ion diamagnetic drift effects. These features were confirmed through the numerical experiments that the stability properties of the QH-mode plasma can be changed to that of the ELMy H-mode one by shifting the plasma profiles inward in the radial direction and halving the ion diamagnetic drift frequency. The reasons of the change in the stability properties are thought as that K/PM is stabilized due to the inward shift of the bootstrap current profile, and PBM is destabilized due to the reduction of the coupled rotation and ion diamagnetic drift stabilizing effect. Importance of these features was validated through numerical experiments with experimental data of other QH-mode plasmas in DIII-D. All the results show that MHD stability properties of QH-mode plasma can be obtained in case that pedestal foot is close to LCFS, ion diamagnetic drift frequency is large due to high ion temperature, and strong rotation shear exists near pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Design of passive and structural conductors for tokamaks using thin-wall eddy current modeling

Abstract A new three-dimensional electromagnetic modeling tool ( ThinCurr ) has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vessel and other conducting structural elements with realistic material resistivities. This model was leveraged to support the design of a passive runaway electron mitigation coil (REMC), studying the effect of various design parameters, including coil resistivity, current quench duration, and plasma vertical position, on the effectiveness of the coil. The REMC is a non-axisymmetric coil designed to passively drive large non-axisymmetric fields during the plasma disruption thereby destroying flux surfaces and deconfining RE seed populations. These studies indicate that current designs should apply substantial 3D fields at the plasma surface during future plasma current disruptions as well as highlight the importance of having the REMC conductors away from the machine midplane in order to ensure they are robust to off-normal disruption scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Overview of the first Wendelstein 7-X long pulse campaign with fully water-cooled plasma facing components

After a long device enhancement phase, scientific operation resumed in 2022. The main new device components are the water cooling of all plasma facing components and the new water-cooled high heat flux divertor units. Water cooling allowed for the first long-pulse operation campaign. A maximum discharge length of 8 min was achieved with a total heating energy of 1.3 GJ. Safe divertor operation was demonstrated in attached and detached mode. Stable detachment is readily achieved in some magnetic configurations but requires impurity seeding in configurations with small magnetic pitch angle within the edge islands. Progress was made in the characterization of transport mechanisms across edge magnetic islands: Measurement of the potential distribution and flow pattern reveals that the islands are associated with a strong poloidal drift, which leads to rapid convection of energy and particles from the last closed flux surface into the scrape-off layer. Using the upgraded plasma heating systems, advanced heating scenarios were developed, which provide improved energy confinement comparable to the scenario, in which the record triple product for stellarators was achieved in the previous operation campaign. However, a magnetic configuration-dependent critical heating power limit of the electron cyclotron resonance heating was observed. Exceeding the respective power limit leads to a degradation of the confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Transport-driven toroidal rotation with general viscosity profile

Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

Predictive modeling of NSTX discharges with the updated multi-mode anomalous transport module

Abstract The objective of this study is twofold: firstly, to demonstrate the consistency between the anomalous transport results produced by updated Multi-Mode Model (MMM) version 9.0.4 and those obtained through gyrokinetic simulations; and secondly, to showcase MMM’s ability to predict electron and ion temperature profiles in low aspect ratio, high beta NSTX discharges. MMM encompasses a range of transport mechanisms driven by electron and ion temperature gradients, trapped electrons, kinetic ballooning, peeling, microtearing, and drift resistive inertial ballooning modes. These modes within MMM are being verified through corresponding gyrokinetic results. The modes that potentially contribute to ion thermal transport are stable in MMM, aligning with both experimental data and findings from linear CGYRO simulations. The isotope effects on these modes are also studied and higher mass is found to be stabilizing, consistent with the experimental trend. The electron thermal power across the flux surface is computed within MMM and compared to experimental measurements and nonlinear CGYRO simulation results. Specifically, the electron temperature gradient modes (ETGM) within MMM account for 2.0 MW of thermal power, consistent with experimental findings. It is noteworthy that the ETGM model requires approximately 5.0 ms of computation time on a standard desktop, while nonlinear CGYRO simulations necessitate 8.0 h on 8 K cores. MMM proves to be highly computationally efficient, a crucial attribute for various applications, including real-time control, tokamak scenario optimization, and uncertainty quantification of experimental data.

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

Successful prediction of tokamak transport in the L-mode regime

A long standing shortfall in the predicted L-mode edge energy transport by reduced quasi-linear models of gyrokinetic turbulent transport has been resolved. The improved model TGLF-SAT2 has higher fidelity to gyrokinetic simulations of the electron-scale contribution to the electron energy transport and the ion-scale flux surface shape dependence of energy transport. The success of TGLF-SAT2 in predicting the L-mode and Ohmic edge profiles is critical to whole pulse simulation and opens the door to prediction of the H-mode power threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Validation of a synthetic fast ion loss detector model for Wendelstein 7-X

Abstract We present the first validated synthetic diagnostic for fast ion loss detectors (FILDs) in the Wendelstein 7-X (W7-X) stellarator. This model has been developed on, and validated against experimental data from, a FILD provided by the National Institute for Fusion Science (NIFS-FILD), with potential future applicability to the existing Faraday Cup FILD (FC-FILD) on W7-X as well as the scintillating FILD (S-FILD) currently under development. A workflow combining Monte Carlo codes BEAMS3D and ASCOT5 is used to track fast ions produced by neutral beam injection from the moment of ionization until they are thermalized or lost from the last closed flux surface, and from there to a virtual plane which serves as a projection of the entrance aperture to the FILD. Simulations in ASCOT5 are analyzed via a geometric method to determine the probability of transmission through the FILD aperture and onto the detector as a function of normalized momentum, pitch angle, gyrophase, and position at the virtual plane. This probability is then applied to the simulated ions arriving from the plasma, producing a simulated signal from a computationally tractable number of simulated fast ions. Simulated signals are presented for two W7-X experiments with neutral beam injection and quantitatively compared with experimental measurements from the NIFS-FILD diagnostic. An estimate of the frequency of charge-exchange with neutral particles in the edge is performed, and it is found that this process may have a significant impact on the measured signals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Quasilinear theory and modelling of gyrokinetic turbulent transport in tokamaks

The theory, development, and validation of reduced quasilinear models of gyrokinetic turbulent transport in the closed flux surface core of tokamaks is reviewed. In combination with neoclassical collisional transport, these models are successful in accurately predicting core tokamak plasma temperature, density, rotation, and impurity profiles in a variety of confinement regimes. Refined experimental tests have been performed to validate the predictions of the quasilinear models, probing changes in the dominant gyrokinetic instabilities, as reflected in fluctuation measurements, cross-phases, and transport properties. These tests continue to produce a deeper understanding of the complex mix of instabilities at both electron and ion gyroradius scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Electron transport barrier and high confinement in configurations with internal islands close to the plasma edge of W7-X

The low magnetic shear in the Wendelstein 7-X (W7-X) stellarator makes it feasible to shape the separatrix by the large islands constituting an island-divertor, and this can be exploited to access various magnetic configurations, including samples of different internal island sizes and locations. To investigate the configuration effects on the plasma confinement, a configuration scan was performed by changing the coil currents to vary the rotational transform between values 5/4 and 5/6 at the plasma boundary with different power levels (2, 4, 6 MW) of electron cyclotron resonance heating (ECRH) at a maximum plasma density of 8 × 10 19 m -2 . neutral beam injection (NBI) heating was also applied during some configurations of the scan to create a density ramp and access high densities beyond the X2 ECRH cutoff. For the magnetic configurations, where the 5/5 and 5/6 island chains were moved closer to separatrix but remaining inside the last closed flux surface, the electron cyclotron emission shows that an electron temperature, T e , pedestal develops already during ECRH heated plasma buildup phase indicating a transport barrier, and the barrier sustains irrespective of changed plasma heating conditions such as NBI in the later part of discharge. The transport barrier is broken by subsequent fast crashes, observed with multiple plasma diagnostics with characteristics such as tokamak edge localized modes, and the corresponding crash amplitude and frequency vary with plasma pressure. The impact of the transport barrier on plasma confinement can be seen through the increased core T e profile, which could be responsible for the overall increase in the stored diamagnetic energy by approximately 10% for these configurations. After the plasma heating is terminated, a backwards transition to a degraded confinement state is also observed. These observations indicate a configuration triggered high confinement mode in low shear W7-X. This work focuses on the occurrence of this transport barrier for different magnetic configurations and its relation to internal magnetic islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The formation of an radial edge electric field due to finite ion orbit width effects is the possible root cause of the H-mode edge

Abstract Full orbit-following simulations of thermal ions show that finite ion-orbit width effects create charge separation near the last closed flux surface (LCFS) which generates a localized radial electric field. Experimentally, edge electric fields are observed in H-mode plasmas and they are necessary for the edge turbulence suppression via the E × B flow shear mechanism. Confined trapped (and to a lesser extent co-passing) ions near the plasma edge form a positive charge distribution outside the LCFS, while thermal electrons are tied more tightly to field lines owing to their small mass and are poorly confined outside the LCFS, hence charge neutrality is violated outside the LCFS. A large number of reported observations from spherical and conventional tokamaks support the results from the simulations although the simulations were not performed fully self consistently. The results suggest ways to lower the H-mode power threshold and optimize the H-mode plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

First experimental confirmation of island SOL geometry effects in a high radiation regime on W7-X

Abstract This work characterizes the detachment behavior and radiation characteristics of the lowiota configuration in the Wendelstein 7-X (W7-X) stellarator. The island scrape-off layer (SOL) of thelow iota has a poloidal mode number of 6 islands surrounding the last closed flux surface (LCFS).The island geometry of the low iota configuration is significantly different to that of the standardmagnetic field configuration, whose detachment characteristics have already been described in previouswork[2, 3, 4]. Experimental results show that the radiation pattern in the low iota configurationis starkly different to that of the standard magnetic field configuration, with radiation concentratedat the island SOL O-points, rather than the X-points. Additionally, this O-point localized radiationpattern is associated with unstable detachment, with both radiation oscillations in experiments and thelack of a self-consistent plasma solution at high radiated power fraction in EMC3-Eirene simulations.EMC3-Eirene simulations are used to understand the radiation distribution. It was found that the O-point localized radiation arises first from local impurity accumulation near the parallel flow stagnation,which is located close to the geometrical center of the island (”O-point”). The local cooling in thisregion leads to plasma condensation in the islands in closest magnetic connection to the divertor targetplates. The heat source to this region of the island, which is thermally isolated from the upstream heatsource in terms of parallel transport, must arise via perpendicular heat transport. This heat source isexpected to be large for the low iota configuration due to its very small internal island field line pitch.This work highlights the importance (complementary to previous work, e. g. [5, 6]) of the internalisland field line pitch not only on the radiation pattern, but also the detachment performance of theisland divertor.

Winters, Victoria (ORCID:0000000181087774)↗

DOE OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Modeling the effect of MHD activity on runaway electron generation during SPARC disruptions

Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 – an extended MHD code with a RE fluid model – we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC – a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium (D 2 ) injection, are considered. Here, our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection (2-5 × 10 21 atoms), while combined D 2 + Ne injection (2 × 10 21 Ne atoms; 1.8 × 10 22 D 2 molecules) produces a lower RE current (<2 MA). With D 2 + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.

Datta, Rishabh [Massachusetts Inst. of Technology ↗