Results for “Large Plasma Device” — Page 8

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

NASA NTRS1994

Temporal study of the current collected by a Langmuir probe biased to high voltages in the wake behind a conducting body

The temporal evolution of the voltage-current characteristic of a spherical Langmuir probe in the near-wake region of a conducting body is investigated in a pulsed plasma device. Measurements of the current collected by a sphere placed on axis behind a conducting disk are presented as a function of the bias voltage and time. The bias voltage on the collecting sphere is large enough, ranging from +500 to -3000 V, to attract ions or electrons across the wake region. In addition, plasma potential profiles as a function of sphere bias voltage and time throughout the nearwake and midwake region are presented. It is found that for a negative sphere bias a threshold voltage exists for current collection. This threshold voltage decreases during the formation of the wake behind the conducting disk. Potential measurements indicate that the sheath region around the negatively biased sphere contracts as the wake region establishes itself. This contraction of the sheath region increases the electric field attracting ions to the sphere and decreases the threshold voltage for current collection. For positive sphere bias voltage the current collected by the sphere does not exhibit a threshold voltage but is monotonically increasing.

Meassick, S.↗

DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗

DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Efficient facemask decontamination via forced ozone convection

The COVID-19 crisis has taken a significant toll on human life and the global economy since its start in early 2020. Healthcare professionals have been particularly vulnerable because of the unprecedented shortage of Facepiece Respirators (FPRs), which act as fundamental tools to protect the medical staff treating the coronavirus patients. In addition, many FPRs are designed to be disposable single-use devices, creating an issue related to the generation of large quantities of non-biodegradable waste. In this contribution, we describe a plasma-based decontamination technique designed to circumvent the shortages of FPRs and alleviate the environmental problems posed by waste generation. The system utilizes a Dielectric Barrier Discharge (DBD) to generate ozone and feed it through the fibers of the FPRs. The flow-through configuration is different than canonical ozone-based sterilization methods, in which the equipment is placed in a sealed ozone-containing enclosure without any flow through the mask polymer fibers. We demonstrate the rapid decontamination of surgical masks using Escherichia coli ( E. coli ) and Vesicular Stomatitis Virus (VSV) as model pathogens, with the flow-through configuration providing a drastic reduction in sterilization time compared to the canonical approach. We also demonstrate that there is no deterioration in mask structure or filtration efficiency resulting from sterilization. Finally, we show that this decontamination approach can be implemented using readily available tools, such as a plastic box, a glass tube, few 3D printed components, and the high-voltage power supply from a plasma globe toy. The prototype assembled for this study is portable and affordable, with effectiveness comparable to that of larger and more expensive equipment.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Simulation of Divertor Performance in ST40 Under Dynamic Double-Null Plasmas

A power fraction model was implemented for the simultaneous prediction of 3-D surface temperature evolution at all four divertor targets in near-double-null (DN) tokamak configurations, which is especially important for compact high-field devices that may not have the ability to dissipate large amounts of power on the high-field side. Evaluating the power-sharing between the four divertor strike points in a disconnected DN configuration is important for understanding the overall power balance, as well as for optimizing the power exhaust performance and prolonging the survivability of the plasma facing components (PFCs). This power-sharing is typically evaluated in terms of the separation between the primary and secondary separatrices at the outboard midplane, $\textit {dR}_{\text {sep}}$. The Heat flux Engineering Analysis Toolkit (HEAT) is coupled with Brunner’s power fraction model to simulate the deposited heat flux and resultant temperature change on 3-D divertor targets in a dynamic DN (DDN) pulse operation in ST40, a high-field spherical tokamak. The simulation results showed that with DDN operation, the operation time has significantly increased compared with single-null geometry configurations.

ST40↗

DOE OSTI2024

Deflection of Plasma by High Magnetic Fields at Variable Pressures (Final Report)

The goal of the proposed work is to study the deflection of plasma by high magnetic fields using the high uniform magnetic field of the Magnetized Dusty Plasma Experiment (MDPX) located at the Magnetized Plasma Research Laboratory (MPRL) at Auburn University. The proposal team consists of faculty from the University of Alabama in Huntsville (UAH). One of the motivations for this work is to study the concept of a magnetic nozzle for pulsed fusion propulsion. Fusion propulsion can provide extremely high specific impulse, >30,000 sec, for fuel-efficient and thus mass-efficient trips to other planets. One such fusion propulsion design of interest in this research is the Pulsed Fission-Fusion (PuFF) concept. PuFF uses a Z-pinch device to create a fusion reaction at the back of the spacecraft inside the magnetic nozzle. As the resulting plasma explosion is isotropic, a large portion of the momentum will be directed toward the ship and needs to be redirected to efficiently produce thrust. The magnetic nozzle is thus necessary to redirect or deflect the plasma from an isotropic expansion into the directed thrust. Thus, the pulsed-fusion propulsion nozzle is herein referred to as a deflection magnetic nozzle, as opposed to the more traditional and well-studied accelerating magnetic nozzle. The purpose of this project was thus to study the behavior of a plasma in the MDPX fields to better understand the plasma-field interactions and how it could lead to deflection and a pulsed fusion magnetic nozzle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Non-dimensional confinement scaling in similar negative triangularity plasmas on the DIII-D and TCV tokamaks

Similarity experiments were performed on the DIII-D and TCV tokamaks to explore the scaling of energy confinement in negative triangularity (NT) plasmas using non-dimensional variables. Near up-down symmetric plasmas with large top-bottom averaged NT were created in a lower single null configuration, with the shape of the separatrix being closely matched between the two devices. The normalized energy confinement is found to weakly improve at increasing collisionality and, between the two devices, shows a machine size scaling behavior between Bohm and gyro-Bohm. Engineering scaling on a large DIII-D dataset is in agreement with the non-dimensional experiment.

confinement scaling↗

NASA NTRS1983

Initial experimental results from a laboratory size beam plasma discharge device

A laboratory beam plasma discharge (BPD) device produced BPD in N2, A, and He. All features of the BPD observed in the device agree with those observed in a large vacuum chamber. The empirical ignition criteria determined in the large chamber apply in the small device but do not fit when used for extrapolation between the large and the small geometry. At some energies and magnetic fields beam currents exist for which the total light output in the BPD state varies by a factor of 2 with a factor of 6 pressure variation. Above 0.0001 torr the BPD width is pressure independent but for lower pressures it expands by as much as a factor of 4 at 0.00002 torr.

Konradi, A.↗

DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Turbulence in electrostatic ion-acoustic shocks

Three types of collisionless electrostatic ion-acoustic shocks are investigated using the University of California, Los Angeles, double plasma device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distribution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.

Means, R. W.↗

NASA NTRS1992

High-flux source of low-energy neutral beams using reflection of ions from metals

Reflection of low-energy ions from surfaces can be applied as a method of producing high-flux beams of low-energy neutral particles, and is an important effect in several areas of plasma technology, such as in the edge region of fusion devices. We have developed a beam source based on acceleration and reflection of ions from a magnetically confined coaxial RF plasma source. The beam provides a large enough flux to allow the energy distribution of the reflected neutrals to be measured despite the inefficiency of detection, by means of an electrostatic cylindrical mirror analyzer coupled with a quadrupole mass spectrometer. Energy distributions have been measured for oxygen, nitrogen, and inert gas ions incident with from 15 to 70 eV reflected from amorphous metal surfaces of several compositions. For ions of lighter atomic mass than the reflecting metal, reflected beams have peaked energy distributions; beams with the peak at 4-32 eV have been measured. The energy and mass dependences of the energy distributions as well as measurements of absolute flux, and angular distribution and divergence are reported. Applications of the neutral beams produced are described.

Cuthbertson, John W.↗

DOE OSTI2026

Predictions of m/n = 2/1 neoclassical tearing mode stabilization via high field side lower hybrid current drive on the DIII-D tokamak

Neoclassical tearing modes (NTMs) are a class of resistive instabilities that arise in tokamaks at rational surfaces and form magnetic islands. These islands flatten the pressure gradient, reducing plasma performance and may lead to disruptions if they grow large enough. Driving current within the island can stabilize the NTM, which has been achieved with electron cyclotron current drive (ECCD) on multiple devices. An alternative to ECCD is lower hybrid current drive (LHCD), which offers the advantages of increased current drive efficiency and reduced system cost. LHCD has been viewed as poorly suited for NTM suppression due to the large spatial extent of the driven current when in the multi-pass absorption regime (as has been the case in all past LHCD experiments). However, the driven current is more localized when in the single pass absorption regime, as is predicted for the DIII-D high field side (HFS) LHCD experiment. This work evaluates the feasibility of NTM suppression with HFS LHCD on DIII-D by predicting the island growth rate for a set of representative DIII-D plasmas via the modified Rutherford equation with and without the application of LHCD. In these plasmas, NTM suppression is achieved at reasonable power levels, even with finite misalignment between LH current and the island. The effect of current condensation was included and found to be most significant at smaller island sizes, assuming an experimentally typical temperature perturbation amplitude of 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Measuring the growth of Alfvén wave parametric decay instability using counter-propagating waves: Theory and simulations

The parametric decay instability (PDI) of Alfvén waves—where a pump Alfvén wave decays into a backward-propagating child Alfvén wave and a forward ion acoustic wave—is a fundamental nonlinear wave-wave interaction and holds significant implications for space and laboratory plasmas. However, to date there has been no direct experimental measurement of PDI. Here, we propose a novel and experimentally viable scheme to quantify the growth of Alfvén wave PDI on a linear device using a large pump Alfvén wave and a small counter-propagating seed Alfvén wave, with the seed-wave frequency tuned to match the backward Alfvén wave generated by standard PDI. Using hybrid simulations, we show that energy transfer from the pump to the seed reduces the latter's spatial damping. By comparing seed-wave amplitudes with and without the pump wave, this damping reduction can be used as a direct and reliable proxy for PDI growth. The method is validated in our simulations across a range of plasma and wave parameters and agrees well with theoretical predictions. Notably, the scheme exhibits no threshold for PDI excitation and is, in principle, readily implementable under current laboratory conditions. Finally, this scheme is a critical step toward solving the challenge of experimentally accessing Alfvén wave PDI and provides an elegant method that may be used to validate fundamental theories of parametric instabilities in controlled laboratory settings.

heliosphere↗

DOE OSTI2024

Glow Discharge Boronization and Real-Time Boronization Using an Impurity Powder Dropper in LHD

In the Large Helical Device (LHD), boronization is utilized for wall conditioning, employing two distinct methods: glow discharge using diborane gas and the injection of boron powder into the plasma via an impurity powder dropper (IPD). Glow discharge boronization is conducted after maintenance when the vacuum vessel is exposed to air and just before the start of an experimental campaign. In contrast, real-time boronization using the IPD is performed during discharges. The distribution of boron on plasma-facing components was assessed through long-term material probes in the case of glow discharge and through numerical simulations for the IPD method. Results indicated non-uniform boron distributions in both cases. The impact of these boronization techniques on impurity levels and the lifetime of their effects were also examined. Both methods effectively reduced oxygen and iron in the plasma, with the iron impurities originating from the first wall. For glow discharge boronization, the reduction in oxygen levels persisted until the end of the experimental campaign, whereas the reduction in iron levels only lasted a few days. In the case of boron injection during long pulse discharges exceeding 200 s, the suppression of oxygen and iron was observed. The effect on iron quickly diminished within several seconds after the injection ended, whereas the effect on oxygen lasted more than 100 s. The difference in the lifetimes of the boronization effects is further discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Turbulence in electrostatic ion acoustic shocks

Three types of collisionless electrostatic ion acoustic shocks are investigated using a double plasma (DP) device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distrubution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.

Means, R. W.↗