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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2022

Hydration-driven stabilization and volume collapse of grain boundaries in Mg2SiO4 forsterite predicted by first-principles simulations

Grain boundaries in mantle minerals are of critical importance to geophysical and geochemical processes of the Earth’s interior. One of the fundamental issues is to understand how the water (H 2 O) component influences the properties of grain boundaries in silicate materials. Here, we report the results of the structure and stability of several tilt grain boundaries in Mg 2 SiO 4 forsterite over the pressure range 0 to 15 GPa using density functional theory-based first-principles simulations. Here, the results suggest greater energetic stability and hydration-driven volume collapse (negative excess volume) at zero pressure for the majority of hydrous grain boundaries relative to the anhydrous (dry) ones. All the hydrous grain boundaries become increasingly favorable at elevated pressures as the calculated hydration enthalpy systematically decreases with increasing pressure. The hydrous components at the interfacial regions are predominantly in the hydroxyl form and, to a lesser extent, in the molecular H 2 O form. Their calculated ratio ranges from 1.6 to 8.7 among the different grain boundary configurations. Our structural analysis also reveals that the hydroxyls are bound to either both Mg and Si or to Mg only. In comparison, the molecular species are bound only to Mg sites. Besides direct oxygen-hydrogen bonding, intermolecular hydrogen bonding becomes important with compression. On the basis of our results, we suggest that local atomic rearrangements caused by dissociative adsorption of water facilitate efficient compaction of the boundary interfaces, which, in turn, results in greater relative stability of hydrous grain boundaries. This means that water prefers to be incorporated within the grain boundaries over the bulk of silicate materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Dry Deposition Methods Based on Turbulence Kinetic Energy: 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single‐Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Organizing bacterial vortex lattices by periodic obstacle arrays

Recent experiments have shown that the complex spatio-temporal vortex structures emerging in active fluids are susceptible to weak geometrical constraints. This observation poses the fundamental question of how boundary effects stabilize a highly ordered pattern from seemingly turbulent motion. Here we show, by a combination of continuum theory and experiments on a bacterial suspension, how artificial obstacles guide the flow profile and reorganize topological defects, which enables the design of bacterial vortex lattices with tunable properties. To this end, the continuum model is extended by appropriate boundary conditions. Beyond the stabilization of square and hexagonal lattices, we also provide a striking example of a chiral, antiferromagnetic lattice exhibiting a net rotational flow, which is induced by arranging the obstacles in a Kagome-like array. Recent experiments showed that weak geometrical constraints can organize topological defects in turbulent bacterial suspensions. Here, the authors use a continuum model to study the connection between symmetry and stability of emergent vortex patterns and the geometry of constraining pillar arrays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Stabilized Nanocrystalline Alloys: The Intersection of Grain Boundary Segregation with Processing Science

Processing science for nanocrystalline metals has largely focused on far-from-equilibrium methods that can generate many grain boundaries with excess defect energy. Conversely, the science of stabilizing nanocrystalline alloys has largely focused on the lowering of that excess defect energy through grain boundary segregation, bringing nanocrystalline structures closer to equilibrium. With increasing technological adoption of stabilized nanocrystalline alloys, there is a substantial need for research at the intersection of these two fields. This review lays out the basic thermodynamic issues of the two subfields and surveys the literature on the most common processing methods, including severe plastic deformation, ball milling, physical vapor deposition, and electrodeposition. We provide an overview of studies that have examined grain boundary segregation through each of these methods and identify general themes. We conclude that there is substantial scope for more systematic work at the intersection of these fields to understand how nonequilibrium processing affects grain boundary segregation.

Materials Science↗