Search NASA⌕ Search

SEARCH · Search NASA

Results for “core-edge”

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

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Core performance predictions in projected SPARC first-campaign plasmas with nonlinear CGYRO

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗

Overview of results from the 2023 DIII-D negative triangularity campaign

Negative triangularity (NT) is a potentially transformative configuration for tokamak-based fusion energy with its high-performance core, edge localized mode (ELM)-free edge, and low-field-side divertors that could readily scale to an integrated reactor solution. Previous NT work on the TCV and DIII-D tokamaks motivated the installation of graphite-tile armor on the low-field-side lower outer wall of DIII-D. A dedicated multiple-week experimental campaign was conducted to qualify the NT scenario for future reactors. During the DIII-D NT campaign, high confinement ($H_\mathrm{98y,2}\gtrsim$ 1), high current ($q_\mathrm{95} \lt$ 3), and high normalized pressure plasmas ($\beta_\mathrm{N} \gt$ 2.5) were simultaneously attained in strongly NT-shaped discharges with average triangularity $\delta_\mathrm{avg}$ = −0.5 that were stably controlled. Experiments covered a wide range of DIII-D operational space (plasma current, toroidal field, electron density and pressure) and did not trigger an ELM in a single discharge as long as sufficiently strong NT was maintained; in contrast, to other high-performance ELM-suppression scenarios that have narrower operating windows. These strong NT plasmas had a lower outer divertor X-point shape and maintained a non-ELMing edge with an electron temperature pedestal, exceeding that of typical L-mode plasmas. Also, the following was achieved during the campaign: high normalized density ($n_\mathrm{e}$/$n_\mathrm{GW}$ of at least 1.7), particle confinement comparable to energy confinement with $Z_\mathrm{eff}\sim2$ , a detached divertor without impurity seeding, and a mantle radiation scenario using extrinsic impurities. These results are promising for a NT fusion pilot plant but further questions on confinement extrapolation and core-edge integration remain, which motivate future NT studies on DIII-D and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor characterization and access to dissipative divertor conditions in negative triangularity discharges in DIII-D

Reduced radial transport, short midplane-to-target parallel connection lengths, and a strong effect of cross-field drifts were responsible for the high densities required for detachment in strong negative triangularity (NT) configurations in DIII-D. Dissipative divertor conditions were achieved in NT discharges at different triangularity, injected power, plasma current, and toroidal field direction. Differences between negative and positive triangularity (PT) discharges are analyzed in this paper to understand the requirements for access to detached divertor conditions: power balance, geometry, radial transport and effect of cross-field particle drifts. Parametric dependencies of access to detachment on plasma current and power flowing into the scrape-off layer remained similar in negative and PT and impurity seeding was observed to reduce the density needed to detach by up to 30% at the expense of core impurity dilution. The impact of triangularity on core-edge integration was tested varying bottom triangularity at fixed top triangularity. The high density needed to detach was not intrinsic to the NT edge as shapes with positive lower triangularity and negative upper triangularity were able to detach at lower upstream densities while maintaining an ELM-free NT edge. Confinement degradation at deeper detachment levels was however observed in all NT shapes, often associated with radiation instabilities.

detachment↗

Dependence of the boundary heat flux width on core and edge profiles in Alcator C-Mod

This work presents new evidence that the heat flux width, λ q , in the Alcator C-Mod tokamak scales with the edge electron pressure, as observed in the ASDEX Upgrade (AUG) tokamak, but the scaling with volume-averaged pressure, $\bar{p}$, from the plasma stored energy, is a better predictor of λ q in Alcator C-Mod than the edge electron pressure. These previous studies, which find that λ q decreases with increasing plasma pressure, imply that a high performance core at high pressure will lead to challenging heat and particle exhaust due to very small λ q . This concern has led to our significant enlargement of the C-Mod database with the electron density, temperature, and pressure profile data from the Thomson scattering and electron cyclotron emission diagnostics. Using the C-Mod database augmented with new profile data, we find that λ q decreases with increasing edge electron pressure as ${\lambda }_{q}\propto {p}_{\mathrm{e},95}^{-0.26}$, similar to results from AUG, and showing the strength of cross-machine comparisons. We also find that ${\lambda }_{q}\propto {p}_{\mathrm{e},\mathrm{c}\mathrm{o}\mathrm{r}\mathrm{e}}^{-0.56}$, consistent with the original finding from C-Mod that the heat flux width scales as ${\bar{p}}^{-0.48}$. The scalings of λ q with separatrix pressure and gradient scale length are found to match the AUG results qualitatively. The C-Mod scalings with edge plasma quantities have more scatter than the $\bar{p}$ scaling, and, importantly, show different trends for H-modes relative to L- and I-mode. Furthermore, investigating the source of this discrepancy presents an opportunity for further study that may improve our ability to predict the heat flux width in different confinement scenarios in the pursuit of optimizing core-edge performance in future reactors.

Alcator C-Mod↗

Feasibility of raised inner strike point equilibria scenario in ITER for detritiation from beryllium co-deposits

Abstract In ITER, tritium retention primarily occurs through co-deposition with beryllium. To avoid exceeding the strict tritium inventory limit, efficient tritium recovery techniques are essential. Baking is the ITER baseline for tritium recovery, but its effectiveness in removing tritium from thick beryllium layers is limited. A raised strike point scenario is considered an alternative method for removing tritium from the ITER inner vertical divertor target by heating components via plasma flux. This paper presents SOLPS-ITER code simulations conducted under various conditions, assessing the divertor performance and tritium outgassing of the raised strike point scenario. As the strike point is raised, recycled neutrals are not efficiently baffled by the dome and scrape-off layer, significantly changing the neutral trajectory and ionization source distribution. This improves detachment accessibility but worsens core-edge compatibility compared to the baseline scenario. However, in the partially detached condition, the impact of raising the strike point, perpendicular transport, and q 95 on target heat flux is not significant, as it primarily scales with the input power. Target heat flux is translated to target surface temperature using a simplified heat transfer model that considers the 3D target monoblock geometry and active cooling condition, excluding Be layer thermal properties. For partially detached divertor conditions, the bulk tungsten monoblock surface temperature remains below the baking temperature, which is insufficient for efficient tritium outgassing under the actively cooled ITER divertor condition. However, considering the potential thermal contact resistance between the beryllium and tungsten layers, which may significantly impact temperature distribution, the temperature of the beryllium layer can be raised to a level sufficient for efficient tritium outgassing. Therefore, the raised strike point scenario can be considered as an alternative in-vessel tritium removal technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

Abstract This paper reports on the refinement (building on (Hatch D.R. et al 2022 Phys. Plasmas 29 062501)) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model’s applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature profile exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high η e = L n e / L T e scenarios but not for low η e scenarios in which microtearing modes (MTMs) have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport-driven toroidal rotation with general viscosity profile

Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simultaneous ELM suppression and divertor detachment via synergistic boron powder and neon injection in EAST

A novel approach for simultaneous power exhaust and edge-localized mode (ELM) control is presented in the Experimental Advanced Superconducting Tokamak discharges, which utilize an ITER-like tungsten divertor. Real-time injection of boron (B) powder and neon (Ne) gas overcomes their limitations encountered when used separately. Pure Ne seeding leads to a narrow operational window constrained by core impurity accumulation and H-mode to L-mode back transitions, while pure solid B injection (SBI) is insufficient for effective divertor cooling. In comparison, their combined use achieves a stable, stationary, ELM-suppressed H-mode with adequate power exhaust. This synergistic scenario features partial energy detachment at the outer divertor while maintaining good plasma confinement (H 98 ∼ 1) with minimal degradation. Two key features of this scenario are: (1) the SBI triggers a persistent Edge Harmonic Mode (EHM), which provides a crucial continuous particle transport channel, preventing Ne and tungsten/molybdenum accumulation without flushing out by ELM, and (2) the B + Ne mixture allows for active optimization of the radiated power profile. Core radiation can be reduced by substituting a portion of the Ne with B, leveraging their complementary non-coronal equilibrium radiation efficiencies. This combined B + Ne injection scheme presents a promising pathway toward integrated core-edge scenarios, offering the potential to minimize total impurity throughput while leveraging an actuator (powder injection) already being considered for ITER.

Boron↗

Description of FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

The urgency to deliver fusion power is growing now more than ever, with increasing pressure for both public programs and private companies to meet milestones timelines and overcome significant remaining technical challenges to ensure growth of a nascent fusion industry in time to meet rapidly growing clean energy demands. With incredible advancements in computation and years of investment in fusion model development and validation, integrated modeling is poised to fill a key role in accelerating the timeline to a fusion pilot plant (FPP). Future fusion pilot plants will operate in regimes far beyond current experience, and device design will rely on physics-based prediction and extrapolation. Many concepts will also rely on simulation to assess safety (shielding, tritium management, materials activation and lifetimes), economics and scalability before the decision to build. Importantly, integrated simulation can be used to reveal and solve the complexities of system integration that may otherwise not be apparent in physical components or models developed in isolation. New experimental test facilities that produce relevant conditions to validate and resolve key technical challenges for various subsystems (materials, blankets, fuel cycle, etc.) have been repeatedly called for by the fusion community but are not yet realized. Integrated modeling has an important role in identifying realistic load conditions (thermal, electromagnetic, plasma, neutron and photon loads, etc.) and defining the components and experiments for these test facilities in order to ensure meaningful validation that sufficiently reduces modeling uncertainties and technical risk for the full integrated reactor. The Fusion REactor Design and Assessment (FREDA) SciDAC project is building a component-based integrated modeling framework & data structure to enable self-consistent, multi-fidelity, iterative optimization workflows for the fusion reactor design process. FREDA aims to shorten the time to viable designs by providing a set of flexible workflows to support the various stages of the design process using an integrated model hierarchy, ranging from the simple analytic descriptions to the highest fidelity, theory-based plasma and engineering modeling developed by the fusion and fission communities. These tools are expected to be needed for timely support of FPP design in the milestone program and in the FIRE collaboratives. The plasma simulation backbone of FREDA is IPS-FASTRAN with newly developed coupled Core-Edge Pedestal-SOL (CESOL) workflows, which is being extended to the far-SOL region up to the plasma facing components. FREDA incorporates the FERMI engineering modeling suite and will enable self-consistent evaluation of the thermal shields, limiters, blanket, magnets, and other surrounding structures with predictions of temperatures, erosion, dpa, activation, tritium generation and transport, creep, corrosion, material degradation, etc. Parametric generation of 3D CAD enables rapid iteration of component geometry in response to plasma and loading specifications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integration of high confinement, high poloidal beta plasma with dual radiated power and detachment controls for divertor protection and ELM suppression

Divertor detachment without serious core confinement quality loss in DIII-D’s high poloidal β scenario has been combined with impurity-induced ELM mitigation without disruption. Use of Ne previously granted access to a detached, non-ELMing regime that retained high confinement quality due to stimulation of Internal Transport Barrier (ITB) growth, but suffered from on-going core fuel dilution and high disruptivity. Excess Ne accumulation in the core plasma has now been avoided by feeding back core radiated power ($P_{rad}$) measurements to control Ne seeding, rather than using attachment fraction ($A_{frac}$) control with Ne; this also reduces disruptivity. At the same time, N 2 seeding is used in a feedback loop with $A_{frac}$ measurements, which previously posed low disruption risk. In this way, the effect of Ne in the core is managed while avoiding excess seeding, and N 2 acts to correct for any excess heat exhaust that might interfere with detachment. The average Ne flow rate was 38% of what was used in pure Ne $A_{frac}$ control, plus average N 2 flow that was 43% of pure N 2 $A_{frac}$ control, even while meeting an even deeper detachment target. Meanwhile, a steep pressure gradient in the core plasma reduces sensitivity to impurity-induced degradation of the pedestal and reasonable confinement quality was maintained despite operational challenges that blocked formation of an ITB in these experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing the negative triangularity reactor core operating space with integrated modeling

Negative triangularity (NT) has received renewed interest as a fusion reactor regime due to its beneficial power-handling properties, including low scrape-off layer power and a larger divertor wetted area that facilitates simple divertor integration. NT experiments have also demonstrated core performance on par with positive triangularity (PT) high confinement mode (H-mode) without edge-localized modes (ELMs), encouraging further study of an NT reactor core. In this work, we use integrated modeling to scope the operating space around two NT reactor strategies. The first is the high-field, compact fusion pilot plant concept Modular, Adjustable, NT ARC (MANTA) and the second is a low field, high aspect ratio concept based on work by Medvedev et al. By integrating equilibrium, core transport, and edge ballooning instability models, we establish a range of operating points with less than 50 MW scrape-off layer power and fusion power comparable to PT H-mode reactor concepts. Heating and seeded impurities are leveraged to accomplish the same fusion performance and scrape-off layer exhaust power for various pressure edge boundary conditions. Scans over these pressure edge conditions accommodate any current uncertainty of the properties of the NT edge and show that the performance of an NT reactor will be extremely dependent on the edge pressure. The high-field case is found to enable lower scrape-off layer power because it is capable of reaching high fusion powers at a relatively compact size, which allows increased separatrix density without exceeding the Greenwald density limit. Adjustments in NT shaping exhibit small changes in fusion power, with an increase in fusion power density seen at weaker NT. Infinite-n ballooning instability models indicate that an NT reactor core can reach fusion powers comparable to leading PT H-mode reactor concepts while remaining ballooning-stable. Seeded krypton is leveraged to further lower scrape-off layer power since NT does not have a requirement to remain in H-mode while still maintaining high confinement. We contextualize the NT reactor operating space by comparing to popular PT H-mode reactor concepts, and find that NT exhibits competitive ELM-free performance with these concepts for a variety of edge conditions while maintaining relatively low scrape-off layer power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fusion pilot plant performance and the role of a sustained high power density tokamak

Here, recent U.S. fusion development strategy reports all recommend that the U.S. should pursue innovative science and technology to enable construction of a fusion pilot plant (FPP) that produces net electricity from fusion at low capital cost. Compact tokamaks have been proposed as a means of potentially reducing the capital cost of a FPP. However, compact steady-state tokamak FPPs face the challenge of integrating a high fraction of self-driven current with high core confinement, plasma pressure, and high divertor parallel heat flux. This integration is sufficiently challenging that a dedicated sustained-high-power-density (SHPD) tokamak facility is proposed by the U.S. community as the optimal way to close this integration gap. Performance projections for the steady-state tokamak FPP regime are presented and a preliminary SHPD device with substantial flexibility in lower aspect ratio (A = 2–2.5), shaping, and divertor configuration to narrow gaps to an FPP is described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height and width by taking the separatrix temperature and density calculated from SOLPS as inputs to EPED. An empirical proxy function based on data from a subset of these discharges is then used to infer the ratio of density at the pedestal top to that at the separatrix based on the electron temperature at the divertor target. This allows the density at the top of the pedestal to be provided to EPED as well, so that knowledge of this quantity is no longer required a priori, while SOLPS radial transport coefficients are needed instead. The gas puff rate provided to SOLPS is then a primary model input for this forward model coupling (along with radial transport coefficients, which are kept fixed through each scan in this study). Comparisons with measurements from DIII-D density scans through a range of target conditions find significant pedestal pressure degradation as the density is increased towards detachment in three different divertor configurations, consistent with the experimental data in these ballooning-limited pedestal regimes.

core-edge integration↗