Capillary instability of a jet carrying an axial current.
Hall current effect on capillary instability of finitely conducting liquid jet carrying axial volume current
SEARCH · Engineering Papers
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.
Hall current effect on capillary instability of finitely conducting liquid jet carrying axial volume current
Electrodeless MGD power generation during entry, using torque on magnet taken as rotating member of generator
FLEXO (Flux-Limited Extended-MHD Ohm's Law) is a production-line multiphysics code developed at Sandia to enable more predictive modeling of target physics on pulsed-power devices. FLEXO uses an extended magnetohydrodynamics (XMHD) model which includes a generalized Ohm's law (GOL), an electron inertia term, and Hall physics. This report describes the code's numerical methods, its computational performance, and test problems of interest.
Magnetohydrodynamic power generation experiments utilizing a cesium-seeded H2-O2 working fluid were carried out using a diverging area Hall duct having an entrance Mach number of 2. The experiments were conducted in a high-field strength cryomagnet facility at field strengths up to 5 tesla. The effects of power takeoff location, axial duct location within the magnetic field, generator loading, B-field strength, and electrode breakdown voltage were investigated. For the operating conditions of these experiments, it is found that the power output increases with the square of the B-field and can be limited by choking of the channel or interelectrode voltage breakdown which occurs at Hall fields greater than 50 volts/insulator. Peak power densities of greater than 100 MW/cu M were achieved.
In a qualification life test of a Hall thruster it was found that the erosion of the acceleration channel practically stopped after approx 5,600 h. Numerical simulations using a two-dimensional axisymmetric plasma solver with a magnetic field-aligned mesh reveal that when the channel receded from its early-in-life to its steady-state configuration the following changes occurred near the wall: (1) reduction of the electric field parallel to the wall that prohibited ions from acquiring significant impact kinetic energy before entering the sheath, (2) reduction of the potential fall in the sheath that further diminished the total energy ions gained before striking the material, and (3) reduction of the ion number density that decreased the flux of ions to the wall. All these changes, found to have been induced by the magnetic field, constituted collectively an effective shielding of the walls from any significant ion bombardment. Thus, we term this process in Hall thrusters "magnetic shielding."
Planetary EDL requires that a spacecraft must reduce its relative velocity from multiple km/s to zero for safe landing. The entry phase requires slowing down through an intense hypersonic flight regime (Mach 25 to Mach 5). Magnetohydrodynamics (MHD) is a next-generation concept with a wide untapped set of technologies that can expand planetary exploration capabilities and create more robust entry vehicles by utilizing in-situ resources, increasing control authority, and reducing heat flux. There is therefore a need for an experimentally validated simulation tool for MHD technologies in hypersonic entry.
The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. Results obtained indicate that the cathode sheath can exert significant resistance to current flow and thus reduce the MHD-generated force by 40% or more. The effect of the sheath on the MHD force is found to be more pronounced at lower applied magnetic field.
Laser and optical emission-based measurements have been developed and implemented for use on low-power hydrogen arcjet thrusters and xenon-propelled electric thrusters. In the case of low power hydrogen arcjets, these laser induce fluorescence measurements constitute the first complete set of data that characterize the velocity and temperature field of such a device. The research performed under the auspices of this NASA grant includes laser-based measurements of atomic hydrogen velocity and translational temperature, ultraviolet absorption measurements of ground state atomic hydrogen, Raman scattering measurements of the electronic ground state of molecular hydrogen, and optical emission based measurements of electronically excited atomic hydrogen, electron number density, and electron temperature. In addition, we have developed a collisional-radiative model of atomic hydrogen for use in conjunction with magnetohydrodynamic models to predict the plasma radiative spectrum, and near-electrode plasma models to better understand current transfer from the electrodes to the plasma. In the final year of the grant, a new program aimed at developing diagnostics for xenon plasma thrusters was initiated, and results on the use of diode lasers for interrogating Hall accelerator plasmas has been presented at recent conferences.
We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.
Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results
The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. The effect of the electrode material on the MHD process is studied. Using thoriated tungsten instead of graphite leads to a four-fold increase in the Lorentz forces and also leads to significantly reduced heat fluxes on the cathode. This is attributed to the much higher electrical conductivity of the thoriated tungsten sheath reducing by orders of magnitude the plasma electrical resistance near the surfaces.
Non-axisymmetric, flow-driven instabilities in the incompressible Hall-MHD model are studied in a differentially rotating cylindrical plasma. It is found that, in the Hall-MHD regime, both whistler waves and ion-cyclotron waves can extract energy from the flow shear, resulting in two distinct branches of global instability. The non-axisymmetric whistler modes grow significantly faster than non-axisymmetric, ideal MHD modes. A discussion of the global whistler instability mechanism is presented in the large-ion-skin-depth, “electron-MHD” limit. When the magnetic field is azimuthal, a subset of the whistler modes having zero axial wave number are uncovered to be destabilized by the “corotation amplifier” mechanism. It is observed that the effect of the Hall term on the non-axisymmetric modes can be appreciable when d i is on the order of a few percent of the width of the cylindrical annulus. Distinct global modes emerge in the strong Hall-MHD regime at significantly stronger magnetic fields than those required for unstable global MHD modes, as the Hall effect weakens the stabilizing “field-line bending” by decoupling ion motion from the magnetic field. These global non-axisymmetric modes may play an important role in weakly ionized accretion disks.
Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground-based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected in auroral, ionospheric, and ground magnetic observations through both magnetohydrodynamic theory and a global coupled magnetosphere-ionosphere simulation of a magnetopause surface eigenmode. These show monochromatic oscillatory field-aligned currents (FACs), due to both the surface mode and its non-resonant Alfvén coupling, are present throughout the magnetosphere. The currents peak in amplitude at the equatorward edge of the magnetopause boundary layer, not the open-closed boundary as previously thought. They also exhibit slow poleward phase motion rather than being purely evanescent. We suggest the upward FAC perturbations may result in periodic auroral brightenings. In the ionosphere, convection vortices circulate the poleward moving FAC structures. Finally, surface mode signals are predicted in the ground magnetic field, with ionospheric Hall currents rotating perturbations by approximately (but not exactly) 90° compared to the magnetosphere. Thus typical dayside magnetopause surface modes should be strongest in the East-West ground magnetic field component. Overall, all ground-based signatures of the magnetopause surface mode are predicted to have the same frequency across L-shells, amplitudes that maximize near the magnetopause's equatorward edge, and larger latitudinal scales than for field line resonance. Implications in terms of ionospheric Joule heating and geomagnetically induced currents are discussed.
We examine the role that ions and electrons play in reconnection using observations from the Magnetospheric Multiscale (MMS) mission on kinetic ion and electron scales, which are much shorter than magnetohydrodynamic scales. This study reports observations with unprecedented high resolution that MMS provides for magnetic field (7.8 ms) and plasma (30 ms for electrons and 150 ms for ions). We analyze and compare approaches to the magnetopause in 2016 November, to the electron diffusion region in the magnetotail in 2017 July followed by a current sheet crossing in 2018 July. Besides magnetic field reversals, changes in the direction of the flow velocity, and ion and electron heating, MMS observed large fluctuations in the electron flow speeds in the magnetotail. As expected from numerical simulations, we have verified that when the field lines and plasma become decoupled a large reconnecting electric field related to the Hall current (1–10 mV/m) is responsible for fast reconnection in the ion diffusion region. Although inertial accelerating forces remain moderate (1–2 mV/m), the electric fields resulting from the divergence of the full electron pressure tensor provide the main contribution to the generalized Ohm’s law at the neutral sheet (as large as 200 mV/m). In our view, this illustrates that when ions decouple electron physics dominates. The results obtained on kinetic scales may be useful for better understanding the physical mechanisms governing reconnection processes in various magnetized laboratory and space plasmas.
Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.
The Magnetized Liner Inertial Fusion (MagLIF) experimental platform at Sandia National Laboratories has realized a number of recent improvements in electrical current and laser preheat coupling but fusion yields are below expectations, based on modeling with clean (that is, simulations with no interfacial instabilities and no mix) two-dimensional (2D) magneto-hydrodynamics (MHD) simulations, by a factor of 5 or more. However, magnetized liner implosions are known to produce helical magneto-Rayleigh–Taylor (MRT) modes and complex stagnation structures that cannot be directly modeled in 2D. This paper presents the results of 3D HYDRA MHD simulations, including helical MRT, showing that degradation from these instabilities can readily reproduce experimental yields, but it is difficult to simultaneously match the yield along with inferred fuel temperature, pressure, and burn history (assessed with x-rays), which are larger in the experiments considered here. Additional analysis methods and future experiments are proposed to help address the discrepancies. The 3D simulations also show, without substantive improvements to stability, changing the applied Bz or preheat is unlikely to increase performance of MagLIF to the same degree as clean 2D simulations. Finally, the first 3D HYDRA simulations with the Hall term show the ab initio production of helical MRT and produce a clear change in stagnation morphology compared to pre-seeded simulations without Hall. In spite of the differences, the simulations still produce comparable fusion performance but with longer, lower power, neutron yield history.
The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.
The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.