Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA LAYER”

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

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Erosion and impurity transport for the edge localized mode suppression window in KSTAR

A 3D analysis of plasma wall interactions and global impurity transport for the edge localized mode suppression window in KSTAR during H-Mode discharges has shown that carbon erosion at the divertor plates is a strong function of the resonant magnetic perturbation (n = 1) coil current and relative phasing. The Generalized Perturbation Equilibrium Code was used to determine a realistic initial perturbed plasma equilibrium, and EMC3-EIRENE was used to calculate the resulting scrape-off layer plasma used in this study as a fixed background for the ERO2.0 plasma–material interaction model. The resulting transport leads to deposition of impurities along the targets positioned at the high-field side of the device. An attempt at calculating the resulting effective charge state has demonstrated a similar dependence on the perturbation coil current and has been able to determine a window for the experimentally observed values of Zeff by including contributions of all ionized carbon charge states and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surrogate models to optimize plasma-assisted atomic layer deposition in high aspect ratio features

In this work, we explore surrogate models to optimize plasma enhanced atomic layer deposition (PEALD) in high aspect ratio features. In plasma-based processes such as PEALD and atomic layer etching (ALE), surface recombination can dominate the reactivity of plasma species with the surface, which can lead to unfeasibly long exposure times to achieve full conformality inside nanostructures like high aspect ratio vias. Using a synthetic dataset based on simulations of PEALD, we train artificial neural networks to predict saturation times based on cross section thickness data obtained for partially coated conditions. The results obtained show that just two experiments in undersaturated conditions contain enough information to predict saturation times within 10% of the ground truth. A surrogate model trained to determine whether surface recombination dominates the plasma–surface interactions in a PEALD process achieves 99% accuracy. This demonstrates that machine learning can provide a new pathway to accelerate the optimization of PEALD processes in areas such as microelectronics. Our approach can be easily extended to ALE and more complex structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of gas injection location on a lithium vapor box divertor in NSTX-U

Abstract The lithium vapor box divertor is a proposed divertor design that will minimize contamination of the upstream plasma in a fusion device, while also ensuring protection of the target. In this design lithium is evaporated near the target by high temperature lithium surfaces, dissipating the plasma heat flux. The lithium vapor box has been predicted via the fluid-kinetic code scrape off layer plasma simulator (SOLPS-ITER) to achieve low ( n L i / n e ∼ 0.05) upstream concentrations of lithium and low target heat fluxes. Here we compare two choices of deuterium gas puff location using SOLPS-ITER, the private flux region (PFR) and the common flux region (CFR), and find significant differences in the contamination level required to reach an acceptable target heat flux (defined here as q tar max ⩽ 10 MW m −2 ). Deuterium gas puffing from the PFR is seen to better reduce upstream lithium contamination. The difference in puffing location is seen to cause changes in the upstream flow of lithium ions. The PFR puff, having better access to the separatrix, can better reduce the upstream-directed flow of lithium near the separatrix which is the primary source of contamination due to a large thermal force in this region. Puffing from the CFR, partially due to inefficacy at reducing separatrix lithium flow, has higher lithium concentration within the plasma. Solutions that reduce the heat flux to below 10 MW m −2 have a range of lithium concentrations between n L i / n e ∼ 0.01–0.12 depending on puff intensity, location, evaporator temperature and recycling at the various plasma facing components. The efficacy of the puffs is tested for sensitivity to deuterium recycling coefficient at the target, evaporators, and main chamber walls. A CFR located puff is found to be more dependent on the recycling coefficients used than a PFR located puff. regardless of the set of recycling coefficients chosen, PFR puffing achieves lower lithium contamination than CFR puffing for a given heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics↗

Divertor detachment and heat exhaust mitigation control in KSTAR with tungsten divertor

KSTAR has recently undergone an upgrade to use a new tungsten divertor to run experiments in ITER-relevant scenarios. Even with a high melting point of tungsten, it is important to control the heat flux impinging on tungsten divertor targets to minimize sputtering and contamination of the core plasma. Heat flux on the divertor is often controlled by increasing the degree of detachment of scrape-off layer plasma from the target plates. In this work, we have demonstrated successful divertor detachment and heat exhaust dissipation control experiments using two different methods. The first method uses attachment fraction as a control variable which is estimated using ion saturation current measurements from embedded Langmuir probes in the divertor. The second method uses a novel machine-learning-based surrogate model of 2D UEDGE simulation database, DivControlNN. We demonstrated running inference operation of DivControlNN in realtime to estimate heat flux at the divertor and use it as the control variable in a feedback loop with impurity gas flow. We present interesting insights from these experiments including a systematic approach to tuning controllers and discuss future improvements in the control infrastructure and control variables for future burning plasma experiments.

KSTAR tungsten divertor operations↗

Observations of wall conditioning by means of boron powder injection in DIII-D H-mode plasmas

In this paper, we report observations from the DIII-D tokamak indicating that boron (B) powder injection in tokamak plasmas improves wall conditions similarly to glow discharge boronization (GDB). Isotopically enriched B powder (B 11 > 95%) was introduced gravitationally in a sequence of H-mode plasma discharges at rates up to ~160 mg s -1 for durations up to 3 s. Boron injection to cumulative amounts ≤0.1 g appeared to improve wall conditions similarly to boronization, with indications of reduced wall fueling, reduced recycling at the outer strike point and reduced impurity content at breakdown. Post-mortem analysis of graphite samples exposed to far scrape-off layer plasma fluxes during boron injection confirm the formation of a B-C layer, with average surface composition B:C ~ 1. The results suggest that injecting boron-rich powders in tokamak plasmas can effectively replenish boron films on carbon plasma facing components to improve wall conditions and extend the duration of the beneficial effects of GDB.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗

Comparison of AlF 3 thin films grown by thermal and plasma enhanced atomic layer deposition

Films of aluminum fluoride (AlF 3 ) deposited by thermal and plasma enhanced atomic layer deposition (PEALD) have been compared using in situ multiwavelength ellipsometry (MWE) and monochromatic x-ray photoelectron spectroscopy (XPS). The AlF 3 films were grown using cyclic exposures of trimethylaluminum, hydrogen fluoride, and H radicals from a remote H 2 inductively coupled plasma. Films were characterized in situ using MWE and XPS for growth rate, film composition, and impurity incorporation. The MWE showed a growth rate of 1.1 and 0.7 Å per cycle, at 100 °C, for thermal and plasma enhanced ALD AlF 3 films, respectively. Carbon incorporation was below the XPS detection limit. The plasma enhanced ALD AlF 3 film showed the presence of Al-Al chemical states, in the Al 2p scans, suggesting the presence of Al-rich clusters with a concentration of 14%. The Al-rich clusters are thought to originate during the hydrogen plasma step of the PEALD process. Finally, the Al-rich clusters were not detected in thermal ALD AlF 3 films using the same precursors and substrate temperature.

36 MATERIALS SCIENCE↗

Growth of nanostructured molybdenum disulfide (MoS 2 ) thin films on a nanohole-patterned substrate using plasma-enhanced atomic layer deposition (ALD)

Nanostructured molybdenum disulfide (MoS 2 ) thin films were grown on a nanohole-patterned silicon substrate using plasma-enhanced atomic layer deposition. A nanoscale hole-patterned silicon substrate was fabricated for the growth of MoS 2 film using the self-assembly-based nanofabrication method. The nanoscale holes can significantly increase the surface area of the substrate while the formation and growth of nanostructures normally start at the surface of the substrate. Hydrogen sulfide (H 2 S) gas was used as the S source in the growth of molybdenum disulfide (MoS 2 ) while molybdenum (V) chloride (MoCl 5 ) powder was used as the Mo source. The MoS 2 film had a stoichiometric ratio of 1 (Mo) to 2 (S), and had peaks of E 1 2g and A 1g , which represent the in-plane and out-plane vibration modes of the Mo–S bond, respectively. It was found that the MoS 2 film grown in the nanoscale hole, especially at the wall of the hole, has more hexagonal-like structures due to the effects of nanoscale space confinement and the nanoscale interface although the film shows an amorphous structure. Post-growth high-temperature annealing ranging from 800 to 900 °C produced local crystalline structures in the film, which are compatible with those reported by other researchers.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Here, plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Baking and plasma pretreatment of sapphire surfaces as a way to facilitate the epitaxial plasma-enhanced atomic layer deposition of GaN thin films

The growth of high-quality epitaxial gallium nitride (GaN) thin films is achieved by using a baking and plasma pretreatment of the substrate prior to the GaN plasma-enhanced atomic layer deposition (PE-ALD). In this work, it is found that such pretreatment makes the GaN films grow coherently on sapphire substrates, following a layer-by-layer growth mechanism. The deposited GaN film shows high crystalline quality, a sharp GaN/sapphire interface, and a flat surface. The possibility of growing high-quality GaN epilayers in this way broadens the range of applications for PE-ALD in GaN-based devices.

36 MATERIALS SCIENCE↗

Effects of radial transport on divertor power and particle flux widths under different operational regimes in EAST

Here, a study of the effects of radial transport on the outer divertor particle and power flux widths (λ js and λ q ) is carried out using plasmas under different operational regimes in the Experimental Advanced Superconducting Tokamak (EAST). In the EAST experiments, the λ js values measured during the grassy intra- and inter-ELM phases are similar, and the averaged value of λ js during grassy ELMy discharges is found to be smaller than that of the intra-ELM phase of type-I ELMy discharges but larger than that of the inter-ELM phase of type-I ELMy discharges. Simulations of scans of the radial particle and heat transport coefficients, D and χ e , performed using the BOUT++ transport code show that both λq and λjs increase with D and χe, especially when D or χ e is larger than a threshold value, indicating that background turbulence starts to matter in the determination of λ q and λ js when D or χ e surpasses the threshold value. A comparison between the simulation of the D scan and the experimental results shows that the different values of λ js obtained under different plasma operational regimes are probably due to the different intensities of background turbulence, which is beyond the scope of Goldston's heuristic drift-based model as reported by Goldston et al (2012 Nucl. Fusion 52 013009), since their model was derived by assuming that turbulent transport does not exceed a certain level. The ratio of λ js to λ q is greatly affected by the radial transport; however, simulations carried out using Scrape-off Layer Plasma Simulation (SOLPS) show that divertor geometry and plasma density have big influences on λ js , and thus could also affect the ratio of λ js to λ q .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.

divertor physics↗

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗

Theoretical Modeling of Reactor Relevant Conditions for Plasma Jet Driven Magneto-Inertial Fusion

The Charger Advanced Power and Propulsion Laboratory (CAPP), a laboratory within the Propulsion Research Center (PRC) at the University of Alabama in Huntsville (UAH) is working with Los Alamos National Laboratory (LANL). to develop models and inform on promising paths for high gain magneto-inertial fusion (MIF) conditions. This report provides a framework for identifying promising conditions for achieving ignition in plasma-jet-driven magneto-inertial fusion (PJMIF)[1]. As proposed, for the first part of the contract, UAH proposes to develop a gain over unity set of stagnation conditions to provide a state of plasma conditions to achieve to set long terms goals for the PJMIF program. Specifically, UAH will model PJMIF stagnation conditions to include radiation, heat transfer, two temperature energy equations, fusion reactivity and nonlocal fusion product deposition, but no hydrodynamics for these purposes. These calculations will use a stationary plasma model to reduce simulation complexity—focusing on a DT target at 10 keV. Subsequent work will include a DD plasma layer acting as an afterburner. UAH will assume an initial magnetic field without any consideration of the topology, just assume a field strength, most likely scaled with consideration of the local hall parameter. This effort will inform the team on the tradeoff between mass, peak target field, etc and the achievable gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma enhanced atomic layer deposition of thin film Li 1+x Mn 2-x O 4 for realization of all solid-state 3D lithium-ion microbatteries

Here, the plasma enhanced atomic layer deposition (PEALD) process for MnO 2 was demonstrated, exhibiting self-limiting growth as well as stable composition and stable growth rate over a temperature window of 205–265 °C. The PEALD process for MnO 2 was combined with the thermal ALD process for LiOH to synthesize Li 1+x Mn 2-x O 4 thin film cathodes, where the stoichiometry was effectively controlled to allow for crystallization in the electrochemically active spinel phase. A 3D nanostructure consisting of an Li 1+x Mn 2-x O 4 cathode layer, ALD Li x Al y Si z O solid electrolyte, and SiGe nanowire anode was fabricated and characterized via TEM. The Li 1+x Mn 2-x O 4 thin films maintained 66% of the areal capacity upon a 100× increase in the rate (4–360 μ A cm -2 ) as well as 97% capacity retention over 100 cycles at ~5C. The ALD Li 1+x Mn 2-x O 4 thin films exhibited a volumetric capacity of 52 μ Ah cm -2 μ m -1 at a C-rate of ~0.5C and, coupled with the high operating voltage (4.0 V), offer some of best areal energy densities for ALD thin film cathodes, making it a viable material for integration with 3D lithium-ion microbatteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗