SEARCH · Engineering Papers
Results for “Kinetic plasmas”
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.
A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges
Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.
VPIC-Kokkos
VPIC is a particle-in-cell (PIC) code that solves the Maxwell-Boltzmann system of equations and can be used to simulate a wide variety of plasmas, including those that are fully kinetic and relativistic. Utilizing the Kokkos performance-portable framework, VPIC achieves high performance on multiple CPU and GPU architectures and is adaptable to future platforms with minimal developer effort. VPIC features very powerful input decks, allowing insertion of arbitrary C++ code for custom diagnostics, boundary conditions, and additional physics models.
Collision-induced high frequency instability in a fully ionized plasma.
Damping of high frequency waves in homogeneous fully ionized plasma with relative electron-ion drift calculated, using fokker-planck kinetic equation
Electric conductivity of weakly turbulent plasmas
Weakly turbulent plasmas static electric conductivity derivation from kinetic equation for linear response to one-particle distribution function
Electromagnetic processes in the atmosphere of pulsars
The work consists of two parts. The first deals with the fine structure of radio pulses. Based on kinetic theory, processes occurring in the plasma shell of a pulsar when external electromagnetic radiation is present are investigated. It is shown that electromagnetic waves cause electrons to drift relative to ions, and initiate longitudinal oscillations. A dispersion equation describing the longitudinal oscillations in magnetized plasma is derived. Conditions for excitation of oscillations are found. Correlation functions of electron density are calculated, along with the coefficients of electromagnetic wave scattering. It is shown that variations in the amplitude of pulsar pulses are associated with scintillations caused by fluctuations in the plasma electron density. The second part of the study presents a mechanism for the radio emission of pulsars. The model of a rotating and a pulsating star, a neutron star with dipolar or more complex magnetic field, is examined.
On ionospheric aerodynamics
This paper presents theoretical methods to determine the gas dynamic and the electrostatic effects due to the interaction caused by a rapidly moving body in the ionosphere. The principles of the methods are derived from the kinetic theory of collision-free plasma. It is shown that the collective behavior of the collision-free plasma makes it possible to use the fluid approach to treat the problems of ionospheric aerodynamics. Various solutions to the system of fluid and field equations that have direct bearing on the ionospheric aerodynamics are presented and discussed. Physical significances of the mathematical results are stressed. Some outstanding unsolved problems in ionospheric aerodynamics are elaborated and discussed.
A theory of local and global processes which affect solar wind electrons. 2: Experimental support
The microscopic characteristics of the Coulomb cross section show that there are three natural subpopulations for plasma electrons: the subthermals with local kinetic energy E kT sub c; the transthermals with kT sub c E 7 kT sub c and the extrathermals E 7 kT sub c. Data from three experimental groups on three different spacecraft in the interplanetary medium over a radial range are presented to support the five interrelations projected between solar wind electron properties and changes in the interplanetary medium: (1) subthermals respond primarily to local changes (compression and rarefactions) in stream dynamics; (2) the extrathermal fraction of the ambient electron density should be anti-correlated with the asymptotic bulk speed; (3) the extrathermal "temperature" should be anti-correlated with the local wind speed at 1 AU; (4) the heat flux carried by electrons should be anti-correlated with the local bulk speed; and (5) the extrathermal differential 'temperature' should be nearly independent of radius within 1 AU.
A predicted energetic ion signature of a neutral line in the geomagnetic tail
Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.
Pulsed plasmoid electric propulsion
A method of electric propulsion is explored where plasmoids such as spheromaks and field reversed configurations (FRC) are formed and then allowed to expand down a diverging conducting shell. The plasmoids contain a toroidal electric current that provides both heating and a confining magnetic field. They are free to translate because there are no externally supplied magnetic fields that would restrict motion. Image currents in the diverging conducting shell keep the plasmoids from contacting the wall. Because these currents translate relative to the wall, losses due to magnetic flux diffusion into the wall are minimized. During the expansion of the plasma in the diverging cone, both the inductive and thermal plasma energy are converted to directed kinetic energy producing thrust. Specific impulses can be in the 4000 to 20000 sec range with thrusts from 0.1 to 1000 Newtons, depending on available power.
Semikinetic and generalized transport models of the polar and solar winds
In order to gain a better understanding of the relative merits of the transport and kinetic approaches to modeling thermal space plasmas, a comparison is presented of a transport (bi-Maxwellian-based 16-moment equations) and a semikinetic description of supersonic flow in the solar wind and also of both supersonic and subsonic flows in the polar wind for 'steady state' conditions. The study shows: remarkable agreement between the two models for supersonic collisionless flows, even for the higher-order moments; the inadequacy of the semikinetic approach for modeling subsonic flows; and the superiority of the 16-moment transport over the semikinetic approach for modeling the solar wind. Further evidence that the bi-Maxwellian-based transport equations are a useful tool for studying 'thermal' space plasmas that develop non-Maxwellian features is provided.
GRO: Black hole models for gamma-ray bursts
The possibility of creating gamma ray bursts (GRB's) from accretion flows on to black holes is investigated. The mechanism of initial energy release in the form of a burst is not understood yet. The typical time scales involved in this energy release and the initial distribution of photons as a function of energy are studied. As a first step the problem is formulated in the Minkowski spacetime for a homogeneous and isotropic burst. For an arbitrary initial distribution of photons, the equations of relativistic kinetic theory are formulated for nonequilibrium plasmas which can take into account various particle creation and annihilation processes and various scattering processes.
Quark and Gluon Relaxation in Quark-Gluon Plasmas
The quasiparticle decay rates for quarks and gluons in quark-gluon plasmas are calculated by solving the kinetic equation. Introducing an infrared cutoff to allow for nonperturbative effects, we evaluate the quasiparticle lifetime at momenta greater than the inverse Debye screening length to leading order in the coupling constant.
Entropy Generation Across Earth's Bow Shock
Earth's bow shock is a transition layer that causes an irreversible change in the state of plasma that is stationary in time. Theories predict entropy increases across the bow shock but entropy has never been directly measured. Cluster and Double Star plasma experiments measure 3D plasma distributions upstream and downstream of the bow shock that allow calculation of Boltzmann's entropy function H and his famous H-theorem, dH/dt ~ O. We present the first direct measurements of entropy density changes across Earth's bow shock. We will show that this entropy generation may be part of the processes that produce the non-thermal plasma distributions is consistent with a kinetic entropy flux model derived from the collisionless Boltzmann equation, giving strong support that solar wind's total entropy across the bow shock remains unchanged. As far as we know, our results are not explained by any existing shock models and should be of interests to theorists.
The Fusion Driven Rocket: Nuclear Propulsion through Direct Conversion of Fusion Energy
The future of manned space exploration and development of space depends critically on the creation of a dramatically more efficient propulsion architecture for in-space transportation. A very persuasive reason for investigating the applicability of nuclear power in rockets is the vast energy density gain of nuclear fuel when compared to chemical combustion energy. The Fusion Driven rocket (FDR) represents a revolutionary approach to fusion propulsion where the power source releases its energy directly into the propellant, not requiring conversion to electricity. It employs a solid lithium propellant that requires no significant tankage mass. The propellant is rapidly heated and accelerated to high exhaust velocity (> 30 km/s), while having no substantial physical interaction with the spacecraft thereby avoiding damage to the rocket and limiting both the thermal heat load and radiator mass. The key to achieving this stems from research at MSNW and the UW on the magnetically driven implosion of metal foils onto a magnetized plasma target to obtain fusion conditions. A logical extension of this work leads to a method that utilizes these metal shells (or liners) to not only achieve fusion conditions, but to serve as the propellant as well. Several low-mass, magnetically driven metal liners are inductively driven to converge radially and axially and form a thick blanket surrounding the target plasmoid and compress the plasmoid to fusion conditions. Virtually all of the radiant, neutron and particle energy from the plasma is absorbed by the encapsulating, thick metal blanket thereby isolating the spacecraft from the fusion. This energy, in addition to the intense Ohmic heating at peak magnetic field compression, is adequate to vaporize and ionize the metal blanket. The expansion of this hot, ionized metal propellant through a magnetically insulated nozzle produces high thrust at the optimal Isp. The energy from the fusion process, along with the waste heat, is thus utilized at very high efficiency. The basic scheme for FDR is illustrated and described in the report (see Fig. 2) The two most critical issues in meeting challenges introduced employing magneto-inertial fusion as the power source is driver efficiency and “stand-off” – the ability to isolate and protect fusion and thruster from the resultant fusion energy. By employing metal shells for compression, it is possible to produce the desired convergent motion inductively by inserting the metal sheets along the inner surface of cylindrical or conically tapered coils. Both stand-off and energy efficiency issues are solved by this arrangement. 3 This two year effort focused on achieving three key criteria for the Fusion Driven Rocket to move forward for technological development: (1) the physics of the FDR must be fully understood and validated, (2) the design and technology development for the FDR required for its implementation in space must be fully characterized, and (3) an in-depth analysis of the rocket design and spacecraft integration as well as mission architectures enabled by the FDR need to be performed. A subscale, laboratory liner compression test facility was assembled at the University of Washington Plasma Dynamics Laboratory with sufficient liner kinetic energy (~ 0.5 MJ) to reach conditions required for fusion breakeven conditions. Detailed experimental studies of the dynamic behavior of the driven liners as well as liner convergence and magnetic compression were performed. The development of both the 1D liner dynamics code and the full 3D ANSYS® liner calculations was achieved. The characterization of both the FDR and spacecraft as well as a design architecture analysis was conducted that included an examination of a wide range of mission architectures and destinations for which this fusion propulsion system would be enabling or critical. In particular a rapid, single launch manned Mars mission was developed.
Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification
A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.
1D EDIPIC
EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.
EDIPIC-2D
EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.