Search NASA⌕ Search

SEARCH · Search NASA

Results for “high gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30

Transition from retrograde to prograde drift instabilities in a magnetron microdischarge

Rotating plasma structures, or “spokes,” in magnetized discharges characterized by perpendicular electric and magnetic fields have been seen in a growing number of studies and are believed to be the result of gradient-driven drift instabilities. Under certain conditions, we have shown the spoke’s rotation to be opposite to the E B direction, i.e., retrograde in its expected direction. Recently, we have linked such counter-intuitive rotation to a local inversion of the electric field. Here, we give further experimental evidence for this inversion and attempt to provide an explanation for a relatively distinct transition seen between retrograde and positive (prograde) drift. In the experiments, a partially magnetized plasma forms inside a magnetron device of mm radius operated with argon. Discharge current–voltage measurements are acquired for a range of argon fill pressure and inter-electrode spacing. We find two branches of operation—a low current branch of negative resistance, coinciding with the retrograde spoke rotation, and a higher current branch of positive resistance, coincident with prograde spoke rotation. We postulate that at low discharge currents, high magnetic field confinement leads to a large density gradient, causing more electron transport to the anode than that demanded by the external circuit. At higher currents, anomalous axial electron transport (across the magnetic field lines) becomes dominant, eliminating the conditions for field inversion. Here, the current thresholds for the field inversion are found to be sensitive to the argon pressure and inter-electrode spacing.

42 ENGINEERING↗

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗

Dynamic Structure and Phase Behavior of a Block Copolymer Electrolyte under dc Polarization

There is an important consideration when designing lithium battery electrolytes for advanced applications is how the electrolyte facilitates ion transport at fast charge and discharge rates. Large current densities are accompanied by large salt concentration gradients across the electrolyte. Nanostructured composite electrolytes have been proposed to enable the use of high energy density lithium metal anodes, but many questions about the interplay between the electrolyte morphology and the salt concentration gradient that forms under dc polarization remain unanswered. To address these questions, we use an in situ small-angle X-ray scattering technique to examine the nanostructure of a polystyrene- block -poly(ethylene oxide) copolymer electrolyte under dc polarization with spatial and temporal resolution. In the quiescent state, the electrolyte exhibits a lamellar morphology. The passage of ionic current in a lithium symmetric cell leads to the formation of concurrent phases: a disordered morphology near the negative electrode, lamellae in the center of the cell, and coexisting lamellae and gyroid near the positive electrode. The most surprising result of this study was obtained after the applied electric field was turned off: a current-induced gyroid phase grows in volume for 6 h in spite of the absence of an obvious driving force. We show that this reflects the formation of localized pockets of salt-dense electrolyte, termed concentration hotspots, under dc polarization. Our methods may be applied to understand the dynamic structure of composite electrolytes at appreciable current densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strong yet ductile nanolamellar high-entropy alloys by additive manufacturing

Additive manufacturing produces net-shaped components layer by layer for engineering applications. The additive manufacture of metal alloys by laser powder bed fusion (L-PBF) involves large temperature gradients and rapid cooling, which enables microstructural refinement at the nanoscale to achieve high strength. However, high-strength nanostructured alloys produced by laser additive manufacturing often have limited ductility. In this work, we use L-PBF to print dual-phase nanolamellar high-entropy alloys (HEAs) of AlCoCrFeNi 2.1 that exhibit a combination of a high yield strength of about 1.3 gigapascals and a large uniform elongation of about 14 percent, which surpasses those of other state-of-the-art additively manufactured metal alloys. The high yield strength stems from the strong strengthening effects of the dual-phase structures that consist of alternating face-centred cubic and body-centred cubic nanolamellae; the body-centred cubic nanolamellae exhibit higher strengths and higher hardening rates than the face-centred cubic nanolamellae. The large tensile ductility arises owing to the high work-hardening capability of the as-printed hierarchical microstructures in the form of dual-phase nanolamellae embedded in microscale eutectic colonies, which have nearly random orientations to promote isotropic mechanical properties. The mechanistic insights into the deformation behaviour of additively manufactured HEAs have broad implications for the development of hierarchical, dual- and multi-phase, nanostructured alloys with exceptional mechanical properties.

36 MATERIALS SCIENCE↗

Development of Transformative Preparation Methods to Push up High Q&G Performance of FRIB Spare HWR Cryomodule Cavities

The FRIB accelerator project construction, a top priority of US nuclear science, was completed in January 2022, and is now moving to user operation. The stable and reliable operation of the accelerating cryomodules is essential in achieving/fulfilling DOE and user expectations. So far, FRIB cryomodules meet all FRIB specifications for cavity performance. However, during the lifetime of machine operation, degradation of cryomodule performance is possible, as reported in similar operating facilities (CEBAF, SNS). If cryomodule degradation is observed at FRIB, the under-performing cryomodule will require replacement/maintenance. In effort to manage operational reliability, FRIB plans to construct a 0.53 half-wave cryomodule to serve as an active spare. In a parallel effort, FRIB will also work toward increasing operational Q and gradient of spare cryomodule cavities to gain an overall performance margin to support future operational reliability. The current FRIB cavity designs have a potential to operate at gradients higher than 8 MV/m, but are currently limited by field emission (FE) and/or high field Q slope (HFQS); known issue in buffered chemical polished (BCP) treated cavities. The proposal looks to develop transformative surface preparation treatments to improve the operational gradient of spare cryomodules higher than 10 MV/m while maintaining high Q. Thus, increasing operational margin by 30 - 50%. With the goal to improve operational reliability set, the proposal will investigate multiple objectives as possible paths forward to achieve an overall increase in cavity performance and gain a better understanding of SRF limiting mechanisms. The proposal will study the application of different chemical surface treatments to 0.53 half-wave cavities, with the addition of low temperature bakes (LTB), and measure their effects on accelerating performance. Proposed chemical treatments to be explored in this proposal include conventional EP acid mixtures, as well as innovated EP and BCP acid mixtures designed to simplify processing paths in migrating FE and HFQS. The proposed transformative treatment wet N-doping also has the potential to replicate recent advancements in SRF technology relating to nitrogen doping and high Q operation without the requirement for an ultra-high vacuum annealing furnace; currently being developed at FNAL and JLAB. In parallel, high Q performance relating to flux trapping will be investigated with the installation of a second layer of magnetic shielding in the vertical test Dewar. The research objectives presented in the proposal, and their corresponding effects on cavity performance, will provide essential knowledge and future guidance to the SRF community and provide possible paths for future SRF based projects and applications.

43 PARTICLE ACCELERATORS↗

High-$Q_0$ treatment of CEBAF 1.5 GHz SRF cavities

The Continuous Electron Beam Accelerator Facility (CEBAF) was the first large-scale accelerator to employ superconducting radiofrequency (SRF) cavities for continuous-wave operation. Ongoing research and development efforts continue to focus on increasing the intrinsic quality factor (Q 0 ) of these cavities in order to reduce cryogenic losses while maintaining operational gradients. Here, in this work, we report on the application of high-Q 0 surface treatments to single-cell and multicell C100 and C75 style 1.5 GHz niobium cavities used in the CEBAF accelerator. Nitrogen infusion and oxygen alloying via medium-temperature baking were applied under heat-treatment constraints relevant to existing cavity hardware. Both processes yielded substantial improvements in Q 0 at moderate accelerating gradients, achieving values of approximately 2 × 10 10 at 2.07 K and 20 MV/m. The effectiveness of nitrogen infusion at reduced annealing temperatures and the successful extension of oxygen alloying to multicell cavities are demonstrated. These results establish viable pathways for implementing high-Q 0 treatments in CEBAF-compatible cavities and support future integration into cryomodules for reduced operational cryogenic load.

Dhakal, Pashupati [Thomas Jefferson National Accel↗

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗

Current advances on Talbot–Lau x-ray imaging diagnostics for high energy density experiments (invited)

Talbot–Lau x-ray interferometry is a refraction-based diagnostic that can map electron density gradients through phase-contrast methods. The Talbot–Lau x-ray deflectometry (TXD) diagnostics have been deployed in several high energy density experiments. To improve diagnostic performance, a monochromatic TXD was implemented on the Multi-Tera Watt (MTW) laser using 8 keV multilayer mirrors (Δθ/θ = 4.5%-5.6%). Copper foil and wire targets were irradiated at 10 14 –10 15 W/cm 2 . Laser pulse length (~10 to 80 ps) and backlighter target configurations were explored in the context of Moiré fringe contrast and spatial resolution. Foil and wire targets delivered increased contrast <30%. The best spatial resolution (<6 μm) was measured for foils irradiated 80° from the surface. Further TXD diagnostic capability enhancement was achieved through the development of advanced data postprocessing tools. The Talbot Interferometry Analysis (TIA) code enabled x-ray refraction measurements from the MTW monochromatic TXD. Additionally, phase, attenuation, and dark-field maps of an ablating x-pinch load were retrieved through TXD. The images show a dense wire core of ~60 μm diameter surrounded by low-density material of ~40 μm thickness with an outer diameter ratio of ~2.3. Attenuation at 8 keV was measured at ~20% for the dense core and ~10% for the low-density material. Instrumental and experimental limitations for monochromatic TXD diagnostics are presented. Enhanced postprocessing capabilities enabled by TIA are demonstrated in the context of high-intensity laser and pulsed power experimental data analysis. Significant advances in TXD diagnostic capabilities are presented. Furthermore, these results inform future diagnostic technique upgrades that will improve the accuracy of plasma characterization through TXD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of precipitation in gradient Inconel 718 superalloy

Inconel 718 (IN718) superalloy is a well-established material for the application of aircraft engine and power generator due to its excellent performance at high temperatures. The outstanding mechanical properties of IN718 are largely determined by the precipitates. Here, we investigated the precipitation behaviors and the corresponding mechanical properties of the heat-treated gradient IN718 alloy. The microscopy studies reveal that high-density hexagonal η phase with chemistry of Ni 3 (Al,Ti,Nb) formed in the nanograined IN718 alloy after annealing at 700 °C. Additionally, the compressive stress introduced during severe plastic deformation may have prohibited the formation of γ” phase while accelerated the precipitation of γ’ phase, consequently favoring the precipitation of η phase. The corresponding evolution of mechanical behavior of gradient IN718 alloy is investigated.

36 MATERIALS SCIENCE↗

Quantitative characterization of gradient microstructures: A study on friction stir spot processing of pure cobalt

Heterogeneous microstructures in polycrystalline metals can enhance the strength and ductility, outperforming homogeneous structures of similar composition. This study investigates deformed cobalt via friction stir spot processing with varying dwell times to uncover the effects of plastic deformation and heat generation on the formation of morphological, phase, and grain boundary character gradients. A new approach to quantify the morphological gradients in materials, which describes grain morphology in terms of density followed by parametric regression, enables direct quantification of processing depth and gradient sharpness. Results show that longer processing times increase the steepness of morphological gradients and reduce the deformation depth for friction stir spot processing with low plunge depths and high tool rotational speeds. The amount of retained FCC is increased in the shorter processing conditions, primarily due to refined grain size, increased defect content, and reduced heat generation. Crystallographic texture analysis of the HCP phase indicated a dominant B-fiber described by (0001) ∥ shear plane normal in the extreme processing conditions and the formation of a P-fiber, shear direction ∥ ⟨11$\bar2$0⟩ for intermediate dwell times. The texture of the FCC phase for low processing times was a C texture {100}⟨011⟩ where longer processing times were dominated by a [001] fiber texture with a main {110}⟨100⟩ orientation and emergence of a slight [111] fiber in the longest processing condition. The approaches outlined in this work give insight into quantifying gradients and improve the understanding of highly deformed cobalt.

36 MATERIALS SCIENCE↗

Efficient Subset Simulation using Hamiltonian Neural Network enhanced Markov Chain Monte Carlo Methods

The Monte Carlo method delivers an unbiased estimate of the probability of failure. However, the variance of the estimate depends on the number of evaluated samples. This number must be very large for estimations of a low probability of failure. If the evaluation of each sample is computationally expensive, the crude Monte Carlo simulation strategy is impracticable. Therefore, subset simulations are used to reduce the required number of evaluations. Subset simulations require a Markov Chain Monte Carlo sampler, such as the random walk Metropolis-Hastings algorithm. The algorithm, however, struggles with sampling in low-probability regions, especially if they are narrow. As a consequence, advanced Markov Chain Monte Carlo simulations have been developed. In particular, the Hamiltonian Monte Carlo method explores the target distribution rapidly. Driven by the idea of Hamiltonian dynamics, this sampler provides a non-random walk through the target distribution. The incorporation of subset simulation and Hamiltonian Monte Carlo methods has shown promising results for reliability analysis. One downside of the Hamiltonian Monte Carlo method is that gradient evaluations are computationally expensive, especially when dealing with high-dimensional problems and evaluating long trajectories. We show that integrating Hamiltonian neural networks in Hamiltonian Monte Carlo simulations significantly speeds up the sampling task. Furthermore, the enhancement of adaptive trajectory length within the Hamiltonian Monte Carlo results in the efficient proposal of the following states. Based on this recent enhancement, we provide a fast sampling strategy for subset simulations using Hamiltonian neural networks to replace the evaluation of the gradient and significantly speed up the Hamiltonian Monte Carlo simulation.

97 MATHEMATICS AND COMPUTING↗

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

The mapping & avoidance of high order cavity modes

In this paper a technique developed in 2016 is presented for mapping the multidimensional, multi-mode high order mode landscape of a beam driven super conducting RF cavity. The initial use case for this technique was intended to allow longitudinal gradient control of a beam driven super conducting cavity while simultaneously avoiding many dangerous high order modes and excessive component stresses. This paper will focus on the stability and component stress considerations of the 56MHz superconducting cavity, located in the 4 o’clock sector interaction region of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. The 56MHz cavity is intended to work in conjunction with the Stochastic Cooling system to improve the longitudinal confinement in the central 197MHz RF bucket and thus increase the integrated luminosity during heavy ion beam runs at RHIC.

43 PARTICLE ACCELERATORS↗

Supernova Hydrodynamics: A Lab-scale Study of the Blast-driven Instability Using High-speed Diagnostics

A novel experimental approach to study the blast-driven instability at a nondiffuse, gaseous interface with a density gradient is presented. Under Euler similarity, this approach enables study of dissipative-scale hydrodynamics relevant to many astrophysical and laboratory high energy density phenomena in a well-resolved manner. The instability is initiated by passing a Taylor–Sedov blast wave originating from a controlled detonation through a perturbed and stably stratified interface between two gases. The facility and driving blast wave are characterized to obtain repeatable conditions and capture large ensembles of time-resolved Mie scattering imaging that show consistent hydrodynamic development. We analyze the instability evolution between different gas pairs to demonstrate the wide range of development and turbulent behavior that may occur between different supernova layers. The mean evolution of the hydrodynamic instability is compared to a buoyancy–drag model that is frequently used to estimate perturbation growth in supernova mixing research. We propose a time delay to this model in order to reproduce the measured instability behavior and demonstrate model robustness in handling flows driven by a time-varying acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗