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 145 records · Page 8

Nonideal Fields Solve the Injection Problem in Relativistic Reconnection

Magnetic reconnection in relativistic plasmas is well established as a fast and efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. With particle-in-cell simulations, we demonstrate the importance of nonideal fields for the early stages (“injection”) of particle acceleration. Most of the particles ending up with high energies (near or above the mean magnetic energy per particle) must have passed through nonideal regions where the assumptions of ideal magnetohydrodynamics are broken (i.e., regions with E > B or nonzero E ∥ = E ∙ B/B), whereas most of the particles that do not experience nonideal fields end up with Lorentz factors of order unity. So, injection by nonideal fields is a necessary prerequisite for further acceleration. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-Dimensional Analysis of Conical Pulsed Inductive Plasma Thruster Performance

A model of the maximum achievable exhaust velocity of a conical theta pinch pulsed inductive thruster is presented. A semi-empirical formula relating coil inductance to both axial and radial current sheet location is developed and incorporated into a circuit model coupled to a momentum equation to evaluate the effect of coil geometry on the axial directed kinetic energy of the exhaust. Inductance measurements as a function of the axial and radial displacement of simulated current sheets from four coils of different geometries are t to a two-dimensional expression to allow the calculation of the Lorentz force at any relevant averaged current sheet location. This relation for two-dimensional inductance, along with an estimate of the maximum possible change in gas-dynamic pressure as the current sheet accelerates into downstream propellant, enables the expansion of a one-dimensional circuit model to two dimensions. The results of this two-dimensional model indicate that radial current sheet motion acts to rapidly decouple the current sheet from the driving coil, leading to losses in axial kinetic energy 10-50 times larger than estimations of the maximum available energy in the compressed propellant. The decreased available energy in the compressed propellant as compared to that of other inductive plasma propulsion concepts suggests that a recovery in the directed axial kinetic energy of the exhaust is unlikely, and that radial compression of the current sheet leads to a loss in exhaust velocity for the operating conditions considered here.

Hallock, A. K.↗

Combined Normal and Disruption, Electromagnetic Transient, Thermal, and Structural Analysis of COMPASS Upgrade

We report that COMPASS Upgrade is a replacement for the COMPASS device at the Institute of Plasma Physics (IPP) of the Czech Academy of Sciences. It is a cryogenic copper machine with a major radius of 0.89 m, with 5 T at the plasma centerline and 2-MA plasma current. It is intended to develop ITER relevant plasma densities and high power fluxes in the divertor region. The entire vacuum vessel is planned to reach 500 °C. The COMPASS-U TF coil is a bitter plate-like design with radial conductor plates that extend the width of the TF inner leg. Currents redistribute radially based on resistive and inductive effects. This necessitates an electromagnetic (EM) transient analysis of the coil current distribution. For copper coils like COMPASS Upgrade, the solution must be coupled with a thermal solution to properly include the temperature-dependent resistive effects. The EM transient solution includes the TF magnetic loads, and a stress pass on the results produces in-plane loads. The EM solution chosen requires inclusion of all the conducting and nonconducting regions surrounding the TF coil. With the TF current distribution solved, the addition of the poloidal coils and resulting background fields allows the determination of the out-of-plane loads. Integration of the Joule heating produced from the EM solution produces temperature distributions throughout the shot, which can be read in to a structural pass to include thermal stresses in the evaluation. This is true of both TF and PF temperature prediction with appropriate packing fractions. Currents in the PF coils, with computed background fields, produce Lorentz loads. Time transients of the PF currents produce startup eddy currents in the structures included in the model; inclusion of the vessel and other passive structures allows these to be included in a stress pass. The EM transient solution includes all the elements needed for a disruption analysis with prescribed motions and quench of the plasma. From a model originally intended to simulate current diffusion in the TF coil, the analysis can be refined to be used to predict the performance of most of the components of the tokamak or provide boundary conditions for more detailed submodels. In this article, the use of the current diffusion model for bounding assessments of the COMPASS Upgrade TF, PF, vessel, external structure, and disruption analyses will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of HTS Twisted Stacked-Tape Cable (TSTC) Conductor for High-Field, High-Current Fusion Magnets

The critical currents of a twisted single-tape were measured at 4.2 K in fields up to 17 T. It was confirmed that the critical current of a twisted-tape was similar to that of a flat tape in c-axis fields. Based on the single-tape critical-currents, a 40-tape (4-mm width) twisted stacked-tape cable (TSTC) conductor was experimentally evaluated at fields up to 17 T at 4.2 K. It was found that the Lorentz load degradation of the TSTC conductor was negligible at the end of a cyclic test. The TSTC cabling method has been considered to be suitable for developing a high-field high-current REBCO conductor for magnet applications as well as power transmission cables. A method to fabricate the inner legs of a D-shape toroidal field coil using a TSTC conductor has been discussed. This method allows the conductor to properly resist the transverse Lorentz load and mitigate ac loss issues by adopting a parallel-HTS-tape flat cable configuration where necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Four no-go theorems on the existence of spin and orbital angular momentum of massless bosons

The past decades have seen substantial interest in the so-called orbital angular momentum (OAM) of light, driven largely by its diverse range of applications. However, there are fundamental theoretical issues with decomposing the angular momentum of massless particles, such as photons, into spin (SAM) and orbital angular momentum parts. While the angular momentum of massive particles has a natural splitting into the Wigner SAM and OAM, there are numerous proposed splittings for photons and no consensus about which is correct. Moreover, it has been shown that most of the proposed SAM and OAM operators do not satisfy the defining commutation relations of angular momentum operators and are thus not legitimate splittings. Here, we prove that it is generally impossible to split the total angular momentum operator of massless bosons, such as photons and gravitons, into spin and orbital parts. We prove two further generalizations of this result, showing that there are no SAM-OAM splittings even if (1) the SAM operator generates non-internal symmetries or (2) if one allows the SAM and OAM operators to generate non-SO(3) symmetries.

Chern numbers↗

A Plasma Instability Theory of GRB Emission

A new theory for gamma-ray burst radiation is presented. In this theory, magnetic fields and relativistic electrons are created through plasma processes arising as a relativistic shell passes through the interstellar medium. The gamma-rays are produced through synchrotron self-Compton emission. It is found that shocks do not arise in this theory, and that efficient gamma-ray emission only occurs for a high Lorentz factor and a high-density interstellar medium. The former explains the absence of gamma-ray bursts with thermal spectra. The latter provides the Compton attenuation theory with an explanation of why the interstellar medium density is always high. The theory predicts the existence of a class of extragalactic optical transient that emit no gamma-rays.

Brainerd, Jerome J.↗

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗

Spin-dependent interactions and heavy-quark transport in the quark-gluon plasma

We extend a previously constructed T-matrix approach to the quark-gluon plasma (QGP) to include the effects of spin-dependent interactions between partons. Following earlier work within the relativistic quark model, the spin-dependent interactions figure as relativistic corrections to the Cornell potential. When applied to the vacuum spectroscopy of quarkonia, in particular their mass splittings in S- and P-wave states, the issue of the Lorentz structure of the confining potential arises. We confirm that a significant admixture of a vector interaction (to the previously assumed scalar interaction) improves the description of the experimental mass splittings. The temperature corrections to the in-medium potential are constrained by results from thermal lattice-QCD for the equation of state (EoS) and heavy-quark (HQ) free energy in a self-consistent set-up for heavy- and light-parton spectral functions in the QGP. We then deploy the refined in-medium heavy-light T-matrix to compute the charm-quark transport coefficients in the QGP. The vector component of the confining potential, through its relativistic corrections, enhances the friction coefficient for charm quarks in the QGP over previous calculations by tens of percent at low momenta and temperatures, and more at higher momenta. Our results are promising for improving the current phenomenology of open heavy-flavor observables at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Static current-sheet models of quiescent prominences

A particular class of theoretical models idealize the prominence to be a discrete flat electric-current sheet suspended vertically in a potential magnetic field. The weight of the prominence is supported by the Lorentz force in the current sheet. These models can be extended to have curved electric-current sheets and to vary three-dimensionally. The equation for force balance is 1 over 4 pi (del times B) times Bdel p- p9 z=zero. Using Cartesian coordinates we take, for simplicity, a uniform gravity with constant acceleration g in the direction -z. If we are interested not in the detailed internal structure of the prominence, but in the global magnetic configuration around the prominence, we may take prominence plasma to be cold. Consideration is given to how such equilibrium states can be constructed. To simplify the mathematical problem, suppose there is no electric current in the atmosphere except for the discrete currents in the cold prominence sheet. Let us take the plane z =0 to be the base of the atmosphere and restrict our attention to the domain z greater than 0. The task we have is to solve for a magnetic field which is everywhere potential except on some free surface S, subject to suit able to boundary conditions. The surface S is determined by requiring that it possesses a discrete electric current density such that the Lorentz force on it is everywhere vertically upward to balance the weight of the material m(S). Since the magnetic field is potential in the external atmosphere, the latter is decoupled from the magnetic field and its plane parallel hydrostatic pressure and density can be prescribed.

Wu, F.↗

The motional Stark effect diagnostic for ITER

An overview of the plans for the motional Stark effect (MSE) diagnostic installation on the International Thermonuclear Experimental Reactor (ITER) is presented. The MSE diagnostic uniquely provides spatially localized magnetic field measurements inside the plasma. These are used to constrain equilibrium reconstructions to determine q(r), the safety factor as a function of minor radius. Meeting the system requirements to deliver q-profiles and related quantities with the specified radial resolution of 20 points over the minor radius, 10 ms time resolution, and better than 10% accuracy is challenging. MSE systems observe the D/H-α emission near 656.3 nm from neutral beams. As the beam atoms traverse the magnetic field, B⃗, at high velocity, v⃗, they experience a Lorentz electric field, v⃗×B⃗, which causes the spectral emission to be split and polarized due to the Stark effect. Traditional MSE-LP (line polarization) measurements determine the direction of the magnetic field in the observation volume using polarimetric analysis of the detected light. The harsh conditions of ITER are expected to deposit thin films of contaminants on the first mirror, which would alter the polarization state of reflected light significantly. On ITER, the combination of high magnetic field strength and high energy beams makes the Stark spectrum resolution suitable for the determination of the magnetic field magnitude from the line shift, so this approach has been selected. Every aspect of the measurement system must be planned for the burning plasma environment and carefully analyzed ahead of time. Current status and plans for the system are presented.

Instruments & Instrumentation↗

Zero curvature is a necessary and sufficient condition for a spin-orbital decomposition

There has been an extended debate regarding the existence of a spin-orbital decomposition of the angular momentum of photons and other massless particles. It was recently shown that there are both geometric and topological obstructions preventing any such decomposition. Here we show that any geometric connection on a particle’s state space induces a splitting of the angular momentum into two operators. These operators are well-defined angular momentum operators if and only if the connection has zero curvature. Massive particles have two canonical curved connections corresponding to boosts and rotations, respectively. Furthermore, these can be uniquely combined to produce a flat connection, and this gives a novel derivation of the Newton-Wigner position operator and the corresponding spin and orbital angular momenta for relativistic massive particles. When the mass is taken to zero, transverse boosts and rotations degenerate, leaving only a single connection for massless particles. This connection produces a commonly proposed splitting of the massless angular momentum into two operators. However, the connection is not flat, explaining why these operators do not satisfy the angular momentum commutation relations and are thus not true spin and orbital angular momentum operators.

Angular momentum↗

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding Large-Scale Dynamos in Unstratified Rotating Shear Flows

For this work, we combine simulations with new analyses that overcome previous pitfalls to explicate how nonhelical mean-field dynamos grow and saturate in unstratified, magnetorotationally driven turbulence. Shear of the mean radial magnetic field amplifies the azimuthal component. Radial fields are regenerated by velocity fluctuations that induce shear of radial magnetic fluctuations, followed by Lorentz and Coriolis forces that source a negative off-diagonal component in the turbulent diffusivity tensor. We present a simple schematic to illustrate this dynamo growth. A different part of the Lorentz force forms a third-order correlator in the mean electromotive force that saturates the dynamo.

accretion disk & black-hole plasma↗

Radiation from Shock-Accelerated Particles

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration, generation of magnetic fields , and associated radiation. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. The shock structure depends on the composition of the jet and ambient plasma (electron-positron or electron-ions). Strong electromagnetic fields are generated in the reverse , jet shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. We have calculated, self-consistently, the radiation from electrons accelerated in the turbulent magnetic fields. We found that the synthetic spectra depend on the Lorentz factor of the jet, its thermal temperature and strength of the generated magnetic fields. The detailed properties of the radiation are important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jet shocks, and supernova remnants

Nishikawa, Ken-ichi↗

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗

Fast magnetization of a high-to-low-beta plasma beam

The magnetization of a high-beta (plasma energy density/magnetic-field energy density) hydrogen-plasma beam injected into a vacuum transverse magnetic field is studied experimentally. Nominal parameters were Ti = 1 eV, Te = 5 eV, n = 3 x 10 to the 13th/cu cm or less, v(i) = 7 x 10 to the 6th cm/sec or less, t(pulse) less than 70 microsec, and Bz = 300 G or less. Plasma characteristics were measured for a wide beam and a downstream distance, x = 300 rho(i) or less, where x is the downstream distance and rho(i) is the ion gyroradius. A brief initial state of diamagnetic propagation is observed, followed by magnetized propagation accompanied by beam compression transverse to B with as much as a factor of 4 increase in density and a slight drift of the beam in the ion Lorentz force direction.

Song, J. J.↗

Pulsed Fission-Fusion (PuFF)

In September 2013 the NASA Innovative Advanced Concept (NIAC) organization awarded a phase I contract to the PuFF team. Our phase 1 proposal discussed a pulsed fission-fusion propulsion system that injected gaseous deuterium (D) and tritium (T) as a mixture in a column, surrounded concentrically by gaseous uranium fluoride (UF6) and then an outer shell of liquid lithium. A high power current would flow down the liquid lithium and the resulting Lorentz force would compress the column by roughly a factor of 10. The compressed column would reach criticality and a combination of fission and fusion reactions would occur. The fission reactions would further energize the fusion center, and the fusion reactions would generate neutrons that promote more complete burnup of the fission fuel. The lithium liner provides some help as a neutron reflector but also acts as a propulsive medium, being converted to plasma which is then expanded against a magnetic nozzle for thrust. The expansion of the (primarily) lithium plasma against the nozzle's magnetic field inducts a current that is used to charge the system for the next pulse. Our concept also included secondary injection of a Field Reversed Configuration (FRC) plasmoid that would provide a secondary compression direction, axially against the column, and push the column away from the injection manifold, increasing the manifold's survivability.Our phase 1 proposal included modeling the above process first under steady state assumptions and second under a time variant integration. We proposed including these results into a Mars concept vehicle and finally proposing promising conditions to be evaluated experimentally in Phase II. In phase I we quickly realized that we needed to modify our approach. Our steady state work was completed as proposed, and the results indicated that one, a two stage compression system was not needed and two, that we wanted to move away from UF6. The steady state model shows much more margin than expected, to the point that we may well reach breakeven with the Charger – 1 facility, a 572 kJ Marx bank currently under refurbishment at UAH. Additionally we found that using gaseous D-T and UF6, provided a relatively simple prospect of using a pulsed injector, made reaching criticality more difficult. The introduction of large amounts of fluorine meant a radiative sink, sapping power from the fusion plasma and was harder to handle. Therefore we moved to a solid uranium target that held D-T under pressure. In so doing we could move our target closer to criticality and remove any material that did not sustain the reaction.

Fusion↗

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗