Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma boundary physics”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Measuring dynamics of differentially rotating, unmagnetized, free-boundary plasma produced by soft x-ray irradiation

We present results from a new experiment, fielded on the MAGPIE pulsed power generator, producing a rotating, boundary-free plasma with no pre-imposed magnetic field. Angular momentum is introduced to the rotating system by the oblique collision of multiple plasma flows, which are driven by x-ray ablation of solid targets using the soft x-rays emitted from the implosion of wire array z-pinches. This produces a hydrodynamically stable plasma which undergoes ~2-3 rotations over the duration of the experiment, significantly more than previous pulsed-power platforms. Estimating the angular frequency from the electron density profile in the inner part of the rotating plasma shows that it has a quasi-Keplerian rotation profile. The system also allows for the addition of a controllable magnetic field. This, combined with the sufficiently large Reynolds (~10 5 ) and magnetic Reynolds (~10) numbers, will enable investigation of the effect of magnetic field on the structure and stability of the rotating plasma.

Physics - Plasma physics↗
DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Observation of quasi-coherent density fluctuation in scrape-off layer enhancing boundary transport in high-βN hybrid plasmas on DIII-D

In this work, we report the observation of a quasi-coherent density fluctuation (QCF) by the Doppler backscattering system in the scrape-off layer (SOL) region of the DIII-D tokamak. This QCF is observed in high power, high performance hybrid plasmas with near double-null divertor (DND) shape during the electron cyclotron heating period. This mode is correlated with a steepened SOL density profile, and leads to significantly elevated particle and heat fluxes between ELMs. The SOL QCF is a long-wavelength ion-scale fluctuation (k θ ρ s ≈0.2-0.4 and k r ρ s ≈0.03), and propagates in the ion diamagnetic direction in the plasma frame. Its radial expanse is about 1.5–2 cm, well beyond the typical width of heat flux λ q on DIII-D. Also, the SOL QCF does not show any clear dependence on the effective SOL collisionality, and thus may raise issues on the control of plasma-material interactions in low collisionality plasmas in which the blob-induced transport is reduced. A linear simulation using BOUT++ with a 5-field reduced model is performed and compared with experimental observations. In simulation results, an interchange-like density perturbation can be driven by the SOL density gradient, and its peak location and the radial width of the density perturbation are in agreement with the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Foreshock Model for Interstellar Shocks of Solar Origin: Voyager 1 and 2 Observations

The Voyager 1 (V1) and Voyager 2 (V2) spacecraft were launched in 1977 on a mission to explore the outer planets and reach the heliopause, the boundary between the hot solar plasma and the relatively cool interstellar plasma. V1 reached the heliopause on 2012 August 25, at 121.6 au, and V2 reached the heliopause on 2018 November 5, at 119.0 au. One of their remarkable discoveries was the detection of shocks propagating into the interstellar plasma from energetic solar events. These shocks are typically preceded by electron plasma oscillations excited by electron beams streaming along interstellar magnetic field lines ahead of the shocks. The frequencies of the plasma oscillations have now provided radial electron density profiles in the outer heliosphere and in the interstellar medium to radial distances of more than 145 au. The oscillations are typically preceded by bursts of high-energy ∼5–100 MeV electrons. These electron bursts are interpreted as being due to the reflection (and acceleration) of cosmic-ray electrons by the shock at the time the shock first contacts the magnetic field line that passes through the spacecraft. Relative timing between the cosmic rays reflected by the shock and the onset of the plasma oscillations allow us, for the first time, to estimate the energy, ∼20–100 eV, of the electron beams responsible for the plasma oscillations. These observations are combined into a self-consistent model called the foreshock model that describes the interaction of shocks of solar origin with the interstellar plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A comparative study of internal kink stability in EU DEMO designs with negative and positive triangularity

Abstract Internal kink (IK) instability is investigated for European demonstration fusion reactor (EU DEMO) plasmas in both negative triangularity (NT) and positive triangularity (PT) configurations. For NT plasmas, the IK becomes more unstable as an ideal conformal wall moves away from the plasma boundary, with the mode growth rate saturating at the wall radial location of about b / a = 1.5 , where a is the plasma minor radius and b the wall radial location. The plasma resistivity destabilizes the IK mode. The effect of sub-sonic toroidal plasma flow is sufficiently weak and can thus be ignored for these EU DEMO equilibria. These results are consistent with those for PT plasmas, albeit with larger mode growth rate in the NT configuration. Both perturbative and self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid calculations predict (partial) stabilization of the IK modes in both NT and PT configurations, with inclusion of various kinetic contributions. Precessional drift motion of trapped fusion-born alphas in EU DEMO produces weak stabilization to the IK mode. Stronger stabilization occurs with the toroidal precession of trapped thermal particles (ions and electrons) and the bounce-transit motion of thermal ions. The stabilization is similar between the NT and PT configurations, due to the similarity of the mode eigenfunction (occupying a nearly circular region in the plasma core) despite the sign difference in the triangularity. The non-perturbative MHD-kinetic hybrid model predicts much less stabilization of the mode than the perturbative model, primarily due to the self-consistent determination of the mode eigenvalue in the former. Generally, no significant difference in the IK mode stability is found between the NT and PT plasmas in EU DEMO.

Physics↗
DOE OSTI2023

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (⁠ ≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. Here, the methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Liquid metal walls

Here, the plasma performance of fusion devices depends strongly on the chosen wall materials. Solid plasma-facing components (PFCs) are predominantly used in present devices, and are the most investigated candidates for fusion designs. High-Z materials such as tungsten (W) are the leading solid PFC material candidates. To date, a material choice that scales to steady-state reactor conditions has not been identified. Moreover, if plasma transient events such as edge-localized modes and disruptions cannot be altogether avoided or sufficiently mitigated, the projected peak heat and particle loads far exceed the power exhaust capabilities of solids. Liquid metal (LM) PFCs represent an intrinsic self-healing boundary that are both resilient to surface damage from transients, and that could handle high steady-state heat and particle fluxes. Flowing LM PFCs can be designed to remove “slag,” the buildup of material erosion due to plasma-material interactions, including charge exchange sputtering in the main chamber. Further, LM offer the prospect to manage hydrogenic species otherwise retained in the PFCs, which is important from a safety and inventory standpoint. The two most promising LMs are lithium (Li) and tin (Sn), although Sn–Li eutectics may be considered. While Sn offers a higher temperature window with low vapor pressure and low hydrogen retention, Li offers the prospect of enhanced energy confinement and higher acceptable core contamination limits, and this section focuses on Li PFCs. An LM PFC development research program developed LM PFC concepts for a nuclear fusion device via engineering design calculations, single-effect experiments, and staged prototypical experiments. A self-consistent design window was identified with liquid Li flow speeds ~5–10 m/s; plasma contamination was negligible for predicted Li evolution rates. While these preconceptual designs hold promise, there is substantial R&D needed to advance the technical readiness levels of LM PFCs for application to fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Strongly intermittent far scrape-off layer fluctuations in Alcator C-Mod plasmas close to the empirical discharge density limit

Abstract Intermittent plasma fluctuations in the boundary region of the Alcator C-Mod device were comprehensively investigated using data time-series from gas puff imaging and mirror Langmuir probe diagnostics. Fluctuations were sampled during stationary plasma conditions in ohmically heated, lower single null diverted configurations with scans in both line-averaged density and plasma current, with Greenwald density fractions up to 0.85. Utilizing a stochastic model, we describe the plasma fluctuations as a super-position of uncorrelated pulses, with large-amplitude events corresponding to blob-like filaments moving through the scrape-off layer (SOL). A deconvolution method is used to estimate the pulse arrival times and amplitudes. The analysis reveals a significant increase of pulse amplitudes and waiting times as the line-averaged density approaches the empirical discharge density limit. Broadened and flattened average radial profiles are thus accompanied by strongly intermittent and large-amplitude fluctuations. Although these filaments are arriving less frequently at high line-averaged densities, we show that there are significant increases in radial far-SOL particle and heat fluxes which will further enhance plasma–wall interactions. The stochastic model has been used as a framework for study of the scalings in the intermittency parameter, flux and mean amplitude and waiting times, and is being used to inform predictive capability for the effects of filamentary transport as a function of Greenwald fraction.

Physics↗
DOE OSTI2020

Experimental investigation of laser ablation of stone polycrystalline targets

We report the results of an experimental investigation of ablation of stone polycrystalline targets of complex multicomponent composition, which imitate the substance of asteroids. The targets were irradiated by nanosecond pulses of a neodymium laser at an energy density Φ{sub L} of up to 5 × 10{sup 4} J cm{sup −2}. The experiments demonstrated the existence of two ablation regimes, with the boundary between them lying at Φ{sub L} ≈ 4000 J cm{sup −2}. The regime change is characterised by a change in the form of the dependence of the surface mass density of removed target material on the laser energy density and by the appearance of a minimum in the dependence of specific energy of destruction on Φ{sub L}. This is supposedly related to the passage from a one-dimensional plasma plume expansion to the three-dimensional one and the corresponding decrease in the efficiency of energy transfer from the laser beam to the target due to a lowering of laser-produced plasma density. Our experiments also showed the existence of a maximum in impulse coupling coefficient C{sub m} as a function of laser energy density (C{sub m} ≈ 6.3 × 10{sup −5} N W{sup −1} for Φ{sub L} = Φ{sub opt} ≈ 25 J cm{sup −2}). Maxima were also recorded in the dependences of the ablation efficiency and average ablation flow velocity on Φ{sub L}. For Φ{sub L} > Φ{sub opt}, the decrease in the function C{sub m}(Φ{sub L}) turns out to be much steeper than for metals and polymer materials. The difference is supposedly due to the lower strength and lower plasticity of the polycrystalline stone targets. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The stability of expanding reactive shocks in a van der Waals fluid

Despite the extensive literature accumulated since the pioneering works of D'yakov and Kontorovich in the 1950s, the stability of steady shocks is still an open question when realistic boundary conditions are accounted. The consideration of a supporting mechanism, which is indeed a necessary condition for shock steadiness, modifies the perturbation shock dynamics in the unstable range. The Noh problem is a suitable example to form steady expanding shocks. This configuration is of great interest to the high-energy-density-physics community because of its direct application to inertial confinement fusion and astrophysics, for which the stagnation of a supersonically converging material via an accretion shock front is ubiquitous. Here in this work, we extend the generalized Noh problem, both base-flow solution and linear stability analysis, to conditions where endothermic or exothermic transformations undergo across the shock. Within the spontaneous acoustic emission conditions found for a van der Waals gas, we find that cylindrical and spherical expanding shocks become literally unstable for sufficiently high mode numbers. Counterintuitively, the effect of exothermicity or endothermicity across the shock is found to be stabilizing or destabilizing, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Synthesis of Large Molecules in Cometary Ice Analogs: Physical Properties Related to Self-Assembly Processes

The combination of realistic laboratory simulations and infrared observations have revolutionized our understanding of interstellar dust and ice-the main component of comets. Since comets and carbonaceous micrometeorites may have been important sources of volatiles and carbon compounds on the early Earth, their organic composition may be related to the origin of life. Ices on grains in molecular clouds contain a variety of simple molecules. The D/H ratios of the comets Hale-Bopp and Hyakutake are consistent with a primarily interstellar ice mixture. Within the cloud and especially in the presolar nebula through the early solar system, these icy grains would have been photoprocessed by the ultraviolet producing more complex species such as hexamethylenetetramine, polyoxymethylenes, and simple keones. We reported at the 1999 Bioastronomy meeting laboratory simulations studied to identify the types of molecules which could have been generated in pre-cometary ices. Experiments were conducted by forming a realistic interstellar mixed-molecular ice (H2O, CH3OH, NH3 and CO) at approximately 10 K under high vacuum irradiated with UV light from a hydrogen plasma lamp. The gas mixture was typically 100:50:1:1, however when different ratios were used material with similar characteristics was still produced. The residue that remained after warming to room temperature was analyzed by HPLC, and by several mass spectrometric methods. This material contains a rich mixture of complex compounds with mass spectral profiles resembling those found in IDPs and meteorites. Surface tension measurements show that an amphiphilic component is also present. These species do not appear in various controls or in unphotolyzed samples. Residues from the simulations were also dispersed in aqueous media for microscopy. The organic material forms 10-40 gm diameter droplets that fluoresce at 300-450 nm under UV excitation. These droplets have a morphology and internal structure which appear strikingly similar to those produced by extracts of the Murchison meteorite. Together, these results suggest a link between organic material photochemically synthesized on the cold grains in dense, interstellar molecular clouds and compounds that may have contributed to the organic inventory of the primitive Earth. For example, the amphiphilic properties of such compounds permit self-assembly into the membranous boundary structures that required for the first forms of cellular life.

Dworkin, Jason P.↗
DOE OSTI2026

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

Experimental study of the influence of gas puff locations on H-mode boundary plasmas with argon seeding on EAST

To investigate the optimal scenario of impurity seeding to obtain divertor plasma detachment for target protection, experiments with Ar&D 2 seeding from two different poloidal locations, the upper outer (UO) divertor target and lower outer (LO) target, were carried out on EAST. Partial energy detachment (the electron temperature near the strike point T e,spt ≤ 10 eV) were obtained with Ar&D 2 mixture puffing from the UO target and, for the first time, from the LO target into H-mode plasmas in the upper single null (USN) configuration. The peak heat flux q t on the UO target was significantly reduced (by ~80%). The rollover of ion flux density j s did not appear probably due to insufficient momentum loss, which is independent of the puff locations. The poloidal asymmetries of particle and heat fluxes on the targets have also been investigated. The UO-dominant asymmetry of particle flux was reversed, while the UO-dominant asymmetry of heat flux was mitigated but not reversed. The plasma confinement dropped by only 14% and 8.2% in the LO-puff case and UO-puff case, respectively, during detachment. The high level of C may contributed greatly to the higher radiation in the bulk plasma region and the greater decline in W MHD in the LO-puff case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗
DOE OSTI2024

Electron transport barrier and high confinement in configurations with internal islands close to the plasma edge of W7-X

The low magnetic shear in the Wendelstein 7-X (W7-X) stellarator makes it feasible to shape the separatrix by the large islands constituting an island-divertor, and this can be exploited to access various magnetic configurations, including samples of different internal island sizes and locations. To investigate the configuration effects on the plasma confinement, a configuration scan was performed by changing the coil currents to vary the rotational transform between values 5/4 and 5/6 at the plasma boundary with different power levels (2, 4, 6 MW) of electron cyclotron resonance heating (ECRH) at a maximum plasma density of 8 × 10 19 m -2 . neutral beam injection (NBI) heating was also applied during some configurations of the scan to create a density ramp and access high densities beyond the X2 ECRH cutoff. For the magnetic configurations, where the 5/5 and 5/6 island chains were moved closer to separatrix but remaining inside the last closed flux surface, the electron cyclotron emission shows that an electron temperature, T e , pedestal develops already during ECRH heated plasma buildup phase indicating a transport barrier, and the barrier sustains irrespective of changed plasma heating conditions such as NBI in the later part of discharge. The transport barrier is broken by subsequent fast crashes, observed with multiple plasma diagnostics with characteristics such as tokamak edge localized modes, and the corresponding crash amplitude and frequency vary with plasma pressure. The impact of the transport barrier on plasma confinement can be seen through the increased core T e profile, which could be responsible for the overall increase in the stored diamagnetic energy by approximately 10% for these configurations. After the plasma heating is terminated, a backwards transition to a degraded confinement state is also observed. These observations indicate a configuration triggered high confinement mode in low shear W7-X. This work focuses on the occurrence of this transport barrier for different magnetic configurations and its relation to internal magnetic islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗