Search NASA⌕ Search

Engineering topics

Seyler, Charles E.

Publications and source records attributed to Seyler, Charles E..

Formulation of 8-moment plasma transport with application to the Nernst effect

8-moment plasma models using two different distribution functions are used to study the Nernst effect and heat transport in dense plasma. These models are presented in hyperbolic form in contrast to traditional parabolic systems derived from perturbing the distribution function, as in Braginskii [Rev. Plasma Phys. 1, 205 (1965)]. The hyperbolic moment formulation can be solved implicitly in time with straightforward and fast local solvers. The numerical implementation of 8-moment models with the relaxation method in the PERSEUS code is also presented. To test 8-moment PERSEUS compared to Braginskii’s transport equations, a verification test for the Nernst thermomagnetic wave by Velikovich et al. [Phys. Plasmas 26, 112702 (2019)] is performed that confirms the presence of the same physics, but with slight differences in the transport coefficients, which are tabulated in the limits of high and low magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Jet Formation Disruption From a Critical Applied Uniform Axial Magnetic Field

We examine the effects of varying the magnitude of an applied, uniform, axial magnetic field (Bz) on the formation of laboratory plasma jets produced by a 1-MA 100-ns rise time pulsed power generator in a radial foil configuration. A Helmholtz coil applies the external magnetic field (0 to 2 T). With a small-enough applied field, the foil surface ablation is relatively azimuthally uniform, and J×B forces drive the ablated plasma inward and upward to form a well-defined azimuthally symmetric jet. Using applied field strengths larger than a critical Bz = 1.1±0.1 T with aluminum foils, the plasma ablation is not azimuthally uniform, discrete bursts of plasma initiate from the foil surface, and the formation of a well-defined jet is disrupted. The critical Bz for this plasma jet formation disruption correlates with the foil material’s electrical resistivity and equation of state (EOS). To better understand the material-dependent phenomenon, we compare the experimental results with 3-D extended magnetohydrodynamics simulations of ablation of a 2 mm × 2 mm by 25-μm slab that represents a section of the radial foil. The simulations initialize the slab in the solid state and include the material resistivity and EOS from the solid to plasma phases. As is observed in the experimental disruption, these simulations show enhanced nonuniform plasma ablation with an applied Bz, which can inhibit the azimuthal uniformity necessary to produce a well-defined plasma jet. Furthermore, the simulations also show a material and resistivity dependence on the ablation process similar to the trend shown by the experiments. These results demonstrate the necessity of accurate, detailed modeling of material properties during the transition from solid to plasma phases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relativistic Modeling Capabilities in Perseus Extended MHD Simulation Code for HED Plasmas

This project is submitted in response to DE-FOA-0001153, and will be carried out by Dr. Charles Seyler and Dr. Nathaniel Hamlin as post doc under the direction of Professor Seyler. The objectives are to enable PERSEUS to self-consistently model relativistic high-energy-density (HED) laboratory and astrophysical phenomena. Electrons and ions in a plasma are generally modeled either as fluids using a two-fluid code, like PERSEUS, or as particles using a Particle-In- Cell (PIC) code. Nonrelativistic PERSEUS solves the two-fluid equations, formulated in terms of a generalized Ohm’s law (GOL), so as to model about nine orders of magnitude in density variation using a local semi-implicit method. Relativistic PERSEUS preserves this structure and therefore retains these advantageous properties, enabling it to model a broader range of relativistic HED phenomena than previous codes. We have also overcome a major technical challenge associated with solving a relativistic system of equations relating conserved quantities (momentum, energy, etc.) to primitive variables (density, velocity, and pressure). We have thus far developed a relativistic version of PERSEUS and demonstrated that it recovers expected nonrelativistic results. We are currently using relativistic PERSEUS to model two HED phenomena, namely laser-plasma interactions and X-pinches, both of which show spectroscopic evidence of relativistic electrons. In an X-pinch, a section of wire carrying sufficiently large current forms a plasma column that pinches down to a point-like source, generating X-rays for use in radiographic imaging. Laser-plasma interactions have been simulated using PIC codes, and we use their results as a benchmark for comparison. However, neither PIC codes nor existing two-fluid schemes can simulate X-pinches, due to the broad range of densities which have thus far made such a simulation computationally prohibitive. In both laser-plasma interaction and Xpinch simulations, we have observed relativistic phenomena that PERSEUS could not previously model. This includes relativistic channeling of a laser into an over-dense deuterium gas (i.e. induced transparency), consistent with PIC simulation results, along with energetic electrons and ions following an X-pinch and likely accompanied by strong X-ray emission, consistent with experimental observations of energetic X-rays. This proposal includes our latest simulation results. We will improve computational efficiency with techniques that allow the use of a less refined numerical grid. This will enable extensions to three-dimensional geometries. We will improve modeling of laser-plasma interactions and X-pinches, including radiative transport processes, in order to improve comparisons with PIC simulations and experimental diagnostics. We will also model power loss through electrode surface plasmas in magnetically insulated transmission lines (MITLs). We will model relativistic astrophysical phenomena in both astrophysical and HED regimes, which could shed light on the validity of laboratory astrophysics for certain problems. These phenomena are currently modeled using PIC codes and existing relativistic two-fluid codes, neither of which can model the range of length scales of which relativistic PERSEUS is capable. There is strong potential for this project to have a broad impact on both the astrophysical and laboratory plasma communities. In the laboratory community, our work is applicable to those concerned with laser-plasma interactions such as laser-driven fusion, and those in the dense Z pinch community concerned with X-pinches and electrode surface plasmas. An improved understanding of the mechanisms of hard (i.e. high-energy) X-ray generation will likely make important contributions toward the use of X-pinches for radiography through the harnessing or mitigation of hard X-rays, depending on the application. The modeling of electrode surface plasmas, including vacuum electron flow, will improve our understanding of power loss to various loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

External Magnetic Field Effects on Ablation of Current-Driven Foils Using an Extended Magnetohydrodynamics Simulation

We numerically model the ablation process of a 25-μm-thick aluminum foil driven by a pulsed-power machine that provides a 1-MA peak current in a 100-ns zero-to-peak rise time. The extended magnetohydrodynamics simulation is a discontinuous Galerkin code with Cartesian coordinates in 3-D and with 25-μm spatial resolution. We investigate here the influence of an external magnetic field normal to the foil surface, B z . During the foil ablation, B z = 1 T causes more nonuniform distributions of density and current compared to B z = 0 T. B z = 4 T delays the generation of surface plasma relative to the 0- and 1-T cases. The understanding of a material's ablation as it undergoes transition from the solid to plasma phases requires detailed knowledge of a material's equation of state and conductivity. This paper of warm dense matter and how instabilities propagate from a solid material to plasma motivates improvements to both numerical simulations and experimental diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A mathematical model of the structure and evolution of small-scale discrete auroral arcs

A three-dimensional fluid model for the structure and evolution of small-scale discrete auroral arcs originating from Alfven waves is developed and used to study the nonlinear macroscopic plasma dynamics of these auroral arcs. The results of simulations show that stationary auroral arcs can be unstable to a collisionless tearing mode which may be responsible for the observed transverse structuring in the form of folds and curls. At late times, the plasma becomes turbulent having transverse electric field power spectra that tend toward a universal k exp -5/3 spectral form.

Seyler, Charles E.↗

Collisional electrostatic ion cyclotron waves as a possible source of energetic heavy ions in the magnetosphere

A new mechanism is proposed for the source of energetic heavy ions (NO/+/, O2/+/, and O/+/) found in the magnetosphere. Simulations using a multispecies particle simulation code for resistive current-driven electrostatic ion cyclotron waves show transverse and parallel bulk heating of bottomside ionospheric heavy ion populations. The dominant mechanism for the transverse bulk heating is resonant ion heating by wave-particle ion trapping. Using a linear kinetic dispersion relation for a magnetized, collisional, homogenous, and multiion plasma, it is found that collisional electrostatic ion cyclotron waves near the NO(+), O2(+), and O(+) gyrofrequencies are unstable to field-aligned currents of 50 microA/sq m for a typical bottomside ionosphere.

Providakes, Jason↗

Particle and fluid simulations of resistive current-driven electrostatic ion cyclotron waves

The results from 1-D numerical simulations of electrostatic ion cyclotron waves (EIC) are presented for a model in which the electrons are a resistive (collisional) fluid. Simulations of both the kinetic and fluid descriptions are performed and compared in order to assess the fundamental limitations of fluid theory for EIC waves. The effect of ion-neutral collisions is also included using a simple Monte Carlo technique. It is found that a small ion-neutral collision frequency destroys the frequency harmonic coupling of kinetic EIC waves and tends to validate the fluid description. The saturation amplitude of the current-driven EIC instability is in agreement with recent laboratory experiments. The coherent nature (extremely narrow spectral width) and phase velocity agree with ground based (coherent backscatter radars) and in situ observations of current-driven EIC waves in the high latitude ionosphere.

Seyler, Charles E.↗