Search NASA⌕ Search

SEARCH · Search NASA

Results for “SOLPS”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Super-X and conventional divertor configurations in MAST-U ohmic L-mode; a comparison facilitated by interpretative modelling

Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,\textrm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16 $\,\pm\,7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e = 5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream n e and T e profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,\textrm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation of the impact of particle recycling on the plasma in MPS-LD device based on the BOUT++ LPD module

A linear plasma device (LPD) module has been developed under the BOUT++ framework to simulate plasma transport in the MPS-LD. However, previously, the LPD module used a simplistic neutral particle model that only includes particle density and velocity, which prevents the full understanding of the plasma-neutrals interactions. In this work, we further optimize the neutral model by using a more complete neutral fluid model containing the continuity equation, momentum equation, and energy equation. The reactions such as charge exchange, excitation, and radiation collisions are included. Since the neutral particle source is mainly provided by particle recycling from the target, a particle recycling model is employed, which includes both fast reflection and slow thermal release. The upgraded LPD module is applied to simulate the argon (Ar) discharge experiment of MPS-LD, and the benchmark against experiment measurement and SOLPS-ITER simulation results are presented. Good agreements are obtained, showing the validation of the upgraded module. After that, the impact of particle recycling on Ar plasma is investigated. It is found that a higher recycling coefficient (R) promotes the achievement of high-density plasma at the target. In conclusion, the recycled Ar atoms change target plasma pressure as well as plasma-neutral collisions, which both contribute to plasma momentum loss, thus promoting the rollover of ion flux to the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermal Response of a Lithium Vapor Divertor to Cyclical Operation

The lithium vapor divertor concept is being developed as a method to achieve detached divertor conditions in a tokamak while minimizing impurity radiation losses from the core plasma. SOLPS-ITER modeling has previously been used to identify some of the geometric constraints and required lithium evaporation rate of a lithium vapor divertor in a medium-sized tokamak during steady-state operation. Here an updated conceptual design based on these operating requirements is introduced and the thermal response of the system is modeled during cyclical operation, consistent with operation in a short-pulse tokamak. Controllability of the temperature of the lithium capillary porous system (CPS) is achieved by adopting a design where there is no line-of-sight for radiation from the plasma to reach the heated CPS surface. Operational strategies to minimize the amount of lithium evaporated between plasma discharges while achieving steady evaporation rates during plasma discharges are discussed and modeled here. The optimal feedforward control strategy demonstrated in this work is to ramp up the temperature of the evaporator as quickly as possible immediately before a plasma discharge and then reduce the heating to match the desired steady-state net evaporation rate just before the plasma discharge begins, allowing the thermal inertia of the system to stabilize the evaporation rate during the first second of the plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Erosion resistance test of SiC mirror sample for ITER divertor VUV spectrometer

A series of tests have been performed to validate the resistance of the Silicon Carbide (SiC) mirror as the first mirror material of ITER VUV spectrometers to all ITER environmental conditions. Here we focused on the erosion (and deposition) of the SiC mirror sample caused by high-energy neutral particles. The flux of neutral particles reaching the first mirror was calculated using the ZEMAX software with a simplified entrance duct model. In the calculation, the particle flux reaching the first wall is necessary and the previously reported values derived from the SOLPS calculations were used. Based on this estimated particle flux, the erosion resistance tests were performed to check the erosion effect due to the high-energy neutral particles on the first mirror of ITER divertor VUV spectrometer. In the experiment, erosion was induced by exposing the SiC mirror sample to hydrogen and deuterium plasmas (and ions). The target fluence of incident ions in experiment is based on the estimation of the flux of neutral particles in ITER. The surface shape, composition, erosion depth, and surface roughness were measured to check the damage of the mirror surface after erosion test. Based on simulations and erosion resistance tests, it was concluded that the SiC mirror can be used as the first mirror of ITER divertor VUV spectrometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling the effect of nitrogen recycling on the erosion and leakage of tungsten impurities from the SAS-VW divertor in DIII-D during nitrogen gas injection

The new SAS-VW divertor in DIII-D has tungsten-coated components to enable the study of tungsten erosion and leakage from a closed, slot-like divertor. A proposed method of actively managing the tungsten impurities is to inject low-Z impurities, such as nitrogen, into the Scrape-off-Layer (SOL) to modify the conditions in the plasma boundary and in turn manipulate both the erosion and subsequent transport of tungsten. Nitrogen injection from the SOL crown has been modeled using SOLPS-ITER to calculate the background plasma including the intrinsic carbon and the injected nitrogen impurities, and subsequently DIVIMP has been used to calculate the tungsten erosion and transport on top of the background plasma solution. This workflow has been used to model scenarios at a variety of nitrogen injection rates with different assumptions about the nitrogen recycling at the target. In the scenario modeled here, an optimal nitrogen injection rate around 3-4 x 10 20 N/s is found to reduce the amount of tungsten reaching the core by a factor of about 2.4. However, when the nitrogen recycling rate at the divertor targets is high, the nitrogen redistributes within the slot leading to increased tungsten sputtering, and the range of injection rates resulting in tungsten mitigation becomes narrower.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Detachment scalings derived from 1D scrape-off-layer simulations

Fusion power plants will require detachment to mitigate sputtering and keep divertor heat fluxes at tolerable levels. Controlling detachment on these devices may require the use of real-time scrape-off-layer modeling to complement the limited set of available diagnostics. In this work, we use the configurable Hermes-3 edge modeling framework to perform time-dependent, fixed-fraction-impurity 1D detachment simulations. Although currently far from real-time, these simulations are used to investigate time-dependent effects and the minimum physics set required for control-relevant modeling. We show that these simulations reproduce the expected rollover of the target ion flux — a typical characteristic of detachment onset. We also perform scans of the input heat flux and impurity concentration and show that the steady-state results closely match the scalings predicted by the 0D time-independent Lengyel–Goedheer model. This allows us to indirectly compare to SOLPS simulations, which find a similar scaling but a lower value for the impurity concentration required for detachment for given upstream conditions. We use this result to suggest a series of improvements for the Hermes simulations, and finally show simulations demonstrating the impact of time-dependence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Particle control via cryopumping and its impact on the edge plasma profiles of Alcator C-Mod

At the high n e proposed for high-field fusion reactors, it is uncertain whether ionization, as opposed to plasma transport, will be most influential in determining n e at the pedestal and separatrix. A database of Alcator C-Mod discharges is analyzed to evaluate the impact of source modification via cryopumping. The database contains similarly-shaped H-modes at fixed I P = 0.8 MA and B t = 5.4 T, spanning a large range in P net and ionization. Measurements from an edge Thomson scattering system are combined with those from a midplane-viewing Ly α camera to evaluate changes to n e and T e in response to changes in ionization rates, S ion ∙ $n^{sep}_e$ and $T^{ped}_e$ are found to be most sensitive to changes to $S^{sep}_{ion}$, as opposed to $n^{ped}_e$ and $T^{sep}_e$. Dimensionless quantities, namely α MHD and v*, are found to regulate attainable pedestal values. Select discharges at different values of P net and in different pumping configurations are analyzed further using SOLPS-ITER. It is determined that changes to plasma transport coefficients are required to self-consistently model both plasma and neutral edge dynamics. Pumping is found to modify the poloidal distribution of atomic neutral density, n 0 , along the separatrix, increasing n 0 at the active X-point. Opaqueness to neutrals from high n e in the divertor is found to play a role in mediating neutral penetration lengths and hence, the poloidal distribution of neutrals along the separatrix. Pumped discharges thus require a larger particle diffusion coefficient than that inferred purely from 1D experimental profiles at the outer midplane.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Improved core-edge compatibility using impurity seeding in the small angle slot (SAS) divertor at DIII-D

Impurity seeding studies in the slot divertor (SAS) at DIII-D have revealed a strong relationship between the detachment onset and pedestal characteristics with both target geometry and impurity species. N 2 seeding in the slot has led to the first simultaneous observation of detachment on the entire suite of boundary diagnostics viewing the SAS without degradation of core confinement. SOLPS-ITER simulations with D+C+N, full cross-field drifts and n-n collisions activated are performed for the first time in DIII-D to interpret the behavior. This illustrates a strong effect of divertor configuration and plasma drifts on the recycling source distribution with significant consequences on plasma flows. Flow reversal is found for both main ions and impurities affecting strongly the impurity transport and providing an explanation for the observed dependence on strike point location of the detachment onset and impurity leakage found in the experiments. Matched discharges with either nitrogen or neon injection show that while nitrogen does not significantly affect the pedestal, neon leads to increased pedestal pressure gradient and improved pedestal stability. Little nitrogen penetrates in the core, but a great amount of neon is found in the pedestal consistent with the different ionization potentials of the two impurities. This research demonstrates that neutral and impurity distributions in the divertor can be controlled through variations in strike point locations in a fixed baffle structure. Divertor geometry combined with impurity seeding enables mitigated divertor heat flux balancing core contamination and thus leading to enhanced divertor dissipation and improved core-edge compatibility which are essential for ITER and for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Volumetric recombination in EMC3-EIRENE: Implementation and first application to the pre-fusion power operation phase in ITER

An upgrade to the 3D plasma boundary model EMC3-EIRENE for detached divertor plasmas is introduced and evaluated. Stabilization of the iterative approximation of a self-consistent plasma–neutral gas solution is required at low divertor plasma temperatures, and this has been achieved by linearization of the energy loss term associated with the electron–neutral gas interaction. An implicit method for volumetric recombination has been implemented, which improves agreement with SOLPS-ITER results for the pre-fusion power operation phase in ITER. Resonant magnetic perturbations for control of edge localized modes are found to result in a partially detached divertor plasma with a significantly broader reduction of particle loads in the traditional strike zone, which is attributed to volumetric recombination. The additional non-axisymmetric far scrape-off layer strike points, on the other hand, remain attached at relatively high temperature so that volumetric recombination is not effective there.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mechanisms behind impurity spreading in a channeled, long-leg divertor configuration

One suggested benefit of long-leg, tightly baffled divertor configurations, a common feature included in the designs of many next-generation fusion devices, is access to a delocalized “spread” of impurity radiation along the leg, allowing for a substantial radiation region while maintaining both sufficient distance from the core to prevent pollution of the plasma fuel and avoiding local overheating of the plasma-facing materials. However, it is unclear whether this increased connection length and divertor volume afforded by the long-leg configuration is compatible with impurity transport, and whether it will enable this desired impurity radiation spread at all. Here, the SOLPS 4.3 code is used to assess the distribution of nitrogen and neon impurity radiation in a long-leg, tightly baffled divertor geometry. Further, there is a strong correspondence between the distribution of impurity radiation in the long leg and impurity recycling on the divertor side walls with both the nitrogen simulation set and the neon simulation set, suggesting that impurity spreading along these long legs is due to the interplay between the cross-field transport and associated neutralization of the impurity on the material surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Edge-localized mode mitigation enabled by active control of pedestal density gradient with new EAST tokamak divertor

Mitigation of large edge-localized modes (ELMs) has been achieved by actively reducing the pedestal density gradient with the EAST new right-angled lower divertor through changing the strike point position from the vertical target to the horizontal target. A series of dedicated experiments in the 2021–2024 EAST campaigns demonstrate that this ELM control solution is highly reproducible in a broad parameter space of edge safety factor q 95 = 4.7–7.1, heating power P total = 2.3–5 MW, and pedestal collisionality $ν_{e,ped}^{*}$ = 1–6, under both favorable and unfavorable magnetic configurations. Higher plasma density could facilitate the achievement of this ELM control solution. Statistical results indicate that the ELM mitigation effect can be observed at relatively larger Greenwald density fraction of f GW > 0.47. In addition, this ELM mitigation effect can be achieved with both lithium-coated and boronized metal walls. The pedestal density gradient is systematically lower in the horizontal target case than that of the vertical target case when the ELM mitigation effect can be observed. SOLPS-ITER simulation results indicate that the pedestal fueling from divertor recycling is significantly lower in the horizontal target case. This could contribute to the formation of a flattened pedestal density profile with small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Development of a tomographic reconstruction method for axisymmetric D α emission profiles in the ITER plasma boundary

A tomographic method for reconstructing the axially symmetric D α emission profile in the ITER scrape-off layer and divertor from D α intensity measurements in the fields of view (FoV) of Vis/IR TV and Divertor Impurity Monitor diagnostics is under development. The method takes into account the strong background signal due to stray light from divertor produced by reflections from the metal wall of the vacuum chamber. The method allows to filter out the reflections in D α signals measured in the main chamber. Here, the accuracy of this method is tested within the framework of a synthetic diagnostic, which uses the Raysect and Cherab numerical codes for ray-tracing simulations. The D α emission profiles, simulated with the SOLPS and OEDGE codes, are used as synthetic experimental data. Uncertainty in light reflection properties of the first wall (FW) is the major factor affecting the accuracy of the D α emission reconstruction. The dependence on the uncertainty in the FW surface roughness and on the ratio of diffusive to specular reflection of the error in recovering useful D α signals (the direct light without reflections) in the FoV of the ITER main chamber H-alpha and Visible Spectroscopy Diagnostic is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of particle transport, neutrals and radiation in magnetically-confined plasmas with Aurora

In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Validating the simulation of beam-ion charge exchange in MAST Upgrade

Abstract Simulation of the impact of charge-exchange (CX) reactions on beam ions in the Mega Amp Spherical Tokamak (MAST) Upgrade was compared to measurements carried out with a fission chamber (neutron fluxes) and a fast ion deuterium-alpha (FIDA) diagnostic. A simple model was developed to reconstruct the outer-midplane neutral density based on measurements of deuterium-alpha emission from edge neutrals, and on Thomson scattering measurements of electron density and temperature. The main computational tools used were the ASCOT orbit-following code and the FIDASIM code for producing synthetic FIDA signals. The neutral density reconstruction agrees qualitatively with SOLPS-ITER modelling and yields a synthetic passive FIDA signal that is consistent with measurement. When CX losses of beam ions are accounted for, predicted neutron emission rates are quantitatively more consistent with measurement. It was necessary to account for CX losses of beam ions in simulations to reproduce the measured passive FIDA signal quantitatively and qualitatively. The results suggest that the neutral density reconstruction is a good approximation, that CX with edge neutrals causes significant beam-ion losses in MAST Upgrade, typically 20% of beam power, and that the ASCOT fast-ion CX model can be used to accurately predict the redistribution and loss of beam ions due to CX.

Physics↗
DOE OSTI2025

1D model of SOL with self consistent calculation of non-coronalradiation and applications to shallow field line angles

A 1D model of the scrape off layer (SOL) is created with a physics based calculation of the impurity radiation and the inclusion of the perpendicular heat flux q ⊥ compared to previous 1D models. The calculation of the impurity radiation includes the addition of a model using parallel and perpendicular motion to calculate the impurity confinement time τ and the inclusion of charge exchange between neutral hydrogen and impurity ions. Reasonable agreement is found between the model and the 2D fluid code SOLPS-ITER. Additionally, an improvement is shown over the constant τ used in previous works The model is used to examine the effects of shallow field line angles at ITER levels of parallel heat flux (q || ). A large increase in radiation results in an increase in the upstream parallel heat flux which can still be in detachment at shallow field line angles for both q ⊥ and the calculated τ for traditional impurities (argon and nitrogen) as well as boron. The physics behind this increase in radiation at shallow field line angles is shown to be a decrease in τ and a change in the temperature profile with the inclusion of q ⊥ .

ITER↗
DOE OSTI2020

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neutral leakage, power dissipation and pedestal fueling in open vs closed divertors

The effect of divertor closure on neutral leakage, power dissipation and pedestal fueling is explored using the SOLPS code for different densities including attached and detached divertor conditions for DIII-D open and closed divertors. Here, a metric for divertor closure in terms of the percentage of neutrals escaping the divertor is presented. At the detachment onset the population of neutrals escaping the divertor is 12% for the closed divertor and 35% for the open divertor. This results in a lower contribution of the convective term to the total heat flux and increased power losses in the closed divertor. The comparison between the two geometries reveals that for the same upstream density, the upstream temperature is 23% higher in the closed divertor. The physical mechanisms for detachment at lower upstream density in the closed divertor are also identified with the reduction of the convective term directly resulting in 40% decrease and the increased power losses contributing with a factor of 2. The baffling also lowers the neutral flux refueling the pedestal resulting in further reduction of the core density for achieving detachment. The results introduced here show that the divertor recycling is the primary fueling mechanism in DIII-D. It is shown that fueling occurring in the region between the X-point and the outer midplane is at least four times lower in the closed divertor due the baffling structure which enhances target neutral trapping and consequently reduced upstream leakage. Due to the openness of the inner leg to the common flux and HFS regions both open and closed divertor experience almost identical fueling at the HFS highlighting the need to also consider the inboard divertor for shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗