SEARCH · Search NASA
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.
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.
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 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.
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.
This project has been a multiyear effort to assess the feasibility of a key process inherent to virtually all fusion propulsion concepts: the expansion of a fusion-grade plasma through a diverging magnetic field. Current fusion energy research touches on this process only indirectly through studies of plasma divertors designed to remove the fusion products from a reactor. This project was aimed at directly addressing propulsion system issues, without the expense of constructing a fusion reactor. Instead, the program designed, constructed, and operated a facility suitable for simulating fusion reactor grade edge plasmas, and to examine their expansion in an expanding magnetic nozzle. The approach was to create and accelerate a dense (up to l0(exp 20)/m) plasma, stagnate it in a converging magnetic field to convert kinetic energy to thermal energy, and examine the subsequent expansion of the hot (100's eV) plasma in a subsequent magnetic nozzle. Throughout the project, there has been a parallel effort between theoretical and numerical design and modelling of the experiment and the experiment itself. In particular, the MACH2 code was used to design and predict the performance of the magnetoplasmadynamic (MPD) plasma accelerator, and to design and predict the design and expected behavior for the magnetic field coils that could be added later. Progress to date includes the theoretical accelerator design and construction, development of the power and vacuum systems to accommodate the powers and mass flow rates of interest to out research, operation of the accelerator and comparison to theoretical predictions, and computational analysis of future magnetic field coils and the expected performance of an integrated source-nozzle experiment.
A device for the production of low energy (3 to 10 eV) neutral atomic beams for surface modification studies is described that reproduces the flux of atomic oxygen in low Earth orbit. The beam is produced by the acceleration of plasma ions onto a negatively biased plate of high-Z metal; the ions are neutralized and reflected by the surface, retaining some fraction of their incident kinetic energy, forming a beam of atoms. The plasma is generated by a coaxial RF exciter which produces a magnetically-confined (4 kG) plasma column. At the end of the column, ions fall through the sheath to the plate, whose bias relative to the plasma can be varied to adjust the beam energy. The source provides a neutral flux approximately equal to 5 x 10(exp 16)/sq cm at a distance of 9 cm and a fluence approximately equal to 10(exp 20)/sq cm in five hours. The composition and energy of inert gas beams was diagnosed using a mass spectometer/energy analyzer. The energy spectra of the beams demonstrate energies in the range 5 to 15 eV, and qualitatively show expected dependences upon incident and reflecting atom species and potential drop. Samples of carbon film, carbon-based paint, Kapton, mylar, and teflon exposed to atomic O beams show erosion quite similar to that observed in orbit on the space shuttle.
The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.
Introduction: An object immersed in an airless plasma environment will experience a natural process of surface charging in order to acheieve current balance, or zero net electric current to the object. It has been shown in recent computer simulations that the small-body plasma environment is very complex [1], considering effects of photoemission, topography, and formation of a plasma wake. For this work we consider an exploration craft (or astronaut) immersed within a plasma environment near an asteroid, which exhibits widely varying solar wind and photoelectric particle fluxes and continuously evolving illumination conditions. Objective: We aim to determine how an explo-ration craft or astronaut suit accumulates charge while located in the "nightside" asteroid wake where the particle fluxes are reduced, and in the dayside near-surface photoelectron sheath, by combining an object charging model [2] with kinetic simulations of a near-asteroid plasma environment [1]. We consider an astronaut floating near the asteroid while not in contact with the surface, as well as an astronaut moving along the surface using their hands/gloves to crawl along. Results: The modeling results suggest that remediation of triboelectric charge via accumulation of plasma currents is an important factor to consider when designing future NEA mission infrastructure, especially if repeated and frequent contact with the surface is planned. In shadowed regions such as the location shown in Fig. 1a, the plasma currents are so low (and the effective charge-remediation timescale so long, e.g. minutes to hours) that repeated contact with the surface tribocharges the glove in an uncontrollable fashion, as shown for two representative electron temperatures in Fig. 2a. The resulting buildup of significant negative charge would eventually initiate some other "current of last resort" [4] such as transport of positively-charged dust, field-emission from the glove, or significant alteration of environmental ion currents within the wake. In contrast, the few-meters-thick dayside photoelectron sheath in which the astronaut of Fig. 1b is immersed in is so rich in electrons (and hence so electrically conductive) that accumulated tribocharge dissipates almost instantaneously (e.g. in less than a ms) as shown in Fig. 2b. As our model astronaut orbits the NEA they would experience plasma currents and associated charge re-mediation times spanning many orders of magnitude, and the fusion between our numerical models provides a detailed understanding of the charging hazards possibly associated with contact-based NEA exploration.
A new experimental approach is discussed that can be adopted for studies in the area of plasma flow interactions with bodies in space. The potential use of the Space Shuttle/Orbiter as a near-earth plasma laboratory for studies in space plasma physics and particularly in solar system plasmas is discussed. This new experimental approach holds great promise for studies in the supersonic and sub-Alfvenic flow regime which has applications to the motion of natural satellites around their mother planets in the solar-system (e.g., the satellite Io around the planet Jupiter). A well conceived experimental and theoretical program can lead to a better physical understanding regarding the validity and range of applicability of using gasdynamic, kinetic, and fluid approaches in describing collisionless plasma flow interactions with bodies in a variety of flow regimes. In addition to the above scientific aspects of the program, significant technological advances can be achieved regarding the interaction of space probes in planetary atmospheres/ionospheres and the reliability of using various plasma diagnostic devices on board spacecraft and large space platforms.
The rms field strengths of electrostatic and electromagnetic turbulence in the earth's bow shock, measured in the frequency range 20 Hz to 200 kHz with the Imp 6 satellite, are examined. The largest strengths of electrostatic turbulence occur when the upstream electron-to-proton temperature ratio is large and the proton temperature is small. No substantial correlation is found among the rms field strengths of electrostatic turbulence, the Alfven Mach number, the ratio of particle to magnetic-field pressure, and the shock normal angle. These results indicate that the strength of electrostatic turbulence in the bow shock is determined by the kinetic properties of the solar-wind plasma. The largest strengths for electromagnetic turbulence occur when the upstream particle density is large and when the shock normal angle is closer to 90 deg, supporting a previous conclusion that whistler waves comprise the electromagnetic turbulence in the bow shock.
Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.
Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in 3D, we conduct a comprehensive particle-in-cell simulation study comparing 2D and 3D evolution of long, thin current sheets in moderately-magnetized, collisionless, relativistically-hot electron- positron plasma, and find dramatic differences. We first systematically characterize this process in 2D, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths, and system sizes. We then show that 3D simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3D current-sheet evolution is not always fundamentally classical reconnection with perturbing 3D effects, but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details, and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in 2D. Intriguingly, nonthermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in 2D and 3D. Though the variety of underlying current-sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.
A number of investigations have been concerned with the kinetic theory and processes associated with a relativistic electron gas. Gould (1981) has considered a condition in which upon the ultimate thermalization the temperature can be such that the electron gas is highly relativistic while the gas of protons and other ions is nonrelativistic. With the nuclear component nonrelativistic but having energies in the MeV range and above, it is necessary to consider the effects of nuclear forces in the scattering of the ions in their thermalization. The effects of nuclear forces in the thermalization of ions in plasmas have been computed, principally in connection with problems of controlle; fusion. The present investigation is concerned with an attempt to express results in analytic form to as great a degree as possible. The p-p problem, which is the fundamental problem in astrophysical plasma, is studied. Attention is given to a low-energy formulation, the s-wave phase shift, the effective stopping number, Fokker-Planck operators, and the interaction with the electron gas.
One approach that has been investigated for trash management aboard long-duration crewed exploration space missions is the so-called Trash-to-Gas approach. This technology uses a thermal process (i.e., combustion, pyrolysis, etc.) to convert crew trash into relatively non-reactive gases (i.e., CO2, CO, etc.) that can be easily vented overboard. One difficulty with this approach is the potential for the contamination of the spacecraft habitable volume from the myriad volatile organic compounds (VOCs) that are concurrently produced from these thermal processes. One possible VOC abatement process is to apply a plasma to the effluent gases. Plasmas provide high energy electrons and ions that are capable of breaking down gaseous organic species generally into smaller compounds. To explore Plasma VOC abatement, the authors constructed a plasma gas system that recirculates a particular VOCladen gas mixture through a plasma torch. This allows the periodic collection of aliquots of gaseous samples for subsequent analyses via gas chromatography mass spectrometry, designed to quantify the remaining VOC materials. The VOCs evaluated included 100 ppm of acetone, benzene, ethanol, and pentane in a background of carbon dioxide. The plasma conversion of all analytes followed first order rate kinetics, and complete VOC elimination was observed after 10 min or less of treatment, for all analytes measured. Additionally, it was observed that the plasma caused the conversion of CO2 into oxygen and CO.
Wave-particle transport in a collisionless plasma is due to particle scattering by enhanced fluctuations associated with the growth of instabilities. In particular, relatively short wave-length kinetic instabilities are frequently invoked to explain many different types of plasma transport in space. Although there is an extensive theoretical and simulation literature describing the potential applications of many such instabilities, there are only a few cases of clear-cut identification of kinetic modes in space. The research described in this paper uses linear Vlasov dispersion theory to study correlation functions and dimensionless correlation function ratios for fluctuations or instabilities in three space plasma regimes. This research shows that both the compressibility and the parallel compressibility are likely to be useful in distinguishing modes in the magnetosheath as well as in the plasma sheet boundary layer and that helicity remains a useful identifier of electromagnetic ion/ion instabilities in the foreshock.
This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.
The outflow of coronal plasma into interplanetary space is a consequence of the coronal heating process. Therefore the formation of the corona and the acceleration of the solar wind should be treated as a single problem. Traditionally the mass or particle flux emanating from the extended corona has been thought of as being determined by the coronal temperature or scale height and the coronal (base) density. This argument follows from considerations of the momentum balance of the corona-wind system from which one obtains models of a close to hydrostatic corona out to the critical point where the flow becomes supersonic. With this approach to the acceleration of the wind is has been difficult to reconcile the relatively small variation observed in the proton flux at 1 AU with the predicted exponential dependence of the proton flux on the coronal temperature. In this talk we would like to emphasize another approach in which coronal energetics play the primary role. The deposition of energy into the corona through some 'mechanical' energy flux is balanced by the various energy sinks available to the corona and the sum of these processes determine the coronal structure, i.e. its temperature and density. The corona loses energy through heat conduction into the transition region, through radiative losses, and through the gravitational potential energy and kinetic energy put into the solar wind itself. We will show from a series of models of the chromosphere transition region-corona-solar wind system that most of the energy deposited in a magnetically open region will go into the solar wind, with roughly half going into kinetic energy and half into lifting the plasma out of the solar gravity field. The coronal base density will adjust itself in such a way that the heat conductive flux flowing into the transition region is radiated away in the upper chromosphere. The coronal temperature is set by the requirements that most of the deposited energy goes into accelerating the solar wind; the coronal scale height will adjust itself so that the solar wind energy losses conform to the amplitude of the input energy. These processes are modified by the 'mode' of energy deposition, and we will show the effects on coronal structure of changing the parameters describing coronal heating as well as the effects of including a helium fluid in the models. However, the location, scale height and/or form of the energy deposition (i.e. heating or direct acceleration) are not too important for the solar wind, the coronal density and temperature structure will vary with the 'mode' of energy deposition, but the solar wind mass flux depends mainly on the amplitude of the energy flux.
The objective is to use the measurements of the equatorial particle distributions to obtain the parallel electric field structure and the evolution of the plasma distribution function along the field line. Appropriate uses of kinetic theory allows us to use the measured ( and inferred) particle distributions to obtain the electric field, and hence the variation on plasma density along the magnetic field line. The approach, here, is to utilize the adiabatic invariants, and assume the plasma distributions are in equilibrium.