Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lorentz 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

PETN spark-gap detonators

The well-developed theory of Lorentz plasma that is dominated by electron–ion interactions is used to calculate the PETN arc characteristics. The spark-gap discharge current is a ramp with 10 to 25 ns rise time to peak and remaining constant subsequently. The approximate formulas for the arc channel conductivity, arc temperature, arc radius, and shock pressure from the arc are obtained from a system of nonlinear ordinary differential equations, which is the similarity solution of hydrodynamic equations similar to the Braginskii approximation. These arc parameters are given for the peak current ranging from 100 A to 1000 A and with different rise times. Representative cases are compared to the nonlinear ordinary differential equation code results. The shock pressures at the peak current are comparable to those from a typical commercial EBW bridgewire burst reported in the literature; the arc radius at the peak current is comparable to a typical bridgewire diameter of 0.0375 mm (e.g., RISI detonators, RP-1, and RP-80). The relevant Pop-Plot for low-density PETN is converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. Finally, this criterion is then used to determine the detonation thresholds, which are comparable with test data obtained by Tucker, et al.

42 ENGINEERING↗
DOE OSTI2021

Confinement of relativistic electrons in a magnetic mirror en route to a magnetized relativistic pair plasma

Creating a magnetized relativistic pair plasma in the laboratory would enable the exploration of unique plasma physics relevant to some of the most energetic events in the universe. As a step toward a laboratory pair plasma, we have demonstrated an effective confinement of multi-MeV electrons inside a pulsed-power-driven 13 T magnetic mirror field with a mirror ratio of 2.6. The confinement is diagnosed by measuring the axial and radial losses with magnetic spectrometers. The loss spectra are consistent with ≤2.5 MeV electrons confined in the mirror for ∼1 ns. With a source of 1012 electron-positron pairs at comparable energies, this magnetic mirror would confine a relativistic pair plasma with Lorentz factor γ∼6 and magnetization σ∼40.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Relativistic Alfvén Waves Entering Charge-starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A /j 0 $\ll$ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Friction in a strongly magnetized neutral plasma

A charged projectile generates an electrostatic wake as it moves through a background plasma. This wake normally gives rise to a drag force that acts antiparallel to the velocity vector of the projectile, and is commonly referred to as the stopping power. A recent theory predicted that an additional component of the friction force that is perpendicular to both the projectile velocity and the Lorentz force arises when a plasma is strongly magnetized. We extend the previous analysis of this transverse friction force, which was based on the one-component plasma model, to treat a two-component charge-neutral electron-ion plasma. The direction and magnitude of the transverse force is found to depend on the projectile speed relative to the electron and ion thermal velocities, and to change sign three times. For projectile speeds below a certain threshold, the transverse force is dominated by collisions with ions, while above this threshold, only electron collisions contribute. Here, the average trajectory of the projectile is significantly altered by the transverse force as it slows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A research program to measure the lifetime of spin polarized fuel

The use of spin polarized fuel could increase the deuterium-tritium (D-T) fusion cross section by a factor of 1.5 and, owing to alpha heating, increase the fusion power by an even larger factor. Issues associated with the use of polarized fuel in a reactor are identified. Theoretically, nuclei remain polarized in a hot fusion plasma. The similarity between the Lorentz force law and the Bloch equations suggests polarization can be preserved despite the rich electromagnetic spectrum present in a magnetic fusion device. The most important depolarization mechanisms can be tested in existing devices. The use of polarized deuterium and 3 He in an experiment avoids the complexities of handling tritium, while encompassing the same nuclear reaction spin-physics, making it a useful proxy to study issues associated with full D-T implementation. 3 He fuel with 65% polarization can be prepared by permeating optically-pumped 3 He into a shell pellet. Dynamically polarized 7 Li-D pellets can achieve 70% vector polarization for the deuterium. Cryogenically-frozen pellets can be injected into fusion facilities by special injectors that minimize depolarizing field gradients. Alternatively, polarized nuclei could be injected as a neutral beam. Once injected, the lifetime of the polarized fuel is monitored through measurements of escaping charged fusion products. Multiple experimental scenarios to measure the polarization lifetime in the DIII-D tokamak and other magnetic-confinement facilities are discussed, followed by outstanding issues that warrant further study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗
DOE OSTI2021

Relativistic flying laser focus by a laser-produced parabolic plasma mirror

Here, the question of electromagnetic field intensification towards the values typical for strong field quantum electrodynamics is of fundamental importance. One of the most promising intensification schemes is based on the relativistic-flying mirror concept, which shows that the electromagnetic radiation reflected by the mirror will be frequency upshifted by a factor of 4γ 2 (γ is the Lorentz factor of the mirror). In laser-plasma interactions, such a mirror travels with relativistic velocities through plasma and typically has a parabolic form, which is advantageous for light intensification. Thus, a relativistic-flying parabolic mirror reflects the counterpropagating radiation in the form of a focused and flying electromagnetic wave with a high frequency. The relativistic-flying motion of the laser focus makes the electric and magnetic field distributions of the focus complicated, and the mathematical expressions describing the field distributions of the focus become of fundamental interest. We present analytical expressions describing the field distribution formed by an ideal flying mirror which has a perfect reflectance over the entire surface and wavelength range. The peak field strength of an incident laser pulse with a center wavelength of λ 0 and an effective beam radius of w e is enhanced by a factor proportional to γ 3 (w e /λ 0 ) in the relativistic limit. Electron-positron pair production is investigated in the context of invariant fields based on the enhanced electromagnetic field. The pair production rate under the relativistic-flying laser focus is modified by the Lorentz γ-factor and the beam radius-wavelength ratio (w e /λ 0 ). We show that the electron-positron pairs can be created by colliding two counterpropagating relativistic-flying laser focuses in vacuum, each of which is formed when a 180 TW laser pulse is reflected by a relativistic-flying parabolic mirror with γ=12.2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Estimate of Lorentz force effects on charged particles

Thermonuclear (TN) fusion reaction between deutrons are D + D → T + P D + D → 3 He + n where P is the proton, D is the deuteron, T is the triton, n is the neutron, and 3 He is the helion. They are born with high energies (0.82 MeV for 3 He, 1.01 MeV for T, 3.02 MeV for P, and 2.45 MeV for n). The charged particles slow down due to the drag force (stopping power). While slowing down, the charged particles deposit their energy, which may affect the TN burn characteristics. When magnetic field is present, the velocity of the charged particle would be affected due to the Lorentz force, which can affect the stopping of the charged particle. As a consequence, TN burn characteristics may change, Our purpose is to asses the importance of the Lorentz force in stopping of the charged particle, based on a simple comparison between the Lorentz force and stopping power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Collective plasma effects of electron–positron pairs in beam-driven QED cascades

Understanding the interplay of strong-field QED and collective plasma effects is important for explaining extreme astrophysical environments like magnetars. It has been shown that QED pair plasma can be produced and observed by passing a relativistic electron beam through an intense laser field. This work presents in detail multiple sets of 3D QED-particle-in-cell simulations to show the creation of pair plasma in the QED cascade. The beam driven method enables a high pair particle density and also a low particle Lorentz factor, which both play equal roles on exhibiting large collective plasma effects. Finite laser frequency upshift is observed with both ideal parameters (24 PW laser colliding with a 300 GeV electron beam) and with existing technologies (3 PW laser colliding with a 30 GeV electron beam).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically dominated (${\sigma }_{0}\gg 1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim {\sigma }_{0}{\gamma }_{\mathrm{th}0}$ (here, ${\gamma }_{\mathrm{th}0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are preferentially oriented in the plane perpendicular to the field. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu {F}_{\nu }\propto {\nu }^{s}$ that is much harder than for isotropic particles. Remarkably, for ${\sigma }_{0}\gg 1$ we find that the angle-integrated spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta {B}_{\mathrm{rms}0}^{2}/{B}_{0}^{2}\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Further, our findings may help explain the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of pulsar wind nebulae.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of Coulomb coupling on friction in strongly magnetized plasmas

The friction force on a test particle traveling through a plasma that is both strongly coupled and strongly magnetized is studied using molecular dynamics simulations. In addition to the usual stopping power component aligned antiparallel to the velocity, a transverse component that is perpendicular to both the velocity and Lorentz force is observed. This component, which was recently discovered in weakly coupled plasmas, is found to increase in both absolute and relative magnitude in the strongly coupled regime. Strong coupling is also observed to induce a third component of the friction force in the direction of the Lorentz force. These first-principles simulations reveal novel physics associated with collisions in strongly coupled, strongly magnetized, plasmas that are not predicted by existing kinetic theories. The effect is expected to influence macroscopic transport in a number of laboratory experiments and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
DOE OSTI2022

An overview of magneto-inertial fusion on the Z Machine at Sandia National Laboratories

We present an overview of the magneto-inertial fusion (MIF) concept MagLIF (Magnetized Liner Inertial Fusion) pursued at Sandia National Laboratories and review some of the most prominent results since the initial experiments in 2013. In MagLIF, a centimeter-scale beryllium tube or "liner" is filled with a fusion fuel, axially pre-magnetized, laser pre-heated, and finally imploded using up to 20 MA from the Z machine. All of these elements are necessary to generate a thermonuclear plasma: laser preheating raises the initial temperature of the fuel, the electrical current implodes the liner and quasi-adiabatically compresses the fuel via the Lorentz force, and the axial magnetic field limits thermal conduction from the hot plasma to the cold liner walls during the implosion. MagLIF is the first MIF concept to demonstrate fusion relevant temperatures, significant fusion production (>10^13 primary DD neutron yield), and magnetic trapping of charged fusion particles. On a 60 MA next-generation pulsed-power machine, two-dimensional simulations suggest that MagLIF has the potential to generate multi-MJ yields with significant self-heating, a long-term goal of the US Stockpile Stewardship Program. At currents exceeding 65 MA, the high gains required for fusion energy could be achievable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Motional Stark Effect Diagnostic for NSTX-U

This report describes the development of the motional Stark Effect diagnostic for NSTX. The MSE concept relies upon observation of the Balmer-alpha (n=3 to n=2) emission from a neutral hydrogen beam traversing a plasma. As a beam at high velocity, $\vec{ν}$, passes through the magnetic field in the plasma, $\vec{B}$, it experiences in its reference frame a Lorentz electric field, $\vec{E}$ = $\vec{ν}$ x $\vec{B}$. This electric field causes the spectral emission to be split and polarized as described by the Stark effect. The transition is separated into a set of lines in which the line spacing is linearly proportional to the electric field magnitude, and the line polarization is dependent on the electric field direction. Traditionally, MSE has used polarimetry to determine the polarization angle of the light, and related that to the magnetic field pitch angle in the plasma. The MSE diagnostic was developed under high magnetic field ($>$ 1~T) conditions, and the lower fields of recent high-beta experiments such as NSTX posed challenges to the technique. At low magnetic fields Doppler broadening effects mask the Stark effect. With the use of a narrow bandwidth filter MSE measurements at very low magnetic fields have been obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma facing components with capillary porous system and liquid metal coolant flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 °C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently, a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [A. Khodak and R. Maingi, Nucl. Mater. Energy 26, 100935 (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In the case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a fusion nuclear science facility [Kessel et al., Fusion Sci. Technol. 75, 886 (2019)]: 10 T toroidal field and 10 MW/m 2 peak incident heat flux. Finally, computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma Facing Components with Capillary Porous System and Liquid Metal Coolant Flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [Khodak and Maingi (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile, as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a Fusion Nuclear Science Facility [Kessel et al. (2019)]: 10T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

High Heat Flux↗
DOE OSTI2021

Design and Analyses of the NSTX-U PF1A Poloidal Field Coil Support

New poloidal field (PF) coils and supports have been designed for the National Spherical Torus Experiment Upgrade (NSTX-U) Recovery project. The PF1A coil structural requirements demand a large preload to pre-compress the coil insulation. Coil support design accounts for all operational conditions and loads of NSTX-U including bakeout, seismic events, and electromagnetic loads. Background magnetic fields produce coil Lorentz forces that react on the structural supports. These Lorentz forces include interactions from the coil current and all plasma events: normal operations, disruptions, and vertical displacement events. This paper will detail the design and design loads and highlight some of the analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗