Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma devices”

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 667 records · Page 37

2-D moving mesh modeling of lithium dryout in open surface liquid metal flow applications

Liquid lithium displays increasing promise as a replacement for solid plasma facing components (PFCs) in fusion device applications. The Liquid Metal Infused Trench (LiMIT) system, developed at the University of Illinois (UIUC), has demonstrated how thermoelectric magnetohydrodynamics (TEMHD) can be harnessed to drive liquid lithium flow in an open surface PFC. However, in the highest heat flux applications, large local acceleration is created via TEMHD, and the sudden increase in velocity can cause the liquid level to expose the underlying solid, eliminating the protective benefits of the lithium. In order to study potential mitigation strategies, a 2-D COMSOL Multiphysics model was developed using the moving mesh module to capture free surface flow. The model depicts the development of the dryout phenomenon for 2 test cases – slow (1 cm/s) and medium (10 cm/s) flow in 5 mm deep trenches – including the liquid level reduction under the high heat flux and the pileup of slower flow both upstream and downstream of the heat stripe. The effectiveness of trench shaping dryout mitigation strategies is examined. For the slow flow case, it is shown that a 1.8 mm ledge placed under the heat stripe will stop dryout, and for the medium flow case, a 2.7 mm ledge is required to mitigate the effect. This model can be used to identify strategies for increasing the viable heat load for open surface liquid lithium PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detonation to Deflagration-Mode Transitions in Pulsed Plasma Accelerators

PULSED plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas. The current in such discharges will ionize the gas and produce a strong magnetic field, which interacts with the flowing current to accelerate the plasma through the Lorentz body force. For the present work, two plasma accelerator types employing this general scheme are of interest: the gas-fed pulsed plasma thruster (PPT) 1 and the quasi-steady magnetoplasmadynamic (MPD) accelerator. 2 The gas-fed pulsed plasma accelerator is generally understood as a completely transient device discharging in ∼1-10 μs. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j × B body force, entraining and accelerating propellant it encounters. This process is sometimes referred to in literature as 'snowplowing' the propellant or accelerating the gas in a detonation-mode because the current sheet representation approximates that of a strong detonation shockwave propagating through the gas. For these devices, acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new 'crowbar' discharge at the breech and depriving the initial sheet of additional acceleration. In general, PPTs typically claim thrust efficiencies (ratio of jet kinetic energy to input electrical energy) registering in the teens or lower. 3 In the quasi-steady MPD accelerator, the pulse is lengthened to ∼1 ms or longer and maintained at an approximately constant level during discharge through the use of a pulse-forming network (PFN) of capacitors. After an initial transient discharge, which is typically short relative to the overall discharge period, the plasma assumes a relatively steady-state configuration, known as 'quasi-steady' MPD operation. 2 In this state, ionized gas flows through a stationary current channel in a manner that is sometimes referred to as deflagration-mode operation owing to the similarities to deflagration waves in gases. The plasma experiences electromagnetic acceleration as the plasma flows through the current channel towards the exit of the device. Quasi-steady MPD thrusters claim efficiencies up to 50% for certain propellants.4 There has been significant and sustained research over several decades on both gas-fed PPTs and quasi-steady MPD thrusters, however there have been pulsed thrusters that do not appear to exactly fit either classification, instead possessing a mixture of operational qualities characteristic of both thruster variants. The Coaxial High ENerGy (CHENG) thruster by Cheng, et al.5 operated on the short 10 μs timescales characteristic of PPTs, but claimed the high thrust densities, high efficiencies, and low electrode erosion rates that are more consistent with the MPD/deflagration mode of plasma acceleration. Gas-fed PPT research by Ziemer, et al. 3, 6 identified two separate regimes of performance in those thrusters. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency as a function of mass bit, while the second regime at very low mass bits (termed Mode II) exhibited an increase in efficiency with decreasing mass bit. Work by Poehlmann et al.7 and by Sitaraman and Raja8 sought to understand the performance of the CHENG thruster and the Mode I/Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes of plasma acceleration representing two distinct sets of solutions to the relevant conservation laws. In these works, it was proposed that the values of the various controllable parameters determined whether the accelerator would operate in detonation or deflagration mode. Our hypothesized view of the acceleration process in the CHENG thruster and in PPTs experiencing a transition from Mode I to Mode II is inspired by observations of the transition from the PPT mode of operation to the quasi-steady MPD mode. Specifically, the quasi-steady MPD was discovered by driving a PPT to extended pulse lengths. Above a certain pulse length threshold the transient plasma current sheet transitions into a stable plasma acceleration mode that 'replicates in every observable detail steady flow self-field magnetoplasmadynamic acceleration.' 9 In the present work, instead of treating the accelerator as if it were only operating in a single mode during a pulse, we consider the initial stage of the discharge in all cases as a current sheet forming at the breach of the accelerator and moving towards the exit as a detonation wave. If the current sheet reaches the exit of the accelerator before the discharge is completed, the view of the acceleration mode transitions to the deflagration mode-type found in quasi-steady MPD thrusters. In previous work10 we presented a modeling framework that first captured the time-evolution of the current sheet (detonation) mode of the thruster and then transitioned into the quasi-steady MPD (deflagration) mode of plasma acceleration. In the present work, variations of the controllable parameters - specifically the pulsed circuit properties, the amount of mass injected into the thruster, and the relative timing between the initial gas injection and the initiation of the plasma current sheet - will be used to explore the thruster performance. A range of parameters are explored to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗

Surrogate approximation of the Grad–Shafranov free boundary problem via stochastic collocation on sparse grids

In magnetic confinement fusion devices, the equilibrium configuration of a plasma is determined by the balance between the hydrostatic pressure in the fluid and the magnetic forces generated by an array of external coils and the plasma itself. The location of the plasma is not known a priori and must be obtained as the solution to a free boundary problem. The partial differential equation that determines the behavior of the combined magnetic field depends on a set of physical parameters (location of the coils, intensity of the electric currents going through them, magnetic permeability, etc.) that are subject to uncertainty and variability. The confinement region is in turn a function of these stochastic parameters as well. In this work, we consider variations on the current intensities running through the external coils as the dominant source of uncertainty. This leads to a parameter space of dimension equal to the number of coils in the reactor. With the aid of a surrogate function built on a sparse grid in parameter space, a Monte Carlo strategy is used to explore the effect that stochasticity in the parameters has on important features of the plasma boundary such as the location of the x-point, the strike points, and shaping attributes such as triangularity and elongation. Furthermore, the use of the surrogate function reduces the time required for the Monte Carlo simulations by factors that range between 7 and over 30.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian inferences of electrical current delivered to shocked transmission lines

Small radius Bdot measurements of electrical current delivered by pulsed power devices are routinely compromised by electrode/convolute plasmas endemic to multiterawatt transmission lines. Inferences of delivered current unaffected by these issues have recently been obtained by numerically optimizing consistency between model predictions and local experimental velocimetry data, but these are only unique for shockless velocity profiles. In the paper we describe a more general Bayesian method capable of inferring current despite the presence of shocks. Additionally, we describe uncertainty estimates and use of the technique on experimental data. This technique is the first to provide uncertainty estimates on the full current trace delivered to an inertial confinement fusion target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Editorial for the 28th Workshop on Magnetohydrodynamic (MHD) Stability Control

This special issue contains four papers which were presented at the 28th Workshop on Magnetohydrodymanic (MHD) Stability Control. The workshop was held at Auburn University and on Zoom from 3–5 October 2024. Furthermore, this annual workshop covers the active and passive control of MHD instabilities such as kink, tearing and edge-localized modes in tokamaks, reversed field pinches, spherical tori or other magnetic confinement devices.

Yang, J. [Princeton Plasma Physics Laboratory (PPP↗

Stable heat and particle flux detachment with efficient particle exhaust in the island divertor of Wendelstein 7-X

The island divertor concept is an innovative and promising idea to handle heat and particle exhaust in stellarators. At the Wendelstein 7-X (W7-X) stellarator, this divertor concept plays a central role in the device mission to demonstrate reactor relevant plasma confinement for steady-state time scales of up to 30 minutes in the high-performance campaign (OP2) starting in 2022. During the recently concluded first campaign with the inertially cooled island divertor, a large step in the experimental qualification of this divertor concept has been made. In discharges heated with electron cylotron resonance heating of 5-6 MW, central densities in the range of 0.7-1.2 × 10 20 m -3 have been reached in combination with full divertor heat flux detachment. Also, significant neutral gas pressures and neutral compression ratios were shown for the first time in combination with reduced divertor particle flux. The divertor heat loads drop by an order of magnitude from >5 MW m -2 to below 0.5 MW m-2 with increasing density, and substantial compression of neutrals reaching neutral pressure in the sub-divertor volume of >6.0 × 10 -4 mbar was seen. These elevated neutral pressure levels can be obtained and maintained with an up to 80% reduction of the particle fluxes onto the divertor target tiles. This discharge scenario was held stably detached for up to 28 seconds, which is equivalent to several hundred energy confinement times τ E and longer than the time scales for current relaxation. No impurity accumulation was seen at constant Z eff ≈ 1.5 and the stored energy stayed constant at levels of W dia >600 kJ. The level of neutral pressure and compression reached in this scenario extrapolates well to the steady-state particle exhaust requirements for high-performance steady-state operation in OP2, in which the fully actively cooled high-heat-flux divertor will be available. Finally, an overview of this recently discovered divertor regime is given and the status of the physics understanding based on modeling of these regimes with the EMC3-EIRENE code is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The dependence of tungsten fuzz layer thickness and porosity on tungsten deposition rate and helium ion fluence

Abstract Fuzz formation on a heated tungsten surface in the presence of a helium-containing plasma and tungsten deposition source was investigated. Tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, W/He ion flux ratio of ∼ 0.4 × 10 − 4 , and varied helium ion fluence from 0.18 to 3.4 × 10 26 m −2 . Fuzz thickness was measured by cross-sectional scanning electron microscopy to increase from 0.22 to 15 µ m with increasing helium ion fluence. No indication of saturation in fuzz thickness at high fluence was observed, in contrast to fuzz produced on a tungsten surface without tungsten deposition. Additional tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, helium ion fluence of ∼ 3.4 × 10 26 m −2 , and varied W/He ion flux ratio from 0.26 to 3.0 × 10 − 4 . Fuzz thickness increased from 7.5 to 120 µ m with increasing W/He ion ratio. A final sample exposed at 1123 K to a mixed helium-deuterium plasma with ion incident energy of 76 eV, helium ion fluence of 0.18 × 10 26 m −2 , and W/He ion flux ratio of 2.2 × 10 − 4 developed nearly identical fuzz structures to that developed in a pure He plasma. As a function of deposited tungsten fluence, all results were found to trace out a single layer-growth curve given by a power law relation, indicating that fuzz thickness is independent of the W/He ion flux ratio in the range investigated and independent of any deuterium present in the plasma. As a result, for tungsten plasma facing walls in magnetic fusion devices at 1000–2000 K with 10 −4 W/He ion flux ratio, fuzz with thicknesses greater than hundreds of microns may form in as little as 10 4 s (in the absence of ELM-induced erosion or annealing), and may more significantly affect its thermophysical properties than fuzz generated without a tungsten deposition source.

Physics↗

Acoustic MEMS Sensor Array for Quench Detection of CICC Superconducting Cables

A novel quench detection method using microelectro- mechanical system (MEMS) sensor technology has been investigated in use for high temperature superconducting (HTS) conductors such REBCO tape cables. The sensor array along a superconducting cable, such as a cable-in-conduit-conductor (CICC), is installed in a cooling channel. It will allow sensitive and quick detection for a local quench of a superconducting cable. This work has confirmed that a quench of a single REBCO tape can be detected in liquid nitrogen by a MEMS piezoelectric microphone sensor. The quench detection design utilizing a MEMS sensor array method is discussed for the case of a toroidal field (TF) magnets of a fusion Tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric propulsion ground testing.

Ground testing of electric propulsion devices and systems, discussing electrostatic and plasma propulsion, testing requirements and facilities

ELECTRIC PROPULSION↗

Large-area soft X-ray imaging system

System consists of large-area focusing collector and position sensitive proportional counter. Device can be used to study plasmas, with X-ray imaging in biological sciences, and crystallography.

Gorenstein, P.↗

ECR plasma thruster research - Preliminary theory and experiments

A preliminary theory of the operation of the electron-cyclotron-resonance (ECR) plasma thruster is described along with an outline of recent experiments. This work is presented to communicate the status of an ongoing research effort directed at developing a unified theory to quantitatively describe the operation of the ECR plasma thruster. The theory is presented as a set of nonlinear ordinary differential equations and boundary conditions which describe the plasma density, velocity, and electron temperature. Diagnostic tools developed to measure plasma conditions in the existing research device are described.

Sercel, Joel C.↗

The Use of Langmuir Probes in Non-Maxwellian Space Plasmas

Disturbance of the Maxwellian plasma may occur in the vicinity of a spacecraft due to photoemission, interactions between the spacecraft and thermospheric gases, or electron emissions from other devices on the spacecraft. Significant non-maxwellian plasma distributions may also occur in nature as a mixture of ionospheric and magnetospheric plasmas or secondaries produced by photoionization in the thermosphere or auroral precipitation. The general formulas for current collection (volt-ampere curves) by planar, cylindrical, and spherical Langmuir probes in isotropic and anisotropic non-maxwellian plasmas are examined. Examples are given of how one may identify and remove the non-maxwellian components in the Langmuir probe current to permit the ionospheric parameters to be determined. Theoretical volt-ampere curves presented for typical examples of non-maxwellian distributions include: two-temperature plasmas and a thermal plasma with an energetic electron beam. If the non-ionospheric electrons are Maxwellian at a temperature distinct from that of the ionosphere electrons, the volt-ampere curves can be fitted directly to obtain the temperatures and densities of both electron components without resorting to differenting the current. For an arbitrary isotropic distribution, the current for retarded particles is shown to be identical for the three geometries. For anisotropic distributions, the three probe geometries are not equally suited for measuring the ionospheric electron temperature and density or for determining the distribution function in the presence of non-maxwellian back-round electrons.

Hoegy, Walter R.↗