Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “divertor geometry”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Updated Thermofluid Performance of the Simplified Flat Variant of the HEMJ

Our group has recently developed and studied “finger”-type divertors that are a simplified version of the helium-cooled modular divertor with multiple jets (HEMJ) using coupled computational fluid dynamics and thermal stress simulations. Such a simplified geometry could reduce complexity and cost given the large number of fingers required to cover the total divertor target area. Previous experimental studies for this simplified flat design reported lower heat transfer coefficients and higher pressure drops than the HEMJ, contrary to numerical predictions. Subsequent measurements determined that the original test section had significant dimensional variations in the jet exit holes. A new test section was therefore manufactured and tested in the Georgia Tech (GT) helium loop. The experimental results presented here for this test section at maximum heat flux of 7.1 MW/m 2 are in good agreement with numerical predictions. Correlations developed from these experimental data are extrapolated to predict the maximum heat flux that can be accommodated by the flat design and the coolant pumping power requirements under prototypical conditions. Lastly, numerical simulations are used to estimate the sensitivity of the flat design to geometric variations typical of manufacturing tolerances and variations in the gap width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2022

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Shadow masks predictions in SPARC tokamak plasma-facing components using HEAT code and machine learning methods

Here, this work uses machine learning (ML) to complement HEAT (Heat flux Engineering Analysis Toolkit) by developing 3-D footprint surrogate models for fast and accurate heat load calculations in the divertor of the SPARC tokamak. The focus is on shadowed regions, or magnetic shadows, caused by the 3-D geometry of plasma-facing components (PFCs). ML classifiers are employed to create a surrogate model for HEAT generated shadow masks, predicting these shadow masks and divertor heat flux profiles based on a diverse range of equilibria and only the plasma current, safety factor(q95) at the edge, and magnetic flux angles as input parameters. The ultimate goal is to integrate the model for real-time control and future operational decisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design considerations for optimizing the transient CHI injector on QUEST

Transient coaxial helicity injection (transient CHI) on the QUEST (Q-shu University Experiment with Steady-State ST) spherical tokamak (ST) has recently validated the floating biased electrode configuration for solenoid-free plasma startup. In support of a significant divertor upgrade on QUEST, the details of the transient CHI injector geometry on QUEST have been examined more closely using tokamak simulation code (TSC) simulations. QUEST uses a HIT-II-like (Helicity Injected Torus) injector configuration in which the injector region is comprised of coaxial injector electrodes located in the lower part of the machine. TSC simulations indicate high injector flux operation may benefit from an increased gap between the coaxial electrodes. This is the first study to examine the differences between an open and closed electrode configuration for solenoid-free plasma startup using transient CHI. Results show that both configurations can generate similar levels of closed flux, but the closed configuration may be easier to implement in some ST designs as a much smaller portion of the injector needs to be close to the injector flux coil. Results show that increasing the electrode gap width from the present 10.8 cm to about 15–20 cm would increase the closed flux fraction by about 40%. The results presented in this paper are generally applicable to the CHI design for other STs.

coaxial helicity injection↗
DOE OSTI2023

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↗
DOE OSTI2025

Active Divertor Heat Flux Control using Impurity Powder Dropper

Divertor plasma-facing components (PFCs) in a tokamak are typically designed to withstand average steady-state heat loads of about 5–10 MW/m 2 , a limit that applies to both solid and liquid lithium (LL) PFCs. Exceeding these design values can result in surface damage to tungsten PFCs or excessive lithium (Li) evaporation in liquid lithium divertor (LLD) PFCs. Since exceeding the divertor heat load limits has serious consequences, it is therefore prudent to develop a tool to reduce the divertor heat load and bring the heat load to within the design limit without affecting the plasma performance. Active low Z impurity injection such as Li has been suggested as a potential solution to mitigate excess heat flux as suggested previously, given that non-coronal radiation can be quite large ~ 20–30 MJ per mole of injected Li. Li is considered desirable for reducing the edge neutral recycling helping to improve plasma energy confinement. In this paper, we model the Impurity Power Dropper (IPD) to investigate its potential of divertor heat flux control. The IPD is typically located at the top of the tokamak device and uses a vertical drift tube of a few meters. In the 2 m drift tube case, the IPD powder is accelerated to ~ 6 m/sec before reaching the plasma with the upper divertor configuration, matching the condition for the in-board side pellet injection case. By modeling the IPD geometry we determined the IPD powder deposition profile, and thus the non-coronal radiation and ionization profiles in time as well. From the enhanced radiation power loss, it is therefore possible to reduce the divertor heat load using the divertor simulation code. In conclusion, the IPD divertor heat flux control can be tested in the facilities with IPD including ST-40, DIII-D, EAST, WEST and NSTX-U.

Active lithium injection↗
DOE OSTI2026

Power and particle exhaust for the ARC fusion power plant

To successfully show that fusion is an attractive energy source, the ARCTM fusion power plant will need to operate with a robust, integrated power and particle exhaust solution. To maximise ARC’s fusion power output while avoiding excessive erosion of the plasma-facing components, we will need to radiatively dissipate most of the power crossing the last-closed flux surface, injecting radiating impurities such as argon or neon to access divertor detachment. Divertor detachment will need to be integrated with a high-performance core plasma, and with efficient impurity pumping to prevent the accumulation of helium ash in the core. To access and control detachment in high-performance plasmas, we have designed ARC with up–down-symmetric divertors, with secondary X-points in long, tightly baffled outer legs. Using a core-edge modelling workflow, we predict that with this divertor design, ARC will access detachment with modest argon seeding in the divertor (c Ar,div ∼0.9%), which should have minimal impact on the core ( Z eff,core <0.5) for reasonable argon enrichment (c Ar,div /c Ar,core =6). Due to the high predicted divertor neutral pressure (p div ∼20 Pa), we predict that ARC will sufficiently pump helium to limit ash accumulation in the core (c He,core <2%) for a helium enrichment of c He,div /c He,core =0.4. ARC’s divertor design is expected to increase the stability of a detachment front in the outer divertor leg, helping to prevent divertor reattachment during smaller heat-flux transients such as scrape-off-layer filaments associated with the quasi-continuous exhaust regime. However, this buffering will not be sufficient to prevent divertor reattachment during large type-I edge-localised modes (ELMs), and as such these will need to be avoided on ARC. Experiments on SPARC will be used to select an integrated scenario which avoids or mitigates type-I-ELMs while maintaining access to detachment, good core fusion performance and sufficient impurity exhaust. SPARC experiments will also be used to finalise ARC’s divertor design, by studying the impact of magnetic and first-wall geometry on detachment stability, impurity enrichment and neutral baffling under conditions similar to those expected for ARC. In conclusion, our analysis finds that ARC will have a viable power and particle exhaust solution which is compatible with high-power operations, and this solution will be validated in experiments on SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2025

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Predicting roughness effects in additively manufactured coolant channels with helical enhancements

Additive manufacturing (AM) is a promising technique for fabrication of complex geometries such as those expected to be utilized in the blanket, first wall, and divertor. In the case of cooling, metallic AM may be exploited to embed geometric enhancements (ribs, rifling, etc.) to improve cooling performance. However, due to the roughness of these unfinished internal AM surfaces, prediction of thermal hydraulic performance in such channels is difficult. In this work, we consider a methodology for predicting pressure drop and heat transfer in AM channels containing helical enhancements (e.g. helical ribs, twisted tapes) that allows the incorporation of roughness data through conventional pipe flow correlations. This methodology is tested using experimental friction factor and heat transfer coefficient data from high-pressure helium coolant flow measurements in AM stainless steel tubes fabricated by laser powder bed fusion. Both a featureless AM tube and one containing helical ribs were considered alongside a conventionally manufactured smooth tube. The AM surface roughness is obtained by profilometry and used to predict an equivalent sand-grain roughness, with this equivalent roughness confirmed through AM featureless tube measurements. Under the proposed methodology, this roughness information is incorporated into predictions of friction factor and Nusselt number for the rifled tube. Furthermore, these predictions agree well with experimental data across a large range of Reynolds numbers, encouraging the use of this methodology for thermal hydraulic analysis of similar systems and design of future coolant channel geometries.

Additive manufacturing↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Heuristic predictions of RMP configurations for ELM suppression in ITER burning plasmas and their impact on divertor performance

Abstract A subspace of resonant magnetic perturbation (RMP) configurations for edge localized mode (ELM) suppression is predicted for H-mode burning plasmas at 15 MA current and 5.3 T magnetic field in ITER. Perturbations to the core plasma can be reduced by a factor of 2 for equivalent edge stability proxies, while the perturbed plasma boundary geometry remains mostly resilient. The striated domain of perturbed field lines connecting from the main plasma (normalized poloidal flux < 1 ) to the divertor targets is found to be significantly larger than the expected heat load width in the absence of RMPs. This facilitates heat load spreading with peak values at an acceptable level below 10 MW m − 2 on the outer target already at moderate gas fueling and low Ne seeding for additional radiative dissipation of the 100 MW of power into the scrape-off layer (SOL). On the inner target, however, re-attachment is predicted away from the equilibrium strike point due to increased upstream heat flux, higher downstream temperature and less efficient impurity radiation.

3D plasma boundary modeling↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗