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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Massive, long-lived electrostatic potentials in a rotating mirror plasma

Abstract Hot plasma is highly conductive in the direction parallel to a magnetic field. This often means that the electrical potential will be nearly constant along any given field line. When this is the case, the cross-field voltage drops in open-field-line magnetic confinement devices are limited by the tolerances of the solid materials wherever the field lines impinge on the plasma-facing components. To circumvent this voltage limitation, it is proposed to arrange large voltage drops in the interior of a device, but coexist with much smaller drops on the boundaries. To avoid prohibitively large dissipation requires both preventing substantial drift-flow shear within flux surfaces and preventing large parallel electric fields from driving large parallel currents. It is demonstrated here that both requirements can be met simultaneously, which opens up the possibility for magnetized plasma tolerating steady-state voltage drops far larger than what might be tolerated in material media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The 3D magnetic topology and plasma dynamics in open stochastic magnetic field lines

We report the thermal quench triggered by locked modes is known to be mainly due to open stochastic magnetic field lines connected to the wall boundary. It is essential to understand the 3D structure of open stochastic field lines since it determines the overall plasma dynamics in the system. In this study, we analyze the 3D magnetic topology for two key concepts, the connection length L c and the effective magnetic mirror ratio M eff , and present a comprehensive picture of electron and ion dynamics related to the magnetic topology. The connection length determines the 3D structure of the ambipolar potential, and a sharp potential drop across distinct L c regions induces the E × B transport and mixing across the field line. The confinement of electrons and ions along the field line is determined by the ambipolar potential and [Formula: see text] configuration. Electron and ion temperatures in magnetic hills (M eff < 1) are lower than in magnetic wells (M eff > 1) because particles in magnetic hills are more likely to escape toward the wall boundary along the field line. The mixing between the magnetic wells and hills by E × B and magnetic drift motions results in collisionless detrapping of electrons and ions, which reduces their temperature efficiently. Numerical simulations of two different magnetic configurations demonstrate the importance of the collisionless detrapping mechanism, which could be the main cause of plasma temperature drop during the thermal quench.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Densities and abundances of hot cometary ions in the coma of P/Halley

On its flight by P/Halley, the Giotto spacecraft carried a High Energy Range Spectrometer (HERS) for measuring the properties of cometary ions picked up by the solar wind in the nearly collisionless regions of the coma. Preliminary estimates of the ion densities observed by HERS were reevaluated and extended; density profiles along the Giotto trajectory are presented for 13 values of ion mass/charge. Comparison with the physical-chemical model of the interaction of sunlight and the solar wind with the comet by other researchers reveals that, with the exception of protons and H2(+), all ion densities were at least an order of magnitude higher than predicted. The high ion densities cannot be explained on the basis of compression of the plasma, but require additional or stronger ionization mechanisms. Ratios of the densities of different ion species reveal an overabundance of carbonaceous material and an underabundance of H2(+) compared to the predictions of the Schmidt. While the densities of solar wind ions (H(+) and He(++)) changed sharply across a magnetic discontinuity located 1.35(10)(exp 5) km from the comet, this feature, which has been called both the 'cometopause' and the 'magnetic pileup boundary' was barely distinguishable in the density profiles of hot cometary ions. This result is consistent with the interpretation that the magnetic pileup boundary detected by Giotto was caused by a discontinuity in the solar wind and is not an intrinsic feature of the interaction of the solar wind with an active comet.

Neugebauer, M.↗
NASA NTRS2014

X-Ray Observations of VY Scl-Type Nova-Like Binaries in the High and Low State

Four VY Scl-type nova-like systems were observed in X-rays during both the low- and the high-optical states. We examined Chandra, ROSAT, Swift and Suzaku archival observations of BZ Cam, MV Lyr, TT Ari and V794 Aql. The X-ray flux of BZ Cam is higher during the low state, but there is no supersoft X-ray source (SSS) as hypothesized in previous articles. No SSS was detected in the low state of the any of the other systems, with the X-ray flux decreasing by a factor between 2 and 50. The best fit to the Swift X-ray spectra is obtained with a multicomponent model of plasma in collisional ionization equilibrium. The high-state high-resolution spectra of TT Ari taken with Chandra Advanced CCD Imaging Spectrometer (ACIS-S) and the Chandra High Energy Transmission Grating (HETG) shows a rich emission line spectrum, with prominent lines of Mg, Si, Ne and S. The complexity of this spectrum seems to have origin in more than one region, or more than one single physical mechanism. While several emission lines are consistent with a cooling flow in an accretion stream, there is at least an additional component. We discuss the origin of this component, which is probably arising in a wind from the system. We also examine the possibility that the VY Scl systems may be intermediate polars, and that while the boundary layer of the accretion disc emits only in the extreme ultraviolet, part of the X-ray flux may be due to magnetically driven accretion.

accretion↗
DOE OSTI2023

X-ray imaging methods for high-energy density physics applications

Large scale high-energy density science facilities continue to grow in scale and complexity worldwide. The increase in driver capabilities, including pulsed-power and lasers, continue to push the boundaries of temperature, pressure, and densities, opening up new physics regimes. X-ray imaging is one of the many diagnostic techniques that are used to probe states of matter in these extreme conditions. Improved fabrication and polishing methods have provided improved x-ray microscope performance, while improving detector and x-ray sources now enable pico-second imaging with few micron resolutions. Here this Review will cover x-ray imaging methods, primarily absorption imaging, and their improvements over the last few decades.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗
NASA NTRS2002

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Khazanov, G. V.↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development of process parameters and post-build conditions for qualification of LPBF 316 SS

To harness the potential of laser powder bed fusion (LPBF) 316H stainless steel (SS) for use in advanced nuclear reactors, extensive research efforts are needed to develop optimal laser processing parameters and appropriate heat treatments, taking into account various manufacturing platforms and locations. Furthermore, it is crucial to assess the material properties in environments relevant to reactor operating conditions. In collaboration with Oak Ridge National Laboratory (ORNL) and Los Alamos National Laboratory (LANL), Argonne National Laboratory (ANL) is dedicated to gaining insights into how the manufacturing process and post-build treatments impact the performance of LPBF 316H SS. This report provides an overview of ANL's research findings for FY23, focusing on LPBF 316H SS produced using a Renishaw AM400 laser system. Specifically, our research has concentrated on three key areas: 1) Optimum laser processing parameters: Leveraging the high-throughput printing technology, we systematically varied the laser power, exposure time and point distance within the same build. By measuring the as-printed porosity, we identified an optimum printing parameter window that produced materials with near-complete density; 2) Post-built treatment development: Materials printed with the optimum laser parameters underwent treatments including stress relief, solution annealing, and hot isostatic pressing. Subsequently, we fabricated tension, creep, and fatigue specimens from materials subjected to these different conditions. Over the next year, a series of scoping tests will be conducted to facilitate the selection of the most suitable post-build treatment condition; 3) Thermal aging effects on microstructure and properties: thermal aging at 550°C, 650°C and 750°C was conducted on LPBF 316H SS specimens up to 2500 h. Microstructural characterization was performed with electron microscopy, and the evolution of the dislocation cell structures and secondary phases was studied. Microhardness tests and tension tests were conducted to quantify changes in mechanical properties. The results revealed that the aged LPBF 316H SS exhibited a high density of Cr-rich and Mo-rich fine precipitates within grains due to the presence of a high density of dislocations and the globally distributed dislocation cell structures, which served as heterogeneous nucleation sites for secondary phases. Along grain boundaries, a Mo-Cr-rich phase, observed in 750°C-aged conditions, displayed significant thermal coarsening. Throughout the aging process, various mechanisms, including stress relaxation, dislocation cell structure recovery, solution hardening, and precipitation hardening, competed to influence materials strength. These results contribute towards establishing the technical basis for qualifying LPBF 316H SS for advanced nuclear reactor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Two-Frequency RF Fields Induced Multipactor in Coaxial Transmission Lines

Multipactor is a nonlinear discharge phenomenon that occurs in vacuum RF systems, potentially leading to signal distortion, power loss, and even permanent damage to high-power components. This study presents a detailed investigation of two-surface multipactor in coaxial transmission lines under two-frequency excitation, using one-dimensional (1D) Monte Carlo simulations validated by three-dimensional (3D) Particle-in-Cell (PIC) results and experimental data. Introducing a second carrier mode is shown to suppress multipactor by reshaping and shrinking the susceptibility region, with the extent and location of suppression strongly dependent on the device aspect ratio and the relative phase of the second mode. Distinct suppression patterns emerge across different frequency–gap distance (fd) regimes, and in certain cases, susceptibility expansion is also observed. The study identifies and distinguishes pure and mixed multipactor modes in coaxial geometry, where analytical mode boundaries are not readily defined. Unlike planar systems, pure-mode regions in coaxial structures overlap with mixed-mode domains, complicating classification. Image charge forces are found to have minimal effect on susceptibility thresholds but do influence electron growth rates. These findings provide new insights into waveform-driven control of multipactor in high-power RF systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Performance Comparison of Helium and sCO 2 as Coolants for Modular Finger-Type Divertors

Several studies at the Georgia Institute of Technology have evaluated the thermal and fluids performance of the helium-cooled modular divertor with multiple jets (HEMJ) over the past decade. This finger-type divertor was studied both experimentally at nearly prototypical conditions and numerically at fully prototypical operating conditions using experimentally validated simulations. Recently, supercritical carbon dioxide (sCO 2 ) has been studied as the primary coolant in power cycles and other applications in various systems, in part because CO 2 achieves the high densities typical of supercritical fluids at relatively low temperatures and pressures, with a critical point of (7.38 MPa, 31°C). This density makes it possible to realize very compact and efficient sCO 2 power cycles. The feasibility of sCO 2 as a coolant for plasma-facing components, specifically the divertor, was therefore evaluated as part of the Fusion Energy System Studies design study activities. This work compares the thermal-fluid performance of helium and sCO 2 in the HEMJ divertor geometry using numerical simulations at prototypical conditions: inlet temperatures T i = 600°C to 700°C, pressures p ≈ 10 MPa, and steady-state incident heat fluxes on the tile q″ < 17 MW/m 2 . The performance is quantified here as the maximum heat flux that can be accommodated by the plasma-facing tile, the pumping power fraction, defined as the ratio of the coolant pumping power to the incident thermal power, and the operating stress limits based on ASME pressure vessel criteria. As expected, helium requires lower mass flow rates and pumping power fractions within imposed maximum temperature limits for the HEMJ pressure boundary. However, it also appears that neither helium nor sCO 2 can remove 10 MW/m 2 of incident heat flux while meeting ASME pressure vessel criteria. Finally, the numerical modeling reveals that sCO 2 may remove slightly higher incident heat fluxes than helium due to the imposed stress limits due to the sCO 2 coolant resulting in smaller local temperature gradients, albeit at a considerably higher pumping power fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

I-mode pedestal relaxation events in the Alcator C-Mod and ASDEX Upgrade tokamaks

Abstract In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light on their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks — Alcator C-Mod and ASDEX Upgrade — is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, a growing oscillating precursor before the PRE onset is observed in the region close to the separatrix in both devices, and a discussion on a possible triggering mechanism is outlined. The PRE relative energy loss from the confined region is found to increase with decreasing pedestal top collisionality ν ped * . Similarly, also the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing ν ped * . Based on these relations, the PRE relative energy loss in future devices such as DEMO and ARC is estimated. Finally, the divertor peak energy fluence due to the PRE is measured on each device. Those values are then compared to the model introduced in Eich et al (2017 Nucl. Mater. Energy 12 84–90) for type-I edge localized modes. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections of the PRE divertor energy fluence to DEMO and ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗