Search NASA⌕ Search

Engineering topics

Moulton, D.

Publications and source records attributed to Moulton, D..

11 records
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Pedestal particle balance studies in JET-ILW H-mode plasmas

JET-ILW type I ELMy H-modes at 2.5 MA/2.8 T with constant NBI heating (23 MW) and gas fuelling rate were performed, utilising edge localised mode (ELM) pacing by vertical kicks and plasma shaping (triangularity, δ) as tools to disentangle the effects of ELMs, inter-ELM transport and edge stability on the pedestal particle balance. In agreement with previous studies, the pedestal confinement improves with increasing δ, mostly due to a significant increase in pedestal density while the ELM frequency ($f_{\mathrm{ELM}}$) is decreased. Improved pedestal confinement with increasing δ was observed even when the pedestal MHD stability was degraded artificially by vertical kicks, implying that increased triangularity may favourably affect the inter-ELM pedestal recovery. The workflow developed to quantify the pedestal particle balance uses high time-resolution profile reflectometry to characterise the inter-ELM evolution of the plasma particle content ($\mathrm{d} N/\mathrm{d}t$), the NEO drift-kinetic solver to evaluate the neoclassical fluxes and interpretative EDGE2D-EIRENE simulations to estimate the edge particle source. The edge particle source is then constrained by deuterium Balmer-α line intensity measurements in the main chamber, which are, however, strongly affected by reflections from the metal walls. The reflections are accounted for by the CHERAB code taking the divertor emission (the brightest light source in the torus) distribution from imaging spectroscopy measurements as input. Our analysis shows that in the second half of the ELM cycle, the volume-integrated particle source is larger than $\mathrm{d} N/\mathrm{d}t$, indicating that transport plays a key role in the inter-ELM pedestal recovery.

fuelling↗
DOE OSTI2023

Novel SOLPS-ITER simulations of X-point target and snowflake divertors

Abstract The design and understanding of alternative divertor configurations may be crucial for achieving acceptable steady-state heat and particle material loads for magnetic confinement fusion reactors. Multiple X-point alternative divertor geometries such as snowflakes and X-point targets have great potential in reducing power loads, but have not yet been simulated widely in codes with kinetic neutrals. This paper discusses recent changes made to the SOLPS-ITER code to allow for the simulation of X-point target and low-field side snowflake divertor geometries. Snowflake simulations using this method are presented, in addition to the first SOLPS-ITER simulation of the X-point target. Analysis of these results show reasonable consistency with the simple modelling and theoretical predictions, supporting the validity of the methodology implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

Abstract MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimizing detachment control using the magnetic configuration of divertors

Abstract As tokamak research moves to reactor conditions, the control of a stable, optimally-detached divertor plasma has become increasingly relevant. Simple predictions of such detachment control have been performed previously using the detachment location sensitivity (DLS) model (Lipschultz et al 2016 Nucl. Fusion 56 056007). In this study the DLS model is extended and combined with SOLPS-ITER simulations of isolated divertor grids to study the effects of alternate divertor magnetic field properties on detachment control. The DLS model predicts that divertors can achieve easier access to detachment through a long connection length, a high total flux expansion, and a high average magnetic field in the divertor compared with that at the x -point. SOLPS-ITER simulations suggest an even stronger impact of total flux expansion and connection length on detachment access than the DLS model predicts. In terms of detachment evolution, both simulation and modelling show a high gradient in the total magnetic field and low B p o l B are able to more easily keep a detachment front at a desired poloidal location. Regions with no stable detachment front locations can arise for high magnetic field gradients directed towards the target. Significant differences have been found between impurity scan simulations and the DLS model. These differences may be attributed to sources and sinks of power and electron pressure, and should be explored further in future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of deuterium and carbon radiation transport in MAST-U tokamak advanced divertors

Modest effects of deuterium and carbon radiation opacity in the Super-X and snowflake divertor plasmas are predicted for MAST Upgrade tokamak with core plasma input power 2.5–5 MW and plasma current 1 MA. The radiation transport modeling is based on the SOLPS-EIRENE and UEDGE code divertor plasma predictions. Two radiation transport models are used: one is based on a full radiation transport equation implemented in the radiation transport and collisional-radiative code CRETIN (without feedback on the background plasma), and another is an internal self-consistent UEDGE model with ionization, recombination, and heating rates corrected for Ly α line trapping based on the escape probability model implemented in CRETIN. In MAST-U, the Super-X and snowflake divertor plasmas are predicted to reach detached regimes at lower upstream densities than the standard divertor, and the conclusion still holds with radiation transport effects included. At neutral densities m -3 , modest Ly α deuterium line trapping with optical depths 10–15 is predicted in the Super-X divertor. Divertor plasmas are optically thin to other Lyman and Balmer lines, as well as to strong C III and C IV lines that are responsible for most of divertor radiated power. Insignificant changes (within a few percent) to divertor deuterium ionization and recombination rates are found. Radiation fluxes on outer divertor target are modified within a factor of 2–3 when the radiation transport is accounted for, and a similar variation is found due to the line shape models that define the absorption and emission line profiles in the radiation transport modeling. The predicted Lyman and Balmer spectral intensities are significantly modified due to radiation trapping. A measurement of divertor radiation transport effects is discussed using the Ly β /Ba α line ratio. In the snowflake divertor configuration, divertor plasmas are found to be optically thin to Lyman series lines within a large range of parameter variations that include upstream density, divertor transport coefficients, and magnetic configurations. Finally, modest radiation transport effects are only found in a few cases with strongest divertor transport and magnetic configurations closest to the ideal snowflake configuration, however, the plasma background models that were used are yet to be validated with an experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Modeling Snowflake Divertors in MAST-U Tokamak

We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Comparison of private flux region instability in conventional and super-X divertor configurations

Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗