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

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗
DOE OSTI2023

Extensive 3D mapping of dislocation structures in bulk aluminum

Abstract Thermomechanical processing such as annealing is one of the main methods to tailor the mechanical properties of materials, however, much is unknown about the reorganization of dislocation structures deep inside macroscopic crystals that give rise to those changes. Here, we demonstrate the self-organization of dislocation structures upon high-temperature annealing in a mm-sized single crystal of aluminum. We map a large embedded 3D volume ( $$100\times 300\times 300$$ 100 × 300 × 300 $$\upmu$$ μ m $$^3$$ 3 ) of dislocation structures using dark field X-ray microscopy (DFXM), a diffraction-based imaging technique. Over the wide field of view, DFXM’s high angular resolution allows us to identify subgrains, separated by dislocation boundaries, which we identify and characterize down to the single-dislocation level using computer-vision methods. We demonstrate how even after long annealing times at high temperatures, the remaining low density of dislocations still pack into well-defined, straight dislocation boundaries (DBs) that lie on specific crystallographic planes. In contrast to conventional grain growth models, our results show that the dihedral angles at the triple junctions are not the predicted 120 $$^{\circ }$$ ∘ , suggesting additional complexities in the boundary stabilization mechanisms. Mapping the local misorientation and lattice strain around these boundaries shows that the observed strain is shear, imparting an average misorientation around the DB of $$\approx$$ ≈ 0.003 to 0.006 $$^{\circ }$$ ∘ .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗
DOE OSTI2025

Superheating gold beyond the predicted entropy catastrophe threshold

In their landmark study, Fecht and Johnson unveiled a phenomenon that they termed the ‘entropy catastrophe’, a critical point where the entropy of superheated crystals equates to that of their liquid counterparts. This point marks the uppermost stability boundary for solids at temperatures typically around three times their melting point. Despite the theoretical prediction of this ultimate stability threshold, its practical exploration has been prevented by numerous intermediate destabilizing events, colloquially known as a hierarchy of catastrophes, which occur at far lower temperatures. Here we experimentally test this limit under ultrafast heating conditions, directly tracking the lattice temperature by using high-resolution inelastic X-ray scattering. Our gold samples are heated to temperatures over 14 times their melting point while retaining their crystalline structure, far surpassing the predicted threshold and suggesting a substantially higher or potentially no limit for superheating. We point to the inability of our samples to expand on these very short timescales as an important difference from previous estimates. These observations provide insights into the dynamics of melting under extreme conditions.

Laser-produced plasmas↗
DOE OSTI2025

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗
DOE OSTI2026

The impact of kinetic and global effects on ideal ballooning 2nd stable pedestals of conventional and low aspect-ratio tokamaks

The EPED model (Snyder et al 2011 Nucl. Fusion 51 103016) has had success in describing the pedestal structure of Type-I ELM and QH-mode plasmas in tokamaks, by combining kinetic ballooning mode (KBM) and peeling–ballooning (PB) constraints. Within EPED, the KBM constraint is often approximated by a technique using calculated ideal ballooning mode (IBM) thresholds, fit to a functional form designed to capture non-local effects. It has been noted that quantitative differences between local ideal MHD and gyro-kinetic (GK) ballooning stability can be larger at low aspect ratio. KBM critical pedestals are consistent with observations in initial studies on conventional and spherical tokamaks. In this work, the application of a reduced model for the calculation of the kinetic ballooning stability boundary is presented based on a novel and newly developed gyro-fluid system (GFS) code (Staebler et al 2023 Phys. Plasmas 30 102501). GFS is observed to capture KBMs in DIII-D as well as the NSTX pedestals, opening a route to integrating this model into EPED. Finally, high but finite $n$ global ballooning modes are observed to limit the access to the local 2nd stability and thus provide a transport mechanism that constrains the pedestal evolution with $β_{p,\textrm{ped}}$. The high $n$ global ballooning stability is approximated by its ideal MHD analog using ELITE. It is shown that nearly-local high $n$ modes with $k_yρ_s$ ~ $0.25 - 0.5$ can provide a proxy for the critical $β_{p,\textrm{ped}}$ when a 2nd stable access exists on DIII-D plasmas. The use of GFS and ELITE scaling in EPED provides improved agreement to EPED1 in an initial comparison with DIII-D pedestal data.

Anastopoulos Tzanis, Michail Savvas [Oak Ridge Nat↗
DOE OSTI2024

Prediction of pellet mass thresholds for ELM triggering in low-collisionality, ITER-like discharges

Abstract In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mechanical and Tribological Performance of Additively Manufactured Nanocrystalline Aluminum via Cryomilling and Cold Spray

In this study, nanocrystalline (NC) aluminum (Al) and magnesium (Mg)-doped Al bulk components were fabricated using a hybrid manufacturing process that combines cryomilling and high-pressure cold spray (HPCS) additive deposition techniques. Yttria-stabilized zirconia (YSZ) was also added during the HPCS process to improve deposition efficiency and build-up thickness via peening. The evolution of morphology, crystallite size, and elemental composition of both cryomilled powders and cold-sprayed (CS’ed) components was systematically characterized using X-ray Diffraction (XRD), Scanning Electron Microscopy (SEM), and Transmission Electron Microscopy (TEM). Mechanical characterization was performed using Vickers microhardness and uniaxial tensile testing, while the tribological behavior was assessed using sliding wear tests under dry/lubricated conditions. XRD analysis revealed that increased cryomilling duration led to significant crystallite refinement, which directly correlated with enhanced hardness and strength. This mechanical strengthening was accompanied by an increase in coefficient of friction (COF) and lower wear rates. The results also showed that the Mg-doped Al exhibited superior hardness, tensile strength, and tribological performance compared to pure Al. The study further explores the underlying mechanisms responsible for these enhancements, highlighting the potential of solute-assisted grain boundary stabilization in tailoring high-performance NC Al alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In-situ radiation response of additively manufactured modified Inconel 718 alloys

Here in this study, a novel alloy of modified Inconel 718 produced by laser powder bed fusion is studied before and after in-situ Kr irradiation up to 3 dpa at 200 and 450 °C. Before irradiation, the microstructure consists of dislocation cells having a misorientation angle less than 5° and with an average size of ~500 nm. There are also second phase particles of MC type carbides, Laves phase and oxides such as Y-O, Y-(Ti)-Al-O. While the microstructure consists of stacking fault tetrahedra, faulted and perfect loops after irradiation at 200 °C, dislocation loops are the primary defects at 450 °C. With increasing dose, the size of the defects remains similar at 200 °C while it increases at 450 °C. This has been attributed to the existence of vacancy type defects at 200 °C and the different defect transport mechanisms at different temperatures. Moreover, matrix and second phase particle compositions seem to be similar after irradiation. The sink strengths of the structures have been calculated and superior radiation resistance of this alloy has been attributed to the existence of fine cell boundaries stabilized by the second phase particles produced by additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Tandem Mass-Selective Cryogenic Digital Ion Traps for Enhanced Cluster Formation

We present the implementation of tandem mass-selective cryogenic ion traps, designed to enhance the range of ion processing capabilities that can be performed prior to spectroscopic interrogation. We show that both the formation of ion clusters and mass filtering steps can be combined in a single cryogenic linear quadrupole ion trap driven by RF square waves. Mass filtering and mass isolation can be achieved by manipulation of the RF frequency and duty cycle. Very importantly, this scheme circumvents the need for high-amplitude RF voltages that can be incompatible with typical cryogenic ion processing conditions. In addition, proper adjustment of the stability boundaries during the clustering process allows for the preferential formation of a specific cluster size rather than a broad distribution of sizes. Lastly, we show that a specific cluster size can be formed, mass-selected, and then transferred to another ion trap for a second completely separate ion processing step. Furthermore, the instrumentation and modular design developed here expand the scope of ionic species and clusters that can be accessed by processing electrosprayed ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

On the new universality class in structurally disordered n -vector model with long-range interactions

We study a stability boundary of a region where nontrivial critical behavior of an n-vector model with long-range power-law decaying interactions is induced by the presence of a structural disorder (e.g., weak quenched dilution). This boundary is given by the marginal dimension of the order parameter nc dependent on space dimension, d, and a control parameter of the interaction decay, σ, below which the model belongs to the new dilution-induced universality class. Exploiting the Harris criterion and recent field theoretical renormalization group results for the pure model with long-range interactions, we get nc as a three loop ɛ = 2σ – d-expansion. We provide numerical values for nc applying series resummation methods. Our results show that not only the Ising systems (n = 1) can belong to the new disorder-induced long-range universality class at d = 2 and 3.

Physics↗
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↗