Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetized plasma”

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 919 records · Page 51

Observation of Fast Current Redistribution in an Imploding Plasma Column

Spectroscopic measurements of the magnetic field evolution in a Z-pinch throughout stagnation and with particularly high spatial resolution reveal a sudden current redistribution from the stagnating plasma (SP) to a low-density plasma (LDP) at larger radii, while the SP continues to implode. Based on the plasma parameters it is shown that the current is transferred to an increasing-conductance LDP outside the stagnation, a process likely to be induced by the large impedance of the SP. In conclusion, since an LDP often exists around imploding plasmas and in various pulsed-power systems, such a fast current redistribution may dramatically affect the behavior and achievable parameters in these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing the plasma-induced thermal loads on a 200 kW light-ion helicon plasma source via infra-red thermography

The light-ion helicon plasma source of the Proto-MPEX linear plasma device has been recently upgraded to enable pulsed (0.5–1 s) operation up to 200 kW. The main objective of this work is to report on the plasma-induced surface heat fluxes incident on the helicon window during high power operation (60–150 kW net power) for the purpose of the design of the upcoming material plasma exposure eXperiment (MPEX). The IR imaging system and associated physics models for the extraction of surface heat fluxes are presented. Furthermore, experimental results demonstrate that the control of the plasma strike point via magnetic flux mapping and the use of dedicated limiters is effective at reducing the heat loads on the dielectric window. Moreover, it is found that the flux mapping must create a gap between the plasma and the dielectric window of at least the plasma radial decay length. Extrapolated to 200 kW net RF power, this mode of operation can reduce the power lost to the dielectric window by 33%. The results presented have direct application to the design of high-powered helicon plasma sources that operate in a steady state. Examples include linear divertor simulators such as MPEX, electric thrusters, negative ion sources for NBI and linear plasma–material interaction test stands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Further rotation reversal studies in C-Mod L-mode plasmas

Studies of core toroidal rotation reversal phenomenology in C-Mod deuterium L-mode plasmas have been expanded to include details of the dependences on plasma current and toroidal magnetic field. Rotation reversal occurs at a critical density, and universal scaling indicates that the product of ncritq95R ∼ BT/2, with ncrit in 1020/m3, R in m, and BT in T. Measurements in H and He plasmas exhibit similar behavior, including a connection with the linear Ohmic confinement/saturated Ohmic confinement transition and the cutoff for non-diffusive heat transport. Electron density and ion cyclotron range of frequencies power modulation experiments suggest that the collisionality ν* is a unifying parameter. Strong impurity puffing causes the critical density to increase, indicating that the situation is more complicated than only collisionality, perhaps involving the details of the effects of dilution on ion temperature gradient mode stability.

Physics↗

Magnetic Field Amplification by a Nonlinear Electron Streaming Instability

Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental inference of flux tunneling between magnetic island chains in tokamaks

We report confinement bifurcations due to flux tunneling in magnetic islands in toroidal plasmas. This occurs in 2/1 islands, triggered by decoupling of non-overlapping adjacent 5/2 islands. While subject to electron cyclotron current drive (ECCD), the temperature at the 2/1 O-point is flat before but peaks after decoupling. Here, this is shown to be due to a bifurcation from stochastic to nested magnetic structure, caused by the removal of flux tunneling through intersecting manifolds of overlapping heteroclinic tangles. 2/1 island stability correlates with coupling events in an ITER baseline scenario plasma where ECCD can not stabilize the coupled islands, showing the critical impact of flux tunneling on disruption avoidance in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Consequences of Rayleigh-Taylor Instabilities on Implosion Dynamics in Gas-puff Z-pinch Experiments (Final Report)

The Laboratory of Plasma Studies (LPS), Cornell University, proposed to carry out a comprehensive series of gas-puff z-pinch experiments on a 1 MA pulsed power machine for the purpose of understanding the dynamics of imploding magnetized high energy density plasmas in the presence of Rayleigh-Taylor Instabilities (RTI). Cornell’s pulsed power facilities, including the associated diagnostic instruments, are largely supported under the auspices of the Center for Pulsed-power-driven High-energy-density Plasmas (the Center), an NNSA-sponsored Center of Excellence. The proposed experiments were carried out predominantly by LPS graduate student, Sophia Rocco (the “primary grad student”), in collaboration with other graduate students and the scientific and technical staff of the Center using Cornell’s 1 MA COBRA pulsed power facility and all of its associated diagnostic instruments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strong suppression of heat conduction in a laboratory replica of galaxy-cluster turbulent plasmas

In conventional gases and plasmas, it is known that heat fluxes are proportional to temperature gradients, with collisions between particles mediating energy flow from hotter to colder regions and the coefficient of thermal conduction given by Spitzer’s theory. However, this theory breaks down in magnetized, turbulent, weakly collisional plasmas, although modifications are difficult to predict from first principles due to the complex, multiscale nature of the problem. Understanding heat transport is important in astrophysical plasmas such as those in galaxy clusters, where observed temperature profiles are explicable only in the presence of a strong suppression of heat conduction compared to Spitzer’s theory. To address this problem, we have created a replica of such a system in a laser laboratory experiment. Our data show a reduction of heat transport by two orders of magnitude or more, leading to large temperature variations on small spatial scales (as is seen in cluster plasmas).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Proton deflectometry with in situ x-ray reference for absolute measurement of electromagnetic fields in high-energy-density plasmas

In this work, we report a technique of proton deflectometry that uses a grid and an in situ reference x-ray grid image for precise measurements of magnetic fields in high-energy-density plasmas. A D 3 He fusion implosion provides a bright point source of both protons and x-rays, which is split into beamlets by a grid. The protons undergo deflections as they propagate through the plasma region of interest, whereas the x-rays travel along straight lines. The x-ray image, therefore, provides a zero-deflection reference image. The line-integrated magnetic fields are inferred from the shifts of beamlets between the deflected (proton) and reference (x-ray) images. We developed a system for analysis of these data, including automatic algorithms to find beamlet locations and to calculate their deflections from the reference image. The technique is verified in an experiment performed at OMEGA to measure a nonuniform magnetic field in vacuum and then applied to observe the interaction of an expanding plasma plume with the magnetic field.

47 OTHER INSTRUMENTATION↗

Scrape-off layer stability, turbulence and transport in MAST-U (Final Report)

U.S. Department of Energy (DOE) supported work at Lodestar Research Corporation under grant DE-SC0019270 is summarized in this final report. The research conducted under this grant explored the properties of the edge plasma in spherical torus magnetic fusion devices. Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The edge and scrape-off layer (SOL) region of fusion plasmas is the region where the hot core plasma transitions to a cooler surrounding boundary plasma. The properties of this boundary plasma determine how interactions with material surfaces will impact both the overall device performance and the lifetime of crucial components such as the walls of the device chamber and other hardware. The research carried out studied filamentary plasma instabilities and turbulent transport in the edge and SOL of the MAST U spherical torus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model

The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.

DIII-D↗

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗

Anticathode effect on multimodal azimuthal oscillations in electron beam generated E×B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (⁠E x B⁠) fields are promising for low-damage (gentle) material processing. However, these plasmas can be subject to the formation of plasma non-uniformities propagating in the ⁠E x B direction. These rotating plasma structures (or “spokes”) enhance the transport of charged species across the magnetic field, which can harm the gentle processing capability of the plasma. In this work, we investigate the role of electrostatically active boundaries on the spoke formation by incorporating a variable bias conducting boundary (known as an anticathode) placed on the axially opposite side of the cathode. Our findings indicate azimuthal mode suppression occurs when the anticathode is electron collecting. Furthermore, we show selective azimuthal mode suppression by biasing the anticathode to an intermediate potential between the cathode and anode potentials. In conclusion, these findings suggest a link between the axial electron confinement in the e-beam generated plasma and azimuthally propagating plasma structure formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamic Phase Transition in Magnetic Reconnection

By examining the entropy production in fully kinetic simulations of collisional plasmas, it is shown that the transition from collisional Sweet-Parker reconnection to collisionless Hall reconnection may be viewed as a thermodynamic phase transition. The phase transition occurs when the reconnection electric field satisfies E = $E_{D}\sqrt{m{e}/{m}_i}$ , where ${m{e}/{m}_i}$ is the electron-to-ion mass ratio and E D is the Dreicer electric field. This condition applies for all ${m{e}/{m}_i}$, including ${m{e}/{m}_i}$ = 1 , where the Hall regime vanishes and a direct phase transition from the collisional to the kinetic regime occurs. In the limit ${m{e}/{m}_i}$ → 0, this condition is equivalent to there being a critical electron temperature $Te ≈ m_iΩ^{2}_{i}δ^{2}$, where $Ω_i$ is the ion cyclotron frequency and δ is the current sheet half-thickness. The heat capacity of the current sheet changes discontinuously across the phase transition, and a critical power law is identified in an effective heat capacity. A model for the time-dependent evolution of an isolated current sheet in the collisional regime is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasilinear modeling of heat flux from microtearing turbulence

Microtearing (MT) turbulence affects plasma confinement and profile evolution in many different magnetic-confinement scenarios, ranging from tokamak core plasmas to the H-mode pedestal and from spherical tokamaks to reversed-field pinches. Thus, an urgent need exists to design and test reduced models of MT turbulence and transport. Here, a heuristic quasilinear model is composed and verified by comparing its predictions against nonlinear MT simulations. It is found to yield good predictions of magnetic flutter transport as key parameters are varied: the collision frequency, the electron temperature gradient, and the normalized plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗