Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW VELOCITY - DISTRIBUTION”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma heating and current sheet structure in anti-parallel magnetic reconnection

A theoretical model and an analytic theory of current sheet structure are presented for understanding anti-parallel driven magnetic reconnection in 2-1/2 dimension in collisionless plasmas. The theoretical model provides formulation to compute the current sheet y-profiles by specifying the profiles of electron and ion flow velocities V ex (x,y) and V ix (x,y). The current sheet solutions depend on the plasma density n in , merging magnetic field B 0 , ion velocity v i , and electron velocity ve in the upstream and the S evz =V ez /V dz parameter where V ez is the electron velocity accelerated by the reconnection electric field E z in the electron orbit meandering region, V dz ≃cE y /B x is the $\overrightarrow{E}$ × $\overrightarrow{B}$ drift velocity as electrons enter the orbit meandering region, B x is the merging magnetic field, and E y is the electrostatic electric field. With simplifying assumptions on the y-profiles of V ex and V ix , we have also developed an analytic theory of the current sheet structure. Analytic expressions for the anomalous resistivity, the electrostatic potential drop, and the maximum E y amplitude E max are obtained. The analytic results agree reasonably well with both the particle-in-cell simulation results and the numerical solutions of the theoretical model. The ions energy gain due to the potential drop is $∝B^2_0/n_{in}$. The electron energy gain is $∝(B^2_0 /8πn_{in})S_{evz}$. The $B^2_0/n_{in}$ scaling of the average ion and electron energy gains are consistent with laboratory experiments and space plasma observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microcalorimeter measurement of x-ray spectra from a high-temperature magnetically confined plasma

A NASA-built x-ray microcalorimeter spectrometer has been installed on the MST facility at the Wisconsin Plasma Physics Laboratory and has recorded x-ray photons emitted by impurity ions of aluminum in a majority deuterium plasma. Much of the x-ray microcalorimeter development has been driven by the needs of astrophysics missions, where imaging arrays with few-eV spectral resolution are required. The goal of our project is to adapt these single-photon-counting microcalorimeters for magnetic fusion energy research and demonstrate the value of such measurements for fusion science. Microcalorimeter spectrometers combine the best characteristics of the x-ray instrumentation currently available on fusion devices: high spectral resolution similar to an x-ray crystal spectrometer and the broadband coverage of an x-ray pulse height analysis system. Fusion experiments are increasingly employing high-Z plasma-facing components and require measurement of the concentration of all impurity ion species in the plasma. This diagnostic has the capability to satisfy this need for multi-species impurity ion data and will also contribute to measurements of impurity ion temperature and flow velocity, Zeff, and electron density. Here, we introduce x-ray microcalorimeter detectors and discuss the diagnostic capability for magnetic fusion energy experiments. We describe our experimental setup and spectrometer operation approach at MST, and we present the results from an initial measurement campaign.

Eckart, M. E. (ORCID:0000000338945889)↗
DOE OSTI2022

Vorticity screening by dense layers

We investigate variable density layers in the presence of vortices, by analogy to line charges in the presence of a dielectric slab. We adapt solutions for dielectric layers from the literature to the variable density fluid case, obtaining exact hypergeometric expressions for flow velocities induced everywhere by a vortex in the vicinity of a finite-thickness slab of different densities. The dimensionless Atwood number, which appears naturally in other variable density phenomena, such as the Rayleigh–Taylor instability, reappears here as the natural ordering parameter for the influence of the slab, which acts primarily to screen the influence of vorticity across the slab to a lower effective circulation. The solution takes the form of scaling correction factors that might be applied to computational models or evaluated directly. Extensions, where vorticity is localized within the dense layer itself and amplified outside of it, are considered, as well as some speculative applications such as inertial confinement fusion capsule designs.

42 ENGINEERING↗
DOE OSTI2023

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Axial confinement of wire array Z-pinch precursor plasmas by a pulsed magnetic mirror field

A novel wire-array Z-pinch platform has been developed to study the effect of a pulsed magnetic mirror field on a collisional high energy density plasma. The mirror is driven in series with the Z-pinch target and, therefore, scales with driver current. Axial confinement is observed based on visible imaging and axial flow velocity measurements. The presence of axial compressing J xB force is determined indirectly based on Thomson scattering and interferometry measurements and corroborated by three-dimensional extended-MHD simulations. Compared to non-magnetized wire array Z-pinch, a modified pulsed mirror configuration is observed to increase precursor plasma thermal energy density by about 30%. If optimized, such a configuration could potentially improve magnetized liner inertial fusion performance by reducing axial plasma end loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The next-generation magnetic recoil spectrometer (MRSnext) on OMEGA and NIF for diagnosing ion temperature, yield, areal density, and alpha heating

The next-generation magnetic recoil spectrometer (MRSnext) is being designed to replace the current MRS at the National Ignition Facility and OMEGA for measurements of the neutron spectrum from an inertial confinement fusion implosion. The MRSnext will provide a far-superior performance and faster data turnaround than the current MRS systems, i.e., a 2× and 6× improvement in energy resolution at the NIF and OMEGA, respectively, and 20× improvement in data turnaround time. The substantially improved performance of the MRSnext is enabled by using electromagnets that provide a short focal plane (12–16 cm) and unprecedented flexibility for a wide range of applications. In addition to being able to measure neutron yield, apparent ion temperature, areal density, and plasma-flow velocity over a wide range of yields, the NIF MRSnext will be able to directly, uniquely assess the alpha heating of the fuel ions through measurements of the alpha knock-on tail in the neutron spectrum. The goal is to implement a radiation-hard electronic detection system capable of providing rapid data acquisition and analysis. The development of the MRSnext will also set the foundation for the more advanced, time-resolving MRSt and serve as a testbed for its implementation on the NIF.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A Q-band frequency tunable Doppler backscattering (DBS) system for pedestal and scrape-off layer density fluctuation and flow measurements in the DIII-D tokamak

We present the design and laboratory tests for a new Q-band frequency tunable Doppler backscattering (DBS) system suitable for probing poloidal wavenumber kñ = 6–8 cm−1 density fluctuations and their flow velocities in the pedestal and scape-off layer (SOL) of the DIII-D tokamak. This system will provide new measurements in the increasingly important and under-diagnosed far pedestal and SOL plasma regions. These results are important for experimental transport studies and necessary for the validation of transport models, both of which are important to fusion energy research. The use of a single tunable frequency reduces the complexity and potential failure points as compared to a multichannel system. This new system utilizes a 33–50 GHz tunable source and will be integrated into the current V-band DBS in DIII-D using a broadband Q- and V-band multiplexer. A full-scale mockup of the quasi-optical system was used to test and optimize the performance. These tests include beam profile measurements at different distances (and angles) from a paraboloidal focusing and steering mirror. The measurements cover the full frequency range 33–75 GHz of the integrated/combined Q–V band DBS system and target a large radial coverage of the low-field side of the plasma from ρ = 1.1 to ρ = 0.5, where ρ is the normalized flux surface radial coordinate.

Instruments & Instrumentation↗
DOE OSTI2024

X-ray sources for in situ wavelength calibration of x-ray imaging crystal spectrometers

X-ray sources for a range of wavelengths are being considered for in situ calibration of X-ray Imaging Crystal Spectrometers (XICSs) and for monitoring line shifts due to changes in the crystal temperature, which can vary during experimental operation over a day [A. Ince-Cushman et al., Rev. Sci. Instrum. 79, 10E302 (2008), L. Delgado-Aparicio et al., Plasma Phys. Control. Fusion 55, 125011 (2013)]. Such crystal temperature dependent shifts, if not accounted for, could be erroneously interpreted as Doppler shifts leading to errors in plasma flow-velocity measurements. The x-ray sources encompass characteristic x-ray lines falling within the wavelength range of 0.9–4.0 Å, relevant for the XICSs on present and future fusion devices. Several technological challenges associated with the development of x-ray sources for in situ calibration are identified and are being addressed in the design of multiple x-ray tubes, which will be installed inside the spectrometer housing of the XICS for the JT-60SA tokamak. These x-ray sources will be especially useful for in situ calibration between plasma discharges. Here, in this paper, laboratory experiments are described that were conducted with a Cu x-ray source, a heated quartz (102) crystal, and a pixelated Pilatus detector to measure the temperature dependent shifts of the Cu Kα 1 and Kα 2 lines at 1.5405 and 1.5443 Å, respectively, and to evaluate the 2d-lattice constant for the Bragg reflecting crystal planes as a function of temperature, which, in the case of in situ wavelength calibration, would have to be used for numerical analysis of the x-ray spectra from the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental investigation of laser ablation of stone polycrystalline targets

We report the results of an experimental investigation of ablation of stone polycrystalline targets of complex multicomponent composition, which imitate the substance of asteroids. The targets were irradiated by nanosecond pulses of a neodymium laser at an energy density Φ{sub L} of up to 5 × 10{sup 4} J cm{sup −2}. The experiments demonstrated the existence of two ablation regimes, with the boundary between them lying at Φ{sub L} ≈ 4000 J cm{sup −2}. The regime change is characterised by a change in the form of the dependence of the surface mass density of removed target material on the laser energy density and by the appearance of a minimum in the dependence of specific energy of destruction on Φ{sub L}. This is supposedly related to the passage from a one-dimensional plasma plume expansion to the three-dimensional one and the corresponding decrease in the efficiency of energy transfer from the laser beam to the target due to a lowering of laser-produced plasma density. Our experiments also showed the existence of a maximum in impulse coupling coefficient C{sub m} as a function of laser energy density (C{sub m} ≈ 6.3 × 10{sup −5} N W{sup −1} for Φ{sub L} = Φ{sub opt} ≈ 25 J cm{sup −2}). Maxima were also recorded in the dependences of the ablation efficiency and average ablation flow velocity on Φ{sub L}. For Φ{sub L} > Φ{sub opt}, the decrease in the function C{sub m}(Φ{sub L}) turns out to be much steeper than for metals and polymer materials. The difference is supposedly due to the lower strength and lower plasticity of the polycrystalline stone targets. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Time-resolved turbulent dynamo in a laser plasma

Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas (Pm < 1). However, the same framework proposes that the fluctuation dynamo should operate differently when Pm ≳ 1, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports an experiment that creates a laboratory Pm ≳ 1 plasma dynamo. We provide a time-resolved characterization of the plasma’s evolution, measuring temperatures, densities, flow velocities, and magnetic fields, which allows us to explore various stages of the fluctuation dynamo’s operation on seed magnetic fields generated by the action of the Biermann-battery mechanism during the initial drive-laser target interaction. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude and saturate dynamically. It is shown that the initial growth of these fields occurs at a much greater rate than the turnover rate of the driving-scale stochastic motions. Here, our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized magnetohydrodynamics (MHD) simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2020

Solitary zonal structures in subcritical drift waves: a minimum model

Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same 'equation of state', which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 &#x00d7; 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E &#x00d7; B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E &#x00d7; B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗