Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas

The pedestal structure, edge transport and linear MHD stability have been analyzed in a series of JET with the ITER-like wall hydrogen (H) and deuterium (D) type I ELMy H-mode plasmas. The pedestal pressure is typically higher in D than in H at the same input power and gas rate, with the difference mainly due to lower density in H than in D (Maggi et al (JET Contributors) 2018 Plasma Phys. Control. Fusion 60 014045). A power balance analysis of the pedestal has shown that higher inter-ELM separatrix loss power is required in H than in D to maintain a similar pedestal top pressure. This is qualitatively consistent with a set of interpretative EDGE2D-EIRENE simulations for H and D plasmas, showing that higher edge particle and heat transport coefficients are needed in H than in D to match the experimental profiles. It has also been concluded that the difference in neutral penetration between H and D leads only to minor changes in the upstream density profiles and with trends opposite to experimental observations. This implies that neutral penetration has a minor role in setting the difference between H and D pedestals, but higher ELM and/or inter-ELM transport are likely to be the main players. The interpretative EDGE2D-EIRENE simulations, with simultaneous upstream and outer divertor target profile constraints, have indicated higher separatrix electron temperature in H than in D for a pair of discharges at low fueling gas rate and similar stored energy (which required higher input power in H than in D at the same gas rate). The isotope dependence of linear MHD pedestal stability has been found to be small, but if a higher separatrix temperature is considered in H than in D, this could lead to destabilization of peeling-ballooning modes and shrinking of the stability boundary, qualitatively consistent with the reduced pedestal confinement in H.

08 HYDROGEN↗
DOE OSTI2021

A novel path to runaway electron mitigation via deuterium injection and current-driven MHD instability

Relativistic electron (RE) beams at high current density (low safety factor, qa) yet very low free-electron density accessed with D2 secondary injection in the DIII-D and JET tokamak are found to exhibit large-scale MHD instabilities that benignly terminate the RE beam. In JET, this technique has enabled termination of MA-level RE currents without measurable first-wall heating. This scenario thus offers an unexpected alternate pathway to achieve RE mitigation without collisional dissipation. Benign termination is explained by two synergistic effects. First, during the MHD-driven RE loss events both experiment and MHD orbit-loss modeling supports a significant increase in the wetted area of the RE loss. Second, as previously identified at JET and DIII-D, the fast kink loss timescale precludes RE beam regeneration and the resulting dangerous conversion of magnetic to RE kinetic energy. During the termination, the RE kinetic energy is lost to the wall, but the current fully transfers to the cold bulk thus enabling benign Ohmic dissipation of the magnetic energy on longer timescales via a conventional current quench. Hydrogenic (D2) secondary injection is found to be the only injected species that enables access to the benign termination. D2 injection: 1) facilitates access to low qa in existing devices (via reduced collisionality & resistivity), 2) minimizes the RE avalanche by ‘purging’ the high-Z atoms from the RE beam, 3) drives recombination of the background plasma, reducing the density and Alfven time, thus accelerating the MHD growth. Furthermore, this phenomenon is found to be accessible when crossing the low qa stability boundary with rising current, falling toroidal field, or contracting minor radius - the latter being the expected scenario for vertically unstable RE beams in ITER. While unexpected, this path scales favorably to fusion-grade tokamaks and offers a novel RE mitigation scenario in principle accessible with the day-one disruption mitigation system (DMS) of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A dimensionality reduction algorithm for mapping tokamak operational regimes using a variational autoencoder (VAE) neural network

Abstract A variational autoencoder (VAE) is a type of unsupervised neural network which is able to learn meaningful data representations in a reduced dimensional space. We present an application of VAE in identifying the operational stability boundary of tokamak plasma discharges. This model was implemented using a dataset of over 3000 discharges from the high beta tokamak-extended pulse (HBT-EP) device. We found the VAE model to be capable of forming a continuous low-dimensional operational space map and identifying the operational boundaries using a specified warning time window. By projecting the operational parameters onto the same reduced space, this provides an intuitive way for the machine operator or an automated control system to perform disruption avoidance using a relevant control actuator as a discharge approaches a boundary. Pre-programmed GPU control experiments were conducted to demonstrate this control technique using HBT-EP’s saddle control coils as a horizontal position actuator, showing the ability to avoid the oncoming disruptive event and extend the duration of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D

Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. Cryogenically frozen deuterium pellets are injected from the low-field side of the DIII-D tokamak at frequencies lower than the natural ELM frequency, and heat flux is measured by infrared cameras. Ideal MHD pedestal stability calculations show that without pellet injection, these low collisionality pedestals were limited by their current density (peeling-limited) rather than their pressure gradient (ballooning-limited). ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. Cameras observing both heat flux and D-alpha emission often find significant toroidally asymmetric striations in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets. Toroidal asymmetries at the inner strike point are similar between natural and pellet-triggered ELMs, indicating that the reduction in peak heat flux and total fluence at that location is robust for the conditions reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Critical role of current-driven instabilities for ELMs in NSTX

The impact of different extended-magnetohydrodynamic (MHD) contributions on the stability thresholds of peeling-ballooning modes in ELMing and ELM-free plasmas in the spherical tokamak NSTX is investigated with the initial value code M3D-C1. We show that ELMing discharges in NSTX are limited by resistive current-driven peeling modes, whereas non-ELMing wide-pedestal H-mode discharges are located near the ideal pressure-driven ballooning threshold. It is demonstrated that extended-MHD can lead to more reliable edge stability predictions than existing ideal-MHD models. Resistive peeling-ballooning modes are found to exist well before the ideal stability threshold is met, and kink-peeling modes exhibit considerable sensitivity to plasma resistivity. Other effects not considered in ideal-MHD models affect PB modes in NSTX in a weaker way. Gyroviscous stress appears stabilizing such that the stability boundary lies closer to the experimental point. Equilibrium rotation can suppress ideal core modes and thus isolate edge modes. These results are important for the development of a predictive pedestal model for low-aspect ratio tokamaks.

pedestal↗
DOE OSTI2024

Decoupling of peeling and ballooning thresholds for pedestal stability and reduction in ELM frequency via enhanced turbulence with edge electron cyclotron heating in DIII-D

Abstract The edge localized mode (ELM) frequency ( f ELM ) decreased by 63% when electron cyclotron heating (ECH) deposition location is shifted from ρ = 0.4 to ρ = 0.8 in DIII-D discharges where the power ratio between neutral beam injection (NBI) and ECH ( P NBI / P ECH ) is kept at ∼1. The performance of the pedestal in the ECH heated case is compared with a pure NBI reference discharge while keeping the total input power constant. All these discharges are performed at balanced input torque conditions. Furthermore, in the pure NBI discharge a strong decoupling of the peeling–ballooning (PB) thresholds is observed. The PB decoupling is preserved when the ECH is deposited at ρ = 0.8 and P NBI / P ECH ∼ 1, while the thresholds manifest a closed stability boundary when the ECH is deposited at ρ = 0.4. The inter-ELM pedestal recovery time is considerably larger for the ECH at ρ = 0.8 case. Increased pedestal turbulence is observed in beam emission spectroscopy (BES), Doppler backscattering and magnetic diagnostics for the ECH at the ρ = 0.8 case. Strong growth of a TEM-like mode is observed in BES and the mode growth is correlated with the decrease in f ELM . In view of these observations, the increased pedestal turbulence seems to be the plausible reason behind the delayed pedestal recovery following an ELM event in the ECH at ρ = 0.8 case, and the preservation of PB decoupling through temperature pedestal profile widening. TRANSP interpretative simulations show that the ECH at the ρ = 0.8 case is more susceptible to ITG/TEM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Characterization of ELM pacing via vertical jogs on DIII-D

Edge localized mode (ELM) pacing via vertical plasma oscillations or jogging has been successfully demonstrated on DIII-D. Rapid vertical movement of the plasma toward the X-point has been shown to effectively trigger ELMs. By vertically oscillating the plasma at a rate of 10 Hz, the ELM frequency increased from ~5 Hz, the natural ELM frequency in similar DIII-D discharges, to 20 Hz. Downward jogs have been observed to trigger multiple ELMs in one cycle. ELMs triggered at higher than natural frequencies lead to smaller decreases in stored energy, from 8% to as little as below 1%. As a consequence, the peak heat flux to the divertor has been observed to be reduced by a factor of ~2. In addition, a reduction in the carbon impurity concentration has been observed. During downward jogs in the lower single null (LSN) configuration, the X-point movement is slower and smaller than the top of the plasma. As a result, a reduction in the plasma cross-section and hence volume has been observed. To understand the mechanism of ELM triggering by jogging, a toy model of the edge toroidal current has been built and tested with DIII-D experiment data. The experimental data and model suggest that when the plasma moves down toward the X-point, a net positive toroidal current is locally induced in the edge region. ELITE stability analysis suggests that this current pushes the plasma state across the peeling side of the peeling–ballooning stability boundary into the unstable region triggering ELMs.

ELM pacing↗
DOE OSTI2025

Doubling fusion power with volumetric optimization in magnetic confinement fusion devices

A technique, volumetric power optimization, is presented for enhancing the power output of magnetic confinement fusion devices. Applied to a tokamak, this approach involves shifting the burning plasma region to a larger plasma volume while introducing minimal perturbations to the plasma boundary shape. This edge perturbation—squareness—is analogous to pinching and stretching the edge boundary. Stability calculations confirm that this edge alteration is compatible with maintaining plasma stability. This optimization method for optimizing fusion power output could improve the performance of magnetic confinement fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗
DOE OSTI2020

Region-Based Stability Analysis of Resilient Distribution Systems with Hybrid Grid-Forming and Grid-Following Inverters

In this paper, distribution system stability considering mixed grid-forming (GFM) and grid-following (GFL) inverters is discussed. A holistic small-signal model of the distribution feeder accommodating both GFM and GFL inverters is derived. In contrast to conventional point-by-point stability analysis, a region-based stability assessment approach is proposed and implemented to evaluate small-signal system stability with both GFM and GFL controls. Furthermore, to enhance the computational efficiency, an Artificial Intelligence (AI) assisted Kernel Ridge Regression (KRR) approach is developed to obtain the stability region boundary. Therefore, with the derived stability region, the impact of critical parameters on system stability is intuitively shown, and more importantly, the relationship between the maximum penetration level and the corresponding feasible range of the selected control parameters can be quantified. The proposed approach has been validated using MATLAB/Simulink.

grid-following↗
DOE OSTI2024

Design and Development of Stable Nanocrystalline High‐Entropy Alloy: Coupling Self‐Stabilization and Solute Grain Boundary Segregation Effects

Abstract Grain growth is prevalent in nanocrystalline (NC) materials at low homologous temperatures. Solute element addition is used to offset excess energy that drives coarsening at grain boundaries (GBs), albeit mostly for simple binary alloys. This thermodynamic approach is considered complicated in multi‐component alloy systems due to complex pairwise interactions among alloying elements. Guided by empirical and GB‐segregation enthalpy considerations for binary‐alloy systems, a novel alloy design strategy, the “ pseudo‐binary thermodynamic ” approach, for stabilizing NC‐high entropy alloys (HEAs) and other multi‐component‐alloy variants is proposed. Using Al 25 Co 25 Cr 25 Fe 25 as a model‐HEA to validate this approach, Zr, Sc, and Hf, are identified as the preferred solutes that would segregate to HEA‐GBs to stabilize it against growth. Using Zr, NC‐Al 25 Co 25 Cr 25 Fe 25 HEAs with minor additions of Zr are synthesized, followed by annealing up to 1123 K. Using advanced characterization techniques— in situ X‐ray diffraction (XRD), scanning/transmission electron microscopy (S/TEM), and atom probe tomography, nanograin stability due to coupling self‐stabilization and solute‐GB segregation effects is reported in HEAs up to substantially high temperatures. The self‐stabilization effect originates from the preferential GB‐segregation of constituent HEA‐elements that stabilizes NC‐Al 25 Co 25 Cr 25 Fe 25 up to 0.5 T m ( T m –melting temperature). Meanwhile, solute‐GB segregation originates from Zr segregation to NC‐Al 25 Co 25 Cr 25 Fe 25 GBs; this results in further stabilization of the phase and grain‐size (≈14 nm) up to ≈0.58 and ≈0.64 T m , respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗